Instruction: Write an abstract for the following acadamic paper.

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

Paper: the properties of quasi - one - dimensional conductors have been subject of intensive studies for half a century . a theory based on the one - dimensional ( 1d ) hubbard model predicts a non - fermi - liquid metallic state.@xcite the situation , however , develops differently when the filling of the system is commensurate , since in this case , electronic interactions may cause a mott insulating state . although this phenomenon is possible in all dimensions , it is particularly strong in the 1d case via umklapp process.@xcite the influence of the interchain coupling on the electronic properties of a half- and quarter - filled quasi-1d interacting electron system was studied theoretically.@xcite it is predicted@xcite that with increasing interchain transfer integral @xmath6 such a system undergoes a deconfinement transition into a metallic state . the schematic phase diagram illustrating this phenomenon is shown in fig . [ fig : deconfine ] . strictly speaking , the deconfinement transition takes place at @xmath7 when @xmath6 becomes larger than a certain value @xmath8 and can be referred to as a quantum phase transition . for @xmath9 the system is in a mott insulating state and the single - particle interchain hopping is completely suppressed . if the interchain transfer integral exceeds a critical value @xmath8 , the insulating state is broken and single - particle interchain hopping becomes possible . thus , the system evolves from a 1d mott - hubbard insulator to a high - dimensional metallic state ( trajectory 1 in fig . [ fig : deconfine ] ) . the theoretical study@xcite of two coupled hubbard chains predicts a close relation of the interchain transfer integral to the charge gap , @xmath10.@xcite while in a simple two - chain system no real deconfinement is observed,@xcite for a large ( infinite ) number of chains@xcite the crossover is expected at some critical value @xmath11.@xcite     when @xmath12 reaches a critical value @xmath13 ( trajectory 1 ) . at high enough temperature , the increase of @xmath12 leads to a transition from a mott insulating to a 1d luttinger liquid ( ll ) conducting state ( trajectory 2 ) and finally to a dimensional crossover into an hdm state ( trajectory 3 ) . the temperature change at a fixed value @xmath12 may also lead to a transition from a mott insulating to a ll state ( trajectory 4 ) , or to a dimensional crossover ( trajectory 5 ) . the horizontal dashed line indicates the path of our pressure study.[fig : deconfine],title="fig : " ] +    if the temperature or energy scale at which the response of the system is probed becomes larger than @xmath12 , the warping of the open fermi surface is masked . therefore , the properties of the system are virtually the same as in the case of a flat fermi surface , i.e. , for a strictly 1d system where the interacting electrons are described by the luttinger liquid ( ll ) model.@xcite thus , the deconfinement transition at high enough temperature leads to a crossover from a mott insulating to a ll metallic state ( see trajectories 2 and 4 in fig .  [ fig : deconfine ] ) . furthermore , a dimensional crossover from a ll to a high - dimensional metallic state occurs when the interchain coupling ( i.e. , the warping of the fermi surface ) becomes larger than the energy of thermal fluctuations ( @xmath14 ) . this can be achieved by increasing the transfer integral @xmath15 ( trajectory 3 ) or by lowering the temperature ( trajectory 5 ) . in the high - dimensional metallic state the fermi - liquid description is valid instead of the ll model . suitable systems to study the predicted dimensional crossover and deconfinement transition are the well - known organic bechgaard - fabre salts@xcite ( tmttf)@xmath16 and ( tmtsf)@xmath16 . they consist of weakly coupled molecular stacks along which the charge transport occurs . hereby , the type of the molecule ( tmtsf or tmttf ) and of the anions @xmath17 determine the coupling between the molecular stacks and thus the dimensionality of the system ( chemical pressure effect ) . in particular , the ( tmttf)@xmath16 salts consist of basically uncoupled chains , which are half - filled due to the dimerization along the chains . thus , they constitute systems of electrons confined to the chains , forming a mott insulating state , and are located on the left side of the diagram given in fig .  [ fig : deconfine ] . at low temperatures the mott insulating state of ( tmttf)@xmath16 salts is transformed into a charge - ordered state characterized by a disproportion of charge of the donor tmttf molecules within the stacks . this ferroelectric transition is strongly pressure - dependent , and the charge - ordered phase can be suppressed by the application of a moderate pressure of a few kbars.@xcite on the other hand , the compounds ( tmtsf)@xmath16 represent examples of quasi-1d metals and they are on the right side of the diagram , since the enhanced interchain coupling leads to a charge deconfinement .    the low - temperature optical conductivity spectra of the tmtsf salts along the stacks are characterized by two distinct features : a narrow mode at zero energy containing a very small part of the spectral weight and a high - energy mode centered at around 200  @xmath18.@xcite the latter feature was interpreted in terms of the mott ( charge ) gap of a 1d insulator.@xcite the frequency dependence of the optical conductivity in the energy range above the effective interchain transfer integral [ about 250  @xmath18  for @xmath19 is well described by a ll model , since the system can be considered as effectively one - dimensional . on the other hand , the zero - energy mode was suggested to originate from the effective doping of the 1d mott insulator due to the deconfined interchain electron hopping . thus , it is related to 2d physics and could be successfully described by the fermi - liquid model.@xcite similar on - chain optical conductivity spectra with the dominating mott gap were also observed for tmttf compounds.@xcite however , the zero - energy mode typical for tmtsf salts is absent , since these systems are in the confined state where interchain hopping is suppressed . thus , optical spectroscopy on the high - energy range ( larger than the interchain transfer integral ) addresses electronic interactions in 1d , while dc resistivity measurements are more susceptible to interchain hopping and thus reveal deviations from a strictly 1d picture . first hints for the deconfinement transition and the dimensional crossover were found by pressure - dependent transport measurements on ( tmttf)@xmath20pf@xmath21 and @xmath22  crystals.@xcite the effect of the dimensional crossover occurring during the temperature decrease was addressed by an optical study of @xmath22.@xcite lateron , the above - mentioned theoretical criterium for the deconfinement transition was verified by comparing the ambient - pressure optical spectra of different bechgaard - fabre salts,@xcite i.e. , making use of the effect of chemical pressure on the interchain coupling . however , this approach has disadvantages , since only a limited number of discrete pressure points corresponding to particular compounds are available . moreover , atomic substitution ( especially for sulfur and selenium atoms ) may appreciably affect not only the interchain transfer integral @xmath12 , but also the stack dimerization , on - chain bandwidth , coulomb repulsion etc .    a more direct and cleaner way to study the deconfinement transition experimentally is the application of external pressure , since the changes in the electrodynamic response can be monitored while the interchain coupling is being continuously tuned . in this paper we present the results of the first pressure - dependent infrared spectroscopic study on the salt @xmath0 which is a prime example of a 1d mott - hubbard insulator at ambient pressure.@xcite the main aim of the present investigation is to determine the pressure - induced changes in the optical response of @xmath0  during the deconfinement transition and the dimensional crossover . during our experiment we thus intended to follow the evolution of the system along the dashed arrow in fig . [ fig : deconfine ] . a quantitative analysis of the measured spectra was performed , in order to extract the value of the interchain transfer integral @xmath6 as a function of pressure . the obtained results are compared to the ambient - pressure optical response of different bechgaard - fabre salts and discussed in terms of the expected crossover from insulating ( tmttf)@xmath16 to metallic ( tmtsf)@xmath16 compounds . single crystals of @xmath0  were grown by standard electrochemical growth procedure . a diamond anvil cell equipped with type iia diamonds suitable for infrared measurements was used to generate pressures up to 6  gpa . finely ground csi powder was chosen as quasi - hydrostatic pressure transmitting medium . the measured samples had as - grown specular surfaces perpendicular to the @xmath23 axis and a thickness of about 50 @xmath24 m . the lateral dimensions were approximately @xmath25  @xmath24m@xmath26 . the pressure in the cell was measured using the ruby luminescence method.@xcite polarized reflectivity spectra were recorded at room temperature utilizing a bruker ifs 66v / s ft - ir spectrometer in the mid - infrared ( mid - ir ) frequency range ( 550 - 8000 @xmath18 ) . to focus the beam on the small sample in the pressure cell , an infrared microscope ( bruker irscope ii ) coupled to the spectrometer was used . the intensity spectrum @xmath27 reflected from the lower air - diamond interface of the diamond anvil served as reference . the obtained reflectivity spectra @xmath28 refer to the absolute reflectivity at the sample - diamond interface calculated according to @xmath29     of @xmath0  at room temperature as a function of pressure for * * e**@xmath30 , * * e**@xmath31 , and * * e**@xmath32 at room temperature . note the different pressure values for * * e**@xmath33 and * * e**@xmath32 . the inset shows the normalized pressure - induced change of the reflectivity @xmath34@xmath35 - 1 at 630  @xmath18.[fig : refl],title="fig : " ] +    where @xmath36 is the intensity reflected from the sample - diamond interface and @xmath37 is the absolute reflectivity of the diamond - air interface . since the refractive index of diamond@xcite @xmath38 is almost frequency independent below @xmath39  @xmath18 , @xmath40 was used in our analysis . furthermore , it was assumed to be pressure - independent.@xcite the reflectivity spectra @xmath28 in the frequency range 1700 - 2700  @xmath18  are affected by multiphonon absorptions of diamond , which can not be perfectly compensated . therefore , this part of the measured spectra is not shown . in order to check the reproducibility of the results , we have performed four high - pressure measurements on different @xmath0crystals . of @xmath0  with the fits , according to eqs . ( [ eq:4p - fit ] ) and ( [ eq : ref - fit ] ) , for * * e**@xmath41 for two pressures ( 0.5 and 3.0  gpa ) ; ( b ) corresponding optical conductivity @xmath42 spectra from the fit compared with the optical conductivity obtained by the kk analysis.[fig : fitting],title="fig : " ] + pressure - dependent reflectivity spectra are presented in fig .  [ fig : refl ] for the polarizations * * e**@xmath30 ( stacking direction ) , * * e**@xmath31 ( direction perpendicular to @xmath1 , in the @xmath1-@xmath43 plane ) , and * * e**@xmath32 ( perpendicular to @xmath1 and @xmath3 ) . since the spectra for * * e**@xmath32 were obtained in a separate experimental run , the pressure values are slightly different from those for * * e**@xmath33 . the * * e**@xmath41 reflectivity spectra include a strong mid - ir band at around 4600  @xmath18  and a number of sharp vibrational modes below 2000  @xmath18 . most of them are totally symmetric electron - molecular - vibration ( emv ) coupled a@xmath44 raman modes of the tmttf molecule.@xcite the application of pressure leads to the gradual increase of the low - frequency part ( @xmath452000  @xmath18 ) of the reflectivity . at the same time , the reflectivity level of the mid - ir band remains almost constant in the whole studied pressure range .    for * * e**@xmath31 the reflectivity level is almost frequency - independent at low pressures ( @xmath452  gpa ) , except for the broad mid - ir absorption band located at around 5300  @xmath18 . the sharp peak at around 700  @xmath18  is related to the f@xmath46 symmetry vibration of the octahedral asf@xmath21 anion.@xcite other weaker vibrational features for * * e**@xmath31 can be assigned to the b@xmath46 and b@xmath47 ir - active modes of the tmttf molecules.@xcite the application of high enough pressure ( @xmath482  gpa ) results in a pronounced increase of the reflectivity in the low - frequency range ( @xmath452000  @xmath18 ) . at the same time , the intensity of the mid - ir band decreases , while its frequency position remains almost unchanged . in contrast to the other studied directions , the reflectivity for * e*@xmath49@xmath50 is very low , almost frequency - independent ( except for the sharp asf@xmath21 vibration at 700  @xmath18 ) , and does not change with the application of pressure up to 6 gpa . the value of the low - frequency reflectivity ( @xmath452000  @xmath18 ) is mostly affected by the drude contribution due to free electrons . therefore , the pressure - induced enhancement of the charge transport should lead to a higher reflectivity at low frequencies . this is illustrated in the inset of fig . [ fig : refl ] , where the change of the reflectivity @xmath51@xmath52 - 1 at 630  @xmath18normalized to the reflectivity at the lowest applied pressure @xmath53 is depicted as a function of applied pressure @xmath54 . the enhancement of the reflectivity along the stacks ( * * e**@xmath55 ) is as large as a factor of two at 5.2  gpa . for * * e**@xmath56 the increase becomes noticeable only above 3  gpa , and reaches one at the highes applied pressure . in contrast , the spectra of * * e**@xmath32 show a very small increase of only @xmath570.2 at 6  gpa , which can not be analyzed quantitatively . in order to extract the optical conductivity for * * e**@xmath30 from the measured pressure - dependent reflectivity spectra , we have carried out a kramers - kronig ( kk ) analysis . the application of this method to our data is not trivial due to two reasons : ( i ) the measured reflectivity spectra cover a limited frequency range ; for the kk analysis they need to be extrapolated to lower and higher frequencies and interpolated within the diamond absorption range ( 1700 - 2700  @xmath18 ) . ( ii ) the standard kk relation between the reflectivity and phase needs to be corrected when it is applied to the sample - diamond interface , @xmath58 , and the necessary correction term contains an _ a priori _ unknown parameter.@xcite    as a first step in the analysis of the reflectivity data @xmath28 , we performed an approximate fitting of the reflectivity data . most of the observed phonon resonances are characterized by an asymmetric lineshape and can not be fitted by a simple three - parameter lorentz model . we therefore used an extended four - parameter factorized oscillator model@xcite for the fitting , which is frequently applied for the analysis of ir phonon modes . the dielectric function in this model is described by @xmath59 where @xmath60 and @xmath61 are transverse and longitudinal frequencies of the @xmath62th polar mode , respectively , @xmath63 and @xmath64 their respective damping constants , and @xmath65 the optical dielectric constant far above the resonance frequencies . the reflectivity is then calculated according to @xmath66 with this procedure we were able to fit the measured reflectivity spectra rather precisely . as examples , we show the results of the fits together with the measured spectra for two pressures ( 0.5 and 3  gpa ) in fig . [ fig : fitting](a ) . spectra of @xmath0  for * * e**@xmath30 as a function of pressure obtained by the kk analysis . inset : pressure dependence of the spectral weight @xmath67 according to eq . ( [ eq : sumrule]).[fig : cond - a],title="fig : " ] +    the obtained fits were used for the extrapolation of the experimental spectra to lower and higher frequencies necessary for the kk analysis . the reflectivity in the diamond absorption range ( 1700 to 2700 @xmath18 ) was linearly interpolated . the kk relation for the phase of the reflectivity @xmath58 has the following form:@xcite @xmath68 ,    \label{eq : kk}\ ] ] where @xmath69 is the position of the reflectivity pole on the imaginary axis in the complex frequency plane . in case of measurements on the sample - air interface , @xmath69 tends towards infinity and the second term vanishes . for the sample - diamond interface the second term must , however , be taken into account . the criterium for the proper value of @xmath69 is the agreement between the optical conductivity obtained by the kk analysis and that from the initial fit . for the illustration of the applied data analysis a comparison of these two spectra is given in fig . [ fig : fitting](b ) for two different pressures ; obviously , the agreement between the spectra is very good . for all studied pressures the best correspondence between the optical conductivity spectra obtained by fitting and the output of the kk analysis is achieved for the parameter @xmath70 @xmath18 . the pressure - dependent optical conductivity spectra @xmath42 of @xmath0  obtained by the kk analysis are presented in fig . [ fig : cond - a ] . for the lowest applied pressure ( 0.5 gpa ) the conductivity spectrum is very similar to the published ambient - pressure data.@xcite it can be viewed as a broad absorption band centered at around 2450  @xmath18  with a fine structure of vibronic modes.@xcite this mid - ir band is usually attributed to the optical transitions across the mott - hubbard gap.@xcite according to our results , the application of pressure causes a gradual shift of the position of the band to lower frequencies , accompanied by the onset of a drude - like conductivity . the spectral weight @xmath67 of the absorption band was estimated by integrating the optical conductivity @xmath42 up to @xmath71  @xmath18@xmath72 here , @xmath73 denotes the on - chain plasma frequency , @xmath74 the electron density , and @xmath75 the electronic band mass in @xmath1 direction . the spectral weight obtained according to eq . ( [ eq : sumrule ] ) is almost pressure - independent ( see inset of fig .  [ fig : cond - a ] ) . the slight monotonic increase of the spectral weight ( about 5% from the lowest to the highest pressure ) can be attributed to the increase of the charge density due to the pressure - induced lattice compression.@xcite     spectra of various bechgaard - fabre salts for * * e**@xmath76 : ( a ) @xmath0  at 5.2 gpa ; ( b ) ( tmtsf)@xmath20clo@xmath77 ; ( c ) ( tmtsf)@xmath20pf@xmath21 ; ( d ) ( tmtsf)@xmath20asf@xmath21 ; ( e ) ( tmttf)@xmath20br ; ( f ) ( tmttf)@xmath20pf@xmath21 ; ( g ) @xmath0  at 0.5 gpa . the spectra ( b)-(f ) are taken from refs .  . [ fig : anion - compare],title="fig : " ] +    based on the spectral weight , we estimated the band mass @xmath75 of @xmath0  for the lowest pressure : with a charge density of @xmath78  @xmath79 , as calculated from the lattice parameters,@xcite we obtain a band mass @xmath75@xmath57 @xmath80 ( @xmath81 being the free electron mass ) according to eq . ( [ eq : sumrule ] ) . this value is larger than the typical band mass in metallic tmtsf salts , where it is almost equal to the free electron mass.@xcite the difference is in good agreement with the calculated values of the room temperature bandwidth for various bechgaard - fabre salts.@xcite    in fig . [ fig : anion - compare ] the room - temperature * * e**@xmath41 optical conductivity spectra of various bechgaard - fabre salts@xcite are compared with the corresponding results of the present @xmath0  study at the lowest and highest applied pressure . for clarity , the curves are offset vertically and arranged from the bottom to the top with increasing chemical pressure , according to the generic phase diagram.@xcite the spectrum of @xmath0  at 5.2  gpa is plotted above all other spectra . obviously , the application of _ chemical _ pressure shifts the strong mid - ir band to lower frequencies , suggesting the decrease of the mott - hubbard gap energy which is mostly related to decrease of the stack dimerization.@xcite as demonstrated in figs .  [ fig : cond - a ] and [ fig : anion - compare ] , the same effects are also induced by the application of _ external _ pressure to @xmath0 . however , a more detailed comparison reveals that the above - mentioned effects induced by chemical pressure appear to be less pronounced in the case of external pressure : even at a relatively high pressure of 5.2  gpa the on - chain optical conductivity of @xmath0  differs from the ambient - pressure spectrum of ( tmtsf)@xmath20asf@xmath21 ; it rather resembles that of ( tmttf)@xmath20br . thus , the spectral weight of the _ on - chain _ optical conductivity strongly depends on the type of cation ( tmtsf or tmttf ) , while external pressure mainly changes the _ inter - chain _ coupling and the energy of the mott - hubbard gap . this issue will be further discussed in section  [ b - direction ] . a measure of the strength of the electronic interaction can be obtained by an analysis of the shape of the mir - ir band , i.e. , the optical transitions across the mott - hubbard gap.@xcite the optical conductivity well above the mott - hubbard gap of a 1d insulator is expected to follow a power - law decay according to@xcite @xmath82 where @xmath83 is the order of commensurability and @xmath84 is the ll exponent characterizing the strength of electronic interaction . @xmath85 corresponds to noninteracting electrons , @xmath86 to repulsive , and @xmath87 to attractive interaction between the electrons . based on the analysis of the low - temperature optical conductivity of several tmtsf compounds , the power - law exponent @xmath88 was determined to be @xmath89 , which corresponds to @xmath90.@xcite since the studied energy range is significantly higher than the thermal energy ( @xmath91 ) , eq . ( [ eq : power - law ] ) should be valid also for the room - temperature optical conductivity . an analysis of the corresponding spectra in ref . gives @xmath88=@xmath92 ( @xmath93 ) . this room - temperature value of @xmath5 is only slightly smaller than the low - temperature value , indicating the weak influence of temperature on the electronic correlations in the bechgaard - fabre salts .     of @xmath0  at different pressures . the straight lines are fits according to eq . ( [ eq : power - law ] ) . inset : power - law exponent @xmath88 and corresponding ll exponent @xmath5 as a function of applied pressure.[fig : slopes],title="fig : " ] +    we applied eq . ( [ eq : power - law ] ) to the analysis of our optical conductivity spectra of @xmath0 . at ambient pressure we obtain @xmath94 ( @xmath95 ) , in agreement with the value obtained by fitting the published spectra of ( tmttf)@xmath20pf@xmath21 and ( tmttf)@xmath20br ( 3.0 and 2.5 , resp.).@xcite this demonstrates the reliability of our analysis . the pressure - induced change of the power - law exponent is illustrated in fig . [ fig : slopes ] , where the measured conductivity spectra are plotted in double - logarithmic scale . one can clearly observe a gradual decrease of the slope with increase of pressure . the inset of fig . [ fig : slopes ] shows the power - law exponent @xmath88 and the respective value of @xmath5 as a function of pressure . at 5.2  gpa we obtain @xmath96 and @xmath97 for @xmath0 , i.e. , the typical value for tmtsf compounds . obviously , the application of external pressure leads to a gradual decrease of the effective repulsive interaction strength in the studied compound as it was generally claimed in the earlier work.@xcite since @xmath84 is close to the value @xmath98 at which quarter - filling umklapp scattering is irrelevant , this means that making the system less interacting decreases strongly the mott gap as it was observed in the previous study.@xcite this could explain the large variation of the mott gap with pressure assuming that the quarter - filled umklapp process dominates . thus , in principle , the effect of the half - filled umklapp due to the stack dimerization is not necessary in order to explain the electronic behavior in the bechgaard - fabre salts , although its influence can not be completely excluded .    furthermore , the interchain hopping is enhanced as the distance between the chains is reduced by external pressure . this effect can be best followed in the transverse optical response as discussed in the next section .      in the reflectivity spectrum for * * e**@xmath99 , i.e. , along the direction perpendicular to @xmath1-@xmath43 plane , no pronounced electronic excitations are found in the whole studied energy range . this indicates an insulating behavior similar to the @xmath100 axis optical conductivity in ( tmtsf)@xmath20clo@xmath77.@xcite the application of pressure up to 5.8  gpa does not induce any appreciable change in the reflectivity spectra .    [ cols="^,^,^,^,^,^,^",options="header " , ]        of @xmath0  for * * e**@xmath31 together with the drude fit according to eq . ( [ eq : drude]).[fig : drude - b],title="fig : " ] +    the optical response along the @xmath3 direction , i.e. , perpendicular to the stacks of tmttf molecules , shows a pronounced pressure - induced increase of the low - frequency reflectivity for pressures above 2  gpa . this increase signals the onset of a drude - like contribution and thus the deconfinement of the charge carriers . the drude - like response can be described quantitatively by estimating the transfer integral in @xmath3 direction . for this purpose , we fitted the low - frequency part of the measured @xmath58 spectra for pressures above 2 gpa with a simple drude model , where the dielectric function is described by @xmath101 with @xmath102 being the background dielectric constant , @xmath103 the plasma frequency , and the optical scattering rate @xmath104 . @xmath104 was considered to be pressure - independent , since at room temperature it is mostly determined by the scattering on thermal phonons . its value @xmath105  @xmath18  agrees with literature data on @xmath22.@xcite the measured reflectivity spectra of @xmath0  together with the fits are shown in fig . [ fig : drude - b ] . the parameters of the fits and the related dc conductivity , @xmath106 , where @xmath107 is the vacuum dielectric constant , are listed in table  [ tab : bandwidth ] .    with the obtained values of @xmath103 we can calculate the interchain transfer integral @xmath4 according to @xmath108 valid for an open fermi surface.@xcite here @xmath109 denotes the unit cell volume , @xmath110 the transfer integral along the stacks , and @xmath43 the interstack separation . in a simple tight - binding model @xmath110 should scale inversely proportional to the effective mass . from the analysis of the spectral weight given in section  [ a - direction ] , we found an effective mass @xmath111 which is almost pressure - independent . for bechgaard salts like @xmath22 , which have an effective mass equal to the free electron mass , @xmath110 is approximately 0.25  ev.@xcite therefore , we estimated the transfer integral to be @xmath112  ev for @xmath0  at all pressures . due to the lack of high - pressure structural data on @xmath0 , we had to make several assumptions , in order to obtain @xmath43 and @xmath109 as a function of pressure : the parameter @xmath43 was assumed to be linearly pressure - dependent in the whole studied pressure range , with the pressure coefficient 2.2  %  gpa@xmath113 taken from the data of gallois et al.@xcite of @xmath22  for pressures up to 1.6  gpa . the unit cell volume @xmath109 above 2  gpa was estimated from the pressure - induced change of the background dielectric constant @xmath114 according to the clausius - mossotti relation@xcite @xmath115 where @xmath116 is the electronic polarizability of the @xmath0unit cell in @xmath3 direction . the value of @xmath116 was chosen to fit the known unit cell volume at pressures below 1.6  gpa@xcite and was assumed to be pressure - independent . the carrier density was calculated assuming one carrier per unit cell . the resulting values of the unit cell volume are listed in table  [ tab : bandwidth ] . we verified that the variation of the lattice parameters affects the calculated charge transfer integral only very weakly compared to the uncertainties of the drude fit .    the pressure dependence of the transfer integral @xmath4 calculated according to eq . ( [ eq : kwak ] ) is shown in fig . [ fig : bandwidth - b ] . above 2.0  gpa , @xmath4 increases approximately linearly with increasing pressure , with a linear pressure coefficient of @xmath575  mev / gpa . at a pressure of 4.5 gpa the value of the transfer integral in @xmath0  becomes equal to the ambient - pressure value in @xmath22(@xmath117  mev).@xcite according to the generic temperature - pressure phase diagram of the bechgaard - fabre salts@xcite the separation on the pressure scale between @xmath0  and @xmath22  is about 3  gpa , i.e. , appreciably smaller than 4.5  gpa . probably , the reason for this discrepancy is the way in which the generic phase diagram was obtained . the various bechgaard - fabre salts are arranged in the diagram mostly according to their ground states ( spin - peierls , superconducting etc . ) , whose nature is not determined solely by interchain coupling . therefore , one would not expect a perfect agreement with our room temperature results related to the equivalence of @xmath4 values .    with increasing pressure the deconfinement transition ( trajectory 2 in fig . [ fig : deconfine ] ) is expected to occur when the mott - hubbard gap @xmath118 equals approximately @xmath119.@xcite in order to verify this condition , we estimate the pressure dependence of the charge gap @xmath120 using the expression for a 1d quarter - filled mott - hubbard insulator@xcite @xmath121 where @xmath122 is the bandwidth and @xmath123 the coupling constant of the corresponding umklapp process . the pressure dependence of @xmath5 was obtained from the analysis of the longitudinal optical response . as already discussed , our experiment indicates that the bandwidth @xmath124 does not strongly depend on pressure and amounts to 560  mev . using the value of @xmath125  mev known from dc resistivity measurements,@xcite the umklapp coupling constant at ambient pressure can be calculated according to eq . ( [ eq : gap ] ) and amounts to @xmath126  mev . since the increase of @xmath5 with increasing pressure ( see inset in fig .  [ fig : slopes ] ) indicates a decrease of the electronic interaction , @xmath127 should decrease as well . the functional relation between @xmath127 and @xmath5 is , however , not known for the case of a strong non - local electron repulsion when @xmath128 . therefore , we calculated the pressure - dependent charge gap @xmath120 assuming that all parameters in eq . ( [ eq : gap ] ) except @xmath5 are pressure - independent . the resulting values of @xmath120 as a function of pressure is depicted in fig . [ fig : bandwidth - b ] . accordingly , for @xmath0  the deconfinement transition occurs at @xmath129  gpa . at the same time , the interchain transfer integral @xmath4 related to single - particle hopping should renormalize to zero in the confined state@xcite this is directly confirmed by our experimental results : the drude - like behavior is not observed below 2  gpa indicating a strong suppression of the interchain charge transfer , as illustrated by the full line in fig . [ fig : bandwidth - b ] .    according to the temperature - pressure phase diagram , which is based on pressure - dependent resistivity measurements on the closely related ( tmttf)@xmath20pf@xmath21 compound,@xcite the deconfinement transition occurs at a pressure of 1  gpa at room temperature . this value is somewhat lower than the pressure of the deconfinement transition in @xmath0  obtained in the present study . possibly , this discrepancy is related to the criteria which were used to specify the transition : in the mentioned work of moser _ _ et al.__@xcite the deconfinement transition was identified based on the resistivity behavior along the @xmath1 axis , whereas in the present work we refer to the onset of the charge transport along the @xmath3 direction . ( open circles ) and the charge gap @xmath118 ( full circles ) as a function of pressure for @xmath0 . the dashed lines are the linear fits of @xmath4 and @xmath119 . @xmath130 is the pressure at which the deconfinement transition takes place . the full line schematically shows the vanishing of @xmath4 at the transition pressure . [ fig : bandwidth - b],title="fig : " ] +    the dimensional crossover ( trajectory 3 in fig . [ fig : deconfine ] ) is expected to reveal itself as an onset of a sharp plasma edge@xcite in the @xmath3 direction , i.e. , the plasma frequency @xmath103 should become significantly larger than the scattering rate @xmath104 ( @xmath573800 cm@xmath113 ) . according to the values of @xmath103 given in table  [ tab : bandwidth ] , the drude response is overdamped ( i.e. , @xmath132 ) in the whole studied pressure range . thus , even at 6  gpa @xmath0  is presumably in a 1d ll state , and the dimensional crossover is expected to occur only at even higher pressure . according to the linear extrapolation ( not shown ) of the pressure dependence of @xmath103 ( which is proportional to @xmath4 ) the condition @xmath133 for the dimensional crossover is fulfilled at a pressure of @xmath134 = 12.8  gpa . the estimated value of the interchain transfer integral at this pressure is @xmath135  mev . the so - obtained value of @xmath134 does not agree with the results of dc resistivity measurements along @xmath50 published by moser et al.@xcite there , a strong increase of the dimensional crossover temperature @xmath136 with application of pressure was demonstrated , such that @xmath136 should reach room temperature at 3  gpa for ( tmttf)@xmath20pf@xmath21 and at 1  gpa for @xmath22 . the former value is smaller by a factor of about four compared to our estimation based on the criterium of the underdamped drude response . therefore , we suggest that this criterium should not be used for the precise determination of the dimensional crossover pressure or temperature . rather , the observation of a well underdamped drude response ( @xmath137 ) could be taken as a criterium that the system can no longer be considered as one - dimensional . in summary , we studied the mid - ir reflectivity of @xmath0  at room temperature with the incident radiation polarized along the @xmath1 , @xmath3 , and @xmath50 direction as a function of quasi - hydrostatic pressure . the pressure dependence of the optical conductivity along the molecular stacks ( * e*@xmath49@xmath1 ) demonstrates a gradual onset of a drude - like response accompanied by a decrease of the charge gap with increasing pressure . we found the band mass along the stacks @xmath138 to be weakly pressure - dependent . assuming that the 1d hubbard model describes the properties of @xmath0  correctly in the high - energy limit ( i.e. , above 2500  @xmath18 ) , where the optical conductivity decays according to a power law @xmath139 , we found the exponent @xmath88 to decrease from 2.9 at ambient pressure to 2.0 at 5.2  gpa . the corresponding ll exponent @xmath5 increases from 0.13 to 0.19 indicating a weakening of the electronic interaction . perpendicular to the stacks the onset of an overdamped drude - like optical conductivity is observed above 2.5  gpa for * e*@xmath49@xmath3 . the analysis of the data suggests a linear pressure - induced change of the interchain transfer integral @xmath4 . the absence of a noticeable drude component below @xmath140  gpa and its appearance above this pressure is interpreted in terms of a pressure - induced deconfinement transition at @xmath130 . for this pressure value the interchain transfer integral @xmath4 is approximately equal to half the charge gap energy @xmath120 . for * e*@xmath49@xmath50 @xmath0  remains insulating up to the highest applied pressure . we thank g. untereiner for crystal growth , as well as m. dumm and t. giamarchi for discussions and comments . financial support by the dfg ( emmy noether program ) is acknowledged . f. nad , p. monceau , c. carcel , and j.  m. fabre , phys . b * 62 * , 1753 ( 2000 ) ; d.  s. chow , f. zamborszky , b. alavi , d.  j. tantillo , a. baur , c.  a. merlic , and s.  e. brown , phys . rev .  lett . * 85 * , 1698 ( 2000 ) ; p. monceau , f. nad , and s. brazovskii , _ ibid _ * 86 * , 4080 ( 2001 ) .
Abstract:

