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: star forming regions and interstellar media have been well studied in all the wavelength regimes from radio to @xmath9-ray . an advantage using far - infrared ( far - ir ) light is that we can trace photo - ionization and dissociation regions without needing to correct for extinction . in the far - ir region , there are important forbidden lines of atoms and ions , such as [ oi ]  63 ,  146 , [ oiii ]  88 , [ nii ]  122 ,  205 , and [ cii ]  158 , which are good tracers of the interstellar media in various gas phases . 30 doradus ( 30 dor ) in the large magellanic cloud ( lmc ) is one of the most famous massive star forming regions ( e.g. @xcite ) . its distance is about 50 kpc and the physical scale of @xmath10 corresponds to 15 pc . in the center , as revealed by the hubble space telescope , there are a lot of massive stars in a compact region , which is named the r136 super star cluster ( @xcite ; @xcite ) . the bubble of @xmath10 in diameter exists around the center of r136 , where gas is almost completely evacuated ( @xcite ; @xcite ) . prominent x - ray and h@xmath11 filamentary and/or shell structures are extended in the north and the west direction from the center @xcite . in the optical emission lines such as [ sii ] @xmath126724 and [ oiii ] @xmath125007 , widely extended components of ionized gas were revealed around the 30 dor region by the university of michigan / ctio magellanic cloud emission line survey ( mcels ; @xcite ) .    from the early era of ir astronomy , the 30 dor region has been studied by various ir observations . using the fifi instrument of the kuiper airborne observatory ( kao ) , @xcite observed the far - ir forbidden lines of [ cii ]  158 and [ oi ]  63 ,  146 in the 30 dor region , and made a line intensity map of [ cii ]  158 with an arc - minute spatial resolution over @xmath13 area . balloon - borne ir telescopes , such as balloon - borne infrared carbon explore ( bice ) , observed the lmc region . bice mapped the entire area of the lmc by the [ cii ] line with a beam of about @xmath14 size . @xcite discussed the photoelectric heating and [ cii ] cooling in the lmc . the first powerful tool for ir spectroscopy was provided by the infrared space observatory ( iso ) operated in the late 1990s . the isocam made line intensity maps of [ neii ]  12.81 and [ neiii ]  15.56 emission lines for the central @xmath15 area of the 30 dor region @xcite . the ionization potentials of ne and ne@xmath16 are 21.56 and 40.96 ev , respectively . therefore , the [ neii ] and [ neiii ] lines are good tracers of highly - ionized gas , similarly to the [ oiii ] line with the ionization potential of 35.12 ev . in the far - ir , the iso / lws observed more than 10 points with different physical conditions for the 30 dor region ( @xcite ; @xcite ) . however , the spatial information provided by the iso / lws was rather poor , because these observations were performed by a single beam .    for technical reasons , far - ir spectral mapping observations with high spatial resolution were difficult . after iso , akari takes a unique position in the far - ir spectroscopy . the japanese ir astronomical satellite , akari @xcite , which was launched on february 21 , 2006 ( ut ) , has two kinds of focal - plane instruments ; infrared camera ( irc ) @xcite and far - infrared surveyor ( fis ) @xcite , to cover the whole ir regime from 2 to 180  . far - ir spectroscopic capability of akari is realized by a fourier transform spectrometer ( fts ) as a function of the fis @xcite . the advantage of fis - fts over the iso / lws is an imaging spectroscopic capability with higher spatial resolution ; each pixel of the fis - fts detector arrays can simultaneously take a spectrum with about 40@xmath17 ( fwhm ) angular resolution , while the iso / lws has the single beam of 80@xmath17 in fwhm . the physical scales corresponding to their resolutions in the lmc are about 10 pc and 20 pc for fis - fts and the iso / lws , respectively . spectral maps with high spatial resolution are essential to investigate the interstellar media , especially those around ionized regions . the structure of an ionized region is determined by the locations of ionizing stars and the distribution of the ambient interstellar media . high spatial resolution maps of the far - ir lines clearly show the structure of star forming regions relating to the environment of star formation . spitzer space telescope made spectral line maps of 30 dor by using the irs and the mips in a sed mode in the spitzer legacy survey of the lmc : surveying the agents of a galaxy s evolution ( sage ; @xcite ) , which revealed the ionization state of the 30 dor region and super star cluster r136 @xcite . akari fis - fts observed many local regions distributed around the 30 dor region and provided their spectral maps . in this paper , we present the large - scale spatial distribution of the [ oiii ]  88 line emission around the 30 dor region .    the advantage of the akari [ oiii ] data over the isocam [ neiii ] and the spitzer data is the spatial coverage of much wider areas , i.e. , those far outside the bubbles around 30 dor . thus our data provide information on the ism influenced by intense radiation from massive stars rather than shocks by stellar winds and supernovae in r136 . the optical [ oiii ] map with mcels also covers a large area . as compared to the far - ir [ oiii ] map , the optical map has a higher spatial resolution , but suffers extinction effects in dense gas regions , which may hamper accurate measurements of the change of the line intensity with the distance from r136 . akari carried out far - ir spectroscopic observations from april 2006 to august 2007 . as a result , about 600 pointed observations were performed with fis - fts . more than two thirds of the fis - fts observations were parallel observations with the irc , in which the irc has priority for observations . since the fields - of - view ( fovs ) of the irc and the fis are separated by about half a degree , fis - fts observed off - primary target skies in the parallel observations . yet these data are useful as a spectroscopic survey of the diffuse interstellar medium , since many irc observations concentrate on the galactic plane and the lmc .    in this paper , we selected 14 pointed observations around the 30 dor regions , which targeted relatively bright regions among the fis - fts observations of the lmc . an observation log is summarized in table [ table : obs_log ] , and observed areas are shown in figure [ fig : obs_area ] , plotted on the spitzer / mips  24 map from the sage program . for comparison , the observational positions of the iso / lws are marked by the small open circles with a size of the beam , and the area of the spectral mapping with the spitzer / mips in a sed mode is indicated by the rectangle in the center of 30 dor . although our observations do not cover the center of 30 dor , they cover wide areas around the 30 dor region and the south ridge . as shown in figure [ fig : obs_area ] , our observational areas are located on the periphery of the area covered with spitzer and thus complementary to the spitzer observations . each observing point has 60 pixels with a fov of 30@xmath18 30@xmath17 for each pixel in the short wavelength band , and 45 pixels with a @xmath19 @xmath20 @xmath19 pixel fov in the long wavelength band . all the observational data analyzed in this paper are taken in a low resolution mode ( named sed mode ) , whose spectral resolution is about 1.2 cm@xmath4 . cccclc no . & date & position ( j2000.0 ) & obs . param . & comment & det . limit + 1 & 2006/10/14 - 21:15 & ( 84.74402 , -69.67834 ) & sed , 0.5s & irc04 ( b;ns ) & 3.9 + 2 & 2006/10/19 - 21:57 & ( 85.20209 , -69.57140 ) & sed , 0.5s & irc04 ( b;ns ) & 1.1 + 3 & 2006/10/21 - 21:54 & ( 85.13556 , -69.52066 ) & sed , 0.5s & irc02 ( b;n ) & 1.7 + 4 & 2006/10/23 - 18:32 & ( 84.85551 , -69.45006 ) & sed , 0.5s & irc02 ( b;n ) & 1.1 + 5 & 2006/10/23 - 20:11 & ( 85.16322 , -69.07068 ) & sed , 0.5s & irc04 ( b;ns ) & 3.4 + 6 & 2006/10/24 - 21:00 & ( 85.09940 , -69.17079 ) & sed , 0.5s & irc02 ( b;n ) & 1.1 + 7 & 2006/10/25 - 18:29 & ( 84.58520 , -69.37970 ) & sed , 0.5s & irc02 ( b;n ) & 1.2 + 8 & 2006/10/25 - 21:48 & ( 84.70576 , -69.23942 ) & sed , 0.5s & irc02 ( b;n ) & 1.0 + 9 & 2006/10/26 - 20:57 & ( 84.30852 , -69.30413 ) & sed , 0.5s & irc02 ( b;n ) & 0.9 + 10 & 2006/10/26 - 22:36 & ( 84.95287 , -68.95736 ) & sed , 0.5s & irc02 ( b;n ) & 1.1 + 11 & 2006/10/27 - 20:06 & ( 84.56963 , -69.02005 ) & sed , 0.5s & irc02 ( b;n ) & 1.4 + 12 & 2006/10/27 - 21:45 & ( 84.03696 , -69.22217 ) & sed , 0.5s & irc02 ( b;n ) & 1.2 + 13 & 2006/10/29 - 20:03 & ( 84.20209 , -69.08269 ) & sed , 0.5s & irc02 ( b;n ) & 1.0 + 14 & 2006/10/30 - 20:51 & ( 83.81157 , -69.14057 ) & sed , 0.5s & irc02 ( b;n ) & 1.3 +   +    ( 120mm,150mm)figure01.eps we used the standard reduction tools for fis - fts to derive spectra . details of the standard tools are described in the akari fts toolkit manual . here , we briefly explain the way to obtain spectra from raw signal and to extract emission lines from the spectra . raw data are processed through the standard pipeline tools , in which a linearity correction , glitch removing and a correction of the optical path difference ( opd ) are performed with an eye - ball check . then , a set of interferograms for each pixel is derived . interferograms are transformed to raw spectra by applying a discrete fourier transformation at the second stage of the pipeline , where the adopting opd range and the apodizing function are selectable . in this paper , we applied a box - type apodizing function and used data in an opd range of @xmath22  0.35  cm for the short wavelength band , a little shorter than the full range ( @xmath22  0.42  cm ) , to avoid a channel fringe @xcite . the resulting spectral resolution of the short wavelength band is reduced to @xmath23  1.4  cm@xmath4 from @xmath23  1.2  cm@xmath4 of the full performance for the flat apodization . in the final step of the pipeline , raw spectra are converted to sky spectra by correcting for instrumental factors . details of the calibration scheme of fis - fts and the accuracy are described in @xcite .    since 30 forward and backward scans are obtained by a single pointed observation , there is a set of 60 spectra for each pixel . from the set of the spectra , we calculated an average spectrum and estimated the error of each spectral element . figure [ fig : sw_spectrum ] shows examples of the averaged spectra obtained in the short wavelength band , exhibiting the clear detection of the [ oiii ]  88 line at 113  cm@xmath4 wavenumber . the structure around the line indicates a sinc function due to the box - type apodizing function . dust continuum emissions are also detected ; the sw26 and sw45 pixels show relatively flat spectra , while the continuum spectrum of the sw28 pixel indicates a cooler dust temperature . from these spectra , we extracted [ oiii ] line intensities by simply subtracting continuum emissions from the signals integrated over 111.0  114.5  cm@xmath4 . the contribution of the continuum emissions at the line position is estimated from signals in both sides of the line , 90.0  111.0  cm@xmath4 and 114.5  130.0  cm@xmath4 . the detection limit depends on the sensitivity of each pixel and also on the environmental condition of observations . the 3@xmath24 detection limits of the [ oiii ] line listed in table [ table : obs_log ] are median values of array pixels for each pointed observation estimated from errors at the [ oiii ] line . the minimal 3@xmath24 detection limit is about 0.6  @xmath20  10@xmath25  w  m@xmath3  sr@xmath4 , and a typical value is about 1  @xmath20  10@xmath25  w  m@xmath3  sr@xmath4 . ( 80mm,75mm)figure02.eps    we also obtain spectra in the long wavelength band , which include the [ cii ]  158 line emission detected in many observations . however , due to the non - uniformity of the detector performance among the pixels , the spectral data are still not well calibrated to discuss the spatial distribution of the [ cii ] line . therefore , in this paper , we treat only the [ oiii ]  88 line , which traces highly - ionized gas . the [ oiii ]  88 intensity map thus obtained is shown in figure [ fig : oiii_map ] together with the mips  24 intensity contour map . the [ oiii ] map is created by re - gridding pixel data on the astronomical coordinates and masking those with signal - to - noise ratios less than two . the mips  24 contour map in the figure is convolved with a gaussian kernel whose width match the spatial resolution of the fis ( @xmath26 ) . the map clearly shows widely distributed [ oiii ]  88 line emission with arc - minute scale structures around 30 dor . from iso archival data , we find that the iso / lws barely detects the [ oiii ]  88 line emission at the position in the fis - fts # 3 area ( see figure [ fig : obs_area ] ) with the intensity of about 2 @xmath20 10@xmath27 w m@xmath3 sr@xmath4 , which is below our typical 3@xmath24 detection limits as shown in table [ table : obs_log ] . the map in figure [ fig : oiii_map ] is consistent with the iso measurement . this line intensity level represents that from quiescent parts of the lmc and the intensities of the spatially extended emission seen in figure [ fig : oiii_map ] are @xmath28 times higher than this level . we significantly detect the [ oiii ]  88 emission at many pixels in each pointed observation . figure [ fig : oiii_map ] shows the global structure of the [ oiii ] emission in about one square degree ; the emission region of the [ oiii ] line is extended around the 30 dor region more widely than the area covered with the spitzer / mips in a sed mode . the structure of the [ oiii ] line emission , the highly - ionized region , is well correlated with that of mips  24 ( contours ) tracing hot dust emission around the 30 dor region . in particular , a filamentary structure at the northern part seen in the mips  24 contours is also recognized well in the [ oiii ] intensity map ( obs.#10 & # 11 ) . spitzer obtained spectra in the 30 dor region with the irs and the mips in a sed mode . the [ oiii ]  88 intensity map taken by the mips covers the center of 30 dor ( rectangle in figure [ fig : obs_area ] ) , partially overlapped with our map at the north edge . the mips measurement shows values about three times larger than ours , but both are roughly consistent by considering large uncertainties of our absolute line intensity calibration @xcite . owing to our clear detection outside the spitzer observing area , we reveal the widely extended structure of the highly - ionized region . in particular , the ionized regions are extended in the north and the west direction along the mips  24 structures , which have spatial correspondence to the x - ray and h@xmath11 filamentary structures .    on the other hand , the south - west region ( obs.#12 & # 14 ) outside the 30 dor does not have a clear correlation between the [ oiii ]  88 and the mips  24 emissions . we also detect the [ oiii ] emission significantly from this region , whose distribution is rather flat as compared to the mips  24 . this region contains the 30 dor c superbubble , which is a large ( 45 pc radius ) shell of material swept up by fast stellar winds and supernovae from tens of massive stars in an ob association . there is a mixture of thermal and nonthermal x - ray emission ( @xcite ; @xcite , @xcite ) ; any photoionization in this region may receive a contribution by hard radiation from the supernova remnants . there is probably also direct collisional ionization to o@xmath1 . the south region far outside 30 dor ( obs.#1  # 4 , # 7 & # 9 ) may have a different relation from the region near 30 dor . at least , the obs.#1 data are apparently influenced by n159 and n160 , other massive star forming regions ( see figure [ fig : obs_area ] ) . the mips  24 emission regions in the south are compact as compared to the 30 dor region . there are extended molecular clouds in this region and ionizing stars are located in dense cores , suggesting that star formation is in an early stage .    for clarity of presentation , we categorize the observational points into four groups according to their positions with respect to 30 dor . the nw group ( obs.#10 , # 11 and # 13 ) is located in the north - west region of 30 dor and has a good correlation with the mips  24 surface brightness . the sw group ( obs.#12 and # 14 ) is located on the far west side of 30 dor and has a rather flat distribution of the [ oiii ] line emission . the se group ( obs.#5 , # 6 and # 8) is located in the south - east region of 30 dor . the sr group ( obs.#1  # 4 , # 7 and # 9 ) is located along the south ridge of molecular clouds . figure [ fig : corr_ha ] shows a scatter plot between the [ oiii ]  88 line intensity and the h@xmath29 intensity , where the h@xmath30 data are taken from the southern h - alpha sky survey atlas @xcite . the different colors correspond to the data from the different group defined above . we can see a global correlation between the [ oiii ] and h@xmath29 emissions for a range of about one and a half orders of magnitude . the correlation , however , is not so tight at the higher line intensity region . a possible reason is difference in ionization degrees of the gases associated with the h@xmath29 and [ oiii ]  88 emission ; much higher energies of photons ( @xmath31 35.12 ev ) are required to ionize o@xmath16 to o@xmath1 than to ionize h atom . ( 80mm,75mm)figure04.eps    we also make a scatter plot between the continuum emission at 88   derived with fis - fts and the mips  24 emission , as shown in figure [ fig : corr_mips ] . the spatial resolution of the mips  24 map is reduced to match the fis - fts map . there is a tight correlation between the two continuum emissions at 24   and 88   for about two orders of magnitude . thus the 88   continuum signal detected with fis - fts reasonably traces dust emission . as shown in the plot , brighter regions exhibit a little bluer colors , indicating hotter dust temperatures , and thus stronger radiation field . all the data points in every group are on almost the same relation , suggesting that the properties of dust are similar throughout the 30 dor region . we reveal the extended distribution of the [ oiii ] emission around the 30 dor region ; the distribution is found to be much wider than the boundary of the previously observed area . even on the peripheral regions @xmath5 away from the center of 30 dor , we obtain [ oiii ] line intensities so high as @xmath2 w m@xmath3 sr@xmath4 . from calculation of a three - level system for the lowest energy levels of o@xmath1 @xcite , we estimate that the gas densities of @xmath6 @xmath7 are required to explain the above line intensities . here we assume that 50 % of oxygen is in the ionization state of o@xmath1 and the o / h abundance ratio is @xmath32 @xcite . for the depth of the emitting region , we adopt 75 pc ( = @xmath33 ) , which is similar to the distance from the center of 30 dor .    to ionize o@xmath16 to o@xmath1 , energetic photons ( @xmath31 35.12 ev ) are required , which can be produced only by very massive stars such as early o type stars and wolf - rayet ( wr ) stars . we estimate the contribution of the massive stars to the ionization of o@xmath16 . r136 has a lot of massive stars in the compact central region , as revealed by the hubble space telescope ( @xcite ; @xcite ) . the ionization structure at 30 dor ( a few arc - minutes around r136 ) is determined by a balance between radiation from massive stars in r136 and the distribution of the interstellar medium ( @xcite ; @xcite ) .    to assess the influence range of radiation field from the r136 super star cluster , we investigate the change of the [ oiii ] intensity with distance from r136 . figure [ fig : r_depend ] shows the radial profiles of the distributions of the [ oiii ] emission and the far - ir continuum emission plotted against the distance from r136 . the upper plots ( filled circles ) show the [ oiii ] line emission profile , while the lower plots ( open circles ) show the far - ir continuum emission profile . the sr group is excluded in this plot . the [ oiii ] emission and far - ir continuum of the nw and se groups show a clear dependence on the distance from r136 as seen in figure [ fig : r_depend ] , while those of the sw group do not ; the latter is likely to have local ionization sources other than r136 . ( 80mm,85mm)figure06.eps    figure [ fig : r_depend ] shows the wide distribution of the highly - ionized region around the central super star cluster r136 . the influence of r136 seems to reach @xmath34 , or 150 pc away from r136 . the intensities of the central [ oiii ] emission detected with the iso / lws are consistent with the extrapolation of the akari / fis - fts data to the inner area of @xmath10 from r136 , as shown in figure [ fig : r_depend ] . the iso / lws data point labeled n157b is located on the supernova remnant ngc2060 ( see figure [ fig : obs_area ] ) . it is notable that the far - ir continuum emission has a radial dependence similar to the [ oiii ] emission . the vertical line in the figure indicates the strmgren radius of hydrogen gas ( about @xmath35 ) determined by the r136 super star cluster . to calculate the value , we assume that the massive stars within a radius of 15 pc ( @xmath10 ) ( @xcite ; @xcite ; @xcite ) concentrate on the center and that ambient gas density and temperature are 100 @xmath7 and @xmath36 k , respectively . the gas density adopted here reflects the above value that is estimated from the [ oiii ] line intensity observed with fis - fts . these gas parameters are also reasonable from the results obtained with the radio @xcite , iso @xcite , and spitzer observations @xcite , which suggested densities around 200 @xmath7 near the center of 30 dor . we consider a conservative case that the gas density is constant at the value estimated in the peripheral region . yet , from figure [ fig : r_depend ] , it is clear that the simple calculation can not explain the widely - extended distribution of the [ oiii ] line emission . the central bubble of 1@xmath37 in diameter would increase the strmgren radius , but its presence does not dramatically change the result . to explain the observed distribution , we have to assume a lot of hidden massive stars , and/or non - uniform distribution of the interstellar media , that is , clumpy distribution and leakage of ionizing photons . it is not likely that a lot of embedded massive stars are hidden in the field . if the gas density is reduced down to 10 @xmath7 , the strmgren radius expands to @xmath38 , which is comparable to the observed size of the distribution of the [ oiii ] emission . therefore a picture of the clumpy distribution of the interstellar media with a volume filling factor of @xmath39 is a plausible explanation . thus it is difficult to explain the ionization in the areas far away from 30 dor solely by radiation from massive stars in r136 . however many other massive stars are also distributed widely in the 30 dor region in addition to those in r136 . @xcite presented a catalog of 1750 massive stars in the lmc with accurate spectral types ; stars in the r136 super star cluster are not included in the catalog . we calculated the strmgren radius of hydrogen gas for each cataloged star @xcite , assuming the physical condition of @xmath40 @xmath7 and @xmath36 k for gas density and temperature , respectively . the gas density again corresponds to the value obtained from the observed [ oiii ] line intensity . the strmgren radius thus obtained is plotted in figure [ fig : stromgren ] with the circle for each massive star . so far the strmgren radii are calculated for h@xmath16 ions , but not for o@xmath1 ions . however our simplified calculation shows that the strmgren radius for o@xmath1 is almost the same as for h@xmath16 as long as ionizing photons have energies significantly higher than the o@xmath16 ionization energy ( 35.12 ev ) . at least , a strmgren radius for o@xmath1 is no larger than that for h@xmath16 ions . thus it is difficult to explain the observed [ oiii ] emission in the outer region of 30 dor by radiation from the distributed massive stars , too . ( 120mm,150mm)figure07.eps    in order to evaluate the contribution of hidden sources of ionization to the observed distribution of the [ oiii ] line emission , in figure [ fig : qi_oiii ] , we plot the [ oiii ] line intensity as a function of the o@xmath16-ionizing photon flux expected from the cataloged stars . we obtain the flux at each position by adding the total absorption - free photon fluxes coming from the cataloged wr stars and o2o7 stars shown in figure [ fig : stromgren ] . here we use the projected distance for each star . as seen in the figure , there is a correlation between the [ oiii ] line intensity and the absorption - free ionizing flux , at least for the nw region and part of the se region . this suggests that the ism is very clumpy in those regions and that there are no significant hidden sources of ionization other than the dense star clusters in figure [ fig : stromgren ] . an exception is the sw region that contains 30 dor c , where the [ oiii ] line emission shows almost no dependence on the ionizing photon flux with line intensities above and below the correlation at smaller and larger photon fluxes , respectively . the lack of the correlation suggests the presence of another diffuse ionizing source , which can be attributed to hard radiation and/or shocks from the supernova remnants 30 dor c. as for the lower line intensities at photon larger fluxes , some areas of the sw region are so closely located near massive stars ( see figure [ fig : stromgren ] ) that the photon fluxes calculated from the projected distances are very likely overestimated . ( 80mm,75mm)figure08.eps    figure [ fig : corr_cnt ] shows a scatter plot between the [ oiii ] line and the 88   continuum emission , both obtained with fis - fts , where we see an overall correlation between them . it is consistent with the similar trend of radial dependence between the [ oiii ] emission and the continuum emission as seen in figure [ fig : r_depend ] . this correlation is interesting because they are likely of different origins ; the far - ir continuum represents the dust emission mostly coming from neutral regions , while the [ oiii ] emission originates in highly - ionized gas regions . thus the correlation suggests that the highly - ionized gas and dust are mixed well and exist in the same region , or have a clumpy structure where the dust is shielded from energetic photons by local dense gases . this interpretation is consistent with the widely extended distribution of the [ oiii ] emission . ( 80mm,75mm)figure09.eps    the ionization of the sw region does not seem to be related to the central sources in the 30 dor region as seen in figure [ fig : r_depend ] . in figure [ fig : corr_cnt ] , all the data points seem to follow the global trend of correlation . by looking locally in detail , however , the sw group has a much weaker dependence of the [ oiii ] intensity on the far - ir continuum emission , i.e. no clear correlation between the [ oiii ] and the dust emission , suggesting that the dust and ionized gases are not mixed well in this region . as described above , most data points of the sw region are located on the supershell of the supernova remnants 30 dor c @xcite . then , a large fraction of dust in the supershell may have been destroyed by the supernovae . we have performed spectroscopic surveys of the lmc with akari / fis - fts with a good spatial resolution of about @xmath10 . we have analyzed spectral data from the 14 observational points distributed widely around the 30 dor region on a degree scale . we find that the [ oiii ]  88 emission is widely extended by more than @xmath0 ( 150 pc ) from the central super star cluster , r136 . since the [ oiii ] line traces highly - ionized gas regions heated by very massive stars , the detection of widely extended [ oiii ] emission suggests that energetic photons are escaping into regions far from massive stars localized at the center of the 30 dor region . it is unlikely that a lot of embedded massive stars are hidden in the outer region of 30 dor . therefore , a plausible interpretation is that the interstellar medium has a patchy or porous structure around the 30 dor region , so that high - energy photons can pervade into wider areas . this is also supported by a good correlation between the [ oiii ] emission and the far - ir continuum emission , suggesting that the gas and dust are well mixed in the highly - ionized region where the dust survives in clumpy dense clouds shielded from the energetic photons . we would express many thanks to an anonymous referee for giving us important comments . we thanks all members of the akari project . particularly , we express our gratitude to the akari data reduction team for their dedicated work in generating the tsd , and developing data analysis pipelines . akari is a jaxa project with the participation of esa . 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 . this work also made use of data product from the southern h - alpha sky survey atlas ( shassa ) , which is supported by the national science foundation . murakami , h. et al . 2007 , , 59 , s369-s376 murakami , n. et al . 2010 , , 62 , 1155 - 1166 onaka , t. et al .  2007 , , 59 , s401-s410 peck , a. b. , goss , w. m. , dickel , h. r. , roelfsema , p. r. , kesteven , m. j. , dickel , j. r. , milne , d. k. , & points , s. d.  1997 , , 486 , 329 - 337          smith , r. c. , leiton , r. , & pizarro , s.  2000 , in stars , gas and dust in galaxies : exploring the links , asp conference proceedings , vol . 221 . alloin , d. , olsen , k. , and galaz , g. , san francisco , 83 - 86 .
Abstract:

