the 2.16-m reflector is an english equatorial mount telescope at xinglong observatory and it obtained the first light in 1989 , which was built / developed by chinese ourselves independently .
the effective aperture of the telescope is 2.16 meter , and the focal ratio of the primary mirror is f/3 @xcite .
the effective aperture of the secondary mirror is 0.717 m. currently , there are two available foci for mounting astronomical observing instruments , the cassegrain focus and the coud focus . for the cassegrain focus , which is a r - c system , the focal ratio is f/9 , and the scale on the focal plane is @xmath3/mm .
while for the coud system , it is f/45 , and the scale on the focal plane is @xmath4/mm .
the most special character of the optical system is that the cassegrain system and the coud system share the same secondary mirror and there is a relay mirror in the coud system , and both systems can sufficiently eliminate the spherical aberration and coma aberration @xcite .    since its first light in 1989 ,
the 2.16-m telescope has been used in various scientific research fields in galactic and extragalactic astrophysics , including the determination of stellar parameters ( e.g. , abundances , surface gravity , temperatures ) of a large sample of stars , the discoveries of substellar and planetary companions of stars , studies of active galaxy nuclei ( agn ) , including the identifications of high - redshift quasars , discoveries and studies of supernovae , e.g. , 1993j @xcite , as well as time - domain science ( e.g. , supernovae , gamma - ray bursts , stellar tidal disruption events , and variable stars ) .
although the 2.16-m telescope is the third largest optical telescope in china right now ( smaller than lamost and the lijiang 2.4-m telescopes , which were installed in 2008 and 2007 , respectively ) , it plays an important role in chinese astronomical observations .
every year there are hundreds of chinese astronomers applying for observing time at the 2.16-m telescope , including significant numbers of new telescope users and graduate students .
therefore , it is of great importance and useful to describe the specific parameters and the observing ability of the telescope and its instruments to the telescope users , which could be very helpful for planning and carrying out observations and for the data reduction .
this paper is organized as follows : the xinglong observatory is introduced in sect.[obs.sec ] ; the telescope and its instruments are described in sect.[tel.sec ] and the efficiencies of the telescope and its instruments also have been estimated in this section .
the relevant science and projects based on observational data obtained by the 2.16-m telescope are described in sect.[sci.sec ] ; finally a brief summary is given and some future facilities are discussed in sect.[sum.sec ] .
the xinglong observatory of national astronomical observatories , cas ( naoc ) ( iau code 327 , coordinates : @xmath5 n , @xmath6 e ) was founded in 1968 . at present , it is one of most primary observing stations of naoc . as the largest optical astronomical observatory site in the continent of asia
, it harbors 9 telescopes with effective aperture greater than 50 cm .
these are lamost , the 2.16-m reflector , a 1.26-m optical & near - infrared reflector , a 1 m alt - az reflector , a 85 cm reflector ( naoc - beijing normal university telescope , nbt ) , a 80 cm reflector ( tsinghua university - naoc , tnt ) , a 60/90 cm schmidt telescope , a 60 cm reflector and a 50 cm reflector .
the average altitude of the xinglong observatory is @xmath7 m and it is located at the south of the main peak of the yanshan mountains , in the xinglong county , hebei province , which is @xmath8 km northeast of beijing . the mean and median seeing values of the xinglong observatory are @xmath9 and @xmath10 , respectively , and on average , there are 117 photometric nights and 230 useful nights per year based on the data of 2007 - 2014 @xcite . for most of the time ,
the wind speed is less than 4 m s@xmath11 ( the mean value is 2 m s@xmath11 ) , and the sky brightness is @xmath1221.1 mag arcsec@xmath13 in @xmath14 band at the zenith @xcite .
each year , more than a hundred astronomers use the telescopes of xinglong observatory to perform the observations for the studies on galactic sciences ( stellar parameters , extinction measurements , galactic structures , exoplanets , etc . ) and extragalactic sciences ( including nearby galaxies , agns , high - redshift quasars ) , as well as time - domain astronomy ( supernovae , gamma - ray bursts , stellar tidal disruption events , and different types of variable stars ) . in recent year , besides the basic daily maintenance of the telescopes , new techniques and methods have been explored by the engineers and technicians of xinglong observatory to improve the efficiency of observations
. meanwhile , the xinglong observatory is also a national popular - science and education base of china for training students from graduate schools , colleges , high schools and other education institutions throughout china , and it has hosted a number of international workshops and summer schools .
as shown in figures  [ fig1 ] and [ fig2 ] , the telescope is a r - c system on an english equatorial mount .
the effective aperture of the primary mirror is 2.16 meter .
currently , there are three primary instruments available for the cassegrain focus : bfosc , omr and the fiber - fed hrs .
previously , a high resolution spectrograph was mounted on the coud focus before the fiber - fed hrs was mounted in 2010 . since then , the coud focus has only been used for a few special experiments , such as adaptive optics tests . as shown in figure  [ fig1 ]
, the diameter of the dome is 22 meter , which is relative large for a 2m - class telescope , due to its design .
the rotation part of the telescope weights 91 tons and the pointing accuracy which is modified with pointing model is @xmath15 . while the tracking precision of the telescope is @xmath16 within 10 minutes , and @xmath17 within 1 hour .
the position precision of the telescope on the sky under guiding is @xmath18 @xcite .
the bfosc is one of the primary instruments of the telescope , for which the design and processing of mechanism and electronics control , the design of grims and assembling and debugging of the whole system were done by university of copenhagen and eso provided consult and optical design of the focal reducer .
it is available for the f/9 cassegrain focus .
the scale of the focal plane is @xmath3/mm . as discussed in sect .
[ intro.sec ] , it can be used for both imaging mode and spectroscopy mode , which is switchable .
figure  [ fig3 ] shows the optical layout of the bfosc instrument , including the aperture wheel , filter wheel , grism / echelle wheel , calibration device , guiding device , collimator , shutter , camera and the ccd detector .
there are eight positions on the aperture wheel : three of them for direct imaging , coronagraph mask and focal adjusting plate , respectively and the other five for long / short slit plates . for the long slits , the lengths are all @xmath19 and there are nine options for the slit widths ( @xmath20 , @xmath21 , @xmath22 , @xmath23 , @xmath10 , @xmath24 , @xmath25 , @xmath26 and @xmath27 ) .
while for the short slits , the slit lengths are various for different slit width : @xmath28 for slit width of @xmath20 , @xmath29 for slit width of @xmath30 , @xmath25 for slit width of @xmath31 and @xmath32 for both slit widths of @xmath24 and @xmath33 .
the coronagraph mask is used to block out the intense light from bright sources near the observing targets , with circular spots of which the diameters are @xmath9 , @xmath34 , @xmath29 , @xmath35 , @xmath36 and @xmath37 .    for the filter wheel ,
there are several sets of filters available .
the johnson - cousins @xmath38 filters are for the broadband photometry .
in addition , there are special filters used for the spectroscopic observing : the @xmath39 band for the transmitting spectral region @xmath40 nm , and the @xmath41 band for removing the 2nd - order spectrum in the wavelength region @xmath42 nm .
the transmission curves are shown in figure  [ fig4 ] .    besides , a series of interference narrow band filters covering [ oiii ] ( rest - frame wavelength of @xmath435007a ) , @xmath44 and heii ( rest - frame wavelength @xmath434686a ) are also available for the bfsoc imaging mode . for series of [ oiii ]
band observations , there are 8 filters of which the central wavelengths are between 500.9 - 536.0 nm with fwhm of 6 nm , corresponding to redshifts of 0 - 20000 km s@xmath11 . while for the series of @xmath45 bands , there are 11 filters of which the central wavelengths are between 656.2 - 706.0 nm with fwhm of 7 nm , corresponding to redshifts of 0 - 22000 km s@xmath11 . for the he bands , there are three filters of which the central wavelengths are between 447.1 - 468.6 nm with fwhm of 6 nm . in addition ,
an [ oiii ] filter with fwhm of 13 nm and an @xmath45 filter with fwhm of 14 nm are also provided for a redshift of z=0 .
these filters can be used for the study of star - forming regions of nearby galaxies at different redshifts .
table  [ t1.tab ] presents the parameters of the grisms / prisms / echelle for the bfosc instrument . from left to right in columns are names , working spectral orders , reciprocal linear dispersions , dispersions and wavelength coverages of the various grisms / prisms / echelle , respectively .
the configurations can be chosen by the users depending on the different requirements of the projects .
for grisms g3/g6/g10 , the blue ends of wavelength coverages are limited by the cutoff of the atmospheric window , while for grisms g5/g8 , the red ends are limited by the size of the ccd .
the low dispersion grisms g10/g11/g12 are also used for the cross - disperser when mounted on the filter wheel .
the echelle e13 is only used for measuring the velocity field of extended sources with the 3rd order spectrum , and the @xmath14 band filter is recommended to be applied together with to remove the other order of the spectrum .    ccccc name & spec.ord . & rec .
disp . & disp . &
range + & & (  / mm ) & (  / pix ) & (  ) + p1 & & @xmath46 & @xmath47 & @xmath48 + g3 & 1 & 139 & 3.12 & @xmath49 + g4 & 1 & 198 & 4.45 & @xmath50 + g5 & 1 & 199 & 4.47 & @xmath51 + g6 & 1 & 88 & 1.98 & @xmath52 + g7 & 1 & 95 & 2.13 & @xmath53 + g8 & 1 & 80 & 1.79 & @xmath54 + g10 & 1 & 392 & 8.80 & @xmath49 + g11 & 1 & 295 & 6.63 & @xmath55 + g12 & 1 & 837 & 18.8 & @xmath56 + e9+g10 & @xmath57 & @xmath58 & @xmath59 & @xmath49 + e9+g11 & @xmath60 & @xmath61 & @xmath62 & @xmath55 + e9+g12 & @xmath63 & @xmath64 & @xmath65 & @xmath56 + e13+v & 3 & 33.1 & 0.666 & @xmath66 +    table  [ t2.tab ] presents the spectral resolutions of some frequently used bfsoc grisms at a minimum slit width @xmath20 and a slit width of @xmath67 used in typical seeing condition , which are estimated with the emission lines of the planetary nebula ( pn ) ic4997 .
the observations were taken on november 24 , 2014 and february 26 , 2016 , with gratings of g4 , g6 , g7 and g8 on bfosc . throughout the nights of the observations ,
the weather were clear and the seeing were between @xmath68 .
the exposure time was 0.2 - 10 seconds depending on the dispersions of the different grisms .
the shortest exposure of 0.2 second is considerably longer than the shutter speed of 1.5 milliseconds .
c|cccc|cccc    wavelength & g4 & g6 & g7 & g8 & g4 & g6 & g7 & g8 + (  ) & & + 4341 & 531 & 1423 & 1204 & & 246 & 463 & 478 & + 4363 & 557 & 1503 & 1252 & & & & & + 4861 & & & & & 262 & 538 & 521 & + 4959 & & & & & 265 & 540 & 524 & + 5007 & 620 & 1413 & 1328 & & 265 & 555 & 532 & + 6563 & 824 & & & 2245 & 321 & & & 820 +    for the calibration device , four lamps can be mounted at most .
one is for the flat - fielding correction and the other lamps can be used for wavelength calibration
. the lights of the lamps illuminate the integrating sphere at first , then modified to be f/9 light beam and finally reflected to the focal plane of the bfosc instrument . the fe lamp and ne lamp
are frequently used for the wavelength calibration .    on the guiding device ,
a mirror reflects the view outside of the observing field of view ( fov ) to an iccd , which is movable in 3d directions to find a proper guide star and also to adjust the focus .
the fov of the guiding device is @xmath69 , which is large enough to find a suitable guide star .    in 2010
an e2v @xmath70 back illuminated ccd , aimo was installed on the spectrograph and the ccd controller was made by the lick observatory .
the size of the ccd is @xmath71 pixels with pixel size of 22.5 @xmath72 m .
the pixel scale is @xmath73 and the fov is @xmath74 according to the size of ccd .
figure  [ fig5 ] shows the quantum efficiency ( q.e . ) of the ccd .
it can be seen that the maximum quantum efficiency is higher than 90% around 5700  of the wavelength .
table  [ t3.tab ] lists the gains , readout noises and readout speed of the bfosc ccd at various readout times .
observers can choose different options for the specific observations and the slow readout speed is applied in most of the time .
ccccc readout speed & mode & gain & readout noise & readout time + & & ( @xmath75/adu ) & ( @xmath75/pix . ) & ( second ) + fast & 0 & 99.58 & 132.5 & 5 + fast & 1 & 49.44 & 103.7 & 5 + fast & 2 & 24.48 & 87.35 & 5 + median & 0 & 2.33 & 6.80 & 8 + median & 1 & 1.22 & 5.93 & 8 + median & 2 & 0.50 & 6.06 & 8 + slow & 0 & 2.43 & 3.25 & 28 + slow & 1 & 1.13 & 2.58 & 28 + slow & 2 & 0.50 & 2.46 & 28 +    there are six observing modes for the bfosc instrument : ( 1 ) direct imaging ; ( 2 ) long slit spectroscopy ; ( 3 ) slitless spectroscopy ; ( 4 ) echelle grism spectroscopy ; ( 5 ) coronograph mask ; and ( 6 ) multiple object spectroscopy .
figure  [ fig6 ] shows the optical layout of the lens and light path for the direct imaging mode and the spectroscopy observing mode , which could be switched in a few minutes .        since 2012
, a multi - object spectroscopy ( mos ) observing mode of bfsoc is available for observers by placing a multiple - aperture mask on the aperture wheel @xcite . in this mode ,
10 - 20 objects can be observed simultaneously , depending on the spatial distribution density of the targets .
the mos improves the observing efficiency of multiple - object observations , such as star - forming activities of hii regions of nearby galaxies , star clusters , or groups and clusters of galaxies .    the total efficiency of the observing system , including the atmosphere , telescope and its instruments , is what the observers are mostly concerned about . in fact
, there are a number of factors to be considered for calculating the total efficiency of the system , such as atmospheric extinction , reflectivity of primary and secondary mirror , transmissions of filters and the quantum efficiency of the ccd .
it can be estimated through observing standard stars in the equation  [ eq1 ] ,    @xmath76    in this formula , @xmath77 is the observed number counts of a standard star per second ( adu s@xmath11 ) ; @xmath78 is the gain of the ccd ( @xmath79 adu@xmath11 ) ; @xmath80 is theoretical photon flux of a standard star derived from its ab mag ( @xmath81 @xmath11 ) ; @xmath82 is the effective band width of the filter in imaging observations or the dispersion of the grating for spectroscopic observations (  ) ; @xmath83 is the effective area of the primary mirror of the telescope ( @xmath84 ) and @xmath43 is the effective wavelength of the filter or the wavelength at which the efficiency is to be computed for the spectroscopy (  ) .    in order to estimate the total efficiency of the 2.16-m telescope , a landolt standard
star pg2336 + 004b was observed on november 21 , 2014 in the broadband @xmath38 bands of the bfosc photometric system by using equation  [ eq1 ] .
the seeing was @xmath85 throughout the whole night of observations and the airmass is @xmath86 .
figure  [ fig7 ] is a plot of the total efficiency in different bands as a function of central wavelength .
the atmospheric extinction , reflectivity of primary and secondary mirrors , transmissions of filters , q.e . of ccd as well as other factors were not corrected in the calculation .
it can be seen that the total efficiency is relative low in the @xmath87 band ( @xmath88% ) and @xmath89 band ( @xmath90% ) , but it is relatively high in @xmath91 bands ( @xmath92% ) .
bands at an airmass of @xmath86 . ]
the limiting magnitude is also an important quantity to evaluate the observing ability of a telescope .
actually at the same night on november 21 , 2014 , a landolt standard star pg2336 + 004 was observed in the @xmath38 bands .
the airmass was @xmath86 .
the limiting magnitude ( signal - to - noise ratio , @xmath93 ) as a function of exposure time in the @xmath38 bands are shown in figure  [ fig8 ] and the same relations but for the limiting magnitudes of @xmath94 are shown in figure  [ fig9 ] .
it is noted that in the 640s exposure the limiting magnitudes of the @xmath95 bands could reach @xmath96 mag in @xmath93 and @xmath97 mag in @xmath94 .
the limiting magnitude in the @xmath87 band is the lowest because the sensitivity of the observing system in the blue band is relatively low , further affected by the effect of significant atmospheric extinction in blue band .
however , although the efficiency in the @xmath98 band is similar to those of the @xmath99 bands , the night - sky background is much brighter ( @xmath88 times ) than in the @xmath99 bands .
in addition , the fringe of the ccd in the @xmath98 band contributes significant noise .
although we have tried the de - fringing method introduced by @xcite , it has not been improved significantly .
therefore , these factors together make the limiting magnitudes @xmath100 mag shallower than those in the @xmath101 bands
.     bands for the signal - to - noise ratio of @xmath93 at an airmass of @xmath86 . ]    , but for a signal - to - noise ratio of @xmath94 . ]    in order to estimate the efficiency of the spectroscopic system of bfosc , the two eso standard stars hr153 and hr9087 were observed on november 24 , 2014 , with several frequently used grisms ( g4 , g5 , g6 , g7 and g8 ) of the bfosc instrument .
the weather was clear and the seeing was @xmath102 .
the exposure time was 0.2 - 10 seconds according to the dispersions of the grisms . to make sure that most of the flux can be obtained and measured for estimating the efficiency , the slit was configured as @xmath26 .
the total efficiency , which is defined above , including atmospheric extinction and instruments , was calculated with equation  [ eq1 ] and shown in figure  [ fig10 ] . for the grisms g4 and g5 ,
the peak efficiencies are @xmath1215% and @xmath1210% respectively , while for the other three grisms , the total efficiencies of the peak are @xmath125% .
the filter @xmath41 was used for removing the 2nd - order spectrum of wavelength @xmath103385 nm .        for the limiting magnitude in the spectroscopic observations ,
the previous observations show that typically it could reach @xmath104 mag of @xmath93 in the grism g6 , with exposure of 1 hour , given by @xcite .
the omr is another low - intermediate resolution spectrograph available on the cassegrain focus of the 2.16-m telescope .
it was made by optomechanics research inc .
( arizona , usa ) at the end of 1994 and was tested on the 2.1 m telescope at kitt peak .
then , it was installed on the 2.16-m telescope of xinglong observatory in 1995 and was available to astronomers in 1996 .    the optical layout of the omr is shown in figure  [ fig11 ] .
the system is composed of slit & decker assembly , filters , spectrograph , calibration system , collimator , gratings , ccd camera , guiding ccd camera , console and data collecting system .
the slit & decker assembly is composed of decker , the slit width adjuster , calibration reflector , and driving device .
the top surface of the slit is aluminized stainless steel with high reflecting power while the slit edge was polished to be very sensitive and should not be touched .
the slit is tilted by @xmath105 with respect to telescope optical axis and the projected slit width can be adjusted remotely through the main console in the range of 0.05 - 1.0 mm ( corresponding to @xmath106 ) , which is displayed on the main console via the voltage value .
the size of the slit jaws is 32.8mm@xmath10738 mm and the effective length of the slit is 28.8 mm ( corresponding to @xmath108 ) .
the slit jaw reflects the light to the guiding system .
there are six positions on the filter wheel for the filters of clear , corning @xmath109 , schott @xmath110 , @xmath111 , @xmath112 and @xmath113 .
all the filters have the same size of 25mm@xmath10725 mm .
the central wavelengths of filters could be matched with the 1st or 2nd order of the spectrum by the softwares automatically .
there are three wavelength calibration lamps ( fe - ne , fe - ar and he - ar ) and a flat - fielding lamp for the calibration device .
all of lamps can be controlled remotely .
for the fe - ne lamp , a standard 1.5-inch istc model # wl-22810 is supplied , and actually other types of lamps can be used as well , such as hamamatsu and starna .
the maximum working current is 20 ma and the normal working current is 10 - 12 ma . the fe - ar lamp istc model is # wl-22611 . for the safety of
he - ar lamp , the dc power should be used for the supply . as for the flat - fielding lamp ,
the power is the standard 12v dc , 1.5a halogen tungsten lamp ( sylvania#808 - 301550 ) and it usually can be used for 1000 - 2000 hours .
meanwhile , there are four condensers for the calibration lamps , which are made of quartz or pyrex glass .
the diameters are 25 mm and the focal lengths are @xmath114 mm .
the field lens is made of quartz , with a diameter of 38 mm and a focal length of 64 mm .
the mirror of the collimator is an off - axis parabolic aluminized reflector , made of pyrex glass .
the aperture of the mirror is @xmath115 mm and the focal length is @xmath116 mm .
the off - axis angle is 8.1@xmath117 and the focus can be adjusted remotely .    at present
there are six reflecting blazed gratings mounted on the omr spectrograph .
the gratings are made of aluminized pyrex glass and can be switched manually .
the parameters of the blazed gratings are given in table  [ t4.tab ] , including the number of the gratings , the grooves , groove areas , reciprocal linear dispersions , dispersions , blaze wavelengths of the 1st order spectra and the blaze angles .
the overall working wavelength coverage of the gratings is @xmath123700 - 10000  . for the specific grating
, the wavelength coverage is adjustable and it can be estimated through parameters of table  [ t4.tab ] and the size of the ccd .
for instance , the wavelength coverage of gratings of 1200 lines / mm , is @xmath118  while for the grating of 300 line / mm , it is @xmath119  .
the spec software , which is the camera controlling software designed and installed with the pi spec-10 ccd camera , can recognize the gratings and the wavelength coverage can be adjusted via spec according to the requirements of the observers .
the camera is a schmidt - cassegrain system .
the useful aperture is @xmath120 mm and the focal length is @xmath121 mm .
the ccd mounted is a pi spec-10 pixis @xmath122 scientific ccd detector , which delivers the highest sensitivity possible and @xmath123bit dynamic range for spectroscopy applications .
the pixel size is 20@xmath10720 @xmath72 m pixels and the ccd size is @xmath124 mm@xmath125 . the deepest cooling temperature is @xmath126 c. the quantum efficiency of the ccd is shown in figure  [ fig12 ] , which is midband .
the guiding system of the omr is composed of a reflecting system and a ccd camera . in the reflecting system ,
the slit reflecting mirror is aluminized and the r.m.s . of the surface flatness
is 1/4 @xmath43 . for the transfer lens
, the diameter is 50 mm , with focal length of @xmath127 mm . while for focal lens
, the diameter is 36 mm , with focal length of @xmath128 mm .
both of the lenses have been coated .
the ccd camera is an alta u42 ccd , made by apogee imaging systems , inc .
the sensor is e2v ccd@xmath129 and the size is 2048@xmath1072048 pixels .
the gain is 1.3e@xmath130/count and the maximum digitized well capacity is 82k e@xmath130 .
the dark current is 0.39e@xmath130/pixel / s .
ccccccc no . & grooves & gro .
area & rec .
disp . & disp . &
+ & ( l / mm)&(mm ) & (  / mm ) & (  / pix ) & ( 1st order ) (  ) & ( @xmath117 ) + 1 & 150 & @xmath131 & 400 & 8.0 & 5000 & 2.2 + 2 & 300 & @xmath131 & 200 & 4.0 & 8000 & 6.5 + 3 & 600 & @xmath132 & 100 & 2.0 & 6500 & 11.3 + 4 & 1200 & @xmath132 & 50 & 1.0 & 8000 & 28.7 + 5 & 300 & @xmath132 & 200 & 4.0 & 4224 & 3.6 + 6 & 1200 & @xmath132 & 50 & 1.0 & 4000 & 13.9 +    similarly , the landolt standard star feige 66 was observed with two blazing gratings 1200b and 1200r of the omr instrument and the total efficiencies ( including the atmosphere extinction , the reflecting rate of mirrors , q.e .
of ccd , etc . )
have been computed through equation  [ eq1 ] . as described above
, the central wavelength and coverage is adjustable , figure  [ fig13 ] shows the total efficiencies of two gratings 1200b and 1200r , in two different central wavelength and wavelength coverages .
it can be seen that the efficiency ranges from @xmath133 to @xmath134 , which is relatively lower than that of the bfosc gratings .
typically , for a star of @xmath135 mag in an exposure time of 1800s , the signal - to - noise ratio is @xmath136 at wavelength of @xmath43=5500  ; for a star of sdss @xmath137 mag in an exposure time of 3600s , the signal - to - noise ratio is @xmath93 on average for the whole band .
previously an echelle spectrograph was mounted on the coud focus of the telescope for high resolution spectroscopic observations .
the special design of the optical system is to allow the cassegrain system and the coud system to share the same secondary mirror and there is a relay mirror in the coud system . when switching between the two systems , the coud system can eliminate spherical aberration and coma aberration sufficiently just by slightly moving the secondary mirror @xcite .
the resolution power was between r=16,000 to 170,000 ( 79 gr / mm ) in the blue beam from 330 nm to 580 nm and r=13,000 to 170,000 ( 31.6 gr / mm ) in the red beam from 520 nm to 1100 nm . for an exposure time of 1 hour
, the limiting magnitude could reach @xmath138 in the red band and @xmath139 mag in the blue band with a signal - to - noise ratio @xmath140 @xcite .
since 2010 , a new fiber - fed high resolution spectrograph ( hrs ) has been developed by the nanjing institute of astronomical optics & technology ( niaot ) to satisfy the scientific requirement of , for instance , exoplanets surveys , the study of stellar abundances , and stellar magnetic activities . the optical layout and light path
are shown in figure [ fig14 ] .
the hrs system is available for the cassegrain focus and composed of four parts : 1 ) the cassegrain connector of the telescope ; 2 ) the fiber connector and the micro - optical module for focal - ratio changing ; 3 ) the main body of the spectrograph ; 4 ) the data collecting part of the ccd camera .
the fiber - fed hrs observing mode could be switched to other instruments ( such as bfosc or omr ) conveniently through the cassegrain focal interface of the telescope within only a few minutes .
the calibration system of hrs , the @xmath141 cell and its heating system , as well as the tip / tilt system are mounted on the interface of the spectrograph .
the whole main body is sealed in a protective box and the temperature variation is less than @xmath142c for a week . at present , a @xmath143 aperture fiber is configured and the tip / tilt system is working to improve the efficiency of the system when the seeing is ideal .
the fiber of @xmath31 is at the commissioning stage and it could be available very soon .
an environmental controlling system can accurately keep the stability of both temperature and humidity for the hrs system .
furthermore , a sub - controlling system of pressure will be installed at this system in the future .
the basic parameters of hrs are listed in table  [ t5.tab ] .
the working wavelength coverage is 360 - 1000 nm and the instrumental efficiency of the spectrograph is @xmath144 for the peak at a wavelength coverage of 640 - 790 nm ( @xmath145 band ) , and @xmath146 for @xmath147  ( the whole working band ) , based on the tests of november 23 , 2010 .
the spectral resolution is 32,000 - 106,000 for the spectrograph and it is @xmath148 at a fixed slit width of 0.19 mm ( corresponding to 18 ) based on the test of april 12 , 2011 .
the stability of the instrument for the velocity measurement is @xmath149 m / s and the temperature is quite stable even for two weeks . the ccd camera is a back - illuminated first - order red - sensitive e2v ccd 203 - 82 .
the size is @xmath150 pixels with a pixel size of 12 @xmath72 m .
the typical quantum efficiency of the ccd is @xmath151 under the temperature of @xmath152c in the wavelength coverage of @xmath153 nm , which is shown as the solid line in figure  [ fig15 ] .
the liquid nitrogen ( ln ) holding time of the system is @xmath154 hours , and the cooling temperature is @xmath155c . for @xmath156 e@xmath130/adu , the readout noise ( ron )
is 2.84 e@xmath130 at the readout speed of 50k , 4.29 e@xmath130 at readout speed of 100k and 7.88 e@xmath130 at readout speed of 200k .
the distance between the nearby orders of the spectrum is @xmath157 pixels and the two fibers ( @xmath143 aperture one and @xmath31 aperture one ) can work at the same time .
the aperture of the fiber is 100 @xmath72 m ( @xmath143 ) and the fiber length is @xmath158 m. the fov for the guiding camera is @xmath159 and in the guiding plate the aperture of pupil is @xmath29 in front of the fiber .
the temperature variation of the system is @xmath160c for a whole night and @xmath161c for a week .
the tip / tilt system is working normally on the cassegrain focus .
the system is currently in its commissioning phase , and its radial velocity precision dose not yet reach its design goal of a few @xmath12cm / s .
it is expected that this goal will be achieved once commissioning is complete .
the guiding ccd camera is a gc1380 , and the software @xmath162 is used for ccd controlling , exposure time configuration and data collection .
the telescope controlling software can monitor the guiding uncertainty of the guiding system .
an astro - frequency comb calibration system , which is developed by peking university ( pku ) , also was installed and is being tested on the 2.16-m telescope .
the full spectral wavelength range of the astro - frequency comb is 160 nm and the central wavelength is @xmath163 nm . the observed spacing of the comb teeth is 29.01 ghz . once working normally , it can greatly improve the measuring accuracy of the stellar radial velocity .
the system has significant advantages over the @xmath141 cell @xcite .
it consists of a series of discrete , equally spaced spectral lines with equal intensity and it is repeatable in long time scale . right now
, the system has been set up and the first spectrum simultaneously from the comb and the flat - fielding lamp has been obtained during the engineering run .
meanwhile , a 25ghz astrocomb optical frequency calibration system manufactured by the menlosystems company , also has been installed on the hrs instrument system recently .
the spectral coverage is @xmath12450 - 720 nm , with the central wavelength of @xmath164 nm . the flatness of the spectrum is r.m.s . @xmath165
db and the observed spacing of the comb teeth is 25 ghz .
the luminous power of the system for the full spectral wavelength range is @xmath166w .
the two new calibration systems are at the commissioning stage .
the two new systems are supposed to be applied in the end of this year .
c , which is shown by the solid line ]    ccc tech .
index & parameters & notes + wave .
3650 - 10,000  & thar lamp + aper . of fiber & @xmath143 / @xmath31 & + spec .
res . & 32,000 - 106,000 & tested on apr . 12 , 2011 + stab .
of inst . &
@xmath167 & + effi . & @xmath168 at peak value & + & @xmath169 % for @xmath147  & tested on nov . 23 , 2010 + & e2v 4k@xmath1074k 12@xmath72 m & + & scientific chip , class 1 & + & on - chip binning & + ccd camera & 50 - 200 khz & 2@xmath1071 binning + & ln holding time : 20 hours & mode suggested + & ( w / ln auto - filling ) & + & ln cooling & + fov of guid . cam . &
@xmath170 & + temp .
/ day & @xmath171 & + temp .
/ week & @xmath172 & +    in order to estimate the total efficiency ( including the atmospheric extinctions and reflectivity of the primary / secondary mirrors , among others ) , a number of standard stars were observed with hrs on june 30 , 2015 . by using equation  [ eq1 ] , the final result shows that it is @xmath1732% at 5500 .
in addition , the limiting magnitudes are estimated .
figure  [ fig16 ] is a plot of the signal - to - noise ratio ( snr ) at @xmath43=5500  as a function of exposure time for the stars from @xmath174 mag to @xmath175 mag for the hrs system , based on the observing data . for a @xmath175 mag star , the typical signal - to - noise ratio @xmath176 for an exposure time of 3000 - 3600 seconds .
a few observing tests give the snr in the typical observing conditions of 2012 , for a star of @xmath177 mag with an exposure time of 3600s , the signal - to - noise ratio is @xmath178 at wavelength of @xmath179  , and for a star of @xmath180 mag in 2400s exposure , the signal - to - noise ratio is @xmath181 at wavelength of @xmath43=6000  ( provided by xiaoling yang and yuqin chen ) .
= 5500  as a function of exposure time for stars from @xmath174 mag to @xmath175 mag by using the fiber - fed hrs spectrograph . ]
since 1989 , when the 2.16-m telescope saw its first light , a great many scientific research projects have been carried out in various fields , including : the study of nearby galaxies ( star formation rates , gas and dust content ) , agns and their supermassive black holes ( smbhs ) , quasars , stellar parameter determinations , exoplanets , supernovae ( sne ) , gamma - ray bursts ( grbs ) , and follow - ups of tidal disruption events .
when lamost was built , the 2.16-m telescope has also been used for the lamost follow - up and stellar library observations .
some of the scientific highlights obtained with the 2.16-m telescope are as follows . in 1993
, @xcite observed the spectra of supernova 1993j with the 2.16-m telescope and found the blue - shifted oxygen lines only four months after the optical discovery of the supernova , which is different from the case of 1987a .
@xmath182 g - type giant stars of @xmath183 mag have been monitored with the coud hrs system @xcite . in a joint planet - search program between china and japan by using the xinglong 2.16-m telescope and okayama astrophysical observatory ( oao ) , a number of substellar companions of intermediate - mass giant stars has been discovered by @xcite .
later , a brown dwarf companion candidate was discovered by @xcite and a long - period giant planet has been discovered by @xcite , which is actually the first exoplanet discovered jointly with the subaru telescope , the oao 1.88 m telescope and the xinglong 2.16-m telescope .    in recent years , in order to improve the efficiency of the telescope and to maximize the scientific output , 7 - 8 key projects with the telescope have been supported for 3 - 5 years .
each project usually owns @xmath184 nights per year , which allows to carry out long - term observations and large - sample surveys . for instance , a number of quasars at intermediate redshift @xmath185 have been identified by @xcite .
an h@xmath186 imaging survey of @xmath121400 nearby alfalfa galaxies is carried out to study the star formation rate and the stellar distributions . from 1997 to 2002 , the spectra of @xmath182 blue compact galaxies ( bcgs ) , which was the largest sample before sdss , have been observed with the 2.16-m telescope , and the metallicity , extinction of dust , star formation rate of the sample galaxies have been derived by @xcite .
a group of peking university has discovered a high redshift quasar ( @xmath187 ) , j2202 + 1509 , in november 2014 with the 2.16-m telescope , which is sdss @xmath188 mag ( the paper is in preparation ) .
the redshift is derived from the emission lines , such as ly@xmath186 at wavelength of 1216  , and nv at wavelength of 1240  .
this demonstrates that the 2.16-m telescope is able to be used to discover and study high redshift quasars .
each year , the call for proposals begins around 20 october for the period of one month .
astronomers who are willing to use the telescope can submit proposals through the website http://astrocloud.china-vo.org before the deadline .
after that , the proposals will be collected and reviewed by the time allocation committee ( tac ) of the 2.16 m telescope .
the probability of observing nights obtained for observers highly depends on the mark ranking .
the too ( target - of - opportunity ) follow - up observations , like transient , sne and grbs , is supported , which are allowed to be applied and observed when the transient sources happen .
the turn - around time ( i.e. the time interval between a too alert and the start of the first observation ) of the 2.16-m telescope depends on the system and configure being used when it happens .
it concludes the following parts : the read - out time of last frame ( @xmath1230 - 50 seconds ) ; the time for changing the grisms and adjusting slit width and sometimes it even takes longer time when the grism is not on the spectragraph ( @xmath125 - 10 minutes ) ; the pointing time of the telescope ( @xmath121 - 5 minutes ) ; the possible time for changing instrument and focusing ( @xmath125 - 10 minutes ) , etc .
therefore , it takes @xmath1210 - 30 minutes in all for starting the too observations .
since the first light of the 2.16-m telescope , there are astronomers from france , japan , taiwan , hongkong , etc . , beyond china mainland , applied and used the telescope and publish papers .
until now , there are more than 150 sci papers published with the data obtained with the 2.16-m telescope since the first light .
all these papers have been peer - reviewed and most of them were published on the high - impact magazines in astrophysics , e.g. , @xmath189 , @xmath190 , @xmath191 and @xmath192 , and the impact factors are around @xmath193 , including the articles published in @xmath194 .
when the key projects are concluded , there will be numerous high impact journals published in the future .
as a 2-meter class optical astronomical telescope in china , the 2.16-m telescope of xinglong observatory of naoc , plays an important role in observational astronomy today .
there are currently three primary instruments , the bfosc , omr and the fiber - fed hrs , available for the telescope users .
when the 2.16-m telescope saw the first light in 1989 , various scientific projects have carried out , based on observations of the telescope in the research areas of nearby galaxies , agn , supernovae , grbs , exoplanets , galactic stars , time - domain astronomy , and many other sciences .
since the lamost spectroscopic survey started , the 2.16-m telescope has also been used for the lamost follow - ups , in order to identify interesting objects , and for stellar library observations .
a great many remarkable studies have been done with the telescope , including the spectroscopic research of the supernova 1993j @xcite . in recent years , in order to improve the utilizing efficiency of the telescope and the scientific output further , 7 - 8 key observing projects have been set up , in which the observers own 20 - 40 nights / year to carry our long - term observations and serial research .    in the future , a number of potential new instruments will be available to the 2.16-m telescope users to improve its observing abilities further , such as an imaging photopolarimeter @xmath195 ( please see , * ? ? ?
* ) , a prototype of which is currently in use at the 1.3-m telescope of the skinakas observatory in crete , greece , where it works well .
a similar photo - polarimeter will probably be installed at the 2.16-m telescope in the next one or two years , which will be the first astronomical photo - polarimeter in china .
it could be used for polarization measurements , for instance of seyfert galaxies , blazars , and grbs .
further , an intermediate resolution spectrograph ( irs ) is being investigated .
if it is installed and committed in the future , the projects of intermediate - resolution spectroscopic observations ( e.g. , x - ray binaries , bright stars , nearby galaxies , and a large number of fainter galactic stars ) could be carried out .
in addition , the testing of the astro - frequency comb calibration system is at the commissioning stage .
when it is finished , the precision of the radial velocity measurements of galactic stars could be improved greatly , even to a few cm s@xmath11 . recently a number of adaptive optics experiments were performed on the 2.16-m telescope by engineers from various chinese institutes , which can improve the spatial resolution and signal - to - noise ratio of targets to a large extent .
the exoplanet detection technology group from the nanjing institute of astronomical optics & technology , naoc and the california state university , northridge collaborate and built the high performance portable adaptive optics ( hpao ) which is mounted on the coud focus of 2.16-m telescope and succeed after the testing observations .
although the limiting magnitude of the current hpao system is faint to 3.8 mag in @xmath196 band , the system will be upgraded in the next step and the results can be improved significantly .
this research was supported by the national natural science foundation of china ( nfsc ) through grants 11003021 , 11373003 , 11273027 , 11303042 and national key basic research program of china ( 973 program ) 2015cb857002 .
z.f . acknowledges a young researcher grant of the national astronomical observatories , chinese academy of sciences .
we thank yuqin chen , liang wang and xiaoying yang for providing the limiting magnitude measurements of the hrs instrument and thank yingwei chen for providing the picture of the 2.16-m telescope .
chen , j. s. , fan , x. m. , tan , x. y.  1987 , acasn , 28 , 303 cui , xiang - qun , zhao , yong - heng , chu , yao - quan , et al .
2012 , raa , 12 , 1197 fan , y .- f . , et al .
2015 , raa , 15 , 918 huang , l. , wei , j .- y . ,
jiang , x .-
j . , lu , x .-
zhou , z .- z .
, li , h .- b .
, 2015 , astronomical research and technology , 12 , 44 huang , l. , wu , h. , & li , h. b. `` bfosc ( beijing faint object spectrograph and camera ) operating manual '' , 2012 , april . king , o. g. , blinov , d. , ramaprakash , a. n. , myserlis , i. , angelakis , e. , balokovic , m. , et al .
, 2014 , , 442 , 1706 kong , x. , cheng , f. z.  2002a , , 389 , 845 kong , x. , cheng , f. z. , weiss , a. , charlot , s. x. , cheng , f. z. , weiss , a. , charlot , s.  2002b , , 396 , 503 kong , x.  2004 , , 425 , 417 li , de - pei  2001 , optical instruments , april , 23 , 2 liu , y .- j . ,
sato , b. , zhao , g. , noguchi , k. , wang , h. , kambe , e. , ando , h. , izumiura , h. , chen , y .- q .
, okada , n. , et al .
2008 , , 672 , 553 liu , y .- j . ,
sato , b. , zhao , g. , ando , h.  2009 , raa , 9 , 1 shi , f. , kong , x. , li , c. , cheng , f .- z .
2005 , , 437 , 849 su , d .- q . , zhou , b .- f . , yu , x .-
1989 , science in china ( a ) , 11 , 1187 wang , l. , sato , b. , zhao , g. , liu , y .- j . ,
noguchi , k. , ando , h. , et al .
2012 , raa , 12 , 84 wang , l. , sato , b. , omiya , m. , harakawa , h. , liu , y .-
j . , et al .
2014 , , 66 , 118 wang , l .- f . ,
& hu , j .- y .  1994 , , 369 , 380 wilken , t. , lovis , c. , manescau , a. , steinmetz , t. , pasquini , l. , lo curto , g. , hansch , t. w. , holzwarth , r. , udem , th .  2010 , , 405 , 16 wu , x .- b . ,
wang , r. , schmidt , k. b. , bian , f .- y . , jiang , l .- h . , fan , x .- h .
2011 , , 142 , 78 wu , x .- b . , zuo , w .- w . , yang , j .- y . , yang , q. , wang , f.  2011 , , 142 , 78 zhang , j .- c . , ge . , l. , lu , x .-
cao , z .- h . ,
chen , x. , mao , y .-
n . , jiang , x .- j . , , 2015 , revised version submitted .
zhao , g. , & li , h .- b .
2001 , chjaa , 1 , 555 .
zhao , g. , qiu , h .-
m . , chen , y .- q . , li , z .- w . , 2000 , , 126 , 461 .
zheng , w .- k . ,
wang , j. , deng , j .- s . , hu , j .- y . , wei , j .- y . , 2009 , astronomical reserach & technology , 6 , 2 zhou , z .- m . , wu , h. , huang , l. , li , h .- b . ,
zhou , z .- z .
, jia , j .- j . ,
zhu , y .- n . , 2014 , raa , 14 , 1393


