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: recent discoveries of red - skewed iron lines in spectra of neutron star ( ns ) sources serpens x-1 @xcite , 4u 1820 - 30 , gx 349 + 2 ( * ? ? ? * c08 hereafter ) and 4u 1636 - 536 @xcite show that the phenomenon of red - skewed lines is not restricted to black hole ( bh ) sources . in this paper we report on the discovery of the asymmetric iron line in the _ suzaku _ spectrum of the ns source cygnus x-2 ( cyg x-2 ) . therefore , cyg x-2 is the fifth ns source which shows strongly asymmetric iron line profile . this indicates that the red - skewed lines in ns sources may be as common as in bh sources . more generally , asymmetric emission lines appear to be abundant in both types of accretion powered x - ray sources . it is crucial to correctly identify the physical origin of the red - skewness of these lines because they can be potentially used to study the properties of the accretion close to accreting objects as well as to constrain the fundamental characteristics of compact objects . c08 interpreted the @xmath1 iron line profiles in terms of relativistically red - shifted emission due to reflection off the accretion disk very close to a compact object . this scenario is commonly accepted as an explanation of strongly red - skewed iron lines in bh sources @xcite . the main motivation for applying the relativistic line formation scenario to the ns case was the fact that the inner radius of the accretion disk , which was predicted by the relativistic line model was consistent with the interpretation of the highest observed kilohertz quasi - periodic oscillation ( khz qpo ) frequency in this sources as a keplerian frequency at this radius . a red - skewed profile of emission lines can be also produced by repeated electron scattering in a diverging outflow as proposed by @xcite , lt07 hereafter , see also references therein . in the framework of the wind model the fluorescent iron line k@xmath2 is formed in the partly ionized wind as a result of illumination by central source . electron scattering of the iron k@xmath0 photons within the ionized expanding flow leads to a decrease of their energy ( redshift ) . this photon redshift is an intrinsic property of any outflow for which divergence is positive . recently @xcite confirmed using multidimensional monte - carlo simulation that for sufficiently high wind densities , moderate fe k@xmath2 emission lines can be formed and that electron scattering in the flow may cause these lines to develop extended red wings . we examine the red - skewed line profile observed in cyg x-2 both in terms of the relativistic paradigm and in the framework of the wind downscattering . the main obstacle in analyzing the high spectral resolution data with the wind model is that analytical solution is not available for the general formulation of this problem . lt07 used monte carlo ( mc ) simulations to model the line profiles produced in the wind environment . in the presented work we provide a consistent analysis of the line profile with the wind model by introducing the lt07 mc code into xspec astrophysical data analysis package . we find that the wind outflow model is able to reproduce the red - skewed line profile with the fit quality similar to that shown by the relativistic reflection models . therefore , in order to distinguish between these two models one has to consider their consistency in a broader phenomenological context . for example , c08 proposed to look for high frequency quasi - periodic oscillations ( hf qpos ) as an additional evidence of a keplerian disk existence close to a ns , which then would necessitate the presence of the reflection components in energy spectra . in contrast , the presence of the opaque wind in the system would result in smearing the signal coming from the central region . this smearing effect leads to a suppression of the fast variability of x - ray emission . unfortunately , high time resolution observations by _ rossi x - ray timing explorer _ ( _ rxte _ ) , simultaneous with the _ suzaku _ observations analyzed in this paper , are not available . thus , we resort to a study based on the _ rxte _ data set with similar spectral characteristics . our comparative analysis indicates that the source was in low variability state during these observations which is more consistent with the wind / outflow scenario . description of _ suzaku _ and supporting _ rxte _ observations as well as details of our spectral fitting are presented are presented in  2 . we discuss implications of the relativistic and wind outflow red - skewed line formation scenarios in cyg x-2 in  3 . conclusions follow in  4 . cyg  x-2 is a low mass x - ray binary ( see * ? ? ? * for review ) which exhibits a z - shape color - color diagram @xcite . the observations of thermonuclear x - ray bursts @xcite identified the nature of the compact object in cyg x-2 as a neutron star . @xcite used the _ rxte _ burst data to estimate the ns mass to be about 1.4 solar masses and radius to be about 9 km . @xcite reported the simultaneous detection of twin khz peaks at 500 and 860 hz and highest single khz qpo at 1007 hz . cyg x-2 was observed by _ @xcite on may 16 , 2006 for a total exposure of 39 ksec ( observation i d : 401049010 ) . however , during more than a half of the observation the satellite were operating in the medium telemetry mode which led to the telemetry saturation and resulted in data unusable for scientific analysis . during intervals when the high telemetry setting were utilized the foreground illuminated xis detectors ( i.e. xis 0 , 2 and 3 ) operated in 3@xmath33 event editing mode with burst clock and 1/4 window settings . the 3@xmath33 event editing mode is not available for the xis 1 for this observation . we therefore used data from xis 0 , 2 , and 3 collected when high telemetry rate was utilized . we reduce the _ suzaku _ xis and hxd data using the _ xselect _ data analysis tool following the guidelines given by _ suzaku _ data analysis guide . the xis images indicate that strong pile - up in the center of the source point spread function ( psf ) leads to a characteristic `` crater '' . in order to remove piled - up data we extracted xis spectra from the annulus regions with the outer radius set to the maximum allowed by the xis detector field of view ( @xmath4 arcsec ) and the inner radius was manually selected to excise the most piled - up inner core of psf ( @xmath5 arcsec ) . the spectra and corresponding responses were extracted by _ xselect _ extractor and _ xisresp _ , _ xisrmfgen _ tools . spectra and responses for individual detectors were then added using _ mathpha _ and _ addrmf _ ftools . we linearly rebin xis spectral and response data to obtain 1024 spectral channels . hxd / pin spectrum was corrected for non - x - ray and cosmic x - ray background . we fit xis and pin spectra jointly in xspec using 0.7 - 9.0 kev energy range for xis data and 15.0 - 40.0 kev range for pin data fixing the cross - normalization factor between xis and pin spectra at unity . due to large calibration uncertainties we also ignore 1.5 - 2.5 kev range for xis spectrum .    for the continuum spectra we choose the sum of thermal ( bbody ) and comptonized ( comptt ; * ? ? ? * ) components , modified by interstellar photoelectric absorption according to @xcite . when we directly fit the data with this model we observe three distinct narrow features in the residuals ( see figure [ line_res ] , left panel a ) . first , we see a line signature at 6 - 7 kev , which is the primary target of our investigation . we also observe a weak excess around 3.2 kev and a prominent line at 1 kev . the feature at 3.2 kev is probably an instrumental artifact . the line at 1 kev was reported previously from cyg x-2 @xcite and presumably belongs to the source spectrum . both features are well represented by gaussian shape ( see residuals on panel b of figure 1 ) . in addition to these narrow lines we observe a building - up excess towards higher energies in xis spectrum . this indicates a presence of a residual pile - up in the regions close to the excised central part of the psf . to mitigate this effect we utilize the xspec pileup convolution model designed to model pile - up effect in ccd detectors . the xspec implementation of the pile - up model was initially designed to describe this effect in _ chandra _ data . however , the model , developed in @xcite , is valid for _ suzaku _ xis detectors also . _ suzaku _ psf of xis ccds is broader and spreads around larger number of pixels . therefore , pile - up model for the xis spectra would require larger number of regions to be considered for pile - up . to model pile - up in xis detectors we used the following parameter settings for the _ pileup _ xspec model : 2 second time frame , maximum of 10 photons to pile up , unity for grade correction and morphing parameter and 5% of psf to consider for pile - up . the number of detector regions to consider independently for pile - up was allowed to change which led to the best fit value of 38 regions . this is consistent with the number of 3@xmath33 pixel regions close to excised central core of psf . this approach successfully removed the hard excess due to the pile - up effect . finally , we excluded 4.5 - 7.5 kev energy range where the line emission is expected to be significant and fit the spectrum again to obtain the fit quality of @xmath6 . the result is shown in panel d of fig . [ line_res ] . the final best - fit parameters for the continuum model are following : hydrogen column @xmath7 @xmath8 for wabs ; seed photons temperature @xmath9 kev , electron temperature @xmath10 kev , optical depth @xmath11 for comptt and _ bbody _ temperature @xmath12 kev . the energy , sigma and equivalent width ( @xmath13 ) for gaussians used to fit lines are @xmath14 kev , @xmath15 kev , @xmath16 ev and @xmath17 kev , @xmath18 kev , @xmath19 ev . when we include the iron k@xmath0 emission region ( 4.5 - 7.5 kev energy band ) and fit the spectrum with the relativistic and wind line models the parameter ranges of the above continuum parameters change insignificantly ( see table 1 ) . this indicates that continuum model weakly depends on the assumed line model . it is worth noting that the values of the continuum parameters imply that during this _ suzaku _ observation cyg x-2 was in the `` high / soft '' state characterized by high opacity of geometrically thin configuration . panel e in fig . [ line_res ] shows the residuals of the best - fit model with iron region noticed in the data . the apparent line profile is broad and red - skewed . in fact , the fit with the gaussian line model ( model 1 ) produces the worst fit quality among the applied models @xmath20 ( see table 1 ) . we apply physically motivated models to describe this red - skewed line profile . namely , the relativistic models diskline @xcite and laor @xcite and the wind outflow model windline ( lt07 mc code ) . we refer to these models as model 2 , 3 and 4 respectively . the parameters of the laor and diskline models are : the line energy @xmath21 , inner and outer disk radius @xmath22 and @xmath23 , emissivity index @xmath24 and inclination angle @xmath25 . during fitting these relativistic models to the data the value of @xmath23 was fixed at its default value of 1000 @xmath26 and 400 @xmath26 for diskline and the laor respectively . to avoid unphysical results we required the value of the inner radius @xmath22 to be greater than @xmath27 as dictated by the ns compactness limit ( see e.g. * ? ? ? the windline model is not a part of the standard set of xspec models and is implemented as a local model during our modeling . the windline model calculates the line profile by means of mc simulations . the main parameters of the model are the input line energy @xmath21 , the optical depth of the wind @xmath28 , the wind electron temperature @xmath29 and the dimensionless wind speed @xmath30 . the additional fixed parameter of the model is the number of individual photons @xmath31 to be used in the mc simulation . we did nt observe any significant dependence of the @xmath32-statistic behavior on the number of photons used in our line simulations for @xmath31 higher than 3@xmath310@xmath33 . for the presented study we used 5@xmath310@xmath33 photons in the windline mc simulations . we summarize the results of our modeling and fit quality in table 1 . the resulting values of @xmath34 for models 2 , 3 and 4 are 1.31 , 1.30 and 1.32 respectively . in figure [ windline_fit ] we present the unfolded view of spectral fit with model 4 where we used 5@xmath310@xmath35 photons in the mc simulation . we also utilize the data from two _ rxte _ observations made on july 25 , 2006 and september 4 , 2004 ( obs . ids 92039 - 01 - 01 - 01 and 90030 - 01 - 39 - 00 correspondingly ) to investigate the evolution of timing properties of cyg x-2 . we extract the _ rxte_/pca energy spectra from standard2 data modes and we use high resolution modes to calculate power density spectra ( pds ) . we fit _ rxte _ energy spectra in xspec with the model wabs(comptt+bbody+gaussian ) , where gaussian is used to model the iron line . @xmath36 column density was fixed at @xmath37 . we obtain the following best - fit values . for obsid 92039 - 01 - 01 - 01 , comptt : @xmath38 kev ( unconstrained ) , @xmath39 kev , @xmath40 , bbody : @xmath41 kev , gaussian : @xmath42 kev , @xmath43 kev , @xmath44 ev . for obsid 90030 - 01 - 39 - 00 , comptt : @xmath45 kev , @xmath46 kev , @xmath47 , bbody : @xmath48 kev , gaussian : @xmath49 kev , @xmath50 kev , @xmath51 ev . _ spectral parameters for observation i d 90030 - 01 - 39 - 00 match closely those for the _ suzaku _ spectrum_. all sky monitor ( asm ) counts rates in 1.5 - 3/3 - 5/5 - 12 kev energy ranges are 14.06/14.85/14.03 cnts / s and 13.8/14.07/13.2 cnts / s during may 15 , 2006 ( date of _ suzaku _ observation ) and during september 4 , 2004 ( date of 90030 - 01 - 39 - 00 _ rxte _ observation ) respectively , showing almost the same fluxes and source hardness . this also strongly indicates that during the _ suzaku _ observation 401049010 and _ rxte _ observation i d 90030 - 01 - 39 - 00 cyg x-2 was almost in the same state . thus , we can arguably assume that the timing characteristics of the source during _ rxte _ observation 90030 - 01 - 39 - 00 should also be similar to those of the _ suzaku _ observation . in figure [ rxte ] we present _ rxte _ energy spectra along with the power spectra related to _ observations 90030 - 01 - 39 - 00 ( red ) , which are expected to be similar to the analyzed _ suzaku _ observations , and for the observation 92039 - 01 - 01 - 01 ( black ) . the presented evidence for the red - skewed iron line in cyg x-2 along with the detection of the asymmetric lines in serp x-1 , 4u 1820 - 30 , gx 349 + 2 ( * ? ? ? * c08 ) and 4u 1636 - 536 @xcite , increases a number of ns showing this effect to five . this indicates that red - skewed line may be as abundant in ns sources as in bhs ( e.g. * ? ? ? c08 interpreted these observed line appearances as the evidence of relativistic distortion due to the reflection from inner disk located close to the ns surface . the main argument in the association of the line skewness with the inner disk is that the highest observed kilohertz qpo in these sources is consistent with the keplerian frequency at the inner disk radius predicted by the line profile . the diskline model , when applied to the _ suzaku _ ( model 2 ) , leads to the best fit value for the disk inner radius equal to the lower limit set by the model , i.e @xmath52 , which formally translates to the radius of @xmath53 km for the ns mass of 1.4@xmath54 . model 3 , which employs the laor component to represent the iron line , exhibits the same statistical performance as model 2 and yields the radius of 9.5 @xmath26 or 20 km . we should note that the laor model was formulated for the extreme kerr bh case [ see @xcite ] . in this model the spin parameter @xmath55 was assumed to be close to unity . a moderately rotating accreting ns is expected to have a spin parameter less than 0.5 . the ns spin in cyg x-2 is not exactly known but if the difference between khz qpo of @xmath56364 hz is taken as a spin estimate then the space - time background near cyg x-2 should satisfy schwarzshild metric within 10% margin of error . c08 found that their values of @xmath22 were consistent with the inner disk radius predicted by the beat - frequency model @xcite from khz qpo values . maximum khz qpo value observed in cyg x-2 is 1007 khz , which translates to 16.7 km radius if interpreted as a keplerian frequency at the inner disk . this inferred radius is in between the @xmath22 values obtained using diskline and laor models . interpretation of the highest khz qpo peak as a keplerian frequency is not unique . in the framework of the transition layer ( tl ) model the _ lower _ khz qpo peak is classified as a keplerian frequency ( see * ? ? ? * and references therein ) . the tl model successfully describes the behavior of khz qpos in ns sources . the tl model qpo classification leads to the inner disk edge radius of of 22.5 km . beyond this radius the accreting gas enters the tl to adjust its motion to the rotation of the central star . while the value of the disk inner radius given by the diskline model for the iron line can be considered to be in satisfactory agreement with the beat - frequency qpo model , it is harder to reconcile with the tl paradigm . the red - skewed line - khz qpo connection is based on the _ highest observed _ khz qpo values ( c08 ) . however , the duty cycle of khz qpos is low . as noted by c08 , the most compelling evidence for the inner disk origin of the line would come from the simultaneous observation of the khz qpos and the red - skewed iron line . at the time of this writing these two effects were not observed simultaneously from the same source . the reason for absence of the simultaneous detection of khz qpos and broad iron line may be a lack of the correlated observation of high spectral resolution x - ray telescopes ( i.e. _ suzaku_/xmm - newton ) and _ rxte_. the cyg x-2 _ suzaku _ observation analyzed in this paper is a striking example how this lack of simultaneous coverage creates an obstacle in conducting scientific investigation . namely , we have to search _ rxte _ archive for observation with similar spectral characteristics to estimate timing properties ( see the previous section and the discussion below ) while simultaneous _ rxte _ data would allow directly test timing properties more reliably .    a general idea about the fast timing properties can be inferred based on the archival _ rxte _ data by matching the spectral parameters shown by _ rxte _ instruments with those observed in the _ suzaku _ spectrum in question . this approach is justified by the firmly established correlations between spectral characteristics and timing properties in ns lmxbs @xcite and particularly in cyg x-2 @xcite . we searched _ rxte _ archive for pointing observations nearest to may 16 , 2006 . despite the fact that _ rxte _ usually conducts frequent monitoring of cyg x-2 , we found that _ observation was made in the middle of 140 day gap in _ rxte _ monitoring of this source . this , and the fact that cyg x-2 is changing its state on a daily basis , does not permit us to test subsecond timing variability directly with _ rxte_. therefore , we have to resort to search for matching spectral properties and rely on spectral - variability correlations . we found that the nearest _ observation with the spectrum similar to the one shown by _ suzaku _ was performed on september 4 , 2004 ( obsid 90030 - 01 - 39 - 00 , see the previous section ) . to compare this observation with the spectral state characterized by harder spectrum and lower opacity we arbitrarily choose the observation 92039 - 01 - 01 - 01 which is fit by comptt model with parameters @xmath40 and @xmath57 . in figure [ rxte ] we show the energy and power spectra for these observations . it is clear that the variability for frequencies higher than 0.1 hz is strongly suppressed and a `` forest''-type pds appears at frequencies above 0.1 hz . the photon spectrum related to this pds closely matches the _ suzaku _ spectral data . this effect of variability suppression can be readily explained by a smearing of a signal in the opaque wind of optical depth @xmath58 because the suppression factor is related to the opacity in the wind as @xmath59 . a higher gaussian equivalent width of 115 ev related to 90030 - 01 - 39 - 00 _ rxte _ observation with respect to that of 86 ev during 92086 - 01 - 01 - 01 observation is another factor to support the low variability - strong line connection .    the broad band variability continuum observed in pds can be produced by diffusion of disk small perturbations propagating towards central object . these disk perturbations result in a modulation of inner mass accretion rate , which then leads to variability in the x - ray flux @xcite . @xcite developed a theory of the diffusive propagation of perturbations in the disk and presented a model for the power spectrum formation . application of the diffusion theory to the fourier power spectra of cyg x-2 and bh source cyg x-1 led to the conclusion that the fast x - ray variability is produced in two configurations : in a large cold outer accretion disk and a compact geometrically thick configuration ( i.e. transition layer ) . long - term variations of the accretion matter supply at the outer accretion disk edge lead to changes of the source spectral state . specifically , in high - soft state the innermost region becomes very compact and relatively cold as a result of the strong mass accretion . presumably , this is the state which we are dealing with in the case of _ rxte _ observation 90030 - 01 - 39 - 00 and the _ suzaku _ observation . however , in this case strong accretion disk presumably extended close to the central object should produce high frequency perturbations leading to a broad band pds continuum or khz qpos in the range @xmath60 hz . none of these is seen in our data . this may be attributed to wind / outflow attenuation . it could present a problem for the relativistic disk reflection model of the iron emission line production . we note however , that the origin of fast variability and physical processes governing its evolution are not yet fully understood . therefore , the above arguments does not provide a solid proof for the wind line model or a dismissal of the relativistic reflection scenario . futures studies on the asymmetric emission lines in x - ray binaries should address the above points using more substantial observational data . observation of millisecond variability in the same data which would require wind opacities of 2@xmath563 would provide grounds for the rejection of the wind model . such cases have not been observed yet . on the other hand , further evidence of the correlation of `` forest''-type pds with the asymmetric line would rule out relativistic reflection paradigm . cyg x-2 is a z source @xcite and presumably accretes at a rate close to the eddington limit . therefore , strong outflows are expected in this source , naturally leading to the production of the red - wing of the iron line in this wind / outflow configuration . this hypothesis is supported by the fact that z sources have strong radio counterparts . @xcite analysis of radio / x - ray correlation in z and bright atoll sources that radio and comptonized emission in these sources originate in same optically thick plasma with the temperature 2.5 - 3 kev near ns . moreover , strong requirement for the disk illumination to be concentrated very close to its inner edge @xcite in the relativistic line formation scenario led @xcite to consider the production of fluorescent emission within the innermost stable orbit region where matter spirals into the compact object . however , as indicated by the models of @xcite and @xcite , the ionization of such a disk by the intense x - ray radiation might further invalidate some of basic assumptions associated with this interpretation . these arguments indicate that the wind / outflow paradigm may indeed be at work as ( or at least contribute to ) an origin of the red - skewness in the iron lines in compact sources . we present the analysis of the _ suzaku _ spectrum from ns cyg x-2 . we discover an k@xmath0 iron line which shows significant red - wing . this is the fifth ns source so far to show asymmetric line profile . we analyze the line in terms of the relativistic emission from the inner accretion disk and in terms of the wind outflow model . we conclude that both model are acceptable according to the statistical performance . however , the wind model appears to give more adequate explanation which does not require the accretion disk inner edge to advance close to ns surface . we also consider the line production scenarios in the context of the timing properties . we identify _ rxte _ observation which has spectrum very similar to that during _ suzaku _ observation which shows the red - skewed line . these _ rxte _ data indicate that the source fast variability is strongly suppressed , which can be attributed to smearing in a strong wind . this lack of high frequency variability weakens the red - skewed line connection with khz qpo and strengthens its connection to a wind / outflow phenomena . ballantyne , d.r . , ross , r.r . & fabian , a.c . 2001 , mnras , 327 , 10 bhattacharyya , s. , & strohmayer , t.  e.  2007 , , 664 , l103 cackett , e.  m. , et al.2008 , , 674 , 415 davis j.e . 2001 , , 562 , 575 56 fabian , a.  c. , rees , m.  j. , stella , l. , & white , n.  e.  1989 , , 238 , 729 fabian , a. c. 1995 , mnras , 277 , l11 hasinger , g. & van der klis , m. 1989 , a&a , 225 , 79      laming , j.m . & titarchuk , l. 2004 , apj , 615 , l121 laor , a.  1991 , , 376 , 90 lattimer , j.  m. , & prakash , m.  2007 , , 442 , 109 laurent , p. & titarchuk , l. 2007 , , 656 , 1056 lewin , w.  h.  g. , van paradijs , j. , & van den heuvel , e.  p.  j.  1995 , cambridge astrophysics series , cambridge , ma : cambridge university press , |c1995 , edited by lewin , walter h.g . ; van paradijs , jan ; van den heuvel , edward p.j . lyubarskii , yu . , e. 1997 , mnras , 292 , 679    miller , j.  m.  2007 , , 45 , 441 miller , m.  c. , lamb , f.  k. , & psaltis , d.  1998 , , 508 , 791 mitsuda , k. , et al . 2007 , , 59 , 1 morrison , r. , & mccammon , d.  1983 , , 270 , 119 nandra , k. et al . 1999 , apj , 523 , l17 nayakshin , s. , kazanas , d. & kallman , t.r . 2000 , apj , 537 , 833 osherovich , v. , & titarchuk , l.  1999 , , 522 , l113 paizis , a. , et al . 2006 , , 459 , 187 pandel , d. , kaaret , p. , & corbel , s.  2008 , , 688 , 1288 reynolds , c.s . & begelman , m.c . 1997 , apj , 488 , 109 sim , s.  a. , long , k.  s. , miller , l. , & turner , t.  j.  2008 , , 388 , 611 smale , a. p. 1998 , apj , 498 , l141 smale , a. p. , et al . 1993 , apj , 410 , 796 sobolev , v. v. 1957 , soviet astron . , 1 , 678 sunyaev , r. a. , & titarchuk , l. g. 1980 , a&a , 86 , 121 titarchuk , l.  1994 , , 434 , 570 titarchuk , l.  2002 , , 578 , l71 . titarchuk , l. , & shaposhnikov , n.  2002 , , 570 , l25    titarchuk , l. , kuznetsov , s. , & shaposhnikov , n.  2007 , , 667 , 404 titarchuk , l.g . , osherovich , v.a . & kuznetsov , s.i . 1999 , apj , 525 , l129 titarchuk , l. , shaposhnikov , n. & arefiev , v. 2007 , apj , 660 , 556 titarchuk , l. & wood , k. 2002 , apj , 577 , l23 wijnands , r. , et al . 1998 , , 493 , l87    lllll @xmath36 , @xmath61 ......... & 0.23@xmath620.03 & 0.23@xmath620.03 & 0.23@xmath620.03 & 0.23@xmath620.03 + @xmath63 , kev ...... & 0.87@xmath620.01 & 0.87@xmath620.01 & 0.87@xmath620.01 & 0.87@xmath620.01 + @xmath64 , kev ........... & 0.14@xmath620.03 & 0.14@xmath620.04 & 0.14@xmath620.03 & 0.14@xmath620.03 + @xmath65 , kev ......... & 2.28@xmath620.02 & 2.32@xmath620.02 & 2.35@xmath620.02 & 2.28@xmath620.02 + @xmath66 .................... & 17.2@xmath620.4 & 17.1@xmath620.4 & 17.0@xmath620.4 & 17.2@xmath620.4 + @xmath21 , kev .......... & 6.65@xmath620.03 & 6.64@xmath620.14 & 6.49@xmath620.10 & 6.78@xmath620.04 + @xmath67 , kev ............. & 0.22@xmath620.04 & - & - & - + @xmath24 ................... & - & 2.3@xmath620.5 & 3.0 & - + @xmath22 , @xmath26 .......... & - & 6.0 & 9.5@xmath621.3 & - + @xmath68 , deg ............. & - & 30.2@xmath626.8 & 34.0@xmath624.1 & - + @xmath28 ..................... & - & - & - & 1.29@xmath620.35 + @xmath30 , @xmath69 .......... & - & - & - & 1.43@xmath620.14 + @xmath29 , kev ......... & - & - & - & 0.16@xmath620.10 + @xmath70 , ev .......... & 28.5@xmath627.1 & 43.2@xmath6211.0 & 50.6@xmath6213.0 & 37@xmath6215 + @xmath71 ............ & 251.2/187 & 243.3/185 & 241.4/186 & 243.5/185 +
Abstract:

Answer: we report on the _ suzaku _ observation of neutron star low - mass x - ray binary cygnus x-2 which reveals a presence of the iron k@xmath0 emission line . the line profile shows a significant red wing . this discovery increases the number of neutron star sources where red - skewed iron lines were observed and strongly suggests that this phenomenon is common not only in black holes but also in other types of accreting compact objects . we examine the line profile in terms of models which attribute its production to the relativistic effects due to reflection of x - ray radiation from a cold accretion disk and also as a result of the line formation in the extended wind / outflow configuration . both models are able to adequately represent the observed line profile . we consider the results of line modeling in the context of subsecond variability . while we were unable to conclusively disqualify one of the models , we find that the wind paradigm has several advantages over the relativistic disk reflection model .


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: within a given taxonomic group , individual biological species are generally considered to be of equal or near - equal biodiversity value . so , for instance , areas with a greater number of species are more valuable than those with fewer ( myers et al . when wild species are ranked by value , this is usually based on their threat of extinction ( see , e.g. sara 2002 ) . however , as pointed out by cousins ( 1991 ) , species are discovered and identified because they are different from other species , which suggests that they may differ in value . in the context of conservation , avise ( 2005 ) has highlighted five different currencies for valuing species : rarity , distribution , ecology , charisma , and phylogeny . here , we consider the value of a species based on its position in a phylogeny . a phylogeny is the directional , acyclic graph depicting relationships between leaves ( species ) , which we define formally in the next section . a phylogeny generally has a root ( which assigns direction ) and edge weights that can represent unique feature diversity ( e.g. as measured by evolutionary time or genetic distance ) . species can be defined by the features they possess , and one measure of their worth is their expected contribution of unique features . in this way , we can use a phylogeny to assign a measure of evolutionary value to a species based on its expected contribution of unique features . because of the highly imbalanced shape of the tree of life , some species in a phylogeny will have far fewer close relatives than others in that phylogeny ( mooers and heard 1997 ) , and these more distantly - related species will be expected to contribute more unique features ( faith 1992 ) . phylogenetic measures of conservation value have a long pedigree ( see , e.g. alschul and lipman 1990 ; may 1990 ) and have begun to be explored in some detail ( haake et al . 2005 ; hartmann and steel 2007 ; pavoine et al . , 2005_b _ ; redding and mooers 2006 ) . so , for example , pavoine and colleagues presented one new measure of originality , a set of sampling weights such that the expected pairwise distance on the tree is maximized . haake and colleagues extended the ` shapley value ' ( shapley 1953 ) from co - operative game theory to the conservation setting to calculate the average distance of a focal species to all possible subsets of taxa . for both measures , more original species are those expected to contribute more to the resulting sets . yet another measure that uniquely apportions the tree to its tips ( isaac et al . 2007 ) and which is the focus of a new international conservation initiative ( the edge initiative , zoological society of london ) scales almost perfectly with the shapley value ( unpublished results ) . one question with these measures concerns the sets that individual species are asked to complement . for instance , given known extinction probabilities for species , some future sets are much more likely than others and so some species will be more valuable because their close relatives are less likely to be included in future sets . here we formalize this idea to extend the shapley value of a species to include pre - assigned extinction probabilities . we then compare our measure with the original shapley value using the prosimian primates as a test case . * definition * let @xmath0 be a rooted or unrooted phylogenetic tree with leaf set @xmath1 , together with an assignment of positive lengths to the edges ( branches ) of @xmath0 . we let @xmath2 denote the length of edge @xmath3 , and let @xmath4 denote the set of edges of @xmath0 . for a subset @xmath5 of @xmath1 , let @xmath6 denote the _ phylogenetic diversity _ of @xmath5 defined as follows . if @xmath0 is unrooted then @xmath6 is the sum of the lengths of the edges ( branches ) of @xmath0 in the minimal subtree that connects @xmath5 . if @xmath0 is rooted , then @xmath6 is the sum of the lengths of the edges of @xmath0 in the minimal subtree of @xmath7 that connects @xmath5 and the root of the tree . figure 1 illustrates these concepts , and includes values at the tips that we will use in the next section . note that although the branch lengths in this example are clock - like , this assumption is not required in any of the results we describe . units for the terminal edges incident with c , d , e , and @xmath8 unit for the other five edges ) . each tip @xmath9 has an associated extinction probability @xmath10 . ( b ) for a subset @xmath11 of taxa that are extant at some future time , the phylogenetic diversity score @xmath6 is the sum of the lengths of the edges indicated in bold . the dashed edges lead to extinct taxa . ] for a leaf @xmath12 , and a subset @xmath13 let @xmath14 the quantity @xmath15 measures how much phylogenetic diversity @xmath16 contributes to the tree that one obtains from @xmath0 once species not in @xmath5 have been pruned out ( for example if they go extinct ) . alternatively , @xmath15 is the marginal increase in phylogenetic diversity of @xmath5 if @xmath16 is added .    now , suppose that each species has an associated extinction probability @xmath17 ( which may vary from species to species )  for example , this may be the probability that the species is extinct in ( say ) 100 years from now ( either globally , or in some specified community ) . we will denote this @xmath17 value for species @xmath9 by @xmath18 . in this paper we consider the simplest model which assumes that the extinction of each species in @xmath1 comprise independent events . given @xmath12 , let @xmath19 denote the random subset of species in @xmath20 which survive ( i.e. do not go extinct ) .    by the independence assumption we have : @xmath21 = \prod_{j \in s } ( 1-\epsilon_j ) \times \prod_{j \in x-\{i\}-s } \epsilon_j.\ ] ]    for @xmath12 , let @xmath22 denote the expected value of @xmath23 . that is , @xmath24 = \sum_{s \subseteq x-\{i\ } } { { \mathbb p}}[{\mathcal s}_i = s ] \delta_{pd}(s , i).\ ] ] we call @xmath22 the _ heightened evolutionary distinctiveness _ of species @xmath16 , and the function @xmath25 the _ heightened evolutionary distinctiveness _ ( he d ) index for @xmath0 . notice that if all the species in @xmath20 were guaranteed to survive , then @xmath22 would be just the length of the pendant edge incident with leaf @xmath16 , however random extinctions mean that @xmath22 will tend to be increased ( ` heightened ' ) over this pendant edge length . a related but different index , based on the shapley value in co - operative game theory , has recently been described by haake et al . this index , denoted here as @xmath26 can be defined ( for unrooted trees ) as follows : for @xmath12 , @xmath27 this index has certain appealing properties . in particular , @xmath28 , and there is a simple formula for quickly computing @xmath29 . the index @xmath26 also has a stochastic interpretation , but this is not based on extinction or survival of species , rather on the expected contribution to @xmath30 of each species under all possible orderings of the total set of species ( for details see haake et al . the index @xmath31 allocates existing @xmath30 ` fairly ' amongst the species , whereas @xmath32 quantifies the expected contribution of each species to future @xmath30 . computing the he d index directly via ( [ defeq ] ) could be problematic as it requires summation over all the subsets of @xmath20 and this grows exponentially with @xmath33 . however we now show that the index can be readily and quickly computed , both for rooted and unrooted trees . this polynomial - time algorithm for computing @xmath32 thus complements ( but is quite different to ) the polynomial - time algorithm described by haake et al . ( 2005 ) for computing @xmath26 .      for a rooted phylogenetic @xmath1tree @xmath0 let @xmath34 denote the set of species in @xmath1 that are descended from @xmath3 ( i.e. the clade that results from deleting @xmath3 from @xmath0 ) . for @xmath12 , let @xmath35 ( @xmath36 ) denote the edges ( branches ) on the path from @xmath16 to the root of @xmath0 , listed in the order they are visited by that path . recall that @xmath2 denotes the length of edge @xmath3 . the proof of the following theorem is given in the appendix . [ th1 ] @xmath37    note that in this ( and the next ) theorem we adopt the convention @xmath38 , which is relevant for the first term ( @xmath39 ) in the sum as @xmath40 is empty . thus the first term in the summation expression for @xmath22 given by theorem  [ th1 ] is simply @xmath41 , the length of the pendant edge of @xmath0 incident with species @xmath16 . we can apply the he d index to the members of the rooted tree depicted in fig . 1 . for example , to compute @xmath42 by using theorem  [ th1 ] we have @xmath43 . by inspection , we can see that the most valuable species will be @xmath44 , since it shares an edge with only one other species above the root , and that this species ( @xmath45 ) has a high @xmath17 . at the other extreme , @xmath46 shares its path to the root with two other species , and one of them ( @xmath47 ) has a low @xmath17 . it should therefore receive a low he d value . the computed values are @xmath48 , @xmath49 , @xmath50 , @xmath51 , and @xmath52 . using the shapley index ( haake et al . 2005 ) , @xmath44 and @xmath45 are ranked first ( with value = 2.63 ) , followed by @xmath47 ( 2.33 ) and then @xmath46 and @xmath53 ( 1.75 ) . pavoine s qe metric ( pavoine et al . 2005 ) returns the same ranking as does the shapley . a portal for computing he d is available at http://www.disconti.nu/-phylo/emd.dpf      we now provide a similar formula for efficiently computing the he d index for unrooted trees . given a leaf @xmath16 of @xmath0 and an edge @xmath3 of @xmath0 , @xmath3 induces a split of @xmath1 into two disjoint subsets , and one of these subsets , which we denote as @xmath54 , contains @xmath16 . the proof of the following theorem is given in the appendix . [ th2 ] @xmath55    notice that the rooted he d index is just a special case of the unrooted he d index ( indeed theorem  [ th1 ] can be deduced from theorem  [ th2 ] ) . to see this , given a rooted tree @xmath0 attach a new leaf @xmath56 to the root via a new edge to obtain an unrooted tree , and assign the new edge length @xmath57 . let @xmath58 . then it is easily seen that the he d index for @xmath0 is just the he d index for the derived unrooted tree .    using theorem  [ th1 ] it can be shown that if @xmath0 is a rooted phylogenetic tree then the condition : @xmath59 holds for all selections of positive branch lengths and @xmath60 s if and only if @xmath0 is a ` star tree ' ( that is , every leaf is adjacent to the root ) . moreover theorem  [ th2 ] shows that there is no unrooted phylogenetic tree @xmath0 for which ( [ pareto ] ) holds for all positive branch lengths and @xmath60 values ( of course ( [ pareto ] ) may hold on phylogenetic trees  either rooted or unrooted  if the branch lengths and @xmath60 values take certain values ) . this contrasts with the index @xmath26 which satisfies @xmath61 on all unrooted phylogenetic trees and choices of branch lengths , a property that is referred to as the pareto efficiency axiom by haake et al . ( 2005 ) . in the setting of this paper we should not be surprised that ( [ pareto ] ) holds for @xmath32 only in very special cases since we are not trying to divide out existing @xmath30 amongst present taxa ( one motivation behind @xmath26 ) but rather quantify the expected contribution each species makes to future @xmath30 . we compared the he d index with the shapley ( haake et al . 2005 ) values for the prosimians ( mammalia : primata ) , a group of approximately 50 species with a broad range of extinction probabilities . this group includes the aye - aye , the lemurs , the lorises and galagos . we made use of a recent dated supertree of the order primates ( vos and mooers 2004 ; vos 2006 ) , see fig . 2 , and red list risk designations from the iucn ( www.iucnredlist.org , accessed february 2006 ) . following isaac et al . ( 2007 ) and redding and mooers ( 2006 ) , we first converted the five categories of risk ( cr , en , vu , nt , and lc ) to probabilities of extinction . under the iucn criteria , the species in the vu category are given a @xmath17=0.1 over the next 100 years . we gave the lowest and highest threat categories very conservative probabilities of extinction over the next 100 years of 0.001 and 0.9 respectively , leaving @xmath62 for en , and @xmath63 for nt : this scale is very similar to that calculated from real population viability analyses for birds ( redding and mooers 2006 ) . we are primarily interested in how the ranking of species changes using different approaches .    : critically endangered , @xmath64 ; @xmath65 : endangered , @xmath62 ; @xmath66 : vulnerable , @xmath67 ; @xmath68 : near threatened , @xmath63 ; @xmath69 : least concern , @xmath70 . edge lengths are on an arbitrary scale that represents time since divergence . ] the bivariate correlation between the metrics is high ( 0.94 ) . both measures chose the aye - aye ( @xmath71 @xmath72 ) as the most important species , followed by @xmath73 @xmath74 . interestingly , the three most highly ranked species under current conservation policy ( the critically endangered lemurs @xmath75 @xmath76 , @xmath77 @xmath78 , @xmath79 @xmath80 ) are nested well up in the tree ( figure 2 ) such that none of them were chosen in the top ten for either sv or he d . if we compare the rest of the rankings for these two metrics , the largest single difference is for the two @xmath81 species : they rank twelfth under sv ( being relatively isolated on the tree ) , but only twenty - sixth under he d : because neither is severely threatened , the chances are good that their common path will persist . both measures are very heavily influenced by the pendant edge ( @xmath82 ) length of the focal species ( with correlations @xmath82 vs. sv=0.94 , and @xmath82 vs. hed=0.98 ) . @xmath83 is always part of the marginal increase to @xmath5 , while interior nodes are most likely represented with high probability , especially for larger and more balanced trees . @xmath82 is , however , a poor predictor of he d for @xmath75 , for @xmath79 @xmath84 , and for @xmath85 @xmath86 ( figure 2 ) . the first two groups contain the three most endangered species , increasing the value of close relatives . @xmath85 is an isolated genus , and @xmath87 s sister species is listed as vulnerable ( @xmath67 ) . likewise for @xmath88 - although it has close relatives and so a short pe , these relatives are at high risk of extinction , which increases its value ; this is what we saw with species @xmath44 in figure 1 . the effect of close relatives risk status on one s own value is precisely the strength of the he d approach . however , the fact that the extinction risk of other species affects a focal species , but its own risk does not is somewhat counter - intuitive . we address this by showing how it is possible to write the he d index as the sum of two terms each of which takes into account the extinction risk of the focal species . to describe this further , let @xmath89 be the random variable which takes the value @xmath16 if the focal species @xmath16 survives ( at the future time under consideration ) and which otherwise takes the value of the emptyset ( i.e. @xmath90 ) if @xmath16 goes extinct . let @xmath91,\ ] ] where , as before @xmath19 is the random subset of species in @xmath20 that survive . in words , @xmath92 is the increase in the expected pd score if we condition on the event that species @xmath16 survives . similarly , let @xmath93.\ ] ] in words , @xmath94 is the decrease in the expected pd score if we condition on the event that species @xmath16 becomes extinct . the following result describes how to compute these two indices easily from the he d index , and verifies that they add together to give the he d index ( its proof is given in the appendix ) . [ th3 ]    * @xmath95 , * @xmath96 , * @xmath97 .    the approach of assigning a value to a species which is a function of its phylogenetic distinctiveness and its extinction probability has been referred to as ` expected loss ' by redding and mooers ( 2006 ) and , more evocatively , an ` edge ' score ( evolutionarily distinct and globally endangered ) by isaac et al . ( 2007 ) .    in the same spirit we will call @xmath98 and @xmath99 ( which extend our he d index @xmath32 ) hedge ( _ heightened evolutionary distinctiveness and globally endangered _ ) scores . the hedge score @xmath92 is more relevant when evaluating actions that might save species , whereas the hedge score @xmath94 is appropriate when evaluating actions that might cause the extinction of species ( such as building a dam ) .    one potential advantage of he d and hedge over previous scores is their flexibility in designing conservation scenarios . so for instance , we can choose iucn - ranked species for which conservation is cheap and/or already partially successful , set their @xmath17 to 0 , and see how rankings of other species change . alternatively , we might want to increase the @xmath17 to 1.0 for certain species to see how others are affected . most generally , he d and hedge could be incorporated in an assessment of species value that included many factors besides risk and future contribution , e.g. the ecological , distributional and aesthetic values enumerated by avise ( 2005 ) , and the costs of recovery and probability of its success . _ proof of theorem  [ th1 ] _    first observe that the only edge lengths that contribute to @xmath15 are those from the set @xmath100 . consequently , for the random set @xmath19 of surviving species of @xmath20 we have @xmath101 where @xmath102 is the @xmath103 indicator random variable that takes the value @xmath8 precisely if @xmath104 is not an edge of the subtree of @xmath0 connecting the taxa in @xmath19 and the root of @xmath0 ; since this is the only situation for which @xmath104 lies in the subtree of @xmath0 connecting @xmath105 but not in the subtree of @xmath0 connecting @xmath19 . thus , by linearity of expectation , @xmath106 = { { \mathbb e}}[\sum_{r\geq 1 } l(e_r ) \cdot i_r({\mathcal s}_i ) ] = \sum_{r\geq 1 } l(e_r ) \cdot { { \mathbb e}}[i_r({\mathcal s}_i)],\ ] ] and since @xmath102 takes the values @xmath57 and @xmath8 , @xmath107 = { { \mathbb p}}[i_r({\mathcal s}_i)=1]$ ] . thus , @xmath108.\ ] ] now , the event ` @xmath109 ' occurs precisely if none of the elements in @xmath40 survive , and this latter event has probability @xmath110 . substituting this into the previous equation establishes the theorem . @xmath111    _ proof of theorem  [ th2 ] _    for @xmath12 and the random subset @xmath112 , we have @xmath113 where @xmath114 is the @xmath103 indicator random variable taking the value @xmath8 precisely if @xmath19 consists of no elements of @xmath115 and at least one element of @xmath116 . thus , @xmath117\ ] ] and by the independence assumption @xmath118 = { { \mathbb p}}[{\mathcal s}_i \cap ( c_i(e)-\{i\ } ) = \emptyset]\cdot   { { \mathbb p}}[{\mathcal s}_i \cap ( x- c_i(e ) ) \neq   \emptyset]\ ] ] and so @xmath118 = ( \prod_{j \in c_i(e)-\{i\}}\epsilon_j ) \cdot ( 1-\prod_{j \in x- c_i(e)}\epsilon_j),\ ] ] as claimed . @xmath111    _ proof of theorem  [ th3 ] _    by definition @xmath119 $ ] , and so @xmath120 = \psi_i + { { \mathbb e}}[pd({\mathcal s}_i)]$ ] . now we can write the unconditional expectation @xmath121 $ ] as the weighted sum of conditional expectations @xmath122 { { \mathbb p}}(i = \{i\ } ) + { { \mathbb e}}[pd({\mathcal s}_i \cup i)|i = \emptyset ] { { \mathbb p}}(i = \emptyset)$ ] and so @xmath123 = ( 1-\epsilon_i){{\mathbb e}}[pd({\mathcal s}_i\cup \{i\ } ) ] + \epsilon_i { { \mathbb e}}[pd({\mathcal s}_i)].\ ] ]    parts ( i ) and ( ii ) now follow by applying these equations ( and the linearity of expectation ) to the definitions of @xmath92 and @xmath94 . part ( iii ) follows directly from parts ( i ) and ( ii ) . @xmath111 this research was supported by a grant from the marsden fund , new zealand to ms and aom . aom received additional support from nserc canada and the institute for advanced study , berlin . we thank dave redding , klaas hartmann , and andy purvis for discussion . 99 altschul , s. f. , and d. j. lipman . equal animals . _ nature _ 348:493494 . avise j.c . phylogenetic units and currencies above and below the species level . pages 77 - 100 _ in _ purvis a , gittleman jl , and brooks t , eds . phylogeny and conservation cambridge university press , cambridge , uk . cousins , s.h . species diversity measurement : choosing the right index . _ trends in ecology and evolution _ 6:190192 . faith , d.p . 1992 . conservation evaluation and phylogenetic diversity . _ biological conservation _ 61:1 - 10 . haake , c - j , kashiwada , a. , and f.e . su . 2005 . the shapley value of phylogenetic trees . imw working paper 363 ( 363 ) . arxiv ( q-bio.qm/0506034 ) hartmann , k. and m. steel . phylogenetic diversity : from combinatorics to ecology . _ in _ o. gascuel and m. steel eds . reconstruction evolution : new mathematical and computational advances . oxford university press . isaac , n.j.b . , turvey , s.t . , collen , b. , waterman , c. , and j.e.m . mammals on the edge : conservation priorities based on threat and phylogeny . _ plos one _ 2(3 ) : e296 . may , r.m . taxonomy as destiny . _ nature _ 347:129130 . mooers , a.o . and s.b . 1997 . inferring process from phylogenetic tree shape . _ quarterly review of biology _ 72:31 - 54 . myers , n.,mittermeier , r. a. , mittermeier , c. g. , da fonseca , g. a. b. , and j. kent . 2000 . biodiversity hotspots for conservation priorities . _ nature _ 403:853858 . pavoine , s. , ollier , s. , and d. pontier . measuring diversity from dissimilarities with rao s quadratic entropy : are any dissimilarities suitable ? _ theoretical population biology _ 67:231239 . pavoine , s. , s. ollier , and a.b . dufour . 2005 . is the originality of a species measurable ? _ ecology letters _ 8:579 - 586 . redding , d.w . and a.o . mooers . 2006 . incorporating evolutionary measures into conservation prioritization . _ conservation biology _ 20:1670 - 1678 . bill c-5 , an act respecting the protection of wildlife species at risk in canada . government of canada , ottawa . shapley , l.s . 1953 . a value for n - person games . _ annals of mathematical studies _ 28:307 - 317 . vos , r.a . , and a.o . mooers . 2004 . reconstructing divergence times for supertrees : a molecular approach . pages 281 - 299 _ in _ bininda - emonds , ed . phylogenetic supertrees : combining information to reveal the tree of life . kluwer academic press , 765 . vos , r. 2006 . inferring large phylogenies : the big tree problem . thesis , simon fraser university .
Abstract:

Answer: if predictions for species extinctions hold , then the ` tree of life ' today may be quite different to that in ( say ) 100 years . we describe a technique to quantify how much each species is likely to contribute to future biodiversity , as measured by its expected contribution to phylogenetic diversity . our approach considers all possible scenarios for the set of species that will be extant at some future time , and weights them according to their likelihood under an independent ( but not identical ) distribution on species extinctions . although the number of extinction scenarios can typically be very large , we show that there is a simple algorithm that will quickly compute this index . the method is implemented and applied to the prosimian primates as a test case , and the associated species ranking is compared to a related measure ( the ` shapley index ' ) . we describe indices for rooted and unrooted trees , and a modification that also includes the focal taxon s probability of extinction , making it directly comparable to some new conservation metrics .    _ keywords : _ phylogenetic diversity , extinction , biodiversity conservation , shapley index    _ short title : _ taxon - specific indices of expected future biodiversity


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: even 60 years after its discovery , the casimir effect continues to fascinate  @xcite . originally calculated as an attractive force between two neutral metallic plates due to the zero - point energy of the electromagnetic field , the casimir effect nowadays plays a role in fields as diverse as gravitation , cosmology , atomic and molecular physics , and quantum field theory  @xcite . from a condensed matter perspective , research has been mainly preoccupied with the inclusion of several contributions that are present in a real experimental setup , e.g. , boundary roughness , shape and geometry of the mirrors , non - zero temperature , and dielectric properties of the boundaries ( a thorough review can be found in ref . @xcite ) . with the theory in place , several recent experiments have verified the existence of the casimir force  @xcite ) .    whereas most publications deal with optically isotropic mirrors , the casimir effect between anisotropic mirrors is a rather recent topic . the van der waals interaction between uniaxial anisotropic media of infinite thickness has been considered in ref . @xcite and formed the starting point for predictions of a torque between two uniaxial birefringent plates  @xcite . the zero - point energy depends on the relative orientation of the plates which gives rise to a torque that seeks to align their principal axes . similar effects are present for plates in which the optical anisotropy is induced by corrugations  @xcite . furthermore , in a search for repulsive instead of attractive casimir forces , and stimulated by the advent of so - called magneto - dielectric metamaterials ( which have a non - negligible magnetic response on top of the electric one ) , the effects of anisotropy of the magnetic permeability on the casimir effect were studied in refs . the casimir force between ferromagnetic mirrors , which of course can not be described by a magnetic permeability , has been treated by us in a recent publication  @xcite . in this case , the boundary conditions for the zero - point electromagnetic field are influenced by the magneto - optical kerr effect , and we have found that this gives rise to a magnetic contribution to the casimir force .    in the present paper , the force between a single ferromagnetic mirror opposite of a non - magnetic one is considered . interestingly , we find that the casimir effect will give rise to a _ magnetic _ anisotropy in the ferromagnet . we will first study a setup where the non - magnetic mirror is optically isotropic , and observe an out - of - plane magnetic anisotropy in the ferromagnetic layer . for this setup , the anisotropy energy arising from the casimir effect is orders of magnitude smaller than the surface anisotropy and magnetocrystalline anisotropy . therefore we will propose another scheme where the non - magnetic mirror is made of a uniaxially birefringent material ( with optical axis in the plane ) . the casimir effect then induces an in - plane magnetic anisotropy that can be detected more easily . in both cases , the magnetic orientation with the lowest energy results from a competition between several contributions to the magneto - optical kerr effect , and will depend on the distance between the plates . the casimir energy ( per unit area ) between two plates @xmath0 and @xmath1 can be expressed in terms of the reflection coefficients of the plates as  @xcite    @xmath2 , \label{eq : casenergy}\ ] ]    where @xmath3 is the electromagnetic field frequency , @xmath4 is the component of the wavevector perpendicular to the mirrors , and @xmath5 is the angle between the incidence plane and the @xmath6-axis . furthermore , @xmath7 are @xmath8 matrices containing the reflection coefficients of the plates ( evaluated at imaginary wavevector and frequency ) : @xmath9 the index @xmath10 ( resp . @xmath11 ) corresponds to a polarization with the electric field perpendicular ( resp . parallel ) to the incidence plane . we adopt the convention that the @xmath11 axis remains unchanged upon reflection .    in this paper , we will consider one of the plates to be ferromagnetic . the reflection coefficients are then influenced by the magneto - optical kerr effect , and depend on the magneto - optical constant @xmath12 . since @xmath12 is small , we will only keep terms up to first order in @xmath12 , in which case the reflection matrix can be written as @xmath13 , with @xmath14\ ] ] independent of @xmath12 , while @xmath15\ ] ] depends on the magnetization . the magnetization orientation has been written out explicitly in terms of the spherical angles @xmath16 and @xmath17 ( the @xmath18-axis points into the ferromagnetic plate ) . the terms @xmath19 and @xmath20 are the fresnel coefficients for a normal metallic mirror , while the terms @xmath21 , @xmath22 , and @xmath23 arise from the polar , longitudinal and transverse kerr effect respectively and are due to the magnetization . all coefficients can be expressed in terms of the diagonal ( @xmath24 ) and off - diagonal ( @xmath25 ) components of the dielectric tensor  @xcite :    [ eq : reflcoefffm ] @xmath26\left[\epsilon_{xx } ( { \rm i}\omega)k_{\perp}c+\xi\right]},\\ r_{sp}^{\parallel}({\rm i}\omega,{\rm i}k_{\perp } ) & = & \frac{-k_{\perp}c \sqrt{\omega^{2}-(k_{\perp}c)^{2 } } \ , \omega\ , \epsilon_{xy}({\rm i}\omega ) } { \left(k_{\perp}c+\xi\right)\left(\epsilon_{xx } ( { \rm i}\omega)k_{\perp}c+\xi\right)\xi } , \\ \delta r_{pp}({\rm i}\omega,{\rm i}k_{\perp } ) & = & \frac{2\sqrt{\omega^{2}-(k_{\perp}c)^{2 } } \ , \epsilon_{xy}({\rm i}\omega ) \ , ( k_{\perp}c ) } { \left(\epsilon_{xx } ( i\omega)k_{\perp}c+\xi\right)^{2}},\end{aligned}\ ] ]    where we have defined : @xmath27 . using the reflection coefficients above , we will now consider the casimir effect between a non - magnetic isotropic mirror @xmath0 and a ferromagnetic mirror @xmath1 . for this particular case , mirror @xmath0 is described by a reflection matrix @xmath28 , and mirror @xmath1 by @xmath29 as defined in eqs . ( [ eq : reflmatrixr0 ] ) ,  ( [ eq : reflmatrixdr ] ) and  ( [ eq : reflcoefffm ] ) . making an expansion in @xmath30 of the logarithm in eq . ( [ eq : casenergy ] ) and keeping only terms in the lowest non - vanishing order , one obtains the following expression for the casimir energy : @xmath31 in this expression , @xmath32 is independent on the magnetization orientation and describes the non - magnetic contribution to the casimir energy , while    [ eq : casenergycontributionsexpressions ] @xmath33    are terms arising respectively from the polar , longitudinal and transverse contributions to the magneto - optical kerr effect . the casimir energy in eq .  ( [ eq : casenergycontributions ] ) is of course independent of the in - plane magnetization , because the considered system is invariant with respect to rotations around the normal to the plates . however , the energy depends on the angle @xmath16 between the magnetization and the normal to the plates , i.e. , the casimir effect gives rise to an out - of - plane magnetic anisotropy . it is important to note that the contribution coming from the polar kerr effect , @xmath34 is proportional to @xmath35 and thus tends to align the magnetization parallel to the plates ( @xmath36 or @xmath37 ) , while the longitudinal and transverse contributions @xmath38 are proportional to @xmath39 and have their energy minimum when the magnetization is perpendicular to the plates ( @xmath40 ) . the physical parameters ( dielectric constants ) of the plates , and as we will see , the interplate distance , will decide which of these competing terms is dominant .    we have numerically calculated the different contributions for a setup where the non - magnetic plate is made of gold , and the magnetic plate is made of iron . the reflection coefficients in eq . ( [ eq : reflcoefffm ] ) will depend on the dielectric tensor of these metals . in order to take the interband transitions properly into account , we have calculated the dielectric tensor from the available optical data , instead of reverting to a simple plasma or drude model . the diagonal part of the dielectric response function @xmath41 of @xmath42 and @xmath43 , evaluated at real frequencies , is tabulated in ref . a significant amount of data points is available in the range @xmath44 for @xmath42 and @xmath45 for @xmath43 . at lower frequencies , where interband transitions do not play a significant role anymore , we have extrapolated the data with a drude model using a plasma frequency @xmath46 and a relaxation time @xmath47 for @xmath42 , while @xmath48 and @xmath49 for @xmath43  @xcite . optical data for the off - diagonal element of the dielectric tensor @xmath50 of @xmath43 is rather scarce . we have used data points between @xmath51 and @xmath52 found in ref . @xcite , and assumed @xmath53 out of this range . the dielectric constants @xmath54 and @xmath55 at imaginary frequencies , which are needed for evaluating the reflection coefficients in eq . ( [ eq : reflcoefffm ] ) , were then obtained by using the kramers - krnig relations . details of the whole procedure , together with plots of the resulting dielectric tensors can be found in ref . @xcite and our earlier work  @xcite . finally , the casimir energy per unit area is calculated by numerical integration of eqs . ( [ eq : casenergycontributionsexpressions ] ) , for interplate distances @xmath56 ranging from @xmath57 to @xmath58 . results for the casimir energy amplitudes @xmath34 , @xmath59 , and @xmath60 are given in fig . [ fig1](a ) . the corresponding force contributions are obtained by deriving the energy expressions with respect to the interplate distance @xmath56 : @xmath61 and are shown in fig . [ fig1](b ) .    for distances larger than approximately @xmath62 , corresponding to the experimentally relevant regime , the contribution @xmath34 arising from the polar kerr effect is clearly dominant . looking back at eq . ( [ eq : casenergycontributions ] ) , this means that the lowest energy state at these distances is one where the magnetization is perpendicular to the plates . however , at smaller distances there is a crossover to a regime where the transverse contribution @xmath60 dominates , and consequently the preferential magnetization direction will interestingly change from perpendicular towards parallel to the plates upon decreasing the interplate distance . looking at the distance dependence of the force ( per unit area ) in fig . [ fig1](b ) , we have found that for long distances the polar contribution varies like @xmath63 , compared to @xmath64 for the other contributions . for very short distances , the transverse force varies like @xmath65 , while the polar and longitudinal terms show ( approximately ) a @xmath66 dependence .    in order to estimate the casimir force in a typical experiment , we consider two disks of radius @xmath67 , at a distance @xmath68 from each other . for this setup , we obtain a force amplitude @xmath69 on the order of a few fn . the anisotropy could thus possibly be detected with state - of - the - art mechanical ( force ) measurement schemes ( @xmath70 resolution is reached in magnetic resonance force microscopy  @xcite ) , but presumably not via some magnetic measurement because the anisotropy contribution due to the casimir effect is considerably smaller ( by several orders of magnitude ) than the shape anisotropy and the magnetocrystalline surface anisotropy of the @xmath43 plate . we have therefore looked at a different configuration where the latter anisotropies are essentially zero . in particular , we considered a setup where the magnetic mirror is now placed opposite of a non - magnetic mirror with uniaxial ( optical ) anisotropy . the latter can , e.g. , be made from a uniaxial birefringent material with the optical axis in the plane . for this kind of system , we have found that the casimir effect will create an in - plane magnetic anisotropy , which could be much easier to detect . our results will be presented in the next section . the reflection coefficients of a mirror with uniaxial optical anisotropy can be expressed in terms of the ordinary ( @xmath71 ) and extraordinary ( @xmath72 ) dielectric constants of the anisotropic material . the corresponding equations can be found , e.g. , in refs . we will consider a mirror with the optical axis in the plane , for which one obtains ( at imaginary frequency and wavevector )  @xcite :    @xmath73   \right ] \nonumber \\ d ' r_{pp}^\text{uni } & = & -\left(k_\perp c + \xi_o\right ) \left(\frac{\xi_o}{\epsilon_o k_\perp c } - 1 \right ) \xi_\alpha^2 \\ & \ & \ \ \ + \left(\xi_e - \xi_o\right ) \left [ \xi_\alpha^2 - \frac{\xi_o}{k_\perp c}\left [ \omega^2 \beta^2 - ( \alpha k_\perp c)^2 \right ]   \right ] \nonumber \\ r_{sp}^\text{uni } & = & -r_{ps}^\text{uni } = \frac{2 \alpha \beta \xi_o \left(\xi_o - \xi_e \right ) \omega}{d'},\end{aligned}\ ] ]    with    @xmath74 , \nonumber \\ \xi_o & = & \sqrt{\omega^2 \left ( \epsilon_o - 1\right ) + \left ( k_\perp c\right)^2 } , \nonumber \\ \xi_e & = & \sqrt{\xi_o^2 + \frac { \epsilon_e - \epsilon_o}{\epsilon_o } \xi_\alpha^2 } , \nonumber \\ \xi_\alpha & = & \sqrt{\left(\epsilon_o - \alpha^2\right ) \omega^2 + ( \alpha k_\perp c)^2 } , \nonumber\end{aligned}\ ] ]    and where @xmath75 and @xmath76 , with @xmath77 the angle between the optical axis and the @xmath6-axis . the dielectric functions @xmath71 and @xmath72 are evaluated at imaginary frequencies . please note that in comparison with the expressions in ref . @xcite , both @xmath20 and @xmath78 have a minus sign . this results from a different convention in ref . @xcite , where the @xmath11 axis reverses sign upon reflection , whereas we keep this axis unchanged  @xcite . furthermore , the definition of @xmath79 and @xmath80 is different because we assume the incidence plane to be the @xmath81 plane as in ref . @xcite , whereas this is the @xmath82 plane in ref . @xcite . it is clear that the angular dependence of these reflection coefficients is much more complicated than what we had for a magnetic mirror in eq . ( [ eq : reflmatrixdr ] ) . therefore it is more difficult to make analytical progress , and we have opted for a direct numerical evaluation of the casimir energy in eq . ( [ eq : casenergy ] ) . we have considered a setup where one ferromagnetic ( @xmath43 ) plate with magnetization in the mirror plane is opposed to a non - magnetic uniaxially anisotropic mirror . the dielectric tensor of the @xmath43 plate is calculated directly from optical data , as has been described in the previous section . the uniaxial plate is assumed to be made of one of the birefringent materials quartz ( @xmath83 ) , calcite ( @xmath84 ) or barium titanate ( @xmath85 ) . since optical data for these materials is not available in a large enough frequency range , we have used a two - oscillator model to describe their dielectric function . although this model is quite simple , it gives a rather good description of the dielectric tensor in inorganic materials  @xcite , and has been used already to estimate the casimir torque between two birefringent plates  @xcite . within this approach , both the ordinary and extraordinary dielectric function of the uniaxial mirror will be modeled by : @xmath86 where @xmath87 and @xmath88 are the characteristic absorption frequencies and absorption strengths in the infrared ( ultraviolet ) range . these constants were calculated in refs . @xcite from the available optical data , and were summarized for the ordinary and extraordinary directions of quartz , calcite and barium titanate in table i of ref . the resulting dielectric functions are plotted in fig . [ fig2 ]  @xcite . from this figure , one sees that quartz shows the weakest birefringence and has @xmath89 for all frequencies . barium titanate is the strongest birefringent material and has @xmath90 , whereas calcite is somewhere in the middle and shows @xmath89 for small frequencies and @xmath90 for larger frequencies . we have calculated the casimir energy @xmath91 as a function of the angle @xmath92 between the magnetization of the magnetic mirror ( assumed to be in - plane ) and the optical axis of the non - magnetic one . results are shown in fig . [ fig3 ] for the three considered birefringent materials , and at distances of @xmath62 and @xmath93 . we have plotted the quantity @xmath94 so that the energy minimum will always lie at zero . all curves could be fitted very well with either a @xmath95 or a @xmath96 dependence , as depicted by the solid line connecting the calculated data points . we observe that at the shortest distance ( @xmath62 in fig . [ fig3](a ) ) , the curves for @xmath85 and @xmath84 show a @xmath96 behavior and thus favor a magnetization that is perpendicular to the optical axis , whereas the curve for @xmath83 varies like @xmath95 and tends to align the magnetization with the optical axis . this difference is due to the different sign of the quantity @xmath97 for the three materials in fig . [ fig2 ] , and we find that the energy minimum corresponds to a situation where the magnetization is oriented perpendicular to the dielectric axis with the largest value of @xmath98 . one should note that although @xmath99 changes sign at some intermediate frequency for calcite , the higher frequencies ( where @xmath100 ) will dominate the integrals in eq . ( [ eq : casenergy ] ) at small distances , and calcite therefore shows the same angular dependence as @xmath85 in fig . [ fig3](a ) . interestingly , at a longer distance ( @xmath93 in fig . [ fig3](b ) ) , the angular dependence is exactly opposite , at least for quartz and barium titanate : we have a @xmath95 behavior for @xmath85 and a @xmath96 for quartz ( see fig .  [ fig3 ] ) so that now the energy is minimized for a magnetization parallel to the dielectric axis with the largest value of @xmath98 . in the previous section , we observed a similar change of the angular dependence for a setup with an optically isotropic mirror , and explained it by a competition between the polar and transverse kerr contributions . although the polar effect is zero in the setup we consider now ( the magnetization lies in the plane , so @xmath36 in eq .  ( [ eq : reflmatrixdr ] ) ) , the two remaining longitudinal and transverse contributions can still compete with each other and will lead to the observed angular crossover . it is however not possible to address both contributions separately , since mixed terms of the form @xmath101 will also appear in the energy expression . calcite is a special case , as it has the same @xmath96 behavior both in the limit of short and long distances ( see fig .  [ fig3 ] ) . in principle , the competition between longitudinal and transverse kerr effect contributions would lead to the same crossover as in the other two considered materials . however , at long distances the low frequencies will determine the integrals in eq . ( [ eq : casenergy ] ) . as depicted in fig . [ fig2](b ) , one has @xmath102 at these low frequencies , compared to @xmath100 for the higher frequencies which determined the integrals at short distances . it is this sign change of @xmath99 that will cancel the angular dependence crossover due to the kerr effect competition . nevertheless , there will be an intermediate distance regime where the higher frequencies with @xmath100 are still the most important , and the competition between the kerr contributions is decisive in changing the angular dependence from @xmath96 to @xmath95 . this becomes clearer in fig . [ fig4 ] , where the amplitude of the @xmath92 oscillations is shown . both the curves for @xmath85 and quartz show a single kink . at the kink , the @xmath95 and @xmath96 contributions to the anisotropy energy are of equal magnitude so the net energy amplitude goes to zero . this kink thus signals a crossover point , and as explained above , one clearly observes two such crossovers for calcite but only a single one for quartz and barium titanate . furthermore we see from fig . [ fig4 ] that , as expected , the anisotropy energy is proportional to the amount of birefringence of the considered material , and thus largest for @xmath85 and smallest for quartz . an experimental detection of these anisotropy effects will be quite challenging since the calculated energies are rather small . however , the studied system has the advantage to lend itself to different measurement setups . firstly , a force measurement could be done by evaporating a droplet of iron on a cantilever and placing it over a birefringent plate . a resonant rf - field can be used to rotate the magnetization , and the casimir force will result in oscillations at the resonance frequency of the cantilever . for an iron droplet with a curvature radius @xmath103 , the force amplitude @xmath104 can be estimated by using the force proximity theorem : @xmath105 . at a distance of @xmath106 , the force is @xmath107 for @xmath85 . force detection with an - resolution has already been shown in similar magnetic resonance force microscopy setups  @xcite . secondly , the anisotropy energy will give rise to a mechanical torque @xmath108 that seeks to rotate the plates towards their minimal energy position . for plates with a radius of @xmath109 , the torque at an interplate distance @xmath110 is @xmath111 . a possible setup to demonstrate such torques can be found in ref . @xcite , where both plates are immersed in a liquid one on top of the other . with a suitable selection of material parameters , a repulsive casimir force will balance gravity and a birefringent disk can be held floating above a hard ferromagnetic plate . the casimir torque will then rotate the plates back to their minimum energy orientation once they are brought out of their equilibrium position . a third measurement scheme would consist of detecting the magnetic anisotropy energy of the ferromagnetic plate directly with a magnetic measurement . by using a cylindrically shaped ferromagnetic mirror cut along the correct axis , both the shape and magnetocrystalline anisotropy can be made small . furthermore , no anisotropy contribution intrinsic to the ferromagnetic plate will depend on the orientation of the birefringent plate , so that parisitic effects can in general be filtered out efficiently . finally , we would like to mention that instead of using a birefringent material such as barium titanate to create the uniaxially anisotropic mirror , the optical anisotropy can also be realized by means of some appropriate nanostructuring , e.g. , by creating parallel stripes on an isotropic gold mirror . the optical anisotropy can then be made much larger , thereby increasing the magnetic anisotropy effects ( a factor of @xmath112 is possible  @xcite ) . a detailed calculation for this case is not at all straightforward , and beyond the scope of the present paper . we have considered the casimir effect between a ferromagnetic and a non - magnetic mirror . if the non - magnetic mirror is optically isotropic , the casimir effect is shown to give rise to an out - of - plane magnetic anisotropy in the ferromagnetic plate . for a non - magnetic mirror with uniaxial anisotropy , we found an in - plane anisotropy . in both cases , the angular position with the lowest energy depends on the distance between the plates and is determined by the competition between different contributions to the magneto - optical kerr effect . different setups for measuring the effect have been proposed . metalidis2002 h.b.g . casimir , proc . akad . wet . * 51 * , 793 ( 1948 ) . m.  bordag , u.  mohideen , and v.m . mostepanenko , phys . rep . * 353 * , 1 ( 2001 ) . s.  k.  lamoreaux , phys . rev . lett . * 78 * , 5 ( 1997 ) . u.  mohideen , and a.  roy , phys . * 81 * , 4549 ( 1998 ) . h.  b.  chan , v.  a.  aksyuk , r.  n.  kleinman , d.  j.  bishop , and f.  capasso , science * 291 * , 1941 ( 2001 ) . g.  bressi , g.  carugno , r.  onofrio , and g.  ruoso , phys . lett . * 88 * , 041804 ( 2002 ) . r.  s.  decca , d.  lopez , e.  fischbach , and d.  e.  krause , phys . * 91 * , 050402 ( 2003 ) . h.  b.  chan , y.  bao , j.  zou , r.  a.  cirelli , f.  klemens , w.  m.  mansfield , and c.  s.  pai , phys . * 101 * , 030401 ( 2008 ) . j.  n.  munday , f.  capasso , and v.  a.  parsegian , nature * 457 * , 170 ( 2009 ) . g.  l.  klimchitskaya , and v.  m.  mostepanenko , contemp . phys . * 47 * , 131 ( 2006 ) . y.  s.  barash , izv . uchebn . zaved . * 16 * , 1227 ( 1973 ) [ sov . * 16 * , 945 ( 1975 ) ] s.  j.  van  enk , phys . a * 52 * , 2569 ( 1995 ) . j.  n.  munday , d.  iannuzzi , y.  barash , and f.  capasso , phys . a * 71 * , 042102 ( 2005 ) . shao , a .- h . tong , and j.  luo , phys . a * 72 * , 022102 ( 2005 ) . r.  b.  rodrigues , p.  a.  maia  neto , a.  lambrecht , and s.  reynaud , europhys . * 76 * , 822 ( 2006 ) . f.  s.  s.  rosa , d.  a.  r.  dalvit , and p.  w.  milonni , phys . a * 78 * , 032117 ( 2008 ) . g.  deng , z .- z . liu , and j.  luo , phys . a * 78 * , 062111 ( 2008 ) . p.  bruno , phys . rev . lett . * 88 * , 240401 ( 2002 ) . g.  metalidis and p.  bruno , phys . a * 66 * , 062102 ( 2002 ) . jaekel and s.  reynaud , j. physique * i 1 * , 1395 ( 1991 ) . g.  metzger , p.  pluvinage , and r.  torguet , ann . 10 * , 5 ( 1965 ) . _ handbook of optical constants of solids _ , edited by e.d . palik ( academic press , new york , 1991 ) . landolt - boernstein , _ numerical data and functional relationships in science and technology _ , new  series  group , iii vol.19 , pt . a ( springer , berlin , 1992 ) . a.  lambrecht and s.  reynaud , eur . j. d * 8 * , 309 ( 2000 ) . d.  rugar , o.  zger , s.  hoen , c.s . yannoni , h .- vieth , and r.d . kendrick , science * 264 * , 1560 ( 1994 ) . stowe , k.  yasumura , t.w . kenny , d.  botkin , k.  wago , and d.  rugar , appl . . lett . * 71 * , 288 ( 1997 ) . d.  rugar , r.  budakian , h.  j.  mamin and b.  w.  chui , nature * 430 * , 329 ( 2004 ) . t.  p.  sosnowski , opt . commun . * 4 * , 408 ( 1972 ) . j.  lekner , j.  phys . : condens . matter * 3 * , 6121 ( 1991 ) . r.  h.  muller , surf . science * 16 * , 14 ( 1969 ) . d.  b.  hough and l.  r.  white , adv . colloid interface sci . * 14 * , 3 ( 1980 ) . l.  bergstrm , adv . colloid interface sci . * 70 * , 125 ( 1997 ) . a similar figure can be found in ref . @xcite , but we have chosen to include it here again since the following discussion will heavily depend on it .
Abstract:

Answer: we consider the casimir interaction between a ferromagnetic and a non - magnetic mirror , and show how the casimir effect gives rise to a magnetic anisotropy in the ferromagnetic layer . the anisotropy is out - of - plane if the non - magnetic plate is optically isotropic . if the non - magnetic plate shows a uniaxial optical anisotropy ( with optical axis in the plate plane ) , we find an in - plane magnetic anisotropy . in both cases , the energetically most favorable magnetization orientation is given by the competition between polar , longitudinal and transverse contributions to the magneto - optical kerr effect , and will therefore depend on the interplate distance . numerical results will be presented for a magnetic plate made out of iron , and non - magnetic plates of gold ( optically isotropic ) , quartz , calcite and barium titanate ( all uniaxially birefringent ) .