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: multilayers built from ferromagnetic layers with various spacer layers include a wide variety of magnetic film systems that have been intensively studied during last years . due to remarkable phenomena as giant magnetoresistance , exchange - spring behaviour and/or exchange bias , and surface enhanced magnetic anisotropy , such nanostructures have already found a number of applications and are considered as promising candidates for nonvolatile magnetic recording media.@xcite on the other hand , nanoscale superlattices and similar structures provide convenient model systems to study different aspects of surface magnetism and magnetic ordering in confining geometry .    in particular , much attention has been given to multilayers composed of antiferromagnetically coupled ferromagnetic nanolayers . @xcite such layered synthetic antiferromagnets can be separated into two classes according to the symmetries ruling their magnetic properties : superlattices with relatively strong uniaxial anisotropies include low - symmetry multilayers with effective uniaxial magnetic anisotropy in the layer planes , e.g. epitaxial systems deposited on ( 110 ) , ( 211 ) faces of cubic substrates.@xcite also multilayer systems with perpendicular anisotropy belong to this class.@xcite magnetization processes in these nanostructures are strongly influenced by the uniaxial anisotropy which is responsible for specific phenomena such as `` surface spin - flop '' @xcite or field - induced metamagnetic jumps.@xcite on the other hand , the uniaxial anisotropy may be absent in layered systems with higher symmetry . these represent another large and intensively investigated class of synthetic antiferromagnetic nanostructures . superlattices with planar magnetization grown on ( 001 ) faces of cubic substrates , e.g. multilayers from materials combinations as co / cu @xcite , fe / si @xcite co / cr , @xcite or fe / cr , @xcite belong to this class . in the case of weak four - fold anisotropy , their magnetic properties are mostly determined by the interlayer exchange interactions which may include an important _ biquadratic _ contribution . @xcite evidence of strong biquadratic exchange interaction has been given in a number of experimental papers for fe / cr two - layers and multilayers @xcite , for fe / si(001 ) multilayers @xcite ) . strong effective four - fold anisotropies have been found in systems such as co / cr(001 ) or fe / cr(001 ) . @xcite the competition between magnetic anisotropy , applied fields and exchange energies may cause complicated magnetic effects and processes . in fact , a great number of novel magnetic configurations and remarkable reorientational effects in external fields have been found in such superlattices . @xcite in particular , recent experimental results using modern depth - resolving techniques reveal spatially inhomogeneous magnetic structures , e.g. in fe / cr(001 ) superlattices,@xcite and specific reorientation effects imposed by four - fold planar ( tetragonal ) anisotropy.@xcite the understanding and interpretation of the complex magnetization processes found in such systems @xcite requires a theoretical underpinning . theoretical activity in this field is largely based on analytical and numerical calculations mostly within phenomenological approaches . @xcite these studies have demonstrated the general validity of the phenomenological models to describe the magnetization processes in antiferromagnetic nanostructures . @xcite for the system under discussion the phenomenological theory has been developed to describe effects of 90  degree couplings and a concomitant complex evolution of domain structures in fe / cr / fe layers.@xcite four - fold anisotropy effects have been theoretically investigated in ref .   for sandwich structures with @xmath2=2 ferromagnetic layers coupled through a spacer . multilayer systems provide also experimental models to study effects of the confining surfaces on antiferromagnetic structures.@xcite in this context , specific inhomogeneous magnetic states described for theoretical models in ref .  , have recently been observed in fe / cr superlattices.@xcite the existing theoretical results , however , are restricted to special cases and are not sufficient for an exhaustive description of the magnetic states and field - induced reorientional effects observed in recent experiment . @xcite in this paper we provide a theoretical analysis of magnetic states and magnetization processes in planar antiferromagnetic superlattices with and without four - fold anisotropy in magnetic fields applied within the multilayer plane . in such magnetic superlattices the exchange interlayer coupling is an oscillatory function of the spacer thickness.@xcite depending on the spacer thickness an alternating sequence of ferromagnetic and antiferromagnetic interlayer couplings is realized , and , by adjusting the spacer thicknesses , very different strengths of antiferromagnetic coupling can be realized . on the other hand , the four - fold anisotropy includes interface - induced contributions which implies a strong dependency of the effective average anisotropy of each ferromagnetic layer on the layer thickness.@xcite thus , the effective magnetic interactions can vary in extremely broad ranges for such multilayers in dependence on the chosen materials combinations and thicknesses , see e.g. ref .  , where for a co / cu(001 ) two - layer system with a wedged spacer layer the ratio between four - fold anisotropy energy and the exchange coupling is changed by orders of magnitude . in contrast to bulk planar antiferromagnets , where an essentially fixed hierarchy for the strengths of the magnetic interactions holds,@xcite in these artificial antiferromagnetic systems the ratios between different magnetic energies , respectively the phenomenological parameters in the magnetic free - energy , may assume practically arbitrary values . moreover , as the interlayer exchange is weak compared to direct exchange interactions , the fields to induce spin - reorientation phenomena are similarly weak and experimentally accessible . the rich phase diagrams for these systems precludes an analysis in all details . the phase space in terms of the phenomenological parameters includes a large variety of different magnetic states with a corresponding multitude of spontaneous and field - induced phase - transformations . in a first step to such an analysis , all laterally homogeneous states in such multilayers must be found . they are the building blocks for a domain theory.@xcite for the case of laterally homogeneous states of each ferromagnetic sublayer one has a system that behaves like an antiferromagnetically coupled chain of stoner - wohlfarth particles . this simplified one - dimensional model for the behaviour across the multilayer stack also yields the limiting case for the magnetization processes with maximum hystereses . again , a direct analysis of all magnetic states even for these one - dimensional models yields an intricate succession of phase diagrams and magnetization curves.@xcite in this paper , we avoid the cumbersome task of listing and classifying _ all _ solutions and transitions . instead , we provide a broad physically intuitive picture of the physical effects due to the different exchange or anisotropic forces , and those imposed by the confining geometry of the system . to this end , we study limiting cases of the model . this includes the case of strictly zero anisotropy with and without biquadratic exchange , and the case of infinite anisotropy with fixed four - fold orientation of the magnetizations in each layer . for the antiferromagnetic two - layer systems ( @xmath2=2 ) , i.e. the sandwich structures ( experimentally realized as ferromagnetic / spacer / ferromagnetic trilayers ) , we provide a detailed investigation of the magnetic phase - diagram for arbitrary orientation of fields in the layer planes . based on this , we can understand the _ basic _ magnetic configurations in the multilayers , and we can give a map of the _ topologically different types _ of magnetic phase diagrams . we use standard methods to analyse magnetic phases and transitions within the phenomenological approach and the theory of phase transitions.@xcite the one - dimensional chain models are considered as composite order - parameters with many components ( @xmath2 components in a multilayer stack composed of @xmath2 ferromagnetic / spacer bilayers ) and a characteristic structure of couplings between the components defined by interlayer exchange and the surfaces . from this point of view , the very rich phase - diagrams and correspondingly complex sequences of magnetic configurations can be understood . for the general cases of the model , the equations for equilibrium and phase stability can be solved only by numerical methods . with the methods and results expounded below , one can extend the analysis to specific experimental cases in all detail . the solutions for the one - dimensional chain models include various field - induced canted and inhomogeneous states with a net magnetization . based on the phase diagrams of these models , the evolution of laterally inhomogeneous ( domain ) states and magnetization processes can be discussed.@xcite in this connection , the coexistence regions of different phases in the vicinity of discontinuous ( first - order ) magnetic phase transitions are important . in external fields , thermodynamically stable domain - configuration from these competing phases can be established in extended multilayers this is crucial for an understanding of the hysteretic magnetization processes under coercivity mechanisms . the paper is organized as follows . after introducing the model and mathematical tools ( sec . ii ) we consider the effects of the bilinear and biquadratic exchange interactions in the next section ( sec . iii ) . in sec . iv we investigate in detail four - fold anisotropy effects in antiferromagnetically coupled two - layers and then discuss the generalization of these finding to the case of multilayers . in sec . iv we discuss domain states and magnetization processes by using qualitative arguments . in the concluding part we discuss possible extensions of the theory , and we suggest some useful experiments to enhance our understanding of magnetization processes in antiferromagnetic superlattices . let us consider a stack of @xmath2 ferromagnetic plates infinite in @xmath3- and @xmath4-directions and with finite thickness along @xmath5-axis . the magnetization of the layers is @xmath6 , and there are indirect interlayer - exchange couplings through spaces between them . the phenomenological energy of the system can be written in the following form    @xmath7 -\mathbf{h}\cdot \sum_{i=1}^{n}\zeta_i \mathbf{m}_i - \frac{1}{2}\sum_{i=1}^{n } k_i\left ( \mathbf{m}_x^4   + \mathbf{m}_y^4 \right )   \label{energy1}\end{aligned}\ ] ]    where @xmath8 ( @xmath9 ) are unity vectors along the @xmath10-th layer magnetization . @xmath11 designate deviations of the magnetization in the @xmath10-th layer from the averaged value @xmath12 . we assume that @xmath13=0 , i.e. the layer magnetizations are restricted to the layer plane . @xmath14 and @xmath15 are constants of bilinear and biquadratic exchange interactions , respectively ; @xmath16 are constants of the in - plane four - fold anisotropy . the functional ( [ energy1 ] ) generalizes similar models considered earlier in a number of studies on exchange @xcite and anisotropy @xcite effects in planar antiferromagnetic systems . within this approach the ferromagnetic layers are considered as homogeneously magnetized blocks with constant values of the magnetic interactions . this assumption deserves some comment . it is well - established that in magnetic nanostructures surface / interface exchange and relativistic interactions strongly modify electronic and magnetic properties within all volume of the magnetic constituents @xcite . this means that the values of the exchange or anisotropy parameters , and the magnetizations include large interface / surface - induced components which may strongly vary across the thickness of the individual layers.@xcite however , the hypothesis of magnetic homogeneity in the ferromagnetic nanolayers in the models of type ( [ energy1 ] ) has a solid physical basis and is justified by successful applications of these models to describe magnetization processes in layered ferro- and antiferromagnetic nanostructures.@xcite this relies on the fact that in ferromagnetic nanolayers the intrinsic ( direct ) exchange coupling are usually very strong and overcome surface / interface induced interactions . thus , they play the dominating role to determine the magnetic order _ within _ the layers , which reacts also very stiffly on all external and induced magnetic forces . furthermore , in these planar systems the stray field forces confine the magnetization of the layers into their plane . as a result , in most of these systems the ferromagnetic layers preserve essentially homogeneous in - plane magnetized states even under the influence of strongly inhomogeneous induced interactions.@xcite in this connection it is important to stress here that the materials parameters in the phenomenological model ( [ energy1 ] ) should be treated as averages over a multilayer _ period_. they comprise ( in integral form ) all intrinsic and induced energy contributions of the magnetic states in the ferromagnetic layers . in particular , for the systems under consideration they may also include a contribution from magnetism of the spacer layers . thus , in contrast to their bulk counterparts which are considered as constants of the magnetic material , in nanoscale multilayer systems the phenomenological parameters strongly depend on many physical and geometrical factors and may considerably vary from sample to sample . parity of @xmath2 also plays an important role . the model ( [ energy1 ] ) with _ even _ number of ferromagnetic layers represents systems with fully compensated magnetization . such superlattices can be treated as analogues to bulk collinear antiferromagnets . superlattices with odd numbers of layers or with unequal thicknesses of layers own a noncompensated magnetization which strongly determines their global magnetic properties . in their response to an applied field , these structures are similar to bulk ferrimagnets . they could be studied by similar methods as used below , but have to be considered as separate class of systems . concentrating on the properties of systems with fully compensated total magnetization , we consider only superlattices with even @xmath2 and equal magnetization in all layers ( @xmath17 ) . in this contribution dedicated to general properties of these antiferromagnetic superlattices , we avoid secondary effects which are related to their micro - structure and interface properties , such as strains , chemical intermixing etc . with that in view we will study the model ( [ energy1 ] ) mainly for the case of _ identical _ ferromagnetic layers assuming that the exchange and anisotropy constants are equal for all the layers , @xmath18 , @xmath19 , and @xmath20 . some of the analytical results can be generalized to the isotropic ( i. e. with @xmath21 ) system with mirror symmetry about the center of the layer stack , i.e.  @xmath22 and @xmath23 . in our problem the magnetization of the _ i_th layer is confined in planar orientations and can be described by the angle @xmath24 between the vector @xmath25 and @xmath3-axis . thus , calculations of the magnetic states for the model ( [ energy1 ] ) can be reduced to the optimization of the function @xmath26 . following the theory of bulk antiferromagnetism it is convenient to introduce here the vectors of the _ total magnetization _ @xmath27 and the _ staggered _ ( or _ antiferromagnetism _ ) vector @xmath28 @xmath29 the energy ( [ energy1 ] ) is invariant under transformation @xmath30 , and , therefore , all magnetic states in this model are degenerate with respect to the sign of the staggered vectors . in the following only solutions with a definite sign of the staggered vector will be discussed . however , one should keep in mind that the full set of solutions includes also those with opposite sign of @xmath28 . the magnetic states with opposite sign of @xmath28 behave identically in an applied magnetic field . thus , for the magnetization processes there is no need to distinguish these solutions .    for the particular case of a two - layer ( @xmath0 ) the energy ( [ energy1 ] ) can be transformed to the following form @xmath31 where @xmath32 , @xmath33 and @xmath34 is the angle between the @xmath35 and the @xmath3-axis . this case of a two - layer is of special interest . the energy ( [ energy2 ] ) functionally coincides with that of a bulk two - sublattice antiferromagnet . this offers the opportunity for useful physical relations and analogies with bulk antiferromagnetism . on the other hand , the two - layers represent the simplest model of antiferromagnetically coupled layered nanostructures . the equations minimizing the function @xmath26 have strongly nonlinear character and no analytical solution are available generally . the body of our results have been obtained by numerical methods . we solve the coupled equations for equilibria @xmath36 by an efficient conjugate gradient minimization.@xcite for systems with large anisotropy the configurations with the magnetizations oriented in the four - fold easy - axis directions are close to mostly metastable minima of the energy . these configurations have been systematically searched for global absolute stability and checked for stability . numerically , this is feasible for @xmath2 not too large with present day computers . stability limits and phase - transitions are determined from the evolution of the smallest eigenvalue @xmath37 of the stability matrix @xmath38 under changing parameters in ( [ energy1 ] ) , i.e. ratios @xmath39 , @xmath40 , or the applied magnetic field . basic spin - configurations in isotropic antiferromagnetically coupled two - layers ( @xmath2=2 ) . at zero field ( a ) the magnetization vectors @xmath25 are antiparallel in adjacent layers but the structure has no fixed orientation in the plane ( infinite degeneracy ) . an in - plane magnetic field lifts the degeneracy and stabilizes states with * l * perpendicular to the field and with * m * along the field - this is a spin - flop ( sf ) phase ( b ) . in a sufficiently strong magnetic field the sf phases transform into the saturated ( _ flip _ ) phase ( c),width=321 ] we start the investigations of the static configurations minimizing energy ( [ energy1 ] ) from the isotropic case ( i.e. with @xmath16 = 0 ) . depending on relations between bilinear and biquadratic constants different collinear and noncollinear configurations are stable.@xcite ) the equations @xmath41 , @xmath42 , @xmath43 determine the region in the parameter space where a collinear _ antiferromagnetic _ ( af ) phase is the zero - field ground - state consisting of blocks of adjacent layers with antiparallel orientation of the magnetization ( fig . we restrict our analysis to this practically important case . the af phase is infinitely degenerate with respect to orientations of the staggered vector * l * in planar directions . the magnetic states for the isotropic two - layer ( @xmath0 ) in the applied field can be obtained by minimization of the energy ( [ energy2 ] ) with @xmath44 . the solutions ( @xmath45 , @xmath46 ) describe the states with the staggered vector perpendicular to the applied field ( so called _ spin - flop _ ( sf ) phase ) . in an increasing magnetic field the magnetization vectors rotate into the field direction ( fig . 1 ( b ) ) , and finally a continuous transition into the ferromagnetic ( _ flip _ ) phase with @xmath47 . this ( _ spin - flip _ transition ) occurs at the exchange - field @xmath48 ( fig . 1 ( c ) ) , which is lower than the corresponding spin - flip field @xmath49 for a bulk system . in superlattices with @xmath1 the magnetic field applied in the plane also lifts the degeneracy and stabilizes the sf phase ( the state with @xmath50 and @xmath51 ) ( figs . however , the magnetic configurations in these sf states and their evolution in the applied field are markedly different from those in bulk antiferromagnets or in the two - layer systems . it turns out that magnetic structures of the multilayers with @xmath2 divisible by four ( @xmath52 called here _ even - even _ ) differ from _ even - odd _ systems with @xmath53 ( @xmath54 ) . in low fields ( @xmath55 ) the solutions for the sf phase are given by a set of linear equations @xmath56 where @xmath57 for even - even and @xmath58 even - odd systems . it is clear that for the isotropic model the direction of the magnetic field in the film plane plays no role . following the definition of the angles @xmath24 we assume here that the magnetic field is applied along the @xmath3-axis . the solutions ( [ sol1 ] ) describe small deviations of the magnetization vectors , from the directions perpendicular to the easy axis ( see the magnetization profiles in figs .  2 , 4 and the configuration in fig . 3 ( b ) ) . in all internal layers ( @xmath59 ) neighbouring pairs retain antiparallel orientations ( e.g. , ( @xmath60 ) and ( @xmath61 ) in fig . 3b ) . for even - even systems the magnetization vectors of the central layers ( @xmath62 , @xmath63 ) preserve perpendicular orientation even in finite ( but weak ) fields ( see layers 6 and 7 in fig . 3 ) . for even - odd systems the magnetization vectors of all layers have finite deviations from perpendicular orientation . towards top and bottom layer @xmath64 or @xmath2 in the stack , these deviations increase . due to antiparallel orientation the pairs have zero net magnetization . the total magnetization arises only due to the top and bottom layers @xmath64 and @xmath2 . note that in these configurations a number of layers rotates against the applied field ( in fig . 2 - 4 the layers with @xmath10 = 2 , 4 and their symmetric counterparts with @xmath10 = 11 , 9 ) . this occurs because , in weak fields , the exchange interactions favouring antiparallel magnetizations in adjacent layers play the dominant role . an increasing magnetic field counteracts and slows down this reverse rotation , when the non - linear evolution for the magnetization structure is reached , and finally the sense of rotation is changed at characteristic fields @xmath65 where @xmath66 . with increasing @xmath2 the number of layers increases which display this reverse rotation , the deviations from orientations perpendicular to the applied low field become larger near the surface layers , and the fields @xmath65 are reached at lower fields ( see inset in fig . 2 for @xmath67 ) . a further set of characteristic fields @xmath68 defines the points where the projection of @xmath69 on the field directions changes the sign ( the field @xmath70 and @xmath71 are indicated in fig . 2 . ) in increasing field these characteristic fields initially are reached for the central layers and at higher fields for those closer to the boundaries . for fields @xmath72 the magnetization directions of all layers have positive projections onto the field direction . in the model with equal exchange constant there is another special field @xmath73 ( independent of @xmath2 ) , where all inner layers have the same projection onto the field direction ( @xmath74 , @xmath75 ) . the parameters of this `` knot '' point are determined from a set of equations @xmath76 the functions @xmath77 , @xmath78 , @xmath79 are plotted in fig . [ knotpoint ] . in particular for superlattices with zero biquadratic exchange @xmath80 one has @xmath81 and @xmath82 . near saturation , the sf phase has only positive projections of the magnetization on the direction of the magnetic field which decreases towards the center similarly to spin configurations obtained numerically in ref .  . evolution of spin - states in multilayers with a field @xmath83 applied in the plane of the layers and for zero anisotropy @xmath84 and biquadratic exchange @xmath80 , hence @xmath85 . main figure : rotation angles of the layers in the plane @xmath24 @xmath86 for superlattice with @xmath87 . at point labeled `` r '' the sense of rotation changes for the magnetization direction @xmath24 ( here , @xmath88 ) ; at point `` @xmath89 '' the magnetization of this layer is perpendicular to the applied field again . point `` @xmath90 '' is the `` knot - point '' ( eq . ( [ knot ] ) ) . inset : @xmath67 only @xmath24 for even @xmath88 and @xmath91 shown . , width=321 ]     spin - states in multilayers with @xmath87 for zero anisotropy @xmath84 and field @xmath83 applied in the plane : ( a ) @xmath92 ( b ) @xmath93 ( c ) @xmath94 ( d ) @xmath95 ( e ) @xmath96 , and ( f ) @xmath97 , i.e. below the transition to the flip - phase ( cf . [ evolutionkzero ] ) . all states have mirror - symmetry about the center of the multilayer stack ( marked by the dashed line ) . , width=321 ]    the spatial inhomogeneity of the sf phases in the multilayers and their remarkable evolution in the applied field are due to the particular _ finite - size _ effects in this type of magnetic nanostructures . in such layered nanostructures all internal layers interact with two adjacent layers but the top and bottom layer have only one neighbouring layer . thus , their exchange coupling is weakened compared to internal layers within the multilayer stack . this strong disbalance of the exchange forces at the boundaries affects the magnetic ordering within the entire multilayer and causes the reorientational processes in the sf phases . this mechanism is also responsible for the occurrence of inhomogeneous states in other noncollinear structures in the multilayers ( see iv.c and v.a ) . note , this important role of the cut exchange - couplings is specific to antiferromagnetic confined systems with noncompensated moments at surfaces . in contrast , the cutting of exchange couplings does not affect _ ferromagnetically _ coupled nanostructures where relative orientation of the magnetization in the layers essentially depend only on magnetic anisotropy , applied fields , and the demagnetizing stray - fields . evolution of spin - flop profiles for the multilayer with finite biquadratic exchange is plotted in fig . 4 . for the systems with collinear ground - states the biquadratic exchange does not induce reorientional transitions but it has rather strong influence on the distribution of the magnetization in the multilayers ( see insets in fig .  4 ) and the value of the characteristic fields ( fig .  5 ) . the effects of a biquadratic exchange @xmath98 , if present in multilayers with collinear antiferromagnetic ground - state , can be easily understood . for magnetic configurations close to the antiparallel orientation of neighbouring layers , it softens the linear exchange forces ; for nearly parallel orientation the system becomes stiffer instead . the evolution of the magnetic states in the external field is distorted accordingly by the presence of the biquadratic exchange @xmath98 ( fig . [ evolutionkappa ] ) . in low fields close to the antiferromagnetic ground - state , the system reacts more strongly on the external field . thus , the fields @xmath99 , where the rotation of the ithe layer reverses , and the field @xmath68 , where perpendicular orientation with respect to field is regained , respectively , are reached at relatively lower fields ( labels `` r '' and `` @xmath89 '' in fig .  [ evolutionkappa ] ) . conversely , the ferromagnetic state is reached only close to the enhanced exchange fields @xmath100 reflecting the stiffening of the effective exchange couplings in nearly ferromagnetic configurations . thus , the magnetization curves acquire a nonlinear character with increasing @xmath98 ( inset ( a ) in fig . [ evolutionkappa ] ) . experimental observation of the fields @xmath99 and the special field @xmath101 ( fig . [ knotpoint ] ) could be used to determine the relative strength of @xmath98 . for large values of @xmath102 the reversal of the rotation of the magnetizations of individual layers takes place at very low fields ( see , inset ( b ) in fig . [ evolutionkappa ] ) . evolution of magnetic states in antiferromagnetic multilayer @xmath87 with zero anisotropy @xmath84 and biquadratic exchange @xmath103 , i.e. with enhanced @xmath49 ( compare with fig . [ evolutionkzero ] for @xmath104 ) . @xmath83 is applied in the layer plane . insets ( a ) : magnetization curves in the range of values of @xmath105 with antiferromagnetic ground - state in zero field . ( b ) : evolution of states ( only shown for layers with even numbers ) in small fields for @xmath106 . , width=321 ]     dependence of `` knot '' point @xmath107 ( eq . ( [ knot ] ) on strength of biquadratic exchange @xmath105 . inset shows the two angles @xmath108 , @xmath109 which characterise this special configuration independent on @xmath2 . , width=321 ]       spin - flip fields for antiferromagnetic multilayers ( zero anisotropy @xmath84 ) in dependence on number of layers @xmath2 . , width=321 ]    near the spin - flip transition from the sf - phase to the ferromagnetic states , the deviations of @xmath25 from the field directions are small ( @xmath110 ) , and the energy of the system ( [ energy1 ] ) can be expressed by the quadratic form @xmath111 , where the matrix @xmath112 has a tridiagonal band structure with nonzero elements only in the main diagonal @xmath113 and with the side diagonal elements @xmath114 . the spin - flip field @xmath115 is determined as the largest root of the equation @xmath116 . in particular for the model with equal parameters these solutions are plotted in fig . the spin - flip field gradually increases from @xmath117 for the two - layer to the `` bulk '' value @xmath118 in the limit of infinite @xmath2 ( fig . this dependence reflects the increase of the average number of exchange bonds from the value 1 for the two - layer and approaching 2 as @xmath2 tends to infinity . effects of the four - fold magnetic anisotropy on the states in the antiferromagnetic superlattices are revealed in detail by applying fields in arbitrary directions . we approach the general case by a detailed investigation on the case of the two - layer . to set the stage , we discuss the highly symmetric phase - diagrams with fields in direction of easy and hard - axes in the plane . this completes earlier work by dieny et al . @xcite on the @xmath2=2 systems . then , we present phase diagrams for arbitrary field directions in the layer plane . as we have seen in the last section , the primary effect of biquadratic exchange is a distortion of phase diagrams in various regions . thus , here we consider only models with @xmath119 to avoid such quantitative effects which will not affect the general topological features of the phase diagrams . optimization of the function @xmath120 eq . ( [ energy2 ] ) yields solutions for the magnetic states and their stability limits . in zero field the tetragonal anisotropy lifts the rotational degeneracy of the af phase and stabilizes two different states with perpendicular orientations of @xmath28 ( fig . 7 ( a ) ) . magnetic configurations in planar antiferromagnetically coupled two - layers with four - fold anisotropy . easy directions are indicated by thin continuous lines , hard directions by dashed lines . in zero field two antiferromagnetic configurations ( i ) ( a ) and ( ii ) ( b ) with the staggered vectors along the easy directions correspond to the ground - state . a magnetic field along the hard directions distorts them into the canted phases ( c ) , ( d ) which transform into the spin - flop state ( e ) . corresponding configurations of net magnetization and staggered vectors are shown for the two canted phases ( c ) ( d ) , and for the spin - flop phase ( e ) . , width=287 ]     magnetic phase diagram of the two - layer system in magnetic fields along _ hard _ directions . the tricritical point ( @xmath121 , @xmath122 ) separates the continuous ( dotted line ) and discontinuous ( thick line ) transitions between the canted and sf phases . a further dotted line gives the critical field of the second - order transition between the sf and flip phases @xmath123 ( [ tran2 ] ) . thin lines give the lower stability limit @xmath124 of the sf phases and the upper limit @xmath125 of the canted phases . the grey area is the region of the phase equilibrium between the two canted phases ( see fig . [ tetragonconfigs ] ( c),(c ) and ( d),(d ) ) . magnetization curves for low and high anisotropy , as indicated in the main panel , are plotted in ( a ) and ( b ) . , width=287 ]      in a field applied along one of the hard axes , the af configurations of the ground - state ( fig . 7(a),(b ) ) are distorted into low symmetry configurations ( so called _ canted _ phases - fig . 7(c),(d ) ) . both these magnetic configurations preserve mirror symmetry with respect to the field direction and have the same energy . the total magnetizations * m * of these canted states deviate from the field directions to different sides at equal angles ( fig . 7(c),(d ) ) . an oblique magnetic field deviating from the hard direction to one or the other side violates the phase balance between these two canted states . it favours the canted phase with the larger projection of * m * onto the field direction . this is the typical situation of a first - order ( or discontinuous ) transition .    in an increasing field the vectors * m * in both canted phases rotate into the field direction and reach it at a certain critical field @xmath126 . at this point both phases transform into the common configuration corresponding to the sf phase ( fig . 7(d ) ) , i.e. a second - order ( or continuous ) phase transition from low symmetry canted phases into the high - symmetry sf - phase occurs in this field . standard analysis yields the following expression for the parameters of the critical point @xmath127 where @xmath128 determines the equilibrium sf configuration at the critical field . the two canted phases are competing phases related by a first - order transition . at the field @xmath124 the discontinuity between these phases disappears , i.e. this is a _ critical end - point _ of the first - order transition . after the transition into the sf phase , the system further evolves by rotation of the magnetizations @xmath25 into the direction of the field . this process is terminated at another critical field @xmath129 ( fig . the transition between the canted and sf phases is continuous only below a certain anisotropy strength , @xmath130 , and becomes discontinuous at higher anisotropy . the parameters of the corresponding tricritical point @xmath121 , @xmath122 have been calculated numerically from the equation @xmath131 together with eqs . ( [ tran2 ] ) . numerically calculated transition fields  @xmath132 and the upper stability fields of the canted phases @xmath125 are plotted in fig . 8 . in this picture the grey area describes the region of the first - order transition , respectively the region of phase - coexistence between the canted phases . a magnetic field applied in one of the easy directions violates the energy balance between the two af states . the state af(i ) with staggered vector parallel to the field does not change its configuration and exists as a metastable collinear state up to the critical field @xmath133 . the phase af(ii ) with the staggered vector perpendicular to the magnetic field transforms into a sf phase which corresponds to the global energy minimum of the system . in an increasing magnetic field this phase evolves similarly to the sf phase along the hard directions , and continuously transforms into the flip states at the critical fields @xmath134 for low anisotropy . the transition between the sf and flip phases changes the order in another tricritical point ( @xmath135 ) . at a critical field @xmath136 ( [ stab1 ] ) the sf phase becomes unstable with respect to flipping distortions , i.e. modes that redirect the magnetization vectors into the field direction but preserve the symmetry of the sf phase . the parameters of the tricritical point ( @xmath137 ) . and the stability limits @xmath138 have been calculated in @xcite . after some algebra the first - order transition field @xmath139 also can be derived in analytical form @xmath140\ , , \label{tran1}\end{aligned}\ ] ] at the transitional field we have @xmath141 . magnetic phase diagram of the two - layer system in magnetic fields along _ easy _ directions . the tricritical point ( @xmath142 , @xmath143 ) separates the continuous @xmath134 ( [ stab1 ] ) ( dashed line ) and discontinuous @xmath134 ( [ tran1])(thick line ) transitions between sf and flip phases . the tricritical point ( @xmath144 , @xmath145 ) separates the continuous @xmath134 ( [ stab1 ] ) ( dashed line ) and discontinuous @xmath134 ( [ tran1])(thick line ) transitions between sf and flip phases . thin lines indicate the lability lines of the individual phases ( [ stab1 ] ) . the grey area designates the stability region of the canted phases . for the anisotropy @xmath146 = 0.47596 ( point @xmath147 ) , the critical field of the sf phase switches from an instability against symmetric `` flipping '' at the line @xmath136 to an instability against canting at @xmath148 given by the two branches @xmath149 from eq . ( [ stabsf2 ] ) . ( inset in main figure shows the details of this process ; see text following eq . ( [ stabsf2 ] ) for further explanation . ) small panels ( a ) and ( b ) show magnetization curves for low and high anisotropies as indicated in the main figure . full thick lines are for the evolution of the equilibrium ground - states , in ( a ) from sf phase to a continuous spin - flip transition into saturation . dotted lines are the evolution starting from the metastable af - phase . in ( b ) the evolution of the saturated state in decreasing fields is given by a thin line . the inner small hysteresis loop shows the transitions from and into the metastable canted phase . , width=287 ]    the ( k , h)-phase diagram for this cases is plotted in fig . 9 .    with increasing strength of the anisotropy , the `` landscape '' described by the energy function @xmath120 acquires additional folds . in particular , metastable canted phases arise in a certain region of the magnetic fields along the easy directions . ( this region is indicated by grey colour in fig . [ kheasy ] ) . for @xmath150 with @xmath151 = 0.28 + 0.08@xmath152 = 0.47596 the sf phase undergoes an instability with respect to a transition into the canted phases . this instability occurs at lability fields which are given by the two branches of the following parametric equation @xmath153 where @xmath154 is given by the the configuration in the sf - state , @xmath155 . both branches of the lability field ( [ stabsf2 ] ) start at the point @xmath156 ( point @xmath147 in fig . [ kheasy ] ) ; @xmath157 meets the lability line @xmath136 of eq . ( [ stab1 ] ) at the `` beak '' ( @xmath158 ) ( point @xmath159 in fig . [ kheasy ] ) . the transition between the metastable sf phase and the metastable canted phases along @xmath160 is continuous between the points @xmath159 and a further tricritical point @xmath161 . for stronger anisotropy the transition between these metastable phases is discontinuous at a first - order transition field @xmath162 . this transition line along with the stability limits of the canted phases have been numerically determined as given in fig . [ kheasy ] . note that such processes between metastable phases may be realized only in systems with high coercivity where the transition into the thermodynamically stable flip phase is suppressed at the field @xmath139 ( [ tran1 ] ) . now we consider the magnetic states in applied fields deviating from the symmetric directions and construct ( @xmath163)-phase diagrams for different values of the parameter @xmath164 . for the ( @xmath163)-diagrams in figs . [ hxhyweak ] to 14 it is assumed that the easy directions coincide with the @xmath3- and @xmath4-axis .     ( a ) schematic ( @xmath163)-phase diagram of the two - layer for @xmath165 . the astroid ( continuous line ) gives the stability limit of the canted metastable phases . thick lines in diagonal directions ( hard axis ) are first order transitions between canted phases . the line for the ( continuous ) transition into the saturated spin - flip state is given by the dotted line . note that the real size of the astroid in the ( @xmath163)-diagram is much smaller than shown here . ( b ) evolution of spin - flip field , respectively lower stability limits of the saturated ( ferromagnetic ) phase ( fm ) . phase fm is stable outside the closed curves for the various values of anisotropy @xmath166 . for @xmath167 the fm - phases are ( meta)stable even in reverse fields and the existence regions overlap for states saturated in different easy directions . , title="fig:",width=287 ]   ( a ) schematic ( @xmath163)-phase diagram of the two - layer for @xmath165 . the astroid ( continuous line ) gives the stability limit of the canted metastable phases . thick lines in diagonal directions ( hard axis ) are first order transitions between canted phases . the line for the ( continuous ) transition into the saturated spin - flip state is given by the dotted line . note that the real size of the astroid in the ( @xmath163)-diagram is much smaller than shown here . ( b ) evolution of spin - flip field , respectively lower stability limits of the saturated ( ferromagnetic ) phase ( fm ) . phase fm is stable outside the closed curves for the various values of anisotropy @xmath166 . for @xmath167 the fm - phases are ( meta)stable even in reverse fields and the existence regions overlap for states saturated in different easy directions . , title="fig:",width=287 ]    this corresponds @xmath168 in the energy eq . ( [ energy2 ] ) . diagrams for @xmath169 can be obtained by rotation of those for @xmath168 through an angle @xmath170 . in the limit of weak anisotropy , @xmath171 , independent minimization with respect to the angle @xmath172 yields @xmath173 and the energy ( [ energy2 ] ) is simplified to the following form @xmath174 \ , . \label{energy3}\end{aligned}\ ] ] the set of equations for the stability limits of solutions , @xmath175 and @xmath176 , yield two closed lines of critical fields in the ( @xmath177)-phase diagram ( fig .  [ hxhyweak ] ) . one of them @xmath178 describes the second - order transition into the flip phase with @xmath179 . the other closed curve can be written in a parametric form @xmath180 where @xmath181 , @xmath182 . it describes an astroid with eight - cusps that confines the region of the canted metastable states ( fig . [ hxhyweak](a ) ) . first quadrant of the ( @xmath163)-phase diagram for antiferromagnetic two - layers with ( a ) @xmath183 ( b ) @xmath184 . easy axes of magnetic anisotropy are along the axes of the diagrams . inset in ( a ) magnifies the region of the `` swallow tail '' for the discontinuous transitions in magnetic fields oriented close to the hard axes . thick lines are lines of first - order transitions . the dotted lines @xmath185  @xmath186 are continuous spin - flips from spin - flop phase to the saturated ( spin - flip ) state . thinner lines are the limits of stability , defining the coexistence regions for the corresponding phases . , title="fig:",width=287 ]   first quadrant of the ( @xmath163)-phase diagram for antiferromagnetic two - layers with ( a ) @xmath183 ( b ) @xmath184 . easy axes of magnetic anisotropy are along the axes of the diagrams . inset in ( a ) magnifies the region of the `` swallow tail '' for the discontinuous transitions in magnetic fields oriented close to the hard axes . thick lines are lines of first - order transitions . the dotted lines @xmath185  @xmath186 are continuous spin - flips from spin - flop phase to the saturated ( spin - flip ) state . thinner lines are the limits of stability , defining the coexistence regions for the corresponding phases . , title="fig:",width=287 ]    the cusps along the easy directions coincide with the stability field of the af phase @xmath133 eq . ( [ stab1 ] ) , and those along the hard directions with @xmath187 eq .  ( [ tran2 ] ) . thick black lines within the astroid and along the hard directions indicate the first - order transitions between the canted phases from fig . [ tetragonconfigs ] . this topology of the ( @xmath163)-phase diagram fig . [ hxhyweak](a ) is preserved up to the first tricritical point @xmath188 . for @xmath189 a first - order transition between sf phase and two canted phases arises along the hard directions ( point @xmath190 in fig . [ hxhyo3](a ) ) . for stronger anisotropies , @xmath191 , discontinuous phase transitions exist also for finite deviations of the field direction from the hard axes . the corresponding line of these first - order transitions have critical end - points , where the difference between the competing phases disappears ( analogously to the end - point of a coexistence - line in gas - liquid phase diagrams ) . for @xmath192 and for fields along the easy directions , the spin - flip into the saturated ( induced ) ferromagnetic state occurs discontinuously . thus , another line of first - order transitions develops for @xmath193 also in oblique fields . correspondingly , the transition lines for spin - flips in ( @xmath163)-phase diagrams consist of continuous and discontinuous sections joined by the tricritical points @xmath185 , @xmath186 . the calculated ( @xmath163)-phase diagrams for @xmath194 in fig . [ hxhyo3 ]  ( a ) and for @xmath195 in fig . [ hxhyo3 ]  ( b ) represent the main topological features of these phase diagrams . point ` @xmath190 ' indicates the field of phase equilibrium between the sf phase and two canted phases and corresponds to the critical field @xmath132 from fig . 8 ; points @xmath196 , @xmath197 are the end - points of the first - order transitions between one of the canted phases and the sf phase . points @xmath198 @xmath199 correspond to the transition field @xmath139 from fig . 9 ; points @xmath185 , @xmath186 are the tricritical points , where the discontinuous transition between the flip phase and the distorted sf phase ends . the lability lines in the vicinity of the transition field @xmath132 have the shape of a `` swallow tail '' . similar phase diagrams arise in uniaxial ferromagnets with strong fourth - order anisotropy.@xcite as the anisotropy constant increases in the region @xmath193 the lines of the first - order transitions extend and near the special value @xmath200=@xmath201 the points @xmath196 , @xmath197 reach the easy directions . for larger @xmath202 canted phases exist as metastable state also for fields in easy directions ( compare fig . however , for increasing fields close to this direction , the discontinuous transition from the spin - flop state into the spin - flip state occurs before the transition into such a canted phase can take place , i.e. the two first - order lines @xmath190  @xmath196 and @xmath185  @xmath198 cross each other close to the easy axis directions in fig . brhl et al . @xcite reported evidence of an intermediate canted state in a co / cu / co / cu(001 ) system with field in an easy direction , following the theoretical prediction of such states by dieny et al . @xcite however , the state was found at lower fields than expected . this may be due to a misorientation of the external field and/or a mosaic of the epitaxial layer system because then the canted state becomes stable already at lower fields , as can be seen from line @xmath190  @xmath196 in phase diagram fig . generally , it is not sufficient to investigate only magnetization behaviour in hard and easy axes directions for a thorough understanding of the magnetization phases in these multilayer systems .    to understand the transformation of ( @xmath163)-diagrams with increasing tetragonal anisotropy the limiting case is useful where @xmath25 are strictly oriented along the easy axes . in this case with _ infinite _ anisotropy , i.e. @xmath203 , our model can be considered as a chain of an antiferromagnetic four - state clock - model or planar potts - model@xcite in a transverse external field . recently , such four - state - clock models were employed to analyse the spin configurations of bulk tetragonal metamagnets with large four - fold anisotropies , such as rare - earth nickel - borocarbides @xcite or rare - earth silver - antimonides.@xcite rich experimental ( @xmath163)-phase diagrams have been obtained and analysed in terms of four - state - clock models , e.g. for honi@xmath204b@xmath204c @xcite or dyagsb@xmath204.@xcite the values of the angles @xmath205 0 , @xmath206 , @xmath207 , or @xmath208 in these four - fold states for @xmath209 are symbolically given by @xmath210 , @xmath211 , @xmath212 , and @xmath213 , e.g. the af and sf phases are ( @xmath210@xmath212 ) and ( @xmath213@xmath211 ) for @xmath0 . in addition to the collinear af and flip states there are phases with perpendicular orientation of the magnetization , or `` 90  degree - folded '' phases ( fig .  [ hxhyinfinite ] ) . the states created by all combinations of these four angles @xmath24 for the magnetization in the multilayer stack exist as metastable states in arbitrary magnetic fields because the energy wells corresponding to these solutions are separated by infinitely high potential barriers . the regions of the absolute stability of these phases are separated by first - order transitions lines shown in fig . [ hxhyinfinite ] . in the first quadrant the transition between the af states and the canted phase ( @xmath214 ) occurs at the line @xmath215 ; the transitions lines from the canted phase ( @xmath214 ) into the flip states ( @xmath216 ) and ( @xmath217 ) are @xmath218 and @xmath219 , correspondingly ( fig . 12 ) . in the points ( 1,0 ) and ( 0,1 ) four phases coexist . thus , under increasing strength of the anisotropy ( @xmath220 ) the ( @xmath163)-phase diagram evolves from that plotted in fig .  [ hxhyweak](a ) for @xmath221 to that in fig . [ hxhyinfinite ] for infinitely large values of @xmath222 . ( @xmath163)-phase diagram of the @xmath0-layer system in the limit of an infinite positive four - fold anisotropy . , width=321 ] the equilibrium magnetic configurations in two - layer systems arise as a results of the competition between the interlayer exchange coupling and tetragonal anisotropy . for multilayers the disbalance of the exchange forces at the boundaries ( sec . ii ) additionally influences the magnetic states . we first describe the structure of the one - dimensional solutions for laterally homogeneous states in finite antiferromagnetic superlattices . these magnetic states are determined by the interplay between cut exchange at the surfaces and the four - fold anisotropy . then we discuss some consequences for multidomain states and magnetization processes , and we discuss the physical nature of other effects which may play a role for the magnetic behaviour of real experimental systems . first , we investigate the ( @xmath163)-diagrams of magnetic states of multilayers with @xmath1 in the limit of infinite four - fold anisotropy . ( @xmath163)-phase diagram of the multilayer in the limit of infinite positive four - fold anisotropy . the phase - region of the l@xmath223x@xmath223-phases coincides with that of the collinear ferrimagnetic f@xmath223 , not shown here . ( for details see text . ) , width=287 ]    the zero - field ground - states are the af phases @xmath224 and @xmath225 we may restrict the field to be oriented in directions @xmath34 in the interval @xmath2260,@xmath207/4@xmath227 $ ] . the configurations for other values of @xmath34 follow from symmetry . then , we may distinguish the four ground - state ( zero - field ) domains ( af1)@xmath228= @xmath229 , ( af2)@xmath228= @xmath230 , ( sf1)@xmath228= @xmath231 , and ( sf2)@xmath228= @xmath232 , where @xmath233 signifies the number of repetitions of the pair in a domain . in external fields @xmath234 configurations with net magnetization will be stabilized . these configurations must have pairs of adjacent moments flipped by 90 degree l=@xmath235 or x=@xmath236 , and by 180 degree f=@xmath237 . the 90-degree - folded configuration `` x '' is less favourable than `` l '' for field orientations in the chosen range , except for fields in easy - axis direction @xmath238 . thus , the magnetized configurations with lowest magnetization and smallest expense of exchange energy are those with only one pair of type `` l '' . in a short - hand , we write l@xmath223 for these configurations of type ( af1)@xmath239l@xmath240(sf2)@xmath241 ( with @xmath242 . ) next , we may form configurations with higher magnetization and smallest expense of exchange energy using one `` l '' and one `` x''-pair : l@xmath223x@xmath223= ( af1)@xmath239l@xmath240(sf2)@xmath243x@xmath240(af1)@xmath244 , ( @xmath245 ) . these configurations have the same energy and , therefore , are degenerate with configurations containing only one `` f''-pair : f@xmath223=(af1)@xmath239f@xmath240(af2)@xmath241 . note that these states are highly degenerate because the `` l''- , `` x''-pair , or `` f''-pair may be placed at arbitrary positions in the stack of @xmath2 layers . in higher fields configurations with various combinations of `` l''- , `` x''- , and `` f''-pairs may be stabilized by an external field . however , the structures with lowest energy , i.e. the absolutely stable states , are rather simple because of the following considerations . if , in external fields , more than one `` l''-pair can be formed starting from the l@xmath223 structure , the formation of the maximum number of l - pairs gives the most favourable configuration . in particular , the structure l@xmath246=(@xmath214)@xmath246 is the lowest energy structure for fields pointing in hard axis - directions @xmath247 in the limit of infinite h. the state with higher magnetization than l@xmath246 for fields pointing closer to the easy axis direction , @xmath248 , and with the smallest expense of exchange energy is in our notation l@xmath249f . the saturated state fm@xmath250f@xmath246 is the most favourable state which is reached whenever energy can be gained in external fields by flipping more moments into position @xmath210 in the states l@xmath249f , f@xmath223 , or l@xmath223x@xmath223 . thus , for general @xmath2 we have only the following phases for infinite positive four - fold anisotropy : af or sf as zero - field ground - states ; for fields with orientation close to the hard direction @xmath247 l@xmath223 and l@xmath246 ; otherwise in intermediate fields two degenerate phases f@xmath223 and l@xmath223x@xmath223 ; at high fields an asymmetric l@xmath249f - phase and the fully saturated ferromagnetic phase fm . we did numerical checks to ascertain that no further energetically stable phases do exist in external fields of arbitrary strength and direction , indeed . thus , we searched for the states of lowest energy by sampling all possible configurations for models with @xmath2=4 , @xmath251 , 12 corroborating our arguments . based on this set of magnetic configurations the resulting ( @xmath163)-phase diagram for general @xmath252 can be calculated analytically ( fig . [ hxhyinfinite2 ] ) . as in the simpler case of fig . [ hxhyinfinite ] for the two - layer system all these states are separated by infinitely high potential barriers and remain metastable for arbitrary fields . the first - order transitions between different phases occur along straight lines as shown in fig . [ hxhyinfinite2 ] .    for finite strength of the four - fold anisotropy and under the influence of the exchange interactions , the basic structures are derived from the phases in fig . [ hxhyinfinite2 ] . under the influence of the field , they are elastically distorted into spatially _ inhomogeneous _ configurations . we have numerically investigated models for such cases with @xmath253 and various values of anisotropy @xmath222 . [ mvshn8 ] displays the general features for the example @xmath254 . the numerically calculated magnetization curves corresponding to the lowest energy states for sufficiently high anisotropies show the sequence of phases present in the infinite - anisotropy phase diagram fig . [ hxhyinfinite2 ] . however , the degeneracy of these phases is lifted because the distortions possible at finite anisotropy yield different gains of energy for the different configurations of the phases l@xmath223 , f@xmath223 , and l@xmath223x@xmath223 . generally , depending on field - orientation , large domains with nearly spin - flop - like configurations @xmath50 are favoured because these configurations can be more easily distorted by the fields yielding a corresponding gain of energy . in particular , phases derived from the ferrimagnetic collinear f@xmath223-phase are disfavoured compared to the l@xmath223x@xmath223-type phases . the nearly collinear ferrimagnetic configurations seem to exist only as metastable states . for fields closer to the hard - axis direction the phases with either one l - pair or @xmath255 l - pairs occur . under the influence of an applied field further discontinuous transitions occur at intermediate anisotropies , as seen e.g. in fig . [ mvshn8 ]  ( a ) for @xmath256 . these transitions are jumps from one energy basin to another which are formed by distorting degenerate configurations mainly of the type l@xmath223 and l@xmath223x@xmath223 .     examples of magnetization in an antiferromagnetic multilayer with @xmath254 . curves corresponding to the evolution of lowest energy states and for various values of anisotropy are shown with field oriented in directions of an easy axis ( a ) , in an oblique field ( b ) and in hard - axis direction ( c ) . @xmath34 is angle between magnetic field and easy axis direction . in each case , various magnetic configurations are shown for the discontinuous evolution of the magnetization curves with highest anisotropy . , title="fig:",width=287 ]   examples of magnetization in an antiferromagnetic multilayer with @xmath254 . curves corresponding to the evolution of lowest energy states and for various values of anisotropy are shown with field oriented in directions of an easy axis ( a ) , in an oblique field ( b ) and in hard - axis direction ( c ) . @xmath34 is angle between magnetic field and easy axis direction . in each case , various magnetic configurations are shown for the discontinuous evolution of the magnetization curves with highest anisotropy . , title="fig:",width=287 ]   examples of magnetization in an antiferromagnetic multilayer with @xmath254 . curves corresponding to the evolution of lowest energy states and for various values of anisotropy are shown with field oriented in directions of an easy axis ( a ) , in an oblique field ( b ) and in hard - axis direction ( c ) . @xmath34 is angle between magnetic field and easy axis direction . in each case , various magnetic configurations are shown for the discontinuous evolution of the magnetization curves with highest anisotropy . , title="fig:",width=287 ]    at low anisotropies the phase diagram attains the behaviour discussed for zero anisotropies in section iii . here the magnetic states are only influenced by the cuts of exchange bonds at the surfaces and spatially inhomogeneous spin - flop states slightly distorted by four - fold anisotropy are realized ( see figs . ( @xmath163)-diagrams for this isotropic case consist of the region of the inhomogeneous sf phase ( fig .  2 , 3 ) separated from the saturated flip state by the critical line @xmath115 ( fig . these global features of the phase - diagrams are similar to the case of antiferromagnetic multilayers with uniaxial anisotropy.@xcite    up to now we have discussed the effects of the competition between the bilinear exchange and four - fold anisotropy . it is clear that finite _ biquadratic coupling _ may substantially change the critical fields and stability regions of the various phases as it affects the elastic stiffness of the system . we note further , that the simple structure of the phase - diagram ruled by the phases present in the infinite anisotropy limit of fig . [ hxhyinfinite2 ] is valid only for the case of equal exchange constants @xmath257 in the multilayer stack . for arbitrary sets of values for @xmath14 in energy ( [ energy1 ] ) the phase diagrams may become considerably more complex and may contain further phases with different combinations of flipped 90 and 180 degree spin - pairs . even in such cases , the outline of the magnetic phase - diagrams , described here for finite equal - constant superlattices , should generally hold : ( i ) the high anisotropy limit is comparably simple with few phases determined by the competing lowest - energy basins of the anisotropy . non - equal constants may lift some degeneracies that are present in superlattices . ( ii ) in intermediate anisotropy range many elastically distorted phases appear , which are derived from stable and metastable high - anisotropy phases . ( iii ) for vanishing anisotropy the phase - diagrams become simple again as only inhomogeneous spin - flop - like phases and the saturated ferromagnetic phase remain in external fields . so far , we have analyzed single domain magnetic configurations . magnetization processes in real systems , however , are usually accompanied by complex reconstruction of the multidomain patterns as those recently observed in fe / cr multilayers.@xcite    there are two main physical mechanism for multidomain states in the systems under discussion . four - fold degeneracy of the ground - state leads to creation of _ antiferromagnetic multidomain structures _ with 90 degree domain walls . unlike the case of magnetic materials with nonzero total magnetization , where multidomain states are caused by demagnetization effects , in antiferromagnets such domains are _ metastable _ and arise during the formation of the ordering state , i.e. they have a kinetic origin.@xcite hence , multidomain patterns observed in antiferromagnetic coupled two - layers and multilayers have irregular morphologies and depend on thermal and magnetic - field histories.@xcite    however , another type of multidomain structures arises in the vicinity of field - induced discontinuous transitions.@xcite such thermodynamically stable _ transitional domain structures _ are formed by domains from states corresponding to the coexisting phases at first - order transitions . these domains are analogous to the domains of a demagnetized ferromagnet . in principle , the equilibrium parameters of such multidomain structures and their boundaries can be calculated by standard methods.@xcite in @xcite such calculations have been carried out for bulk easy - plane tetragonal antiferromagnets . _ magnetoelastic interactions _ lead to modification of the inhomogeneous states and decrease the regions of the multidomain states up to their complete suppression.@xcite this and other coercivity mechanism partly block the development of the equilibrium states . as a result , in real systems the evolution of multidomain states @xcite is accompanied by rather strong hystereses.@xcite for experiments , there are two important consequences related to the starting states and history dependence of magnetization processes . the `` texture '' of a real antiferromagnetic state in zero - field depends on the detailed kinetics imposed by e.g. cooling rates or deposition conditions . on the other hand , field cycling by inner loops for an antiferromagnetic multilayer may yield various metastable configurations and multidomain structures owing to the very wide coexistence region of the domains whenever sizable magnetic anisotropies are present in the multilayer stack within a phenomenological approach we develop the theory of reorientation transition in antiferromagnetically coupled superlattices with in - plane magnetization . detailed investigations of the surface effects in the isotropic multilayers ( sec . iii ) and four - fold anisotropy effects in two - layer systems ( sec . iv ) reveal the most important features of the system : ( i ) complex evolution of the inhomogeneous states ( figs .  2 - 4 ) is imposed by the strong disbalance of the exchange coupling , i.e. by the _ cut of the exchange bonds_. ( ii ) remarkable field - induced reorientational processes occur due to enhanced four - fold anisotropy ( figs . 8 to [ mvshn8 ] ) . the model used here , corresponding solutions and phase - diagrams include as special cases earlier theoretical studies on surface @xcite and four - fold anisotropy @xcite effects . the results are in essential accordance with existing experimental observations on inhomogeneous distribution in the spin - flop phase near the saturation field @xcite , and some effects of four - fold anisotropy in co / cu(001 ) wedged two - layers @xcite and fe / cr superlattices.@xcite our approach and results enable a _ qualitative _ analysis of the magnetization processes in the multilayered systems ( sec . v ) . in spite of the rather complex phase - diagrams of these systems , the analysis can be extended towards a quantitative description of real systems belonging to the class of artificial layered antiferromagnets described by eq . ( [ energy1 ] ) .    so far most experimental work is carried out only for special conditions , often data are collected only with fields along easy - axes . these results only cover small regions of the ( @xmath258,@xmath259)-phase diagrams ( figs .  [ hxhyweak ] to [ hxhyinfinite2 ] ) and do not capture the rich varieties of magnetizations processes in such systems . it is remarkable that many interesting effects , such as 90-degree folded phases , transitions into asymmetric canted phases etc . , are present already in antiferromagnetically coupled two - layers.@xcite as we have seen , phase diagrams for multilayer systems with @xmath260 become very complex . hence , systems with few layers are probably a better starting point for detailed investigations on magnetization processes . such experiments could be used to assess magnetic parameters and quality of such systems . generally , magnetization processes and checks for their dependence on magnetic and thermal pre - history should be made by applying fields in oblique directions and/or under rotating fields . only then , the behaviour of the ( @xmath258,@xmath259)-phase diagrams can be usefully compared with detailed theoretical investigations . further theoretical work may address models with non - equal constants . also effectively ferrimagnetic systems with odd numbers @xmath2 of layers and with different layer thicknesses may be interesting . some experimental data for magnetization processes in such systems exist , e.g. for co / cr two - layers with different ferromagnetic layers thicknesses,@xcite and for odd - numbered multilayers.@xcite interesting reorientational effects should arise also in magnetic fields , which are applied perpendicular or under arbitrary angle to the layer plane . antiferromagnetically coupled superlattices may also undergo transitions into perpendicularly magnetized states for certain thicknesses of the individual ferromagnetic layers , as observed for co / cr(001).@xcite in such cases the stray field must be taken into account already for laterally homogeneous states . concluding we state that the magnetic effects and phenomena discussed in this paper can be used for detailed investigations on such aspects of nanomagnetism as interlayer - exchange interactions , reorientational , and multidomain processes . g. j. strijkers , s. m. zhou , f. y. yang , and c. l. chien , phys . b * 62 * , 13896 ( 2000 ) ; k. y. kim , s. h. jang , k. h. shin , h. j. kim , and t. kang , j.  appl . phys .  * 89 * , 7612 ( 2001 ) ; a. moser _ et al . _ , j.  phys . d * 35 * , r157 ( 2002 ) . m. a. howson _ et al . _ , b * 49 * , 9560 ( 1994 ) . r. w. wang , d. l. mills , e. e. fullerton , j. e. mattson , and s. d. bader , phys . lett .  * 72 * , 920 ( 1994 ) ; s. g. e. te velthuis , j. s. jiang , s. d. bader , and g.p . felcher , phys . lett .  * 89 * , 127203 ( 2002 ) . p. steadman _ et al . _ , lett .  * 89 * , 077201 ( 2002 ) ; s. g. e. te velthuis , j. s. jiang , g. p. felcher , appl . .  lett .   * 77 * , 2222 ( 2000 ) . o. hellwig , t. l. kirk , j. b. kortright , a. berger , and e. e. fullerton , nature materials * 2 * , 112 ( 2003 ) . k. ounadjela _ _ , phys . b * 45 * , 7768 ( 1992 ) ; s. hamada , k. himi , t. okuno , and k. takanashi , j.  magn .  magn . * 240 * , 539 ( 2002 ) . h. itoh , h. yanagihara , k. suzuki , and e. kita , j.  magn . mater .  * 257 * , 184 ( 2003 ) . a. n. bogdanov and u. k. rler , phys . b * 68 * , 012407 ( 2003 ) ; u. k. rler and a. n. bogdanov , cond - mat/0307494 ( unpublished ) . k. brhl , s. di nunzio , f. schreiber , th . zeidler , and h. zabel , j.  appl . * 75 * , 6184 ( 1994 ) ; k. brhl , th . zeidler , f. schreiber , a. schreyer , and h. zabel , j.  magn . mater .  * 130 * , l1 ( 1994 ) . a. van der graaf , m. valkier , j. kohlhepp , and f. j. a. den broeder , j.  magn . * 165 * , 157 ( 1997 ) . r. w. e. van de kruijs _ et al . _ , phys . b * 65 * , 064440 ( 2002 ) . t. zeidler , f. schreiber , and h. zabel , j.  appl . phys .  * 79 * , 4793 ( 1996 ) . j. j. picconatto , m. j. pechan , and e. e. fullerton , j.  appl . phys .  * 81 * , 5058 ( 1997 ) . j. z. hilt , j. j. picconatto , a. obrien , m. j. pechan , and e. e. fullerton , j.  magn . magn .  mater .  * 198 * -*1999 * , 387 ( 1999 ) . v. v. ustinov , m. a. milyaev , l. n. romashev , t. p. krinitsina , e. a. kravtsov , j.  magn . * 226 * -*230 * , 1811 ( 2001 ) . c. chesman _ _ , phys .  rev . b * 58 * , 101 ( 1998 ) .    v. lauter - pasyuk , h. j. lauter , b. p. toperverg , l. romashev , and v. ustinov , phys . 89 * , 167203 ( 2002 ) ; v. lauter - pasyuk _ et al . _ , j.  magn . magn .  mater . * 258 * -*259 * , 382 ( 2003 ) . k. temst , physica b  * 276 - 278 * , 684 ( 2000 ) . d. l. nagy _ et al . _ , lett .  * 88 * , 157202 ( 2002 ) ; l. bottan _ et al . _ , j.  magn . mater .  * 240 * , 514 ( 2002 ) . h. lauter _ et al . _ , j.  magn . * 258 * -*259 * , 338 ( 2003 ) . s. o. demokritov , j.  phys . d : appl . * 31 * , 925 ( 1998 ) . b. dieny , j. p. gavigan , and j. p. rebouillat , j.  phys . : condens . matter * 2 * , 159 ( 1990 ) ; b. dieny and j. p. gavigan , j.  phys . : condens . matter * 2 * , 187 ( 1990 ) . f. c. nrtemann , r. l. stamps , a. s. carrio , and r. e. camley , phys . b * 46 * , 10847 ( 1992 ) ; a. l. dantas and a. s. carrio , phys .  rev . b * 59 * , 1223 ( 1999 ) . a. hubert , r. schfer , _ magnetic domains _ ( springer - verlag , berlin , 1998 ) . d. l. mills , j.  magn . mater .  * 198 * -*199 * , 334 ( 1999 ) . p. grnberg , r. schreiber , y. pang , m. b. brodsky , and h. sowers , phys . lett .  * 57 * , 2442 ( 1986 ) ; s. s. p. parkin , n. more , and k. p. roche , phys .  rev .  lett .  * 64 * , 2304 ( 1990 ) . m. d. stiles , j.  magn . mater .  * 200 * , 322(1999 ) . v. g. baryakhtar , a. n. bogdanov , and d. a. yablonskii , usp . * 156 * , 47 ( 1988 ) [ sov . usp .  * 31 * , 810 ( 1988 ) ] . l. d. landau , e. m. lifshitz , _ statistical physics _ part i. course of theoretical physics ; v.  5 ( butterworth - heinemann , oxford 1980 ) . f. y. wu , rev .  mod . phys .  * 54 * , 235 ( 1982 ) . mller and v. n. narozhnyi , rep . phys .  * 64 * , 943 ( 2001 ) ; p. c. canfield and s. l. budko , j.  alloys compd .  * 262 * , 169 ( 1997 ) ; g. andr _ et al . _ , physica b * 292 * , 176 ( 2000 ) . p. c. canfield _ et al . _ , b * 55 * , 970 ( 1997 ) ; v. a. kalatsky and v. l. pokrovsky , phys . b * 57 * , 5485 ( 1998 ) . k. d. myers , p. c. canfield , v. a. kalatsky , and v. l. pokrovsky , phys . b * 59 * , 1121 ( 1999 ) . a. hubert , _ theorie der domnenwnde in geordneten medien _ ( springer , berlin , 1974 ) .
Abstract:

Answer: in antiferromagnetically coupled superlattices grown on ( 001 ) faces of cubic substrates , e.g. based on materials combinations as co / cu , fe / si , co / cr , or fe / cr , the magnetic states evolve under competing influence of bilinear and biquadratic exchange interactions , surface - enhanced four - fold in - plane anisotropy , and specific finite - size effects . using phenomenological ( micromagnetic ) theory , a comprehensive survey of the magnetic states and reorientation transitions has been carried out for multilayer systems with even number of ferromagnetic sub - layers and magnetizations in the plane . in two - layer systems ( @xmath0 ) the phase diagrams in dependence on components of the applied field in the plane include `` swallow - tail '' type regions of ( metastable ) multistate co - existence and a number of continuous and discontinuous reorientation transitions induced by radial and transversal components of the applied field . in multilayers ( @xmath1 ) noncollinear states are spatially inhomogeneous with magnetization varying across the multilayer stack . for weak four - fold anisotropy the magnetic states under influence of an applied field evolve by a complex continuous reorientation into the saturated state . at higher anisotropy they transform into various inhomogeneous and asymmetric structures . the discontinuous transitions between the magnetic states in these two - layers and multilayers are characterized by broad ranges of multi - phase coexistence of the ( metastable ) states and give rise to specific transitional domain structures . donetsk institute for physics and technology , 340114 donetsk , ukraine