Answer: we studied the pressure dependence of the room - temperature infrared reflectivity of @xmath0  along all three optical axes . this anisotropic organic compound consists of molecular stacks with orbital overlap along the @xmath1 direction ; due to electronic correlations the system is a quasi - one - dimensional mott insulator with a charge gap @xmath2  mev . the gap is gradually reduced with increasing external pressure , accompanied by the onset of a drude contribution along the stacking direction . in the perpendicular @xmath3 direction a drude - like optical response is observed for pressures above 2  gpa . this behavior is interpreted in terms of a deconfinement of the electrons in a one - dimensional mott insulator , i.e.  an insulator - to - metal transition which occurs when the interchain transfer integral @xmath4 is approximately equal to half of the charge gap energy . we estimate the values of @xmath4 and the luttinger liquid parameter @xmath5 as a function of pressure .


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: interstellar shocks are seen in association with various astrophysical objects , from wind - blowing massive stars @xcite , mass - ejecting giant stars ( e.g. , @xcite ) , expanding supernova remnants ( snrs ; e.g. , @xcite ) , to energetic pulsars @xcite . in these objects , fast moving ejecta collide with the surrounding interstellar medium ( ism ) , and the strong interaction drives a shock front . depending on the kinetic energy carried into the shock and local ism conditions , the postshock hot material can be bright at different wavelengths . observations of shocks not only help study shock physics , but also provide probes of the associated objects and ism . depending on temperature @xmath3 , the sound speed @xmath4 in the ism [ @xmath5 km s@xmath6 is in the range of 1100 km s@xmath7 . as a result , when stellar objects have high space velocities @xmath8 , they can be moving supersonically in the ism . in such cases , a shock will appear as a cometary structure ( so - called bow shock ) , as interstellar material is swept up into a dense cone - type shell ( e.g. , @xcite ) . one type of bow shock is known to be associated with pulsars . having space velocities of @xmath9100 km s@xmath7 ( e.g. , @xcite ) and strong winds ( which carry much of the spin - down energy @xmath10 of pulsars ) , pulsars moving in the ism are candidates to drive bow shocks . thus far , six pulsars have been found to have a typical bow shock , revealed by the shocks optical balmer line emission @xcite . the detected line nebulae are understood to arise from de - excitations of neutral h atoms in the shocked ambient gas , following the processes of collisional excitation or charge - exchange ( e.g. , @xcite ) . in addition , a pulsar bow shock may also be revealed by the detection of a termination shock of the pulsar wind . the termination shock lies inside a bow shock , and can appear bright at x - ray or radio energies due to synchrotron radiation ( @xcite ; @xcite ) . in our _ spitzer _  infrared array camera ( irac ) survey of 7 relatively young pulsars , for the purpose of probing the general existence of debris disks around isolated neutron stars ( wang et al . 2013 , in preparation ) , a mid - infrared ( mir ) cometary nebula was serendipitously detected near one of the pulsars , psr j1549@xmath04848  ( hereafter j1549 ) . naively thinking that the mir structure is associated with the pulsar , we would have found the first mir bow shock that arises from the interaction between the pulsar wind and ism . we thus made follow - up 24 and 70  @xmath1 m imaging and 7.514.4  @xmath1 m spectroscopic observations with _ spitzer _ and optical observations with the magellan telescopes after the initial detection of the nebula at 8.0 @xmath1 m . however , our detailed data analyses raise questions about the association between the cometary nebula and the pulsar . on the basis of our observational results , the origin of the nebula is not certain . here we report our multi - wavelength observations , data analyses , and results . the observations and data reduction are described in  [ sec : obs ] , and the data analyses and results are presented in   [ sec : ana ] and   [ sec : res ] , respectively . we discuss the possible origin of the cometary nebula in  [ sec : disc ] and summarize the results in   [ sec : sum ] . our _ spitzer _   and ground - based observations of the cometary nebula are summarized in table  [ tab : obs ] . below we describe them and the related data reductions .        as part of our survey program , the initial _ spitzer _  observations of j1549  were carried out on 2007 september 12 . the imaging instrument used was the infrared array camera ( irac ; @xcite ) . the field of our target was imaged in the pair of 4.5 and 8.0 @xmath1 m channels of the irac simultaneously . the detectors at the short and long wavelengths were insb and si : as devices , respectively , with 256@xmath11256 pixels and a plate scale of 12/pixel . the field of view ( fov ) was 52@xmath1152 . the frame time was 30 s , with 26.8 effective exposure time per frame . the total exposure time at each channel was 17.87 min . a cometary nebula was detected in the 8.0 @xmath1 m image ( figure  [ fig : field ] ) around j1549 , but not in the 4.5 @xmath1 m image .    to study the nebula and derive its properties , we requested _ spitzer _  observations at the longer wavelengths . using multiband imaging photometer for _ spitzer_(mips ; @xcite ) , 24 and 70 @xmath1 m broad - band imaging of the source field was carried out on 2008 august 30 . the 24 @xmath1 m detector was a 128@xmath11128 pixel si : as array , with a pixel size of 255 and a 54@xmath1154 fov . the 70 @xmath1 m detector was a 32@xmath1132 pixel ge : ga array . we chose its wide fov observing mode for imaging , which provided a pixel scale of 985@xmath111006 and a 52@xmath1126 fov . at 24 and 70 @xmath1 m , the exposure times were 11.8 and 10.5 min , respectively ( frame times @xmath12 3 sec ) . we started with the basic calibrated data ( bcd ) , provided by _ spitzer _ science center ( ssc ) . the data were produced from raw images through the irac and mips data pipelines ( irac software version was s16.1.10 and mips software version was s18.1.0 ) in ssc . the detailed reduction in the pipelines can be found in the irac data handbook ( version 3.0 ) and mips data handbook ( version 3.3.1 ) .    for irac data , we first used the irac artifact mitigation program , distributed as contributed software by ssc , to clean the bcd images . several bright stars in the target field caused column pull - down , row pull - up , or optical banding artifacts . we then corrected each bcd image for array location dependence , using the array location correction images provided by ssc . to combine bcd images into a final post - bcd ( pbcd ) mosaic at each channel , we used mopex , an ssc s package for reducing and analyzing imaging data .    for mips 24 @xmath1 m data , the bcd images do not have any artifacts that need corrections . we used mopex to combine the bcd images into a final pbcd mosaic . the mips 70 @xmath1 m bcd images contain horizontal and vertical strips , which were due to response variations and stim latents ( stims were mips s built - in calibration light sources , which flashed during an observation sequence ; for details , see mips data handbook ) , respectively . we used the gert package , which is provided for mips 70 & 160 @xmath1 m data reductions , to clean up the bcd images . during the reduction , the region that contains our source was masked out . the bcd images were combined into one final pbcd mosaic using mopex . our _ spitzer _  spectroscopic observations were made on 2009 september 14 . the instrument was the infrared spectrograph ( irs ; @xcite ) . we used the short - low ( sl ) module of the irs for the observations , which provided a wavelength coverage of 7.514.4 @xmath1 m ( first order ) and a resolving power of 61120 . the detector was a 128@xmath11128 arsenic doped silicon ( si : as ) array , with a pixel scale of 1.8/pixel . the slit in the first order had a length@xmath11width of 57@xmath1137 . we obtained a set of 16 spectra at each of two default nod positions along the slit . the integration time for each individual spectrum was 60.95 s , and the total on - source time was 32.5 min . in addition , since our target is an extended source , we used the exact same observing mode and obtained spectra of a nearby sky region for background subtraction . we started with the bcd images , produced from the ssc irsx channel-0 software ( version s18.7.0 ) . we first combined all the sky bcd images into one master sky image and subtracted this master image from each of the target bcd images . we then used the irsclean_mask program , provided by ssc , to construct a super mask by combining the rogue pixel mask for the campaign irsx009500 ( in which our irs observations were made ) and additional bad pixels identified by irsclean_mask . the super mask was used in the irsclean program for cleaning up bad pixels in the bcd images . the bcd images were combined into one final image by using the coa2d.pro program . we used spice , an irs spectrum extraction software provided by ssc , to extract the spectrum of the nebula . the width of the extraction box was 15 pixels ( 27 ) , which covers the 2-d spectrum image well . the spectra at two nod positions were extracted and average - combined into the final source spectrum .      launched on 2009 december 14 , the wide - field infrared survey explorer ( wise ) mapped the entire sky at 3.4 , 4.6 , 12 , and 22 @xmath1 m ( called w1 , w2 , w3 , and w4 bands , respectively ) in 2010 with fwhms of 6.1 , 6.4 , 6.5 , and 12.0   in the four bands , respectively ( see @xcite for details ) . the wise all - sky images and source catalogue were released in 2012 march . we downloaded the wise image data of the source field from the infrared processing and analysis center ( ipac ) , and found that the cometary nebula was detected in the w3 and w4 bands but not in the first two bands . the wise observations were made between 2010 feb . 22 - 25 and the depth of coverage was @xmath1313.7 pixels ( corresponding to 110  s on - source integration time ) . we also observed the target field in @xmath14 broad band and h@xmath15 narrow band on 2008 june 15 , using inamori - magellan areal camera and spectrograph ( imacs ) on the 6.5-m baade magellan telescope at las campanas observatory in chile . the long camera ( known as f/4 ) of imacs was used . the detector of the camera consists of 8 2048@xmath114096 pixel@xmath16 site st-002a ccds . under the f/4 imaging mode , the camera provides a field coverage of 154@xmath11154 and a pixel scale of 0111 pixel@xmath7 . the exposure times of the images were 4 min in @xmath14 and 18.3 min in h@xmath15 . the observing conditions were good , with seeing being 065 . in addition , we observed standard stars pg1633 + 099 @xcite in @xmath14 band and ltt7379 @xcite in h@xmath15  for flux calibration . we included near - infrared ( nir ) @xmath17 imaging data of the target field , obtained on 2006 may 17 as part of our ir survey of a few young pulsars ( like j1549 ) . the observation was also made with the baade magellan telescope . the nir camera was persson s auxiliary nasmyth infrared camera ( panic ; @xcite).the detector was a rockwell hawaii 1024@xmath111024 hgcdte array , having a field of view ( fov ) of 2@xmath112   and a pixel scale of 0125 pixel@xmath7 . the total on - source exposure time was 22.5 min . during the exposure , the telescope was dithered in a 3@xmath113 grid with offsets of 10  to obtain a measurement of the sky background . the observing conditions were good , with 04 seeing in @xmath17 . we used the iraf data analysis package for data reduction . the images were bias - subtracted and flat - fielded . in addition for the nir data , a sky image was made by filtering out stars from each set of dithered images in one observation . the sky image was subtracted from the set of images , and then the sky - subtracted images were shifted and combined into one final image of the target field . the pointing of the irac frames is typically accurate to 05 . in order to determine more accurately the position of the cometary nebula and locate j1549  on the images , we astrometrically calibrated the irac 8.0 @xmath1 m image by matching stars detected in this image to two micron all - sky survey ( 2mass ; @xcite ) stars . for this calibration , 190 2mass stars were used . the nominal uncertainty of the calibrated image is dominated by the 2mass systematic uncertainty ( @xmath12015 , with respect to the international celestial reference system ) . the other images were also positionally calibrated by matching them to the 8.0 @xmath1 m image . we studied the cometary structure by semi - quantitatively determining its shape . the 8.0 @xmath1 m image was used for the detailed study because it has the best resolution among the _ spitzer _  and wise detections . we first removed all point sources near the cometary structure . using the software package apex ( multi - frame extraction ) provided by ssc , point sources in the field were identified and subtracted from each bcd image . the residual frames were then combined into a pbcd mosaic , which is shown in figure  [ fig : bow ] . we then manually selected small regions along the three parts of the nebula : the two wings and the head region . narrow box regions set parallel to the three nebula parts were defined , with each box having a length@xmath11width of 5@xmath111 pixel@xmath16 ( 61@xmath1112 ) . photons within each box were counted and the uncertainty on the photon counts was estimated as the square - root of the total counts . this way a brightness curve that was perpendicular to the narrow box regions ( or to the nebula parts ) was obtained . using a gaussian function , we performed least - squares fitting to the brightness curve and thus determined the location of a brightness peak point . the uncertainty for each location was also estimated from the fitting . in total 19 such brightness peak points were determined ( see figure  [ fig : bow ] ) . we obtained the flux of the head region of the cometary nebula at each band . a circular region with a radius of 21 , which well covers the head region , was used as the photometry aperture . a box region of 80@xmath1150  near the head was used for estimating the sky background brightness . for photometry at 8.0 @xmath1 m , aperture correction was required because irac imaging uses point sources for flux calibration . we used the function form provided by ssc to calculate the correction factor and found a value of 0.83 for a radius of 21 . the flux value at 8.0 @xmath1 m is the measured flux value times the correction factor . we also derived the flux upper limit for the 4.5 @xmath1 m non - detection , because the deep image provides an additional constraint on mir emission from the cometary nebula . following the guidelines given for aperture photometry on the wise atlas images , wise w3 and w4 fluxes were also obtained . the broad - band flux measurements and 4.5 @xmath1 m upper limit are given in table  1 . the shape of the ir cometary nebula indicates that it is likely a typical bow shock . the location of the associated object that drives the ir nebula can thus be estimated by fitting the nebula with a standard bow - shock shape , since the distance from an associated object to the apex of the associated bow shock ( so - called  standoff distance " ) sets the length scale of the bow shock . we considered an analytic solution derived by @xcite , @xmath18 , where @xmath19 is the distance between the bow shock shell and the associated object , @xmath20 is the polar angle from the axis of symmetry ( i.e. , the associated object s moving direction ) , and @xmath21 is the standoff distance . with this solution , we fit the positions of the brightness peak points obtained in  [ subsec : mor ] . the free parameters in the fitting were @xmath22 , @xmath23 location of the associated object , @xmath21 , and the rotation angle of the axis of the bow shock in the plane of the sky . to simplify the fitting , we assumed that the axis of symmetry of the bow shock lies nearly in the plane of the sky ( note that an inclination angle such as 60  will result in a smaller projected distance between the associated object and the ` apparent ' apex of the bow shock ) . the obtained ( 3@xmath24 ) region in which the associated object is located is indicated as a dotted contour in figure  [ fig : bow ] . it is a small region very close to the head of the cometary nebula , and approximately 15  away from j1549  ( whose position is r.a . = 15@xmath2549@xmath262115 , decl . = @xmath04848374 , equinox j2000.0 , and has an uncertainty of @xmath121 ; @xcite ) . because the data points of the brightness peaks do not provide tight constraints along the axis of symmetry , the region is more extended along this direction . the best - fit location and bow shock solution are shown as a plus sign and a dashed curve , respectively , in figure  [ fig : bow ] . as can be seen , the standard bow shock shape generally fits the data points , although the tail of the nebula is wider , deviating away from the model . from the best fit ( @xmath27=14 for 15 degrees of freedom ) , we found @xmath28 .      in figure [ fig : spec ] , we show the broad - band and _ spitzer _  irs spectrum of the head region of the cometary nebula , which is dereddened with @xmath29 @xcite using the interstellar extinction laws given by @xcite for wavelengths @xmath308 @xmath1 m and by @xcite for wavelengths @xmath318 @xmath1 m spitzer_/irac 4.5 @xmath1 m flux upper limit is also included . the irs spectrum was flux calibrated by equaling the average value of its 7.58.5 @xmath1 m region to the irac 8.0 @xmath1 m broad - band flux value . a few emission lines were detected in the irs spectrum , which are polycyclic aromatic hydrocarbon ( pah ) 11.3 , 8.6 , and likely 7.6 @xmath1 m . these pah lines are commonly seen in many emission nebulae and arise from vibrational modes of pah in dust ( e.g. , @xcite ) . two h@xmath2 rotational lines were also detected , although weakly . they are h@xmath2s(3 ) 9.7 @xmath1 m and h@xmath2s(2 ) 12.3  @xmath1 m . in general , h@xmath2 lines are associated with photodissociation regions ( pdrs ; e.g. , @xcite ) or shock interaction regions ( e.g. , @xcite ) . we note that the broad dip at @xmath139.7 @xmath1 m is not the well - known silicate absorption feature . rather , its appearance is caused by the nearby strong pah emission features @xcite . the total galactic reddening along the line of sight to the source is only approximately @xmath29 ( for @xmath32=3305 , @xmath33=43 ; @xcite ) . if the dip is caused by interstellar reddening , @xmath34 should be @xmath35 , which is too large to be consistent with those estimated from the all - sky dust @xcite or neutral hydrogen density @xcite maps . also , as shown below in this section , the local column density of the cometary nebula is much lower than that implied by the reddening @xmath36 . we first tested using a simple dust grain model @xcite to reproduce the spectrum we obtained . this model was developed to explain the observed ir emission from supernova remnants ( snrs ) , in which dust grains are heated both by collisions with ambient thermal gas and by radiation flux , and could work for the ir nebula in our case . considering the general case that the grains absorption / emission efficiency decreases as @xmath37 , we found that two components with temperatures of @xmath38 k and @xmath39 k ( assuming a simplified , uniform grain size of radius 0.1 @xmath1 m ) , can generally describe the spectrum ( figure  [ fig : spec ] ) . the dust masses thus simply estimated are 1.0@xmath40  g and 1.9@xmath41  g , respectively , where @xmath42 is the source distance in units of kpc .    to better reproduce the irs spectrum , a detailed model including pah particles is needed . we considered the model spectra derived and provided by @xcite , although we note that in their model , dust grains are heated by a field of starlight . according to their calculations , for the case of dust grains being heated by a single radiation intensity , the flux ratio between 24 and 70 @xmath1 m is mainly sensitive to radiation with energy density per unit frequency @xmath43 , where @xmath44 is the interstellar radiation field ( isrf ) estimated by @xcite for the solar neighborhood and @xmath45 is a dimensionless scaling factor . we thus found that @xmath46 is required for our data points . however , the @xmath46 model spectra have much stronger pah features than those in our irs spectrum ( even with the lowest pah mass fraction @xmath47 ; @xmath48 is found to be applicable to the dust in the milky way ) . in order to have a good fit to the irs spectrum , we found that the dust should be heated by a single radiation intensity , @xmath49 , plus by a power - law distribution of starlight intensities ranging from @xmath4910@xmath50 ( where @xmath48 , the power - law index @xmath51 is fixed , @xmath52 , and the fraction of the dust mass that is exposed to the distribution of starlight intensities is @xmath53 ; cf . equation ( 23 ) in @xcite ) . this model spectrum is displayed in figure  [ fig : spec ] , which generally describes the pah features in our spectrum . as a result , the dust mass @xmath54 can be estimated from the total value of the 24 , 70 , and 160 @xmath1 m fluxes ( cf . equation ( 34 ) in @xcite ) . since no flux measurement at 160 @xmath1 m was obtained , @xmath54 should be @xmath55  g and @xmath56  g when assuming the flux at 160 @xmath1 m is equal to that at 70@xmath1 m . the value range is consistent with that derived above from the simple dust grain model . we further estimated the properties of the head region of the bow shock . considering again the dust size of 0.1 @xmath1 m ( and density of 3 g  @xmath57 ; @xcite ) , the total luminosities in the region were 5.7@xmath58 and 1.1@xmath59 erg  s@xmath7 for the 250  k and 55  k dust grains , respectively . simply assuming that the head region is a 21  radius ( or 0.1@xmath42 pc ) sphere , the number density ( dominated by the cold dust grains ) can be estimated to be @xmath1390 @xmath57 , where a mass ratio of 1 to 100 between the dust to gas is used @xcite .      in order to search for possible line emission from the bow shock , we made both @xmath14 and @xmath60 imaging observations of the target field with the magellan telescopes ; the wide - band image was used to serve as the continuum flux . we did not detect any @xmath60 nebula in our observations . both optical @xmath14- and @xmath60-band images of the target field are shown in figure  [ fig : bowrh ] . the @xmath14-band detection limit was @xmath1323 mag , and the @xmath60 flux upper limit ( 3@xmath24 ) was 2.9@xmath61 erg  s@xmath7  @xmath62arcsec@xmath63 . provided with the results we have obtained from the above analyses , we discuss the possible origin of this bow shock by first considering a stellar object and then j1549  as the associated star .          as shown by our analysis in  [ subsec : loo ] , considering an isotropic wind / outflow from an object @xcite , the location of the associated object would likely be within a small area near the head region of the bow shock , which is more than 15  away from psr  j1549 . such ir bow shocks can be driven by either massive stars with a strong wind @xcite or giant stars ejecting mass @xcite , as long as the associated stars are moving supersonically in the ism . the stand - off distance of such a bow shock is given by @xmath64 ( e.g. , @xcite ) , where @xmath65 is the mass loss rate , @xmath66 is the wind velocity , and @xmath67 is the density of the ism . for a massive o / b star with typical values of @xmath68  yr@xmath7 , @xmath69  km  s@xmath7 , and @xmath70  km  s@xmath7 ( see , e.g. , @xcite ) , the required ambient ism hydrogen density to produce the observed bow shock is @xmath71 @xmath57 , where @xmath72 is the mass loss rate in units of @xmath73  yr@xmath7 , @xmath74 is the wind velocity in units of 1000  km  s@xmath7 , and @xmath75 is the star s space velocity in units of 30  km  s@xmath7 . here we have used @xmath76  cm and @xmath77 with @xmath78 ( @xmath79 is the mean nucleus number per hydrogen atom and @xmath80 is the mass of the hydrogen nucleus ) . using the low end value of @xmath81 and @xmath82 ( see below for a discussion of the distance ) , along with the other mass - loss parameters given above , we find @xmath83 , broadly consistent with values expected for the ism density . alternatively , for a giant star with mass ejection , @xmath84  km  s@xmath7 and other parameters similar to those above ( e.g. , @xcite ) , the required ambient density is @xmath85 ; reasonable adjustments to the parameters can easily lead to values more consistent with the ism density . these calculations thus show that it is possible to have either a massive star or a giant star with mass ejection as the object that drives the bow shock . however we evaluated the possible association of the bow shock with any in - field stars and obtained negative results . first there were no stars detected in the associated object region in our optical images down to @xmath86 ( see figure  [ fig : bowrh ] ) . we also analyzed the nearby stars by constructing a color - magnitude diagram of @xmath17 versus @xmath87 ( figure  [ fig : rk ] ) . as can be seen , the stars are consistent with being low - mass main sequence stars ( or white dwarfs ) at distances of @xmath882.5 kpc . massive or giant stars would have to be at distances of @xmath136 kpc and @xmath1310 kpc , respectively , and approximately 0.50.8 kpc away from the galactic plane , a place not likely to form bow shocks due to the low density of the ism . the detection of h@xmath2 lines in the irs spectrum suggests that the associated object should probably be @xmath89 kpc away from the galatic plane , since molecular gas generally inhabits within the mid - plane of the galactic disk @xcite . given their high space velocities ( @xmath9100 km s@xmath7 ) and strong winds , pulsars can naturally drive bow shocks in the ism . detections of bow shocks around pulsars at multiple wavelength ranges have provided sufficient evidence for it @xcite . the existence of ir bow shocks , arising due to irradiation or collisional heating , can be expected . the pulsar j1549  has a spin - down energy of @xmath90 erg s@xmath7 and a distance of 1.5 kpc @xcite , implying that if it were associated with the bow shock , approximately 1.6% of the energy is used to illuminate the dust , which is a reasonable fraction to consider ( note the total luminosity from the head of the bow shock is @xmath91 erg s@xmath7 at the distance of 1.5 kpc ) . on the other hand , the pulsar s kinematic energy can also contribute to dust heating . considering a radius @xmath92  pc ( the head region we have defined when the distance is 1.5 kpc ) interaction area , the energy received by the ambient medium is @xmath93 erg  s@xmath7 , where @xmath94 @xmath57 and @xmath95 km  s@xmath7 , the pulsar s velocity , are used . therefore there is sufficient energy to illuminate the bow shock . the obvious challenge to this interpretation is j1549 s location , which is at least 15   ( or 0.11 pc if at the pulsar s distance ) away from the determined area for the associated object of the bow shock . however , we note that a pulsar wind can be highly anisotropic and for peculiar cases , long jet - like outflows have been detected from pulsars @xcite . calculations for such anisotropic wind cases to produce bow shocks have been considered ( e.g. , @xcite ) . one feature that has been revealed is the off - axis location of an associated star if the star s wind points away from the direction of the stellar motion , for which the star s location would otherwise be determined to lie in the apparent axis of symmetry when simply fitting the stand - off region @xcite . @xcite did further detailed simulations by considering the anisotropy of pulsar winds for pulsar bow shocks , and have shown that a variety of shape changes with respect to the ` standard ' shape of bow shocks . in particular , the broadening of the tails of the bow shocks in their simulations sa and sb ( see figure 13 in their paper ) , when a pulsar runs through a wall of the high - density ism with finite width , is possibly observed in our case ( see figure  [ fig : bow ] ) . we thus did test fitting to our cometary nebula by using equation a11 given in @xcite , an analytic solution for jet - like outflows ( @xmath96 , @xmath97 , and @xmath98 were simply set to be 0 , 0 , and 5 , where @xmath96 sets the angle between the directions of the pulsar s wind and its motion ; for details see @xcite ) . the best - fit solution ( @xmath99 for 15 degrees of freedom ) is shown in figure  [ fig : bow ] , which is nearly the same as that found in section  [ subsec : loo ] . the allowed 3@xmath24 region for the associated object is however substantially elongated toward j1549 . while the pulsar is still 5.8   away from the region , the fitting results suggest that it may be possible for the pulsar to drive the bow shock if it had a highly - collimated wind and the wind were tilted away from the direction of the pulsar s motion . in order to fully explore this possibility , detailed hydrodynamic simulations , such as those carried out by @xcite , are needed . in addition , if j1549  did have a highly collimated pulsar wind , the interaction of its wind with the ambient medium could be detectable at x - ray energies . x - ray imaging should be carried out to search for the interaction , which would lie from j1549  to the stand - off region of the ir nebula , thus verifying the association of the pulsar with the nebula . the bow shock has a direction towards the galactic plane ( figure  [ fig : field ] ) , implying that the associated object is moving in the same direction . high - spatial resolution radio imaging of the pulsar , to detect its proper motion , should also be carried out . if j1549  is also moving in this direction , its chance of being associated with the bow shock would be substantially increased . we have discovered an ir bow shock around the pulsar j1549 , and from our follow - up imaging and spectroscopic observations and wise all - sky survey data , a composite spectrum of its head region from 8 to 70 @xmath1 m , which includes a _ spitzer _ irs 7.514.4 @xmath1 m spectrum , was obtained . in the irs spectrum , a few pah and h@xmath2 emission features were detected . considering the overall spectrum with these features , the emission likely arises from the dust grains in the interaction region heated by collisions and/or radiation flux from the associated object . the dust mass in the region was estimated to be 2@xmath1110@xmath100  g. we analyzed the shape of the bow shock by fitting with a standard analytic solution for such objects , and the derived region in which the associated object would be located is 15  away from j1549 . this raises the question of which object drives this bow shock . no stars were detected in the associated object region down to @xmath101 , and moreover our studies of the in - field stars indicate that none of them are likely bright o / b stars or red giants located near the galactic plane , which are able to drive such bow shocks due to their strong winds or mass ejection . given the considered properties of pulsar winds , including a few peculiar pulsar jet - like outflows , we suggest that it is possible for j1549  to be associated with the ir bow shock if it had a highly collimated wind . hydrodynamic simulations of the interaction of such a wind with the ism can help explore this possibility . in addition , x - ray imaging to detect the interaction of its wind with the ambient medium and high - spatial resolution radio imaging to detect the proper motion of the pulsar will help verify their association . this work is based in part on observations made with the _ spitzer _ space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa . the publication makes use of data products from the wide - field infrared survey explorer , which is a joint project of the university of california , los angeles , and the jet propulsion laboratory / california institute of technology , funded by nasa . this research was supported by national basic research program of china ( 973 project 2009cb824800 ) , and national natural science foundation of china ( 11073042 ) . zw is a research fellow of the one - hundred - talents project of chinese academy of sciences . ps acknowledges partial support from nasa contract nas8 - 03060 . vmk holds a canada research chair and the lorne trottier chair in astrophysics & cosmology , and is a fellow of the royal society of canada . llccc 2006 may 17 & magellan / panic & @xmath17 & 22.5 & + 2007 sept 12 & spitzer / irac & 4.5 & 26.8 & @xmath1020.8 + & & 8.0 & 26.8 & 10.8@xmath1031.6 + 2008 june 15 & magellan / imacs & @xmath14 & 4.0 & + & & h@xmath15 & 18.3 & + 2008 aug 30 & spitzer / mips & 24 & 11.8 & 20.9@xmath1030.9 + & & 70 & 10.5 & 82@xmath10332 + 2009 sept 14 & spitzer / irs & 7.514.4 & 32.5 & + 2010 feb 2225 & wise & 3.4 & 1.8 & + & & 4.6 & 1.8 & + & & 12 & 1.8 & 11.4@xmath1030.6 + & & 22 & 1.8 & 21.1@xmath1031.6 +
Abstract:

Answer: we report on the discovery of an infrared cometary nebula around psr j1549@xmath04848 in our spitzer survey of a few middle - aged radio pulsars . following the discovery , multi - wavelength imaging and spectroscopic observations of the nebula were carried out . we detected the nebula in spitzer irac 8.0 , mips 24 and 70 @xmath1 m imaging and in spitzer irs 7.514.4 @xmath1 m spectroscopic observations , and also in the wise all - sky survey at 12 and 22 @xmath1 m . these data were analyzed in detail , and we find that the nebula can be described with a standard bow - shock shape , and that its spectrum contains polycyclic aromatic hydrocarbon and h@xmath2 emission features . however , it is not certain which object drives the nebula . we analyze the field stars and conclude that none of them can be the associated object because stars with a strong wind or mass ejection that usually produce bow shocks are much brighter than the field stars . the pulsar is approximately 15  away from the region in which the associated object is expected to be located . in order to resolve the discrepancy , we suggest that a highly collimated wind could be emitted from the pulsar and produce the bow shock . x - ray imaging to detect the interaction of the wind with the ambient medium and high - spatial resolution radio imaging to determine the proper motion of the pulsar should be carried out , which will help verify the association of the pulsar with the bow shock nebula .