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: torque oscillator ( sto ) is a nanoscaled device , utilizing a dc current through a magnetized magnetic multi - layer structure to yield a steady - state voltage oscillation @xcite-@xcite . the operating frequency of this voltage oscillation typically lies in the microwave range @xcite and can be widely tuned by altering the magnetic field and the dc current . the typical structure of an sto is presented in fig . it is composed of two magnetic layers , the `` free '' layer ( fl ) and the `` polarizing '' layer ( pl , or so - called `` fixed '' layer ) , being decoupled by a non - magnetic ( nm ) `` spacer '' . in this multi - layer structure , a spin - polarized current transfers angular momentum from the pl to the fl , so as to exert a torque on the local magnetization of the fl . by using this transferred torque to compensate the magnetic damping of the fl , magnetization dynamics can be sustained at microwave frequencies . the magnetization dynamics of the fl in the presence of spin - polarized current can be described by the landau - lifshitz - gilbert equation with a slonczewski spin - transfer term ( llgs ) @xcite , @xcite @xmath0+\frac{\alpha(\xi)}{m_\text{0}}[{{\mathbf{m}}}\times{\frac{{{\mathrm{d}}}{{\mathbf{m}}}}{{{\mathrm{d}}}t}}]+{\gamma}{{\boldsymbol{\tau_\text{stt } } } } \label{llgeq}\end{aligned}\ ] ] where @xmath1 is the gyromagnetic ratio , @xmath2 is the damping parameter , @xmath3 is the saturation magnetization , @xmath4 is the magnetization of the fl , @xmath5 is the effective magnetic field acting on the fl , and @xmath6 is the spin transfer torque ( stt ) , which is used to cancel out the damping term so as to achieve a steady precession of the fl magnetization . @xmath6 can be expanded as @xmath7 + b_\text{j}{{\mathbf{m}}}\times{{\mathbf{m_\text{p}}}}$ ] , where @xmath8 is the magnetization of the pl , @xmath9 and @xmath10 are the bias - dependent coefficients of the in - plane and perpendicular torque , respectively . the stt has been used in magnetic tunnel junction ( mtj ) devices and modeled @xcite , @xcite to develop several novel spintronic devices . for example , the stt - based mtj has been employed in logic circuit designs @xcite , @xcite , and in the spin transfer torque magnetoresistive random access memory ( stt - mram ) @xcite , @xcite , which outperforms sram in terms of power consumption and cost , as well as being non - volatile . the stt mtj models @xcite , @xcite used in logic circuits or memory applications , predict the switching between `` 0 '' and `` 1 '' , then set either the dc resistance or the dc voltage of the mtj sto accordingly . they do not contain the rf dynamics of stos and hence are not suitable to present stos behavior . in stos , the magnetization dynamics are detected by means of the magnetoresistance ( mr ) effect . this effect predicts a change in the resistance of the multi - layer depending on the relative orientation between the fl and pl . when stt stabilizes a precession of the fl magnetization , an oscillatory resistance is also established leading to an rf voltage generation by virtue of ohm s law . this voltage oscillation generated by the sto can be expressed as @xmath11 where @xmath12 is the dc resistance of the sto under a specific biasing condition , @xmath13 is the applied current for driving the sto , @xmath14 is the amplitude of the resistance oscillation , @xmath15 is the frequency generated by the sto and @xmath16 represents the random phase fluctuation ( i.e. phase noise ) of the sto . stos exhibit a unique blend of features : high operating frequency , extremely wide tunability , high integration level with cmos technology , and fast turn - on time ( < 1 ns ) @xcite . thanks to these features , stos show potential to be used as up - converting mixers without the need for a local oscillator @xcite , frequency detectors @xcite , magnetic field sensors @xcite , and oscillators @xcite . currently , stos using an insulator as the spacer , in contrast to those using thin metallic spacers , offer higher output power and hence are more suitable for applications . the sto implemented with the thin insulator forms a tunnel junction between the magnetic layers @xcite , and it is the so - called mtj sto . the mtj sto is used as the base in this work .    to enable the use of mtj stos in applications , a model is necessary to capture their characteristics and to be used in system- and circuit - level designs . the core of sto modeling is to solve the llgs . one possible approach to solve the llgs is based on micromagnetics @xcite , which describe the magnetization dynamics of the sto on a microscopic scale . however , this approach is not possible to implement using a hardware description language , such as verilog - a , due to the required complicated numeric computation . another approach is based on the macrospin approximation , which assumes that only a spatially - uniform magnetization precession is excited and the spin - polarized current is uniform across the area of the free layer @xcite . the macrospin - based sto model is able to provide an analytic solution with acceptable accuracy , while allowing implementation in a hardware description language . as a result , the macrospin - based analytical model can be employed in the design of sto - based systems . two such models , which provide the rf characteristics of the mtj sto , have been proposed in @xcite . however , these models contain equations that are usable for matching one specific device and are not fully verified by either experiments or theory . additionally , the dc operating point of these models has not been analyzed , limiting the accuracy of the model for circuit- or system - level design . furthermore , only one specific device has been used to verify the models , which is not sufficient . this paper presents a new comprehensive and compact mtj sto model , which can overcome the issues of the existing mtj sto models and be applied to mtj stos with arbitrary parameters . our model is based on the hamiltonian formalism presented in @xcite as well as extensive analysis of the literature , and it is completed by estimations including but not limited to the dc operating point and electrical rf power . we verify our model by three different mtj sto measurements under different possible biasing conditions , published by different groups . all the characteristics of the proposed model follow the measurements closely . furthermore , the proposed analytical model provides a comparable accuracy with the micromagnetics - based model . thereafter , it is implemented in verilog - a , encapsulating the characteristics of mtj stos , and is ready for being used by the device and circuit community to implement sto - based systems . this paper is organized as follows . in part i , the theoretical analyses of the effective magnetic field and the characteristics of the proposed mtj sto model are provided . part ii describes the implementation of the proposed model in verilog - a and presents the simulation results of the stand - alone mtj sto model . finally , the proposed model is further validated in an sto - based system , where the sto verilog - a model is simulated together with cmos rf circuits . the actual structure of mtj stos is complex and can vary in different mtj stos @xcite-@xcite , as it is given in . however , the generalized stack structure of the mtj sto , as it is shown in fig . 1(a ) and considered in @xcite , is identical for different mtj stos . based on this stack structure , an analytical model of the mtj sto using the actual cross section size , will be proposed and then verified by comparing it against three different mtj stos @xcite-@xcite provided by different research groups . a summary of the important parameters used in this analytical model to obtain the characteristics of these three mtj stos @xcite-@xcite is presented in table ii . the values of most of the parameters in table ii are available in @xcite-@xcite as either measured or suggested values . most of these values are directly employed in this work . nevertheless , some values that are provided by @xcite-@xcite are slightly adjusted so as to match the analytical model to experiments under different biasing conditions . the adjustments are less than 25@xmath17 of the measured or suggested values . besides , the values of some parameters are not given in @xcite-@xcite , as noted in table ii , so that the empirical values from literature are used . @xmath18 and @xmath19 are the only phenomenological parameters and will be detailed later . it should be mentioned that the notations of parameters and sto output characteristics , as well as the notations of external and effective magnetic fields , used in different publications @xcite-@xcite are not identical . to obtain the characteristics of the sto , such as operating frequency , output power , and linewidth ( the full width at half - maximum ) , according to the macrospin - based analytical model in @xcite , the first step is to calculate the effective field . .actual structures of different mtj stos [ cols="<,<,<",options="header " , ]     values without notes are either measured or suggested values in @xcite-@xcite . these values are slightly adjusted based on the measured or suggested values given in @xcite-@xcite . * empirical values .    * * phenomenological parameters .    the effective field , based on @xcite , can be solved by the boundary conditions derived from the integral formulation of maxwell s equations for magnetic fields . however , the boundary conditions ( eq.(102a-102c ) in @xcite ) take only the magnetic anisotropy field @xmath20 into account for simplicity , which is not sufficient to achieve an accurate mtj sto model since the perpendicular torque and the inter - layer coupling field @xmath21 @xcite-@xcite are considerable in the mtj sto . the perpendicular torque introduces an additional field with the amplitude of @xmath22 ( see eq.1 ) . consequently , the boundary conditions including @xmath22 and @xmath21 ( @xmath22 and @xmath21 are both along @xmath23 axis ) , along @xmath23 , @xmath24 , @xmath25 axes can be derived as    @xmath26    @xmath27    @xmath28    where @xmath29 , @xmath30 and @xmath31 are the magnitude , out - of - plane and in - plane angles of the effective field ; @xmath32 , @xmath33 and @xmath34 are the magnitude , out - of - plane and in - plane angles of the external applied field , as illustrated in fig . the perpendicular torque coefficient @xmath35 , as detailedly studied in @xcite is a function of the voltage @xmath36 across the sto . it is usually determined from the comparison of the theoretical and experimental characteristics . in @xcite , the expression @xmath35=37@xmath37 with the unit of oe , is used to estimate the additional field introduced by the perpendicular torque for the free - layer excitation mode , where the edge mode is neglected . although an additional quadratic term has been previously reported @xcite , its magnitude is relatively small and does not show any significant impact when included in our model . because of this and the uncertainty in the exact value for this quadratic coefficient , only the linear dependence is implemented . this linear dependence from @xcite is found to accurately match the characteristics of the mtj stos also in @xcite and @xcite . for the sake of generality , the proposed model allows users changing the coefficients of @xmath35 easily , for the cases where the above expression can not be successfully applied .    for mtj stos , the external field is typically applied in - plane ( @xmath33 = 0 ) @xcite so that the boundary conditions ( eq.(3a-3c ) ) can be reduced to    @xmath38    @xmath39    where eq.(4a ) gives the solution to @xmath31 , from which @xmath40 can be obtained from eq.(4b ) . these simplified equations enable easy implementation of the equation solver in verilog - a and rapid simulations , so that they are employed to obtain the effective field of mtj stos . the electrical signal generated by an mtj sto comprises a dc and an ac ( oscillating ) component . the dc component of the mtj sto model can be expressed by either the dc voltage across the mtj sto or its dc resistance , since the dc current is applied externally . the ac ( oscillating ) component is characterized by the operating frequency , peak power ( or amplitude ) and linewidth .      in the existing models @xcite , the mtj stos are assumed to be in the anti - parallel state ( with a fixed dc resistance of @xmath41 ) , and the dc operating point is not analyzed . nonetheless , in order to compute the other characteristics of the mtj sto and enable the design of on - chip biasing circuits for it , the dc operating point of the mtj sto is of importance . the dc voltage across the mtj sto can be simply calculated by @xmath42 , where @xmath12 is given by @xcite @xmath43 where @xmath41 and @xmath44 are the resistances of the mtj sto in anti - parallel and parallel states respectively . as it can be noticed from eq.(5 ) , the dc operating point ( or dc resistance ) of the mtj sto can vary greatly as @xmath31 changes . the operating frequency @xmath15 , according to @xcite , is given by @xmath45 where @xmath46 and @xmath47 are the ferromagnetic resonance ( fmr ) frequency and coefficient of the nonlinear frequency shift respectively , which can be obtained from eq.(103 ) , eq.(104a ) and eq.(105a ) in @xcite once the effective field is determined ; @xmath48 is the dimensionless power , which is a coefficient that is proportional to the experimentally measured rf power of the sto @xcite and is involved in calculating all the characteristics of mtj stos , hence it is a critical quantity . as detailed in @xcite ( eq.(84b ) ) , @xmath48 is a function of the supercriticality @xmath49 , nonlinear damping coefficient @xmath50 and noise power @xmath19 . @xmath49 is defined as the ratio between @xmath13 and the threshold current @xmath51 of the mtj sto . @xmath51 ( see eq.(22 ) in @xcite ) is a function of @xmath31 , the nonlinear damping rate @xmath52 calculated based on the effective field ( see eq.(104b ) in @xcite ) , and the coefficient @xmath53 expressed as @xcite @xmath54 where @xmath55 is the dimensionless spin polarization efficiency , @xmath56 is the spectroscopic land@xmath57 factor , @xmath58 is the bohr magneton , @xmath59 is the modulus of the electron charge , @xmath60 is the thickness of fl , and @xmath61 is the area of the current - carrying region . for mtj stos with cofe and cofeb as the pl , @xmath55 is @xmath62 @xmath63 @xmath64 and @xmath65 @xmath63 @xmath66 , respectively @xcite . to obtain the other important parameter @xmath50 that is required to solve @xmath48 , @xmath18 involved in the expression of @xmath50 ( see eq.(105b ) ) needs to be determined . @xmath18 is the first coefficient in the expansion of @xmath2 and is usually considered as a phenomenological parameter @xcite . the condition for @xmath18 is that the resulting @xmath50 falls within its typical range @xmath67 @xcite . likewise , the noise power @xmath19 is typically treated as a constant value between 0 and 0.2 @xcite . to determine @xmath18 and @xmath19 , they need to be roughly adjusted and tested together in order to enable a good agreement of the mtj sto characteristics between the proposed model and the experimental results .    to validate the operating frequencies of the proposed mtj sto model , we compare the numerical results to the experimental measurements provided by @xcite-@xcite under different biasing conditions , and the comparison is depicted in fig . 2 .    fig . 2(a ) shows the comparison between the model and the measured @xmath15 ( @xcite ) at @xmath68 = 300 oe , with respect to @xmath34 as well as @xmath13 . for @xmath31 < 70@xmath69 , as @xmath13 increases , a decreasing @xmath15 is observed in both the modeled and the measured results as presented in fig . at @xmath31 = 70@xmath69 and 90@xmath69 , @xmath15 are similar for different @xmath13 . this is mainly because the nonlinear frequency shift @xmath47 in eq.(6 ) is close to zero at these angles , so that the operating frequencies for different @xmath13 are mostly determined by the fmr . this observation is in agreement with the measured @xmath47 @xcite . generally , the operating frequency of the proposed model matches the measured data in a large region . the discrepancy between the model and the measured data may be due to the simplification of the macrospin - based model , which is not able to represent the complexities of a real device such as microscopic dynamics , roughness and the coupling between layers .    to further verify the dependence of the external magnetic field on @xmath15 of the proposed model , @xmath15 of the model as a function of the external field and the biasing voltage @xmath70 , is compared with the measured @xmath15 given in @xcite , and illustrated in fig . the results obtained from the model matches well with the measured data . for low magnetic fields , a small discrepancy between the model and measured results can be found . this could be due to the fact that the layers are not fully saturated under low magnetic fields . the measured @xmath15 in @xcite and @xmath15 obtained from the micromagnetics simulation results @xcite , are also used to compare with @xmath15 of the model . in @xcite , the mtj sto was measured under a fixed external field of 200  oe , so that the comparison is made under this condition . the comparison is shown in fig . 2(c ) . within 0@xmath69 < 50@xmath69 , the proposed model offers good agreement with the results obtained from the micromagnetic - based model . however , similar discrepancies are identified between the measured results and their respective models . this may be caused by the multiple modes of the sample , where the frequencies of the first and second modes appear close to each other at low field magnitudes and low field angles @xcite . the multi - mode behavior is neither captured by the macrospin - based model , nor correctly estimated by the micromagnetics - based model . regarding the large field angles ( 50@xmath69 < @xmath34 < 90@xmath69 ) , inconsistency of measured results is found between @xmath13- and field - swept measurements . this might be explained by electromigration , which degrades the mtj sto during the measurement . in this case , no model can completely predict the operation of the mtj sto ; yet , both of the macrospin- and micromagnetics - based models provide reasonable agreement with @xmath13- and field - swept measurements respectively . the comparison in fig . 2(c ) indicates that the proposed analytical model is able to capture the core behaviors of the mtj sto in spite of the considerable complexities in the mtj sto , which are not yet totally identified and fully understood . in addition , the model provides the operating frequencies of mtj stos , which have similar accuracy as the ones obtained from the micromagnetic simulations . moreover , the operating frequencies measured by sweeping @xmath13 at a fixed external field angle @xmath34 = 40@xmath69 that are available in @xcite , are compared with the modeled results . for different @xmath13 , the proposed model matches the measured data with a slight frequency offset of about 0.15 ghz ( fig . 2(c ) ) , which is less than 5% of the operating frequency .      for estimating the peak power generated by the mtj sto , expressions determining the magnetoresistance and precession angle were introduced and used in the existing mtj sto models @xcite . these expressions , however , are not yet validated by either the theory or experiments . in this work , the expressions that have been typically employed to estimate the precession angle @xmath71 based on the measured power @xcite and validated by experiments @xcite , the peak power of the fundamental signal of the mtj sto is given by @xcite @xmath72 where @xmath73 is the rf power transfer efficiency at frequency @xmath74 and assumed to be 1 ( no loss due to parasitics ) , @xmath75 is the bessel function of the first kind , and @xmath76 is defined as @xmath77 to model the output power of the mtj sto , an expression for estimating @xmath71 is required to determine the peak power generated by the mtj sto using eq.(8 ) . @xmath71 can be estimated based on the equation provided by @xcite , which is re - written as @xmath78 eq.(10 ) gives the dependence of @xmath71 on the dimensionless power @xmath48 converted from the experimental results . here , the dimensionless power @xmath48 has already been obtained in the previous subsection based on the analytical theory . eq.(8 ) and eq.(9 ) contain bessel functions of the first kind , which are difficult to implement in verilog - a . since mtj stos have limited precession angles @xcite , based on the expansions of bessel functions , two approximations , including @xmath79 and @xmath80 , can be used to simplify eq.(8 ) and eq.(9 ) . in addition , in order to reproduce the voltage oscillation as described by eq.(2 ) in the verilog - a model , @xmath14 in eq.(2 ) can be derived and simplified based on eq.(8 ) as @xmath81 eq.(9 ) can be also reduced to @xmath82 the peak power and the voltage amplitude of the fundamental signal of the mtj sto can be determined based on eq.(11 ) and eq.(12 ) . the power generated by the second harmonic is not considered in this work , yet can be computed in a similar way and added to the proposed model . figure 3 presents the comparison of the peak power ( a ) and the peak power spectral density ( psd ) ( b ) between the proposed model and the available experimental results ( @xcite ) .    the output peak power of the proposed model , presented in fig . 3(a ) , is derived as @xmath83 and compared with the measured peak power @xcite . it follows the measured data @xcite under different biasing conditions , where the @xmath13 and @xmath31 are swept . the minor discrepancies between the modeled and measured peak power might be due to the imperfections of loss de - embedding during the measurement .    to further validate the peak power ( or voltage amplitude ) of the proposed model as a function of @xmath13 , the time domain signals generated by the proposed model are converted to the frequency domain to obtain their psds , so as to compare the peak psds against the measured ones in @xcite . to generate these signals for evaluating the amplitudes of the spectrum , the required information of linewidth is taken from the measured data @xcite . in addition , as it can be calculated from @xcite , the sum of peak psds of the second and third modes takes approximately 1/3 of the total peak psds of all the modes . accordingly , a coefficient of 0.67 is used to downscale the modeled peak psds , in order to compare the modeled peak psds of the first mode with the measured ones when multiple modes exist . the comparison of the modeled and measured peak psds is shown in fig . 3(b ) . for different @xmath13 , the peak psds obtained from the proposed mtj sto model are very close to those given by the experiments @xcite .      in the existing mtj sto models @xcite , the equation used to obtain the linewidth @xmath84 is given by eq.(95 ) in @xcite @xmath85 where @xmath86 is the normalized dimensionless nonlinear frequency shift , @xmath87 is the positive damping rate , @xmath88 is the product of the boltzmann constant and the temperature , and @xmath89 is the oscillation energy . @xmath86 , @xmath87 and @xmath89 can be computed based on the known parameters and the solved nonlinear coefficients ( see eq.(33 ) , eq.(19b ) and eq.(77 ) in @xcite ) . nonetheless , eq.(13 ) is only valid for the above - threshold regime , where @xmath13 is larger than @xmath51 @xcite . thus , it should not be applied to all the biasing conditions . regarding the below - threshold regime ( @xmath90 ) , @xmath91 is found to be @xcite @xmath92 noticed from eq.(14 ) , as the mtj sto operates in the near - threshold regime ( @xmath93 ) , @xmath84 goes to zero , which can not be achieved in practical cases . our calculations suggest that eq.(14 ) is not valid for the regime @xmath94 , since the computed linewidth is much lower than practical values . hence , the regime @xmath94 is considered as the near - threshold regime in this work . however , as explained in @xcite , the linewidth of an mtj sto operating in near - threshold regime , is difficult to approximate theoretically . fortunately , our calculations show that in this near - threshold regime , using eq.(13 ) can offer reasonable agreement between the modeled and the measured linewidth . since this work is targeting mtj sto modeling rather than exploring the linewidth of the mtj sto , eq.(13 ) is simply applied to approximate the linewidth in the near - threshold regime , even though , it is not theoretically validated .    to verify the linewidth of the proposed model , the available measurement results of two different mtj stos @xcite are employed . figure 4(a ) shows the comparison between the measured @xmath84 given in @xcite and the @xmath84 obtained from the proposed model , as a function of @xmath34 and @xmath13 . for different angles , the general trend of the modeled @xmath84 follows the trend of the measured @xmath84 , owing to the different equations used for computing the linewidth in different regimes . besides , for low currents ( @xmath13 = 0.5 ma ) , the linewidth of the proposed model shows an increasing trend as a function of the field angle , indicating that the mtj sto is operating in the below - threshold regime . the threshold point shifts towards higher field angles as @xmath13 increases . figure 4(a ) demonstrates that the proposed model is able to offer the linewidth with the correct order of magnitude and proper trend regardless of the biasing conditions , despite the huge simplification of linewidth estimation provided by the macrospin approach . the linewidth of the proposed model is also compared with the measured data provided in @xcite , and depicted in fig . 4(b ) . as @xmath13 increases , the measured @xmath84 remains almost flat , while the @xmath84 of the model is decreasing . the trend of the linewidth obtained from the model , yet , agrees with the theoretical analysis in @xcite . the disagreement is likely due to the fact that current induces stronger dynamics , which physically translates to the generation of higher order spin wave modes and possibly mode - hopping events that broaden the linewidth @xcite . for mtj stos with sufficiently large cross section , such modes can only be represented micromagnetically , affecting the accuracy of the macrospin approach . moreover , mtj stos have been reported to exhibit not only the analytically treated white frequency noise but also colored , @xmath95 type noise for low fluctuation frequencies @xcite . frequency noise of the @xmath95 type has been shown to lead to an increased measured linewidth @xcite , making it another plausible explanation of the discrepancy between the measured data and the proposed model . a comprehensive and compact analytical mtj sto model based on the macrospin approximation as well as the physics - based equations of the sto , has been proposed . the perpendicular torque and inter - layer coupling field that significantly affect the characteristics , have both been considered in the proposed model to calculate all the characteristics of mtj stos . the model has been compared under different biasing conditions against the experimental data obtained from three different mtj stos . despite its simplicity , the model can reproduce the experimental data with an acceptable degree of accuracy . therefore , it is suitable for being implemented in a hardware description language , which enables the evaluation and utilization of mtj stos in real and extensive applications . s. i. kiselev , j. c. sankey , i. n. krivorotov , n. c. emley , r. j. schoelkopf , r. a. buhrman , and d. c. ralph , `` microwave oscillations of a nanomagnet driven by a spin - polarized current , '' _ nature _ , vol . 425 , no . 6956 , pp . 380383 , 2003 . g. d. panagopoulos , c. augustine , and k. roy , `` physics - based spice - compatible compact model for simulating hybrid mtj / cmos circuits , '' _ ieee trans . electron devices _ , vol . 9 , pp . 28082814 , 2013 . y. zhang , w. zhao , y. lakys , j. o. klein , j. v. kim , d. ravelosona , and c. chappert , `` compact modeling of perpendicular - anisotropy cofeb / mgo magnetic tunnel junctions , '' _ ieee trans . electron devices _ , vol . 3 , pp . 819826 , 2012 . x. yao , j. harms , a. lyle , f. ebrahimi , y. zhang , and j .- p . wang , `` magnetic tunnel junction - based spintronic logic units operated by spin transfer torque , '' _ ieee trans . 1 , pp . 120126 , 2012 . c. augustine , n. n. mojumder , x. fong , s. h. choday , s. p. park , and k. roy , `` spin - transfer torque mrams for low power memories : perspective and prospective , '' _ ieee sens . j. _ , vol . 12 , no . 4 , pp . 756766 , 2012 . s. bonetti , p. muduli , f. mancoff , and j. @xmath96kerman , `` spin torque oscillator frequency versus magnetic field angle : the prospect of operation beyond 65 ghz , '' _ appl . 10 , p. 102507 , 2009 . i. n. krivorotov , n. c. emley , j. c. sankey , s. i. kiselev , d. c. ralph , and r. a. buhrman , `` time - domain measurements of nanomagnet dynamics driven by spin - transfer torques , '' _ science _ , 5707 , pp . 228231 , 2005 . p. k. muduli , y. pogoryelov , s. bonetti , g. consolo , f. mancoff , and j. kerman , `` nonlinear frequency and amplitude modulation of a nanocontact - based spin - torque oscillator , '' _ phys . b _ , vol . 81 , no . 14 , p. 140408 , 2010 . a. a. tulapurkar , y. suzuki , a. fukushima , h. kubota , h. maehara , k. tsunekawa , d. d. djayaprawira , n. watanabe , and s. yuasa , `` spin - torque diode effect in magnetic tunnel junctions , '' _ nature _ , vol . 7066 , pp . 339342 , 2005 . p. m. braganca , b. a. gurney , b. a. wilson , j. a. katine , s. maat , and j. r. childress , `` nanoscale magnetic field detection using a spin torque oscillator , '' _ nanotechnology _ , vol . 23 , p. 235202 , 2010 . t. chen , p. drrenfeld , s. rodriguez , j. kerman , and a. rusu , `` a highly tunable microwave oscillator based on mtj sto technology , '' _ microw . 9 , pp . 20922095 , 2014 . i. firastrau , d. gusakova , d. houssameddine , u. ebels , m .- c . cyrille , b. delaet , b. dieny , o. redon , j .- ch . toussaint , and l. d. buda - prejbeanu , `` modeling of the perpendicular polarizer - planar free layer spin torque oscillator : micromagnetic simulations , '' _ phys . b _ , vol . 78 , no . 2 , p. 024437 , 2008 . h. lim , s. ahn , m. kim , s. lee , and h. shin , `` a new circuit model for spin - torque oscillator including perpendicular torque of magnetic tunnel junction , '' _ adv . condens . matter phys . 2013 , p. 169312 , 2013 . m. kim , h lim , s. ahn , s. lee , and h. shin , `` advanced circuit - level model of magnetic tunnel junction - based spin - torque oscillator with perpendicular anisotropy field , '' _ j. semicond . _ , vol . 13 , no . 6 , pp 556561 , 2013 .    p. k. muduli , o. g. heinonen , and j. kerman , `` bias dependence of perpendicular spin torque and of free - and fixed - layer eigenmodes in mgo - based nanopillars , '' _ phys . b _ , vol . 18 , p. 184410 , 2011 . z. m. zeng , p. upadhyaya , p. k. amiri , k. h. cheung , j. a. katine , j. langer , k. l. wang , and h. w. jiang , `` enhancement of microwave emission in magnetic tunnel junction oscillators through in - plane field orientation , '' _ appl . _ , vol . 99 . 3 , p. 032503 , 2011 . p. gangmei , a. neudert , m. k. marcham , r. j. hicken , m. a. gubbins , x. cao , r. r. lamberton , and a. b. johnston , `` thermal and spin - transfer - torque excitation of precessional modes in magnetic tunnel junction nanopillars with symmetric interfaces and a thick free layer , '' _ phys . b _ , vol . 88 , no . 13 , p. 134415 , 2013 . h. kubota , a. fukushima , k. yakushiji , t. nagahama , s. yuasa , k. ando , h. maehara , y. nagamine , k. tsunekawa , d. d. djayaprawira , n. watanabe and y. suzuki , `` quantitative measurement of voltage dependence of spin - transfer torque in mgo - based magnetic tunnel junctions , '' _ nature phys . 1 , pp . 3741 , 2007 . a. m. deac , a. fukushima , h. kubota , h. maehara , y. suzuki , s. yuasa , y. nagamine , k. tsunekawa , d. d. djayaprawira , and n. watanabe , `` bias - driven high - power microwave emission from mgo - based tunnel magnetoresistance devices , '' _ nature phys . 4 , no . 10 , pp . 803809 , 2008 . j. c. sankey , y .- t . cui , j. z. sun , j. c. slonczewski , r. a. buhrman , and d. c. ralph , `` measurement of the spin - transfer - torque vector in magnetic tunnel junctions , '' _ nature phys . 1 , pp . 6771 , 2007 . t. chen , a. eklund , e , iacocca , s. rodriguez , g. malm , j. kerman and a. rusu , `` comprehensive and macrospin - based magnetic tunnel junction spin torque oscillator model  part ii : verilog - a model implementation , '' submitted to _ ieee trans . electron devices _ , 2014 . a. c. reilly , w. park , r. slater , b. ouaglal , r. loloee , w.p pratt jr . , and j. bass,``perpendicular giant magnetoresistance of co91fe9/cu exchange - biased spin - valves : further evidence for a unified picture , '' _ j. magn l269l274 , 1999 . h. oshima , k. nagasaki , y. seyama , y. shimizu , s. eguchi , and a. tanaka , `` perpendicular giant magnetoresistance of cofeb / cu single and dual spin - valve films , '' _ j. appl 10 , pp . 81058107 , 2002 . c. boone , j. a. katine , j. r. childress , j. zhu , x. cheng , and i. n. krivorotov , `` experimental test of an analytical theory of spin - torque - oscillator dynamics , '' _ phys . b _ , vol . 79 , no . 14 , p. 140404 , 2009 . m. quinsat , d. gusakova , j. f. sierra , j. p. michel , d. houssameddine , b. delaet , m. -c . cyrille , u. ebels , b. dieny , l. d. buda - prejbeanu , j. a. katine , d. mauri , a. zeltser , m. prigent , j. -c . nallatamby and r. sommet , `` amplitude and phase noise of magnetic tunnel junction oscillators , '' _ appl . 18 , p. 182507 , 2010 . m. w. keller , m. r. pufall , w. h. rippard , and t. j. silva , `` nonwhite frequency noise in spin torque oscillators and its effect on spectral linewidth , '' _ phys . b _ , vol . 82 , no . 5 , p. 054416 , 2010 . tingsu chen ( s11 ) received the b.sc . degree in communication engineering from the nanjing university of information science and technology , china , and the m.sc . degree in system - on - chip design from the kth royal institute of technology , sweden , in 2009 and 2011 , respectively . she is currently working toward the ph.d . degree at kth with the research area of high frequency circuit design for spin torque oscillator technology . anders eklund ( s13 ) received the m.sc . degree in engineering physics from kth royal institute of technology , sweden , in 2011 . he is currently working towards a ph.d . degree in physics at kth , experimentally investigating the frequency stability of spin torque oscillators by means of electrical characterization and synchrotron x - ray measurements . ezio iacocca ( s08 ) received the b.sc . degree in electronic engineering from the simn bolvar university , caracas , venezuela ( 08 ) , the m.sc . in nanotechnology from the royal institute of technology , stockholm , sweden ( 10 ) , and the ph.d . in physics from the university of gothenburg , gothenburg , sweden ( 14 ) . his research focuses on the magnetodynamical modes of spin transfer torque driven nano oscillators and their applications in communication and storage technology . saul rodriguez ( m06 ) received the b.sc . degree in electrical engineering from the army polytechnic school ( espe ) , quito , ecuador , and the m.sc . degree in system - on - chip design and the ph.d . degree in electronic and computer systems from the kth royal institute of technology , stockholm , sweden . in 2001 , 2005 , and 2009 , respectively . his current research interests include rf cmos circuit design for wideband frond - ends , ultralow - power circuits for medical applications and graphene - based rf , and ams circuits . b. gunnar malm ( m98 - sm10 ) was born in stockholm , sweden , in 1972 . he received the m.s . from uppsala university , sweden , in 1997 , the phd in solid - state electronics 2002 , from royal institute of technology ( kth ) , stockholm . he is an associate professor at the school of ict , kth since 2011 . his recent work includes silicon photonics , silicon - carbide technology for extreme environments and spintronics . he also serves on the kth sustainability council . johan @xmath96kerman ( m06 ) has an ing . degree ( 94 ) from epfl , switzerland , a m.sc . in physics ( 96 ) from lth , sweden , and a ph.d . in materials physics ( 00 ) from kth royal institute of technology , stockholm . in 2008 he was appointed full professor at university of gothenburg and is a guest professor at kth royal institute of technology . he is also the founder of nanosc ab and nanosc instruments ab . ana rusu ( m92 ) received the m.sc . degree in electronics and telecommunications and ph.d . degree in electronics in 1983 and 1998 , respectively . she has been with kth royal institute of technology , stockholm , sweden , since 2001 , where she is professor in electronic circuits for integrated systems . her research interests include low / ultralow power high performance cmos circuits and systems , sto - based systems , rf graphene and high temperature sic circuits .
Abstract:

