the present knowledge of solar and stellar magnetism derives essentially from the interpretation of polarization in spectral lines due to the zeeman effect . in the presence of a weak magnetic field ,
the radiative transfer equation for polarized light can be solved analytically and it is known as the `` weak field approximation '' .
although very simple , it relies on a set of strong assumptions that we discuss later .
however , the weak field approximation is applicable in many scenarios and gives very good results for the inference of magnetic fields as compared with more elaborate methods . in solar physics , it is at the heart of the success of a large number of synoptic magnetographs , like those of big bear @xcite .
it has been also used recently to produce vector magnetograms with the imaging magnetograph experiment ( imax ) instrument @xcite onboard the sunrise balloon @xcite . in stellar spectropolarimetry ,
the weak field approximation is at the base of the least - squares deconvolution ( lsd ; * ? ? ?
* ) , the most successful technique used to detect and measure magnetic fields in solar - type stars , or in some recent works on central stars of planetary nebulae @xcite , white dwarfs @xcite , pulsating stars @xcite , hot subdwarfs @xcite and ap and bp stars @xcite .    most of the scientific cases commented so far deal with unresolved structures which means that the weak field approximation has to be understood in terms of fluxes of the stokes parameters instead of specific intensities .
recently , @xcite derived the weak field expressions in terms of fluxes , particularizing to the case of the stellar dipole .
they neglected rotation to simplify the equations and to derive analytical expressions for the inference of the magnetic field .
when rotation is taken into account ( in other words , there is a correlation between the magnetic field and the line - of - sight velocity ) , the problem becomes much more challenging and it is usually solved numerically ( e.g. , * ? ? ?
* for a recent effort ) . however , we show in this paper that it is possible to obtain analytical expressions for the stokes vector fluxes in the weak field approximation .
there are two interests on this effort .
first , having an analytical expression for the stokes flux might be of help to investigate the interplay between rotation and the magnetic field .
second , it can be used to speed up inversion codes based on the weak - field approximation because no numerical integration is needed .
this will facilitate the application of bayesian inference codes , like the one developed by @xcite , that are based on markov chain monte carlo methods that need to carry out the synthesis thousands of times to sample the posterior distribution .
this is the first paper of a series dealing with analytical expressions for the fluxes of the stokes parameters in the weak field limit
. this paper deals with the expressions for a rotating star with a magnetic field on its surface .
the results depend on the correlation between the magnetic field and the velocity and are general for any rotation profile or any magnetic field configuration .
we particularize the equations to the case of a rotating magnetic stellar dipole .
let us consider a star with arbitrary velocity and magnetic fields on its surface .
we define both fields for an arbitrary point @xmath0 in the surface using the reference frame represented in fig .
[ fig : fig_geom ] .
the line - of - sight ( los ) is along the @xmath1 axis , while we choose the axis @xmath2 as the one defining the reference direction for positive stokes @xmath3
.     axis .
the direction @xmath2 defines the reference for @xmath4 .
the vector @xmath0 defines any position on the stellar surface.,scaledwidth=50.0% ]    our aim is to compute the local stokes profiles emerging from every point of the visible surface and compute its integral . as explained in the introduction , a simple ( although subject to approximations ) way to relate the emergent stokes profiles with the magnetic field vector
is the weak - field approximation .
this approximation ( combined with the weak - line approximation for linear polarization ) leads to an analytical solution to the radiative transfer equation .
it is based on the following three fundamental assumptions : i ) the magnetic field vector is constant along the line of sight on the region of the stellar atmosphere where the spectral line is formed , ii ) the magnetic field intensity is sufficiently weak so that the zeeman splitting is smaller than any other line broadening mechanism and iii ) the los velocity and any broadening mechanism are constant with height in the line formation region . under these simplifying conditions ,
the magnetic field can be considered as a perturbation to the zero - field case .
this approximation is not strictly valid in general , but it is a good approximation whenever the magnetic field strength is weak in the line formation region in the absence of strong gradients along the los .
quantitatively , the approximation holds whenever @xmath5 , where @xmath6 ( see eq .
( 3.14 ) in * ? ? ?
* ) is the zeeman splitting that is proportional to the field strength , @xmath7 is the dominant broadening mechanism and @xmath8 represents the effective land factor of the line , which quantifies the magnetic sensitivity of the line and that only depends on the quantum numbers of the transition ( see , e.g. , * ? ? ? * ) . from its definition , it is clear that the weak field regime occurs at different field strengths for different spectral lines , for different stellar objects and for different spectral resolutions . for instance , for a magnetically sensitive iron spectral line with @xmath9 at 5000   with a standard thermal broadening at @xmath10 k , the weak - field regime is applicable for fields below @xmath11600 g if the spectral resolution of the spectrograph is smaller than the line broadening . under the previous simplifying assumptions , and to first order in @xmath6 , the intensity profile @xmath12 of a spectral line formed in a weak magnetic field is fully insensitive to the magnetic field . in other words , it fulfills the standard radiative transfer equation in the absence of a magnetic field ( e.g. , * ? ? ?
. likewise , at first order in @xmath6 , the circular polarization profile , i.e. , the stokes @xmath13 profile for a given spectral line , has the following expression : @xmath14 where the symbol @xmath15 stands for the longitudinal component of the magnetic field ( the projection along the los in the reference system shown in fig .
[ fig : fig_geom ] ) .
the symbol @xmath16 stands for the wavelength while @xmath17 is the central wavelength of the transition . the constant @xmath18 .    at first order in @xmath6 ,
linear polarization is zero . in order to obtain an expression for
the stokes profiles characterizing linear polarization we have to expand the radiative transfer equation to second order in @xmath6 and assume that the spectral line is not saturated and weak .
under these assumptions , a good approximation to the intensity profile emerging from the stellar atmosphere is to consider it , in general , proportional to a voigt function ( see eqs .
( 9.84 ) and ( 9.86 ) in * ? ? ?
for the sake of simplicity , we consider the case of a gaussian line profile .
all our results can be derived for a voigt profile but the final expressions are more complicated without any fundamental difference .
this assumption somehow reduces the scope of applicability of our results for linear polarization because the line shape self - similarity is only valid for lines whose line depth is smaller than @xmath1140% according to the results of @xcite .    with this assumption , the weak - field approximation for the stokes @xmath19 and @xmath20 , defined in the @xmath21 reference system shown in fig .
[ fig : fig_geom ] , are ( see * ? ? ?
* ; * ? ? ?
* ) : @xmath22 where the symbol @xmath23 plays the role of the effective land factor for linear polarization and quantifies the sensitivity of linear polarization to the magnetic field .
again , it is only a function of the quantum numbers of the transition @xcite .
the symbols @xmath24 and @xmath25 stand for the components of the magnetic field vector along the @xmath2 and @xmath26 axes , respectively .
the stokes profiles in any other reference system that is defined as a rotation around the @xmath1 axis by an angle @xmath27 are given by the usual rotation @xcite : @xmath28 where the quantities with the subindex `` new '' refer in this case to the new reference frame .
when observing unresolved objects , the detected polarized radiation is obtained as an integration in the plane of the sky of the individual stokes parameters at each point of the stellar surface ( @xmath29 ) , all of them referred to a common reference system .
we describe the stellar disk using the polar coordinates @xmath30 and @xmath31 , where @xmath30 is normalized to the stellar radius , so that its range is from 0 to 1 ( see * ? ? ?
therefore , the average of any arbitrary function @xmath32 on the surface @xmath29 will be given by : @xmath33 as a consequence , it is obvious that the observed values of the stokes flux vector depends on the surface distribution of the magnetic field , on the centre - to - limb variation ( clv ) of the radiation , and on the doppler effect due to the rotation of the star .
following standard approximations and assuming that the star has azimuthal symmetry , @xcite parameterized the intensity at each point of the star as a product of two functions :
@xmath34 , that only depends on the wavelength ( the intensity at the disc centre since they neglect rotation ) , and @xmath35 , the clv that depends on the position on the disk and is wavelength independent ( e.g. , * ? ? ?
* ; * ? ? ?
* ) : @xmath36 the clv is given in terms of @xmath37 , where @xmath38 is the astrocentric angle between the normal to a point in the stellar surface and the line of sight ( see fig . [
fig : fig_geom ] ) . the parameters @xmath39 and @xmath40 have values between 0 and 1 and
are supposed to be constant along the spectral line ( although they can vary from line to line ) . with this assumption
, we also force the line depth to be @xmath41-independent .
note that the values of @xmath39 and @xmath40 have to fulfill the condition that the intensity has to be positive .
this simple parameterization allowed them to express the weak field equations in terms of the observed fluxes , thus mimicking the expressions for the specific stokes vector .
this was possible because the intensity flux is assumed to be the product of the wavelength - dependent intensity at disk center and the wavelength independent clv . using this relationship ,
the weak - field and weak - line approximation for the integrated polarized flux remains formally the same than for the specific stokes vector but the components of the magnetic field appear weighted by the clv law .
when the velocity field is included , the situation turns out to be much more complicated because the intensity itself depends on the local projection of the velocity along the @xmath1 axis ( @xmath42 ) , which subsequently depends on the position along the stellar disk in a potentially complicated way . in this case , the fluxes of the stokes vector are expressed as : @xmath43\frac{\partial^2 { i}(\lambda , v_z(\rho,\theta),\mu)}{\partial \lambda^2 } \nonumber \\
f_{u}(\lambda)&=&-\frac{{c}^2}{2}\lambda_0 ^ 4 g \int_0 ^ 1 \rho d\rho \int_0^{2\pi } d\theta~b_x(\rho,\theta)b_y(\rho,\theta ) \frac{\partial^2 { i}(\lambda , v_z(\rho,\theta),\mu)}{\partial \lambda^2 } \nonumber \\
f_{v}(\lambda)&=&-{c}\lambda_0 ^ 2 g \int_0 ^ 1 \rho d\rho \int_0^{2\pi }
d\theta~b_z(\rho,\theta ) \frac{\partial { i}(\lambda , v_z(\rho,\theta),\mu)}{\partial \lambda},\end{aligned}\ ] ] the immediate consequence of the dependence on @xmath44 is that the integrals can not be separated .
to proceed , it is crucial to decouple the rotation velocity and the wavelength to end up with expressions for the weak field approximation in terms of the observed fluxes . in order to do that
, we assume that the shape of the local absorption profile is independent of the position on the stellar surface , apart from the doppler shift introduced by @xmath42 , and that the clv is wavelength independent . furthermore ,
if we assume that the absorption profile is gaussian , we end up with : @xmath45\right ] , \label{eq : intensity_gaussian}\ ] ] where @xmath46 stands for the continuum intensity , @xmath47 for the line depth , and @xmath48 for the line width .
this approximation is not appropriate for stars with large radial temperature and density gradients because the lines close to the limb will differ strongly from the lines at the disk center .
apart from the gaussian shape , it is possible to work with voigt functions but the subsequent expressions turn out to be more complicated and no new physics is gained , so we prefer to stick with the gaussian shape .
the gaussian shape is also consistent with the assumption of an unsaturated and weak spectral line that is necessary for eqs .
( [ eq : q_u ] ) to hold .
we can carry out the taylor expansion of @xmath49 around @xmath50 ( we drop now the @xmath41-dependence ) to obtain the following expression for the intensity at each point of the stellar surface : @xmath51 defining the emission line profile at zero velocity as @xmath52,\ ] ] the absorption profile at an arbitrary velocity shift can be written as : @xmath53    taking into account that the gaussian function is the generating function of the hermite polynomials : @xmath54=(-1)^n\exp\left[-\frac{x^2}{2 } \right ] \textit{he}_n(x),\ ] ] where @xmath55 is the hermite polynomial of order @xmath56 with weight function @xmath57 , also known as chebyshev - hermite polynomials ( see * ? ? ?
if a voigt function is used instead for the local profile , the derivatives are more complex and the recurrence relations of @xcite can be used . taking the previous definition into account , the derivatives ( @xmath58 ) of the intensity profile can be expressed as : @xmath59 so that the intensity profile of eq .
( [ eq : intensity_gaussian ] ) is given by the following gram - charlier series : @xmath60,\ ] ] where we have extended the summation to @xmath61 to include the profile at zero velocity , since @xmath62 .
the previous formula shows that the intensity is given as the addition of infinite terms ( potentially a large number of them for obtaining convergence ) given by the product of a function that depends on the wavelength and another one that depends on the velocity along the line of sight .
the latter is the only one that changes on the stellar surface , so that this variable separation allows us to carry out the surface integration . in order to simplify the following computations , it is interesting to define the function @xmath63 so that the intensity profile and its first and second derivatives with respect to the wavelength are given by : @xmath64 \nonumber \\
{ i}'(\lambda , v_z , \mu)&= & -\psi'(\lambda ) f(\mu ) \mathcal{h}(\lambda , v_z ) - \psi(\lambda ) f(\mu ) \mathcal{h}'(\lambda , v_z ) \nonumber \\ { i}''(\lambda , v_z , \mu)&= & -\psi''(\lambda ) f(\mu ) \mathcal{h}(\lambda , v_z ) - 2\psi'(\lambda ) f(\mu ) \mathcal{h}'(\lambda , v_z ) -\psi(\lambda ) f(\mu ) \mathcal{h}''(\lambda , v_z),\end{aligned}\ ] ] where the @xmath65 and @xmath66 symbols represent the first and second derivatives with respect to the wavelength .
applying the previous developments and dropping the dependences of the variables to avoid crowding , the polarized flux in the weak - field weak - line approximation can be written as : @xmath67 \nonumber \\
f_{u}(\lambda ) & = & \frac{{c}^2}{2}\lambda_0 ^ 4 g   \left [ \psi''(\lambda ) \overline{b_xb_yf\mathcal{h } } + 2\psi'(\lambda ) \overline{b_xb_yf\mathcal{h } ' } + \psi(\lambda ) \overline{b_xb_yf\mathcal{h } '' } \right ] \nonumber \\
f_{v}(\lambda)&=&{c}\lambda_0 ^ 2 g \left [ \psi'(\lambda ) \overline{b_zf\mathcal{h } } + \psi(\lambda ) \overline{b_zf\mathcal{h } ' } \right ] , \label{eq : fi_star}\end{aligned}\ ] ] where we remind that the overlines represent integration in the plane of the sky , as indicated in eq .
( [ eq : integration_sky ] ) .
it is important to note that the area asymmetry of the polarized flux is zero as can be verified by computing the wavelength integral of the previous equations ( it is easy to verify that the integral is zero by using the properties of the hermite polynomials ) .
the symbol @xmath68 stands for the continuum intensity flux and it is trivially given by : @xmath69 as expected , the stokes parameters depend on quantities averaged over the stellar surface . using the expansion on hermite polynomials of the absorption spectral line
, the terms appearing in the previous equations can be written as the following summations : @xmath70 the functions @xmath71 are given by : @xmath72    it is important to note that eqs .
( [ eq : fi_star ] ) and ( [ eq : averages_surface ] ) are formally invariant for any magnetic field configuration , rotation profile and center - to - limb variation once the assumptions stated before ( sufficiently weak magnetic field , wavelength independent clv , unsaturated and weak gaussian line ) are fulfilled .
the information about the magnetic configuration , rotation and the clv is encoded in the @xmath71 functions , which may depend on additional external parameters .
this is the case of the dipolar field with solid body rotation ( see sec . [
sec : dipole ] ) , in which these functions depend on the orientation of the dipole and the angular velocity vector .
the integrated stokes parameters also depend on terms that include the wavelength derivative of the @xmath73 function .
thanks to the effective separation between wavelength and velocity , these are easily obtained from the previous ones by using known properties of the hermite polynomials .
for instance : @xmath74 and the rest of terms are computed following the same strategy .    when the star is not rotating , the only terms contributing to the summations in eq .
( [ eq : averages_surface ] ) are those with @xmath61 .
likewise , the terms with derivatives of eqs .
( [ eq : averages_surface_d ] ) are always zero .
the polarized fluxes of eqs .
( [ eq : fi_star ] ) can be used to explain stellar observations using a standard optimization method to fit the functional forms to the observations .
however , we note that the expressions for the polarized fluxes depend on @xmath47 and @xmath48 , properties of the local line profiles to which we do not have direct access .
it is possible , though , to drop the dependence on @xmath47 from the equations if we use the fact that @xmath75 , @xmath76 and @xmath77 can be written in terms of @xmath78 , @xmath79 and @xmath80 , which are , otherwise , observed quantities .
this is the strategy followed by @xcite to write the weak - field weak - line approximation for a non - rotating star in terms of the observed fluxes ( and the parameters of the clv and the gaussian profile ) .
so : @xmath81 which gives , after isolation of @xmath75 , @xmath76 and @xmath77 : @xmath82 \nonumber \\
f_{u}(\lambda ) & = & -\frac{{c}^2}{2 } \lambda_0 ^ 4 g \left [ \frac{\overline{b_xb_y f \mathcal{h}}}{\overline{f\mathcal{h } } } f''_\mathcal{i } + 2\left ( \frac{\overline{b_xb_y f \mathcal{h}'}}{\overline{f\mathcal{h } } } - \frac{\overline{b_xb_y f \mathcal{h}}~\overline{f\mathcal{h}'}}{\overline{f\mathcal{h}}^2 } \right ) f'_{i}\right .
\nonumber \\   & + & \left .
\left ( \frac{\overline{b_xb_y f \mathcal{h}''}}{\overline{f\mathcal{h } } } - \frac{\overline{b_xb_yf\mathcal{h}}~\overline{f\mathcal{h}''}}{\overline { f\mathcal{h}}^2 }   + 2\frac{\overline{b_xb_yf\mathcal{h}}~\overline{f\mathcal{h}'}^2}{\overline { f\mathcal{h}}^3 } -2\frac{\overline{b_xb_y\mathcal{h}'}~\overline{f\mathcal{h}'}}{\overline { f\mathcal{h}}^2 } \right ) ( f_c - f_{i})\right ] \nonumber \\
f_{v}(\lambda ) & = & -{c } \lambda_0 ^ 2 g \left [ \frac{\overline{b_z f \mathcal{h}}}{\overline{f\mathcal{h } } } f'_{i }   +   \left ( \frac{\overline{b_z f \mathcal{h}}~\overline{f\mathcal{h}'}}{\overline{f\mathcal{h}}^2 } - \frac{\overline{b_z f \mathcal{h}'}}{\overline{f\mathcal{h}}}\right ) ( f_c - f_{i } ) \right ] \label{eq : fi_star_observer}\end{aligned}\ ] ] these expressions are the generalization of the weak - field weak - line approximation to an unresolved star with an arbitrary velocity and magnetic field on the surface .
the dependency on the local absorption @xmath47 has disappeared once the observed @xmath83 is used . however , the local broadening is still present and has to be either inferred from the observations or fixed to a reasonable value .
the set of eqs .
( [ eq : fi_star_observer ] ) represent a description of the stokes parameters in terms of observed quantities and also in terms of wavelength - dependent unknowns that are explicit functions of the field and velocity distribution on the star ( e.g. , @xmath84 ) .
both contributions are conveniently separated . in principle
, one could use these expressions to obtain information about the magnetism of the star following two strategies .
the first one is to parameterize these terms by fixing a model for the magnetic structure and infer the values of the parameters .
this is the approach we follow in the next section .
the second one is to non - parametrically infer these wavelength - dependent unknowns directly using a suitable inversion procedure .
since this would require additional regularization techniques , a full treatment of the problem would have to be done in a bayesian framework .
an interesting property of the integrated stokes parameters when rotation is taken into account is that the circular polarization flux depends on the derivative with wavelength of the flux ( @xmath85 ) but also on the flux itself ( @xmath83 ) .
likewise , the linear polarization fluxes depend on the flux and its first and second order derivative with respect to the wavelength . when particularizing the polarized flux of eq .
( [ eq : fi_star_observer ] ) to the non - rotating star , we recover the equations developed by @xcite .
this is a consistency check of our equations .
let us consider a star that is rotating with an angular velocity vector @xmath86 , as shown in fig .
[ fig : fig_geom_dipole ] . in this section
we use two different reference systems for our computations .
the arbitrary reference system @xmath87 , shown in the left panel of fig . [ fig : fig_geom_dipole ] , is the one associated with an observer .
the line - of - sight ( los ) is along the @xmath1 axis and the axis @xmath88 ( arbitrarily chosen by the observer ) defines the reference direction for positive stokes @xmath3 .
the angular velocity vector is defined in this reference system by its modulus ( @xmath89 ) , the inclination angle ( @xmath90 ) and the azimuth ( @xmath27 ) with respect to the @xmath88 axis .
the second reference system , @xmath21 , is obtained after a rotation of an angle @xmath27 around the @xmath1 axis , so that the angular velocity vector is now contained in the @xmath91 plane .
this second reference system will help us in obtaining the final analytical expressions .
one of the simplest non - trivial configuration of a stellar magnetic field is that of a dipole , whose axis is determined by a unit vector @xmath92 along the dipole moment and an intrinsic dipolar field .
the three components of the unit vector @xmath92 on the cartesian system @xmath21 are : @xmath93 where @xmath94 is the inclination of the dipolar moment with respect to the los and @xmath95 is its azimuth with respect to the @xmath2 axis . in the arbitrary reference system @xmath87 , the azimuth is transformed to @xmath96 .
the emergent stokes parameters @xmath97 and @xmath98 are then transformed to the arbitrary reference system following eq .
( [ eq : rotation_stokes ] ) .
note that our notation can be translated to that used by @xcite making the transformations ( in format ours@xmath99theirs ) : @xmath100 , @xmath101 , @xmath102 , @xmath103 .
therefore , the cosine of the inclination of the dipole axis with respect to the rotational axis , usually termed @xmath104 , is plainly given by the dot product @xmath105 , i.e. , @xmath106 .
the magnetic field at the surface of the star of a dipole placed at the center of the star can be expressed in cartesian coordinates in the @xmath21 system as : @xmath107 \nonumber \\ b_y&=&-\frac{b_p}{2}\left [ e_y - 3y(x e_x+y e_y+z e_z ) \right ] \nonumber \\ b_z&=&-\frac{b_p}{2}\left [ e_z - 3z(x e_x+y e_y+z e_z ) \right],\end{aligned}\ ] ] where @xmath108 is the intensity of the magnetic field at the poles .
an arbitrary point @xmath109 in the visible surface of the star can be written in terms of the polar coordinates on the plane of the sky @xmath29 as : @xmath110 where we have selected the positive solution for @xmath111 to isolate the visible half of the star .        for simplicity , we assume that the star is rotating as a solid body because it greatly simplifies the equations ( the introduction of a non - rigid rotation is formally identical , except for the computation of the @xmath71 integrals ) , so that : @xmath112 in terms of the polar coordinates on the stellar surface , the component of the velocity along the los in the @xmath21 reference system is given by : @xmath113    taking all these definitions into account and after some tedious algebra , it is possible to find relatively simple expressions for the @xmath71 integrals , which now depend on the set of parameters @xmath114 .
the procedure is to plug the functional form of the dipolar magnetic field and carry out the angular and radial integrals on the visible disk .
we end up with : @xmath115 \right.\nonumber \\ & + & b_1(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_3(n ) + a f_4(n ) + b f_5(n)\right ] \nonumber \\ & + & c_1(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_2(n+1 ) + a f_3(n+1 ) + b \left .
f_4(n+1)\right ] \right\ } \nonumber \\ x_2(i , j,\alpha,\omega , n , a , b ) & = & ( \omega \sin \alpha)^n \left\{a_2(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n+2 ) + a f_2(n+2 ) + b f_3(n+2)\right ] \right .
+ & b_2(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n+4 ) + a f_2(n+4 ) + b f_3(n+4)\right ] \nonumber \\ & + & c_2(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_2(n+1 ) + a f_3(n+1 ) + b f_4(n+1)\right ] \nonumber \\ & + & d_2(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_2(n+3 ) + a f_3(n+3 ) + b f_4(n+3)\right ] \nonumber \\ & + & \left .
e_2(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n ) + a f_2(n ) + b f_3(n)\right ] \right\ } \nonumber \\ x_3(i , j,\alpha,\omega , n , a , b ) & = & ( \omega \sin \alpha)^n \left\ { a_3(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n+2 ) + a f_2(n+2 ) + b f_3(n+2)\right ] \right . \nonumber \\ & + & b_3(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n+4 ) + a f_2(n+4 ) + b f_3(n+4)\right ] \nonumber \\ & + & c_3(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_2(n+1 ) + a f_3(n+1 ) + b f_4(n+1)\right ] \nonumber \\ & + & d_3(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_2(n+3 ) + a f_3(n+3 ) + b f_4(n+3)\right ] \nonumber \\ & + & \left .
e_3(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n ) + a f_2(n ) + b f_3(n)\right ] \right\ } \nonumber \\ x_4(i , j,\alpha,\omega , n , a , b ) & = & ( \omega \sin \alpha)^n a_4(i , j,\alpha,\omega , n ) \left [ ( 1-a - b ) f_1(n ) + a f_2(n ) + b f_3(n)\right ] .
\label{eq : x_functions}\end{aligned}\ ] ] these functions depend on the orientation of the dipole given by the angles @xmath94 and @xmath95 , the angular velocity parameters and the parameters of the clv and the index @xmath56 of the summation .
the following functions take into account the dependence on the orientation of the dipole and only depend on angular integrals : @xmath116 -6\sin^2 i \cos^2 j s(n ) + 6 s(n+2 ) \nonumber \\ b_2(i , j,\alpha,\omega , n ) & = & 9 \cos^2
i \left [ s(n+2 ) - s(n ) + 2s(n+4 ) \right ] \nonumber \\ & + & 9 \sin^2 i \cos^2 j \left [ s(n ) - 4 s(n+2 ) + 4s(n+4)\right ] \nonumber \\ & + & 9s(n+2 ) - 18s(n+4 )   \nonumber \\ c_2(i , j,\alpha,\omega , n ) & = & 6 \cos i \sin i \sin j s(n+1 ) \nonumber \\ d_2(i , j,\alpha,\omega , n ) & = & 18   \cos i \sin i \sin j \left [ s(n+1 ) - 2s(n+3 ) \right ]   \nonumber \\ e_2(i , j,\alpha,\omega , n ) & = & \sin^2 i \cos ( 2j ) s(n ) \nonumber \\ a_3(i , j,\alpha,\omega , n ) & = & -3   \sin^2 i \cos j \sin j s(n ) \nonumber \\ b_3(i , j,\alpha,\omega , n ) & = & 18   \sin^2 i \cos j \sin j \left [ s(n+2 ) - s(n+4 ) \right ]   \nonumber \\ c_3(i , j,\alpha,\omega , n ) & = & -3 \cos i \sin
i \cos j s(n+1 ) \nonumber \\ d_3(i , j,\alpha,\omega , n ) & = & 18 \cos i \sin i \cos j \left [ s(n+1 ) - s(n+3 ) \right ] \nonumber \\ e_3(i , j,\alpha,\omega , n ) & = & \sin^2 i \cos j \sin j s(n ) \nonumber \\ a_4(i , j,\alpha,\omega , n ) & = &   s(n ) , \label{eq : abcde_functions}\end{aligned}\ ] ] with the angular integral @xmath117 given by : @xmath118 \sqrt{\pi } \gamma \left ( \frac{1+n}{2 } \right ) \gamma^{-1 } \left ( 1+\frac{n}{2 } \right ) , \label{eq : s_functions}\ ] ] which is valid for @xmath119 .
likewise , the radial integrals @xmath120 are independent on the orientation of the dipole and are computed from : @xmath121 which are valid if @xmath122 .
once the expressions for @xmath71 have been particularized to the case of a rotating dipole , it is possible to analyze in detail the shape of the profiles according to eq .
( [ eq : fi_star ] ) .
we analyze the symmetry properties of the circular polarization profile , although this can be extended easily to the linear polarization profiles .
recall that the profile @xmath75 is a symmetric function of wavelength with respect to the rest wavelength of the line , while @xmath76 is antisymmetric .
according to the properties of the hermite polynomials , the terms with @xmath123 in the expression for @xmath124 in eq .
( [ eq : averages ] ) are symmetric , while those with @xmath125 are always antisymmetric .
likewise , for @xmath126 , the terms with @xmath123 are antisymmetric , while those with @xmath125 are symmetric . consequently , the polarized flux @xmath127 can be written as the addition of two functions , one symmetric and the other one antisymmetric : @xmath128 where @xmath129 \nonumber \\ f_{v}^\mathrm{antisym}(\lambda ) & = &   { c}\lambda_0 ^ 2 g \left [ \psi'(\lambda ) \overline{b_zf\mathcal{h}}_\mathrm{sym } + \psi(\lambda ) \overline{b_zf\mathcal{h}'}_\mathrm{antisym } \right].\end{aligned}\ ] ] taking into account the definition of the @xmath130 function of eq .
( [ eq : x_functions ] ) , it is easy to verify that the symmetric contribution to @xmath127 depends only on @xmath131 , while the antisymmetric contribution depends only on @xmath132
. consequently , we will find standard antisymmetric circular polarization profiles whenever @xmath133 , while they will be fully symmetric when @xmath134 .
stellar rotation is taken into account by keeping fixed the los and rotating the vector @xmath92 around the angular rotation velocity @xmath86 .
this is accomplished by computing the new rotated dipolar moment unit vector as : @xmath135 where the rotation matrix is given by : @xmath136,\end{aligned}\ ] ] with @xmath137 indicating the stellar rotation phase and taking values between @xmath138 and @xmath139 .
note that our initial of phases is an arbitrary position of the vectors , while using the standard notation , the initial phase coincides with the @xmath86 and @xmath92 vectors lying in the same plane parallel to the @xmath1 axis .
this implies a phase shift between our @xmath137 angle and the standard definition that is given by : @xmath140 finally , the inclination and azimuth angles with respect to the @xmath21 reference system are obtained using the standard relations , taking into account that @xmath141 is a unit vector : @xmath142 , \qquad j_\mathrm{rot } = \arctan \frac{(e_\mathrm{rot})_y}{(e_\mathrm{rot})_x}\ ] ]
all the functions needed in the previous equations are well behaved and can be computed numerically without much complexity .
the main difficulty resides on the fact that , if the rotation velocity @xmath143 is large , the convergence of the expansion in eqs .
( [ eq : averages_surface ] ) is slow and oscillating , thus making it hard to compute .
the main problem is that the hermite polynomials grow in amplitude while the coefficients decrease slow .
only after many terms ( the number of which depend on the ratio between the rotation velocity and the broadening ) , the coefficients start to compensate for the increase in the amplitude of the hermite polynomials and convergence is eventually obtained .
this is a well - known defect of the gram - charlier expansion ( e.g. , * ? ? ?
when the rotation velocity is much larger than the doppler width of the local spectral line ( larger than a factor @xmath144 ) , the compensation between the amplitude increase of the hermite polynomials and the coefficients of the series expansion occurs for very large values of @xmath56 [ see eqs .
( [ eq : averages_surface ] ) ] . at this moment
, the double precision computations start to become unprecise and the terms in the expansion do not compensate adequately .
this purely numerical drawback can be overcome if arbitrary precision methods are used or if more elaborate methods like the edgeworth expansion are adapted to our problem @xcite . in any case , we have been able to numerically compute expansions up to @xmath145 . in other words , our method with a double precision summation of the expansions shown in eqs .
( [ eq : averages_surface ] ) is presently applicable for rotation velocities not larger than 8 times the doppler width of the local spectral line .    in any case , it is advantageous to use the following set of modified hermite polynomials : @xmath146 which fulfill the following recursion formula : @xmath147 computing the hermite polynomials with this recursion formula , the expansions of eqs .
( [ eq : averages_surface ] ) can be extended to much larger values of @xmath56 .
the advantage of using the modified @xmath148 hermite polynomials is that their amplitude grow much slower with @xmath56 than @xmath55 due to the inclusion of the term @xmath149 and @xmath150 .
this facilitates the computation of each term in the expansions without numerical overflows .
we note that the hermite polynomials are computed using their recursion relations @xcite .    from the computational point of view , once the rotation angular velocity of the star is set , the wavelength dependence of eqs .
( [ eq : averages_surface ] ) and their derivatives in terms of the hermite polynomials can be precomputed once and used in all subsequent calculations .
this greatly accelerates the process because only the @xmath151 functions have to be computed for a given orientation of the dipolar moment and added until convergence .
a free computer program in idl is available to compute the observed stokes fluxes for a rotating dipole using the expressions in this paper . in this non - optimized version of the code , we synthesize profiles evaluated at 500 wavelength points with a step of 8 m  in @xmath1135 - 40 ms .
given our assumptions ( weak - field , weak - line , fixed center - to - limb variation , global dipolar field ) , one of the most direct applications of our equations is to study the magnetism of hot stars .
for this reason , we use typical values of the broadening and rotation velocity of such stars ( e.g. , * ? ? ? * ) .
we show in fig .
[ fig : examples ] two examples for a dipolar magnetic field with different field configurations .
both computations share the parameters for the clv ( @xmath152 and @xmath153 ) , the magnetic field strength in the pole ( @xmath154 g ) and the rotation velocity on the surface ( @xmath155 km s@xmath156 ) .
the artificial selected spectral line is centered at @xmath157   and has an effective land factor of @xmath9 and its equivalent for linear polarization of @xmath158 . since the magnetic field takes its maximum value at the poles , the spectral line can be safely considered in the weak - field regime .
the line broadening is set to 0.1  , typical of an observation at a spectral resolution of @xmath159 of a line with a microturbulent width of @xmath115 km s@xmath156 .
these values of the broadening and rotation velocity are close to the numerical limit that we presently have for computing the summations of eqs .
( [ eq : averages ] ) . given this limitation ,
our scheme is presently not applicable to fast rotators with narrow local absorption profiles .
the computed intensity flux normalized to the continuum flux is shown in the upper panel of each figure .
it is independent of the phase since , under the weak field regime , the intensity does not depend on the magnetic field and we have assumed that the star is structureless .
however , it takes into account the rotation broadening and the presence of a clv .
concerning the linear polarization fluxes , we see maximum amplitudes around 0.02% , while this increases to values around 100 times larger for circular polarization .
this is a direct consequence of the fact that linear polarization appears only in second order .
we have generalized the expressions for the stokes parameters in the weak - field weak - line regime to the stellar case in which rotation is taken into account .
this produces an intensity flux profile that is broadened by rotation while it is insensitive to the magnetic field . as a consequence , our formalism
can only be applied to stars in which the intensity profile is not observed to be modulated by rotation .
however , the circular and linear polarized fluxes are sensitive to the topology of the magnetic field .
we have demonstrated that the inclusion of doppler shifts generates a dependence of circular polarization of the intensity flux itself , and not only on its wavelength derivative . likewise , the linear polarization profiles depend now on the first and second wavelength derivative and on the intensity flux itself .
we have verified also that the expressions developed in this paper simplify to those obtained by @xcite in the non - rotating dipole case .
the general weak - field expressions have been particularized to a rotating dipole .
the extension to a more general multipolar magnetic field is relatively easy and can be done fast if analytical expressions for the angular and radial integrals exist . in subsequent papers of this series
we will deal with thin magnetized disks in keplerian rotation and unresolved structures in the quiet sun .
likewise , we will also explore the possibility of using more elaborate expansions to improve the convergence .
we thank r. manso sainz for fruitful discussions .
financial support by the spanish ministry of economy and competitiveness through projects aya2010 - 18029 ( solar magnetism and astrophysical spectropolarimetry ) and consolider - ingenio 2010 csd2009 - 00038 is gratefully acknowledged .                                      ,
v. , del toro iniesta , j.  c. , lvarez - herrero , a. , domingo , v. , bonet , j.  a. , gonzlez fernndez , l. , lpez jimnez , a. , pastor , c. , gasent blesa , j.  l. , mellado , p. , piqueras , j. , aparicio , b. , balaguer , m. , ballesteros , e. , belenguer , t. , bellot rubio , l.  r. , berkefeld , t. , collados , m. , deutsch , w. , feller , a. , girela , f. , grauf , b. , heredero , r.  l. , herranz , m. , jernimo , j.  m. , laguna , h. , meller , r. , menndez , m. , morales , r. , orozco surez , d. , ramos , g. , reina , m. , ramos , j.  l. , rodrguez , p. , snchez , a. , uribe - patarroyo , n. , barthol , p. , gandorfer , a. , knoelker , m. , schmidt , w. , solanki , s.  k. , & vargas domnguez , s. 2011 , sol .
phys . , 268 , 57            , s.  k. , barthol , p. , danilovic , s. , feller , a. , gandorfer , a. , hirzberger , j. , riethmller , t.  l. , schssler , m. , bonet , j.  a. , martnez pillet , v. , del toro iniesta , j.  c. , domingo , v. , palacios , j. , knlker , m. , bello gonzlez , n. , berkefeld , t. , franz , m. , schmidt , w. , & title , a.  m. 2010 , apjl , 723 , l127      , j. 2011 , in astronomical society of the pacific conference series , vol . 437 , solar polarization 6 , ed . j.  r.  kuhn , d.  m.  harrington , h.  lin , s.  v.  berdyugina , j.  trujillo - bueno , s.  l.  keil , & t.  rimmele , 83


in this paper we are concerned with the study of the three - dimensional map @xmath6 defined by @xmath7 with three parameters , @xmath8 , @xmath9 , and @xmath10 .
this map is a natural generalization of the famous two - dimensional hnon map @xcite ; indeed , like the latter , the map ( [ eq:3hm1 ] ) is quadratic , has constant jacobian @xmath11 and moreover , when @xmath12 , reduces to the two - dimensional hnon map .
therefore , it is natural to call ( [ eq:3hm1 ] ) the  _
3d hnon map_. " we will study a number of properties of ( [ eq:3hm1 ] ) , including :    1 .
its nontrivial dynamics and , most importantly , the existence of _ wild - hyperbolic _ strange attractors ; 2 .   its origination in homoclinic bifurcations ; 3 .
its structure as a member of the family of three - dimensional quadratic maps with constant jacobian .
three - dimensional maps have not been as widely studied as one- and two - dimensional maps , though in recent years there has been a considerable increase in their study .
one of important reasons for this fact is that multi - dimensional dynamical systems ( the dimension of the phase space is at least four for flows and three for maps ) can exhibit complicated dynamics that is decidedly distinct from lower - dimensional cases . in particular they can possess a new variety of strange attractor called a _ wild - hyperbolic attractor _ @xcite .
wild - hyperbolic attractors have no uniform hyperbolic structure ; moreover , they may contain homoclinic tangencies and coexisting periodic orbits of different types .
however , these attractors contain no stable invariant subsets such as attracting periodic orbits .
most crucially , the latter property holds for open regions in the space of smooth dynamical systems .
thus , wild - hyperbolic attractors are very interesting objects and their detection in concrete models is a fundamental problem of nonlinear dynamics .
most of the strange attractors occurring in two - dimensional , smooth maps ( or three - dimensional flows ) are not structurally stable since they possess homoclinic tangencies .
for such systems , even though the orbits are observed to be quite chaotic for some parameter values , very small changes in the parameters may destroy the aperiodicity and give rise to periodic behavior .
such parameter regions are called _ stability windows_. even if these windows are not  observable `` in numerical computations , their appearance ( under arbitrarily small perturbations ) is inevitable because of the prevalence of homoclinic tangencies .
such attractors were called _
quasi - attractors _ in @xcite .
one of the most challenging problems in dynamical systems is to distinguish aperiodic motion from extremely long periodic motion .
moreover , since most of the chaotic attractors we see in applications are definitely quasi - attractors , it is unclear in which sense they are chaotic , or how one can define and measure the ' ' probability " for a quasi - attractor to be chaotic .
however , it is well - known that there are strange attractors that are free of these problems ; for example , hyperbolic attractors and the lorenz attractor . though the latter is not structurally stable @xcite,@xcite ( in contrast to hyperbolic attractors ) it contains no stable periodic orbits , every orbit in it has positive maximal lyapunov exponent , etc .
moreover , these properties are robust @xcite,@xcite .
the reason is that lorenz attractor possesses a pseudo - hyperbolic structure . in terms of turaev and shilnikov
@xcite this means that the following properties hold :    1 .
there is a direction in which the flow is strongly contracting ( `` strongly '' means that any possible contraction in transverse directions is always strictly weaker ) ; and 2 .
transverse to this direction the flow expands areas .
the robustness of pseudo - hyperbolicity ensures the robustness of the chaotic behavior of orbits in the lorenz attractor .
moreover , pseudo - hyperbolicity is maintained even for ( small , indeed ) time - periodic perturbations so that a periodically forced lorenz attractor is an example of a genuinely  strange " attractor since it contains no stable invariant subsets ( for more discussion see @xcite,@xcite,@xcite ) .
however it is necessary to note that homoclinic tangencies arise inevitably under such perturbations .
consequently systems with periodically forced lorenz attractors fall into newhouse regions and demonstrate an extremely rich dynamics @xcite,@xcite .
such attractors were called wild - hyperbolic attractors in @xcite ; in this paper the theory of wild - pseudo - hyperbolic attractor was discussed and an example of a wild - hyperbolic spiral attractor was constructed .
generalizing these ideas , we say that wild - hyperbolic attractors possess the following two main properties that distinguish them from other attractors : ( 1 ) wild - hyperbolic attractors allow homoclinic tangencies ( hence , they belong to newhouse regions ) ; ( 2 ) every such attractor and all nearby attractors ( in the @xmath13-topology with @xmath14 ) have no stable periodic orbits .
one of the most felicitous properties of wild - hyperbolic attractors is that they can be created in _
local bifurcations _ of periodic orbits ; in particular , when the orbit has three or more floquet multipliers on the unit circle @xcite .
consequently , the construction of wild - hyperbolic attractors in concrete models is possible when the model contains sufficiently many parameters to permit such a degeneracy .
indeed , the map ( [ eq:3hm1 ] ) is sufficiently complex and has been shown to possess wild - hyperbolic attractors @xcite . moreover , numerical simulations have shown that these attractors persist in quite big parameter domains and can be order - one in size @xcite .    in  [ sec :
attractors ] we will show that ( [ eq:3hm1 ] ) has a small , wild - hyperbolic lorenz attractor in a ( small ) parameter domain near the point @xmath15 at which a fixed point with multipliers @xmath16 exists . the general case of fixed points with such a triplet of multipliers was considered in @xcite , where it was shown that two distinct cases occur .
one of these cases has a local flow normal form that coincides with the shimitsu - marioka system ( see ( [ = sm ] ) ) .
this system has a lorenz attractor in some parameter domains @xcite,@xcite .
as was noted in @xcite a periodically forced lorenz system ( when the perturbation is rather small ) can give rise to a wild - hyperbolic attractor .
moreover , this does not destroy the pseudo - hyperbolicity of the attractor @xcite,@xcite .
note that the flow normal form only approximates the map and , since the residual terms are exponentially small in time , it is reasonable that the map will have a wild - hyperbolic attractor ( which is small , in fact , because the approximation must be very good ) , see more discussions in @xcite .
the second property of ( [ eq:3hm1 ] ) that we mentioned above is that this map can arise as a local normal form for multi - dimensional systems with a homoclinic tangency .
this shows that such maps are not exotic and gives a reason for an extensive study of their properties .
it is well - known that quadratic maps appear naturally when studying bifurcations of quadratic homoclinic tangencies @xcite,@xcite,@xcite,@xcite,@xcite .
they appear as normal forms of rescaled first return maps to the neighborhood of a point on the homoclinic orbit .
there are four standard forms for such one- and two - dimensional  homoclinic maps : "    1 .
@xmath17 ( the logistic map ) ; 2 .
@xmath18 ( the standard hnon map ) ; 3 .
@xmath19 ( the mira map ) ; 4 .
@xmath20 ( a generalized hnon map ) .
the third case is a two - dimensional endomorphism and is one of the well - studied polynomial maps ( see , e.g. @xcite ) ; we call it the  mira map " following @xcite )
. the fourth case can be obtained as a small perturbation ( due to @xmath21-terms ) of the standard hnon map .
this map was obtained in a study of bifurcations of two - dimensional diffeomorphisms with a homoclinic tangency to a saddle fixed point ( with unit jacobian ) @xcite,@xcite .
the 3d hnon map ( [ eq:3hm1 ] ) is absent from this list , and in  [ sec : local][sec : rescaling ] we will fill this gap .    in  [ sec : local ] we study bifurcations of a three - dimensional diffeomorphism with quadratic homoclinic tangency to a saddle - focus fixed point . for a three - dimensional map , such points can be of two types : @xmath22 or @xmath23 ( sometimes called _ type a _ and _ type b _ , respectively ) .
recall that @xmath24 saddle is a fixed point with @xmath25 , @xmath26 , and @xmath27 . in our case the two - dimensional manifold will be assumed to correspond to a complex conjugate pair of multipliers .
we will consider the case where the the product the absolute value of the multipliers , @xmath28 is one .
the bifurcations in the general case , where @xmath29 , were studied in @xcite,@xcite .
we will use their  rescaling method " to study our case as well .
the essence of this method is that a suitable change of coordinates and parameters in a neighborhood of the homoclinic orbit is found to bring the map into a standard form .
if the homoclinic tangency is assumed to be quadratic the resulting map is also typically quadratic .
for example , for a type @xmath22 saddle - focus with @xmath30 , the rescaled first return maps coincide with the standard hnon map up to asymptotically small terms ( when the time of return tend to infinity ) . for a saddle - focus of type @xmath23 ,
the corresponding rescaled first return map is the mira map .
note the condition @xmath31 leads to maps that are effectively two - dimensional ; consequently all three - dimensional volumes near the fixed point are collapsed .
this is not the case when @xmath32 for the unperturbed map when three - dimensional rescaled maps can appear a priori .
indeed , we will show that when @xmath33 the map ( [ eq:3hm1 ] ) occurs as the normal form of the rescaled map for type @xmath22 saddle - focus ; and for the type @xmath23 case the corresponding normal form is ) is also the rescaled first return map for a four - dimensional diffeomorphism that has a quadratic homoclinic tangency to a saddle - focus fixed point with two pairs of complex conjugate leading multipliers @xcite,@xcite,@xcite . ]
@xmath34    both ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) are quadratic three - dimensional maps with constant jacobian @xmath35 .
it is natural to call them the three - dimensional hnon maps of the _ first _ and _ second kind _ , respectively .
the quadratic 3d hnon maps are elements of the group of polynomial automorphisms .
the study of the dynamics of polynomial mappings has a long history applied dynamics ; for example , they often used in the study of particle accelerators @xcite .
the study of _ cremona maps _
, i.e. , polynomial maps with constant jacobian , also has intrinsic mathematical interest @xcite .
an interesting mathematical problem is to obtain a normal form for an arbitrary degree-@xmath36 cremona map .
this problem is unsolved ; the major obstruction is that it is not known in general that such a map has a polynomial inverse .
this is the content of the much studied ,  jacobian conjecture " :    let @xmath37 be a cremona map
. then @xmath38 is bijective and has a polynomial inverse .
this conjecture is still open .
the set of polynomial maps with polynomial inverses is called the _ affine cremona group_. for the case of the plane , the structure of this group was obtained by jung @xcite .
jung s theorem states that planar polynomial diffeomorphisms are `` tame ; '' that is , they can be written as a finite composition of `` elementary '' ( or triangular ) and affine maps .
more recently friedland and milnor @xcite showed that any map in this group is either conjugate to a composition of generalized hnon maps , @xmath39 where @xmath40 is a polynomial and @xmath41 , or else the map is dynamically trivial .
not much is known about the structure of the affine cremona group for higher dimensions , though it was recently proved that the three - dimensional cremona map @xmath42 is not tame @xcite .
thus the normal form for the 3d affine cremona group must consist of more than tame maps .
however , it is still possible that the nontame maps are dynamically trivial , like the above example .
another interesting problem concerns the approximation of smooth diffeomorphisms by polynomial diffeomorphisms .
one result along this direction is for the class of @xmath43 symplectic diffeomorphisms , which , as turaev has shown @xcite , can be approximated by compositions of the symplectic cremona maps @xmath44 where @xmath45 and @xmath46 is a polynomial .
related approximation methods are used in the study of particle accelerators @xcite , and results have also been obtained for analytic symplectic maps @xcite .    some work has been done on finding normal forms for quadratic symplectic and quadratic volume - preserving maps .
moser showed that every quadratic symplectic map in @xmath47 can be written as the composition of an affine symplectic map and a symplectic shear @xcite .
for the four - dimensional case , he obtained a normal form .
for the case of volume - preserving maps , it was shown in @xcite,@xcite that any quadratic three - dimensional , volume - preserving map whose inverse is also quadratic and whose dynamics is non - trivial can be transformed ( by affine transformations ) to the form @xmath48 where @xmath49 , @xmath50 and one of the parameters @xmath51 or @xmath52 can be eliminated ( by a coordinate shift ) . by analogy with the planar results , there are two additional cases but they have trivial dynamics ( see  [ sec : quadratic ] ) .
we also note that it is not hard to construct three - dimensional quadratic diffeomorphisms that do not fall into this classification because they have an inverse of higher degree .
for example , the diffeomorphism @xmath53 has a quartic inverse .
note that both of the maps ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) fall into the form ( [ eq:3hm3 ] ) .
thus the results of @xcite show that when @xmath54 both of the 3d hnon maps occur as normal forms for quadratic elements of the affine cremona group . moreover , even when @xmath55 the map ( [ eq:3hm2 ] ) is a normal form for a certain class of quadratic maps @xcite .
tatjer showed that every quadratic map , @xmath38 , for which @xmath56 is also quadratic , but has no invariant linear foliation is linearly conjugate to ( [ eq:3hm2 ] ) . in
[ sec : statement ] , we will generalize the results of @xcite to show that ( [ eq:3hm3 ] ) is the the normal form for quadratic diffeomorphisms with quadratic inverses when @xmath55 as well .
consider a three - dimensional diffeomorphism @xmath57 for some @xmath58 that satisfies the following conditions ( see fig .
[ fig : type21 ] ) :    1 .
@xmath0 has a type @xmath22 saddle - focus fixed point @xmath59 with multipliers @xmath60 where @xmath61 ( w.l.o.g .
@xmath62 ) , ; 2 .   @xmath63 .
3 .   the stable , @xmath64 , and unstable , @xmath65 , invariant manifolds of the fixed point @xmath59 have a quadratic homoclinic tangency at the points of a homoclinic orbit @xmath66 .
saddle - focus for the map @xmath0 under assumputions ( a)-(c ) .
[ fig : type21],width=288 ]    diffeomorphisms that are close to @xmath0 and satisfy conditions ( b ) and ( c ) comprise , in the space of @xmath13 diffeomorphisms , a codimension - two , locally - connected surface @xmath67the bifurcation surface .
our main problem is to study bifurcations of periodic orbits in families of diffeomorphisms that transversely cross @xmath67 .
we denote such a family by @xmath4 so that @xmath0 satisfies ( a)-(c ) .
the first question that arises is the optimal choice of the governing parameters @xmath68 .
it is evident that the set of such parameters must include both a parameter @xmath69 of the splitting of the initial quadratic homoclinic tangency and a parameter @xmath70 , similar to @xmath71 , that varies the jacobian of the saddle - focus fixed point .
however , it was shown in @xcite that diffeomorphisms on @xmath67 possess continuous ( topological ) conjugacy invariants on the set of non - wandering orbits , so - called @xmath72-moduli .
the most important modulus is the angular argument @xmath73 of the complex multipliers of the saddle - focus ( see  [ sec : rescaling ] ) . by definition ,
@xmath72-moduli are natural governing parameters .
therefore , we include the @xmath72-modulus @xmath73 into the set of governing parameters .
thus , we take as the third parameter @xmath74 .
more generally , we must consider a three - parameter family @xmath75 , @xmath76 , that is transversal to @xmath67 at @xmath5 and such that @xmath77    consider a sufficiently small and fixed neighbourhood @xmath78 of the orbit @xmath66 .
it can be represented as a union of a small neighbourhood @xmath79 of @xmath59 with a number of small neighbourhoods of those points of @xmath66 which do not belong to @xmath79 . in the family @xmath75
we study bifurcations of _ single - round periodic orbits _ from @xmath78that is orbits periodic that enter @xmath79 exactly once .
every point of such an orbit of period @xmath80 can be considered as a fixed point of the corresponding first return map @xmath81 along orbits lying in @xmath78 . in  [ sec : rescaling ] , we will prove the following result about the form of this first return map in some domains of the parameters .
[ tm:1 ] let @xmath75 be the three parameter family defined above .
then , in any neighborhood of the origin , @xmath5 , of the space of parameters , there exist infinitely many open domains , @xmath82 , such that the first return map @xmath81 takes the form @xmath83 where @xmath84 are rescaled coordinates and the domains of definition of the new coordinates and the parameters @xmath85 are arbitrarily large and cover , as @xmath86 , all finite values . furthermore
, the jacobian , @xmath10 , is positive or negative depending on the orientability of the map @xmath81 .
the terms denoted @xmath87 are functions ( of the coordinates and parameters ) that tend to zero as @xmath88 along with all derivatives up to order @xmath89 .
theorem [ tm:1 ] shows that the 3d hnon map ( [ eq:3hm1 ] ) is a normal form for this homoclinic bifurcation , since it appears when we omit the @xmath87-terms in ( [ eq : t07 + 1 ] ) .
we will see that the new parameters @xmath90 and @xmath10 can be given as functions of @xmath91 by @xmath92 where @xmath93 is some coefficient ( an invariant of @xmath0 ) and @xmath94 .
the exact relations will be given in formulas ( [ mpar * ] ) , ( [ eq : t04 ] ) and ( [ eq : t06 ] ) below .
thus , we can study the bifurcations of single round periodic orbits in the family @xmath75 , when @xmath95 , by studying the bifurcations of fixed points of the 3d hnon map ( [ eq:3hm1 ] ) .
certainly , the problem of studying bifurcations in the 3d hnon map is much simpler than the corresponding problem for the original first return maps .
moreover , the 3d hnon map provides a simple , standard model for which various well - tried methods ( including numerical ones ) can be applied .
however , the bifurcation problem as whole for the 3d hnon map is very complicated and here we focus only on two aspects . in  [ sec : attractors ] we first give the simple construction of bifurcation surfaces , lines and points corresponding , respectively , to codimension one , two and three bifurcations of the fixed points .
then we consider the problem of the existence of wild - hyperbolic strange attractors .
we will obtain the following result .
[ tm:2 ] let @xmath75 be the three parameter general family defined as before .
then , in any neighborhood of the origin @xmath5 of the space of parameters , there exist infinitely many open domains @xmath96 such that the diffeomorphism @xmath75 of the family at @xmath97 has a wild - hyperbolic strange attractor of lorenz - type .
the proof of this theorem is given in
[ sec : attractors ] following the analysis of @xcite
. the main idea of the proof is to apply the results of paper @xcite where it was shown that maps having a fixed point with the triplet @xmath98 of multipliers can , in some cases , be reduced to a flow normal form that coincides with the shimitsu - marioka system ( see ( [ fl1 ] ) ) which possesses a lorenz attractor for some parameter domains @xcite,@xcite .
our verification of these facts for the corresponding fixed point of the 3d hnon map allows us to obtain th .
[ tm:2 ] . a justification of the fact that the attractor in the 3d hnon map is actually the wild - hyperbolic lorenz - type attractor is given in @xcite .
consequently , we will omit these details .
solving the bifurcation problem for the family @xmath75 allows us to solve almost automatically an analogous problem for a diffeomorphism @xmath99 ( @xmath100 with a type @xmath23 saddle focus .
specifically , assume that @xmath101 satisfies ( see fig .
[ fig : type12 ] ) :    1 .
@xmath101 has a type @xmath23 fixed point @xmath102 with multipliers @xmath103 where @xmath104 ( w.l.o.g .
@xmath105 ) ; 2 .   @xmath106 ; and 3 .
the stable and unstable invariant manifolds of @xmath102 have a quadratic homoclinic tangency at the points of some homoclinic orbit @xmath107 .
saddle - focus for the map @xmath101 under assumptions ( a@xmath108)-(c@xmath108 ) .
[ fig : type12],width=288 ]    as in the case of the type @xmath22 saddle - focus , we consider a three - parameter family @xmath109 where @xmath110 .
we suppose here that @xmath111 is a splitting parameter , @xmath112 .
[ tm:3 ] let @xmath109 be the three - parameter general family defined above .
then , in any neighborhood of @xmath113 there exist infinitely many open domains , @xmath114 , such that the first return map @xmath115 takes the form @xmath116 where @xmath84 are rescaled coordinates and the domains of definition of the new coordinates and the parameters @xmath85 are arbitrarily large and cover , as @xmath86 , all finite values .
furthermore , the jacobian , @xmath10 , is positive or negative depending on the orientability of the map @xmath117.the terms denoted @xmath87 are functions ( of the coordinates and parameters ) that tend to zero as @xmath88 along with all derivatives up to order @xmath89 .
the proof immediate from th .
indeed , the inverse @xmath118 is , evidently , a diffeomorphism satisfying the original conditions ( a(c ) .
thus , the map ( [ eq : t07 + 1 ] ) is , in this case , the rescaled map for the inverse first return map @xmath119 .
the last can be be easily written in form ( [ eq : t07 + 2 ] ) , as we show explicitly below .
the inverse of ( [ eq:3hm1 ] ) is of interest in any case . since ( [ eq:3hm1 ] ) has the constant jacobian @xmath10 it is invertible whenever @xmath120 .
this inverse is @xmath121 this map can be rewritten in the `` standard '' form of @xcite , if we define new variables @xmath122 and new parameters @xmath123 then the resulting map has precisely the form ( [ eq:3hm2 ] ) . as we stated above map
is well - known in homoclinic dynamics : it was found to be the rescaled first return map for four - dimensional diffeomorphisms that have a quadratic homoclinic tangency to a saddle - focus fixed point with two pairs of complex conjugate leading multipliers @xcite,@xcite,@xcite .
the dynamics of each of the maps ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) is extremely rich .
however , for the volume - contracting case , we note that these maps demonstrate distinct types of stable dynamics . for example , ( [ eq:3hm1 ] ) gives rise to attractors similar to those of the lorenz system , as we will see in  [ sec : attractors ] .
however , ( [ eq:3hm2 ] ) need not have attractors of this  wild " type , although certainly , a `` lorenz repeller '' will be present .
nevertheless , we can expect existence of  extremely small " lorenz - type attractors in ( [ eq:3hm2 ] ) by the following argument . as we will see in  [ sec : attractors ] , the map ( [ eq:3hm2 ] ) can have a fixed point with multipliers @xmath124 at certain values of its parameters .
the flow normal form near this point can have a homoclinic loop to a saddle - focus and , thus , ( [ eq:3hm2 ] ) can have a fixed point @xmath59 with multipliers @xmath125 , where @xmath126 and @xmath127 , such that the manifolds @xmath128 and @xmath129 will have a ( quadratic ) homoclinic tangency . consequently , we can apply th .  [ tm:2 ] to this situation giving immediately the promised  extremely small " attractors .
nevertheless , as it is shown in @xcite the lorenz - type wild - hyperbolic attractor for ( [ eq:3hm1 ] ) can be quite big and were observed to have a large basin .
the proof of existence of big , globally attracting wild - hyperbolic attractors in map ( [ eq:3hm2 ] ) is still open problem .    finally , we consider the more general structure of polynomial diffeomorphisms of @xmath130 . as we remarked in the introduction ,
the 3d hnon maps ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) have the form of ( [ eq:3hm3 ] ) when @xmath54 . as we will shown in  [ sec : quadratic ] , this result has the following generalization ;    [ tm:4 ] let @xmath131 be a quadratic diffeomorphism with a quadratic inverse and constant jacobian @xmath132 .
then , after an affine coordinate transformation @xmath38 can be written in the form ( [ eq:3hm3 ] ) or else the dynamics of @xmath38 is trivial .
note that both of the maps ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) fall into this classification .
the classification of the dynamics of the trivial cases can be done from the normal forms ( [ eq : qforms ] ) that will be obtained in  [ sec : quadratic ] .
we begin the proof of th .
[ tm:1 ] by constructing a local map @xmath133 , valid in a neighborhood @xmath79 of the fixed point @xmath134 , that simplifies the description of the local unstable @xmath135 and local stable @xmath136 manifolds of @xmath59 .
then we obtain a global map @xmath137 that takes a neighborhood of a point on @xmath138 to a neighborhood of a point of @xmath139 .    by @xcite,@xcite
, one can introduce @xmath140 coordinates @xmath141 in a neighborhood @xmath78 of the fixed point @xmath59 such that the local map @xmath142 can be written , for small enough @xmath68 , in the form ) holds for general case , including @xmath33 .
for the volume - preserving case the analogous form was established in @xcite ] @xmath143 where @xmath144 is the rotation by angle @xmath73 , and @xmath145 here the coefficients , @xmath146 , and @xmath147 as well as @xmath148 and @xmath149 depend ( smoothly ) on the parameters @xmath68 . the form ( [ eq : t0 ] ) will be called the _ main normal form _ ; it is very convenient for calculations . in particular ,
the local stable and unstable manifolds of @xmath59 are @xmath150 and @xmath151 in the neighborhood @xmath78 .
when @xmath152 , choose a homoclinic point @xmath153 on the homoclinic orbit @xmath66 and a corresponding point @xmath154 that is an image of @xmath155 : @xmath156 for some @xmath157 , recall fig .
[ fig : type21 ] .
let @xmath158 and @xmath159 be some small neighborhoods of the points @xmath160 and @xmath155 .
the global map @xmath161 is well defined for small enough @xmath68 , and can be written in the coordinates of ( [ eq : t0 ] ) as @xmath162 where the @xmath163 matrix @xmath164 , the vectors @xmath165 , @xmath166 , and @xmath167 , the coefficients @xmath168 and @xmath36 , and the higher order terms all depend smoothly on @xmath68 .
we noted especially the dependence of the coefficient @xmath169 that will become the splitting parameter for the manifolds @xmath64 and @xmath65 near the homoclinic point @xmath160 .
as discussed in the previous section , we therefore denote it by @xmath170 and consider as the first of the governing parameters .    the assumption ( c ) that the tangency between @xmath65 and @xmath64 is quadratic implies that @xmath171 in addition , since @xmath0 is assumed to be a diffeomorphism , and since the jacobian @xmath172then we must have @xmath173 .
indeed if @xmath174 is not identically zero , then one can use the linear coordinate change , @xmath175 where @xmath176 , ( [ eq : rotation ] ) , is a rotation by an angle @xmath177 , to transform the map ( [ eq : t1 ] ) to a new map of the same form such that the vector @xmath166 becomes @xmath178 here the new matrix @xmath164 , point @xmath165 , and vector @xmath167 are modified from the corresponding coefficients in the original ( [ eq : t1 ] ) .
now the nonvanishing of the determinant of the jacobian ( [ eq : j1 ] ) implies that @xmath179    one of merits of the main normal form ( [ eq : t0 ] ) is that the iterations @xmath180 for any @xmath80 can be calculated quite easily , namely , in a form close to the form when @xmath133 is linear .
in particular for the linear case we have @xmath181 this is the so - called _ cross - form _ of map @xmath182 .
we can write @xmath182 in an analogous cross - form when @xmath133 is in the form ( [ eq : t0 ] )    [ lem2 ] for any positive integer @xmath80 and for any sufficiently small @xmath68 the map @xmath183 can be written in the cross - form @xmath184 where the vector @xmath185 and the function @xmath186 are uniformly bounded along with all derivatives up to order @xmath187 . in addition , all of the derivatives of order @xmath188 for the functions @xmath189 and @xmath190 tend to @xmath191 as @xmath192  .
let @xmath193 denote the orbit : @xmath194 .
then it is easy to see that the @xmath195 iterate of ( [ eq : t0 ] ) is @xmath196    the solution of ( [ eq:3 * ] ) for any given sufficiently small @xmath197 and large @xmath80 ( the so - called boundary value problem ) can be found by successive approximations @xcite if the zeroth approximation is taken as the corresponding solution of the linear problem . using the fact that @xmath198 , we obtain the following estimates @xmath199 where @xmath200 and @xmath201 are some positive constants estimating the norms for derivatives of the functions @xmath148 and @xmath149 in the domain @xmath78 , from ( [ eq : t0 ] ) .
these inequalities evidently imply the desired result for the map .
estimates for the derivatives can be found analogously .
general results of this form for near saddle fixed point were given in @xcite,@xcite .
here we give sharper estimates that apply in this particular case .
the main task of this section is to obtain an analytical form for the first return map for diffeomorphisms of the family @xmath75 to a neighborhood of the homoclinic point @xmath160 . to construct this map we first find an image of the local map @xmath133 for which @xmath202
this will be true for each @xmath203 where @xmath204 is some sufficiently large integer .
for each such @xmath80 there are strips @xmath205 and @xmath206 such that @xmath207 .
a first return map is obtained by composing @xmath182 with the global map @xmath208 , to obtain a map on a neighborhood of @xmath160 , recall fig .  [
fig : type21 ] .
thus the first return map will be of the form @xmath209 under the global map @xmath137 the strip @xmath210 is transformed into three - dimensional horse - shaped region @xmath211 that has non - empty intersection with @xmath210 when @xmath68 is small enough .
thus , the map @xmath212 is a generalized , three - dimensional , horseshoe map .
we will show that these horseshoe maps can possess rich dynamics ; for example , they can lorenz - like , wild - hyperbolic strange attractors .
the goal is the prove the following lemma :    [ rescaling lemma . ]
[ lm : rescaling ] for any sufficiently large @xmath80 such that @xmath213 the first return map @xmath212 can be brought , by means of an affine transformation of coordinates , to the following form @xmath214 where @xmath215 are some coefficients , @xmath216 .
\end{array}\ ] ] also , in ( [ eq : t07 ] ) the @xmath217-terms have the indicated asymptotic behavior for all derivatives up to order @xmath187 .    using formulas ( [ eq : t1 ] ) , ( [ eq : t0kk ] ) and conditions ( [ eq : van3 ] ) and ( [ eq : b1 ] ) we can write the first return map @xmath218 in the form @xmath219 where the @xmath163 matrix @xmath220 and @xmath221 introduce new coordinates @xmath222 where @xmath223 , in order to eliminate the affine terms in the first two equations and the term linear in @xmath224 in the third equation of ( [ eq : van8 ] ) . then , the map ( [ eq : van8 ] ) is rewritten as @xmath225 where @xmath226 and @xmath227 .
\end{array}\ ] ]    next , we rescale the map ( [ van9 ] ) introducing new coordinates ( x , y , z ) defined through @xmath228 where we suppose that values of @xmath73 are such that ( [ eq : t02 ] ) holds . then , map ( [ van9 ] ) is rewritten in the new coordinates in form ( [ eq : t07 ] ) .
now theorem  [ tm:1 ] easily follows from this lemma .
indeed , we have only to indicate those domains @xmath82 of the initial parameters where the coefficients in the third equation of ( [ eq : t07 ] ) are finite . denote the coefficient of the third equation of ( [ eq : t07 ] ) in front of @xmath229 as @xmath9 , i.e. @xmath230 since @xmath231 and @xmath232 , the coefficients
@xmath8 and @xmath9 can take arbitrary finite values when varying parameters @xmath233 ( for @xmath8 ) and @xmath73 ( for @xmath9 ) , see formulas ( [ mpar * ] ) and ( [ eq : t04 ] ) .
moreover , the finiteness of @xmath9 means that the corresponding values of the trigonometrical term @xmath234 are asymptotically small and , in any case , do not exceed @xmath235 in order . for such values of @xmath73
, we obtain @xmath236 thus , by ( [ van7 ] ) , @xmath237 in this case .
denote as @xmath10 the coefficient in front of @xmath238 from the third equation of ( [ eq : t07 ] ) . by virtue of ( [ eq : t05 ] ) and ( [ van7 ] ) we have @xmath239 since @xmath240 at @xmath5
, it follows from ( [ eq : t06 ] ) that @xmath10 can take ( when @xmath80 is large ) any arbitrary finite value ( positive or negative depending on the sign of @xmath241 ) as the value of the jacobian @xmath242 of the local map @xmath133 varies .
this completes the proof of theorem  [ tm:1 ] .      1 .
we noted that at least three parameters must be used to study bifurcations of a homoclinic orbit for a three - dimensional saddle - focus .
the governing parameters that we use are @xmath243 corresponding to the splitting distance between the stable and unstable manifolds , the phase @xmath244 of the focus multiplier , and the jacobian @xmath245 of the saddle - focus .
if the map @xmath0 is exactly volume preserving , i.e. , @xmath246 , then the resulting map ( [ eq : t07 ] ) is as well , i.e. @xmath247 .
such maps were studied in @xcite . moreover , th .
[ tm:2 ] implies that there are nearby maps ( by unfolding the parameter @xmath10 to make the map volume contracting ) that have strange attractors .
in this section we will study the dynamics of the 3d hnon map ( [ eq:3hm1 ] ) and prove th .
the dynamics can be complex , especially when @xmath248 is close to one , essentially because when @xmath249 the map has , for certain values of the parameters , fixed points with three multipliers are on the unit circle .
it is well known that bifurcations of such points can lead to appearance of homoclinic tangencies , invariant circles , and even strange attractors ( as @xmath248 decreases below one ) or repellers ( as @xmath248 increases above one ) @xcite,@xcite .
we are interested primarily in the stable dynamics of the map ( [ eq:3hm1 ] ) , i.e. , in characterizing its attractors .
one way to find attractors with complex dynamics is to study the local normal forms that arise near fixed points that have three multipliers that are @xmath250 .
we follow here the analysis of @xcite where local normal forms of the equilibrium of a flow with three zero eigenvalues and with additional symmetries were analyzed and conditions for the existence of strange attractors were given .
these results can be applied to the analysis of bifurcations of codimension - three fixed points of three - dimensional maps . in this case
we can often construct some approximated flow normal form whose time - one shift map coincides ( up to terms of a certain order ) with some power of the original map restricted on a neighborhood of the fixed point .
usually , the presence of multipliers @xmath251 or @xmath252 implies that the corresponding approximated flow possesses some symmetry , such as a reflection @xmath253 , or a rotation , @xmath254 .
one of the most interesting of their results is that lorenz attractors can appear by local bifurcations of an equilibrium with three zero eigenvalues when there is a reflection symmetry through the origin @xcite,@xcite,@xcite .
these results can be applied to the local bifurcation analysis of ( [ eq:3hm1 ] ) when the fixed points has multipliers @xmath255 .    as is true for the general volume - preserving quadratic diffeomorphism of lomel and meiss @xcite , the 3d hnon map ( [ eq:3hm1 ] ) has at most two fixed points on the diagonal : @xmath256 providing that the discriminant , @xmath257 , is nonnegative , @xmath258    the primary bifurcations of these fixed points occur when they have multipliers on the unit circle .
these bifurcations are easily determined using the characteristic polynomial of the jacobian of ( [ eq:3hm1 ] ) , @xmath259 the first codimension - one case corresponds to a single multiplier one , where @xmath260 , or equivalently to the surface where @xmath257 vanishes : @xmath261 this is simply a saddle - node bifurcation that creates the pair of fixed points ( [ eq : fp ] ) in the domain @xmath262 .
similarly , a period - doubling bifurcation occurs when there is a multiplier @xmath251 , i.e. , where @xmath263 , or equivalently for @xmath264 , which gives @xmath265 note that this bifurcation occurs for the point @xmath266 when @xmath267 , and for the point @xmath268 when @xmath269 .
the final codimension - one bifurcation corresponds to the hopf bifurcation when pair of multipliers is on the unit circle , @xmath270 , so that @xmath271 or equivalently on the surface @xmath272 here we have written the equation for the surface @xmath273 in parametric form using the phase @xmath274 .
when @xmath275 , this bifurcation only occurs for the fixed point @xmath266 .
there are three codimension - two curves in the @xmath85 parameter space . the hopf surface
, @xmath276 terminates at @xmath277 when the multipliers are @xmath278 along the curve @xmath279 .
it also terminates on @xmath280 where the multipliers are @xmath281 along the curve @xmath282 .
finally , the period - doubling and saddle - node curves are tangent when the multipliers are @xmath283 on the curve @xmath284 .
the bifurcation surfaces are most easily viewed in two - dimensional slices through constant @xmath10 planes .
an example of the resulting bifurcation curves are shown when @xmath285 in fig .
[ fig : type21bifb06fp1 ] and fig .
[ fig : type21bifb06fp2 ] .
note that @xmath266 is stable inside the curvilinear triangle bounded by @xmath277 , @xmath280 and @xmath273 , when @xmath286 .
the second fixed point , @xmath268 , is never stable when @xmath275 .    ) for the fixed point @xmath266 when @xmath287 .
fixed points do not exist in the shaded region .
the solid curve is the saddle - node @xmath277 , the dotted curve is the period doubling , @xmath280 , and dot - dash curve is the hopf bifurcation , @xmath273 .
the fourth curve ( dashed ) corresponds to the existence of a double eigenvalue .
there are five corresponding regions , and a representative configuration of the multipliers in the complex-@xmath288-plane is shown for each .
the codimension - two bifurcations are shown as the  square " ( period - doubling - hopf ) ,  triangle " ( saddle - node - period doubling ) and   circle " ( saddle - node - hopf ) points .
[ fig : type21bifb06fp1],width=432 ]    ) in for the fixed point @xmath268 when @xmath287 .
the curves have the same meaning as in fig .
[ fig : type21bifb06fp1 ] . for this fixed point
there is no hopf bifurcation .
[ fig : type21bifb06fp2],width=384 ]    a similar analysis can be used to derive the bifurcation surfaces for the map ( [ eq:3hm2 ] ) .
indeed , the formula for the fixed points , and consequently the saddle - node surface is the same as for ( [ eq:3hm1 ] ) .
the period - doubling and hopf surfaces are slightly modified .
an example of a slice through these surfaces for @xmath285 is shown in fig .
[ fig : type12bifb06fp1 ] and fig .
[ fig : type12bifb06fp2 ] .    ) for the fixed point @xmath266 when @xmath287 .
the curves have the same meaning as in fig .
[ fig : type21bifb06fp1 ] .
[ fig : type12bifb06fp1],width=432 ]    ) in for the fixed point @xmath268 when @xmath287 .
the curves have the same meaning as in fig .
[ fig : type21bifb06fp1 ] .
for this fixed point there is no hopf bifurcation .
[ fig : type12bifb06fp2],width=384 ]
we will now focus on the point @xmath266 and expand around the codimension - three bifurcation point where the multipliers are @xmath289 , that is about the point @xmath290 where @xmath291 .
note that on the parameter slice @xmath54 , the hopf curve ( [ eq : hopf ] ) coincides with the saddle - node curve ( [ eq : sn ] ) so that there is no hopf bifurcation .
moreover the codimension - two points in fig .
[ fig : type21bifb06fp1 ] corresponding to the multipliers @xmath292 ( square ) and @xmath293 ( triangle ) collide at the codimension - three point .
otherwise the bifurcation diagram for @xmath266 looks similar the slice in the figure .    our goal is to prove th .
this can be done with the aid of the following lemma .
[ smlem ] in a region of the parameters @xmath294 near the point @xmath295 the second iterate of ( [ eq:3hm1 ] ) in some neighborhood of @xmath266 is approximately the time - one map of the the three - dimensional flow @xmath296 up to arbitrary small correction terms . moreover the coordinates @xmath297 and parameters @xmath288 and @xmath298 can take arbitrary finite values .
shifting coordinates @xmath299 brings map ( [ eq:3hm1 ] ) to the form @xmath300 introduce new ( small ) parameters @xmath301 by the following formulas @xmath302 note that @xmath303 the first inequality means that @xmath275 , and the second is fulfilled automatically when the fixed points exists , e.g. to the right of the @xmath277 curve in fig .
[ fig : type21bifb06fp1 ] .
we have also the following relation @xmath304 between the new parameters and the parameter @xmath8 .
next , we bring map ( [ hen1 ] ) to the following form @xmath305 where @xmath306 note that the jacobian of the map ( [ res1 ] ) is close to the jordan normal form for the multipliers @xmath289 when @xmath307 are small
. the first equation of ( [ res1 ] ) does not contain any high order terms , and in the second and third components , we take into account only quadratic terms .
moreover , these quadratic terms have been further simplified by standard normal form analysis , so that only the resonant ones remain .
the second power of the map ( [ res1 ] ) can be embedded into a flow .
namely , using the picard iteration procedure ( of the second order ) , we can see that the map ( [ res1 ] ) coincides , up to terms of order @xmath308 , with the time - one map of the following flow @xmath309    after a smooth change of coordinates of the form @xmath310 equations ( [ fl1 ] ) are written as follows @xmath311 where @xmath312 flow ( [ flow ] ) is the same that was considered @xcite and , in our case we have also that @xmath313 .
then , we can introduce new time @xmath314 and new rescaled coordinates @xmath297 , by formulas @xmath315 @xmath316 such that system ( [ flow ] ) is reduced to the following one @xmath317 where @xmath318 since values of @xmath319 can be arbitrarily small , it implies that the rescaled coordinates @xmath297 and new parameters @xmath288 and @xmath298 can take arbitrarily finite values .    note that if we omit the @xmath320 terms in the system ( [ sm ] ) , we obtain system ( [ = sm ] ) , which is the well - known , shimizu - morioka model .
this system was intensively studied in @xcite,@xcite where , in particular , the existence of the lorenz attractor was proved for some domain of the parameters .
for example , the domain of the existence of the lorenz attractor in shimizu - morioka model is known to be contained in a rectangle @xmath321 ( @xcite , fig . 13 ) .
using this fact we can pointed out the following boundaries of values of the parameters @xmath322 : @xmath323 ( and , respectively , of the parameters @xmath85 ) such that in the corresponding domain of the initial parameters @xmath85 our 3d hnon map has a wild - hyperbolic lorenz - type attractor .
this completes the proof of the main theorem .
suppose that @xmath324 is a quadratic cremona map , i.e. , is a diffeomorphism with constant jacobian @xmath325 .
then it is easy to see that without loss of generality we can write it in the form @xmath326 where @xmath327 is the diagonal scaling map with @xmath328 , @xmath200 is an affine volume - preserving map , @xmath329 is a vector of quadratic forms so that the map @xmath330 is a volume - preserving _
quadratic shear _
@xcite .    as in th .
[ tm:4 ] we now assume that the map @xmath38 has a polynomial inverse ( i.e. , it is an element of the affine cremona group ) and that the inverse is also a quadratic map .
the same property then will also hold for the volume - preserving map @xmath331 .
it was shown in @xcite that maps of the form @xmath332 can be transformed , by an affine change of variables , to one of the three following forms : @xmath333 where @xmath334 is a scalar quadratic form .
consequently @xmath335 for some affine transformation @xmath164 , and @xmath336 , @xmath337 , or @xmath338 .
using the same transformation on @xmath38 gives @xmath339 but the scaling map @xmath340 commutes with @xmath164 , so @xmath341 finally we can apply an additional scaling transformation to the variables and parameters to transform @xmath342 to one of the three forms @xmath343 where the parameters have been rescaled .
the dynamics of @xmath344 and @xmath345 are trivial .
indeed for @xmath344 the @xmath346 dynamics decouples from the remaining dynamics , and for @xmath345 the @xmath347 dynamics is linear and decouples from the @xmath348 dynamics .
thus the only dynamically interesting case is @xmath349 .
since this map is the same as ( [ eq:3hm3 ] ) , we have proved th .
[ tm:4 ] .
an example of this analysis consider the act - map ( arneodo - coullet - tresser map ) which was introduced by arneodo , coullet and tresser in their study of strange attractors in a family of differential equations on @xmath130 with homoclinic points of shilnikov type @xcite .
@xmath350 where @xmath351 are real parameters with @xmath352 , @xmath353 , and @xmath354 an integer .
the act - map always has a fixed point at the origin . nevertheless , when @xmath355 it is a quadratic family with a quadratic inverse , so it follows from th .
[ tm:4 ] that it is affinely conjugate to ( [ eq:3hm3 ] ) with @xmath356 .    to see this
explicitly define new coordinates using the transformation @xmath357 then , the map is rewritten as @xmath358 where where @xmath359 . for the case
@xmath355 we can rescale the coordinates by a factor @xmath360 to set the coefficient of @xmath361 to one . thus the map is written in the promised form ( [ eq:3hm3 ] ) .    note that if we do one additional transformation and shift the variables by @xmath362 we can also obtain the standard form ( [ eq:3hm2 ] ) . here
we define @xmath363 , and @xmath364 .
the inequality @xmath365 holds guaranteeing the existence of a fixed point for the map .
we have seen that many phenomena in the dynamics of three - dimensional hnon maps are markedly different from the possible phenomena predicted by our knowledge of two - dimensional dynamics .
primarily , we have seen these new phenomena when the jacobian , @xmath10 , is close to @xmath366 .
two natural problems appear here , regarding the understanding of bifurcation scenarios that lead to chaotic three - dimensional dynamics .
the first is to understand  how 2d dynamics is transformed into 3d dynamics " upon the variation of @xmath10 from values near zero to those close to @xmath367 .
the second problem  which is , in a sense , the inverse of the first one  is to understand
how 3d volume - preserving dynamics is transformed into dissipative ( and chaotic ) dynamics " as @xmath10 moves away from @xmath366 . in the current paper ( and also in @xcite ) , some partial results related to the second problem have been obtained ; however , the first problem is  as of yet  unstudied .
however , our analysis implies that these global problems have no unified , qualitative solution .
one general conclusion is that the stable dynamics of the hnon maps of the first and second kind , ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) , have unmistakably different characters .
the bifurcations in ( [ eq:3hm1 ] ) and ( [ eq:3hm2 ] ) both lead eventually to the formation of three - dimensional versions of smale horseshoes .
however , these have different characters : the horseshoes of ( [ eq:3hm1 ] ) contain type @xmath22 saddles of , while those of ( [ eq:3hm2 ] ) contain type @xmath23 saddles .
thus , a third global problem appears : to understand  how the character of the attractors changes as one changes type of nonlinearity in 3d hnon families . " in this case , it seems that using the family ( [ eq:3hm3 ] ) with the parameters @xmath368 determining the type of nonlinearity may be very profitable
.    while the problems posed above are of general interest , there are also many partial problems relating directly to the dynamics of quadratic maps .
some of these may also be solved in the the near future :    1 .   does the hnon map of the second kind , ( [ eq:3hm2 ] ) , have a large scale , wild - hyperbolic attractor ? 2 .
what are the dynamical properties of the local normal form for ( [ eq:3hm2 ] ) near fixed points with three unit multipliers ?
+ this is related to the `` reverse - shimitsu - marioka '' system since ( [ eq:3hm2 ] ) is the inverse of ( [ eq:3hm1 ] ) .
thus we will obtain the model ( [ = sm ] ) , with @xmath369 replaced by @xmath370 .
consequently , the question is to determine the properties of the stable dynamics of this system .
what are main bifurcations scenarios the 3d hnon map ( [ eq:3hm1 ] ) that lead from simple dynamics to a wild - hyperbolic , lorenz attractor ?
+ there are two distinct , but useful approaches : ( a ) to fix @xmath10 and consider the corresponding two - parameter family ; ( b ) to consider @xmath10 as the main governing parameter and to move for given values of @xmath371 from small values of @xmath10 to values close to @xmath366 .
what role is played by the codimension - two bifurcations ( recall  [ sec : bifurcations ] ) in the formation of attracting invariant sets ? + the analysis of codimension - one bifurcations alone gives us only a representation of the main hyperbolic elements of the dynamics .
a detail analysis of codimension - three bifurcations ( as in
[ sec : bifurcations ] ) allows us to find in local normal forms ( using the flow embedding approach ) and verify the existence of small , wild - hyperbolic attractors .
analogous methods can be applied to the study of codimension - two bifurcations .
we believe that this approach will provide an understanding of principal mechanisms for the creation of wild - hyperbolic attractors .    in summary , three - dimensional hnon maps have an abundance of new dynamical behavior which , it seems , has an important place in the theory of multidimensional dynamics as whole .
we believe that they form a natural direction for a program of investigation similar to the well developed programs in one and two - dimensional dynamics .
this is especially true since these maps cover ( up to the affine conjugacy ) all possible quadratic three - dimensional automorphisms .
dragt and d.t .
symplectic maps and computation of orbits in particle accelerators . in _ integration algorithms and classical mechanics ( toronto , on , 1993 ) _ , pages 5985 .
soc . , providence , ri , 1996 .
dynamics and moduli of @xmath375-conjugacy of 4d - diffeomorphisms with a structurally unstable homoclinic orbit to a saddle - focus fixed point . , 200:107134 , 2000 .
providence , ri .
bifurcation and chaos in the marioka - shimizu system . in _ methods of qualitative theory of differential equations _ , pages 180193 .
gorky , 1986 .
english translation in selecta math .
soviet . 10
( 1991 ) , 105 - 117 .