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: in the last decades point contact spectroscopy has been proved to be a very powerful method to study electron - phonon interaction in metals ( for reviews see , and ) . in the recent years it has been realized that in mesoscopic metallic systems the electron scattering by some dynamical defects may play a major role in understanding the energy relaxation and dephasing processes of the electrons at low temperature ( for reviews see and ) . such defects can also be studied by point contact spectroscopy provided that they are situated in the contact region . these defects can either be magnetic impurities or structural defects where an atom or a group of atoms have a metastable , almost double degenerate groundstate forming a two - level system ( tls ) . the tls spectroscopy is very similar to the phonon spectroscopy , but as the tlss can not propagate like phonons , they are not in thermal equilibrium with the bath . the tls can be formed as defects in the bulk region of crystalline materials , or at dislocations , at surfaces due to impurities and also in amorphous materials . in most of the cases the internal structure of the tlss are unknown and they can behave very differently .    the dynamics of magnetic impurities shows the kondo effect@xcite where by lowering the temperature a kondo resonance is developed at the range of the characteristic kondo temperature , @xmath0 and approaching zero temperature a magnetic singlet is formed by the impurity spin and the electrons . in that case near zero temperature the magnetic singlet acts just like a static impurity and the dynamical processes are frozen out . the subject of slow defects can be described by the conventional theory of tlss and its application to the point contact is the present review . there are , however , many experimental facts which can not be described in that way , and show resemblance to the magnetic kondo anomalies even if their magnetic origin is very unlikely . we make an effort to point at those experiments by comparing the observed behavior with the predictions of the slow tls model . several of these are believed to be due to dynamical defects , different from slow tls . it must be emphasized that there is no theory which is generally accepted and can be applied . earlier it has been proposed that they are fast tlss , which can be described by the so called two - channel kondo model.@xcite that proposal has recently been criticized and debated.@xcite these scatterers now appear rather as `` phantoms '' and the present review does not cover that issue , we , however , present the earlier interpretation just to help in reading the earlier experimental publications .    * . in that case an atom or a group of atoms can have a transition between the two states of lowest energies by direct tunneling or at higher temperature by thermal excitation over the barrier . in the former case those systems are also called as tunneling centers ( for review see @xcite ) . the transition rate can be on the scale of seconds and even on longer scales . in that case e.g. the specific heat may depend on the speed of the measurement as to reach thermal equilibrium takes longer time . the resistivity of such mesoscopic system may jump in time between two or several values on the above time scale , as in small mesoscopic systems the resistivity depends on the position of a single atom even if the two possible atomic positions are much smaller than an atomic distance . the study of point contacts containing tlss in the contact region turned out to be very powerful as the tls may result in zero bias anomalies shown by the @xmath1 characteristics . the electron passing through the point contact can interact dynamically with the tls or it can be elastically scattered , which scattering depends on the actual position of the atom forming the tls . the inelastic scattering may result in a back - scattering of the electron in which case the electron does not proceed from one electrode to the other when the applied voltage is larger than the splitting of the tls and the process shows up as a resistivity minimum at zero bias . in the case of elastic scattering the transmission rate may depend on the position of a single atom and the anomalies can have either sign . the characteristic relaxation times from the tls causing slow telegraph fluctuation of resistance to those being responsible for zero bias anomalies in the @xmath1 curve varies over several orders of magnitude , but for all these so - called slow tlss the average transition time is much longer than the electron - tls interaction time . * the fast tlss are strongly debated whether they can show kondo anomaly.@xcite here we present the previous ideas just to make contact with the extended literature , and the present status of the problem will be discussed at the end of the paper . if a conduction electron experiences several switching occasions one after the other we speak about fast tls . these fast tlss may show kondo anomaly as the tls s two states can be described by a quasi - spin replacing the real impurity spin of the magnetic kondo problem and the scattered electron can have different spherical momenta e.g. @xmath2- and @xmath3-states playing similar role to the real spin of the conduction electron . the spin - flip process of the magnetic kondo problem is replaced in the latter case by electron assisted transition between the two states of the tls . that trasnition was is originally believed to be electron assisted tunneling , which is unlikely according to the present theories,@xcite however an electron assisted transition without tunneling is still possible.@xcite in point contacts containing magnetic impurities the kondo resonance scattering contributes to the back - scattering rate , thus the kondo impurities result in zero bias peak in the dynamical resistivity , which has a width characterized by the kondo temperature , @xmath0 . the sign of the peak is just opposite to those caused by slow tlss . i.k . yanson and his co - workers@xcite and also other groups@xcite studied the zero bias anomalies due to intentionally placed magnetic impurities . considering the kondo effect there is a significant difference between the magnetic case and the tls . in the latter case the conduction electrons , additionally to their angular momenta occurring in the hamiltonian and the coupling , have a further degree of freedom , the real spin , which makes the electronic sea double degenerate although it does not appear in the coupling constants . referring to that extra degeneracy the phenomenon is called two - channel kondo problem ( 2-ck ) in contrast to the magnetic one - channel case ( 1-ck ) ( for earlier review see ) . that extra degeneracy prevents the system to have a singlet ground state at low temperature , it has finite entropy and , therefore , the fermi liquid behavior does not show up as far as a low energy cutoff ( e.g. the energy difference between the two lowest states of the tls , known as splitting ) is reached as the temperature is lowered . + the condition for the formation of a 2-ck ground state is that the energy splitting is negligible compared to the kondo temperature , which is generally not expected . the zero bias anomalies have been believed to be powerful tools to make difference between 1-ck and 2-ck problems ( for a review see ) . for example annealing or electromigration could result in disappearance or modification of the zero bias anomalies in the case of structural defects while such effects are not expected for magnetic impurities.@xcite similar issues also appear in the electron dephasing time which has been intensively debated recently@xcite . another version of the two - channel kondo effect was proposed to explain the zero bias anomalies in uru@xmath4si@xmath4 point contacts@xcite where two degenerate localized electron orbitals of the u atom play the role of the tls which interact with the conduction electrons.@xcite    the issue of the 2-ck problem was responsible for the great interest in the zero bias anomalies due to tls and , therefore , it has a great importance to distinguish between anomalies due to slow tlss and fast scatterers . even if a considerable amount of the discussions are carried out considering fast tlss there is no doubt that in many cases the slow ones are dominant e.g. where the zero bias anomaly in the point contact spectrum has a positive sign.@xcite working on the issues listed above we learned that the literature on the theory of slow tlss in point contacts is very much scattered over many papers and journals , the predictions should be collected from a large number of papers , especially , published in russian and ukrainian physical journals some of them not available at smaller libraries . that experience encouraged us to collect and at several occasions to complete those results . the present review is written in a self - contained manner , thus no further reading is necessary to follow the theory . we want to emphasize that the authors did not consider as a task to give the historically proper complete list of references but they picked up those which are the most appropriate ones to get further information . the model of two level tunneling systems was at first introduced to explain the unexpected low - temperature behavior of specific heat and thermal conductivity in disordered solids which can not be described by phonon excitations.@xcite the model based on two level systems , which was able to explain the linear temperature dependence of specific heat of amorphous dielectrics at low temperatures was suggested by anderson , halperin and varma@xcite and phillips.@xcite this so - called tunneling model assumes that in an amorphous state some atoms or groups of atoms may switch between two , energetically nearly equivalent configurations . this situation is modeled as a double well potential with two stable states differing by energy @xmath5 , and a potential barrier between them with a tunneling rate @xmath6 ( see fig . [ tls ] ) .     between the two positions , and a tunneling probability @xmath6 for crossing the barrier between the two metastable states . @xmath7 and @xmath8 denote the width of the barrier and the distance between the minima , respectively . _ _ ]    the tunneling probability can be approximated in terms of the gamow parameter ( @xmath9 ) : @xmath10 , where @xmath9 is defined by the characteristic height , @xmath11 and width , @xmath7 of the potential barrier , and the mass , @xmath12 of the tunneling particle : @xmath13 furthermore it is assumed that the tlss in a solid are uniformly distributed in terms of the parameters @xmath5 and @xmath9 , that is the density of states @xmath14 for the first excited states of the tlss is constant : @xmath15 experimental studies do indeed support this assumption.@xcite    the two positions of the system are described by a quasispin , and the hamiltonian of the tls is expressed in terms of the quasispin s pauli - operators : h_tls=12 _ , ( b^+_^z_b_+ _ 0b^+_^x_b _ ) , [ h_tls ] where the indexes @xmath16 and @xmath17 correspond to the two states of the tls ( @xmath18 and @xmath19 ) , and @xmath20 and @xmath21 are the creation and annihilation operators for the state @xmath16 . this hamiltonian can be diagonalized by writing the eigenstates as:|e_+&=&|1+|2 +    between the two eigenstates is : e=.the coefficients @xmath22 and @xmath23 can be expressed as : = + = .[munu]for a highly asymmetric tls , where @xmath24 the energy eigenstates are the two positions of the tls , that is @xmath25 and @xmath26 . in the opposite case , where the tunneling energy is much larger than the splitting ( @xmath27 ) the energy eigenstates are delocalized between the two positions , and @xmath28 .      in a disordered material two level systems are not isolated objects , but they interact with their neighborhood . for instance a strain induced by an elastic wave may change the parameters of the two level systems . it is assumed that only the asymmetry parameter @xmath5 is influenced by an elastic wave , thus the interaction between two level systems and phonons can be treated as a perturbation in @xmath5 : @xmath29 where @xmath30 is the coupling constant , and @xmath31 is the appropriate local elastic strain tensor .    in metallic samples the coupling to conduction electrons has to be taken into account as well . the general form of the hamiltonian describing the interaction between the electrons and the tls is:@xcite h_el - tls=__1,_2 , + , , b^+_a^+__2v^__1_2a__1b _ , [ gen_int_tls ] where @xmath32 and @xmath33 are the electron creation and annihilation operators , @xmath34 is the spin index of the electron and @xmath11 is the interaction matrix element , which can be decomposed in terms of pauli operators:@xcite v^__1_2=_i = x , y , z v^i__1_2^i_+v^0__1_2_. [ v_pauli ] the last term stands for potential scattering on the average positioned tls , the term with @xmath35 describes the difference between scattering amplitudes for the two positions of the tls , while @xmath36 and @xmath37 correspond to the electron assisted transitions between the two states of the tls . the interaction hamiltonian is further simplified if the electron - tls interaction time is much shorter than the transition time of the tls , thus the tls stays in a given position during the electron scattering . for these so - called slow two level systems the electron assisted transitions described by @xmath38 and @xmath39 can be neglected , and the hamiltonian is : h_el - tls=__1,_2 , v^z__1_2a^+__2a__1 _ , ^z_b^+_b_+__1,_2,v^0__1_2a^+__2a__1 _ , _ b^+_b _ , [ int_slowtls]the interaction matrix elements between the delocalized energy eigenstates are : e_-|h_el - tls|e_+=__1,_2,a__2^a__1,[matrix_elements_inel]which gives the amplitude of the electron s inelastic scattering on the tls , which is associated with an energy change , @xmath40 . the matrix elements corresponding to the elastic processes are written as : e_+|h_el - tls|e_+&=&__1,_2,a__2^a__1 + e_-|h_el - tls|e_-&=&__1,_2,a__2^a__1.[matrix_elements_elastic]in both cases the energy of the scattered electron is conserved , however the scattering cross section is different if the tls is in the @xmath41 or @xmath42 state .    by simple argumentation the interaction strength , @xmath43 can be estimated as : @xmath44,@xcite as the difference between the electron scattering amplitudes in the two positions must be proportional to the distance of those positions on the scale of the fermi wavelength , and the strength of the scattering , @xmath11 . @xmath45 is very model dependent , e.g.  the tunneling atom can be different from the most of the atoms . the electron assisted transition processes are only important in the case of fastly relaxing two level fluctuators . in this case the hamiltonian given by eq . ( [ gen_int_tls ] ) can be scaled to the two channel kondo problem ( 2ck ) with the hamilton operator : h_int=_tlsv_p_1,p_2 , , , a^+_p_1_a_p_2 , [ 2ck ] where the the momentum dependence of @xmath46 is decomposed into appropriately choosen spherical waves with indices @xmath47 , and @xmath48 ( for details see ) . @xmath49 is a pseudo - spin corresponding to the spherical indices of the electrons , whereas @xmath50 is a pseudo - spin corresponding to the states of the tls , and @xmath11 is the interaction strength . in this model the electron spin is an extra degree of freedom compared to the single channel magnetic impurity kondo problem , and this quantum number is conserved during the interaction and the coupling is independent of that . in contrast to the single - channel kondo problem , which has a non - degenerated , @xmath51 ground state because of the spin compensation cloud formed by the conduction electrons , the ground state of the two channel kondo problem has a non - zero entropy . in the two channel kondo model the conduction electrons overscreen the impurity spin by forming a non - trivial spin state . the low temperature regime of the two channel kondo model was found to have non - fermi - liquid behavior . this paper is devoted to review the literature of the slow tlss and point contacts . the kondo regime of the tlss is out of present consideration , but one can find an up - to - date review of this field in the work of cox and zawadowski.@xcite      in the previous part of this section the scattering rates due to tlss are calculated in born approximation . considering only the screening interaction @xmath43 and ignoring the assisted tunneling rates ( @xmath52 ) kondo@xcite has found power law correction to the tunneling rate which is due to the building up of a screening cloud or friedel oscillation around the changed position of the tunneling atom . that is closely related to the anderson s orthogonality catastrophe@xcite and the x - ray absorption problem ( see e.g.  ) . that phenomenon is associated with creating a large number of electron - hole pairs as the position of the tunneling atom is changed . the renormalized dimensionless coupling constant ( usually denoted by @xmath53 ) can be expressed in terms of the phase shift @xmath54 describing the scattering of the electrons by the atom in the s - wave channel ( see e.g.  ) and the separation distance , @xmath8 between the two positions of the tls @xmath55 ^ 2 , \label{k}\ ] ] where @xmath56 with the spherical bessel function , @xmath57 and @xmath58 for only @xmath2-wave scattering . in the weak coupling limit @xmath59 and @xmath60 that expression is proportional to @xmath61 . in many other publications the notation @xmath16 is used for @xmath53 ( see e.g.  ) . the behavior of a tls coupled to electronic heat bath ( ohmic region ) by the coupling @xmath53 has been studied in great detail ( see e.g.   and ) . in general there are two different regions regarding the coupling strength @xmath53 . here the asymmetry parameter @xmath5 is disregarded and @xmath62 is taken .    * . at @xmath63 the atom has a periodic motion between the two positions . as @xmath53 is increased the oscillation is more damped and the renormalized tunneling frequency @xmath64 is given by @xmath65 where @xmath66 is the high frequency cutoff , which has first believed to be the electronic bandwidth ( see e.g. , but according to the adiabatic renormalization it is rather the energy scale of the lower excitations in the potential well where several excited levels are considered.@xcite * . the motion is incoherent and the atom is spending longer times in a position before tunneling as @xmath53 is increased . many experimental facts indicate ( see next subsection ) that the tlss are in the first region ( @xmath58 ) and even can be in the weak coupling limit ( @xmath67 ) . the main issue is how the quasielastic and inelastic scattering rates are modified by the coupling @xmath53 . that problem was investigated by h. grabert , s. linkwitz , s. dattagupta , u. weiss@xcite considering the neutron scattering by tls . those results are very instructive for the present case , there are , however , essential differences :    * korringa type of relaxations are very important in the relaxation of the tls coupled to the ohmic electronic heat bath . if energy @xmath68 is transferred to the tls the spectrum of the tls is broadened by a relaxation rate proportional to @xmath68 . in the case of point contacts that energy may arise from the energy of the electron passing through the point contact . at finite bias , @xmath11 the system is , however , out of equilibrium and in the stationary state of the tls the averaged energy of the tls is proportional to @xmath11 also . thus , the nonequlibrium situation should be treated by using e.g. the keldysh technique.@xcite * the energy transfer is not determined by the applied voltage , in contrast to the neutron scattering experiment where the energy of the incoming and outgoing neutrons are simultaneously measured . in the present case by small change of the voltage the energy of the extra electrons passing through the point contact are given ( @xmath69 ) but the energy of the scattered electrons is not determined , thus the energy of the outgoing electron shows some distribution . that results in an additional smearing in the spectrum .    in the following we give a brief summary of the neutron scattering case.@xcite very accurate results are obtained for low temperature and energies except the range @xmath70 . that is certainly the interesting range considering especially for the long range tails of the characteristics which are crucial for the interpretations of the experiments.@xcite the calculation was carried out by expanding in the tunneling events and summing up . more comprehensive studies were performed by using numerical renormalization group@xcite for the entire energy range of energy and even including the asymmetry , @xmath5 and later even for nonequilibrium.@xcite . the previous results@xcite have a very good fit , by the formula @xmath71 where @xmath33 , @xmath72 , and @xmath73 are fitted parameters . that clearly shows a double peak structure with energies @xmath74 . the lines are broadened by a korringa type energy relaxation rate proportional to @xmath68 . in the nearly lorentzian lineshape no direct role of coupling @xmath53 was discovered . those results are strongly indicating that in weak coupling region that would be very hard to discover any sign of anomalous relaxation rate in the tunneling characteristics . that situation is further complicated by the averaging over the distribution of the energy of the scattered electrons and the out of equilibrium situation . we find , that the experimental studies of the tunneling characteristics are not appropriate tool to look for the interaction effects , except that @xmath6 is renormalized . that situation is very different from the cases of one and two channel kondo impurities , where the characteristics are determined by the kondo resonances . there are several experimental methods to investigate the kinetics of tlss . the first two methods to be discussed are studying the tlss by their interaction with phonons . in heat release measurements the time dependence of heat release or specific heat is measured : after a sudden cool down of the sample due to the energy relaxation of the tlss it takes a long period of time to reach thermal equilibrium . the second group of measurements investigates the propagation of sound in amorphous materials measuring sound attenuation and sound velocity.@xcite    heat release measurements basically focus on amorphous dielectrics and superconductors , as in metals the contribution of conduction electrons to the heat capacity usually dominates that of tlss , thus the estimation of @xmath75 is more difficult . according to the heat release measurements of kol et al.@xcite in amorphous metal fe@xmath76b@xmath77si@xmath78 the density of the tlss is about @xmath79 . though in acoustic measurements basically the phonons are investigated , if the electron - tls coupling is strong , the results may significantly differ from that on amorphous dielectrics.@xcite a simple description for ultrasound absorption by two level systems is as follows . the occupation numbers for the upper and the lower states of the tls are @xmath80 and @xmath81 , respectively . the probability for ultrasound absorption is @xmath82 , where @xmath14 is the ultrasound intensity , and @xmath16 is the absorption coefficient . the probability of relaxation is @xmath83 , where @xmath84 is the relaxation time of the tls . in balance these two quantities are equal , thus one obtains : @xmath85 the absorbed energy is : @xmath86 it is easy to see , that at large enough intensities ( @xmath87 ) , the energy absorption saturates to a value inversely proportional to the relaxation time ( @xmath88 ) . in insulators and semiconductors the available maximal ultrasound power is enough to drive the absorption into saturation at @xmath89k , but in metals such saturation is only obtained at temperatures as low as @xmath90mk . it implies that the relaxation times in metallic samples are much shorter than in disordered insulators . the interpretation for this experimental observation is , that the electron - tls scattering processes are dominating in the relaxation of tlss . assuming a standard korringa - like electron - tls relaxation , where an electron - hole pair is created , at low enough temperatures one obtains : @xmath91 that is , the ultrasound absorption measurements give the opportunity to make an estimation for the electron - tls coupling parameter ( @xmath53 ) , and the density of the tlss ( @xmath75 ) . in fact the acoustic properties of disordered metals do not show universal behavior as in the case of insulating glasses . there are several theoretical approaches ( the mentioned standard korringa - like relaxation,@xcite the strong coupling theory of vladr and zawadowski,@xcite the electron - polaron effect considered by kagan and prokofev@xcite ) that may explain the aspects of the relaxation of tlss due to conduction electrons.@xcite still , taking into account the simple electron - tls coupling given by eq . ( [ gen_int_tls ] ) , the acoustic measurements can give estimations for the density of two level systems ( @xmath75 ) or the electron - tls coupling constant ( @xmath53 ) . in case of superconductors the korringa relaxation is suppressed by the presence of the superconducting gap , but the gap disappears if a high enough magnetic field is applied on the sample , switching on the tls - electron interaction.@xcite according to the measurements of esquinazi et al.@xcite in the normal and superconducting state of @xmath92 compared with the theory of kagan and prokofev@xcite the coupling parameter is approximately @xmath93 . the acoustic measurements of coppersmith and golding@xcite on the normal conducting amorphous metal pd@xmath94si@xmath95cu@xmath96 estimated the coupling constant as @xmath97 . for the density of the two level systems a lot of results are available in the literature , and it can be stated generally , that it is the same within a factor of less than 10 over a large variety of disordered material , being metallic glasses or dielectric amorphous material . ( for tlss with energy splitting less than 1 kelvin it is approximately 1 - 10ppm . ) a detailed review of heat release and acoustic measurements in disordered material can be found in , and for earlier data see . electron - tls interaction can be studied directly by measuring electric transport in disordered metals . point contact spectroscopy offers the possibility to investigate the properties of even a single two level system centered in the vicinity of the contact . the most spectacular sign of two level systems in point contacts is the so - called telegraph noise : the resistance of the contact is fluctuating between two ( or more ) discrete values on the timescale of seconds . one can estimate the average lifetimes in the excited ( @xmath98 ) and ground state ( @xmath99 ) by recording hundreds of transitions , and fitting the resulting histograms of lifetimes ( @xmath100 and @xmath101 ) with exponential decay functions ( @xmath102 and @xmath103 , respectively ) . the two life times are related to each other by the detailed balance , @xmath104 , where @xmath5 is the asymmetry parameter . if the tls jumps to the other state , the electron screening cloud also needs to rearrange , which makes the jumps of the tls slower . ( the building up of the electronic screening cloud is related to a process which is similar to the x - ray absorption in metallic systems@xcite , and which is also called as electron - polaron effect.@xcite ) from the theoretical point of view this slow - down of the tls motion is treated as a renormalization of the tunneling amplitude:@xcite @xmath105 where @xmath66 is the bath cutoff frequency for which originally the electron bandwidth had been taken,@xcite but later it was proposed that it is replaced by a typical energy of the next higher excitation in the potential well.@xcite golding et al.@xcite studied the two level fluctuation in mesoscopic disordered bi samples . they argued that their samples were in the strong coupling limit , where the coupling of the tls to conduction electrons is very strong compared to the tunneling matrix element , thus the latter can be treated as a small perturbation . in this limit the tunneling of the tls is incoherent , because of the dephasing due to the fast oscillations of the electron bath . here a theoretical scaling function can be set up : @xmath106 , which does indeed agree with the experimental result . ( the energy splitting of the the tls can be tuned by changing the electron density , which can be reached by applying weak magnetic field.@xcite ) in the limit @xmath107 the scaling function is constant , i.e. @xmath108 , that is the coupling parameter can be determined from a single fit , giving @xmath109 for the [ couplingconst ] particular tls measured . the effect of fastly relaxing tls in point contacts can not be resolved as a telegraph fluctuation of the resistance , it causes an anomalous behavior in the voltage dependence of the differential resistance around zero bias , a so - called zero bias anomaly ( zba ) , which will be a subject of detailed discussion in this review . metallic systems , where two macroscopic electrodes are connected via a contact with small cross section are called `` point contacts '' ( pc ) in general , regardless of the actual size of the contact area . it can be generally stated , that the resistance of a point contact is mostly determined by the narrow neighborhood of the junction ; therefore , a pc acts as a `` microscope '' magnifying all kinds of phenomena occurring in the small contact region . + ( a ) ( b ) ( c ) ( d )    several methods have been worked out to produce extremely small contacts between two conducting leads . figure  [ 4methods ] presents the four most important experimental techniques . the first one ( fig . [ 4methods]a ; ) is referred to as mechanically controllable break junction technique ( mcbj ) . here the sample ( practically a piece of metallic wire ) is fixed on the top of a flexible beam , and a small notch is established between the anchoring points . the contact is created in situ at low temperature by breaking the sample with bending the beam , thus one obtains clear and adjustable junction on atomic length scale . the second method ( fig . [ 4methods]b ; ) uses nanolithography to establish a small hole in a silicon nitride insulating membrane by etching . if the etching is stopped just when the hole breaks through , the diameter at the bottom edge remains extremely small ( @xmath110 3 nm ) , well below the usual resolution of lithography ( @xmath110 40 nm ) . in the next step metal is evaporated on both sides , creating a high - quality point contact device . this method can provide extremely stable and clear contacts on atomic length scale , however the contact diameter can not be varied during the measurement . the third drawing ( fig . [ 4methods]c ) shows a similar arrangement to a scanning - tunneling microscope , a so - called spear - anvil geometry : a vertically movable , sharply tapered needle is pressed onto a flat surface . finally , fig . [ 4methods]d shows a simple technique , where the edges of two wires are brought into contact . in arrangements ( a ) , ( c ) and ( d ) usually a differential screw mechanism is used to adjust the contact supplemented with a piezo - crystal for fine tuning . the first application of point contacts was carried out by igor yanson@xcite to investigate electron - phonon scattering in nanojunctions . he found , that the point contact spectrum , obtained as the second derivative of current with respect to the bias voltage ( @xmath111 ) contains structure due to the electron - phonon interaction described by the eliashberg function @xmath112.@xcite this simple method for measuring electron - phonon interaction spectrum became a popular application of pc spectroscopy , however , it can be used to probe other electron scattering processes as well , like electron - tls interaction , which is the central topic of this paper . theoretically , point - contacts are considered as two bulk electrodes connected through a narrow constriction . the simplest , and most commonly used pc model is presented in fig . [ pcmodels.fig]a . this so - called opening - like point contact is an orifice with diameter @xmath8 in an infinite isolating plane between the two electrodes . another extreme limit is the channel - like pc : a long , narrow neck between the bulk regions with the length being much larger than the diameter , @xmath113 ( fig .  [ pcmodels.fig]b ) . the crossover between the both cases can be obtained by considering the point - contacts as rotational hyperboloids with different opening angles ( fig . [ pcmodels.fig]c ) . in most of the cases the shape of the pc does not influence the character of physical processes in the constriction radically , and the main parameter is the ratio of the contact diameter ( @xmath8 ) and other characteristic length scales in the system . three fundamental length scales are the mean free paths connected to different scattering processes ( @xmath114 ) ; the fermi wavelength of electrons ( @xmath115 ) ; and the atomic diameter ( @xmath116 ) .    '' in an insulating screen between two metallic half - spaces . ( b ) two bulk regions are connected with a long , narrow conducting neck . ( c ) rotational hyperboloid.__,title="fig:",height=170 ] '' in an insulating screen between two metallic half - spaces . ( b ) two bulk regions are connected with a long , narrow conducting neck . ( c ) rotational hyperboloid.__,title="fig:",height=170 ] '' in an insulating screen between two metallic half - spaces . ( b ) two bulk regions are connected with a long , narrow conducting neck . ( c ) rotational hyperboloid.__,title="fig:",height=170 ] + ( a ) ( b ) ( c )    if the fermi wavelength or the atomic diameter becomes comparable to the contact size , we speak about _ quantum point contacts _ and _ atomic - sized point contacts _ , respectively . these systems are reviewed in detail in . in this paper we consider contacts , that are neither _ atomic _ nor _ quantum _ , but that are small enough compared to certain mean free paths , and the calculations are basically performed for an opening type geometry . concerning the mean free paths , we have to make difference between the elastic and inelastic scatterings . under usual experimental conditions the elastic mean free path ( @xmath117 ) is smaller than the inelastic one ( @xmath118 ) . here the inelastic mean free path is the _ length of the path _ an electron travels between two inelastic scatterings ( @xmath119 ) . mostly the important parameter is not @xmath118 , but the inelastic diffusive length , that is the average _ distance _ an electron can diffuse between two inelastic scatterings : @xmath120 . if the contact diameter is much smaller than any of the mean free paths @xmath121 , we speak about a _ ballistic _ contact . in this case the electron travels through the constriction without any scattering ( except for the reflection on the walls ) . on the other hand , if @xmath122 , the electron makes a diffusive motion in the contact , and accordingly we speak about the _ diffusive _ regime . at contact diameters exceeding the inelastic diffusive length ( @xmath123 ) the excess energy of the electrons is dissipated inside the constriction , which causes a considerable joule heating in the contact . this limit is called _ thermal _ regime . in the following subsections these different limits are discussed . in many cases the system is characterized by the knudsen number , @xmath124 , which was first introduced for the problem of the gas outflowing from a tank through a hole,@xcite but in the recent decades it has been used to characterize point contacts as well .      first the diffusive contacts are treated , for which the electric potential ( @xmath125 ) can be determined by classical equations . ( for a general discussion see ) if the mean free path of the electrons is much shorter than the dimension of the contact ( @xmath126 ) then the current density , @xmath127 is given by ohm s law in terms of the electric field , @xmath128 or the electric potential , @xmath125 : @xmath129 where @xmath34 is the conductivity of the metal . furthermore due to the charge neutrality in metals the continuity equation holds : @xmath130 if the conductivity is considered to be constant , these two equations result in the laplace equation for the electric potential : @xmath131    in this phenomenological approach the scattering processes are included in the conductance , which is inversely proportional to the mean free path , @xmath132 . this treatment does not distinguish between elastic and inelastic electron scattering . in the further treatment of the structures in the dynamical conductance ( sec .  [ scatteringinpc ] ) the inelastic scatterings play crucial role , therefore , instead of the phenomenological theory the kinetic equations must be used . here the solution of the laplace equation is presented in order to determine the linear resistance of a diffusive contact , and to determine the position dependence of the electric potential . furthermore , we note that the laplace equation is only applicable if @xmath34 is constant is space . this assumption is certainly not valid if the temperature of the contact neighborhood is inhomogeneous , which happens in the thermal regime . this situation is treated in the next subsection . here we assume that the inelastic diffusive length is much larger than the dimension of the contact , thus the temperature of the contact region is constant . the laplace equation should be solved with the boundary conditions : ( z ) . [ bound ] this problem was first solved by maxwell.@xcite due to the geometry of the problem the laplace equation is most easily solved in an elliptical coordinate system demonstrated in fig . [ elliptic_system ] ( see e.g. , ) . = 8.0 cm 0.7truecm    the transformation between a traditional cylindrical coordinate system @xmath133 and the elliptical coordinates @xmath134 is determined by the following equations of the ellipses and hyperbolas : + z^2&=&1 + r^2-z^2&=&1 , [ ell_hyp ] where @xmath135 is the radius of the contact ; @xmath136 ; and @xmath137 . using eq . ( [ ell_hyp ] ) the old coordinates can be expressed in terms of the new ones : r = a , z = a . in this coordinate system the orifice is given by the @xmath138 surface , while the isolating layer is defined by @xmath139 . the next task is to determine the components of the metric tensor in the new coordinates : s^2&=&dr^2+dz^2+r^2d^2 + or , using @xmath9 and @xmath22 , s^2&=&a^2 ^ 2+^2d^2 + a^2 ^ 2+^2d^2 + a^2(1+^2)(1-^2)d^2 . the prefactors are the components of the metric tensor in the new ( also orthogonal ) coordinate system : g_=a^2 ^ 2+^2 , g_=a^2 ^ 2+^2 , g_=a^2(1+^2)(1-^2 ) . substituting these into the form of the laplace operator expressed by the general orthogonal coordinates : = 1 , one can get easily = 1 . taking the geometry of the problem into account it is obvious that = 0 . the boundary condition eq . ( [ bound ] ) can only be satisfied if at @xmath140 the potential is independent of @xmath22 , that is @xmath141 . as an ansatz we extend this condition to arbitrary @xmath9 values which means that the equipotential surfaces are the @xmath142 surfaces . as these surfaces are orthogonal to the @xmath143 plane it is obvious that electric field has no normal component at the insulating layer . after this consideration the laplace equation takes the following simple form = 0 . [ laplace_apert ] the solution of this equation with the boundary conditions is : ( ) = -v ( ) . [ phi_lambda ] along the @xmath144-axis ( @xmath145 ) , the potential is changing as ( z)=-v , [ eq31 ] which is plotted in fig . [ potent_fig ] . = 0.6 0.7truecm    inside the opening ( @xmath146 ) the electric field has only @xmath144 component : e_z=-=-1|_=0=v1.[ez]it should be noted that the field is the largest at the edge of the opening ( @xmath147 ) and most of the current flows through that region . the total current flowing through the contact is : i=_0^a 2rdre_z(r)=2av , and thus the resistance of a diffusive point contact is : r=1 . this calculation was performed for an opening type contact . in case of a channel type geometry with @xmath113 the potential drops in the channel , and the resistance is easily written using ohm s law : @xmath148      if the dimension of the contact is larger than the inelastic mean free path , then the electron s excess energy is dissipated in the contact region , and joule heating takes place . in this case the contact center can be considerably overheated compared to the bath temperature , and thus the conductivity is also position dependent : @xmath149 . this phenomenon can be treated by considering the equations both for the electrical and the heat conduction . the electrical and thermal current densities are : @xmath150 where @xmath151 is the heat conductivity and @xmath152 is the position dependent temperature . the continuity equations for the electrical and thermal current are determined by the charge neutrality and the joule heating respectively : @xmath153 furthermore , we assume the validity of the the wiedemann - franz law : = , [ wiedemann ] where @xmath154 is the lorentz number . combining these equations we get : @xmath155 from these equations it can be easily seen that the temperature is generally related to the electric potential as : @xmath156 the constant term is determined by the boundary condition : @xmath157 thus : @xmath158 in an opening type contact the potential is zero at the contact surface , so the local temperature is written as : @xmath159 according to this relation a bias voltage of @xmath160mv already heats up a contact from cryogenic temperatures ( @xmath161k ) to room temperature . this feature has crucial importance in experimental studies . if the measurements are performed on a large contact or at high bias voltages , it can easily happen that the @xmath1 curve presents the temperature dependence of the conductance instead of the spectroscopic features,@xcite which will be discussed in sec . [ scatteringinpc ] . if the linear dimension of the contact is much smaller than the minimum of the mean free paths the electrons are travelling along straight trajectories , and in the first approximation no scatterings are taken into consideration . this system is described by the equations of semiclassical dynamics , so the space and momentum dependent distribution function , @xmath162 is to be determined beside the electric potential , @xmath163 . in this limit the electrons `` remember '' which side of the contact they are coming from , thus the distribution function at position @xmath164 can be expressed as a sum of the terms corresponding to electrons coming from the left or the right side of the contact.@xcite the distribution function can be determined by solving the homogeneous boltzmann equation which contains no collision integral:@xcite @xmath165 where @xmath128 is the electronic field , @xmath166 , and @xmath167 is the velocity of the electron with momentum @xmath168 . note that the electronic charge @xmath169 is _ positive _ , so the charge of an electron is @xmath170 . far from the constriction the boundary condition for the distribution function is that it has to tend to the equilibrium distribution function . @xmath171 where @xmath172 . if we apply @xmath11 voltage on the sample , the boundary condition for the potential @xmath125 takes the form @xmath173 the solution of the boltzmann equation can be found by using the trajectory method.@xcite ( the trajectories of the electrons are considered to be straight because we are interested in the limit @xmath174 and the current we are searching for , depends linearly on the electrical field . it can be shown that the bending of trajectory only contributes to the current in the higher order of the field . ) in this way @xmath175 holds combined with the charge neutrality condition @xmath176=0 , \label{charge}\ ] ] where @xmath177 denotes the integration over the trajectory of the electron coming from the distant reservoir to the contact region at point @xmath164 . the distribution function of an electron satisfying the conditions ( [ charge ] ) , ( [ trajectory ] ) and ( [ limit1 ] ) takes the form in the linear order of the field @xmath178 which is constant along the trajectory as the energy @xmath179 is conserved . the expression @xmath180 describes the different value of the chemical potential in the right and left hand sides of the constriction ( @xmath181 ) : @xmath182 ( for an illustration see fig . [ eta.fig ] . ) that can be also expressed in terms of the solid angle , @xmath183 at which the opening of the contact is seen from point @xmath164 : @xmath184 a simple visualization of the distribution function is presented in fig . [ distr_ball ] . = 0.25 0.7truecm    = 0.7 0.7truecm    the potential can be derived by substituting the form ( [ distrfunc ] ) of the distribution function into the neutrality condition ( [ charge]):@xcite 2e^3=0 , where the integral for equienergetic surfaces is easily expressed by the solid angle @xmath185 , and just an energy integral remains : = 0 , [ neutr2 ] where the upper and lower signs correspond to the cases @xmath186 and @xmath187 , respectively . the equation can be solved easily by using the identity ( f_0(+a)-f_0())d =- a[identity]for an arbitrary energy shift @xmath33 . we can get the potential in the whole space as:@xcite @xmath188 \mathrm{sign}(z ) . \label{potent_op}\ ] ] here @xmath185 is the solid angle at which the contact is seen from the position @xmath164 . the potential depending on the distance measured from the contact along the @xmath144 axis is shown in fig . [ potent_fig ] ( solid curve ) . the current density , @xmath127 is determined by the distribution function as follows : ( ) = -2e . [ currdens ] the current through the contact is calculated by integrating the @xmath144-component of the current density over the area of the contact , @xmath189 : @xmath190 at the contact surface half of the electrons go to the left with a distribution function @xmath191 and half of the electrons goes to the right with a distribution function @xmath192 , thus the current is written as i =- e_a da__=(f_0(+)-f_0(- ) ) , where the integral @xmath193 is taken over an equienergetic surface in the @xmath168-space . at low enough temperature and voltage ( @xmath194 and @xmath195 ) the expression @xmath196 equals @xmath197 according to the identity ( [ identity ] ) , so the current is written as : i=()_fsv,[sharvin]where @xmath189 is the area of the contact , @xmath198 is the area of the fermi surface , @xmath199 and @xmath200 denotes the average taken over the fermi surface . the conductance of a ballistic point - contact given by the formula ( [ sharvin ] ) is known as the sharvin conductance.@xcite for free electron gas the sharvin conductance is simplified as : g_s= ( ) ^2 . [ sharvinconductance ]      in the intermediate region between the diffusive and ballistic regime an interpolating formula can be set up by numerically solving the boltzmann equation for arbitrary ratio of the contact diameter and mean free path , @xmath114:@xcite @xmath201 where @xmath202 is a numerically determinable monotonous function , @xmath203 ; @xmath204 . note that the first term is exactly the sharvin resistance by putting the drude resistivity into the formula , @xmath205 , thus it is actually independent of @xmath114 . this formula provides the possibility to estimate the contact diameter from the contact resistance . the voltage applied on a point contact results in a nonequilibrium distribution of the conduction electrons in the contact region . an electron coming from the right reservoir has an energy larger by @xmath197 than those coming from the left reservoir . this energy can be released through inelastic scattering processes , which can happen in such a way that an electron that has already crossed the contact is scattered back through the opening . this so - called backscattering reduces the current , thus the energy dependence of the scattering probability can be traced by measuring the nonlinearity of the current as the bias voltage is varied . this phenomenon was first used by igor yanson@xcite ( for a review see ) to study the phonon spectra and the electron - phonon interaction and since that it is widely applied .    in this section we discuss the nonlinearities in the current voltage characteristics due to the scattering on slow tlss , which show strong similarities to localized phonons . there is , however , an essential difference between the two cases , while the phonons contribute to the voltage region up to @xmath206mev , the tlss manifest themselves below or even well below @xmath207mev due to their characteristic energies . the main contribution to the @xmath1 characteristics comes from the close region of the contact , therefore a microscopic process like a transition of one atom from one position to another in the contact region occurs as a significant measurable change in the current . considering the tlss the main advantage of the point contact spectroscopy is that in case of small enough contacts one can investigate just one or few scattering centers . the current through the contact can be derived by solving the classical stationary boltzmann equation @xmath208 where @xmath209 is the collision integral for the electron with momentum @xmath168 at position @xmath164 , which is assumed to be local in real space .    in the case of bulk phonons the collision integral is @xmath210 \delta(\varepsilon_{{\z{p}}+{\z{q}}}-\varepsilon_{\z{p}}-\omega_{\z{q}}^\alpha)+ \nonumber\\ & { } & \left[f_{{\z{p}}-{\z{q}}}(1-f_{\z{p}})n_{\z{q}}^\alpha - f_{\z{p}}(1-f_{{\z{p}}-{\z{q}}})(n_{\z{q}}^\alpha+1)\right ] \delta(\varepsilon_{{\z{p}}-{\z{q}}}-\varepsilon_{\z{p}}-\omega_{\z{q}}^\alpha)\bigg\ } , \label{eq70}\end{aligned}\ ] ] where @xmath211 is proportional to the squared matrix element of the electron - phonon coupling , @xmath212 is the phonon occupation number for momentum @xmath213 and energy @xmath214 in phonon branch @xmath16 . there is another term of the collision integral which is due to the elastic ( impurity ) scattering but that is not given here . the collision integral vanishes far from the point contact because well inside the electrodes thermal equilibrium exists . the situation is further simplified when the following assumptions are valid : ( i ) the electron - phonon interaction is homogeneous in the real space ; ( ii ) the phonon distribution corresponds to thermal equilibrium . in that case the information on the contact geometry and also the elastic scattering due to impurities in the dirty limit can be expressed by a geometrical factor @xmath215 in the expression of the current which depends only on the momenta @xmath168 and @xmath216 of the incoming and outgoing electrons . that factor is frequently called as the @xmath53-factor ( see e.g.  ) . in the case of phonons the @xmath53-factor for a ballistic opening like contact is written as : @xmath217 where @xmath218 @xmath219 are the unit vector parallel with @xmath168 and @xmath216 , while @xmath220 and @xmath221 are the @xmath144 components of these unit vectors , where the @xmath144 direction is the axis of the contact . the @xmath53-factor can be also calculated for a diffusive contact : @xmath222 where @xmath117 is the elastic mean free path , and @xmath8 is the contact diameter .    at zero temperature the logarithmic derivative of the resistance can be expressed by a function @xmath223 , where @xmath224 is the density of states of phonons and @xmath225 is the modified electron - phonon coupling which contains the @xmath53-factor . without that factor it is just the electron - phonon coupling strength and in that case @xmath226 is the eliashberg function known in the theory of superconductivity . the main result is @xmath227 in case of single crystals @xmath226 depends also on the orientation of the crystal , which makes it possible to study the anisotropy of the phonon spectrum.@xcite measuring @xmath228 the phonon spectra were determined for many different cases and the structures in them were identified as the details of the spectra.@xcite those structures , however , are superimposed on a continuous background .    until now , it is assumed that the phonons are in thermal equilibrium , which is the correct assumption as far as the phonons at the contact region arrive from a large distance without any collision where the thermal bulk distribution is realized . that is , however , not the case when in the contact region the phonon mean free path is comparable with the contact size . then the phonons generated by the non - equilibrium electrons in the contact region can not escape thus the phonons are also out of equilibrium . the measure of non - equilibrium is limited by the energy relaxation time of the phonons @xmath229 . the non - equilibrium localized phonons contribute to the background resistance.@xcite that can be studied by measuring the amplitude of the background as a function of the frequency of the applied voltage@xcite and at high frequency ( @xmath230 ) that is decreasing . these phenomena will be studied in the context of scattering on tlss as well .      in the following the boltzmann equation will be solved considering only the electron - tls interaction . it will be assumed that the density of the tlss is low thus an electron scatters only once in the orifice region and the contribution of the different tlss are additive . first the ballistic limit is treated , where both the elastic and inelastic mean free paths are larger than the size of the contact . the collision term related to a single tls is very similar to the case of phonons , the momentum is , however , not conserved . furthermore , the tls has only two states with population @xmath231 and @xmath232 ( @xmath233 ) corresponding to the tls energy levels @xmath234 and @xmath235 ( see sec . [ tlsinsolids ] ) . the situation is different from those cases of phonons , where the phonon mean free path is large compared to the size of orifice , @xmath236 and the phonon distribution is the bulk one which is in thermal equilibrium with the contact regions . in the case of tlss the interaction is localized in space , thus the tlss in the contact region are decoupled from the bulk region and they can be considerably out of thermal equilibrium due to the voltage drop at the contact . therefore , the situation is similar to the special case of localized phonons with @xmath237 shortly discussed in the introductory remarks of this section .    the collision term similar to eq . ( [ eq70 ] ) can be written as the sum of inelastic and elastic terms as @xmath238 . in the inelastic case the tls jumps to its other state due to the scattering , and thus the energy of the electron changes by @xmath239 . the strength of the interaction is characterized by @xmath240 , which is obtained from the interaction matrix elements ( eq . [ matrix_elements_inel ] ) by fermi s golden rule as @xmath241 the interaction strength is symmetric in the incoming and outgoing momentum , and in the neighborhood of the fermi surface it depends only on the unit vectors @xmath218 and @xmath219 , i.e.@xmath242 . the collision integral for the inelastic scattering is written as : @xmath243 where @xmath244 is the place of the tls . the current correction due to the elastic scattering processes can be written similarly with two remarks : ( i ) in an elastic process the energy of the tls does not change , thus @xmath245 must be inserted ; ( ii ) the scattering cross section is different for the tls being in the `` @xmath246 '' and the `` @xmath247 '' state . according to the matrix elements in eq . ( [ matrix_elements_elastic ] ) , the two scattering strengths are given as : @xmath248 after introducing these notations the correction to the collision integral for the elastic scattering is : @xmath249\delta(\varepsilon_{{\z{p}}}-\varepsilon_{{\z{p}}'})\label{irtlsel}.\end{aligned}\ ] ]    the electron distribution function @xmath250 and the electric potential is expanded in terms of @xmath251 . at a large distance measured from the orifice the potential is constant and @xmath252 . thus , at @xmath253 the electrons are in thermal equilibrium . the distribution function @xmath250 and the potential @xmath254 can be expanded as @xmath255 and similarly the electric field @xmath256 where the upper indices @xmath257 , @xmath258 label the order in the strength of the electron - tls coupling . the zero order terms have been previously calculated and the following treatment is restricted to the first order terms . the previous results are @xmath259 and according to eq . ( [ potent_op ] ) @xmath260 for a round orifice perpendicular to the @xmath144-axis , where @xmath261 is the fermi function . the first order terms in the boltzmann equation given by eq . ( [ boltzmann_new ] ) can be written as @xmath262 where the label @xmath257 of the collision term indicates that the collisions are calculated with the distribution functions of zeroth order , @xmath263 . the change in the potential is determined again by the neutrality condition given by eq . ( [ charge ] ) for @xmath263 , which is @xmath264 the boundary conditions are @xmath265 and @xmath266 in order to determine @xmath267 and @xmath268 the trajectory method is used ( see e.g.  ) dealing with phonons . electron trajectories are considered where the electron comes from the left or right and arrives to the plane of the orifice at time @xmath269 at point @xmath164 with momentum @xmath168 to calculate the current flowing through the orifice . these trajectories correspond to zeroth order . moving along the trajectory @xmath270 the distribution function varies as @xmath271 where @xmath272 , and @xmath273 . this equation can be integrated and using the boltzmann equation one finds for @xmath267 at @xmath269 @xmath274 here the collision term describes the adding or taking off electrons to the trajectory arriving at the contact at @xmath269 . in the expression of @xmath267 only the terms linear in the voltage are kept , thus @xmath275 in the above equation can be approximated by the zeroth order term @xmath276 . taking into account that @xmath277 the first term of the integral can be performed , thus @xmath278 the equation of electrical neutrality taken at @xmath269 with @xmath279 and @xmath280 combined with the equation above determines @xmath281 @xmath282 the general expression of the electrical current flowing through the surface of the contact is @xmath283 where @xmath284 is the surface element of the contact and the integral is taken over the orifice , the factor @xmath285 is due to the electron spin . the change @xmath286 in the total current due to the presence of the tls can be separated as @xmath287 and @xmath288 can be further split , whether the electron scattering is elastic or inelastic @xmath289 making use of the expression ( [ eq : f1.2 ] ) for @xmath267 only the collision term contributes to the current because the other term is even in the momentum . the integral according to the time can be transformed to the one along the path using @xmath290 where @xmath291 is the element of the path . the expression for the change of the current due to the presence of collision on tls is @xmath292 the path of the integral is changed due to the collisions but it may contain reflections by the insulating surface of the contact . the next step is to introduce the collision term due to a single tls at position @xmath244 , given by eq . ( [ irtlsin ] ) . as the paths arriving at the surface element of the orifice @xmath293 are straight lines , the volume element for the scattering event is @xmath294 where @xmath295 is the angle between the path and @xmath144 direction . now the integration over different paths and position of the collision gives @xmath296 where @xmath297 is defined by @xmath298 and the momentum integral is restricted to the solid angle @xmath299 at which the contact can be seen from the tls at position @xmath244 . furthermore , @xmath300 as the electron passes the contact from the tls ( for an illustration see fig . [ pctls1 ] ) . the final expression is obtained using @xmath301 where @xmath302 is the solid angle element in the momentum space and @xmath303 is the conduction electron density of states for one spin direction . assuming a spherical fermi surface @xmath304     and the angle @xmath295 is also indicated . _ _ ]      using the expression of the collisison term @xmath305 ( [ irtlsin ] ) the change in the current due to the inelastic scattering can be written as @xmath306 \delta(\varepsilon+e-\varepsilon ' ) + \left[f_{\z{p}'}^0(1-f_{\z{p}}^0)n_+-f_{\z{p}}^0(1-f_{\z{p}'}^0)n_-\right ] \delta(\varepsilon - e-\varepsilon')\big\},\end{aligned}\ ] ] where the distribution functions @xmath307 are given by eq . ( [ eq : f0 ] ) and taken at @xmath308 , and the dimensionless notation @xmath309 is introduced . the potential @xmath310 appearing in these expressions can be dropped as the energy integral variables can be shifted . the essential contribution is given by @xmath311 , which has a positive sign if the electron at @xmath308 arrives from the left contact and a negative sign if it is arriving from the right contact .    the final expression is @xmath312 \delta(\varepsilon+e-\varepsilon')+\nonumber\\ & & \big[f_0(\varepsilon'+{ev\over2}\eta(\z{p}'))(1-f_0(\varepsilon+{ev\over2}\eta(\z{p})))n_+- f_0(\varepsilon+{ev\over2}\eta(\z{p}))(1-f_0(\varepsilon'+{ev\over2}\eta(\z{p}')))n_- \big ] \delta(\varepsilon - e-\varepsilon')\bigg\}.\nonumber\\\end{aligned}\ ] ]    this expression can be divided into two parts . if the electron with momentum @xmath216 arrives from the same side of the contact ( @xmath313 ) as the unscattered electron with momentum @xmath168 , thus @xmath314 , then the electron is scattered forward , while for @xmath315 it is scattered backward ( see fig . [ fig : paths ] ) . the total current can be divided according to that , thus @xmath316     signs are also indicated.,scaledwidth=70.0% ]    the forward scattering cancels out . that can be seen by looking e.g.  the first and the last term in the previous expression and changing @xmath317 , @xmath318    in the calculation of the backscattering the case @xmath319 is taken first . then the electron contributing to the current originally comes from the left ( @xmath320 ) and after the scattering goes from right to the left ( @xmath321 ) . @xmath322 the factor due to the angular integrals @xmath323 plays the role of the geometrical factor for the bulk phonons . only the backscattering contributes like in the case of phonons as it can be seen from eq . ( [ kfactorball ] ) for ballistic regions where the @xmath53-factor contains the factor @xmath324 .    in the following the geometrical factor is introduced which depends on the position of the tls , @xmath244 , and it is independent of the strength of the interaction @xmath325 where @xmath326    if the momentum dependence of the interaction is ignored then @xmath327 and @xmath328 is a pure geometrical factor . further on this simplification is assumed . if the tls is positioned on the geometrical axis of a circular contact then the above integral can be evaluated as @xmath329    in the second part of the integral in eq . ( [ diin.right ] ) the variables are changed as @xmath317 and the integral with respect @xmath330 is performed , furthermore , in some of the integrals the variable is shifted as @xmath331 . the result is @xmath332n_-\nonumber\\ & + & \left[f_0(\varepsilon+e+ev)-f_0(\varepsilon+e - ev)\right]n_+ \big\}.\label{eq:112}\end{aligned}\ ] ]    now the tls on the left is considered ( sign@xmath333 ) . then for backscattering the signs of @xmath334 and @xmath335 are reversed ( see fig .  [ fig : paths]b ) that is equivalent of changing the sign of @xmath11 ( @xmath336 ) in eq . ( [ diin.right ] ) and then the expression in the curly bracket also changes its sign as it can be seen by inserting @xmath317 . that negative sign is cancelled by sign@xmath333 , thus the total current is unchanged . therefore , the current does not depend on whether the tls is on the right or left as it should be . after inserting the fermi functions and performing the energy integral , the following result is obtained for the inelastic current correction : @xmath337 where @xmath338- { e - ev\over { 2}}\bigg[\coth({e - ev\over { 2k_bt}})-1\bigg]\right\},\nonumber\\ d=\frac{2e}{\hbar}wk_{\z{r}}\left\ { { ev - e\over { 2}}\bigg[\coth({ev - e\over { 2k_bt}})-1\bigg]- { -e - ev\over { 2}}\bigg[\coth({-e - ev\over { 2k_bt}})-1\bigg]\right\}.\end{aligned}\ ] ] these expressions are simplified at @xmath62 as @xmath339    if the tls was in thermal equilibrium with the bath , at zero temperature @xmath340 and @xmath341 would hold . in that case the voltage dependence of the conductance is easily evaluated , it contains a jump - like decrease at the excitation energy of the tls : @xmath342 where @xmath343 . that result means if the tls is in the ground state @xmath341 then the tls can be excited only by a voltage @xmath344 , where the backscattering reduces the current . the second derivative has a dirac delta peak at the excitation energy @xmath345 this result is similar to the phonon result , that is the @xmath346 shows the structure of the excitation spectrum of a scatterer in the contact region . however , for a tls situated in the contact region the occupation of the upper level is not zero , as the system is driven out of equilibrium . due to this feature a nontrivial contribution occurs in the background resistance beside the spectroscopic feature at the excitation energy . this can be evaluated after calculating the voltage dependence of the occupation number , @xmath347 , which is performed later in this section . a similar calculation can be performed for the current correction due to the elastic scattering on tlss by using the collision integral @xmath348 ( [ irtlsel ] ) instead of @xmath349 ( [ irtlsin ] ) . the formulas for the elastic scattering can be derived easily by inserting @xmath245 and modifying the scattering strength in the results for the inelastic case . for elastic scattering the strength of the interaction can be different for the tls being in the `` @xmath246 '' and `` @xmath247 '' state , thus in eq . ( [ cd ] ) @xmath350 must be replaced by @xmath351 ^ -$ ] and @xmath351^+$ ] in the expression for @xmath352 and @xmath353 , respectively . after these considerations the result is : @xmath354 ^ -+n_+\cdot[wk_\z{r}]^+\right).\ ] ] note , that in the elastic case the voltage dependent `` @xmath355 '' terms cancel out , thus the only voltage dependence comes from the voltage dependence of the occupation number . if the tls is in thermal equilibrium with the bath , then @xmath356 holds , thus the elastic current correction is a voltage independent constant term . similarly , if the scattering cross sections for the two states of the tls are equal , then the voltage dependence cancels out due to the condition @xmath233 . in the following we use a simplified notation for the elastic current correction : @xmath357 where @xmath358 and @xmath359 are considered as the reduction of the conductance due to the tls being in the `` @xmath246 '' and the `` @xmath247 '' state , respectively.@xcite the average change in the conductance due to the elastic scattering is characterized by @xmath360 , which is written as follows according to formulas ( [ diel ] , [ kr ] , [ kisw ] , [ wel ] ) : @xmath361 ^ 2+\left[(\mu^2-\nu^2)v^z\right]^2\right),\ ] ] where @xmath45 and @xmath43 are the the interaction matrix elements from eq . ( [ int_slowtls ] ) assuming isotropic scattering . the difference between the conductances corresponding to the two states of the tls are written as : @xmath362    again , the voltage dependence of the current correction can only be determined after calculating the occupation number as a function of bias voltage , which is done in the following part . it is important to note that in case @xmath363 the scattering strength is symmetric for the two states , thus @xmath364 .      in the following @xmath231 will be calculated in the nonequilibrium case , where it is determined by the interaction with the nonequilibrium electrons . the transition probability from the `` @xmath247 '' to the `` @xmath246 '' state due to an electron scattered from momentum state @xmath168 to @xmath216 has been considered in the collision term ( [ irtlsin ] ) and it is : @xmath365 and similarly : @xmath366 the transition probabilities for the tls are obtained after integrating for the electron momenta : @xmath367 where the expression ( [ eq : f0 ] ) for @xmath368 is used . the kinetic equation for @xmath231 is : @xmath369 in the fermi functions the electron momentum is only involved in @xmath334 , which tells which reservoir the electron at position @xmath244 with momentum @xmath168 is coming from . therefore the integral can be performed by separating the regions of integration into four cases with respect to the possible values of @xmath334 and @xmath335 . in all of these four regions the factor containing the fermi functions is independent of @xmath370 and @xmath371 . there is an essential simplification if @xmath372 is replaced by a momentum independent averaged one , @xmath373 . take e.g.  the case where the tls is on the right hand side of the contact ( @xmath374 ) and introduce the solid angle @xmath299 at which the opening can be seen from point @xmath244 . in this case the integrals for the different regions are written as : @xmath375 this simplification is justified only after taking an average over large number of tlss but that simplifies the calculation and makes the result more transparent . in this simplified case the kinetic equation for the tls is obtained after calculating the energy integrals for the different cases . ( as before , @xmath376 is eliminated by shifting the integral variables . ) the final result is independent whether the tls is on the right or left , and it can be written in the form : @xmath377 where the coefficients @xmath189 and @xmath72 are : @xmath378+\nonumber\\ + \left ( \left(1-\frac{\omega_\z{r}}{4\pi}\right)^2 + \left ( \frac{\omega_\z{r}}{4\pi}\right)^2\right ) \frac{e}{2 } \left[\coth(\frac{e}{2k_bt})-1\right ] \bigg\ } , \nonumber\\ b={w\over\hbar}\bigg\{\left({1-{\omega_\z{r}\over { 4\pi}}}\right){\omega_\z{r}\over { 4\pi}}\left[{{ev - e\over{2}}\left({\coth ( { ev - e\over{2k_bt}})-1}\right)+ { -ev - e\over{2}}\left({\coth ( { -ev - e\over{2k_bt}})-1}\right)}\right]+\nonumber\\ + \left({\left({1-{\omega_\z{r}\over { 4\pi}}}\right)^2 + \left({{\omega_\z{r}\over{4\pi}}}\right)^2}\right ) { -e\over 2}\left[{\coth ( { -e\over{2k_bt}})-1}\right]\bigg\},\end{aligned}\ ] ] respectively . the relaxations are of korringa types , which depend on @xmath11 only if the two electrons involved have different chareacters in terms of @xmath334 . the stationary value of @xmath231 is obtained as @xmath379 which is an even function of @xmath11 . at @xmath62 this expression simplifies essentially : @xmath380 where @xmath381 . @xmath231 is a continuous function of @xmath11 . far from the opening @xmath382 thus @xmath383 at voltages for which @xmath384 as the electron gas is in thermal equilibrium . for @xmath385 and @xmath386 @xmath387 as for large enough voltage the levels are equally occupied @xmath388 . if the tls is in the middle of the contact , then @xmath389 . if the tls is far away from the contact then it is in thermal equilibrium , which can be obtained by taking @xmath390 and the deviation from equilibrium is large when the tls is in the middle of the contact region @xmath391 . here other relaxation mechanism different from the scattering of electrons is not taken into account as the generation of bulk phonons is very weak as the relevant phase space is very small at low energies . at large concentration of tlss the tlss are interacting and the collective effects may modify the stationary values of the occupation numbers . now the inelastic and elastic contributions to the current and the conductance are calculated using the stationary occupation numbers obtained . the expressions will be given when the tls is just in the middle of the round opening , thus @xmath391 ( @xmath392 ) , but the calculation for arbitrary @xmath393 can be easily performed as well .    at zero temperature the inelastic current correction is obtained by inserting the stationary value of @xmath231 given by eq . ( [ nt0 ] ) into ( [ iinel2 ] ) : @xmath394 which is a continuous function of @xmath11 . after differentiation with respect to the voltage the correction to the conductance is obtained as : @xmath395 and the second derivative of the current is : @xmath396 at positive voltages the second derivative of the current shows a negative dirac delta peak at @xmath397 , which reflects the energy spectrum of a single tls . above the excitation energy a positive background is seen due to the nonequilibrium distribution of the tls occupation number .    as a next step the elastic contribution is determined at @xmath62 for a single tls positioned in the middle of the contact . the change in the conductance due to the elastic scattering can be calculated by differentiating eq . ( [ iel ] ) . the elastic current correction contains a linear term , which causes a constant , voltage independent reduction of the conductance . experimentally it is hard to separate this constant contribution , thus we calculate only the voltage dependent part by subtracting the zero bias conductance : @xmath398 the second derivative of the @xmath1 curve is @xmath399 again , a dirac delta peak reflects the spectrum of the tls , and a continuous background arises at @xmath344 due to the nonequilibrium distribution . contrary to the inelastic case , the dirac delta peak and the background have the same sign ; furthermore , in the elastic case the sign of the peak can either be positive or negative depending on the sign of @xmath400 . the amplitude @xmath400 is related to the universal conductance fluctuation @xmath401 ( see e.g.  ) .    for arbitrary position of the tls the occupation number @xmath231 is zero at @xmath402 and @xmath403 at @xmath404 . therefore , the total amplitudes for the change in the conductance in the elastic and inelastic case can be generally written as : @xmath405 where isotropic scattering was assumed , and the formulas ( [ scatteringstrength ] , [ kisw ] , [ kr ] , [ dgamma ] ) were used . the results are given in the unit of the universal conductance quantum , @xmath401 . according to eq . ( [ munu ] ) for a highly asymmetric tls , where the energy splitting @xmath5 is much larger than the transition term @xmath6 , the equations @xmath406 and @xmath407 hold . in this case the inelastic term is suppressed . in the opposite case , where @xmath408 the elastic term is suppressed as @xmath409 . it must be also noted that the inelastic term depends only on the matrix element @xmath43 , whereas the elastic term is influenced both by @xmath43 and @xmath45 . the inelastic term can be roughly estimated as @xmath410 using @xmath411 , @xmath412 , @xmath413 ( see ) . the amplitude sharply drops by moving the tls further from the orifice than its diameter because of the geometrical factor . the estimation of the elastic term is more difficult as @xmath45 is strongly model dependent .    in the following the effect of a tls is considered at finite temperature . the formulas for the conductance and the second derivative of the @xmath1 curve can be explicitly calculated using equations ( [ iinel ] , [ cd ] , [ diel ] , [ ab ] , [ nplus ] ) ; however , these equations are very complicated , thus the results are demonstrated by figures . figure [ ball_cent ] shows the voltage dependence of the conductance for the inelastic and elastic case respectively . in fig . [ secder_ball ] the second derivative of the @xmath1 curve , the so - called point contact spectrum is presented at various temperatures both for the elastic and inelastic case . in fig . [ secder_ballq ] both the inelastic and elastic contributions are compared for a tls positioned in the contact center ( @xmath391 ) and for a tls being farther away ( @xmath414 ) .     and @xmath415 , respectively ) . the temperature is @xmath416 . the sign of the elastic contribution can be either positive or negative depending on the sign of ( @xmath400 ) . on the other hand , the inelastic scattering always causes a decrease in the conductance.__,title="fig:",scaledwidth=50.0% ] and @xmath415 , respectively ) . the temperature is @xmath416 . the sign of the elastic contribution can be either positive or negative depending on the sign of ( @xmath400 ) . on the other hand , the inelastic scattering always causes a decrease in the conductance.__,title="fig:",scaledwidth=50.0% ]     curve for inelastic and elastic scattering at various temperatures for a ballistic point contact ( @xmath417 , and @xmath418 , respectively ) . the inelastic contribution shows a negative peak at the excitation energy of the tls and a positive background tail at higher voltages . the elastic contributions similarly shows a peak at @xmath397 and a background tail at higher voltages , but here the sign of the peak and the tail is the same . in the elastic case the sign of the signal can either be positive or negative depending on the sign of ( @xmath400).__,title="fig:",scaledwidth=50.0% ] curve for inelastic and elastic scattering at various temperatures for a ballistic point contact ( @xmath417 , and @xmath418 , respectively ) . the inelastic contribution shows a negative peak at the excitation energy of the tls and a positive background tail at higher voltages . the elastic contributions similarly shows a peak at @xmath397 and a background tail at higher voltages , but here the sign of the peak and the tail is the same . in the elastic case the sign of the signal can either be positive or negative depending on the sign of ( @xmath400).__,title="fig:",scaledwidth=50.0% ]     curve for inelastic and elastic scattering for two different positions of the tls : in one case the tls is in the contact center , in the other case it is farther away on the contact axis , so that the opening is seen in a solid angle @xmath419 . the tls being farther away from the contact is closer to thermal equilibrium , which is reflected by the reduced amplitude of the background tail with respect to the spectroscopic peak in the inelastic contribution.__,title="fig:",scaledwidth=50.0% ] curve for inelastic and elastic scattering for two different positions of the tls : in one case the tls is in the contact center , in the other case it is farther away on the contact axis , so that the opening is seen in a solid angle @xmath419 . the tls being farther away from the contact is closer to thermal equilibrium , which is reflected by the reduced amplitude of the background tail with respect to the spectroscopic peak in the inelastic contribution.__,title="fig:",scaledwidth=50.0% ]      as it has been discussed in sec . [ maxwell ] the phenomenological treatment must be replaced by a theory based on kinetic equation , where effect of the inelastic scatterings are taken into account in the distribution function . in the following the elastic and inelastic scatterings are treated in equal footing . as a first step the distribution function is determined in the presence of elastic scattering and next the contributions of the tls are treated as a perturbation . consider first the elastic scattering in the limit @xmath420 , where the resistance is very large . because of the very strong elastic scattering , the electrons are immediately redistributed concerning the direction of their momenta . therefore , in the limit @xmath421 the distribution function depends only on the energy of the electrons , @xmath422 . as the electron arriving at the contact is either coming from the left or the right reservoirs and the electron energies are changed due to the external potential @xmath254 the distribution function must be the superposition of the distribution of the electrons coming from the left or right with amplitudes @xmath423 and @xmath424 , respectively : f_()=_0()f_0(_-e ( ) + ev2)+(1-_0())f_0(_- e()-ev2 ) . [ diffdistr ] using the charge neutrality condition ( @xmath425 @xmath423 can be determined by a similar treatment used earlier in case of eq . ( [ potent_op ] ) : _ 0()=12+()v . [ alpha ] this formula is valid in case of arbitrary geometry because only the shape of the potential function contains information about the details of the geometry . for an opening type circular contact the potential was determined in sec . [ maxwell ] by solving the laplace equation using a hyperbolic coordinate system . = 0.6 0.7truecm    the distribution function contains two sharp steps due to the two different fermi energies in the two reservoirs ( see fig . [ distrmax ] ) . such steps are measured by the saclay group studying short metallic wires of type fig . [ pcmodels.fig]b by attaching an extra tunneling diode on the side of the wire . ( see ) .    in the ballistic limit the distribution function is very similar only the factors @xmath423 and @xmath424 are replaced by the geometrical factors @xmath426 and @xmath427 determining the solid angles in which the ballistic electron arrives from the left or right reservoirs to the point @xmath164 . the electron distribution function ( [ diffdistr ] ) obtained for @xmath421 does not result in any current in agreement that the resistivity is infinite in this limit . for finite @xmath428 the factor @xmath423 must be replaced by a momentum dependent one @xmath429 , and the current is due to @xmath430 , which has a strong asymmetric momentum dependence and @xmath431 as @xmath421 . the @xmath432 will be determined by the very elegant theory of kulik , shekter and shkorbatov,@xcite who have extended the treatment of the ballistic limit to the diffusive one treating the scattering on phonons and defects which are large in space . the strong elastic scattering on defects is combined with a weak inelastic scattering . the limit considered is specified as @xmath433 where @xmath8 is the size of the contact and @xmath117 and @xmath118 are the elastic and inelastic mean free paths of the electrons , respectively . introducing the corresponding relaxation times , @xmath434 and @xmath435 the inequality can be rewritten as @xmath436 . that means that an electron can diffuse out from the contact region of size @xmath8 with small probability of inelastic scattering , thus double inelastic scatterings can be neglected . in the following @xmath435 is related to the relaxation time due to tlss . also in the diffusive limit the current correction due to the tlss is derived by solving the boltzman equation ( [ boltzmann_new ] ) with an appropriate collision term . in this case the collision integral @xmath209 has two contributions as the elastic impurity part @xmath437 and the contribution of individual tlss at position @xmath244 , @xmath438 thus @xmath439 where the sum is due to different tlss , but in the following it is assumed that the concentration of the tlss is so low that first a single tls at position @xmath244 is considered and in the final result the contributions of the different tlss are additive . the kinetic equation can be arranged as @xmath440 first the left hand side is treated in the diffusive limit and the right hand side is considered as a weak perturbation . the impurity part of the collision integral is @xmath441 where @xmath442 is the elastic transition probability , @xmath443 is the velocity at momentum @xmath216 perpendicular to the equienergetic surface for which the integral is performed . the electric field can be expressed by the electric potential @xmath254 as @xmath444 , where @xmath125 is determined by the neutrality condition as in the ballistic case . the distribution function @xmath445 satisfies the following boundary condition : at the boundary , @xmath446 between the insulator and the metal we assume mirror reflection @xmath447 where an incoming electron with momentum @xmath168 is reflected with the momentum @xmath448 . furthermore , very far from the contact @xmath313 the equilibrium distribution is recovered with chemical potential @xmath22 , thus @xmath449 and @xmath22 satisfies the charge neutrality . the distribution function in the absence of the collision term due to tls is denoted by @xmath450 and the first order correction due to the tls is @xmath451 , thus @xmath452 where the higher order corrections are neglected . the kinetic equation for @xmath450 contains the static impurity contributions , thus @xmath453 and @xmath451 satisfies the equation @xmath454 where the electric field is also expanded as @xmath455 . the term @xmath456 is linearized in the collision , thus @xmath445 is replaced by @xmath450 . similarly , the current is @xmath457 , where @xmath458 and the integral with respect @xmath459 is taken on a dividing surface representing the point contact .    according to the introductory remarks of this section , at finite elastic mean free path the distribution function can be expressed with a momentum dependent parameter , @xmath429 : @xmath460 where the electrons arriving from far left ( right ) at the contact have distribution function @xmath461 ( @xmath462 ) , thus @xmath463 as the collision term ( [ eq : i_i ] ) is linear in the distribution function the equation for @xmath450 is satisfied with @xmath429 as well : @xmath464 in the @xmath421 limit the solution of this equation is @xmath423 , given by ( [ alpha ] ) . at finite elastic mean free path @xmath429 can be expanded as : @xmath465 where @xmath432 is the first order term in the small parameter , @xmath466 . the momentum dependent correction , @xmath467 can be determined by a simple argumentation . at non - zero relaxation time the probability @xmath468 can be considered as the momentum independent probability @xmath469 taken at the position of the last elastic collision , that is : @xmath470 the expansion of this formula in @xmath434 gives : @xmath471 the value of @xmath469 is determined by the potential , which drops in the contact region , thus @xmath472 , and @xmath473 . the value of @xmath467 can also be obtained by inserting ( [ deltaalpha ] ) into the boltzmann equation for @xmath468 ( [ eq : alpha_kineq ] ) . after neglecting the higher order terms in the small parameters @xmath474 and @xmath475 the following simple formula is achieved : @xmath476 for isotropic scattering @xmath477 holds , the relaxation time approximation is appropriate , thus the collision integral is expressed as : @xmath478 which provides the relaxation time , @xmath434 due to the elastic scattering by the impurities . inserting this result into ( [ boltzalpha ] ) the formula ( [ deltaalpha2 ] ) for @xmath467 is regained . the next step is to determine @xmath479 using the equation ( [ boltzalpha ] ) for @xmath468 . the boundary condition for @xmath480 as very left the electron comes from the left contact where the electrons are in equilibrium thus @xmath481 furthermore , the collision term contributes to the distribution function only if the system is out of equilibrium , thus @xmath479 is proportional to the electric field . therefore , the second term @xmath482 in equation ( [ eq : kineq_sep ] ) can be dropped as it is proportional to @xmath483 . the correction to the kinetic equation is determined by eq . ( [ eq : kineq_sep ] ) where after linearizing each term in the field the ansatz with the function @xmath484 @xmath485 can be used . the terms proportional to @xmath486 cancel out in eq . ( [ eq : kineq_sep ] ) . thus @xmath487 with the boundary condition @xmath488 and @xmath489 . as it will be shown below the change in the distribution due to the collision term containing the interaction with the tls can be expressed as linear response with the aid of a green s function @xmath490 : @xmath491 where @xmath492 and @xmath493 with the boundary condition @xmath494 and @xmath495 . as @xmath490 is closely related to the correlation function @xmath496 the time reversal symmetry can be applied as @xmath497    the function @xmath498 is defined as @xmath499 where the integral with respect to @xmath500 is a surface integral of variable @xmath501 taken on the orifice . this function satisfies the equation @xmath502 which follows from the equation ( [ chipr2 ] ) defining @xmath490 . the boundary conditions are the usual ones @xmath503 and @xmath504 . this equation for @xmath498 is closely related to eq . ( [ boltzalpha ] ) for @xmath429 introduced earlier , but without the inhomogeneous term . ( the sign in front of the collision integral can be changed by replacing @xmath505 and taking the relation @xmath506 into account . ) by comparing the two equations the expression @xmath507 is obtained , which automatically satisfies the boundary condition for @xmath498 .    using the expression ( [ eq : galpha ] ) for @xmath498 the current can be calculated and only the @xmath168-dependent part is contributing , as the total number of electrons at position @xmath164 is not changed by the collision . therefore : @xmath508    the kinetic equation ( [ eq : chi_kineq ] ) for @xmath509 can be checked using equations ( [ chipr1 ] ) , ( [ chipr2 ] ) , ( [ chipr3 ] ) and inserting twice the identity ( [ gpp ] ) . the next task is to calculate the correction to the current due to the inelastic scattering on the tls . using the expression ( [ chi ] ) for @xmath479 ( dropping the first term which is even in momentum @xmath168 in leading order ) : @xmath510 where the variable @xmath164 is taken over the contact area , @xmath511 . inserting the expression of @xmath512 the result can be rewritten as @xmath513 that result can be further rewritten in terms of @xmath498 as @xmath514 using eq . ( [ eq : galpha ] ) for @xmath498 the momentum independent terms do not contribute due to eq . ( [ qwerty ] ) . the remaining term contains @xmath467 which was determined earlier and @xmath515 holds . furthermore , @xmath516 and in this way the final expression for the current is @xmath517 where @xmath467 is taken at the position of the tls , @xmath244 . because of the odd parity of @xmath467 only the odd part of the collision integral contributes . using the equation ( [ falpha ] ) for @xmath450 the collision integral for the inelastic scattering can be written similarly to the ballistic case ( eq .  [ irtlsin ] ) : @xmath518 \left[\alpha_{\hat{\z{p}}}(1-f_0(\varepsilon+{ev\over2 } ) ) + ( 1-\alpha_{\hat{\z{p}}})(1-f_0(\varepsilon-{ev\over2}))\right ] \right.\cdot\nonumber\\ & { } & \;\;\;\bigg[n_-\delta(\varepsilon'-\varepsilon - e)+ n_+\delta(\varepsilon'-\varepsilon+e)\bigg]\nonumber\\ & { } & \;-\bigg[\quad\hat{\z{p}}'\leftrightarrow\hat{\z{p}}\quad \varepsilon'\leftrightarrow\varepsilon\quad\bigg ] \bigg[\quad\hat{\z{p}}'\leftrightarrow\hat{\z{p}}\quad \varepsilon'\leftrightarrow\varepsilon\quad\bigg ] \bigg[\quad \varepsilon'\leftrightarrow\varepsilon\quad\bigg ] \bigg\}\end{aligned}\ ] ] is obtained , where in the second set of brackets the variables are changed as @xmath317 , @xmath519 and all the energy variables @xmath520 and @xmath330 are shifted by @xmath521 thus @xmath254 drops out . now using eq . ( [ deltaalpha ] ) @xmath522 and keeping only the odd part of the collision integral , after changing the variables the current is @xmath523 \left[\alpha_0(1-f_0(\varepsilon+{ev\over2}))+ ( 1-\alpha_0)(1-f_0(\varepsilon-{ev\over2}))\right ] \right.\right.\nonumber\\ & { } & \;\;\;-\delta\alpha_{\hat{\z{p } } } \left[\alpha_0f_0(\varepsilon'+{ev\over2}))+ ( 1-\alpha_0)f_0(\varepsilon'-{ev\over2})\right ] \left[f_0(\varepsilon+{ev\over2})-f_0(\varepsilon-{ev\over2})\right ] \bigg]\times\nonumber\\ & { } & \;\;\;\left[n_-\delta(\varepsilon'-\varepsilon - e)+ n_+\delta(\varepsilon'-\varepsilon+e)\right]\nonumber\\ & { } & \;\;-\bigg[\quad\hat{\z{p}}'\leftrightarrow\hat{\z{p}}\quad \varepsilon'\leftrightarrow\varepsilon\quad\bigg]\bigg\},\end{aligned}\ ] ] where @xmath469 and @xmath467 are taken at the site of the tls ( @xmath308 ) .    in this expression two squared averaged matrix elements occur @xmath524 where the expression ( [ deltaalpha2 ] ) for @xmath467 was used . furthermore @xmath525 holds as @xmath526 and @xmath527 integrating over the variable @xmath330 the following expression is obtained @xmath528\right\},\nonumber\end{aligned}\ ] ] where the geometrical factor @xmath529 is introduced for the diffusive limit as @xmath530 the expression ( [ idiff ] ) can be brought to a form which is comparable to the result ( [ eq:112 ] ) in the ballistic limit . the key of the algebra applied is in the integral @xmath531 the variable can be shifted by an arbitrary energy @xmath532 as the contribution comes from the neighborhood of the fermi energy and the bandwidth cutoff does not play any role . the result of a lengthy algebra is i^(1)=-wk_^diff f_0()\{(f_0(-e - ev)- f_0(-e+ev))n_-+(f_0(+e - ev)- f_0(+e+ev))n_+}.this expression coincides with the one obtained in the ballistic limit ( [ eq:112 ] ) after replacing the geometrical factor @xmath533 . therefore , the final results in the ballistic limit ( [ iinel ] , [ cd ] , [ iinel2 ] ) are also valid in the diffusive limit , just the geometrical factor , @xmath534 is different . the elastic contribution can be obtained by inserting @xmath245 and considering different scattering strengths for the two states of the tls in the diffusive limit as well . therefore , the ballistic result for the elastic case ( [ diel ] ) is also valid in the diffusive regime with the replaced geometrical factor . the next step is to calculate @xmath231 . in the present case whether the electrons come from the left or right reservoir is determined by @xmath469 or @xmath535 , thus the previous results in the ballistic regime ( [ nplus ] , [ ab ] , [ nt0 ] ) can be adopted after changing @xmath536 the momentum dependent correction @xmath467 to @xmath469 contributes only if the momentum dependence of @xmath537 is kept but that is small in case @xmath538 . concluding , in the diffusive limit the shape of the current correction due to the scattering on a tls is the same as in the ballistic case , the only essential difference is in the geometrical factors , which are independent of the strength of the interaction between the electrons and the tlss . in the ballistic case @xmath534 is in the range of unity if the tls is situated in the contact region . however , @xmath529 contains other quantities proportional to @xmath539 which indicates how much @xmath16 is changed in space on the scale of the impurity scattering mean free path . we have seen that @xmath540 scales with the inverse of the contact diameter , @xmath8 , thus @xmath541 in case where the tls is nearby the orifice . far from the orifice ( @xmath313 ) @xmath542 . the amplitudes of the change in the conductance for the important tls ( being close to the contact ) are scaled down by @xmath543 which can be a very small factor in the dirty limit . the largest contribution arises from the tlss for which @xmath544 is the largest , thus for the tls just nearby the edge of the orifice . ( see eq . [ ez ] )    the change of the geometrical factor due to the impurity scattering length is similar to the case of electron - phonon interaction , where integration is taken according to location of the interaction.@xcite the @xmath53-factors for electron - phonon scattering are given in ( [ kfactorball ] , [ kfactordiff ] ) for the ballistic and diffusive case , respectively . in the phonon case the amplitude of the signal in the diffusive limit is reduced only by @xmath545 .      up to this point the scattering on tlss in point contacts was investigated as a function of bias voltage , which provides information about the excitation spectrum of the tls , and also about the nonequilibrium distribution of the tls states . a further possibility to study the dynamical properties of the tlss is the study of the response to a harmonic excitation . let us consider that a voltage of @xmath546 is applied on the point contact , where the dc bias voltage , @xmath547 is varied on the scale of the excitation energy of the tls to cover the spectroscopic peak at @xmath548 , and also the background tail at @xmath344 . the amplitude of the harmonic excitation is considered to be much smaller than the excitation energy , @xmath549 . furthermore , it is assumed that the energy @xmath550 is much smaller than the excitation energy of the tls , thus the quantum nature of the radiation should not be taken into account . in this case the current flowing through the contact can be expanded in second order in terms of the time dependent voltage contribution : @xmath551 where the equation @xmath552 is used . the voltage dependence of the harmonic response , @xmath553 provides the differential resistance curve , @xmath554 . the second harmonic response , @xmath555 directly provides the point contact spectrum , i.e.  the second derivative of the @xmath1 curve . the dc component , which can be measured as a time average of the response signal : @xmath556 also shows the second derivative of the @xmath1 curve after subtracting the current without harmonic excitation , @xmath557 . further on this term is referred to as the dc shift signal , denoted by @xmath558 . the above considerations are , however , only valid if the system is fast enough to follow the time dependent voltage . in the case of two level systems , the dynamical behavior is determined by the kinetic equation ( [ dndt ] ) : @xmath559 where the coefficients @xmath189 and @xmath72 are given by eq . ( [ ab ] ) . both of the coefficients depend on the actual voltage @xmath560 at given @xmath561 . for simplicity the notation @xmath562 is introduced , which is the inverse of the relaxation rate of the tls ( @xmath563 ) . for slow two level systems the value of the relaxation time can be estimated using eq . ( [ dndt],[kisw],[scatteringstrength ] ) as follows : |b ( 2)^2(_0 v^z)^2[relaxtime]inserting @xmath412 , @xmath564mev , and @xmath413@xcite into this equation a relaxation time of @xmath565s is obtained . if the alternating voltage is faster than the relaxation of the tls ( @xmath566 ) then a more complicated procedure is necessary to determine the time dependent response of the system . this phenomenon , however can be used to trace the typical relaxation time of the tls by studying the response signal as a function of the frequency @xmath68 . the first calculation with time dependent applied voltage was performed by kozub and kulik.@xcite in the next part a general solution is given for the time dependent occupation number @xmath567 at arbitrary frequency and amplitude of the harmonic excitation signal . after that various limits are treated , where the formulas are essentially simplified . the kinetic equation ( [ mw1 ] ) is a linear inhomogeneous differential equation which has a general solution @xmath568 where the coefficient @xmath569 is determined by the initial condition at time @xmath570 . according to the definitions given by eq . ( [ ab ] ) @xmath571 and @xmath572 hold , thus the exponent in the last term of eq . ( [ mw2 ] ) monotonically increases with time . this means that @xmath569 is essential only for the transient behavior and for long time it can be dropped . for periodic voltage @xmath189 and @xmath573 are also a periodic function of time : @xmath574 , @xmath575 . the time integrals can be divided according to the periods as : @xmath576 where @xmath577 , @xmath578 is an integer and @xmath579 . in this way : @xmath580 where @xmath581 and @xmath582 using these identities , one gets : @xmath583 the final result is : @xmath584 which gives the following stationary solution in the long time @xmath585 limit : @xmath586 turning to dimensionless variables like @xmath587 , and @xmath588 and also new notations , @xmath589 , @xmath590 , @xmath591 and @xmath592 one gets : @xmath593 it is important to note that @xmath594 is proportional to the relaxation rate @xmath595 and @xmath596 to the excitation rate of the upper level , which are independent of the frequency @xmath68 .    in the following two simple limits are treated : @xmath597 ( @xmath598 ) ; and @xmath599 ( @xmath600 ) . [ [ low - frequency - limit - omegataull-1 ] ] low frequency limit ( @xmath601 ) + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    assuming that @xmath602 ( @xmath603 ) the first part of the right hand side of eq . ( [ mw10 ] ) disappears as @xmath604 . ( note that @xmath594 is a monotonous function . ) the contribution to the second part also disappears as far as @xmath605 is finite . in the vicinity of the upper limit of the integral ( @xmath606 ) , however , @xmath607 can be expanded as @xmath608 , where @xmath609 , thus @xmath610 where the term arising from the lower limit is dropped as it is exponentionally decreasing with @xmath611 . the final result for the limit @xmath611 is : @xmath612 that result also holds for @xmath613 , where the condition @xmath614 is not satisfied any more . without making a rigorous proof , it can be argued , that @xmath613 is at the beginning of the time period @xmath615 . in the long time limit ( @xmath616 ) the time can be shifted , therefore , that time is not special , thus the result ( [ mw11b ] ) must be generally valid .    in this way the static result is recovered with @xmath189 and @xmath617 taken with the actual voltage at time @xmath561 . therefore , at arbitrary time the system is in the stationary state , and the the expansion of the current in ( [ iexpansion ] ) is valid . [ [ high - frequency - limit - omegataugg-1 ] ] high frequency limit ( @xmath566 ) + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    taking @xmath618 only the first part of the right hand side of eq . ( [ mw10 ] ) contributes . the denominator can be expanded as @xmath619 and the result is @xmath620 where @xmath621 stands for an average for a time period and the definition of @xmath622 is used . as the tls reacts slowly on that time scale , thus it reacts on the average value of the parameters @xmath189 and @xmath617 . note that @xmath623 . ( a)(b)(c )    in order to calculate the time dependent current first the occupation number @xmath567 must be calculated , which is not an easy task in the general case as the formula given by eq . ( [ mw10 ] ) contains double integrals . the next step is to insert the obtained @xmath567 into the equation for the current ( [ iinel ] ) . if one is interested in the average current ( dc ) and the harmonics of the current , it is not enough just to determine the the averaged @xmath231 and its harmonics . that can be demonstrated easily at @xmath62 by eq . ( [ iinel2 ] ) . at @xmath624 the expression contains a term , which is independent of the occupation number . this causes a nonlinearity in the @xmath1 curve in every case , regardless of the voltage dependence of the occupation number . for instance , if the tls was in thermal equilibrium with the bath then @xmath340 would hold at any bias voltage . in this case the @xmath1 curve is linear both at @xmath625 and at @xmath344 , but it has a breaking point at @xmath548 . if the alternating voltage is smaller than @xmath40 in the whole time period ( @xmath626 ) , then the response to the excitation is completely linear , and no higher harmonics are generated . this situation is demonstrated in fig . [ rectification1.fig]a . similarly , the response is linear if @xmath627 at any time . however , if the condition @xmath626 is only valid in a part of the time period , the system gives a rectified response to the harmonic excitation , thus higher harmonics are also generated . this case is shown in fig . [ rectification1.fig]b . the tls can be excited only in the part of the time that is demonstrated by the shadowed areas , where the conductance is reduced due to the possibility of backscattering . in fig . [ rectification1.fig]c a third case is presented , where the amplitude of the alternating voltage is so large , that rectification occurs both at @xmath397 and at @xmath628 . in this case the contributions from the both polarities work against each other , and the higher harmonic generation is reduced .    in a general case both the nonlinearity of the occupation number @xmath347 and the above discussed intrinsic nonlinearity of the @xmath1 curve cause higher harmonic generation . however , if the period time of the excitation is much shorter than the characteristic relaxation time of the tls then the occupation number can not follow the alternating voltage , thus it sets at a constant steady state value . in this situation the higher harmonic generation is still present due to the intrinsic nonlinearity of the @xmath1 curve .    in the following section these features are discussed in two special cases . in sec . [ rectification.sec ] the @xmath62 the limit is treated in the high frequency limit ( @xmath629 ) , where the rectified current and the harmonics can be easily calculated at arbitrary value of the amplitude of the excitation , @xmath630 . these results are also valid at finite temperature if @xmath631 . in the second case ( sec . [ expansion.sec ] ) an expansion with respect to @xmath630 is applied . this approach is valid if the amplitude of the excitation is small ( @xmath632 ) , however in this limit the frequency dependence can be investigated in a wide range of frequencies . the rectification due to a tls in the contact can be nicely demonstrated at @xmath62 in the @xmath629 limit . we have seen that at @xmath629 the occupation number @xmath231 becomes time independent which makes it possible to perform the calculation analytically . according to eq . [ mw11c ] @xmath633 , where @xmath634 and @xmath635 are the time average of these coefficients . inserting zero temperature into the formula ( [ ab ] ) the voltage dependence of @xmath189 and @xmath573 can be written as : @xmath636 therefore , in case of amplitudes for which @xmath637 holds for the total period of time the occupation number is zero . if that is not satisfied , then the actual voltage at a certain part of the time can excite the tls out of the groundstate . two cases must be distinguished : ( i ) where @xmath638 but @xmath639 ; ( ii ) where @xmath638 and @xmath640 . case ( i)case ( ii )    figure [ rectification2.fig ] show the actual voltage as a function of time using @xmath641 in these two cases . the shadowed areas indicate the regions where @xmath642 . the phases , where the shadowed areas are ending are denoted by @xmath643 ( in case ( i ) ) and @xmath644 ( in case ( ii ) ) . all of the following results for case ( i ) can be applied also for case ( ii ) by inserting @xmath645 . these phases can be expressed by @xmath547 , @xmath630 and @xmath40 as follows :    @xmath646    the time average of @xmath189 is calculated as : @xmath647{\rm d}\varphi+ \int\limits_{\bar\varphi_0}^{\pi}\left[-e(v_0+v_1\cos\varphi)-e\right]{\rm d}\varphi\right\}=\\ = \frac{\kappa}{\pi}e\left [ v_0(\varphi_0-\varphi_0^{\star})+v_1(\sin\varphi_0+\sin\varphi_0^{\star})\right]-\frac{\kappa}{\pi}e(\varphi_0+\varphi_0^{\star}),\end{aligned}\ ] ] and the occupation number is written as @xmath648 after having the result for @xmath231 the inelastic current correction can be calculated using eq . [ iinel2 ] . at the points , where @xmath649 the time derivative of the current is discontinuous , which results in higher harmonic generation . the second harmonic generation with @xmath650 and @xmath651 and also the shift of the dc signal can be easily calculated : @xmath652 the results for the @xmath650 are plotted for different values of the excitation amplitude , @xmath630 . in this case the width of the peak at @xmath653 is determined by the amplitude @xmath630 , and these results are also valid at finite temperature if @xmath654 . the curves are plotted for different amplitudes of the excitation signal . _ _ ]      at finite frequencies the occupation number is already time dependent ( @xmath567 ) , thus the above considerations can not be applied . a calculation for intermediate frequencies can be performed in the @xmath655 limit . in this case the integrals in the general solution for @xmath231 ( eq .  [ mw2 ] ) are calculated after expanding the coefficients @xmath656 and @xmath657 with respect to @xmath630 . the expansion for @xmath189 is : @xmath658 and similarly for @xmath617 . the derivatives are taken at @xmath659 . in the general solution for @xmath231 first the quantity @xmath660 must be calculated . the non - oscillating part of the exponent is linearly increasing with time , thus the expansion can be carried out only for the oscillating part assuming that @xmath661 and @xmath662 : @xmath663 where @xmath664 the shift of the relaxation rate by a term proportional to @xmath665 is due to the nonlinear behavior . the next step is to calculate the integrals @xmath666 . all the terms obtained are proportional to the exponential , @xmath667 , which is cancelled out by the factor @xmath668 in the expression ( [ mw2 ] ) . in the remaining part the higher order terms @xmath669 are dropped . the next tedious task is to collect all the terms in the coefficients of the expansion . the occupation number can be written in the following form : @xmath670 where @xmath671 is the stationary occupation number at @xmath672 ( @xmath673 ) , @xmath674 and @xmath675 define the amplitudes of the harmonic response to the excitation , @xmath676 is the shift of the stationary value due to the nonlinear behavior , while @xmath677 and @xmath678 define the amplitudes of the second harmonic generation . the values of these coefficients are given by the following equations : @xmath679 the conditions for the expansion have excluded the narrow frequency region around @xmath680 . accordingly , the above results contain singular terms in @xmath68 , which diverges as @xmath681 . on the other hand , it has been seen that the @xmath681 limit is well defined for @xmath231 , and indeed , it can be shown that the singular terms cancel out in the low frequency limit .    ) and the second harmonic terms ( @xmath682 , @xmath683 ) as a function of bias voltage at different excitation frequencies in the inelastic case . the right panels presents the frequency dependence of the corresponding terms at a fixed voltage of @xmath397 . the temperature is @xmath416.__,scaledwidth=60.0% ]    the expression for the current correction can be expressed similarly to ( [ mw20 ] ) : @xmath684 in the inelastic case the current correction given by eq . ( [ iinel ] ) can be written as @xmath685 where @xmath686 . these must be expanded similarly to equation ( [ mw17 ] ) given for @xmath189 . after this step the coefficients in ( [ mw23 ] ) can be given as follows : @xmath687 these results can be analytically calculated , and for demonstration the experimentally important quantities ( @xmath688 , @xmath689 , and @xmath690 ) are plotted in fig . [ mwin.fig ] . the left panels show the voltage dependence at different @xmath68 values , while the right panels demonstrate the @xmath68 dependence at the fixed bias voltage value of @xmath397 . at low frequencies , the tls occupation number can follow the alternating voltage , thus both @xmath690 and @xmath558 are proportional to the second derivative of the @xmath1 curve ( see eq . [ iexpansion ] ) . in this limit there is no phase shift , and thus @xmath691 . in accordance with the static results , the curves show a negative peak at @xmath397 , and a positive tail at higher voltages . in the high frequency limit the tls can not follow the excitation , thus the occupation number sets at a constant value . due to the intrinsic nonlinearity of the @xmath1 curve , even in the high frequency limit a pronounced negative peak is observed at @xmath397 , but the positive background disappears . ( at zero temperature the amplitude of the negative dirac delta peak at @xmath397 in @xmath690 is exactly @xmath692 times smaller in the @xmath693 limit than in the @xmath681 limit . this factor of @xmath692 is well demonstrated by the @xmath68 dependence of @xmath694 in fig . [ mwin.fig ] as well . ) around the characteristic frequency of the tls ( @xmath695 ) a significant phase shift is observed , as demonstrated by the curves for @xmath689 .    ) and the second harmonic terms ( @xmath696 , @xmath697 ) as a function of bias voltage at different excitation frequencies in the elastic case . the right panels presents the frequency dependence of the corresponding terms at a fixed voltage of @xmath397 . the temperature is @xmath416.__,scaledwidth=60.0% ]    in the elastic case the different current components can be expressed by using eq . ( [ iel ] ) : @xmath698 in fig . [ mwel.fig ] the second order terms are plotted for the elastic case . the nonlinearity of the @xmath1 curve of the elastic contribution is basically determined by the nonlinearity of @xmath347 . above the characteristic frequency of the tls the occupation number can not follow the excitation , thus both @xmath699 and @xmath699 are suppressed . the shift of the dc component ( @xmath700 ) shows the spectroscopic peak at @xmath397 even at @xmath701 , just the amplitude is reduced by @xmath702 . around the characteristic frequency of the tls a significant @xmath699 contribution is observed . the dynamical properties of the system are also demonstrated by the phaseshift in the second harmonic contributions : @xmath703 . in fig . [ phaseshift.fig ] the phaseshift @xmath704 is plotted at @xmath397 both for the inelastic and elastic case . in the elastic contribution the phaseshift grows from @xmath705 to @xmath706 between the @xmath707 and the @xmath701 limit . contrary , in the inelastic case the phaseshift is zero in both limits , and even at @xmath708 a only minor phaseshift of @xmath709 is observed , as the nonlinearity is basically related to the intrinsic nonlinearity of @xmath1 curve and the nonlinearity of @xmath347 plays a less important role .     for the elastic and the inelastic contribution , respectively . the temperature is @xmath416.__,scaledwidth=35.0% ]      in this section the theory of the scattering on slow tlss in point contacts was reviewed . in sec . [ ballistic ] it was seen that the electron distribution functions and the electric potential around a ballistic contact can be determined by solving the homogeneous boltzmann equation for a point contact geometry . the solution yields a voltage independent conductance , which is known as the `` sharvin conductance '' ( eq . [ sharvinconductance ] ) . the effect of the scattering on a slow tls near the contact can be considered by perturbation in the collision term of the homogeneous boltzmann equation . the electron - tls interaction may result in backscattering processes , where an electron that has already crossed the contact is scattered back through it . these processes give a strongly voltage dependent correction to the current , for instance at zero temperature backscattering can only occur if the voltage bias exceeds the excitation energy of the tls . the nonlinear current correction due to the inelastic scattering of the electrons on a single tls can be explicitly given , if the occupation number of the upper level @xmath231 , the electron - tls coupling strength @xmath373 , and the geometrical factor @xmath710 are known . it has crucial importance that the interaction of a tls with the neighborhood is dominated by the electron - tls scattering , and the phonon contribution is significantly smaller because of the relatively small phonon phase space at energies much smaller than the debye energy . therefore , the excited tls mainly relaxes its energy through the electrons , which have a nonequilibrium distribution near the contact . as a consequence , the tls occupation number is also not the equilibrium one , and it is strongly voltage dependent . the value of @xmath347 can be determined by a simple dynamical equation using fermi s golden rule . due to voltage dependence of @xmath231 a nonlinear current correction arises from the elastic scattering of the electrons on the tls as well , provided that the scattering cross sections are different for the two states of the tls . both for the elastic and inelastic scattering the point contact spectrum ( i.e.  the second derivative of the @xmath711 curve ) shows a sharp peak at the excitation energy of the tls , and a background tail at higher voltages due to the nonequilibrium occupations . the current correction due to the scattering on a single tls can also be determined for a diffusive point contact . in this case the equations are more difficult , as the collision integral due to the elastic scattering on defects should also be included in the boltzmann equation . surprisingly , the lengthy calculations give exactly the same results as in ballistic contacts , just the geometrical prefactors are different . the response to a time dependent bias voltage was also reviewed . due to the nonlinearity of the @xmath711 curve the application of a small ac modulation beside the dc bias causes two important features in the response signal : ( i ) a second harmonic signal is generated , ( ii ) the average value of the current is shifted , which is referred to as the dc shift signal . if the frequency of the ac modulation is much smaller than the relaxation frequency of the tls and its amplitude is small compared to the energy splitting , then both of these signals are proportional to the second derivative of the @xmath711 curve . on the other hand , if the modulating signal is fast compared to the relaxation of the tls then the occupation number @xmath231 can not follow the variation of the bias voltage , which has important consequences . previously it was believed that in the high frequency limit both the elastic and inelastic contributions from the scattering on the tls are completely suppressed in the response signal . our more detailed calculations have shown that the dc shift signal contains a frequency independent contribution due to the rectification at voltages in the range of the tls energies thus it is not decisive for the relaxation dynamics of the tls . we have demonstrated that both the elastic and inelastic contributions are almost independent of the irradiation frequency in the dc shift signal . a more appropriate quantity for these studies is the second harmonic generation signal . in this case the elastic contribution is completely suppressed and the inelastic contribution shrinks by a factor of @xmath692 if the irradiation frequency exceeds the relaxation frequency of the tls . we have also investigated the phase shift in the second harmonic signal which might be examined in subsequent experimental studies . there is another suggestion by kozub and rudin@xcite to explain zero bias anomalies based on the adiabatic renormalization of the splitting of the tls by electrons . that gives , however , a weak temperature dependence because the main contribution comes from the electrons far from the fermi energy , thus that is not related to infrared singularities . the interaction between the tls due to the friedel oscillation was also considered and magnetic field dependence is suggested . however , it is hard to compare the theory with experiment in its present form . it is also interesting to note that tlss and magnetic impurities interacting with metallic electrons show several similar features . both could lead to logarithmic corrections known as the kondo effect . the main similarity is that the spin @xmath712 has also two states . without external magnetic field or internal magnetic field characteristic to the spin glass systems the two states @xmath713 are degenerate , just like a tls without splitting ( @xmath714 ) . in that case the occupations of the two spin states @xmath715 and that can not be affected by the applied voltage on the point contact because of spin degeneracy . in born approximation that gives a single peak at zero bias in the point contact spectra @xmath716 , and no background occurs . the logarithmic kondo corrections , however , broaden the peak which looks very similar to background in case of phonons and tls , but that is not due to a nonequilibrium distribution.@xcite the situation is different , where the spin degeneracy is lifted by the magnetic field . in that case @xmath717 , thus the occupation numbers depend on the applied voltage and background occurs even without kondo corrections .    in the following section the former experimental results are discussed and confronted with the model of slow two level systems . this review is aimed at collecting and summarizing the main results concerning the theory of the scattering on slow two - level systems in point contacts , and to compare those with the experiments . as a point contact probes the restricted area around the contact region , it can be used to investigate a few or even a single tls located in the contact area . therefore , some measurements can be directly used to study the nature of tlss in the contact region , and to check the relevance of the model presented .    in disordered metallic point contacts tlss with a wide variety of relaxation times can be present . a spectacular sign of the presence of tlss in a point contact is the telegraph fluctuation of resistance , caused by defects fluctuating between metastable configurations on an experimentally resolvable timescale.@xcite examples for such so - called two level fluctuators are shown in fig . [ buhrmantls ] . panel ( a ) shows the simplest telegraph fluctuation : the resistance is switching between two values , as discussed by considering the elastic processes ( see sec .  [ elasticscattering ] ) , while the switching itself is an inelastic process . in panel ( b ) two independent telegraph fluctuations are superimposed upon each other , while on panel ( c ) and ( d ) one two - level fluctuator modulates the amplitude or the average switching time of the other . the screening of these two - level systems by electrons causes a slow - down of the tls motion which can be described by the renormalization of the tunneling amplitude.@xcite this phenomenon can be used to estimate the electron - tls coupling constant , as it is discussed in and . this review concentrates on slow tlss where the relaxation time is much larger than the electron - tls interaction time , but still the tls is fast enough not to be able to resolve it as a telegraph fluctuation of the resistance . it was seen , that such `` slow '' tlss can be studied by point - contact spectroscopy , a singularity in the pc spectrum ( @xmath718 ) is expected if the voltage bias coincides with the energy splitting of the tls close to the contact . figures [ ball_cent ] and [ secder_ball ] demonstrated that the contributions of both the elastic and inelastic scattering exhibit this anomaly at @xmath397 . in the inelastic case this anomaly always corresponds to a maximum of the differential conductance , @xmath719 , which is reflected by a negative peak in the point contact spectrum @xmath720 at @xmath397 . in the elastic case a similar peak is observed , however its sign can either be negative or positive , depending on the sign of @xmath400 ( see eq . [ iel ] ) . a similar peak is expected in the pc spectrum for fast kondo - like tlss , for which the electron assisted transition processes are not negligible ( see sec .  [ fasttls ] ) . for the fast tlss the anomaly appears as a minimum in the differential conductance around zero bias , and the width of the zero bias minimum is determined by the kondo temperature @xmath0 . this feature is shown as a positive peak in the @xmath720 curve at @xmath721.@xcite that peak shows very strong similarity to the kondo peak in the presence of magnetic impurities,@xcite but no magnetic field dependence is observed , and the peak is changing by annealing . the sign of the anomaly has crucial importance in distinguishing the slow tlss from fast tlss . along this discussion the anomaly exhibiting a conductance maximum at zero bias is called `` positive '' zero bias anomaly ( zba ) , while the conductance minimum is called `` negative '' zba . we note that the labelling of the sign of the anomaly can be confusing , as in some cases the voltage dependence of either the differential conductance or the differential resistance is plotted , in other studies the point contact spectrum is presented . furthermore , concerning the pc spectrum , in theory the @xmath720 curves have direct physical meaning ( showing the excitation spectrum ) , whereas in experimental studies usually the second derivative of the inverse curve , @xmath720 is plotted . regardless of the way how the pc spectrum is plotted , the sign of the anomaly can be determined by comparing it with the sign of the contribution of the phonons . for the inelastic scattering on tlss the increase of the bias voltages causes an enhancement of the back - scattering , thus this anomaly has the same sign as the phonon part of the spectrum . in contrast , for kondo - like anomalies a resonant scattering is observed , which is restricted to the energy range of the kondo temperature . in this case the back - scattering is enhanced at zero bias , and it decreases as the bias is elevated above @xmath722 , thus the anomaly has negative sign compared to the phonon contribution . in terms of this labelling , the fast tlss should always show a negative anomaly , the inelastic contribution from slow tlss shows up as a positive anomaly , while the elastic contribution has an indefinite sign .    in the experiments different zero bias anomalies were observed . a basic distinction should be made between the experiments performed on pure metals , and those , performed on disordered systems , like metallic glasses . in the first case the temperature dependence of the bulk resistance does not show any anomalous behavior , but due to the magnifying effect of point contacts , the scattering on a few defects near the contact might be seen in the point contact spectrum . in this case the appearance of the anomaly depends on the way the contact is created , and usually the anomaly can be destroyed by annealing . in metallic glasses , however , a high degree of disorder is initially present , and the anomalous feature seen in the pc spectrum is also reflected by the anomalous low - temperature behavior of the bulk resistivity@xcite ( see also ) . several metallic contacts show anomalies with different signs depending on the way the contact is prepared . in the following we concentrate on a few typical and most interesting cases . the zero bias anomalies were widely studied in copper nanocontacts created with different methods . the first experiments by touching two cu electrodes@xcite reported a pronounced positive anomaly at @xmath723mev [ fig .  [ akimenko.fig]a ] . the zba has a highly asymmetric shape , which can be well fitted with the theoretical prediction for the elastic scattering on a single slow tls . the inelastic scattering could be excluded , as in this case the anomaly should show a basically symmetric peak with a background of opposite sign [ see fig .  [ secder_ballq ] ] . furthermore , as the zba has consequently shown a positive sign , a kondo - like explanation could also be excluded . in some occasions further low - energy peaks were also observed [ fig . [ akimenko.fig]b ] , which were argued to be due to the higher excitations of the tls .    , title="fig : " ] , title="fig : " ]    similar zbas were seen in the experiments of keijsers et al.@xcite on copper thin film mechanically controllable break junctions [ fig . [ posneganomaly]a ] . these samples were created by vapor deposition , and judging from the residual resistance ratio and the intensity of the phonon peaks in the pc spectrum they were of low crystalline quality implying that the contacts were in the diffusive regime . however , the authors claim that the zba is not due to the intrinsic disorder of the sample , but it correlates with the amount of bending and stretching of the film . in these thin film samples the zba has consequently shown a positive sign , and it could be well fitted with the elastic scattering on a slow tls , similarly to the above discussed measurements of akimentko et al . the same work of keijsers et al.@xcite has also presented measurements on single crystal cu break junctions [ fig .  [ posneganomaly]b ] . in this case the well resolved electron - phonon interaction spectrum implies a large elastic mean free path . in these measurements the zba has consequently shown a negative sign . here , the elastic scattering on a slow tls is still a reasonable explanation , however , the kondo - like behavior also has to be considered . ralph et al.@xcite have investigated nanofabricated cu contacts . in this case the size of the contact naturally could not be varied during the measurements , but as a great advantage , the contact was `` fixed '' , and thus a detailed study of the temperature and magnetic field dependence could be performed . in these experiments a well defined negative anomaly was observed if the sample was quenched [ fig . [ ralphzba.fig ] ] , and the anomaly disappeared after annealing . the detailed theoretical analysis of the curves ( including scaling arguments)@xcite has shown very good agreement with the theory of fast kondo - like tlss , even if their existence is still an open question.@xcite the relatively high kondo temperature ( @xmath724k ) can be obtained in the framework of recent theory@xcite only if the center is light ( e.g.  h or he ) as the energy scale is the kinetic energy.@xcite    mv is attributed to the scattering on fast two level systems , while the structures at higher voltages originate from the scattering on phonons . [ taken from with kind permission of the american institute of physics . ] ( b ) differential conductance of a cu nanoconstriction at @xmath725mk . the curve measured in @xmath726 t magnetic field is shifted down by @xmath727 for clarity , and only a slight dependence on the field is observed , which supports the nonmagnetic origin of the anomaly . for comparison the inset shows the conductance of cu constriction with @xmath728ppm mn at @xmath160mk , showing a pronounced zeeman splitting in magnetic field . [ taken from with kind permission of the american institute of physics . ] similar zero bias anomaly was also observed in ti nanojunctions.@xcite , title="fig : " ] mv is attributed to the scattering on fast two level systems , while the structures at higher voltages originate from the scattering on phonons . [ taken from with kind permission of the american institute of physics . ] ( b ) differential conductance of a cu nanoconstriction at @xmath725mk . the curve measured in @xmath726 t magnetic field is shifted down by @xmath727 for clarity , and only a slight dependence on the field is observed , which supports the nonmagnetic origin of the anomaly . for comparison the inset shows the conductance of cu constriction with @xmath728ppm mn at @xmath160mk , showing a pronounced zeeman splitting in magnetic field . [ taken from with kind permission of the american institute of physics . ] similar zero bias anomaly was also observed in ti nanojunctions.@xcite , title="fig : " ]    in all of the above experiments the zbas were attributed to a single or a few tlss . in the measurements , where consequently a positive zba was observed both the inelastic scattering on a slow tls and the kondo - like behavior can be ruled out due to shape and the sign of the anomaly , respectively . therefore , the elastic scattering on a slow tls remains as an explanation , which indeed gives good fit for the curves . however , it should be noted that for elastic scattering the sign of the anomaly is expected to be contact dependent . a consequent positive sign might be explained by the following argument : the atom in the higher level of the tls can be considered to be more far away from the ordered arrangement of the atoms , thus its scattering cross section is larger than that of the low energy state , which results in an enhanced backscattering as the bias is increased . it is interesting to note that the positive zbas in exhibit a pronounced contact size dependence : as the diameter of the junction is decreased the zba shifts to higher energy values [ fig . [ posneganomaly]a ] . the zbas showing kondo effect in magnetically doped contacts have exhibited a similar size effect.@xcite the latter phenomenon was explained by zarnd and udvardi as a result of density of states oscillations nearby the surface.@xcite similar size effect is expected for the kondo - like scattering on fast tlss , but we stress that this explanation does not apply to the zbas in fig .  [ posneganomaly]a due to the positive sign of the anomaly . the zero bias anomalies are also characteristic for metallic glass samples , which were always found to show negative anomaly@xcite [ for examples see fig . [ telegraphzba ] ] in accordance with the anomalous low - temperature increase of the bulk resistance.@xcite in metallic glasses tlss with a wide distribution of parameters can be present , thus the interpretation of the results is even more difficult . usually the zba is attributed to the scattering on tlss , however the magnetic origin also can not be excluded . the anomalies can not be suppressed by magnatic field in contrast to the magnetic kondo effect , however the anomalous part of the temperature dependence of the resistance shows a small but well defined magnetoresitance , which is not present at higher temperatures.@xcite    another interesting feature is the telegraph fluctuation of the zba demonstrated by the upper curves in fig . [ telegraphzba ] . this telegraph noise is only present in the narrow voltage region of the zba , that is the amplitude of the fluctuation vanishes as the voltage is increased ( fig . [ diffzba ] ) . this extraordinary behavior may be understood by assuming that a slow two level fluctuator is not directly changing the resistance ( like in fig . [ buhrmantls ] ) , but it is modulating the parameters of a fast tls in the neighborhood . therefore , the amplitude of the fluctuation only depends on the parameters of a _ single _ or a few tlss with the slowly moving defect being in one or the other metastable position , and the contributions of electron - phonon interaction , or other tlss being more far away from this particular defect are less significant . comparing the voltage dependence of the fluctuation with theory , it turns out that the model of fast kondo - like tlss is in good agreement with the measurement , while the theory of slow tlss gives a poor fit because of the long tail of the pc spectrum.@xcite      the common feature in the yet reported zbas is that the anomaly appears in the voltage region @xmath729mv . usually the anomaly shows a negative sign , however even in these cases the elastic contribution from slow tlss should be considered . a direct distinction between the fast and slow tlss could be made by studying the relaxation dynamics of the system . in terms of the two models the position of the singularity reflects the kondo temperature ( @xmath0 ) or the energy splitting ( @xmath40 ) , respectively . inserting the typical voltage range of the anomaly into the models , the standard korringa - like relaxation rate of slow tlss due to the creation of an electron and hole pair is @xmath730s ( see eq . [ relaxtime ] ) , while in the two channel kondo model the relaxation time is in the range of @xmath731s.@xcite therefore the boundary between the two cases can be investigated with microwave irradiation measurements in the ghz frequency range . experimentally the high frequency measurements were usually performed by irradiating the contact with a chopped microwave signal , and measuring the change of the dc current as the microwave is on and off . according to eq . [ iexpansion ] the voltage dependence of the dc shift signal is proportional to the second derivative of the @xmath711 curve if the frequency of the irradiation signal is slow compared to the relaxation of the tls , and the perturbation due to the alternating voltage is weak . previously , it was believed that above the characteristic frequency of the tls the dc shift signal should vanish , however the calculations were only performed for the second generation signal in the case of elastic scattering.@xcite measurements on various systems have shown that the peak in the dc shift signal can not be suppressed in the frequency range of the measurements ( @xmath732ghz).@xcite therefore , these studies concluded that the anomaly comes from fast tlss with a characteristic frequency @xmath733s@xmath734 . however , according to our results in sec . [ expansion.sec ] these measurements are not decisive as the dc shift signal _ should not vanish for slow tlss either _ , it should be only slightly reduced . ( see upper panels in fig . [ mwin.fig ] and fig . [ mwel.fig ] ) . the arbitrariness of the substraction of the background signal makes it even more difficult to draw a clear conclusion concerning the frequency dependence of the peak . more conclusive results can be obtained by measuring the second harmonic response to the high frequency irradiation , which is , however , experimentally more difficult . according to our calculations , in the second harmonic signal the peak due to the elastic scattering should disappear at high frequency , while the peak due to inelastic scattering only decreases by a factor of  2/3 ( bottom panels in fig . [ mwel.fig ] and [ mwin.fig ] , respectively ) . recent second harmonic generation measurements on ni@xmath735nb@xmath736 metallic glasses have shown,@xcite that the zba can not be suppressed with frequencies up to @xmath737ghz . this measurement excludes the elastic contribution of slow tlss due to the lack of frequency dependence , and it also excludes the inelastic contribution due to slow tlss as the sign of the peak is negative . however , the authors propose that the frequency dispersion of the smooth background signal might come from the scattering on slow tlss . an interesting result has been presented by balkashin et al . on annealed ni@xmath735nb@xmath736 metallic glass samples.@xcite the dc shift response signal was measured at various microwave irradiation frequencies both before and after annealing . the unannealed samples have shown insignificant frequency dispersion up the maximum frequency of the measurements ( @xmath737ghz ) , in agreement with the previous studies . in the annealed samples , however , the zba was completely suppressed at high microwave frequencies [ fig . [ annealing.fig ] ] . as we have discussed , the zba in the dc shift response signal should not dissappear for slow tlss either , thus the experimental observations on the annealed samples can not be interpreted in terms of the presented theory for slow tlss . this observation questions whether the zba originates from the scattering on tlss in the studied ni@xmath738nb@xmath739 metallic glasses .    up to now , the high frequency measurements were only performed on metallic glass samples . unfortunately in metallic glasses various interactions coexist , which makes the interpretation of the results difficult . especially serious problem is the normalization of the measured response signals . in pure samples the curves can be normalized to the phonon peaks , as the electron - phonon relaxation time is very small ( @xmath740s ) . in metallic glasses , however the phonon peaks can not be identified , and a smooth background is observed , which can contain various contributions . in disordered systems the disorder induced electron - electron interaction ( see e.g. ) and the scattering on nonequilibrium phonons@xcite are also important , furthermore the scattering on a variety of tlss with different excitation energies can also contribute to the background . additionally , the materials under study contain magnetic ions , thus the interactions of magnetic origin also can not be excluded . due to the complexity of the system the background signal can also have frequency dispersion in the frequency range of the measurements , which hampers the comparison of the response signals measured at different frequencies . a more clear picture could be obtained by investigating the high frequency behavior of the zbas observed in clean contacts of pure metals . in this case the zba is attributed to the scattering on a single or a few tlss , and due to the ballistic nature of the contact clear phonon peaks are observed at higher energies . in this case the microwave response signals are easily compared by normalizing the curves to the phonon peaks , and thus the measurement of the second harmonic response signal provides explicit information about the relaxation dynamics of the system . especially interesting would be the study of positive zbas in copper nanojunctions ( or in other metals like ag , au and pt ) , which are explained by the elastic scattering on slow tlss . nb@xmath736 with ag in measurements at different frequencies : a sonic frequency ( 1 ) and microwave frequencies @xmath741 0.48ghz ( 2 ) , 2.036ghz ( 3 ) , and 4.72ghz ( 4 ) . the anomaly is strongly suppressed with the microwave irradiation , whereas the zba in unannealed samples is unaffected by the same irradiation . [ taken from with kind permission of the american institute of physics . ] ]      a large variety of point contact spectroscopy measurements show zbas , which may be attributed to the scattering of electrons on tlss . certain measurements on cu nanocontacts show positive zbas , which can be well explained by the elastic scattering on slow tls . in this case , however , the reason for the consequently positive sign of the anomaly is not well understood . another group of measurements show negative zbas , and among these the experiments on cu nanocinstrictions@xcite are well interpreted in the framework of the two channel kondo model . in these studies the zbas are attributed to a few or even a single tls , which allows a direct comparison of the experimental observations with the theoretical predictions . zbas are frequently observed in metallic glass samples , however , in this case the complexity of the system makes the interpretation of the results more difficult . interesting attempts were made by investigating the frequency dispersion of the zbas in metallic glasses , however some of the results were misinterpreted due to the lack of the theoretical calculation for the frequency dependence of the dc shift signal . novel second harmonic generation measurements on ni@xmath738nb@xmath739 metallic glasses has proved that in this system the zba can not be explained by the scattering on slow tlss , however slow tlss can give contribution to the background . the study of the high frequency behavior of a few tlss near clean metallic contacts might provide a more clear information about the relaxation dynamics of tlss .    in the cases cited above the origins of the tlss are not well known . tlss can occur , however , in crystalline material with deviation from the stochiometric composition or disorder as well . the tls is formed in crystalline neighborhood where the size or the shape of the cavity in which the atom is sitting can be modified by the surroundings . in such cases clear indication of orbital kondo effect was found in pbte if some pb atoms are replaced by smaller ge atoms.@xcite however , this material is not so easy to produce in a controlled way . recently a new family of material was discovered , which are metallic and they are composed of u , th , as , se and s.@xcite the systematic study of these well - controlled materials may also provide the possibility to understand the tls - electron interaction in more detail . the point contact spectroscopy is a very powerful method to investigate the dynamics of slow tlss and magnetic kondo impurities and resonances of unknown , nonmagnetic origin , presumably due to some structured defects at zero bias . that spectroscopy is complementary of tunneling spectroscopy with oxide barrier@xcite or break junction in the tunneling limit,@xcite but there are similarities and differences as well . the local excitations in the barrier as molecular vibrations or magnetic spin flip processes open new tunneling channels resulting in extra current at voltages above the excitation energy . that gives just the opposite sign compared to the point contact case , as those anomalies are the negative ones . the magnetic impurities in the electrodes nearby the barrier can also depress the local tunneling density of states around the fermi level in the energy range of the kondo temperature.@xcite that gives a negative zero bias anomaly just like in the point contact case .    in the future it would be very important to make systematic studies concerning the tls , starting with the slow tls with contributions of arbitrary signs due to dominating elastic processes . the detailed studies of the frequency dependence based on the theoretical part of the present review would be especially useful and that would clarify some of the contradictions of the present time . it is very important to study further the resonances of unknown origin which have strong resemblance to nonmagnetic kondo behavior@xcite and to make attempts to determine their origin in at least a few cases . the origin of those resonances has been the subject of theoretical debates in the recent years . originally it was proposed that kondo effect could arise from tunneling tlss with non - commutative couplings e.g. due to the screening by electrons and electron assisted tunneling ( see e.g.  ) . it was shown , however , that the adiabatic screening of the moving atom by electrons depresses the tunneling rate in such a strong manner,@xcite that the actual kondo temperature is negligible . recently , it has been suggested that the situation can be changed if the atom in the double potential well has the ground - state below the potential barrier , but the first excited state is above it.@xcite then this objection could be avoided , but the value of the splitting raises many questions . that theory for an intermediate heavy atom gives too small @xmath0 , not more than @xmath742k and that can not explain the resonances of typical width of @xmath743k . there is , however , the possibility of light atoms like hydrogen,@xcite as the dominating energy scale is the kinetic energy of the atom in the potential well . therefore , it would be very important in the future to study hydrogen defect either in atomic states or water contamination , where the hydrogen could still change its position . summarizing , the origin of slow tls and the dubious resonances deserve extended further studies . we acknowledge to all the authors laying down the principles of slow tlss in point contacts . we are especially thankful to i.k . yanson and sz . csonka for the careful reading of the manuscript , and also to k.e . nagaev , j. von delft , j.m . van ruitenbeek , o.p . balkashin , a.n . omelyanchouk , i.o . kulik , g. zarnd , d. c. ralph , g. mihly , r. a. buhrman and h. van kempen for their useful remarks . part of this work was supported by the hungarian research funds otka t026327 , ts040878 , t037451 , t038162 , t034243 .
Abstract:

Answer: a great variety of experiments , like heat release measurements , acoustic measurements , and transport measurements on mesoscopic samples have proved that two level systems ( tlss ) have a crucial role in the low temperature thermal and electric properties of disordered systems . this paper is aimed at reviewing the role of slow tlss in point contacts . first the theory of point contacts is summarized , concentrating on the discussion of different point contact models , and on the different regimes of electron flow in the contact , mainly focusing on the ballistic and diffusive limit . the boltzmann equation is solved in both regimes , and the position dependence of the electrical potential is determined . then the scattering processes in point contacts are investigated , particularly concentrating on the scattering on slow tlss . if the the electron assisted transitions between the two states are negligible the electron - two level system interaction can be treated with a simplified hamiltonian . the scattering on such slow tlss causes nonlinearity in the current - voltage characteristics of the point contact , which can be determined using fermi s golden role . these calculations are presented showing both the contribution of elastic and inelastic scattering , and including the dependence on the position of the tls , and on the effect of high frequency irradiation . these results are used to discuss the differences between these slow tlss and the fast centers which may be described by the two channel kondo model . the available experimental results are analyzed , distinguishing between the effects due to the different types of tlss .