Answer: magnetic tunnel junction ( mtj ) spin torque oscillators ( sto ) have shown the potential to be used in a wide range of microwave and sensing applications . to evaluate potential uses of mtj sto technology in various applications , an analytical model that can capture mtj sto s characteristics , while enabling system- and circuit - level designs , is of great importance . an analytical model based on macrospin approximation is necessary for these designs since it allows implementation in hardware description languages . this paper presents a new macrospin - based , comprehensive and compact mtj sto model , which can be used for various mtj stos to estimate the performance of mtj stos together with their application - specific integrated circuits . to adequately present the complete model , this paper is divided into two parts . in part i , the analytical model is introduced and verified by comparing it against measured data of three different mtj stos , varying the angle and magnitude of the magnetic field , as well as the dc biasing current . the proposed analytical model is suitable for being implemented in verilog - a and used for efficient simulations at , circuit- and system - levels . in part ii , the full verilog - a implementation of the analytical model with accurate phase noise generation is presented and verified by simulations . = 1    spin torque oscillator , magnetic tunnel junction , macrospin , analytical 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: in an hour s talk such as this it is impossible to cover all that is known on the subject , especially since in some respects detail is of the essence in dealing with exact solutions . the three major recent reviews of rather general character @xcite have a total of over 1200 pages and that is before moving to more specialized reviews such as @xcite . ( the short summary of what is new in @xcite as compared with the first edition , given below as an appendix , illustrates some ways in which the field has changed in the last 25 years . ) thus i shall be selective , dwelling at some length on particular solutions rather than trying to cram in as many as possible in the time . for example , i give particular attention to the schwarzschild solution , the first solution known that did not have constant curvature . i apologize to the authors of the many excellent pieces of work i do not touch on .    before learning the techniques for finding solutions of the einstein equations and discovering their properties , one should ask `` why this is a worthwhile endeavour ? '' . some colleagues do seem to regard it as something of a backwater in the theory . the reason it is still important stems from the nonlinearity of general relativity , one of its essential features . to understand the meaning of the theory , there are really three approaches . one can seek to prove global results , such as are described in @xcite and were the subject of the isaac newton institute programme in progress at the time of writing ( see http://www.newton.cam.ac.uk/webseminars/pg+ws/2005/gmr/ ) . one can try to use approximation , either in the form of iterated perturbation methods or numerical solutions . finally , one can use exact solutions : as mason and woodhouse said , `` they combine tractability with nonlinearity , so they make it possible to explore nonlinear phenomena while working with explicit solutions '' @xcite , and as we shall see below , they have had considerable impact on the theory . we also need to ask : `` what is a solution ? '' . we could assume a form for the metric , calculate its einstein tensor , and so obtain a form for the energy - momentum through einstein s equations . the pointlessness of such ` solutions ' was made clear by synge @xcite . more generally , choosing a more complicated form of energy - momentum with a simple form of metric usually reduces the number of equations to be solved and makes the task of finding solutions easier . the vacuum case , and cases with an equivalent set of equations to actually solve , are more difficult .    of course , exact solutions are very special cases . however , as bicak @xcite noted , feynman said `` the physicist is always interested in the special case . he is talking about something , he is not talking abstractly about anything . he wants to discuss the gravity law in three dimensions : he never wants the arbitrary force case in @xmath0 dimensions . so a certain amount of reducing is necessary ... '' . the special cases known exactly can be very useful as examples in guiding the global or approximative approaches . indeed without such uses finding exact solutions would be more like stamp - collecting than science . it requires ingenuity and the objects found may be beautiful , which is fine if you think relativity is a branch of pure mathematics , but if you think it is physics , that is not so good , despite feynman s point . it is not the formulae for the solutions that are really of interest : `` at present the main problem concerning solutions , in our opinion , is not to construct more but rather to understand more completely the known solutions with respect to _ their local geometry , symmetries , singularities , sources , extensions , completeness , topology and stability _ '' ( my emphasis ) . this remark by ehlers and kundt in 1962 @xcite is as true now as when it was written , perhaps more so as a result of the large number of further solutions found since then . part of my purpose here is to illustrate the crucial role such understanding of exact solutions has played in the development of our understanding of general relativity itself . the illustrations of uses of specific solutions given below naturally focus on some of the best - known solutions . one might imagine that the very large number of other solutions known have little use . for many that may be true , at least so far , but many others have helped the exploration of the physics . the effort required to find solutions is not trivial . compared with newtonian gravity , general relativity has one more independent variable , 9 more dependent ones ( taking the metric approach ) , and equations of degree 8 in these variables rather than one , the result of these changes being that the general form of the field equations expands to 10 partial differential equations each with thousands of terms ( in terms of the metric and coordinates ) . this complexity is the reason only some special solutions can be found , principally those with some special symmetry or algebraic property . one might nevertheless hope for a general solution . the closest to this which i know of ( the formalism of @xcite ) is however largely intractable . so i feel hopes for a useful general solution are slim and instead one has to make the simplifying assumptions already mentioned . here the choices have developed over the last 90 years . one disadvantage of those choices is that different ones may lead us to the same solution . i repeat a jeu desprit i first gave in a review some years ago @xcite , by giving a list of the most commonly - rediscovered solutions :    1 . flat space 2 . the schwarzschild solution 3 .   the kasner solutions 4 . plane waves 5 . conformally flat perfect fluids 6 . the taub - nut family of solutions 7 . static spherically symmetric perfect fluids 8 . cylindrically symmetric stationary electrovac solutions 9 . plane symmetric fluid and electrovac solutions 10 . spherically symmetric shearfree fluids    the 7th and subsequent places in this list are actually rather open to debate , for example the harris zund class are strong contenders . i will say something later about how to recognize known solutions . i should add that in this review i am going to stick firmly to the number of dimensions that i know i exist in , i.e.  4 , and avoid trespassing on the 5- and higher - dimensional work described by others . the development of exact solutions in such theories seems to me to be for the most part still at the stage of using only very simple metric forms . since the title of this meeting is ` a century of relativity ' , not the 90 of general relativity , i am obliged to begin at the beginning , meaning the first solution , the spacetime of special relativity , minkowski space : @xmath1 it is flat , empty and is usually taken to have @xmath2 topology .    the role of this solution has been considerable , for example :    * it provides the prototype for asymptotic flatness ( see ehlers contribution to this volume ) * it is ideal for ` cutting and pasting ' , which was a technique used to great effect while the concepts of causal structure were being developed * it was the setting for the first work on acceleration horizons ( the unruh effect ) * in suitable coordinates , it is the milne universe which is often used as the extreme robertson - walker model with @xmath3 . apart from the use of this example in astrophysical predictions , this second choice of slicing of flat space , together with the three foliations of de sitter space , helps to illustrate the fact that the open , flat or closed nature of spacetimes can depend on the slicing . * it is a special case of many other models , whence its frequent rediscovery * it is the main background for quantum field theory ( qft ) , and plays a role even for qft in curved spaces . the other solutions one could regard as in some sense known before gr was discovered were the other spaces of constant curvature , the de sitter and anti - de sitter spaces : @xmath4 ^ 2},\ ] ] where @xmath5 . these have played a role in , for example ,    * inflation ( for de sitter ) * the ads / cft correspondence ( anti de sitter ) * the `` no hair '' theorems ( de sitter ) * the development of our understanding of particle horizons and event horizons . bicak @xcite notes that stability against general non - linear vacuum perturbations has been proved for these three solutions . this is a ` robustness ' result : it raises the question of how much of what we do is robust in this sense . up to the 1950s the main method of finding new solutions started by postulating a coordinate form of the metric . authors generally assumed one of :    * spherical symmetry , * cylindrical symmetry , * staticity or stationarity and axisymmetry , or * a plane wave form .    in the large bibliography amassed as background for writing @xcite , only a part of which appears in the book itself , i found that nearly all papers before 1950 belonged to one of these groups , by far the largest group being the spherically symmetric cases . it should be noted that the metric forms were usually just written down , whereas now we would derive them from group - theoretic or other invariant assumptions .    as successful examples of this method , one can cite the work of schwarzschild , droste , levi - civita , kasner , chazy , curzon , brinkmann , baldwin and jeffrey , lewis , van stockum , papapetrou , and gdel , refs.@xcite  @xcite , and of course the robertson - walker metrics . many of these solutions have been important in later discussions , for example in establishing the reality of gravitational waves or elucidating the nature of directional singularities , but i will focus only on the two best - known .      the original form of the schwarzschild solution @xcite used a radial coordinate @xmath6 which had its origin at the horizon : for clarity i denote this @xmath7 below . schwarzschild also gave ( contrary to some statements in the literature ) the form most often quoted now : @xmath8 where @xmath9 . here i have renamed schwarzschild s @xmath10 as @xmath6 . schwarzschild regarded this as an auxiliary form because it did not fulfil the condition @xmath11 which einstein had imposed in the initial formulation of general relativity , a condition of course later discarded . this solution initiated a discussion ( see e.g.   @xcite ) on the meaning of the surface @xmath12 and the absence of a solution clearly analogous to a point mass in newtonian theory . there is of course no alternative point mass solution to consider due to the uniqueness of ( [ schw ] ) as the spherically symmetric vacuum solution ( birkhoff s theorem ) . there are still authors who argue that schwarzschild s original @xmath13 should be regarded as a singularity representing a point mass . the horizon clearly is not a regular point since the area of surrounding spheres has a limit @xmath14 ( which is part of the reason @xmath12 is now understood as a sphere , not a point ) . the argument can be made at various levels of sophistication , the simplest being that the metric components are singular ( similar things could be said of the axis of spherical polars , but nobody argues this is singular ! ) . to counter such arguments , it helps to note first that the horizon @xmath12 ( @xmath13 ) is not in the coordinate patch for ( [ schw ] ) . ( the waters here have been somewhat muddied by hilbert s arguments that these coordinates could be continued to the interior , in which he overlooked that at the horizon the three coordinates @xmath15 do not parametrize a three - dimensional manifold , as becomes immediately apparent on passing to the kruskal - szekeres picture . ) the whole schwarzschild patch @xmath16 is isometric to a region of the kruskal - szekeres solution ( region i in the conformal diagram given as fig . 1 ) , @xmath17 where @xmath6 is defined implicitly by the equations giving @xmath18 and @xmath19 in terms of the previous @xmath20 and @xmath6 : @xmath21 by inspection of this form of the solution we see that considering the bounding surface @xmath12 at a given time as a point amounts to topologically identifying all the points on a sphere . one can produce an analogous effect by cutting and pasting flat space ( in @xmath7-like coordinates ) . ( 8,5 ) ( 2.5,4)(1,0)3 ( 2.5,1)(1,0)3 ( 2.5,4)(1,-1)3 ( 2.5,1)(1,1)3 ( 1,2.5)(1,-1)1.5 ( 1,2.5)(1,1)1.5 ( 5.5,4)(1,-1)1.5 ( 5.5,1)(1,1)1.5 ( 5.6,2.6)i ( 2.5,2.5)ii ( 3.8,3.25)iii ( 3.8,1.5)iv ( 7.25,2.5)@xmath22 ( 5.6,4.1)@xmath23 ( 6.25,3.5)@xmath24 ( 6.25,1.5)@xmath25 ( 3.5,4.25)@xmath26 ( 3.5,0.5)@xmath26 ( 4,2.5)(5.5,3.5)(7,2.5 ) ( 4,2.5)(5.5,2.5)(7,2.5 ) ( 4,2.5)(5.5,1.5)(7,2.5 ) ( 5.5,1)(4.5,2.5)(5.5,4 ) ( 5.5,1)(5.5,2.5)(5.5,4 ) ( 5.5,1)(6.5,2.5)(5.5,4 )    one thing to note in these arguments is the ambivalence with which we treat coordinates . introductions to general relativity always emphasize the covariance of the equations , but practical examples often implicitly communicate the importance of particular coordinates . the problems this causes are seen at their worst when refereeing weak papers on exact solutions , where authors often refuse to accept that their solution is not new on the grounds its coordinate form is different from the known ones . the full understanding of the schwarzschild horizon was an important element in the general understanding of global properties of spacetimes which developed in the 1960s . conformal pictures drawn on the same principles as fig . 1 were also developed for the schwarzschild solution s generalizations given by @xmath27 the reissner - nordstrm solution being given by @xmath28 and the kttler solution by @xmath29 . incidentally , these provide two illustrations of a ` metatheorem ' that all named results have the wrong name : the kttler solution is commonly called schwarzschild - anti - de sitter , though as far as i know neither schwarzschild nor de sitter gave this solution , and weyl gave the general class to which the reissner - nordstrm metric belongs before nordstrm s paper appeared . other examples of this are the `` tolman - bondi '' solutions , due to lematre @xcite , whose paper tolman @xcite cited , and the `` bertotti - robinson '' solution , which robinson himself noted was first found by levi - civita @xcite . the schwarzschild solution had a pivotal role in two other developments . approximations to it were used in predictions of the classical tests of general relativity : those approximations were later generalized to the ppn formalism which provides the basis of analysis of solar system tests of relativity . it was also the setting in which hawking made the link between the laws of black hole mechanics , quantum field theory and thermodynamics in his discovery of hawking radiation . the thermodynamic identification of black hole surface area with entropy is in general not fully understood in terms of microscopic states , unlike entropy in normal statistical mechanics ( see @xcite ) . however , for the ` extreme ' charged case ( which is supersymmetric ) or near - extreme charged black holes , a microscopic basis in string theory has been given @xcite . yet again the schwarzschild solution s generalizations are playing a major role in advancing our understanding .    to summarize , the schwarzschild solution assisted the development of our understanding of the following :    * there is no `` point mass '' solution in relativity , since there is no point centre ; * the role of coordinates has to be properly understood ; * global concepts such as + black holes , + event horizons , + apparent horizons , + trapped surfaces and the singularity theorems , + cosmic censorship and naked singularities ; * ppn expansions ; and * hawking radiation , qft in curved spaces and black hole entropy . these are the solutions with the metric @xmath30.\ ] ] where @xmath31    it would be impossible to overemphasize the importance of these solutions in cosmology . they are fundamental to the inferences from observation of the presence on the cosmological scale of dark matter and `` dark energy '' ( above the densities required by dynamics of galaxy clusters ) . these matters were discussed by other speakers at the meeting so i will not give details here . there are many specific solutions due to einstein himself , de sitter , friedman , eddington , lematre , and so on ( see ch .  14 of @xcite ) . the pivotal role of the ones due to friedman and lematre led to those authors names being coupled with the names of robertson and walker in the short name flrw . the role of robertson and walker was that they independently elucidated the geometrical basis of the metric form in the 1930s . one may also note the wide use of the lematre - tolman - bondi models , spherically symmetric inhomogeneous dust models , to describe collapse and voids , primordial black holes , and other inhomogeneities in cosmology ( see @xcite ) . for example , it has been shown that the number distribution of galaxies can be modelled by a non - evolving population in an ltb spacetime rather than an evolving population in flrw . two very big steps forward were made in the early 1950s , which provided the framework for many of the new solutions found up to now . the first of these was taub s 1951 paper @xcite which contained the introduction of group - theoretic and differential geometric methods used in an essential way . ( there was also some interaction with parallel work of gdel . ) this paper also discovered the taub portion of the important taub - nut solution , and hinted at tetrad methods , though these latter were not fully developed until the 1960s , by ellis @xcite , estabrook and wahlquist @xcite , newman and penrose @xcite and others ( a development which i helped to codify @xcite as far as the orthonormal tetrad version was concerned ) . the second was petrov s classification of the weyl tensor ( @xcite ) . this , together with the solutions of the kundt class , played an important part in the development of gravitational radiation theory and formed a big step on the road to invariant classification of solutions , discussed later . one can briefly describe the petrov classification as following from the fact that for any non - zero weyl tensor @xmath32 , as defined by @xmath33}^b+2\delta _ { [ c}^{[a}r_{\;\;d]}^{b]},\ ] ] there are four null vectors @xmath34 , the principal null directions , which satisfy the equation @xmath35bc[d}k_{f]}k^bk^c=0.\ ] ] if these four vectors are distinct , one has the general case , known as petrov type i. the remaining cases , where two or more coincide , are called algebraically special . if just one pair coincide , we have type ii , if three coincide , type iii , and if all four coincide , type n : the other possible degeneracy , two coinciding pairs , is called type d. the case of zero weyl tensor , where the spacetime is conformally flat , is sometimes called type o.    these two developments led to the two main themes of the organization of the exact solutions book @xcite , i.e.  classification by symmetry groups and algebraically special metrics .    work on the latter area was much assisted by the development of the already - mentioned newman - penrose formalism @xcite , a calculation technique not based on coordinates but on a tetrad of null vectors , and therefore well adapted to the study of algebraically special spacetimes . these ideas led to the discovery of a number of important solutions , for example the following .    * the taub - nut solution ( or rather its nut part ) : @xcite . * the kerr solution , the rotating black hole : @xcite . * the robinson - trautman class : @xcite . * ( later ) the kerr - schild ansatz and solution class : @xcite .    the excellent review of ehlers and kundt @xcite , perhaps the first modern general review of exact solutions , summarized the first decade of the resulting developments . i will now pause in describing methods for finding solutions and discuss the uses of three of these solutions or solution classes . the metric in this case can be given as @xmath36 where @xmath37 , which is positive in the taub region and negative in the nut region , is given by @xmath38 this played such a role that misner described it as a `` counterexample to almost anything '' @xcite . papers by misner @xcite and misner and taub @xcite established that the taub and nut regions can be joined , that the nut region contains closed timelike lines and no sensible cauchy surfaces , that there are two inequivalent maximal analytic extensions of the taub region ( or one non - hausdorff manifold with both extensions ) , that taub - nut space is nonsingular in the sense of a curvature singularity , and that there are geodesics of finite affine parameter length .    to summarize , it has the following properties :    * the topology of group orbits changes at the horizon ; * there are closed timelike lines in the nut region ; * the boundary of the taub region has closed null geodesics ; * there is geodesic incompleteness at finite affine parameter without a curvature singularity ; and * there are inequivalent extensions , or a non - hausdorff one .    the solutions also have applications and generalizations outside strict general relativity , e.g.  in string theory or euclidean quantum gravity . the solution has thus had a great influence on studies of exact solutions and cosmological models which are spatially - homogeneous , and more generally on those which are hypersurface - homogeneous and self - similar ( see e.g.  the discussion in @xcite ) , on cosmology in general , and on our understanding of global analysis and singularities in space - times . the plane waves were first found by brinkmann @xcite but their significance , showing among other things that gravitational waves were definitely not a coordinate effect , was not appreciated until the 1950s , in work of bondi , pirani and robinson @xcite , peres @xcite , hely @xcite and others . these are the metrics @xmath40 with @xmath41 . they are members of the more general @xmath39-wave class , which are given , for electrovacuum cases , by @xmath42},\ ] ] where the functions @xmath43 and @xmath44 are arbitrary functions analytic in @xmath45 and dependent on the retarded time coordinate @xmath18 . here @xmath46 is the constant in the einstein equations .    among their properties are :    * the wave front speed is the speed of light ; * they have a transverse character ; * there are focussing effects leading to caustics ; * there is no global cauchy surface ( i.e.  no initial value problem ) . these solutions are comprehensively discussed in the review of ehlers and kundt @xcite , so the treatment in the exact solutions book @xcite adds virtually nothing . they have interesting singularity and horizon behaviour , and gave rise to a number of special cases of interest , such as sandwich waves and impulsive waves . one can also study collisions between them @xcite and the resulting metrics have been studied in the monograph of griffiths @xcite : see also ch . 25 of @xcite . they have more recently been used in understanding higher - dimensional theories and as examples in quantum field theory ( they have the property of having no quantum corrections ) . together with the robertson - walker and schwarzschild solutions , this is probably the best - known exact solution , because it represents the unique rotating vacuum black hole . it has a number of interesting mathematical properties , having indeed been found as the complement to the nut investigation of petrov type d vacuum solutions ( the nut solutions had initially been thought to be the only such solutions ) . the metric is @xmath47   \\ &   & -\left ( 1-\dfrac{2mr}{r^2+a^2\cos ^2\vartheta } \right ) \bigg ( { { \mbox{\rm d}}}t+\dfrac{2mar\sin ^2\vartheta { { \mbox{\rm d}}}\varphi } { r^2 - 2mr+a^2\cos ^2\vartheta } \bigg)^2 . \nonumber\end{aligned}\ ] ]    when the global structure of this metric was worked out , it emerged that it has a ring singularity ( for @xmath48 ) through which further exterior regions can be connected : see e.g.  the summary in @xcite . moreover , the hamilton - jacobi and klein - gordon equations are separable in this metric , which is related to the fact that it has a non - trivial killing tensor . it exhibits the phenomenon of an ergosphere , a region outside the black hole horizon but within which any particle has to corotate around the hole . there is a relation to work on characterization of stationary axisymmetric spacetimes by multipole moments : the kerr solution has very specific relations between its moments which do not appear to be found in physical bodies of rotating fluid , posing a question about possible sources or the process of approach to the kerr solution as the eventual black hole outcome of a collapse . work on its mathematical properties , however , is likely to be exceeded by the many papers on its astrophysical implications . for example , the ergosphere classically allows the penrose process , in which a part of a body dividing within the ergosphere can emerge with more energy than the original body entered with . the wave version of this is superradiance , where the scattered wave has more energy than the incident wave , and this phenomenon was one of the stimuli to the laws of black hole mechanics later explained by hawking . it is also related to the blandford - znajek mechanism @xcite , in which a magnetic field threading the black hole can extract rotational energy from the hole . the most important of all the astrophysical uses of the kerr solution is probably as a basis for accretion disk physics , thought , for example , to be responsible for the x - ray emission of x - ray binaries in the sky , and used in explanations of larger objects such as jets in active galactic nuclei . observations of astronomical accretion disks are now suggesting the objects at their centres really are kerr black holes , as the only way to explain both their short periods and other orbital data @xcite . chandrasekhar remarked `` in my entire scientific life , extending over forty - five years , the most shattering experience has been the realization that an exact solution of einstein s equations , discovered by the new zealand mathematician roy kerr , provides an absolutely exact representation of untold numbers of massive black holes that populate the universe ... ''      a third set of novel techniques appearing for the first time in papers in the 50s and early 60s ( by buchdahl , ehlers , and bonnor , for instance @xcite ) took rather longer to grow to maturity than the first two methods described in this section : this third topic is that of the generating techniques . although they exist for metrics with one killing vector , they are most used for stationary axisymmetric solutions , and the other classes with two commuting killing vectors : cylindrical waves and colliding plane waves , boost - rotation symmetric spacetimes and cosmologies with two commuting spacelike killing vectors . they work not only for vacuum , but for other forms of matter with characteristic propagation speed equal to the speed of light : massless scalar fields ( or stiff fluid in the case of a timelike gradient of the field ) , ( massless ) neutrinos , and electromagnetism . this area exploded in the 1970s and 1980s in work of many people . the number of methods proliferated , and created a secondary industry of understanding the relations between them . at the same time the number of applications grew hugely . there is so much work that i shall not attempt to give references here but instead refer the reader to @xcite ( see also @xcite and @xcite ) . the best known formulation for the basic equations is that due to ernst . for the stationary axisymmetric metrics @xmath49 where the metric functions @xmath37 , @xmath50 , @xmath51 , @xmath52 depend only on the coordinates @xmath53 which label the points on the 2-surfaces @xmath54 orthogonal to the orbits , and the electromagnetic field is given by a complex potential @xmath55 , it reads @xmath56 to recover the metric one needs the equations @xmath57 where @xmath58 is the levi - civita tensor in the @xmath54 . although i will not attempt a review here of the methods , the works of geroch , neugebauer and kramer , hoenselaers , kinnersley and xanthopoulos ( hkx ) , kinnersley and chitre , harrison , belinski and zakharov , hauser and ernst , yamazaki , and cosgrove were so influential their names should be mentioned . these methods themselves can now be seen as embedded in an even more general context of symmetries of differential equations and integrable systems , involving concepts such as inverse scattering and lax pairs , bcklund transformations , riemann - hilbert problems , prolongation and so on ( see @xcite or ch . 10 of @xcite ) . thus the solutions have offshoots in mathematics and in other physical theories . they provide a unification of results on known solutions ( for example , all stationary axisymmetric electrovacuum solutions in which a portion of the rotation axis is regular can be generated from flat space ) . they also enable an infinite number of solutions to be obtained . at one stage there were a significant number of papers which exploited this possibility by exhibiting specific solutions , but it quickly became apparent that merely obtaining a new solution , when it can be done infinitely often ( at least in principle ) , is pointless . attention is now normally directed to ways to generate solutions with predetermined characteristics : for example one can ask what class of axis data gives a certain feature to the solution ?    of the solutions obtained or obtainable by these methods , a number have been of importance or interest , for example the tomimatsu - sato family , the double kerr solution , the neugebauer - meinel dust disk , and a number of colliding wave and cosmological solutions . brevity precludes detailed discussion and the reader is referred to the literature already cited . i now return to the fundamental question raised earlier : how can we compare solutions ? in particular how can we test for local isometric equivalence ? this is called the equivalence problem , and belongs to the general class of recognition problems . it has largely been answered , in theory and in practice , using ideas due to cartan , brans , karlhede and others , with practical implementation and development by man , and later myself and my group , although in a formal sense its final step is undecideable . i will briefly describe the procedure now : a full description would be a lecture in itself . for a review see ch . 9 of @xcite . the key point is that scalar polynomial invariants , i.e.  polynomials in the riemann tensor and its derivatives in which all indices are contracted over , are insufficient . instead we need curvature invariants of a more general sort , the ` cartan invariants ' . these are , for example , given by the components of the weyl tensor referred to a tetrad chosen using the principal null directions . these enable local characterization of the metric . the basic idea is to use invariantly defined tetrads , like the one from the principal null directions , and take components of the riemann tensor and its derivatives in this frame . thus the method has links with the petrov classification and tetrad methods which were introduced in the 1950s and 1960s . counting functionally independent invariants gives the dimension of the symmetry group , thus linking to the group theoretic ideas of that period , and additional information available can give the group structure . these characterization methods can be used to check if solutions found are really new , or to search among known solutions for examples with desired local properties . they could in principle be used to find solutions , as well as classify known ones , and first examples of this method have been developed by bradley , karlhede and marklund . other uses of this approach , i.e. the direct use of invariants , have been in understanding the limits of families of solutions without needing trial and error for the appropriate coordinate transformations ( for example , studying the limits of the schwarzschild family as @xmath59 ) ; proving ( non ) existence of matchings by characterizing the geometries of the proposed matching surfaces ( e.g.   in work of my student daniel cox @xcite ) ; and in providing a method to ` unravel ' directional singularities ( in work of my student john taylor @xcite ) . the classifying quantities might also be used to give a topology on the space of solutions , which could even be of interest in numerical relativity . summarizing , the main methods for finding solutions in current use are still those outlined in section 4 , but the understanding and classification of these solutions can now be done in an invariant manner which should enable better use of the solutions found and may also provide a fresh and more invariant way to find more . the second edition of the exact solutions book contains about 400 pages of new material , covering hundreds of new solutions and references . in its preparation the authors read about 4000 new papers ( as well as the 3000 read for the first edition ) . so there is far more material than i could talk about . most of the new material is integrated into and expands existing chapters and sections : additional sections were added on the ghp formalism and other calculi , junction conditions and so on , and the chapter on solutions obtained by generating techniques was almost completely re - written .      * homotheties * characterization by invariants * generating techniques themselves * dynamical systems methods * inhomogeneous solutions with two spacelike killing vectors * colliding plane waves * special vector and tensor fields . j.  bicak , `` the role of exact solutions of einstein s equations in the developments of general relativity and astrophysics : selected themes , '' in _ einstein s field equations and their physical implications _ , springer verlag , heidelberg , 2000 , vol . 540 of _ lecture notes in physics_.              j.  ehlers , and w.  kundt , `` exact solutions of the gravitational field equations , '' in _ gravitation : an introduction to current research _ , edited by l.  witten , wiley , new york and london , 1962 , pp . 49101 . m.  maccallum , `` an overview of exact solutions of einstein s equations and their classification , '' in _ highlights in gravitation and cosmology ( proceedings of the goa conference 1987 ) _ , edited by b.  iyer , a.  kembhavi , j.  narlikar , and c.  vishveshwara , cambridge university press , london and new york , 1989 , pp . 314 . a.  petrov , _ scientific proceedings of kazan state university ( named after v.i . ulyanov - lenin ) , jubilee ( 1804 - 1954 ) collection _ * 114 * , 5569 ( 1954 ) , translation by j. jezierski and m.a.h . maccallum , with introduction by m.a.h . maccallum , gen . . grav . * 32 * , 1661 - 1685 ( 2000 ) . r.  kerr , and a.  schild , `` a new class of vacuum solutions of the einstein field equations , '' in _ atti del convegno sulla relativit generale ; problemi dellenergia e ondi gravitationali . _ , barbra , firenze , 1965 , p. 222 . c.  misner , `` taub - nut space as a counterexample to almost anything , '' in _ relativity theory and astrophysics , vol . 1 : relativity and cosmology _ , edited by j.  ehlers , lectures in applied mathematics , volume 8 , american mathematical society , providence , r.i . , 1963 ,
Abstract:

Answer: the evolution of the methods used to find solutions of einstein s field equations during the last 100 years is described . early papers used assumptions on the coordinate forms of the metrics . since the 1950s more invariant methods have been deployed in most new papers . the uses to which the solutions found have been put are discussed , and it is shown that they have played an important role in the development of many aspects , both mathematical and physical , of general relativity . address = school of mathematical sciences , queen mary , university of london , mile end road , london e1 4ns , u.k . + email : m.a.h.maccallum@qmul.ac.uk