in many applications one faces the challenge to model objects , or parts of objects , in a mathematical framework . as an example , one important task is to extract an object from a given data set and manipulate it in a post - processing step in order to obtain further information .
typical applications include medical imaging , object tracking in a sequence of images , but also object recognition , where the post - processing step consists of the comparison of the extracted object with a database of reference objects .
similarly , such a comparison can be necessary in medical imaging in order to distinguish between healthy and diseased organs . to that end , however , one has to be able to decide whether two given objects are similar or not .
this requires a representation of the objects that makes the application of standard similarity measures possible .    finding a suitable representation of the object of interest , depending on the type of application , is crucial as a first step . for simplicity
, it is often assumed that the object is a simply connected bounded domain , allowing for the identification of the domain with its boundary . from a mathematical point of view , this assumption reduces the complexity of the representation . in addition , there exists a larger number of descriptions of boundaries than of domains , and , consequently , more mathematical tools to analyze the geometry of the underlying objects .    in 2d
a common approach is to encode the contour of an object by the curvature function of its boundary curve .
this approach has , for instance , been used in @xcite , where the authors set up a shape space of planar curves , where the shapes are implicitly encoded by the curvature function .
the main advantage of using a differential invariant  the most prominent representative being the curvature  to represent an object is the well investigated mathematical framework of this type of invariants ( see @xcite ) .
since all kinds of differential invariants are based on derivatives , they suffer from the shortcoming of being sensitive with respect to small perturbations . to bypass this shortcoming manay et al .
@xcite proposed to use integral invariants instead of their differential counterparts ( see also @xcite ) .
integral invariants have similar invariance properties as differential invariants , but have proven to be considerably more robust with respect to noise . their theory ,
however , is not that well investigated as opposed to the theory of their differential counterparts .    beside the classical approach of differential invariants and the novel approach of integral invariants , there exist several other concepts for encoding an object .
for instance , in @xcite the authors use the zero level set of a harmonic function , which is uniquely determined by prescribing two functions on the boundary of an annulus , to encode the boundary of a 2d object ( see @xcite for a generalization to compact surfaces in 3d ) .
a similar encoding of the object by a function is given in the article of sharon and mumford @xcite . here ,
the authors first map the 2d object , which is supposed to be a smooth and simply closed curve , to the interior of the unit disc in the complex plane via the riemann mapping theorem .
this conformal mapping is composed with a second one , generated out of the exterior of the original object , and the composition is restricted to the boundary of the unit disc .
thus , the final mapping , which the authors call the fingerprint of the object , is a diffeomorphism from the unit circle onto itself .
one of the challenges in object encoding is the question of uniqueness of the encoding .
more precisely , in many applications , e.g.  object matching , the correspondence between the object and its encoding should be one - to - one .
thus , a thorough investigation of the operator that maps an object to its encoding is needed . in case of the encoding by a harmonic or conformal mapping
 if possible  uniqueness is well known . also for the encoding of an arc length parameterized curve by its curvature function
