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: general theories of open quantum dynamics as introduced in @xcite provide the mathematical pathway to the characterization of real - world quantum mechanical systems , subject to dissipation and dephasing by environmental interactions . such effects are crucial across a multitude of fields ranging from solid - state to chemical physics , quantum optics , and mesoscopic physics .    in the context of a classical environment , linear dissipation can be modeled by the caldeira - leggett oscillator approach @xcite , closely linked to langevin equations @xcite which offer a concise formalism based on retarded friction kernels and gaussian random forces ( thermal noise ) . the quantum analogue of friction , however , needs a much more subtle treatment since it typically creates substantial correlations between system and environment . moreover , quantum fluctuations of a thermal reservoir are non - zero at any temperature , leading to interesting phenomena and non - trivial ground states .    within the quantum regime the dynamical properties of the reduced density matrix are paramount . by tracing out reservoir degrees of freedom from the global dynamics , the focus is narrowed towards a relevant subsystem . any systematic treatment of dynamical features like decoherence of quantum states , dissipation of energy , relaxation to equilibrium or non - equilibrium steady states requires a consistent procedure to distill a dynamical map or an equation of motion for the reduced density matrix from the unitary evolution of system and reservoir . the markovian approximation typically generalizes the classical probabilistic technique of a dynamical semigroup in analogy to the differential chapman - kolmogorov equation @xcite . the finite - dimensional mathematical framework of quantum dynamical semigroups traces back to the seminal work by gorini , kossakowski , sudarshan @xcite and simultaneously by lindblad @xcite . while the resulting quantum master equations of lindblad form provide an easy - to - use set of tools for many applications , their perturbative nature fails in the presence of strong environment coupling , long correlation time scales or entanglement in the initial state . the degree of non - markovianity which is inherent to the density matrix evolution and which causes pronounced retardation effects in reservoir - mediated self - interactions constitutes an active field of research @xcite . whereas considerable advances have been made in the characterization of dynamical generators as non - markovian @xcite , much less attention is paid to the question of how non - markovian behavior arises from the hamiltonian description of a system - reservoir model .    beyond perturbative dynamics of memoryless master equations and related strategies like quantum jumps @xcite or quantum state diffusion @xcite , a pool of numerically exact simulation methods has been developed , each with specific strengths and specific weaknesses . one can distinguish between methods set - up in the full hilbert space of system and bath degrees of freedom and those considering the dynamics of the reduced density operator of the system alone . the former include approaches based on e.g.  the numerical renormalization group ( nrg ) @xcite , the multiconfiguration hartree ( mctdh ) @xcite , and the density matrix renormalization group ( dmrg ) @xcite . the latter can all be derived from the path integral formulation pioneered by feynman and vernon @xcite . they treat the functional integration either directly such as the path integral quantum monte carlo ( pimc ) @xcite and the quasi - adiabatic propagator ( quapi ) @xcite or cast it in some form of time evolution equations . this is by no means straightforward due to the bath induced time retardation , a problem that always appears at lower temperatures . equivalence to the path integral expression is then only guaranteed for a nested hierarchy of those equations @xcite or time evolution equations carrying stochastic forces @xcite from which the reduced density follows after a proper averaging . these stochastic approaches exploit the intimate connection between the description of a quantum reservoir in terms of an influence functional and stochastic processes . in fact , influence functionals do not only arise when a partial trace is taken over environmental degrees of freedom , they are also representations of random forces sampled from a classical probability space @xcite . this stochastic construction can be reversed , leading to an unraveling of quantum mechanical influence functionals into time - local stochastic action terms @xcite ; we thus obtain the dynamics of the reduced system through statistical averaging of random state samples generated by numerically solving a single time - local stochastic liouville - von neumann equation ( sln @xcite ) . compared to other methods this provides a very transparent formulation of non - makrovian quantum dynamics with the particular benefit that the consistent inclusion of external time dependent fields is straightforward @xcite . since the random forces resulting from the exact mapping of the quantum reservoir to a probability space are not purely real , the resulting stochastic propagation is non - unitary . therefore , the signal - to - noise ratio of empirical statistics based on sln propagation deteriorates for long - time . this issue is reminiscent of the sign problem in real - time path integrals @xcite and represents a major hurdle in solving the most general c - number noise stochastic liouville - von neumann equation for time intervals much longer than the timescales of relaxation and dephasing . here we modify a strategy , recently presented by one of us @xcite , which uses a projection operator @xcite based finite - memory scheme to solve the complex - noise sln numerically . finite - memory stochastic propagation ( fmsp ) significantly lowers the effect of statistical fluctuations and leads to a much faster convergence of long - time sample trajectories , with a gain in computational efficiency by several orders of magnitude . while the relevant memory timescales used before were reservoir correlation times @xcite , we use a different projector here , adapted to the case of finite dephasing times . if either the dissipative coupling or the reservoir temperature exceed certain thresholds , this results in a shorter memory time and better statistics . we apply this framework to two- and three - level systems and compare numerical data with perturbative predictions . particular emphasis is put on the transition from coherent to incoherent population dynamics . since the new method allows for converged simulations also in the long time domain , transfer rates for sequential hopping and super - exchange @xcite can be extracted which are of relevance for charge or energy transfer in molecular aggregates and arrays of artificial atoms , e.g. quantum dot structures . the paper is organized as follows : we start in sec . [ sln ] with a concise discussion of the sln and its simplified version for ohmic spectral densities . the new scheme based in projection operator techniques is introduced in sec . [ tcbd ] , before the two - level system in sec . [ twolevel ] and the three - level structure in sec . [ threelevel ] are analyzed . we consider a distinguished system which is embedded in an environment with a large number of degrees of freedom . the hamiltonian of such a model comprises a system , a reservoir and an interaction term @xmath0 for bosonic elementary excitations of the reservoir we assume @xmath1 together with a bilinear coupling part @xmath2 that links the system coordinate @xmath3 to the bath force @xmath4 . from the unitary time evolution of the global density matrix @xmath5 that belongs to the product space @xmath6 , we recover the reduced density matrix @xmath7 by a partial trace over the reservoir s degrees of freedom @xmath8 and a factorizing initial condition @xmath9 . we thereby assume the reservoir to be initially in thermal equilibrium , @xmath10 . while traditional open system techniques focus on a perturbative treatment in the interaction hamiltonian @xmath11 , we derive our stochastic approach from influence functionals , a path integral concept introduced by feynman and vernon @xcite , applying to reservoirs with gaussian fluctuations of the force field @xmath12 , @xmath13\pazocal a^{\ast}[q'(\tau ' ) ] \pazocal f [ q(\tau)q'(\tau ' ) ] . \label{fvdoubint}\ ] ] here @xmath14 is the probability amplitude for paths governed by the system action alone ( a pure phase factor ) , and @xmath15 & = & - \frac{1}{\hbar^2 } \int_0^t d\tau \int_0^\tau d\tau ' ( q(\tau)-q'(\tau))\nonumber\\ & \times&[l(\tau-\tau')q(\tau ' ) - l^*(\tau-\tau')q'(\tau ' ) ] . \label{eq : lnf}\end{aligned}\ ] ] all effects of the dissipative environment on the propagation of the system can thus be expressed in terms of the _ free _ fluctuations of the reservoir , characterized by the complex force correlation function @xmath16}{\sinh(\hbar\beta\omega/2 ) } . \label{lcorr}\end{aligned}\ ] ] here we have introduced a spectral density @xmath17 , which can be determined microscopically from the frequencies and couplings of a quasi - continuum of reservoir modes , or , alternatively , from the fourier transform of the imaginary part of eq . ( [ lcorr ] ) when the correlation function has been obtained by other means . the latter definition is somewhat more general than the model of an oscillator bath . it is to be noted that the thermal timescale @xmath18 must be considered long in a quantum system with thermal energy @xmath19 smaller than the level spacing . in the opposite limit @xmath20 the classical version of the fluctuation - dissipation theorem @xcite can be recovered from eq . ( [ lcorr ] ) .    due to the non - local nature of the influence functional ( [ eq : lnf ] ) , there is no simple way to recover an equation of motion from eq . ( [ fvdoubint ] ) . however , since @xmath21 $ ] is a gaussian functional of the path variables , it shows great formal similarity to generating functionals of classical gaussian noise . a classical process formally equivalent to a quantum reservoir can thus be constructed by means of a stochastic decomposition @xcite . the resulting stochastic liouville - von neumann equation ( sln ) contains two stochastic processes @xmath22 and @xmath23 , corresponding to two independent functions of the functional @xmath24 , @xmath25\nonumber \\ & + & \frac{i}{\hbar}\xi(t)[q,\rho_{z}(t)]+\frac{i}{2}\nu(t)\{q,\rho_{z}(t)\}. \label{sln}\end{aligned}\ ] ] we thus map the reduced system evolution to stochastic propagation in probability space of gaussian noise forces with zero bias and correlations which match the quantum mechanical correlation function @xmath26 , @xmath27 it is obvious that there are no real - valued processes with these correlations , however , complex - valued stochastic processes which obey eqs . ( [ corrxixi])([corrnunu ] ) do exist .    even though eq . ( [ sln ] ) contains no term recognizable as a damping term ( in a mathematical sense ) , yet it provides an exact numerical approach to open quantum systems with any type of gaussian reservoir : stochastic samples @xmath28 obtained by propagating ( [ sln ] ) with specific noise samples have no obvious physical meaning . the _ physical _ density matrix , whose evolution is damped , is obtained by averaging the samples , @xmath29 $ ] . there is also a version of the sln which is designed for reservoirs with ohmic characteristic , i.e.  for spectral densities of the form @xmath30 up to a high frequency cutoff @xmath31 being significantly larger than any other frequency of the problem ( including the thermal time @xmath18 ) . this class of reservoirs is of particular relevance due to numerous realizations ranging from atomic to condensed matter physics . in this case , the imaginary part of the reservoir correlation function can be considered as the time derivative of a dirac @xmath32function , @xmath33 accordingly , memory effects arise only from the real part @xmath34 , while the imaginary part can be represented by a time - local damping operator acting on the reduced density matrix . this results in the simplified so - called sled dynamics ( stochastic liouville equation with dissipation ) @xcite with one real - valued noise force @xmath22 , @xmath35-\xi(t)[q,\rho_{\xi}]\right)\nonumber\\ & + & \frac{\gamma}{2i\hbar}[q,\{p,\rho_{\xi}\ } ] \label{sled}\end{aligned}\ ] ] where @xmath36 . ( [ sled ] ) has been derived for potential models with canonical variables @xmath3 and @xmath37 with @xmath38 = i\hbar$ ] . in a wider context , it is still valid for moderate dissipative strength @xcite with the substitutions @xmath39 and @xmath40 $ ] .    both stochastic formulations for open quantum dynamics eq . ( [ sln ] ) and eq . ( [ sled ] ) have been successfully used to solve problems in a variety of areas . they apply to systems with discrete hilbert space as well as continuous degrees of freedom and have the particular benefit that they allow for a natural inclusion of external time dependent fields irrespective of amplitude and frequency . specific applications comprise spin - boson dynamics @xcite , optimal control of open systems @xcite , semiclassical dynamics @xcite , molecular energy transfer @xcite , generation of entanglement @xcite , and heat and work fluctuations @xcite , to name but a few . however , there is a price to be paid for the generality and simplicity of eqs . ( [ sln ] ) , ( [ sled ] ) . the correlation function eq . ( [ corrnunu ] ) requires the complex - valued process @xmath23 to have a random phase ; hence eq . ( [ sln ] ) describes non - unitary propagation . as in the paradigmatic case of multiplicative noise , geometric brownian motion @xcite , the stochastic variance of observables [ taken with respect to the probability measure of @xmath22 and @xmath23 ] grows rapidly with increasing time @xmath42 , making the computational approach prohibitively expensive in the limit of very long times . this problem has recently been solved by one of us @xcite by formally identifying the stochastic averaging as a projection operation . this allows the identification of relevant and irrelevant projections of the ensemble of state samples , with beneficial simplifications for finite memory times of the reservoir . the finite - memory stochastic propagation ( fmsp ) approach has solved the problem of deteriorating long - term statistics . however , the finite asymptotic value of the sampling variance can still be high in the case of strong coupling . here , we address this case using a different projection operator , based on finite decoherence timescales instead of reservoir correlation times . we restrict ourselves to reservoirs with ohmic - type spectral densities for which the sled in eq . ( [ sled ] ) is the appropriate starting point ; generalizations will be presented elsewhere .    in this section we first describe the general strategy and will then turn to specific applications in the remainder . since the latter focus on systems with discrete hilbert space , for convenience we make use of the language developed within the path integral description of discrete open quantum systems . there , it is customary to label path segments where the path labels @xmath43 and @xmath44 differ as `` blips '' , contributing to off - diagonal matrix elements ( coherences ) of @xmath7 and intervening periods with equal labels as `` sojourns '' , contributing to diagonal elements .    in an open quantum system , dephasing by the environment sets an effective upper limit on the duration of a blip . this observation can be transferred to our stochastic propagation methods using projection operators . the operator @xmath45 projects on sojourn - type intermediate states , while its complement @xmath46 projects on blip - type states .    with these projectors , eq . ( [ sled ] ) can be rewritten as a set of coupled equations of motion for diagonal and off - diagonal elements of @xmath47 , @xmath48\mathcal q \end{pmatrix } \begin{pmatrix }      \rho_{\mathcal p } \\ \rho_{\mathcal q } \end{pmatrix } \label{sys3}\ ] ] with a deterministic liouvillian superoperator @xmath49+\frac{\gamma}{2 i \hbar }         [ q,\{p,\cdot\}]$ ] and a stochastic superoperator @xmath50\xi(t)$ ] . the propagation of the sled equation ( [ sled ] ) can now be re - cast in the form of a nakajima - zwanzig equation @xcite @xmath51 where @xmath52 now is the projected ( diagonal ) part of the density matrix , @xmath53 is equal to the r.h.s . of eq . ( [ sled ] ) , and the symbol `` @xmath54 '' denotes time ordering . the time - ordered exponential now represents propagation during `` blip '' periods . generally , the lower integration boundary @xmath55 needs to be set to zero if full equivalence to equations of type ( [ sys3 ] ) is required . however , when it is known that dephasing sets an upper limit to blip times , @xmath55 can be raised , provided that @xmath56 remains large compared to the dephasing time .    in order to translate this approach into an algorithm , it is advantageous not to choose @xmath56 constant : the number of different initial conditions for the irrelevant part must be kept manageable . [ fig1 ] provides an illustration of the resulting algorithm , time - correlated blip dynamics ( tcbd ) . the propagation of coherences is realized in multiple overlapping segments on the time axis . an arbitrary , but fixed number of segments @xmath57 , i.e. @xmath58 with @xmath59 is used throughout the propagation . these segments are initialized with staggered starting times . with an inter - segment spacing @xmath60 chosen larger than the numerical timestep @xmath61 . the individual segments comprise a maximum memory window of @xmath62 . for each propagation step at a given time @xmath42 , the segment with the longest history among all @xmath58 is used in the propagation of @xmath63 . whenever one segment has `` aged '' beyond @xmath62 , it is reset and replaced by its next best follower , thus starting a continuous recycling of memory trails .    formally limiting the maximum blip length leads to improved statistics of our stochastic simulations , since the stochastic part of @xmath64 applies only to blip periods . the choice of a pre - defined number of memory segments ( independent of the propagation timestep ) reduces the required numerical operations significantly . compared to a naive solution of eq . ( [ eq : zwanzig ] ) with fixed difference @xmath56 , resulting in an unfavorable scaling of complexity with @xmath61 as @xmath65 ^ 2\right)$ ] , we lower the complexity to @xmath66 , typically an order of magnitude lower in absolute terms .    in the remainder of this work , we apply the new non - markovian propagation method tcbd to both a two - level system ( tls ) immersed in a heat bath ( spin - boson model ) and a quantum system that is effectively restricted by a three - dimensional hilbert space . in the transparent context of these discrete systems , we compare equilibrium properties to an analytical theory of dissipative two - state evolution , the non - interacting blip approximation ( niba ) and uncover super - exchange phenomena due to virtual particle transfer . the spin - boson model is a generic two - state system coupled linearly to a dissipative environment @xmath67 the coupling to the environment is conventionally taken as @xmath68 . the ohmic environment is described by a spectral density of the generic form @xmath69 where the constant @xmath70 denotes a coupling constant , and @xmath31 is a uv cutoff . in the context of the spin - boson problem one conventionally works with the so - called kondo parameter @xmath71 as a dimensionless coupling parameter . the sled and the tcbd are then applicable for moderate damping @xmath72 . subsequently , one sets @xmath73 , @xmath74 and @xmath75 so that the deterministic and stochastic superoperators introduced in eq . ( [ sys3 ] ) read @xmath76-\frac{\gamma}{2 i    \hbar } \delta[\sigma_z,\{\sigma_y,\cdot\}]\nonumber \\ \pazocal l_{\xi } & = & \frac{i}{\hbar}[\sigma_z,\cdot]\xi(t).\end{aligned}\ ] ] parameterizing the reduced density matrix through pseudospin expectation values leads to an intuitive picture of the dynamics @xmath77 we will later make use of the fact that the diagonal part is determined by the single parameter @xmath78 . it is instructive to compare our projection technique with the analytic niba approach @xcite . path integral techniques were first used to derive the niba ; with discrete path variables @xmath79 , eq . ( [ fvdoubint ] ) reads @xmath80\pazocal a^{\ast}[\sigma'(\tau ' ) ] \pazocal f [ \sigma(\tau)\sigma'(\tau ' ) ] . \label{nidoubint}\ ] ] since the path functions @xmath81 are piecewise constant , their derivatives are sums of delta functions at isolated points ( [ eq : lnf ] ) can then be integrated by parts ; this results in a double sum over `` interactions '' between discrete `` charges '' ( jumps in path variables ) , which depend on the twice - integrated reservoir correlation function @xmath82 \nonumber\\ & & + i \sin(\omega t)\}. \label{qfun}\end{aligned}\ ] ] the niba approach assumes that the resulting form of eq . ( [ eq : lnf ] ) can be approximated by omitting interactions between `` charges '' which are not part of the same blip . the laplace transform of the path sum can then be given in analytic form . for the present purpose , it is more instructive to note that this result can also be cast in the form of an equation of motion @xcite @xmath83 \label{niba}\ ] ] with integral kernels that depend on the relative time @xmath84 @xmath85 the function @xmath86 is given by @xmath87 for large @xmath31 and moderate damping ( @xmath88 ) . the integro - differential equation ( [ niba ] ) has a mathematical structure similar to the nakajima - zwanzig equation ( [ eq : zwanzig ] ) . this is not a coincidence . the dynamical variable @xmath78 of eq . ( [ niba ] ) is representative of the entire diagonal part of the density matrix , after averaging over reservoir degrees of freedom ( or , in our case , noise representing the reservoir ) . niba thus is related to the application of a different projection operator @xmath89 which sets off - diagonal elements to zero _ and _ takes the expectation value of the diagonal elements ( the relation @xmath90 is obvious ) .    however , the integral kernels of niba are equivalent to propagating with @xmath91 rather than @xmath92 . therefore eq . ( [ niba ] ) is not quite a nakajima - zwanzig equation . the difference between niba and tcbd could thus succinctly be stated in the following manner : niba tacitly omits the projection @xmath93 from the dynamics . apart from providing reference data and being a useful conceptual reference point , the niba theory provides us with a quantitative model for the system s decoherence time based on the reservoir s dissipative properties . when the function @xmath94 increases with @xmath95 at long times , long blips are suppressed as the dissipative factor @xmath96 in eq . ( [ difac ] ) becomes smaller . an estimate for a memory window @xmath62 of the tcbd s @xmath97 segments can be obtained for @xmath98 and sufficiently strong suppression of intra - blip interactions @xmath99 . considerably small errors are achieved by memory lengths of the order of @xmath100 longer memory time frames @xmath62 reduce potential errors due to a truncation of the off - diagonal trajectories . at the same time , they increase the amount of random noise that accumulates along the stochastic propagation of @xmath97 and smears out the original transfer signal . while the tcbd method allows to manually access intra - blip correlation lengths through @xmath62 , it puts no restrictions on inter - blip interactions , thereby extending the niba theory .    the efficiency gain of a numerical spin - boson simulation based on the sled eq . ( [ sled ] ) , and the proposed tcbd method , proves to be especially striking in the case of strong dephasing , i.e. for large coupling parameters @xmath101 . beside the dynamical observables @xmath102 and @xmath103 , fig . [ fig2 ] and fig . [ fig3 ] compare the variance of the stochastic sample trajectories @xmath104_t$ ] of a straightforward sled solution to the tcbd propagation of the reduced system . apparently , both first and second order statistics confirm substantial faster convergence of the tcbd scheme to the thermal equilibrium @xmath105 as well as a significant reduction in sample noise . + despite its sound quantitative applicability in a wide range of fields , the niba flaws in describing long time dynamics of @xmath106 ( @xmath107 ) correctly at low temperatures @xmath108 and finite energy bias @xmath109 . considering the equilibrium values for the @xmath103 component @xcite , the niba result comes down to a loss in symmetry and strict localization for zero temperature @xmath110 this confinement to one of the wells stands in contrast to the weak - damping equilibrium with respect to thermally occupied eigenstates of the tls @xmath111 with effective tunneling matrix element @xmath112^{1/2(1-k)}\left(\frac{\delta}{\omega_c}\right)^{k/(1-k)}\delta . \label{tuneff}\ ] ] while the equilibrium prediction of the niba theory clearly fails in the presence of an even infinitesimal bias energy , the tcbd provides correct equilibration in the scaling limit . note that due to the finite frequency cutoff in eq . ( [ osd ] ) , the tcbd approaches a steady state value determined by the effective tunnel splitting @xmath113 for coupling parameters @xmath114 . + results in fig . [ fig4 ] illustrate the relaxation process of @xmath103 for a finite bias and low temperatures .    while both the niba and the tcbd provide nearly indistinguishable data up to times @xmath115 , significant deviations appear for longer times . we note in passing the numerical stability of the tcbd which allows to access also typical equilibration time scales . while in the tcbd approach all time - nonlocal correlations induced by the reservoir are consistently taken into account as a vital ingredient for the reduced dynamics , the niba neglects long - ranged interactions such as inter - blip correlations . in this sense , the tcbd method constitutes a systematic extension of the niba . we will now demonstrate the adaptability of the tcbd method for multilevel systems by investigating population transfer in a three - state structure . for this purpose , we consider ( cf .  fig . [ fig5a ] ) a symmetric donor - bridge - acceptor ( dba ) system @xcite , with degenerate donor state @xmath116 , acceptor state @xmath117 , and a bridge state @xmath118 being energetically lifted by a bias of height @xmath109 . despite its simplicity the corresponding open quantum dynamics even in presence of an ohmic environment ( [ osd ] ) is rather complex and has been explored in depth as a model to access fundamental processes such as coherent / incoherent dynamics and thermally activated sequential hopping versus long - range quantum tunneling ( super - exchange ) .    within the spin-1 basis @xmath119 the hamiltonian of the system with site basis eigenvectors @xmath116 , @xmath118 and @xmath117 reads @xmath120 in eq . ( [ sled ] ) , the position coordinate @xmath3 is then represented by the @xmath121 , while the momentum operator @xmath37 follows from the respective heisenberg equation of motion @xmath122 = \delta s_y$ ] . accordingly , @xmath123 so that the deterministic and stochastic superoperators in eq . ( [ sys3 ] ) are derived as @xmath76+\frac{\gamma}{2 i \hbar } \delta[s_z,\{s_y,\cdot\}]\nonumber \\ \pazocal l_{\xi } & = & \frac{i}{\hbar } [ s_z,\cdot]\xi(t).\end{aligned}\ ] ]    as a first result , we show in fig . [ fig5 ] - [ fig7 ] the population dynamics of the site occupations @xmath124 and @xmath125 with @xmath126 for three different bridge energies @xmath127 , @xmath128 and @xmath129 ( in units of @xmath130 ) . for moderate bridge energies one observes coherent transfer of populations towards thermal equilibrium for the chosen coupling parameter @xmath131 and inverse temperature @xmath132 , while incoherent ( monotonous ) decay appears for high - lying bridges . the tcbd captures these qualitatively different dynamical regimes accurately and in domains of parameters space which are notoriously challenging , namely , stronger coupling and very low temperatures . for the remainder we will demonstrate how the tcbd approach can be used to extract transfer rates in case the population dynamics appears to be incoherent as in fig . these rates are of particular relevance for charge and energy transfer in molecular complexes or quantum dot structures , where the three - state system is the simplest model to exhibit sequential hopping from site to site as well as non - local tunneling between donor @xmath133 and acceptor @xmath134 . according to fig . [ fig5a ] , one has two different transfer channels , a sequential channel with transfer rate @xmath135 and a super - exchange channel @xmath136 . the respective dominance of these transfer mechanisms allows to classify the reduced system evolution as predominantly classical ( hopping ) or quantum ( tunneling ) . in the classical limit of high temperatures @xmath137 and low bridge state energies @xmath109 , the former channel is expected to prevail , while for lower temperatures and higher energy barriers , quantum non - locality in both the system and the reservoir degrees of freedom become increasingly important . assuming that these two rates define the only relevant time scales leads to a simple population dynamics of the form @xmath138 the eigenvalues of the rate matrix are obtained as @xmath139 , @xmath140 and @xmath141 with eigenvectors @xmath142 so that the general solution to ( [ dba ] ) is given by @xmath143 note that due to symmetries and detailed balance one has @xmath144 with boltzmann distributed equilibrium occupation probabilities @xmath145 and @xmath146 . asymptotically , the dynamics ( [ dba ] ) tends towards @xmath147 , @xmath148 . the expression ( [ eigsol ] ) constitutes the basis for a numerical extraction of the transfer rates @xmath149 and @xmath150 from simulation data . for this purpose , one introduces auxiliary functions @xmath151 and @xmath152 to cast eq . ( [ eigsol ] ) in the compact form @xmath153 such that the dynamical behaviour of @xmath151 and @xmath152 expressed by the population differences @xmath154 is obtained as @xmath155 \\ b(t)&=&\frac{1}{2}[\tilde{p}_3(t)-\tilde{p}_1(t)]\ , .\end{aligned}\ ] ] the numerical transfer rates @xmath149 and @xmath150 are now extracted by linear least square fits to the logarithm of the auxiliary functions @xmath151 and @xmath152 . choosing two parameter fit functions @xmath156 and @xmath157 , the transition rates can be computed from the fitting parameters as @xmath158      to compare the benchmark results obtained from the tcbd approach with approximate predictions , we come back to the niba discussed already in the previous section . as one of the most powerful perturbative treatments , the niba has also been the basis to derive analytic expressions for both sequential as well as super - exchange rates in the the domains @xmath159 and @xmath160 . this way , one arrives at the sequential forward rate @xmath161 \label{db}\ ] ] which can also be obtained from a fermi s golden rule calculation . upon expanding the dissipative function @xmath162 in eq . ( [ qfun ] ) to lowest order in @xmath163 and @xmath164 , a simple analytical expression is gained , i.e. , @xmath165 valid for coupling strength @xmath114 and with the effective tunneling matrix element eq . ( [ tuneff ] ) . going beyond the second order perturbative treatment to include also fourth order terms in @xmath130 , then leads to an approximate expression for the super - exchange rate @xmath166 \ , . \label{sqm}\ ] ]    while within the classical , high temperature domain @xmath167 thermally activated processes should be dominant ( [ db ] ) and manifest themselves in the typical arrhenius behavior @xmath168 , an increase in the bridge energy @xmath109 and lower temperature @xmath169 should lead to decay rates according to ( [ sqm ] ) . this is indeed seen in fig . [ dbsqmcr ] , where we depict the rates eq . ( [ db ] ) and ( [ sqm ] ) for varying bridge energies . the regime of sequential transfer crosses over to a regime where quantum tunneling dominates for sufficiently large @xmath109 , a behaviour that can now be compared to exact numerical results from the tcbd based on the two rate model of eq . ( [ dba ] ) .    a comparison between the theoretically predicted classical - quantum crossover in fig . [ dbsqmcr ] and the numerical results in fig . [ fig10 ] reveals at least good qualitative agreement while quantitatively the niba rates are not reliable .    in fig . [ fig11 ] we compare more specifically numerical results for @xmath170 with the niba golden rule predictions . the numerical data are in good agreement with an expected arrhenius behavior , i.e.  @xmath171 , but exceed the niba rates quite substantially .    for increasing barrier heights , we expect the super - exchange rate to control the population decay if @xmath172 . the perturbative treatment leads to a characteristic algebraic dependence @xmath173 with @xmath174 , in contrast to the exponential one for sequential hopping . the results in fig . [ fig12 ] for moderate to large bridge energies verify the expected power law dependence with an exponent @xmath175 , close to the perturbative expectation .      in case of weak system - reservoir couplings , the master equations build standard perturbative approaches to open system dynamics @xcite . the corresponding time evolution equations for the reduced density matrix are formulated in terms of the eigenstates of the bare systems with the dissipator inducing transitions between these states . while these methods clearly fail for stronger couplings , it is nevertheless instructive to analyze their deficiencies in the context of rate dynamics .    here , we use a standard master equation ( me ) for which the dynamics of the diagonal elements of the reduced density matrix decouples from that of the off - diagonal elements . in the eigenstate representation of the hamiltonian eq . ( [ 3ham ] ) , @xmath176 , one then has @xmath177 \label{pma}\ ] ] where @xmath178 . the transition rates are obtained as @xmath179 with @xmath180 and the thermal occupation @xmath181 $ ] . the stationary solution of ( [ pma ] ) reproduces the gibbs state distribution @xmath182 with partition function @xmath183 . the dynamics of the off - diagonal elements can simply be solved @xmath184 with decay rates @xmath185 \nonumber\\ & & - \frac{1}{\hbar^2 } { \ensuremath{\left\langlem\right| } } s_z { \ensuremath{\left|m\right\rangle } } { \ensuremath{\left\langlen\right| } } s_z { \ensuremath{\left|n\right\rangle } } d(0 ) . \label{oddec}\end{aligned}\ ] ] now , the population dynamics @xmath186 in the site representation is obtained from this time evolution by a simple unitary transformation yielding @xmath187 note that the site populations are also determined by the dynamics of the off - diagonal elements in the eigenstate representation ; the respective timescales are entirely determined through eqs . ( [ trama ] ) , ( [ oddec ] ) . we start in fig . [ fig15 ] with the population dynamics in the coherent regime ( weak coupling ) . while the master results capture the transient oscillatory pattern qualitatively correctly , quantitatively deviations are clearly apparent . the time scales for relaxation towards thermal equilibrium are quite different with the master results approaching a steady state on a much faster time scale than the tcbd data . this may be attributed to the relatively low temperature which is beyond the validity of the born - markov approximation on which the master equation is based . the steady state values for the populations are basically identical though . discrepancies in the relaxation dynamics substantially increase for stronger system bath couplings , see fig . [ fig16 ] , when the dynamics tends to become a simple decay in time . based on the rate extraction method presented before , we can now also extract respective rates from the master equation dynamics providing us with corresponding rates @xmath188 and super - exchange tunneling @xmath189 . the sequential transfer rate in fig . [ fig18 ] reveals indeed an exponential decay @xmath190 , however , with substantial deviations @xmath191 from the thermal value @xmath192 . the situation for the super - exchange rates is even worse as the extracted rates lack a physical interpretation , see fig . [ fig19 ] : the numerical @xmath189 saturate for growing bridge energies to eventually become independent of @xmath109 at all , in contradiction to an algebraic decay . this result is not so astonishing since , as already mentioned above , the parameter domain where incoherent population dynamics exhibits substantial quantum effects lies at the very edge or even beyond the range of validity of master equations . what is interesting nevertheless , is the fact that the explicit time scales appearing in eq . ( [ sitra ] ) can be directly related to the extracted super - exchange rate @xmath189 . this is shown in fig . [ fig20 ] . it turns out that indeed non - local processes in the site representation ( super - exchange ) correspond to the decay rate of off - diagonal elements of the density matrix in the energy representation implying that already for moderate bridge energies @xmath193 . in this work we combined a stochastic description of open quantum dynamics with projection operator techniques to improve convergence properties . while the general strategy has been outlined recently by one of us @xcite , we here restricted ourselves to a broad class of systems for which the dephasing time is finite . for the dynamics of the coherences of the reduced density one then keeps track of bath induced retardation effects only within a time window @xmath62 which must be tuned until convergence is achieved . while for @xmath194 with @xmath195 being the full propagation time one recovers the full stochastic description , the new tcbd is superior if @xmath62 is sufficiently shorter than @xmath195 . in parameter space this applies to the non - perturbative regime beyond weak coupling , a domain which is notoriously difficult to tackle , particularly at low temperatures . the new tcbd scheme captures systems with continuous degree of freedom for arbitrary coupling strengths and to those with discrete hilbert space up to moderate couplings . it allows to approach also the long time domain , where equilibration sets in . for two- and three - level systems we have shown explicitly that the new scheme covers both coherent as well as incoherent dynamics . in this latter regime the method is sufficiently accurate to extract quantum relaxation rates in the long time limit . it now allows for a variety of applications , for example , the dynamics of nonlinear quantum oscillators also in presence of external time dependent driving , the heat transfer through spin chains , or the quantum dynamics under optimal control protocols .
Abstract:

Answer: the non - markovian dynamics of open quantum systems is still a challenging task , particularly in the non - perturbative regime at low temperatures . while the stochastic liouville - von neumann equation ( sln ) provides a formally exact tool to tackle this problem for both discrete and continuous degrees of freedom , its performance deteriorates for long times due to an inherently non - unitary propagator . here we present a scheme which combines the sln with projector operator techniques based on finite dephasing times , gaining substantial improvements in terms of memory storage and statistics . the approach allows for systematic convergence and is applicable in regions of parameter space where perturbative methods fail , up to the long time domain . findings are applied to the coherent and incoherent quantum dynamics of two- and three - level systems . in the long time domain sequential and super - exchange transfer rates are extracted and compared to perturbative predictions .


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: a practical single photon source should meet the demands of high fidelity , brightness and efficiency . as a source of photon pairs , parametric down - conversion ( pdc ) enables conditional state preparation  heralding on one of the twin photons ( idler ) gives exact timing information about the other ( signal ) . in addition to this , single pass pdc sources have the advantage of operational simplicity and the quality of the quantum state preparation is mainly limited by the available detectors . recent experiments have demonstrated a full tomographic characterization of fock states or photon added / subtracted states prepared from the pdc emission in the low gain regime  @xcite . nevertheless , since the photon - number distribution of signal and idler is thermal , higher photon - number components degrade the fidelity of the prepared target state . this effect becomes prominent at higher gains i.e. if high preparation rates are needed . in addition to photon - number correlation , spectral correlation between the pdc photon pairs can not be avoided in conventional sources  @xcite . although the spectral correlation is a useful resource for the generation of highly entangled states  @xcite , the preparation of non - gaussian states in one of the twin modes leads to mixedness whenever the other field mode is traced out  @xcite .    in the past , the manipulation of the spectral properties has been in the interest of broadband and efficient non - linear classical frequency conversion  @xcite . more recently , tailoring these properties has been suggested as a method of source engineering for quantum optics  @xcite . several sophisticated methods have been proposed for state decorrelation , which unfortunately often require means that are available only under specific conditions  @xcite . nevertheless , the decorrelation of the pdc photons is always possible via spectral filtering , which drastically changes the shape of the spectral correlation function . the influence of the spectral correlation on the purity of a single photon state prepared by filtering was first experimentally demonstrated by wasilewsky _ et al . _ @xcite . in addition to the optimal state preparation , a precise characterization of the prepared photonic states both in photon - number and spectral degrees of freedom is required . the latest developments in photon counting technologies allows to characterize photonic states via the direct detection of the photon - number statistics . one of these techniques applicable for pulsed light is time - multiplexed detection ( tmd ) , which enables a loss - tolerant reconstruction of photon statistics  @xcite . regarding the spectral properties , their impact on the state preparation can be studied via hong - ou - mandel ( hom ) type interference , originally applied to observe the temporal length of pdc photons  @xcite . this experiment can be implemented with a symmetric beam splitter and a coincidence counting device . two indistinguishable photons impinging simultaneously on the different ports of the beam splitter bunch together and no coincidences are detected . thus , the measurement of hom interference dip between signal and idler is a measure of their spectral indistinguishability , or in other words of their equity  @xcite . however , for characterizing the fidelity of the prepared single photon , one must employ an independent , spectrally pure reference field , which can be a highly attenuated coherent state . the measurement of hom interference between these two independent sources allows to quantify their spectral overlap , as evidenced in several studies  @xcite . considering the direct probing of more evolved non - gaussian states , such as displaced single photons , the determination of this overlap is essential @xcite .    lately , pdc sources with extreme high brightness have been accomplished by employing waveguided structures rather than bulk crystals  @xcite . as filtering happens at the cost of photon flux , waveguides ( wgs ) can offer benefits for preparing high fidelity single photons . in this paper we show how to optimally adapt the bandwidth of the filter such that the trade - off between the fidelity and brightness of the source is balanced . in sec . [ sec_ellipse ] we provide easy applicable tools for analyzing and estimating the spectral correlation of a waveguided pdc source and give a recipe of how to design a filter for optimal state decorrelation . the spectral properties and the performance of the employed filters are verified in sec . [ sec_decorrelation ] . our main focus , the characteristics of the filtered state is presented in sections [ sec_statistics ] and [ sec_hom ] ; we study the statistics of the heralded state and finally determine the spectral overlap between signal and reference . this allows us , to our knowledge for the first time , to directly extract the overall fidelity of a heralded single photon by utilizing the hom interference . the quantum state of photon pairs generated via pdc is usually described as @xmath0 where @xmath1 is the pair creation probability proportional to the pump power , and the labels @xmath2 refer to signal and idler . the spectral correlation function @xmath3 arises as a consequence of energy and momentum conservation . the former is constrained by the spectrum of the pump pulse @xmath4 and the latter , being governed by the dispersion of the non - linear optical medium of the length @xmath5 , is characterized by the phase - matching ( pm ) function @xmath6 . we study the spectral correlation in the vicinity of the central frequencies of pump , signal and idler @xmath7 ( @xmath8 ) , at which @xmath9 and @xmath10 . the grating period @xmath11 of the quasi - phasematched structure can be used for adapting the desired central wavelengths . ( see e.g.  @xcite ) .    for simplicity , we describe the spectrum of the pump by a gaussian envelope @xmath12 where @xmath13 is the detuning from the central wavelength . in addition , we expand the phase - mismatch @xmath14 appearing in the pm function in first order taylor series around @xmath7 and use the gaussian approximation for the sinc - function . thus , the pm function can be estimated by @xmath15 where @xmath16 is determined by the group - velocity mismatch of the non - linear medium and the parameter @xmath17 adapts the width of the gaussian to the width of the original sinc - function . finally , the amplitude of the spectral correlation function is written as @xmath18    \underbrace{\left[\begin{array}{cc }    \frac{1}{\sigma^{2}}+ \gamma \frac{l^{2}}{4}\kappa_{s}^{2 } &   \frac{1}{\sigma^{2}}+ \gamma \frac{l^{2}}{4}\kappa_{s}\kappa_{i } \\   \frac{1}{\sigma^{2}}+ \gamma \frac{l^{2}}{4}\kappa_{s}\kappa_{i } &    \frac{1}{\sigma^{2}}+ \gamma \frac{l^{2}}{4}\kappa_{i}^{2 } \end{array } \right]}_{\mathcal{m } }    \left[\begin{array}{c } \nu_{s}\\ \nu_{i }   \end{array } \right ] } \bigg ) , \label{eq_ellipse}\end{aligned}\ ] ] which represents an ellipse in the @xmath19space as illustrated in fig .  [ fig_correlation](a ) . for a decorrelated and separable photonic state the ellipse is oriented parallel to the @xmath19coordinates , i.e. the matrix @xmath20 in eq . ( [ eq_ellipse ] ) is diagonal . non - vanishing off - diagonal terms correspond to a correlated photonic state . we stress that these properties can be tailored by modifying the pump bandwidth or the pm properties , e.g. by adjusting the waveguide length . the width of the ellipse axes are determined by the eigenvalues of @xmath20 and are equal to @xmath21 whereas the tilt can be extracted from the eigenvectors . in the limit of spectrally broadband pump ( @xmath22 ) the tilt @xmath23 is given by @xmath24 and the minor axis of the ellipse is purely determined by the pm properties , @xmath25 in fig . [ fig_correlation](b ) we show the behavior of this width as a function of @xmath5 considering the expected material properties of our wg . even though this theoretical prediction gives an indication about the pm properties , the actual spectral correlation function needs to be accurately characterized . this is especially important for waveguided pdc as the non - ideal properties of the guide typically influence the pm function significantly .    in terms of fock state preparation , the joint amplitude of signal and idler in eq . ( [ eq_ellipse ] ) has to be decoupled . to optimally achieve this goal we employ a short wg , thus broadening the pm bandwidth , and a femtosecond pump , which allows the state decorrelation via a proper design of the spectral filter . filtering reduces the width of the ellipse s major axis and modifies the correlation such that the ellipse evolves into a circular shape or has a very small remaining tilt . thence , the bandwidth of the filter has to be matched with the width of the ellipse s minor axis . one of the most direct ways to characterize the tilt of the spectral correlation ellipse is to measure the unconditional spectral marginal distributions of signal and idler . the tilt is determined by the simple relation @xmath26 , where @xmath27 is the fwhm of the marginal . however , several spatial modes can propagate in the wg and an accurate determination of the tilt via this method becomes difficult  @xcite . another possibility relies on detuning the pump wavelength and observing the respective change in the central wavelengths of signal and idler . similarly , we can also employ narrowband filters in signal and idler arms , and record the wavelengths at which photon - number correlation appears , thus mapping the spectral correlation function .    in all pdc measurements presented in this paper we employed a long , type - ii , periodically poled ktp wg , which was pumped by a frequency doubled ti : sapphire laser operating at ultrafast regime ( , fwhm , autocorrelation length and repetition rate ) . for the estimation of the ellipse s tilt we separated signal and idler photons in a polarizing beam splitter and launched them into single mode ( sm ) fibers . the single photon spectra was measured with a sensitive spectrograph ( andor ) and the photon - number correlation with avalanche photo diodes ( apds ) as illustrated in fig . [ fig_marginals](a ) . first of all , we measured the spectral marginal distributions of signal and idler as shown in fig . [ fig_marginals](b ) . after that , two narrowband interference filters were set on the beam paths of signal and idler and the photon - number correlation was mapped by tuning the filters . for the accurate determination of the ellipse s tilt the filter wavelengths were recorded with the spectrograph . from the results shown in fig . [ fig_marginals](c ) we extracted the following value for the tilt : @xmath28 .    regarding the major axis of the correlation ellipse , note that its width is highly dependent on the bandwidth of the pump as illustrated in fig . [ fig_correlation](a ) . the fwhm of the major axis can be roughly estimated from the marginal distributions by @xmath29 . our results indicate that pumping the wg with broad pulse leads to an approximately broad major axis , which increases to values between , if the pump is broad . according to fig . [ fig_correlation](b ) the exact determination of the pm bandwidth of our wg , i.e. the minor axis of the ellipse , requires sub - nanometer resolution . instead of employing sensitive monochromators , we studied the second harmonic ( sh ) response of the wg . in this process two photons at the fundamental frequency give rise to a single up - converted one . therefore , the fundamental probe always lies on the line @xmath30 of the @xmath31space , as illustrated in fig . [ fig_shg_response](a ) . in other words , tuning the frequency of the fundamental probe and measuring the respective sh response allows us to characterize the width of the pm function . if the fundamental probe is narrow with respect to the pm bandwidth @xmath32 , and the tilt of the pm only slightly deviates from 45@xmath33 , the envelope of the sh response is approximately equal to @xmath34    we recorded the sh response with a grating spectrometer by tuning the fundamental probe ( , fwhm , and autocorrelation length ) in steps of approximately over the phase matching curve [ setup shown in fig . [ fig_shg_response](b ) ] . in our case , the measurement is limited by both the resolution of the spectrograph ( at ) and the bandwidth of the fundamental probe beam . nevertheless , our results shown in fig . [ fig_shg_response](c ) indicate that we can probe the shape of the pm function and that its bandwidth is on the order of the probe bandwidth .    by combining our findings we are able to reconstruct the correlation ellipse @xmath35 of our photon source . using the estimated values of the tilt and width of the pm , and taking into account a pump with fwhm of we obtain the ellipse shown in fig . [ fig_ellipse_rec](a ) . the result implies that signal and idler are highly frequency correlated , which is a consequence of the long interaction length of the pdc photons in the wg . further , our result suggests that the optimal filter bandwidth is on the order of . unfortunately , spectral filters at with high transmission and narrow bandwidth are rarely available . in consequence , we study in the next section the decorrelation of the generated photonic state by employing spectral filters with and bandwidths . we expect slight remaining correlations after sending the prepared state through these filters as shown in figs  [ fig_ellipse_rec](b ) and [ fig_ellipse_rec](c ) , respectively . the indistinguishability between signal and idler can be studied via a hom interference , where signal and idler are superimposed in a 50/50 beam splitter . the full visibility is achieved either if the correlation function is circular , or if it is symmetric with respect to the axis of energy preservation , i.e. if the correlation ellipse is oriented at @xmath36  @xcite . here we employ this measurement to test and corroborate the acquired knowledge about the spectral ellipse . we adapt the analysis of avenhaus et al . @xcite to the study of the hom visibility with spectral filter . the interference between signal and idler can be characterized by the overlap factor @xmath37 using the gaussian approximation of the spectral correlation function , the overlap factor can be expressed as a product of a spectral ( @xmath38 ) and temporal ( @xmath39 ) overlap i.e. , @xmath40 . while spectral mismatch results from an asymmetric spectral correlation function , temporal mismatch originates from the phase terms , which are ultimately related to the properties of the birefringent media . these can be inferred from the properties of the correlation ellipse via eq . ( [ eq_pm ] ) . if perfect overlap was necessary , the temporal mismatch could be compensated with an optical delay line . the defined overlap factors depend on the ellipse s tilt @xmath23 and aspect ratio @xmath41 , defined as the ratio between the major and minor axes . they can be expressed as @xmath42 @xmath43\,\mathcal{a}^{2 } + [ 1-\sin({2\theta})]^{2 } } { [ 1+\mathcal{a}^{4}][1-\sin^{2}({2\theta})]+ 2\,\mathcal{a}^{2}[1+\sin^{2}({2\theta } ) ] } } \right)}\ ; , \label{eq_o_temporal}\ ] ] where we take use of an estimation that the pm width and the ellipse s minor axis are equal . since the twin - photons of a type - ii process have orthogonal polarizations , we measured hom interference fringes between signal and idler by rotating a half - wave plate in front of a polarizing beam splitter and recording the coincidence counts between the two output arms . the experiment was realized in three different cases : ( i ) without spectral filter , ( ii ) including in the setup a spectral filter with fwhm of and , ( iii ) decreasing the width of the filter down to . the filters were inserted in the joint beam path of signal and idler as shown in fig . [ fig_setups](a ) . according to auxiliary experimental observations , the time delay between p- and s - polarized modes is on the order of . since the temporal width of signal and idler , estimated by the filter bandwidth , is on the order of and no delay line is inserted in the setup , the temporal overlap can not be neglected in our measurements . the visibility is defined by @xmath44 , where @xmath45 indicates the amount of coincident counts . in terms of the overlap factor , the visibility can also be expressed as @xmath46 . the classical case , in which signal and idler are completely distinguishable , corresponds to the situation of zero overlap , i.e. @xmath47 . we measured , without any background subtraction , the following visibilities : @xmath48 , @xmath49 , and @xmath50 , corresponding to each of the above mentioned cases . the corresponding overlap factors must be compared with the theoretical prediction given in eqs  ( [ eq_o_spectral ] ) and ( [ eq_o_temporal ] ) . for this we approximate the aspect ratios in a straightforward manner . when no filtering is applied , @xmath41 is estimated as the ratio of the experimentally extracted major and minor axes , i.e , @xmath51 . when spectral filtering is applied , we determine @xmath41 as the ratio of the filter bandwidth to the pm width i.e. , @xmath52 and @xmath53 . as expected , the spectral overlap alone does not explain our results , as can be seen in fig . [ fig_aspect ] . by including the temporal overlap in the analysis , the employed model predicts the right tendency , even though in the experiment we do not have any compensation for dispersion or chirp . we conclude that our model , which includes only experimentally extracted characteristics of the pm , explains with high confidence results of an independent experiment . in addition , it predicts how the spectral correlation between signal and idler can be manipulated with the selected filter . in the next sections we study how the heralded state is affected by the higher photon - number contributions and remaining spectral correlation . in the low gain regime heralding of pure signal states is normally realized by considering a single spectral filter at the trigger . however , at higher gains multiple modes are excited simultaneously . thus , filtering only the trigger results in the detection of uncorrelated signal modes . in addition to this , the higher photon number contributions decrease the fidelity of the heralded single - photon state . we studied these effects by measuring , how the heralded statistics depends on whether or not a narrowband spectral filter is added in the signal arm . as before , we pumped our pdc source with a frequency doubled ti : sapphire laser . the pump beam was sent through a simple spatial mode cleaner before being launched into the wg . after that , the pump was filtered out and the twins were divided into two independent paths . the idler was sent through broad interference filter and , finally , coupled into sm fiber . we measured the heralded statistics in two different configurations : either without filter or with a 2.5  nm interference filter inserted in the signal arm . the setup shown in fig . [ fig_setups](b ) was arranged such that the idler was used as a trigger and the statistics of the signal was measured with a `` two - bin tmd '' . this detector consists of a single 50/50 beam splitter , constructed by using a half - wave plate in front of a polarizing beam splitter . the loss tolerance of the tmd was achieved by accurate calibration of loss @xmath54 and convolution matrices @xmath55 @xcite . the signal statistics @xmath56 was then deduced by inverting the click statistics @xmath57 with a maximum - likekihood technique @xcite .    from our results one can observe that the heralded statistics presents a large two - photon contribution [ fig . [ fig_statistics_filtered](a ) ] , which is reduced by placing a narrow band filter in the signal arm [ fig . [ fig_statistics_filtered](b ) ] . we understand this result as an indication that , at higher gains , a large amount of distinct spectral modes are excited simultaneously . if two idler photons are generated there is a high probability that just one of them has a good overlap with the filter mode . therefore , by using a filter in the signal arm it is possible to reject the signal photon created together with the spectrally mismatched idler . the click statistics were inverted using the calibrated efficiencies , which in both cases were @xmath58 , corresponding to the ratio of coincident counts to singles . the highest trigger rate was , which is equivalent to a heralding probability of 0.16% . regarding the single - photon state preparation , the one - photon contribution of the filtered signal was 95% at this trigger rate . a comparison between the filtered and unfiltered statistics allows us to extract the degree of mode reduction given by the filter . for that , we model the multimode statistics as the convolution of several uniformly occupied spectral modes  @xcite . we begin the analysis with the single - mode case . the heralded statistics , being conditioned on a trigger event , has the form @xmath59 moreover , the joint probability of having @xmath60 photons in signal and detecting a click from idler is given by @xmath61 $ ] , where @xmath62 is the efficiency at which the filtered mode is detected and @xmath63 the probability of creating @xmath60 signal ( idler ) photons , which always occur in pairs for pdc . since in the experimental realization losses are very high , the joint probability can be approximated by @xmath64 . the overall click probability is estimated as @xmath65 , where @xmath66 indicates the mean photon number of the marginal pdc statistics . in this loss regime , the heralded statistics become independent of the trigger efficiency @xmath62 and its mean photon number can be evaluated as @xmath67 , with @xmath68 being the second moment of the marginal pdc statistics .    in the case of a multimode pdc with several uniformly excited modes , the statistics @xmath69 is given by the convolution of the modes , each having the mean photon number @xmath70 and the second moment @xmath71 . for @xmath72 excited modes we find that @xmath73 and @xmath74 , from which we can calculate the multimode mean photon number of the heralded statistics . considering a thermal single mode marginal distribution , @xmath75 with @xmath76 being the gain of the pdc mode , the mean photon number of the heralded statistics can be approximated by @xmath77 we conclude that , in the limit of low pump powers , such that @xmath78 , the slope of the mean photon number depends on the number of modes . in fig .  [ fig_statistics_filtered](c ) we show the mean photon number of the measured heralded statistics with respect to the pump power . by comparing the slopes , we estimate that the reduction in the number of modes due to the filter at signal is on the order of 30 . we emphasize that the model used is an approximation , where all modes are treated similarly , even though the trigger mode is a special one . in this section we study the hom interference between two independent sources : the heralded signal and a strongly attenuated coherent field , here called reference . the visibility of the interference dip depends , on the one hand , on the presence of higher photon - number components and , on the other hand , on the spectral overlap between the target state and the reference . we experimentally investigate the cases of two- and three - fold hom interferences . the former case corresponds to the interference between signal and reference neglecting the pdc photon - number correlation , while in the latter case the signal state is heralded by trigger events . we begin our study by considering a spectrally pure signal state . further , we decompose the signal only up to the second order in the photon - number basis , which is justified under the experimental conditions . after that , we complete our description by adding the effect of spectral correlations between signal and idler . in the experiment , signal and reference are launched into the ports @xmath79 and @xmath80 of a 50/50 beam splitter . the density matrix of the input state @xmath81 can be written as @xmath82 in which @xmath83 ( @xmath84 ) corresponds to the statistics of the photon - number mixed signal mode with @xmath85 . the reference field has amplitude @xmath86 and spectral distribution @xmath87 , such that @xmath88 . since mode `` @xmath79 '' is spectrally pure , its one- and two - photon components can be respectively described by @xmath89 and @xmath90 , in which the spectral distribution @xmath91 is normalized , i.e. @xmath92 . interference does not take place when the two input pulses arrive at the beam splitter at different times . to model that , we explicitly write the evolution of the fields considering a time delay @xmath93 . in the frequency domain , the beam splitter relations are given by @xmath94e^{i\omega \tau } , \nonumber \\ \hat{b}^{\dagger}(\omega ) & = & \frac{1}{\sqrt{2}}\left [ \hat{c}^{\dagger}(\omega)+\hat{d}^{\dagger}(\omega)\right]\,,\end{aligned}\ ] ] in which @xmath95 and @xmath96 are the two output modes . these equations allow us to calculate the evolved state @xmath97 . finally , the probability of a coincident click can be calculated via @xmath98 , where @xmath99 describes the povm of two simultaneous clicks in the apds . using the bosonic commutation relation @xmath100= \delta(\omega-\tilde{\omega})$ ] and the fact that @xmath101 , the probability of detecting a coincident click is given by @xmath102 with @xmath103 . the functions @xmath104 and @xmath105 have the form @xmath106 where @xmath107 and @xmath108 . at origin ( @xmath109 ) these factors determine the spectral overlap of the one- and two - photon wave - packets with respect to the reference . finally , the visibility of the interference dip is defined by @xmath110    in order to simplify our model we compare the independent impact of the one- and two - photon components of signal on the visibility . considering perfect spectral overlap between these components and the reference field , and appliying values of @xmath111 compatible with the experimental realization , we estimate that the interference of the one - photon component would have a visibility of above 98% , whereas the one of the two - photon component would be less than 5% . therefore , in the regime of power considered here , the two - photon component of signal gives only a background contribution to the coincidences . thus , we consider the spectral degree of freedom exclusively for the one - photon component of signal and reference . further , we approximate @xmath112 , which is well valid within the experimental frame , and neglect the terms arriving from the interference of more than two photons ( total number in mode `` @xmath79 '' plus `` @xmath80 '' ) . the remaining contributions are shown in fig . [ fig_ball_model ] and they can be retrieved from eq . ( [ eq_p ] ) as follows . the first term corresponds to accidental counts due to the two - photon component of the reference and it can be re - expressed as @xmath113 [ fig . [ fig_ball_model](a ) ] . the second and third terms account for the events where both signal and reference contain one single photon [ fig . [ fig_ball_model](b ) ] . the probability of simultaneous accidental clicks due to one photon from signal and another from reference is given in the second term as @xmath114 . considering the first - order photon - number contribution , the third term describing the quantum interference between the one - photon components of signal and reference can be rewritten as @xmath115 . the fourth term takes into account accidental counts coming from signal , @xmath116 [ fig . [ fig_ball_model](c ) ] . the last term , corresponding to interference of at least three photons , is an order of magnitude smaller the other terms and thence it is neglected .    with the described simplifications , eq . ( [ eq_p ] ) can be rewritten as @xmath117 after substituting eq . ( [ eq_pt ] ) into eq . ( [ eq_visibility_def ] ) , the visibility can be written with respect to the mean photon number of the reference @xmath111 as @xmath118 which reaches its maximum at @xmath119 .    in the experimental situation the joint - spectral distribution of signal and idler is not separable , and thus the prepared state is spectrally mixed . our model can be adapted to this condition by replacing the function @xmath120 by a new one , in which a trace operation over the idler mode is realized . the spectral density of the pdc state after filtering is given by @xmath121 where the @xmath122 is the intensity transmission profile of the signal and idler filters employed in the state generation and @xmath123 accounts for the normalization . apart from the visibility , the width of the hom interference dip offers information about the spectral properties of the prepared state . using the gaussian approximation for the spectral correlation function , reference , and filter profiles , the width of the hom interference dip can be estimated via eq . ( [ eq_pt ] ) . disregarding the background contributions , it can be written as @xmath124 where the temporal width of the interference dip , @xmath125 , is determined by the inverse widths of the pump @xmath126 , reference @xmath127 , signal filter @xmath128 , and phase matching @xmath129 , as well as by the ellipse s tilt @xmath23 . the amplitude @xmath130 describes the maximal expected value of the overlap , ranging from zero to one . we emphasize that the two- and three - fold interference dips have the same temporal width . for the measurement of the hom interference we used the setup shown in fig . [ fig_setups](c ) . a portion of the original laser beam was heavily attenuated and launched through a polarization maintaining sm fiber - coupler in s - polarization . the signal was launched to the fiber in p - polarization and overlapped with the reference at a 50/50 beam splitter , consisting of a half - wave plate and a polarizing beam splitter . two 1  nm broad interference filters were inserted in the setup : one in the path of the trigger and the other in the joint arm of the signal and reference . we measured the hom interference dip in the two- and three - fold cases for different amplitudes of the reference . by blocking the reference arm we measured the statistics of the signal state and obtained the results summarized in table [ table_statistics ] . .measured signal state statistics @xmath83 @xmath131 without loss inversion . [ cols="^,^,^,^",options="header " , ]     in figs  [ fig_3-fold](a ) and [ fig_3-fold](b ) we show example measurements of the two- and three - fold interference dips that were recorded with reference powers close to the optimized values @xmath132 in fig . [ fig_overlap ] . the measured visibilities were @xmath133 and @xmath134 , respectively . the temporal width of the interference dips was @xmath135  ps . following eq . ( [ eq_tau ] ) we theoretically study how this width changes as a function of the pm bandwidth [ fig .  [ fig_3-fold](c ) ] . the comparison of the measurement with our model suggests a fwhm of for the pm , which is in good agreement with the value extracted from sh measurements . further exploring our knowledge about the spectral correlation function , we use eqs . ( [ eq_overlapfacators ] ) and ( [ eq_g_pdc ] ) in order to estimate the maximal spectral overlap @xmath130 in the two- and three - fold cases . in fig . [ fig_3-fold](c ) we illustrate the theoretical predictions and compare them with the measured overlap values . in the three - fold case we find a good agreement with the model . however , the measured value of the two - fold overlap deviates largely from the expected one . as no selection of the detection events by triggering is applied in the two - fold case , we believe that the properties of non - ideal wg , such as fluorescence and higher order spatial modes with complex spectral structure , degrade the result . in addition to this , the expected values have not been achieved due to other experimental imperfections , in particular mechanical instabilities and spatial mode mismatch . we have studied the ability to tailor the spectral properties of a conventional pp - ktp waveguide . information about the spectral correlation between signal and idler was retrieved from simple experiments , enabling to model our system with a standard theory of parametric down - conversion . the viability of this model was tested via the hong - ou - mandel interference between the twin beams . since the production of more evolved non - gaussian states , as displaced single - photon fock states , requires an overlap between the signal and an independent reference field , we studied the quantum interference between the heralded state and a coherent reference field . the spectral overlap was directly extracted from the visibility of the hong - ou - mandel interference . further , we experimentally characterized the statistics of the heralded state . considering both the spectral overlap and the one - photon contribution , we obtained by a direct detection the fidelity of 78% between the prepared state and a single - photon fock state . in addition , our observations show that heralding decreases the number of modes only if both signal and idler are spectrally filtered . although we found that in the current experiment the spectral overlap between the reference and the signal state is not yet optimized , we provide a clear and simple recipe for tailoring the characteristics of the spectral correlation and for testing the degree of decoupling , i.e. , the fidelity of the heralded state . we thank m.  avenhaus , a.  christ and p.  j.  mosley for valuable and helpful discussions . this work was supported by the ec under the grant agreement corner ( fp7-ict-213681 ) . k.n.c .  acknowledges the support from the alexander von humboldt foundation . h.  s. poh , c.  y. lum , i.  marcikic , a.  lamas - linares , and c.  kurtsiefer , `` joint spectrum mapping of polarization entanglement in spontaneous parametric down - conversion , '' phys . a * 75 * , 043816 ( 2007 ) .        m.  v. fedorov , m.  a. efremov , p.  a. volkov , e.  v. moreva , s.  s. straupe , and s.  p. kulik , `` anisotropically and high entanglement of biphoton states generated in spontaneous parametric down - conversion , '' phys . lett . * 99 * , 063901 ( 2007 ) . m.  avenhaus , m.  v. chekhova , l.  a. krivitsky , g.  leuchs , and ch . silberhorn , `` experimental verification of high spectral entanglement for pulsed waveguided spontaneous parametric down - conversion , '' phys .  rev . a * 79 * , 043836 ( 2009 ) . a.  b. uren , ch . silberhorn , r.  erdmann , k.  banaszek , w.  p. grice , i.  a. walmsley , and m.  g. raymer , `` genreration of pure single photon wavepackets by conditional preparation based on spontaneous parametric downconversion , '' laser  phys . * 15 * , 146161 ( 2005 ) . p.  j. mosley , j.  s. lundeen , b.  j. smith , p.  wasylczyk , a.  b. uren , ch . silberhorn , and i.  a. walmsley , `` heralded generation of ultrafast single photons in pure quantum states , '' phys . 100 * , 133601 ( 2008 ) . m.  halder , j.  fulconis , b.  cemlyn , a.  clark , c.  xiong , w.  j. wadsworth , and j.  g. rarity , `` nonclassical 2-photon interference with separate intrinsically narrowband fibre sources , '' opt .  express * 17 * , 46704676 ( 2009 ) . m.  avenhaus , h.  b. coldenstrodt - ronge , k.  laiho , w.  mauerer , i.  a. walmsley , and ch . silberhorn , `` photon number statistics of multimode parametric down - conversion , '' phys .  rev .  lett . * 101 * , 053601 ( 2008 ) . p.  kolenderski and k.  banaszek , `` testing single - photon wave packets by hong - ou - mandel interference , '' in _ proceedings of ieee 2008 6th national conference on telecommunication technologies and ieee 2008 2nd malaysia conference on photonics . nctt - mcp 2008 . _ , pp . 9497 ( 2008 ) . m.  fiorentino , s.  m. spillane , r.  g. beausoleil , t.  d. roberts , p.  battle , and m.  w. munro , `` spontaneous parametric down - conversion in periodically poled ktp waveguides and bulk crystals , '' opt . express * 15 * , 74797488 ( 2007 ) . d.  achilles , ch . silberhorn , c.  sliwa , k.  banaszek , i.  a. walmsley , m.  j. fitch , b.  c. jacobs , t.  b. pittman , and j.  d. franson , `` photon - number - resolving detection using time - multiplexing , '' j.  mod .  opt . * 51 * , 14991515 ( 2004 ) .
Abstract:

Answer: parametric down - conversion ( pdc ) offers the possibility to control the fabrication of non - gaussian states such as fock states . however , in conventional pdc sources energy and momentum conservation introduce strict frequency and photon number correlations , which impact the fidelity of the prepared state . in our work we optimize the preparation of single - photon fock states from the emission of waveguided pdc via spectral filtering . we study the effect of correlations via photon number resolving detection and quantum interference . our measurements show how the reduction of mixedness due to filtering can be evaluated . interfering the prepared photon with a coherent state we establish an experimentally measured fidelity of the produced target state of 78% .