Answer: we present the distribution map of the far - infrared [ oiii ]  88 line emission around the 30 doradus ( 30 dor ) region in the large magellanic cloud obtained with the fourier transform spectrometer of the far - infrared surveyor onboard akari . the map reveals that the [ oiii ] emission is widely distributed by more than @xmath0 around the super star cluster r136 , implying that the 30 dor region is affluent with interstellar radiation field hard enough to ionize o@xmath1 . the observed [ oiii ] line intensities are as high as @xmath2 w m@xmath3 sr@xmath4 on the peripheral regions @xmath5 away from the center of 30 dor , which requires gas densities of @xmath6 @xmath7 . however the observed size of the distribution of the [ oiii ] emission is too large to be explained by massive stars in the 30 dor region enshrouded by clouds with the constant gas density of @xmath8 @xmath7 . therefore the surrounding structure is likely to be highly clumpy . we also find a global correlation between the [ oiii ] and the far - infrared continuum emission , suggesting that the gas and dust are well mixed in the highly - ionized region where the dust survives in clumpy dense clouds shielded from the energetic photons .


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 rapid pace of extrasolar planet discovery in recent years has provided detailed statistical data about the distribution of planets in our galaxy . one of the many surprises in these data is the large number of planets with significant eccentricity measured through the radial velocity ( rv ) method @xcite , particularly in those planets orbiting beyond @xmath0 au , where tidal damping is inefficient . in contrast , the planets in our own solar system ( aside from mercury ) have low - eccentricity orbits . while eccentricity has been primarily measured through rv , eccentricity can also be inferred for the many exoplanets discovered through the transit method by utilizing transit timing variations ( ttvs ) @xcite and through the so - called ` photoeccentric effect ' on the transit lightcurve @xcite . @xcite recently showed through a comparison of orbital and transit timescales that the statistical eccentricity distribution of the planets detected by the _ kepler _ spacecraft may be similar to the distributions of those that have direct rv measurements of their eccentricity , however @xcite point out that this depends on the uncertain stellar radii given for the kepler sample , particularly for smaller stars . @xcite also showed that the mean eccentricity of _ kepler _ planet candidates increases with increasing planet size plateauing above neptune - sized planets ( @xmath2 ) . the planets are thought to have migrated into their present orbits . migration processes fall into two broad classes , planet - disk interactions @xcite , where planets migrate smoothly through the disk ; and few - body interactions , where planets are moved inwards through interaction with perturbing or scattering bodies ( see e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? multi - body dynamical processes have been studied intensely , in part to explain the presence of hot jupiter systems that were found by early rv observations . however , nearly a third of the likely planet candidates found by the _ kepler _ observatory are thought to be in co - planar multi - planet systems @xcite , indicating that their evolution was likely dominated by disk interaction . the results of disk migration also provide the initial conditions of the planet locations and eccentricity distribution for the multi - body interactions that take place after the disk has dissipated . given an initial distribution of sufficiently non - zero eccentricities ( @xmath3 ) , @xcite showed that dynamical relaxation of the planetary systems could explain the eccentricity distribution of higher eccentricity systems . migration driven by planet - disk interactions has been investigated as a means to excite eccentricity in migrating planets . disk resonances exchange both energy and angular momentum between a planet and the disk , modifying the semi - major axis and the eccentricity of the planet . certain resonances tend to pump eccentricity , while others tend to damp it . it is generally thought that the net effect of disk - planet interactions is to damp eccentricity . for small planets , this is almost certainly the case @xcite . for giant planets , @xcite and @xcite suggested that saturation of the eccentricity - damping first - order corotation resonances could lead to eccentricity growth of planets able to clear a gap in the disk . however , recent numerical studies @xcite have not seen such saturation in simulations , finding that eccentricity is damped for giant planets except for extreme cases where the gap is wide enough to avoid the corotation resonances entirely . entropy gradients in the disk can modify corotation resonances through baroclinic effects . previous studies on the effects of entropy gradients in non - barotropic disks have focussed on the primary ( co - orbital ) co - rotation resonance , numerically solving for the response of the disk @xcite , and the effect of the non - barotropic non - linear horseshoe torques on type i migration of small planets embedded in the protoplanetary disk . these non - barotropic effects can allow halting or even reversal of migration for certain systems . @xcite recently utilized the numerical fitting formulae from @xcite to predict the effect of stellar illumination on migration of small embedded planets .    in a companion paper @xcite , we go beyond these previous studies by deriving a fully analytic form for the linear non - barotropic corotation torques , for both primary and first - order corotation resonances . here , we apply these results to eccentricity excitation in disks in which the gap formed by a migrating planet is heated by stellar irradiation . the general picture is as follows . as a gap clearing planet perturbs the disk , lindblad resonances excite the planet s eccentricity while non co - orbital corotation torques damp its eccentricity . in a barotropic disk , the damping dominates unless the gap is very large . we show here that the entropy gradients at the edge of an irradiated gap act to reduce the non - barotropic corotation torque , allowing the eccentricity - exciting lindblad resonances to dominate . the requirement for the gap to be irradiated in order to enable eccentricity growth in the disk offers a potential explanation for the lack of eccentric planets at intermediate orbital periods around low metallicity stars found by @xcite . in low metallicity disks , heating of the disk at the dust sublimation radius may lead to shadowing of the disk in the region @xmath0@xmath4 @xcite , causing eccentricity damping in this `` eccentricity valley '' region . an outline of the paper is as follows . in section 2 we discuss the nature of non - barotropic corotation torques and their dependence on background entropy gradient . in section 3 we review the mechanism of eccentricity damping for giant planets in barotropic disks . in section 4 we show that for reasonable assumptions about the disk structure and insolation the resulting non - barotropic torques can lead to eccentricity excitation . in section 5 , we demonstrate how the `` eccentricity valley '' may be explained by eccentricity damping and excitation in self - shadowed and fully irradiated disks respectively . finally we discuss the implication of our results in section 6 . a quantity that plays an important role in the determination of the corotation torque is the vortensity , @xmath5 , ( also known as the potential or specific vorticity ) where @xmath6 is the z - component of a fluid s vorticity . @xcite showed that the torque associated with the corotation resonance has the sign of the radial vortensity gradient in the disk . for disks with orbital angular velocity @xmath7 , @xmath8 can be written as = = , where @xmath9 is the square of the radial epicyclic frequency , @xmath10 is the local oort constant , and @xmath11 is the disk surface density .    the continuity equation and the curl of the momentum equation for ideal fluids can be combined to obtain the lagrangian ( co - moving ) time derivative for vortensity ( the vorticity equation , see e.g. * ? ? ? [ vorticityeqn ] the right hand side of the equation above is the baroclinic term , which contributes when the lines of constant ( surface ) density @xmath11 and ( vertically integrated ) pressure @xmath12 are not parallel in a fluid .    in most previous works studying eccentricity evolution due to disk - planet interaction , the disk is assumed to be barotropic , such that pressure is only a function of density , @xmath13 . this causes the lines of constant @xmath12 and @xmath11 to be parallel for any perturbations , such that the baroclinic term is zero . thus for barotropic disks , the vortensity is conserved for all perturbations .    in a non - barotropic disk , where entropy gradients exist , the baroclinic term is non - zero . for an ideal gas with adiabatic index @xmath14 where @xmath15 , we can write equation as = , [ vorticity2 ] which shows that density perturbations may generate vortensity in the presence of an entropy gradient . a planet orbiting with non - zero eccentricity has a time varying gravitational potential which can be split into the fourier components _ p(r ) = _ l , m _ l , m , where @xmath16 is the pattern frequency for the ( @xmath17 ) component , and @xmath18 and @xmath19 are the orbital and epicyclic frequencies of the planet . only the @xmath20 ( primary ) and @xmath21 ( first - order ) , components are non - zero to first order in eccentricity .    at lindblad resonances , where the pattern frequency of the gravitational perturbation in the frame co - moving with the disk material is equal to the local epicyclic frequency ( @xmath22 ) , the disk material is able to respond strongly to such perturbations and waves are launched , resulting in lindblad torques .    at a corotation resonance , the disk material moves with the same speed as the pattern speed for a particular potential component . the primary corotation resonances ( where @xmath23 ) are located at the same orbital radius as the planet , while the first - order corotation resonances ( where @xmath24 ) are located both interior and exterior to the orbit of the planet ( see figure [ resonancelocation2 ] ) . in the remainder of this paper we will use the terms `` first - order corotation resonance '' and `` non co - orbital corotation resonance '' interchangeably . a full calculation of the linear non - barotropic corotation torque must also consider the disk enthalpy response to the perturbing potential . this has been previously studied for co - orbital ( primary ) corotation resonances by numerically calculating this response @xcite . in the companion paper to this work , @xcite , we employ an approach similar to that of @xcite to derive a fully analytic form for the linear non - barotropic corotation torque for both primary and first - order corotation resonances . while the full expression from @xcite should be used to evaluate such torques in the general case , for the purposes of this discussion we will adopt the simplifying limit of a cold ideal gas disk , such that the adiabatic sound speed , @xmath25 , is small compared to the orbital velocity ( such that the disk scale height , @xmath26 , is small ) . we allow the disk to have entropy gradients such that there can be a non - zero radial brunt - visla frequency , @xmath27 , defined to be n_r^2 - . for simplicity we also assume that @xmath28 such that the disk is stably stratified radially , though this may not be the case for some disks , particularly those where the baroclinic instability only occurs slowly @xcite , or where magnetic fields have prevented instability from occurring @xcite . for the gaps examined later in the paper @xmath28 near the gap edges , as both the pressure and entropy gradients change sign .    for @xmath29 the expression for the corotation torque of the disk acting on the planet can be given as @xcite , which is not unreasonable for the models we consider later . for sharper entropy profiles or thicker disks the full expression of equation 47 in @xcite should be used . in sharp gaps where @xmath30 the 3-dimensional torques should be considered ( see e.g. * ? ? ? * ; * ? ? ? ] ^(cr , cold)_l , m -_r_c _ eff + o(n_r^4/^4 ) [ torqueeqn ] , where the `` effective '' vortensity is given by _ eff = /(s^1/ ) = , and @xmath31 is the radial location of the corotation resonance . the sign of the corotation torque for a cold disk is thus controlled by the combination of the radial gradients of the vortensity and the disk entropy . in the limit where @xmath32 and @xmath33 constant , we recover the original expression for barotropic the corotation torque from @xcite ( see their eq . [ 56 ] ) , where the sign of the torque depends only on the gradient of the vortensity . note that if the characteristic timescale of energy transport in the disk is shorter than the perturbation ( orbital ) timescale , then the perturbations can not be treated as adiabatic , but will be isothermal , and the barotropic torque will also be recovered .      ) , inner lindblad ( @xmath34 ) and outer lindblad ( @xmath35 ) resonances are shown relative to the gap - clearing planet for the @xmath36 primary ( @xmath20 ) and first order ( @xmath21 ) potential components . also shown are the dominant contributions to the eccentricity evolution assuming a cleared gap . the eccentricity damping by the non co - orbital corotation torques tends to slightly dominate ( by @xmath37 ) the eccentricity pumping by the first - order external lindblad resonances @xcite leading to a net eccentricity damping by the disk - planet interaction . [ resonancelocation2 ] ]    the nature of the corotation torque can be understood qualitatively by considering the effect of the potential components on collisionless particles near the corotation point @xcite . just interior to the corotation resonance a particle feels a more slowly moving potential which slows it down and moves it towards the corotation radius . just outside the corotation a particle feels a more quickly moving potential , which speeds it up and moves it towards corotation . there is therefore a net particle flux towards the corotation , with the particle starting interior to the corotation gaining angular momentum and those starting exterior to the corotation losing it . the net torque on disk ( and conversely on the planet ) then depends on the relative flux of material into corotation from either side .    for barotropic perturbations the vortensity is conserved by equation , such that the vortensity is carried with fluid elements regardless of the torques exerted by the planet . in order to have the net vortensity of the incoming material match the vortensity at the corotation the relative flux towards the corotation from either side must scale with the background vortensity gradient @xmath38 , and thus the net torque at a corotation resonance must be proportional to this value for barotropic disks . however , for disks where an entropy gradient exists , such density and pressure perturbations can generate vortensity . the flux of material towards the corotation from either side needed to match the net vortensity at the corotation must be altered by the baroclinic term of the vortensity equation , thus the net torque is modified to include the effect of the entropy gradient .    ) , inner lindblad ( @xmath34 ) and outer lindblad ( @xmath35 ) resonances are shown relative to the gap - clearing planet for the @xmath36 primary ( @xmath20 ) and first order ( @xmath21 ) potential components . also shown are the dominant contributions to the eccentricity evolution assuming a cleared gap . the presence of an entropy gradient at the non co - orbital ( first order ) corotation resonances can cause a reduction ( or even a reversal ) of the corotation torque compared to the barotropic case ( see figure [ resonancelocation2 ] ) , allowing the first - order external lindblad resonances to dominate and thus a net eccentricity excitation by the disk - planet interaction . [ resonanceentropy2 ] ] planets exchange both energy and angular momentum with the disk through the resonances at different rates which allows eccentricity evolution of the orbit to occur @xcite . the eccentricity of a planet , @xmath39 , can be given in terms of its energy @xmath40 and its angular momentum @xmath41 , e_p^2 = 1 + . the transfer of energy and angular momentum between the planet and the disk at any resonance occurs due to interaction of the potential component @xmath36 , with the disk material . this potential component is constant in frame corotating with @xmath42 , so that the jacobi constant @xmath43 of the disk material is conserved , where @xmath44 and @xmath45 are the energy and angular momentum of the disk . consequently the resonance transfers energy and angular momentum between the disk and the planet such that = _ l , m = _ l , m _ l , m , where @xmath46 is the torque acting on the planet . combining these equations to lowest order in eccentricity we find @xcite , e_p = ( _ l , m - _ p ) [ ecceqn ] .    external lindblad torques ( those not coincident with the planet location , see fig 1 ) tend to drive eccentricity growth ( and migration ) , while co - orbital lindblad torques damp the eccentricity of the planet @xcite . the co - orbital lindblad resonances tend to dominate for small planets by a factor of @xmath47 @xcite , thus disk - interaction tends to damp eccentricity of small embedded planets . if a planet is massive enough to clear a gap , such that the density of the co - orbital material falls dramatically , the damping by the co - orbital lindblad resonances is reduced . @xcite then found that in a barotropic disk , for small eccentricity , the dominant terms are the eccentricity excitation by the first - order external lindblad resonances , and the eccentricity damping by the non co - orbital corotation torques ( see figure [ resonancelocation2 ] ) , due to the vortensity ( density ) gradients in the gap .    for first - order corotation resonances located at the gap edge _ interior _ to the planet orbit , the positive vortensity gradient ( negative density gradient ) leads to a negative corotation torque on the planet ( @xmath48 ) , in the barotropic limit of equation . since @xmath49 interior to the planet location , equation yields eccentricity damping .    for first - order corotation resonances located at the gap edge _ exterior _ to the planet orbit , the negative vortensity gradient ( positive density gradient ) leads to a positive torque ( @xmath50 ) . however , @xmath51 exterior to the planet s orbit , and eccentricity is again damped . @xcite found that the full unsaturated corotation torques tended to dominate by a small amount ( @xmath37 ) leading to a net eccentricity damping , but along with @xcite noted that only a small amount of saturation of the corotation was necessary to allow net eccentricity excitation to occur . whether this condition for eccentricity growth is met depends in detail on the non - linear interactions that dictate the corotation saturation . numerical simulations ( see e.g. * ? ? ? * and references therein ) have shown that for some disk and planet parameters eccentricity growth can be seen . in particular for extremely large planets where a very large gap has been cleared , such that no material is present to interact at the non co - orbital corotation resonances , eccentricity growth driven by the first - order lindblad resonances is seen @xcite . however , @xcite recently performed a series of high - resolution 3-d smooth particle hydrodynamic simulations for various planet and disk parameters , and found that for realistic planet mass and disk properties planet eccentricity tends to be damped ( aside from the systems with extremely large gaps discussed above ) , implying that insufficient saturation of the non co - orbital corotation resonances occurs in realistic systems to excite eccentricity . this led @xcite and others ( see e.g. * ? ? ? * ) to conclude that disk interactions can not be a significant source of the eccentricity currently observed in many exoplanets . @xcite recently performed several simulations calculating the damping rates for eccentricity and inclination . we note , however , that all of these simulations were performed using a locally isothermal equation of state and without considering the effect of stellar illumination on the thermal structure of gap , i.e.  the gaps were effectively barotropic . stellar illumination of the disk can significantly alter its geometrical structure and temperature profile ( see e.g. * ? ? ? in particular , @xcite , @xcite and @xcite have used monte carlo and semi - analytic radiative transfer codes to show how insolation modifies the thermal structure around a gap opened by a planet . all three of these works treat radiative transfer in an axisymmetric model disk . in the latter two , the hydrostatic vertical structure is also iterated to consistency with the radiation field . if the disk has a flaring shape and the gap is optically - thin , the star directly illuminates and heats a strip of material along the top of the gap s outer edge . the outer edge in turn emits infrared photons , some of which travel down into the gap , lifting temperatures there above the equilibrium levels found with no gap . the gap s temperature is nearly uniform and determined by the balance between the received infrared photons and the radiative cooling , proportional to the planck - mean opacity of the material in the gap . the temperature profile of a gap depends on the geometry of the disk , in particular the angle of incident stellar illumination , and the amount of luminosity intercepted by the outer gap edge . as a specific example , figure [ tsigmaplot ] shows the surface density and midplane temperature profiles for two models with a jupiter - mass planet 5 au from a solar - mass star . the parameters are identical to the viscous model disk and 4600  k star from @xcite , but we repeated the calculations increasing the number of photon packets ten - fold to @xmath52 on the fifth iteration , to ensure that the temperature gradients at the gap edges are adequately sampled . the gap structure was obtained using the prescription from @xcite , with viscous stress parameter @xmath53 , planetary accretion efficiency @xmath54 , gap width parameter @xmath55 and planet / star mass ratio @xmath56 . the luminosities in the two models span the range expected between the beginning ( @xmath57 ) and end ( @xmath58 ) of the planet - forming epoch . ( blue ; in units of g @xmath59 ) and temperature @xmath60 ( red ; in units of k ) profiles for a gap - clearing @xmath61 planet located 5 au from a @xmath62 star with luminosity @xmath58 ( solid ) or @xmath57 ( dashed ) from the monte carlo radiative transfer calculations of @xcite . the green curves show the ratio of the characteristic density and temperature lengthscales , @xmath63 . dotted lines are also shown corresponding to torque ratio @xmath64 of various values , assuming adiabatic index @xmath65 , for a diatomic ideal gas . we see that while the entropy gradients at the gap edges are insufficient to reverse the sign of the corotation torques ( @xmath66 remains greater than zero ) , they are enough to significantly reduce the torques , with @xmath67 in the region of the gap edges . the vertical dashed lines show the locations of the @xmath68 non co - orbital corotation resonances , assuming a keplerian rotation profile in the disk . the fluctuations in the logarithmic derivative ratio at and inside the interior m=2 resonance seemingly imply @xmath69 ( and thus @xmath70 ) , however , these fluctuations are due to monte carlo noise from the low packet counts in the previous steps of the model iteration , and likely to shrink if more packets were used . [ tsigmaplot ] ]      in a non - barotropic disk , the presence of an entropy gradient at the gap edge can modify the balance of eccentricity excitation and damping ( see figure [ resonanceentropy2 ] ) . in equation the corotation torques are modified by the presence of an entropy gradient with the torque proportional to the logarithmic derivative of the effective vortensity , = ( ) - - . the torque due to an isentropic ( @xmath71 ) disk background is identical to the torque for a barotropic disk . for an ideal gas disk with adiabatic index @xmath14 and temperature @xmath60 , we can write the logarithmic entropy gradient of the background as = - ( - 1 ) , giving = - - . for a gap edge we can expect @xmath72 , which allows us to estimate the ratio of the non - barotropic and barotropic torques as , [ eq : chi ] + . in an isentropic disk , @xmath73 is constant such that the temperature and density gradients must be related by ( ) _ s = - 1 , and @xmath74 , and the torque is the same as the barotropic value . however as we now show , in the case of a heated gap @xmath66 can be significantly lower than 1 . @xcite estimated that in a barotropic disk the eccentricity damping by the non co - orbital corotation torques exceeds the eccentricity pumping by first - order lindblad torques by @xmath37 . the lindblad resonances are not significantly affected by the presence of entropy gradients as their locations are only perturbed slightly so that the perturbation frequency in the co - moving frame , @xmath75 , matches the epicyclic frequency modified to include the effect of the radial buoyancy . if we take @xmath76 to be the corotation torque reduction factor necessary to allow eccentricity driving by the lindblad resonances to win , then we arrive at the condition for eccentricity driving by the disk , < _ crit - 1 . [ critheat ] when the radial temperature scale height is much larger than the radial density scale height ( @xmath77 ) , as is the case in the almost isothermal heated gap , we expect a reduction in the torque by a factor of @xmath78 . an even smaller torque can result when the temperature gradient is of the opposite sign to the density gradient ( if there is a strong temperature inversion going into the gap ) . this result assumes that there is sufficient thermal diffusion to prevent the torque from thermally saturating . this is discussed further in @xcite , where it is shown that the dimensionless diffusivity required at the location of the first - order corotation resonances is small . the minimum dimensionless diffusivity required to maintain linearity is d_e > 10 ^ -7 ( ) ^3/2 , for the model gap width and location used in @xcite , where @xmath79 is the thermal conductivity , @xmath80 is the specific heat at constant pressure , @xmath81 is the volume density , and all quantities are evaluated at the resonance location . if the diffusivity is too small , then the gradient of the entropy perturbation becomes comparable to the background entropy gradient , and the torque thermally saturates , returning to its barotropic value . the minimum thermal diffusivity @xmath82 is exceeded by several orders of magnitude near the gap we consider here .    to show a specific example , we have plotted the logarithmic derivative ratio for the gap models discussed in  4.1 in figure [ tsigmaplot ] , where we see that the insolation into the gap heated it enough to provide @xmath67 at the gap edge . in fact , the heating of the gap has reduced the @xmath83 term in equation ( [ eq : chi ] ) so much that the value of @xmath66 is close to @xmath84 . this level of gap heating can be achieved for systems where the gap width and disk flaring angle allow the stellar luminosity absorbed by the outer edge of the gap to be at least comparable to the heating due to accretion . _ thus we expect that many giant planets experience eccentricity excitation by the disk when they are relatively close to the star and have their gaps exposed to insolation _ , in contrast to the results of previous barotropic simulations ( e.g. * ? ? ? if the light from the central star is blocked , or the angle of incident starlight is sufficiently small ( in , for example , a disk where the dust has significantly settled ) , then the entropy gradients will be much less steep , and the corotation torque should be closer to the barotropic level for which the planets will experience eccentricity damping due to the disk .    the gaps involved must be sufficiently clean to suppress the damping of eccentricity by the co - orbital lindblad resonances . defining @xmath85 , @xmath86 and @xmath87 to be the relative strengths of the non co - orbital corotation , external lindblad , and co - orbital lindblad resonances respectively ( c.f . * ) , we have @xmath88 @xcite . assuming a @xmath37 net eccentricity driving from the external lindblad and corotation resonances , such that @xmath89 , material in the gap will lead to eccentricity damping unless the gap is sufficiently clear @xmath90 . it should be noted that the final value of a planet s eccentricity can not be determined through the linear analyses described above if the planet experiences eccentricity driving through the disk . non - linear effects quickly become important as the eccentricity grows large and the planet nears the gap edge . however we note that even in barotropic simulations , once very wide gaps made the first - order corotation torques ineffective , the eccentricity quickly reached large values @xcite . also , we assumed the disk is more massive than the planet . movement of material within the disk could limit the planet s final eccentricity in cases where disk and planet masses are comparable . @xcite recently showed that there is a puzzling deficit of eccentric planets with semi - major axis between 0.1 and 1 au around metal - poor host stars ( see figure [ eccdata1 ] for more details ) . by assuming that eccentricity can only be a result of a strong dynamical interaction from a third body , and noting the well known relationship between giant planet occurrence and metallicity ( see e.g. * ? ? ? * ) , they conclude that this metallicity dependent `` eccentricity valley '' is evidence that only systems with many giant planets can produce eccentric planets between 0.1 and 1 au . however , as we have shown above , while a giant planet is embedded in the protoplanetary disk its eccentricity can be driven by interaction with the disk , provided its gap is sufficiently heated by the stellar illumination . this provides a means of exciting eccentricities even during the disk phase of the planet s orbital evolution , particularly for planets close to their host stars , as is the case for the bulk of the rv sample , and could provide the source of eccentricity for planets around metal rich stars . below we suggest a possible mechanism for preventing this from occurring in metal poor stars located in the `` eccentricity valley '' region , through disk self - shadowing . the progenitors of the solar mass stars studied by the rv sample in figure [ eccdata1 ] , are thought to be the classical t tauri stars , which possess accretion disks from which planets will likely emerge . @xcite showed that in a significant fraction of such stars a near - infrared excess exists , which is fit well by a single blackbody at the dust sublimation temperature . this was found to be consistent with a dust rim @xcite located at the dust sublimation radius . fitting this model , @xcite found that if the gas interior to the dust sublimation radius is optically thin @xcite , then the disk at the sublimation radius , located at @xmath91 , where @xmath92 is the dust sublimation temperature , becomes puffed up as it receives insolation directly from the star , heating it far more than would be the case in a normal flared disk . this raised optically thick rim blocks the starlight , and causes a shadowed region to form , where there is no direct stellar illumination of the disk , as shown schematically in figure [ diskschematic ] . the disk is again illuminated by the star at @xmath93 @xcite where it begins to flare , yielding a range for the shadowed region of @xmath94 au , where embedded planets would have their eccentricity damped . the gas interior to the dust sublimation radius has opacity dominated by molecular lines ( e.g. h@xmath95o ) between @xmath96k - @xmath97 k @xcite . thus for higher metallicities ( or accretion rates ) the gas interior to the dust sublimation radius can become optically thick , preventing the direct insolation and puffing - up of the dust rim . the disk should then have a normal flared shape @xcite , which `` sees '' the star at all radii and has no shadowed region . a sufficiently flared disk would drive eccentricity growth for giant gap - forming planets embedded in the `` eccentricity valley '' region .    ) becomes puffed up as it receives insolation directly from the star , heating it far more than would be the case in a normal flared disk . this raised rim causes a shadowed region to form , where there is no direct stellar illumination of the disk . a giant planet located in this region does not have its gap heated by the central star . @xcite showed that for t tauri disks , the dust rim radius should be @xmath98 au . the shadowed region should extend to the radius where insolation causes the disk to flare , @xmath99au @xcite . _ below : _ when the metallicity ( or accretion rate ) is sufficiently high , the gas interior to the dust sublimation region is optically thick , and the disk has a ` normal ' flared shape @xcite , which has gap heating at all radii.[diskschematic ] ] we have shown that entropy gradients and a non - barotropic equation of state modify the non co - orbital corotation torques for planets embedded in their natal disks . we have provided a simplified form for the torque valid in the limit of a cold disk , with @xmath100 . a more general analytic formula for the corotation torque is provided in a companion paper @xcite . if starlight heats the gap around a giant planet enough to satisfy equation , then the corotation torques are reduced to the point where the planet s eccentricity is excited rather than damped . this threshold is reached if enough of the stellar luminosity is absorbed on the gap s outer rim . we have shown that the threshold is easily exceeded in the flared structure obtained for a standard viscous disk model by placing the material surrounding the gap in radiative balance with the starlight . our starlight transfer calculations are axisymmetric and thus do not include the planet s wakes , but do capture the gross shape of the disk surface under vertical hydrostatic balance . thus we expect that starlight heating could have enabled many planets to gain eccentricity through disk interaction , rather than solely through later post - disk dynamical interactions as previously thought . note that eccentricity driving requires the gap around the giant planet to be sufficiently cleared that the co - orbital lindblad torques are suppressed , since these would dominate the eccentricity evolution for an embedded planet ( see e.g. * ? ? ? * for further discussion ) . the eccentricity driving has important implications for the inferred dynamical histories of giant planets in the currently known rv sample , and for the relative importance of orbital evolution due to the disk versus strong dynamical interactions , such as planet - planet scattering or the kozai mechanism . such disk interactions may be the source of the isolated eccentric giant planets . for the _ kepler _ sample , the inferred mean eccentricity dependence on the planet size @xcite , rising with @xmath101 and plateauing above @xmath102 , may be consistent with an eccentricity growth process being tied to gap clearing conditions . the eccentricity distribution at the end of the protoplanetary disk phase also provides the initial conditions for multi - planet interactions as discussed by @xcite , and further may help to explain the excess of planets with low , but non - zero eccentricity that are not explained by these interactions . we have suggested that the `` eccentricity valley '' for planets between @xmath0 and @xmath1 au around low - metallicity host stars @xcite may be explained by noting that this range also corresponds to the region shadowed by the dust rim that is formed when an optically thin gas interior to the dust sublimation point allows direct insolation of the dust wall at the sublimation radius . in the self - shadowed region between @xmath0 and @xmath103 au the disk interactions should damp eccentricity , resulting in circular orbits . for higher metallicities , and higher accretion rates , the flow interior to the dust sublimation radius can become optically thick , preventing direct insolation from puffing up the dust rim , and causing the shadowed region to disappear . such fully - illuminated disks may produce eccentric planets in the `` eccentricity valley '' region , as is seen for systems with [ fe / h ] @xmath104 . one might expect the extent of the shadowed region to decrease smoothly with increasing metallicity , as the height of the disk rim varies . intriguingly , a hint of this may be present in the lower envelope of the planet distribution shown in figure [ eccdata1 ] . we have only studied the linear growth of eccentricity in a system with a heated gap . to fully understand the eccentricity evolution of planetary systems and follow the eccentricity growth into the non - linear regime , detailed numerical simulations are needed that include non - barotropic equations of state and the effect of stellar illumination on the gap . additional modeling of the dust rim and the shadowed region of protoplanetary disks , and the geometric effects of a planetary gap in these regions is also needed to confirm if the `` eccentricity valley '' arises due to self - shadowing and non - barotropic torques in fully - illuminated disks . note that accretion heating is neglected in the two radiative transfer and hydrostatic balance calculations shown in figure [ tsigmaplot ] , as is shock heating associated with the planet s wake . starlight will modify the corotation torques only when the gap is heated by illumination at least as much as the gap edges are heated by accretion and shocks . this would depend on the accretion rate at the end of the disk phase , and the gap opening dynamics and resulting structure . it will also be important to study the expected entropy gradients in disks that are not irradiated or in regions of irradiated disks in which accretion heating dominates . dt was supported by funding from the lorne trottier chair in astrophysics and cosmology , and the canadian institute for advanced research . dt would like to thank kostas gourgouliatos , philip muirhead , jason wright , john johnson , and particularly rebekah dawson for helpful advice and useful discussions during the course of this work . njt carried out his part of the research at the jet propulsion laboratory , california institute of technology , under a contract with the national aeronautics and space administration and with the support of the nasa origins of solar systems program under grant 11-oss11 - 0074 . ac was supported by an nserc discovery grant and is an associate member of the cifar cosmology and gravity program .
Abstract:

Answer: we show that the first order ( non co - orbital ) corotation torques are significantly modified by entropy gradients in a non - barotropic protoplanetary disk . such non - barotropic torques can dramatically alter the balance that , for barotropic cases , results in the net eccentricity damping for giant gap - clearing planets embedded in the disk . we demonstrate that stellar illumination can heat the gap enough for the planet s orbital eccentricity to instead be excited . we also discuss the `` eccentricity valley '' noted in the known exoplanet population , where low - metallicity stars have a deficit of eccentric planets between @xmath0 and @xmath1 au compared to metal - rich systems @xcite . we show that this feature in the planet distribution may be due to the self - shadowing of the disk by a rim located at the dust sublimation radius @xmath0 au , which is known to exist for several t tauri systems . in the shadowed region between @xmath0 and @xmath1 au lack of gap insolation allows disk interactions to damp eccentricity . outside such shadowed regions stellar illumination can heat the planetary gaps and drive eccentricity growth for giant planets . we suggest that the self - shadowing does not arise at higher metallicity due to the increased optical depth of the gas interior to the dust sublimation radius .