, it is known that one obtains a one - to - one correspondence between the curve and its encoding ( up to rigid body motions ) .
one even has a complete characterization of the set of functions that arise as curvature functions of a class of sufficiently regular curves ( see @xcite ) . for integral invariants
the situation is different ; the cone area invariant , first introduced in @xcite , is an injective mapping independent of the space dimension , but its application is limited to star - shaped objects .
in contrast , for the circular integral invariant , which is the integral invariant most common in use , there exists no proof for the uniqueness conjecture so far .
this article is a contribution towards this goal : we first prove that the integral invariant is @xmath0-times continuously frchet differentiable in a neighborhood of the circle , seen as a mapping from @xmath1 to @xmath2 , @xmath3 .
then we show that the frchet differential is injective on some @xmath1-neighborhood of the circle , @xmath4 , @xmath5 .
the proofs of these results are based on the implicit function theorem on banach spaces and an application of riesz  schauder theory . using the injectivity result , a taylor series expansion , and an interpolation inequality for @xmath6 norms , we obtain the injectivity of the integral invariant on a @xmath7 neighborhood @xmath8 , @xmath4 .
more precisely , we show that , in @xmath8 , the @xmath2-norm of the difference of the integral invariants of two curves can be estimated from below by their @xmath2-distance .
let @xmath9 be the space of all continuous embeddings from @xmath10 to @xmath11 .
then every curve @xmath12 has a unique interior , denoted by @xmath13 .
following @xcite , this allows us to introduce the circular integral invariant :    for given @xmath14 we define the _ circular integral invariant _
@xmath15 \colon s^{1 } \to \mathbb{r}\ ] ] of a curve @xmath12 as @xmath16(\varphi ) : = \operatorname{area}\bigl(b_r\bigl(\gamma(\varphi)\bigr)\cap\operatorname{int}(\gamma)\bigr)\,,\ ] ] where @xmath17 denotes the ball of radius @xmath18 centered at @xmath19 .
the circular integral invariant behaves well under several group actions :    * @xmath20 is invariant with respect to euclidean motions : for @xmath21 $ ] we have @xmath22 = i_r[\gamma ] \;.\ ] ] * @xmath20 is equivariant with respect to reparametrizations : for every homeomorphism @xmath23 we have @xmath24 = i_r[\gamma]\circ\phi \;.\ ] ] * for every scalar @xmath25 we have @xmath26 = t^2
i_{r / t}[\gamma ] \;.\ ] ]    the observations above suggest to consider the integral invariant on the space @xmath27 of all curves modulo euclidean motions and reparametrizations .
moreover , we assume as an additional smoothness property that the considered curves are of class @xmath28 , @xmath29
. then it makes sense to use the following representation of @xmath27 , as it avoids working with equivalence classes of curves .
denote by @xmath30 , @xmath29 , the space of all curves @xmath31 satisfying the following conditions :    * @xmath32 has constant speed , i.e. , there exists a constant @xmath33 such that @xmath34 for all @xmath35 . *
@xmath36 and @xmath37 , where we identify the circle @xmath10 with the interval @xmath38 .
* @xmath32 is an embedding , that is , @xmath39 for all @xmath40 .    for the proof of our main theorem
we need the following result from differential geometry concerning the manifold structure of @xmath41 :    [ thm : smooth_submanifold ] for @xmath29 the space @xmath41 is a smooth submanifold of the banach space of all @xmath2-curves from @xmath10 to @xmath11 .
its tangent space @xmath42 at a curve @xmath43 consists of all @xmath2-curves @xmath44 with @xmath45    the proof of the submanifold result is similar to ( * ? ? ? * thm .
2.2 ) . in our case
the situation is less complicated , as we only deal with @xmath28-curves instead of sobolev curves .
the constant speed parameterization yields the condition @xmath46 the remaining constraints follow directly from the initial conditions .    under additional smoothness assumptions on @xmath32
we obtain the following characterization of the tangent space @xmath42 .
[ lem : tangent_space_characterization ] let @xmath4 and @xmath47 with curvature function @xmath48 then the tangent space @xmath49 consists of all @xmath2-curves @xmath50 satisfying :    * @xmath51 . * @xmath52 . * @xmath53 .
theorem [ thm : smooth_submanifold ] and the fact that @xmath54 imply that there exists a constant @xmath55 such that @xmath56 using the initial conditions for @xmath44 , we obtain the initial conditions for @xmath57 and @xmath58 and the value of @xmath59 .
we are now able to formulate the main result of this article . here and in the following
we denote by @xmath60 the @xmath2 norm on the space of curves .
similarly , if @xmath61 is a bounded linear operator , we denote its operator norm by @xmath62 , and we use the same notation for norms of multi - linear operators .
the circular integral invariant @xmath63 , @xmath29 , @xmath5 is @xmath0-times continuously frchet differentiable on a neighborhood @xmath64 of the circle of radius @xmath65 and its tangential mapping is injective on this neighborhood
.    moreover there exists a neighborhood @xmath66 of the circle with respect to the topology induced by the @xmath7-norm and a constant @xmath67 such that for every @xmath32 , @xmath68 the stability estimate @xmath69 - i_r[\tilde{\gamma}]\rvert_k    \ge c \lvert \gamma - \tilde{\gamma}\rvert_k\ ] ] holds . in particular ,
the mapping @xmath20 is injective on @xmath8 .
the condition on @xmath18 to be smaller than @xmath70 is necessary , because otherwise the circular integral invariant in each point @xmath71 is constant equal to @xmath72 , the area of the circle , and the same holds for any sufficiently small deformation of the circle which preserves the area .
in the following , we discuss the differentiability of @xmath20 and derive an analytic formula for @xmath20 and , under certain smoothness assumptions , its derivative @xmath73 valid in a neighborhood of the circle . as a first step , we recall the following result on the differentiability of the composition mapping . to that end , we need the following definitions of differentiability of mappings on banach spaces .
let @xmath74 , @xmath75 be banach spaces , @xmath76 , and @xmath77 .
the mapping @xmath78 is called @xmath0-times _ weakly differentiable _ , if it is @xmath0-times gteaux differentiable and its gteaux differential @xmath79 is continuous as a mapping @xmath80 in contrast , it is called @xmath0-times _ continuously frchet differentiable _ or _ of class @xmath81 _ , if it is @xmath0-times gteaux differentiable and @xmath79 is continuous as a mapping @xmath82 here @xmath83 is the banach space of @xmath0-linear mappings from @xmath84 to @xmath75 equipped with the operator norm .
note that the continuity requirement for weak differentiability is weaker than for continuous frchet differentiability , and thus a weakly differentiable mapping need not be frchet differentiable of the same order .
the following lemma shows , however , that it is frchet differentiable of lower order .    [ le : weakfrechet ] let @xmath74 , @xmath75 be banach spaces and @xmath85 @xmath86-times weakly differentiable with @xmath5
. then @xmath78 is @xmath0-times continuously frchet differentiable .
let @xmath87 and @xmath88 .
we have to show that there exists @xmath89 such that for every @xmath90 with @xmath91 the inequality @xmath92 holds .    using the continuity of @xmath93 at the point @xmath94 and the fact that @xmath95
, we obtain the existence of @xmath96 such that @xmath97 now let @xmath98 , let @xmath90 with @xmath99 and @xmath100 with @xmath101
. then we have @xmath102    [ composition ] for every @xmath103 , @xmath5 the composition mapping @xmath104 is @xmath86-times weakly differentiable and therefore @xmath0-times continuously frchet differentiable .
the result on weak differentiability has been shown in  ( * ? ? ?
* section  6.9 ) ( note that the result in @xcite has been shown for the space @xmath105 , which , however , is equivalent to @xmath106 in the case of a compact manifold ) .
the statement concerning the continuous frchet differentiability follows from lemma  [ le : weakfrechet ] .    in order to simplify the notation
, we will sometimes omit the domain and range of the function spaces in expressions like @xmath107 and write @xmath2 instead , if no confusion is possible .
we will need two different types of derivatives for the formulation of our results : first , frchet derivatives in the function space @xmath2 , and , second , derivatives of functions @xmath108 with respect to their argument @xmath109 . in order to highlight this difference ,
we use the following notation : for a function @xmath110 , we denote by @xmath111 its frchet derivative . in contrast , if @xmath112 , then @xmath113 denotes its derivative in the parameter space .    in order to make the notation less cumbersome , we omit the argument @xmath71 in @xmath114 , @xmath115 and similar expressions if the argument is clear from the context .
[ intersection_parameters ] for each @xmath103 , @xmath5 there exists a neighborhood @xmath116 of the constant speed parameterized circle of radius @xmath65 such that the following hold :    1 .   for each @xmath117 and @xmath35 the circle @xmath118 intersects the curve @xmath32 in exactly two points , denoted by @xmath119 and @xmath120 . here
@xmath121 denotes the previous intersection parameter and @xmath122 the next one ( see figure [ fig : sketch ] ) .
the mappings @xmath123 are @xmath0-times continuously frchet differentiable .
the first derivatives in direction @xmath124 are given by @xmath125 here @xmath126 denotes the group of @xmath28-diffeomorphisms on the unit circle .
denote by @xmath127 the constant speed parameterized circle , that is , @xmath128 .
then , for every @xmath35 , the circle @xmath129 intersects @xmath130 precisely at the two points @xmath131 , where @xmath132 now define the mapping @xmath133 obviously , the mappings @xmath134 and @xmath135 we are searching for satisfy @xmath136 and @xmath137 .
in particular , the equation @xmath138 has the two solutions @xmath139 and @xmath140 .
now , lemma  [ composition ] implies that the mapping @xmath141 , and consequently also @xmath78 , is @xmath0-times continuously frchet differentiable .
moreover , it is easy to see that the derivative of @xmath78 at @xmath142 in direction @xmath143 is given as @xmath144 which is obviously an isomorphism of @xmath107 , as the assumption @xmath145 implies that @xmath146 .
thus the implicit function theorem on banach spaces ( see  ( * ? ? ?
i.5 ) ) implies the existence of a neighborhood @xmath116 of @xmath130 and unique @xmath0-times continuously frchet differentiable mappings @xmath147 ,  @xmath148 satisfying the equations @xmath149 .
the formula for the directional derivative of @xmath150 at @xmath32 in direction @xmath44 now follows from the fact that @xmath151 which is a simple application of the chain rule .
the formula for @xmath152 can be derived analogously .
[ variation_general ] for each @xmath103 and @xmath153 there exists a neighborhood @xmath116 of the circle of radius @xmath65 such that the following hold :    1 .   for each @xmath117 the circular integral invariant @xmath154 $ ] can be written as @xmath155(\varphi )        =        \frac{1}{2 } \int_{m_{\gamma}}^{p_{\gamma } } \langle \gamma(\psi)-\gamma,\dot{\gamma}(\psi)^{\bot } \rangle d\psi        + \frac{r^2}{2 } \arccos \biggl(\frac{\langle \gamma(p_{\gamma})-\gamma,\gamma(m_{\gamma})-\gamma \rangle}{r^2 } \biggr ) \;.\ ] ] 2 .
the circular integral invariant @xmath156\ , ,        \end{aligned}\ ] ] is @xmath0-times continuously frchet differentiable .
its derivative in direction @xmath157 is given by @xmath158(\sigma)\\        & =        2 \int_{m_{\gamma}}^{p_{\gamma } } \langle \sigma(\psi),\dot{\gamma}(\psi)^{\bot } \rangle d\psi \\        & \qquad        - \langle \sigma,\gamma(p_{\gamma})^{\bot}-\gamma(m_{\gamma})^{\bot } \rangle        + \langle \gamma(p_{\gamma})-\gamma,\sigma(p_{\gamma})^{\bot } \rangle        - \langle \gamma(m_{\gamma})-\gamma,\sigma(m_{\gamma})^{\bot } \rangle   \\        & \qquad        + \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(p_{\gamma})^{\bot } \rangle        \frac{\langle \sigma-\sigma(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle}{\langle \dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle}\\        & \qquad        - \langle \gamma(m_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})^{\bot}\rangle        \frac{\langle \sigma-\sigma(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle}{\langle \dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle}\\        & \qquad        - \frac{r^2}{\sqrt{r^4 - \langle \gamma(p_{\gamma})-\gamma,\gamma(m_{\gamma})-\gamma\rangle^2}}\cdot\\        & \qquad \qquad        \cdot\bigg ( \ !        \frac{\langle
\sigma-\sigma(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle}{\langle \dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle }        \langle \dot{\gamma}(p_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle \\        & \qquad \qquad \qquad        + \langle \sigma(p_{\gamma})-\sigma,\gamma(m_{\gamma})-\gamma \rangle        + \langle \gamma(p_{\gamma})-\gamma,\sigma(m_{\gamma})-\sigma \rangle\\        & \qquad \qquad \qquad        + \frac{\langle \sigma-\sigma(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle}{\langle \dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle }        \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})\rangle        \ ! \bigg)\;.      \end{aligned}\ ] ]    under the given assumptions for @xmath159 , @xmath18 and @xmath8 , formula   can be easily deduced from figure  [ fig : sketch ] .
a term by term investigation of formula  , using lemma  [ composition ] and lemma  [ intersection_parameters ] , shows that @xmath20 is of class @xmath81 on @xmath8 .    to calculate the differential of @xmath20 we treat the two terms of formula   separately .
for the first term we obtain @xmath160 a simple application of the chain rule yields for the second term @xmath161 using the formulas for the intersection parameters , we obtain the desired result .
in the special case where @xmath32 equals the unit circle the lemma above reduces to :    [ formelkreis ] let @xmath162 , @xmath103 , @xmath5 , be the constant speed parameterized unit circle , that is , @xmath163 and let @xmath164 .
then the derivative of @xmath154 $ ] in direction @xmath124 with @xmath165 is given by @xmath166(\sigma)(\varphi )      =      \int_{\varphi-\vartheta}^{\varphi+\vartheta } a(\psi ) \
, d\psi - 2 \sin(\vartheta ) a(\varphi )      =      \big ( \chi_{[-\vartheta,\vartheta ] } \ast a \big)(\varphi ) - 2 \sin(\vartheta ) a(\varphi)\ ] ] with @xmath167    the proof of this lemma is postponed to the appendix .
we have already shown in section  [ se : diff ] that @xmath20 is of class @xmath81 .
we now show the local injectivity of @xmath73 . without loss of generality
we may assume that @xmath168 .
let @xmath8 be the neighborhood of the unit circle defined in theorem  [ variation_general ] .
the formula for @xmath73 implies that for every @xmath169 the mapping @xmath170 $ ] is bounded as a mapping from @xmath171 to @xmath107 .
thus , @xmath170 $ ] has a unique bounded extension @xmath172 \colon c^k(s^1;\mathbb{r}^2 ) \to c^k(s^1;\mathbb{r})$ ] .
moreover , the mapping @xmath173 is continuous with respect to the @xmath1-topology seen as a mapping from @xmath8 to @xmath174 .
in addition , we denote for @xmath117 by @xmath175 $ ] the restriction of @xmath172 $ ] to @xmath49 .
denote now by @xmath176 the constant speed parameterized unit circle .
define for given @xmath177 the function @xmath178 by @xmath179 because @xmath130 is a @xmath180-curve , it follows that @xmath181 is @xmath28 . using lemma [ lem : tangent_space_characterization ]
it follows that @xmath182 and @xmath183 .
define now the space @xmath184 then it follows that @xmath185 is a bounded linear mapping from @xmath186 to @xmath187 .
in addition , it follows from lemma [ lem : tangent_space_characterization ] that @xmath185 is boundedly invertible with @xmath188 given by @xmath189 with @xmath190 the expression for @xmath191 is due to the periodicity of @xmath58 , which implies that @xmath192 therefore @xmath185 is in fact an isomorphism between @xmath186 and @xmath187 .    according to lemma
[ formelkreis ] , the mapping @xmath193 $ ] evaluated at @xmath194 can be written as @xmath195(\sigma ) = \chi_{[-\vartheta,\vartheta ] } \ast a - 2 \sin(\vartheta ) a \ ; .\ ] ] thus @xmath193 $ ] can be decomposed into @xmath195 = b \circ \imath \circ a \;,\ ] ] where the operator @xmath196 is given by @xmath197 } \ast a - 2 \sin(\vartheta ) a\ ] ] and @xmath198 is the embedding from @xmath187 into @xmath199 .
lemma  [ le : compact ] ( see appendix ) implies that the mapping @xmath200}\ast a$ ] is compact and thus @xmath201 is a compact perturbation of the identity .
therefore the riesz  schauder theory ( see  ( * ? ? ?
x.5 ) ) implies that @xmath201 has a closed range .
next we compute the kernel of @xmath201 . to that end
we consider the mapping in the fourier basis .
a short calculation shows that in this basis the operator @xmath201 is the diagonal operator that maps a sequence of ( complex ) fourier coefficients @xmath202 to the sequence @xmath203 , where @xmath204 & \text { else . }
\end{cases}\ ] ] because @xmath205 and @xmath206 whenever @xmath207 ( see lemma [ lem : sine_inequality ] in the appendix ) , it follows that the kernel of @xmath201 consists of the functions @xmath57 of the form @xmath208 for some @xmath209 , @xmath210 .    in the next step
we show that the kernel of @xmath193 = b \circ \imath \circ a$ ] is trivial .
therefore assume that @xmath211 . because @xmath212 , it follows that @xmath213 ; because @xmath214 , it follows that @xmath215 .
together , this shows that @xmath216 , implying that the intersection of @xmath217 with @xmath187 is trivial .
since @xmath185 is an isomorphism this proves the injectivity of @xmath218 $ ] .
we have thus shown that @xmath218 = b \circ \imath \circ a$ ] is strongly closed and injective .
now note that the set of strongly closed and injective , bounded linear functionals between two banach spaces @xmath74 and @xmath75 is open with respect to the norm topology on @xmath219 . because of the continuity of @xmath173 and therefore @xmath220 , this proves the existence of a neighborhood @xmath221 such that @xmath175 $ ] is injective for every @xmath222 .
in particular , this proves the injectivity of @xmath170 $ ] for every @xmath223 .    for proving the local injectivity ,
let @xmath32 , @xmath224 .
because the mapping @xmath225 is of class @xmath226 , it has a taylor expansion of the form @xmath227 = i_r[\gamma ] + i_r'[\gamma](\tilde{\gamma}-\gamma )    + \int_0 ^ 1(1-t)i_r''[\gamma+t(\tilde{\gamma}-\gamma)](\tilde{\gamma}-\gamma,\tilde{\gamma}-\gamma)\,dt.\ ] ] thus @xmath228-i_r[\gamma]\rvert_k      & \ge \lvert i_r'[\gamma](\tilde{\gamma}-\gamma)\rvert_k      - \bigl\rvert\int_0 ^ 1(1-t)i_r''[\gamma+t(\tilde{\gamma}-\gamma)](\tilde{\gamma}-\gamma,\tilde{\gamma}-\gamma)\,dt\bigr\rvert_k\\      & \ge \lvert i_r'[\gamma](\tilde{\gamma}-\gamma)\rvert_k      - \int_0 ^ 1\bigl\rvert i_r''[\gamma+t(\tilde{\gamma}-\gamma)](\tilde{\gamma}-\gamma,\tilde{\gamma}-\gamma)\bigr\rvert_k\,dt    \end{aligned}\ ] ] because @xmath229 is of class @xmath226 , it follows that @xmath230 is continuous .
thus there exists a convex neighborhood @xmath231 of the circle with respect to the @xmath232-norm and a constant @xmath233 such that @xmath234\rvert_{l^2(c^{k+3},c^k ) } \le c_1 $ ] for every @xmath235 .
consequently , @xmath236-i_r[\gamma]\rvert_k      \ge \lvert i_r'[\gamma](\tilde{\gamma}-\gamma)\rvert_k - c_1 \lvert \tilde{\gamma}-\gamma\rvert_{k+3}^2\ ] ] for every @xmath32 , @xmath237 .
since @xmath238 is continuous , for every @xmath88 there exists a neighbourhood @xmath239 of the constant speed parameterized unit circle @xmath130 such that @xmath240 - i_r'[\gamma_0]\rvert < \varepsilon$ ] for every @xmath241 .
in particular , we have for @xmath241 @xmath242(\tilde{\gamma}-\gamma)\rvert_k    \ge \lvert i_r'[\gamma_0](\tilde{\gamma}-\gamma)\rvert_k    - \varepsilon\lvert \tilde{\gamma}-\gamma\rvert_k.\ ] ]    since @xmath41 is a smooth submanifold of @xmath243 , there exists a neighborhood @xmath244 of @xmath130 and a smooth diffeomorphism @xmath245 such that @xmath246 = 0 ,    \qquad
\phi(\mathcal{v}_2 \cap \mathcal{c}_k ) \subset t_{\gamma_0}\mathcal{c}_k ,    \quad\text { and } \quad    \phi'[\gamma_0 ] = \operatorname{id}.\ ] ] in order to show that such a map exists , let @xmath247 be any submanifold chart centered at @xmath130 .
that is , there exists a closed linear subspace @xmath248 such that @xmath249
. then the mapping @xmath250^{-1 } \circ \psi$ ] has the desired properties .
thus the continuous invertibility of @xmath251 $ ] on @xmath252 seen as a mapping from @xmath2 to @xmath2 implies that there exists @xmath253 such that for every @xmath32 , @xmath254 we have @xmath255(\tilde{\gamma}-\gamma)\rvert_k      & \ge \lvert i_r'[\gamma_0](\phi(\tilde{\gamma})-\phi(\gamma))\rvert_k-      \lvert i_r'[\gamma_0](\tilde{\gamma}-\gamma-\phi(\tilde{\gamma})+\phi(\gamma))\rvert_k\\      & \ge c_2 \lvert \phi(\tilde{\gamma})-\phi(\gamma)\rvert_k      - c_3 \lvert \tilde{\gamma}-\gamma-\phi(\tilde{\gamma})+\phi(\gamma)\rvert_k    \end{aligned}\ ] ] with @xmath256\rvert_{k+3,k}$ ] .
developing @xmath257 in a taylor expansion at @xmath32 , we obtain after possibly choosing a smaller neighbourhood @xmath258(\tilde{\gamma}-\gamma ) + r(\gamma,\tilde{\gamma})\ ] ] with @xmath259 for some @xmath260 independent of @xmath32 . inserting this taylor expansion into   yields
@xmath261(\tilde{\gamma}-\gamma)\rvert_k      & \ge c_2\lvert \phi'[\gamma](\tilde{\gamma}-\gamma)\rvert_k -         c_3\lvert ( \operatorname{id}-\phi'[\gamma])(\tilde{\gamma}-\gamma)\rvert_k\\      & \qquad\qquad        - ( c_2+c_3)\lvert r(\gamma,\tilde{\gamma})\rvert_k\\      & \ge c_2\lvert \tilde{\gamma}-\gamma\rvert_k        - ( c_2+c_3)\lvert ( \operatorname{id}-\phi'[\gamma])(\tilde{\gamma}-\gamma)\rvert_k\\      & \qquad\qquad        - c_4(c_2+c_3)\lvert \gamma-\tilde{\gamma}\rvert_k^2\\      & \ge \bigl(c_2- c_5\lvert \operatorname{id}-\phi'[\gamma]\rvert_{k , k } - c_6\lvert \gamma-\tilde{\gamma}\rvert_k\bigr )        \lvert \gamma-\tilde{\gamma}\rvert_k .
\end{aligned}\ ] ] because @xmath262 is smooth and @xmath263 = \operatorname{id}$ ] , it follows that there exists a neighbourhood @xmath264 of @xmath130 such that @xmath265(\tilde{\gamma}-\gamma)\rvert_k \ge c_7 \lvert \gamma-\tilde{\gamma}\rvert_k\ ] ] for every @xmath32 , @xmath266 .
collecting the inequalities  , , and  , we obtain @xmath267-i_r[\gamma]\rvert_k    \ge ( c_7-\varepsilon)\lvert \gamma-\tilde{\gamma}\rvert_k - c_1\lvert\gamma-\tilde{\gamma}\rvert_{k+3}^2\ ] ] for every @xmath32 , @xmath268 .    in order to obtain the desired result
, we use the interpolation inequality ( see  ( * ? ? ?
* theorem  2.2.1 , p.  143 ) ) @xmath269    choosing @xmath88 in   sufficiently small , we obtain the estimate @xmath270-i_r[\gamma]\rvert_k    \ge ( c_9-c_8\lvert\gamma-\tilde{\gamma}\rvert_{k+6})\lvert\gamma-\tilde{\gamma}\rvert_k\ ] ] for every @xmath32 , @xmath271 .
thus there exists a neighbourhood @xmath272 of the circle @xmath130 and a constant @xmath273 such that @xmath270-i_r[\gamma]\rvert_k    \ge c\lvert\gamma-\tilde{\gamma}\rvert_k\ ] ] for every @xmath32 , @xmath68 .
in this article , we have shown an injectivity result for the circular integral invariant on a @xmath7 neighborhood @xmath66 of the circle .
note , however , that the derived result does not prove the continuous invertibility of the invariant on @xmath66
.    the classical approach for proving such a result would be the usage of the inverse function theorem . to that end , however , we would require that the mapping @xmath20 was continuously frchet differentiable and its derivative @xmath73 an isomorphism of the corresponding tangent spaces .
although our results prove that @xmath73 can be extended to an isomorphism @xmath220 of @xmath252 , we can not use the inverse function theorem , as the mapping @xmath20 is not frchet differentiable ( and not even gteaux differentiable ) from @xmath2 to @xmath2  in fact , our results do not even prove that @xmath20 maps @xmath2 curves into @xmath2 integral invariants .
conversely , seen as a mapping from @xmath274 to @xmath2 , the tangent mapping , though injective , can not be a surjection near the circle .
another issue are the rather stringent smoothness assumptions .
we believe that it is possible to relax these assumptions , as it seems probable that @xmath20 is twice weakly differentiable as a mapping from @xmath274 to @xmath2 , which would indicate a possible uniqueness result on a @xmath275-neighborhood of the circle .
the authors want to thank peter michor and gnther hrmann for their comments and suggestions , which helped to improve the article .
we thank the referee for the careful reading of the article and for pointing out a mistake in the proof of the original main theorem .
[ le : compact ] for every @xmath3 , the mapping @xmath276 } \ast a$ ] is compact as a mapping from @xmath199 to @xmath199 .
denoting by @xmath277 the unit ball in @xmath199 , we have to show that the image of @xmath277 under @xmath278 is precompact in @xmath199 . applying the arzel  ascoli theorem ( see  ( * ? ? ?
iii.3 ) ) , we have to show that @xmath279 is bounded and the first @xmath280 derivatives of the functions in @xmath279 are equicontinuous .
the boundedness of @xmath279 is obvious , @xmath278 being a bounded linear mapping ( of norm @xmath281 ) .
now assume that @xmath282 .
then @xmath283 . because @xmath284 for all @xmath285 , it follows that @xmath286 is lipschitz continuous with lipschitz constant at most 2 .
hence @xmath279 is a precompact set .
assume first that @xmath289 .
we show that in this case the equation @xmath290 has the only solutions @xmath291 and @xmath292 .
first note that the strict concavity of the sine function on the interval @xmath293 $ ] implies that on this interval we only have two solutions , namely @xmath294 and @xmath295 .
moreover , the concavity of the sine implies that @xmath296 , and therefore @xmath297 .
this , however , implies that the equation @xmath298 can not have any solutions for @xmath299 , as the right hand side is strictly larger than @xmath300 . the fact that @xmath301 is the only negative solution follows by symmetry . in particular , setting @xmath302
, this proves the assertion in the case @xmath289 .
now assume that @xmath303 and let @xmath304 .
then @xmath305 now , if @xmath306 is odd , then @xmath307 thus @xmath308 , if and only if @xmath309 . because @xmath310 , the first part of the proof can be applied , showing that @xmath311 unless @xmath312 .
on the other hand , if @xmath306 is even , we have @xmath313 thus , @xmath308 , if and only if @xmath314 .
now note that the equation @xmath315 has only the trivial solution @xmath291 , because @xmath316 and the left hand side is negative for @xmath317 . as a consequence , the equation @xmath314 only holds for @xmath318 , which concludes the proof .
let @xmath319 be the constant speed parameterized unit circle .
then @xmath320 because @xmath32 is the unit circle , there exists @xmath321 such that @xmath322 obviously , all assumptions of lemma  [ variation_general ] are satisfied . it remains to calculate all the terms that appear in the expression of the frchet derivative in lemma  [ variation_general ] for the special case of the unit circle .
in particular , we obtain @xmath323 in the following we treat the derivative in direction @xmath324 and @xmath325 separately . for the derivative in normal direction
@xmath324 we get @xmath326(a\dot{\gamma}^{\bot } )      =      2 \int_{m_{\gamma}}^{p_{\gamma } } a(\psi ) d\psi - a \langle \dot{\gamma}^{\bot},\gamma(p_{\gamma})^{\bot}-\gamma(m_{\gamma})^{\bot } \rangle\\      & \qquad      - a(p_{\gamma } ) \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(p_{\gamma } ) \rangle      + a(m_{\gamma } ) \langle \gamma(m_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})\rangle   \\      & \qquad      + \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(p_{\gamma})^{\bot } \rangle      \frac{\langle a\dot{\gamma}^{\bot}-a(p_{\gamma})\dot{\gamma}(p_{\gamma})^{\bot } , \gamma(p_{\gamma})-\gamma \rangle }           { \langle \dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle}\\           & \qquad           - \langle \gamma(m_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})^{\bot } \rangle           \frac{\langle a\dot{\gamma}^{\bot}-a(m_{\gamma})\dot{\gamma}(m_{\gamma})^{\bot},\gamma(m_{\gamma})-\gamma \rangle }                { \langle \dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle}\\                & \qquad                - \frac{r^2}{\sqrt{r^4 - \langle \gamma(p_{\gamma})-\gamma,\gamma(m_{\gamma})-\gamma \rangle^2}}\cdot\\                & \qquad\quad \cdot                \bigg ( \ !                \frac{\langle a\dot{\gamma}^{\bot}-a(p_{\gamma})\dot{\gamma}(p_{\gamma})^{\bot},\gamma(p_{\gamma})-\gamma \rangle }                     { \langle \dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle }                     \langle \dot{\gamma}(p_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle \\                     & \qquad \qquad                     + \langle a(p_{\gamma})\dot{\gamma}(p_{\gamma})^{\bot}-a\dot{\gamma}^{\bot},\gamma(m_{\gamma})-\gamma \rangle                     + \langle \gamma(p_{\gamma})-\gamma , a(m_{\gamma})\dot{\gamma}(m_{\gamma})^{\bot}-a\dot{\gamma}^{\bot } \rangle\\                     & \qquad \qquad                     + \frac{\langle a\dot{\gamma}^{\bot}-a(m_{\gamma})\dot{\gamma}(m_{\gamma})^{\bot},\gamma(m_{\gamma})-\gamma \rangle }                     { \langle \dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle }                     \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})\rangle                     \ ! \bigg)\;.\\    \end{aligned}\ ] ] inserting the expressions for @xmath32 equal to the unit circle that have been calculated previously , we obtain @xmath326(a\dot{\gamma}^{\bot})\\      & = 2 \int_{\varphi-\vartheta}^{\varphi+\vartheta } a(\psi ) d\psi \\      & \qquad      - \sin(\vartheta ) \big ( a(p_{\gamma } ) + 2a + a(m_{\gamma } ) \big )      - \frac{(\cos(\vartheta)-1)^2 } { \sin(\vartheta ) } \big ( a(p_{\gamma } ) + 2a + a(m_{\gamma } ) \big)\\      & \qquad      -\frac{1}{\sin(\vartheta ) }      \big ( \ !
\big ( a + a(p_{\gamma } ) \big ) \big ( \ ! \cos(\vartheta)-1 \big ) \big ( 1 - 2\cos(\vartheta ) \big)-2a(\cos(\vartheta)-1)\\      & \qquad \qquad      + \big ( a(p_{\gamma } ) + a(m_{\gamma } ) \big ) \big ( \ !
\cos^2(\vartheta)-\cos(\vartheta)-\sin^2(\vartheta ) \big)\\      & \qquad \qquad      + \big ( a + a(m_{\gamma } ) \big ) \big ( \ ! \cos(\vartheta)-1
\big ) \big ( 1 - 2\cos(\vartheta ) \big )      \ ! \big)\\      & = 2 \int_{\varphi-\vartheta}^{\varphi+\vartheta } a(\psi ) d\psi \\      & \qquad      + \frac{\cos(\vartheta)-1}{\sin(\vartheta ) } \big ( 2a(p_{\gamma } ) + 4a + 2a(m_{\gamma } ) - 2a(p_{\gamma } ) + 4\cos(\vartheta)a - 2a(m_{\gamma } ) \big)\\      & = 2 \int_{\varphi-\vartheta}^{\varphi+\vartheta } a(\psi ) d\psi      + 4a \frac{(\cos(\vartheta)-1)(\cos(\vartheta)+1)}{\sin(\vartheta)}\\      & = 2 \int_{\varphi-\vartheta}^{\varphi+\vartheta } a(\psi ) d\psi - 4a\sin(\vartheta ) \ ; .
\end{aligned}\ ] ] inserting the formulas above in the expression for @xmath327(b \dot{\gamma})$ ] yields @xmath328(b\dot{\gamma } )      =      2 \int_{m_{\gamma}}^{p_{\gamma } } b(\psi ) \langle \dot{\gamma}(\psi),\dot{\gamma}(\psi)^{\bot } \rangle d\psi \\      & \quad      - \langle b \dot{\gamma},\gamma(p_{\gamma})^{\bot}-\gamma(m_{\gamma})^{\bot } \rangle      + \langle \gamma(p_{\gamma})-\gamma , b(p_{\gamma})\dot{\gamma}(p_{\gamma})^{\bot } \rangle\\      & \quad      - \langle \gamma(m_{\gamma})-\gamma , b(m_{\gamma})\dot{\gamma}(m_{\gamma})^{\bot } \rangle      + \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(p_{\gamma})^{\bot }
\rangle      \frac{\langle b\dot{\gamma}-b(p_{\gamma})\dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle }       { \langle \dot{\gamma}(p_{\gamma } ) , \gamma(p_{\gamma})-\gamma \rangle}\\      & \quad       - \langle \gamma(m_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma})^{\bot } \rangle       \frac{\langle b\dot{\gamma}-b(m_{\gamma})\dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle }            { \langle \dot{\gamma}(m_{\gamma } ) , \gamma(m_{\gamma})-\gamma \rangle}\\            & \quad            - \frac{r^2}{\sqrt{r^4 - \langle \gamma(p_{\gamma})-\gamma,\gamma(m_{\gamma})-\gamma \rangle^2}}\cdot\\                & \qquad\quad            \cdot\bigg ( \ !            \frac{\langle b\dot{\gamma}-b(p_{\gamma})\dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle }             { \langle \dot{\gamma}(p_{\gamma}),\gamma(p_{\gamma})-\gamma \rangle }             \langle \dot{\gamma}(p_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle \\             & \qquad \qquad             + \langle b(p_{\gamma})\dot{\gamma}(p_{\gamma})-b\dot{\gamma},\gamma(m_{\gamma})-\gamma \rangle             + \langle \gamma(p_{\gamma})-\gamma , b(m_{\gamma})\dot{\gamma}(m_{\gamma})-b\dot{\gamma } \rangle\\             & \qquad \qquad             + \frac{\langle b\dot{\gamma}-b(m_{\gamma})\dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle }             { \langle \dot{\gamma}(m_{\gamma}),\gamma(m_{\gamma})-\gamma \rangle }             \langle \gamma(p_{\gamma})-\gamma,\dot{\gamma}(m_{\gamma } ) \rangle             \ ! \bigg)\\             & = 0 - 0 + \big ( 1-\cos(\vartheta ) \big ) \big ( b(p_{\gamma})-b(m_{\gamma})+b - b(p_{\gamma})-b+b(m_{\gamma } ) \big)\\             & \quad -\big ( \ !
\big ( 1 - 2\cos(\vartheta)\big ) \big ( b - b(p_{\gamma})+b(p_{\gamma})-b(m_{\gamma})+b(m_{\gamma})-b \big ) + ( b - b ) \ ! \big)\\             & = 0 \ ; .
\end{aligned}\ ] ] therefore , @xmath166(a\dot{\gamma}^{\bot } + b\dot{\gamma } ) = \chi_{[-\vartheta,\vartheta ] } \ast a - 2\sin(\vartheta ) a\;.      \qedhere\ ] ]          a.  duci , a.  j. yezzi , jr .
, s.  k. mitter , and s.  soatto .
shape representation via harmonic embedding . in _ computer vision , 2003 .
ninth ieee international conference on _ , volume  1 , pages 656662 .
ieee , october 2003 .
huang , s.  flry , n.  gelfand , m.  hofer , and h.  pottmann .
reassembling fractured objects by geometric matching . in _
acm siggraph 2006 papers _ , siggraph 06 , pages 569578 , new york , ny , usa , 2006 .
acm .
s.  manay , b .- w .
hong , a.  j. yezzi , jr . , and s.  soatto .
integral invariant signatures . in t.
pajdla and j.  matas , editors , _ computer vision - eccv 2004 _ , volume 3024 of _ lecture notes in computer science _ , pages 8799 .
springer berlin / heidelberg , 2004 .
w. michor , some geometric evolution equations arising as geodesic equations on groups of diffeomorphisms including the hamiltonian approach .
in a.  bove , f.  colombini , d.  del santo , editors , _ phase space analysis of partial differential equations _ ,
volume 69 of _ progr .
nonlinear differential equations appl .
_ , pages 133215 .
birkhuser boston ( ma ) , 2006 .
yang , y .- k .
hu , and h.  pottmann .
robust principal curvatures on multiple scales . in k.
polthier and a.  sheffer , editors , _ proceedings of the fourth eurographics symposium on geometry processing _ ,
pages 223226 , aire - la - ville , switzerland , switzerland , 2006 .
eurographics association .