after the effect of spin transfer torque ( stt ) in thin magnetic films had been predicted  @xcite and then experimentally proven  @xcite , it was generally believed that current - controlled spin valve devices would replace soon the memory cells operated by external magnetic field .
such a technological progress , if realized , would certainly offer higher data storage density and faster manipulation with the information stored on a hard drive memory . however , it became clear soon that some important issues must be solved before devices based on spin torque could be used in practice .
the most important is the reduction of current density needed for magnetic excitation ( switching ) in thin films , as well as enhancement of switching efficiency and thermal stability .
some progress has been made by using more complex spin valve structures and/or various subtle switching schemes based on optimized current and field pulses  @xcite .
a significant enhancement of thermal stability can be achieved by replacing a simple free layer ( single homogeneous layer ) with a system of two magnetic films separated by a thin nonmagnetic spacer , known as composite free layer ( cfl ) .
the spacer layer is usually thin so there is a strong rkky exchange coupling between magnetic layers  @xcite . in practice
, antiferromagnetic configuration is preferred as it reduces the overall magnetic moment of the cfl structure and makes the system less vulnerable to external magnetic fields and thermal agitation . when the antiferromagnetically coupled magnetic layers are identical , we call the structure synthetic antiferromagnet ( syaf ) . if they are different , then the cfl has uncompensated magnetic moment and such a system will be referred to as synthetic ferrimagnet ( syf ) .
current and/or field induced dynamics of cfls is currently a subject of both experimental and theoretical investigations  @xcite .
switching scheme of syaf by magnetic field pulses has been proposed in ref .
[ ] , and then the possibility of current - induced switching of syaf was demonstrated experimentally  @xcite . in turn , the possibility of critical current reduction has been shown for a cfl with ferromagnetically coupled magnetic layers  @xcite .
however , the reduction of critical current in the case of antiferromagnetically coupled cfls still remains an open problem . in a recent numerical study on switching a syaf free layer  @xcite it has been shown that the corresponding critical current in most cases is higher than the current required for switching of a simple free layer , and only in a narrow range of relevant parameters ( exchange coupling , layer thickness , etc . ) the critical current is reduced .
hence , proper understanding of current - induced dynamics of cfls is highly required .
we also note , that cfl can be used as a polarizer , too .
indeed , it has been shown recently  @xcite that syaf used as a reference layer ( with magnetic moment fixed fixed to adjacent antiferromangentic layer due to exchange anisotropy ) might be excited due to dynamical coupling  @xcite with a simple sensing layer  @xcite .
the main objective of this paper is to study current - induced dynamics of a cfl with antiferomagnetic rkky coupling in metallic spin valve pillars .
we consider a system @xmath0/@xmath1/@xmath2/@xmath3/@xmath4/@xmath5 , shown in fig .
[ fig : scheme ] , where @xmath0 is an antiferromagnetic layer ( used to bias magnetization of the reference magnetic layer @xmath1 ) , @xmath3 and @xmath5 are two magnetic layers , while @xmath2 and @xmath4 are non - magnetic spacers .
the part @xmath3/@xmath4/@xmath5 constitutes the cfl structure with antiferromagnetic interlayer exchange coupling .
we examine current - induced dynamics of both syaf and syf free layers .
these two structures differ only in the thickness of @xmath3 layer , while rkky coupling and other pillar parameters remain the same .
we assume that spin - dependent electron transport is diffusive in nature , and employ the model described in refs .
[ ] and  [ ] . an important advantage of this model is the fact that it enables calculating spin current components and spin accumulation consistently in all magnetic and nonmagnetic layers , as well as current - induced torques exerted on all magnetic components .
the torques acting at the internal interfaces of cfl introduce additional dynamical coupling between the corresponding magnetic layers .
consequently , the magnetic dynamics of cfl has been modelled by two coupled macrospins and described in terms of landau - lifshitz - gilbert ( llg ) equation .
in addition , we derive some analytical expressions for critical currents from the stability conditions of linearized llg equation in the static points  @xcite , and discuss results in the context of numerical simulations .    the paper is organized as follows . in section 2
we describe the assumed models for spin dynamics and stt calculations . in section
3 we analyze stt acting on cfl and present results from numerical simulations on current - induced switching and magnetic dynamics . some additional information on stt calculation can be found in the appendix .
critical currents are derived and discussed in section  4 .
finally , summary and general conclusions are in section  5 .
scheme of the spin valve pillar structure with a composite free layer . ]    in the macrospin approximation , magnetization dynamics of the cfl is described by two coupled landau - lifshitz - gilbert ( llg ) equations , @xmath6 for @xmath7 , where @xmath8 stands for a unit vector along the net spin moment of the @xmath9-th layer , whereas @xmath10 and @xmath11 are the effective field and current - induced torque , respectively , both acting on @xmath8 .
the damping parameter @xmath12 and the saturation magnetization @xmath13 are assumed the same for both magnetic components of the cfl .
furthermore , @xmath14 is the gyromagnetic ratio , @xmath15 is the vacuum permeability , and @xmath16 stands for thickness of the @xmath17 layer .    the effective magnetic field for the @xmath17 layer is @xmath18 where @xmath19 and @xmath20 . in the latter equation , @xmath21 is the external magnetic field applied along the easy axis in the layers plane ( and oriented opposite to the axis @xmath22 ) , @xmath23 is the uniaxial anisotropy field ( the same for both magnetic layers ) , and @xmath24 is the self - demagnetization field of the @xmath17 layer with @xmath25 being the corresponding demagnetization tensor .
similarly , @xmath26 describes the magnetostatic influence of the layers @xmath1 and @xmath27 on the layer @xmath17 , respectively . here , @xmath28 is the saturation magnetization of the layer @xmath29 , which might be generally different from @xmath13 .
components of the tensors @xmath25 , @xmath30 and @xmath31 used in our simulations have been determined by the numerical method introduced by newell _
et al . _
this method was originally developed for magnetic systems with non - uniform magnetization . to implement it into a macrospin model we considered discretized magnetic layers with uniform magnetizations , calculated tensors in each cell of the layer , and then averaged them along the whole layer .
since these tensors are diagonal , the demagnetization and magnetostatic fields can be expressed as @xmath32 , and @xmath33 , with @xmath34 , @xmath35 , and @xmath36 denoting the components of the vector @xmath8 ( @xmath37 ) in the coordinate system shown in fig.1 .
finally , @xmath38 stands for the rkky exchange field acting on @xmath8 , which is related to the rkky coupling constant as @xmath39  @xcite .    to include thermal effects
we add to the effective field  ( [ eq : heff ] ) a stochastic thermal field @xmath40 .
for both spins its components obey the rules for gaussian random processes @xmath41 and @xmath42 , where @xmath43 and @xmath44 . here
, @xmath45 is the noise amplitude , which is related to the effective temperature , @xmath46 , as @xmath47 where @xmath48 is the boltzmann constant , and @xmath49 is the volume of @xmath17 layer .
in general , the current - induced torques acting on @xmath50 and @xmath51 can be expressed as a sum of their in - plane and out - of - plane components @xmath52 and @xmath53 , respectively . in a cfl structure , the layer @xmath3 is influenced by stt induced by the polarizer @xmath1 , as well as by stt due to the layer @xmath5 . in turn , the layer @xmath5 is influenced by the torques from the layer @xmath3 .
hence we can write    [ eqs : stt ] @xmath54\ , , \\        { { \bm \tau}_{1 \perp } } & = i\bhs_1 \times \left ( b_1^{(0 ) } \bhs_0 + b_1^{(2 ) } \bhs_2 \right)\ , , \\        { { \bm \tau}_{2 \parallel } } & = ia_2^{(1)}\ , \bhs_2 \times \left ( \bhs_2 \times \bhs_1 \right)\ , , \\        { { \bm \tau}_{2 \perp } } & = ib_2^{(1)}\ , \bhs_2 \times \bhs_1\ , ,      \end{aligned}\ ] ]    where @xmath55 is the charge current density , which is positive when electrons flow from the layer @xmath5 towards @xmath1 ( see fig .
[ fig : scheme ] ) , while the parameters @xmath56 and @xmath57 ( @xmath58 ) are independent of current @xmath55 , but generally depend on magnetic configuration .
we write the current density in the spin space as @xmath59 , where @xmath60 is the particle current density ( @xmath61 ) , @xmath62 is the spin current density ( in the units of @xmath63 ) , @xmath64 is the vector of pauli matrices , and @xmath65 is a @xmath66 unit matrix . in frame of the diffusive transport model  @xcite , the parameters @xmath56 and @xmath57 are given by the formulas @xmath67 where the angles @xmath68 and @xmath69 are given by @xmath70 and @xmath71 . here , @xmath72 and @xmath73 are transversal ( to @xmath50 ) components of spin current in the layer @xmath2 ( taken at the @xmath2/@xmath3 interface ) written in the local coordinate system of @xmath50 .
thus , @xmath72 is parallel to the vector @xmath74 ( lying in the plane defined by @xmath75 and @xmath50 ) and @xmath73 is aligned with @xmath76 ( perpendicular to the plane defined by @xmath75 and @xmath50 ) .
note , that the @xmath22-component of spin current is aligned along @xmath50 and does not contribute to stt  @xcite .
similarly , we define the torque amplitudes acting inside the cfl . here ,
@xmath77 and @xmath78 are the components of spin current in the layer @xmath4 ( taken at the n2/f1 interface ) , transversal to @xmath50 and lying respectively in the plane and perpendicularly to the plane defined by @xmath50 and @xmath51 .
analogically , @xmath79 and @xmath80 are spin current components in @xmath4 transversal to @xmath51 .
the relevant spin current transformations are given in appendix  [ sec : transformations ] .
note , that the spin current components depend linearly on @xmath55 , so the parameters @xmath56 and @xmath57 are independent of the current density @xmath55 .
in this section we present results on our numerical simulations of current - induced dynamics for two metallic pillar structures including cfl with antiferromagnetic interlayer exchange coupling .
as described in the introduction , the considered pillars have the general structure @xmath0/@xmath1/@xmath2/@xmath3/@xmath4/@xmath5 ( see fig.(1 ) ) . more specifically , we consider spin valves cu - irmn(10)/py(8)/cu(8)/co(@xmath81)/ru(1)/co(@xmath82 ) - cu , where the numbers in brackets stand for the layer thicknesses in nanometers .
the layer py(8 ) is the permalloy polarizing layer with its magnetization fixed due to exchange coupling to irmn .
in turn , co(@xmath81)/ru(1)/co(@xmath82 ) is the cfl ( @xmath3/@xmath4/@xmath5 structure ) with antiferromagnetic rkky exchange coupling _ via _ the thin ruthenium layer .
the coupling constant between co layers has been assumed as @xmath83 , which is close to experimentally observed values  @xcite .
here , we shall analyze two different geometries of cfl .
the first one is a syaf structure with @xmath84 , while the second one is a synthetic ferrimagnet ( syf ) with @xmath85 .
simulations have been based on numerical integration of the two coupled llg equations  ( [ eq : llg ] ) with simultaneous calculations of stt , see eq.([eqs : stt ] ) .
we have assumed typical values of the relevant parameters , i.e. , the damping parameter has been set to @xmath86 , while the uniaxial anisotropy field @xmath87 in both magnetic layers of the cfl . in turn
, saturation magnetization of cobalt has been assumed as @xmath88 , and for permalloy @xmath89 .
the demagnetization field and magnetostatic interaction of magnetic layers have been calculated for layers of elliptical cross - section , with the major and minor axes equal to @xmath90 and @xmath91 , respectively .    for both structures under consideration
we have analyzed the current - induced dynamics as a function of current density and external magnetic field .
the results have been presented in the form of diagrams displaying time - averaged values of the pillar resistance .
numerical integration of  eq.([eq : llg ] ) has been performed using corrector  predictor heun scheme , and the results have been verified for integration steps in the range from @xmath92 up to @xmath93 .
the stt components acting on cfl spins have been calculated at each integration step from the spin currents , which have been numerically calculated from the appropriate boundary conditions  @xcite .
similarly , resistance of the studied pillars has been calculated from spin accumulation in frame of the model used also for the stt description ( for details see also ref .  ) .
let us analyze first the angular dependence of stt components in the structures under consideration .
although the thicknesses of magnetic layers in the studied syaf and syf structures are different , the angular dependence of stt components as well as their amplitudes are very similar . thus , the analysis of stt in syaf applies also qualitatively to the studied syf free layer .
( color online ) angular dependence of the cartesian components of stt , in the units of @xmath94 , acting at n / f interfaces , when magnetization of the syaf structure is rotated rigidly with both @xmath50 and @xmath51 remaining in the corresponding layer planes . here
, @xmath95 is an angle between @xmath50 and @xmath96 .
@xmath51 is tilted away from the antiparallel configuration with @xmath50 by an angle of @xmath97.,width=302 ]    first , we analyze stt components in the case when syaf is rotated as a rigid structure , i.e. the antiparallel configuration of @xmath50 and @xmath51 is maintained . to have a nonzero torque between @xmath3 and @xmath5 layers , @xmath51 has been tilted away from the antiparallel configuration by an angle of @xmath97 .
figure  [ fig : stt1 ] shows all three cartesian components ( see fig.1 for definition of the coordinate system ) of stt acting at n / f interfaces as a function of the angle @xmath95 between @xmath96 and @xmath50 . while the @xmath98 and @xmath22-components are in the plane of the layers ( the spins of cfl are rotated in the layer plane ) , the component @xmath99 is normal to the layer plane .
however , @xmath100 remains negligible at all interfaces of the cfl .
the stt acting at @xmath2/@xmath3 reveals a standard ( non - wavy  @xcite ) angular dependence , and vanishes when @xmath50 is collinear with @xmath75 .
its amplitude is comparable to stt in standard spin valves with a simple free layer .
the stt at @xmath3/@xmath4 and @xmath4/@xmath5 interfaces also depends on the angle @xmath95 . however , they are about two order of magnitude smaller , which is a consequence of a small angle ( @xmath97 ) assumed between @xmath50 and @xmath51 .    as will be shown in the following ,
cfl is usually not switched as a rigid structure , but generally forms a configuration which deviates from the antiparallel one .
figure  [ fig : stt2 ] shows how the stt components at the @xmath3/@xmath4 and @xmath4/@xmath5 interfaces vary when @xmath51 is rotated from @xmath96 by an angle @xmath101 , while @xmath50 remains fixed and is parallel to @xmath102 .
in such a case , the torque acting at @xmath2/@xmath3 interface remains zero , as @xmath50 stays collinear to @xmath75 .
as before , the out - of - plane components are also negligible in comparison to the in - plane ones .
the in - plane components of stt reveal standard angular dependence at both interfaces .
the amplitude of stt at the internal interfaces of cfl is comparable to that acting at the @xmath2/@xmath3 interface in the case of noncollinear configuration of @xmath75 and @xmath50 , when @xmath51 is fixed in the direction antiparallel to @xmath75 .
( color online ) angular dependence of the cartesian stt components in the units of @xmath103 acting on f1/n2 and n2/f1 interfaces when @xmath104 and @xmath51 is rotated from @xmath96 in the layer s plane by an angle @xmath101.,width=302 ]    when @xmath50 is noncollinear to @xmath75 , the spin accumulation in @xmath2 layer increases and consequently the amplitude of stt at @xmath3/@xmath4 and @xmath4/@xmath5 decreases . in turn , when @xmath50 is antiparallel to @xmath75 , the stt inside the cfl structure is reduced by more than a factor of 2 .
nonetheless , the stt acting at the internal interfaces of the studied cfl layers might have a significant effect on their current - induced dynamics and switching process , provided the magnetic configuration of cfl might deviate remarkably from its initial antiparallel configuration .
first , we examine dynamics of the syaf free layer . from symmetry
we have @xmath105 , and we have set @xmath106 , which corresponds to @xmath107 .
we have performed a number of independent numerical simulations modelling syaf dynamics induced by constant current and constant in - plane external magnetic field .
the latter is assumed to be smaller than the critical field for transition to spin - flop phase of syaf .
accordingly , each simulation started from an initial state close to @xmath108 . to have a non - zero initial stt for @xmath50 ,
both spins of the syaf have been tilted by @xmath97 in the layer plane so that they remained collinear .    from the results of numerical simulations
we have constructed a map of time - averaged resistance , shown in fig .  [
fig : saf_dynam1](a ) .
the resistance has been averaged in the time interval of @xmath109 following initial @xmath110 equilibration time of the dynamics .
the diagram shows only that part of the resistance , which depends on magnetic configuration , and hence varies with cfl dynamics  @xcite .
the constant part of resistance , due to bulk and interfacial resistances of the studied structure , has been calculated to be as large as @xmath111 .
( color online ) ( a ) averaged resistance of cu - irmn(10)/py(8)/cu(8)/co(2)/cu(1)/co(2)-cu spin valve pillar with a saf free layer as a function of current density and applied magnetic field .
examples of switching processes at @xmath112 and @xmath113 ( b  e ) and @xmath114 ( f  i ) . figures ( b ) and ( f ) show dynamics of @xmath22-components of both spin moments , ( c ) and ( g ) present the overal magnetization of the free layer , ( d ) and ( h ) show the corresponding variation of pillar resistance , ( e ) and ( i ) show the spin trajectories of @xmath50 ( red solid line ) and @xmath51 ( black dashed line ) in the time interval from @xmath115 to @xmath116 , where switching takes place .
, width=302 ]    for the assumed initial configuration , magnetic dynamics has been observed only for negative current density . when the current is small , no dynamics is observed since the spin motion is damped into the closest collinear state ( @xmath108 , marked as @xmath117 ) of high resistance .
after exceeding a certain threshold value of current density , there is a drop in the averaged resistance , which indicates current induced dynamics of the syaf free layer .
figures  [ fig : saf_dynam1](b ) and ( f ) show that this drop is associated with switching of the whole syaf structure into an oposite state ( @xmath118 , marked as @xmath119 )    from fig.4(a ) follows that the threshold current for dynamics onset markedly depends on the applied field and reaches maximum at a certain value of @xmath21 , @xmath120 .
furthermore , it appears that mechanisms of the switching process for @xmath121 and @xmath122 are qualitatively different . to distinguish these two mechanisms , we present in figs .
[ fig : saf_dynam1](b  i ) basic characteristics of switching , calculated for @xmath123 and for fields @xmath113 , which is below @xmath124 [ figs .  [ fig : saf_dynam1](b  e ) ] , and @xmath125 , which lies above @xmath124 [ figs .  [
fig : saf_dynam1](f  i ) ] . figs .
[ fig : saf_dynam1](b ) and ( f ) present time evolution of the z - components of both spins . to better understand the syaf dynamics , in figs .
[ fig : saf_dynam1](c ) and ( g ) we plotted the amplitude of overal syaf magnetization , defined as @xmath126 .
this parameter vanishes for antiparallel alignment of both spins of cfl , but becomes nonzero when the configuration deviates from the antiparallel one .
magnetization of syaf is also a measure of the cfl coupling to external magnetic field .
futhermore , figs .
[ fig : saf_dynam1](d ) and ( h ) show the corresponding time variation of the resistance , @xmath127 , which might be directly extracted from experimental measurements as well .
in addition , in figs .  [
fig : saf_dynam1](e ) and ( i ) we show the trajectories of @xmath50 and @xmath51 in the real space taken from the time interval from @xmath128 to @xmath129 . in addition , from figs .
[ fig : saf_dynam1](a ) it has been found that the point where the threshold current reaches its maximum is located at @xmath130 , which indicates its relation to magnetostatic interaction of f2 and fixed polarizer .
this also has been confirmed by analogical simulations disregarding the magnetostatic coupling between magnetic layers , which resulted in similar diagram , but with @xmath131 ( not shown ) .
this fact significantly facilitates understanding the mechanism of syaf switching .
the initial configuration assumed above was @xmath132 with @xmath102 ( @xmath133 ) .
when the magnitude of current density is large enough and @xmath134 , orientation of @xmath50 becomes unstable and @xmath50 starts to precess with small angle around @xmath135 .
initial precession of @xmath50 induces precession of @xmath51  mainly via the rkky coupling .
generally , response to the exchange field is slower than current - induced dynamics .
therefore , a difference in precession phase of @xmath51 and @xmath50 appears , and configuration of syaf deviates from the initial antiparallel one .
this in turn enhances the stt acting on f2 , which tends to switch @xmath51 .
its amplitude , however , is small in comparison to the strong rkky coupling .
further scenario of the dynamics depends then on the external magnetic field . when @xmath121 [ figs .
[ fig : saf_dynam1](b  e ) ] the zeeman energy of @xmath51 has a maximum in the initial state and external magnetic field tends to switch @xmath51 to the opposite orientation .
competition between the torques acting on syaf results in out - of - plane precessions of both spins .
after several precessions @xmath50 reaches the opposite static state , which is stable due to stt . in turn
, @xmath51 is only slightly affected by stt , and its dynamics is damped in the external magnetic and rkky exchange fields .
in contrast , when @xmath122 [ figs .
[ fig : saf_dynam1](f  i ) ] , zeeman energy of f2 has a local minimum in the initial state , which stabilizes @xmath51 .
therefore , in a certain range of current density , syaf does not switch but remains in self - sustained coherent in - plane precessions ( red area in the upper part of fig .
[ fig : saf_dynam1](a ) ) . for a sufficient current density ,
the syaf structure becomes destabilized and the precessional angle increases until the spins pass the @xmath136-plane .
consequently , the precessional angle decreases and spins of the syaf are stabilized in the opposite state ( @xmath137 ) . moreover , as shown in fig .
[ fig : saf_dynam1](c ) , the switching process for @xmath121 is connected with high distortion of syaf configuration , where @xmath138 in a certain point reaches its maximum value ( corresponding to parallel orientation of both spins ) .
contrary , the @xmath138 remains small for @xmath122 [ figs .  [
fig : saf_dynam1](g ) ] , and the effective magnetic moment of the free layer stays smaller than magnetic moment of a single layer . this might play an important role in applications of spin - torque devices based on cfls .
the two switching mechanisms described above dominate the current - induced dynamics when the current density is close to the dynamics threshold . for higher current densities ,
the nonlinearities in syaf dynamics become more pronounced , which results in bistable behavior of the dynamics , especially for @xmath121 and @xmath139 . in that region , the number of out - of - plane precessions before syaf switching increases with the current density .
however , their precessional angle increases in time and consequently @xmath50 might reach an out - of - plane static point slightly tilted away from the @xmath140 direction while @xmath141 remains in the layers plane .
the out - of - plane static states ( marked as @xmath142 ) have small resistance and appear as dark red spots in the diagram shown in fig .
[ fig : saf_dynam1](a ) .
in addition , from the analysis of the dispersion of pillar resistance ( not shown ) one finds that except of a narrow region close to the dynamics threshold with persistent in - plane precessions , no significant steady - state dynamics of syaf element appears . as will be shown below , such a behavior
might be observed when cfl becomes asymmetric ( syf free layer ) .
let us study now spin valve with syf as a free layer , assuming @xmath143 and @xmath144 .
accordingly , @xmath145 remains @xmath146 while @xmath147 is reduced to @xmath148 .
as in the case of syaf , from the averaged time - dependent part of the pillar resistance we have constructed a diagram presenting current - induced dynamics , see fig .
[ fig : sf_dynam1](a ) .
the static part of resistance is now @xmath149 .
( color online ) ( a ) averaged resistance of cu - irmn(10)/py(8)/cu(8)/co(4)/cu(1)/co(2)-cu spin valve pillar with a syf free layer , presented as a function of current density and applied magnetic field .
examples of current - induced dynamics for @xmath150 and @xmath151 ( b  e ) and @xmath113 ( f
i ) . panels ( b ) and ( f ) show dynamics of @xmath22-components of @xmath50 and @xmath51 , ( c ) and ( g ) present the overal magnetization of the free layer , ( d ) and ( h ) show the corresponding variation of pillar resistance , ( e ) and ( i ) show spin trajectories of @xmath50 ( red solid line ) and @xmath51 ( black dashed line ) taken from a time interval as large as @xmath152 after @xmath153 of initial equilibration.,width=302 ]    the diagram has some features similar to those studied in the previous subsection .
however , the maximum critical current is shifted towards negative values of @xmath21 , even if magnetostatic interaction between magnetic layers is neglected .
this asymmetry is caused by the difference in exchange and demagnetization fields acting on layers @xmath3 and @xmath5 .
moreover , this difference leads to more complex dynamics of the cfl s spins than that in the case of syaf .
generally , there are several dynamic regimes to be distinguished .
the first one is the region of switching from @xmath117 configuration to the opposite one , @xmath119 , which is located at largest values of @xmath21 in the diagram .
mechanism of the switching is similar to that of syaf shown in figs .
[ fig : saf_dynam1](f  i ) , where cfl changes its configuration just _ via _ in - plane precessional states with a small value of @xmath138 ( weak distortion of the antiparallel alignment of @xmath50 and @xmath51 )
. furthermore , the darker area above @xmath154 indicates one of the possible self - sustained dynamic regimes of syf , i.e. the in - plane precessions ( ipp ) ; see figs .
[ fig : sf_dynam1](b  e ) .
this precessional regime starts directly after the syf switching , and @xmath50 and @xmath51 precess around @xmath96 and @xmath135 , respectively . due to different effective fields in @xmath3 and @xmath5 , and energy gains due to stt , the spins precess with different precessional angles [ fig .
[ fig : sf_dynam1](e ) ] and consequently different frequencies .
because of the strong interlayer coupling and spin transfer between the layers , amplitudes of their precessions are periodically modulated in time .
this modulation appears also in the time dependence of pillar resistance .
conversely , below @xmath154 the dynamic is dominated by large angle out - of - plane precessions ( opp ) of both spins , as shown in figs .
[ fig : sf_dynam1](f  i ) .
this dynamic state is connected with a strong distortion of the antiparallel cfl configuration , i.e. large value of @xmath138 , and large variation of the resistance . from fig .
[ fig : sf_dynam1](i ) one can see that trajectories of @xmath50 and @xmath51 are rather complicated including both ipp and opp regimes with dominant opp .      from the analysis of current - induced dynamics we found that self - sustained dynamics in structures with syf free layer
is much richer than that in systems with syaf free layer [ see figs .  [
fig : sf_dynam1](b  i ) ] .
therefore , in this section we restrict ourselves to dynamic regimes of the syf free layer only .
more specifically , we shall examine the power spectral density ( psd ) as a function of current density and external magnetic field .    in the simulation we started from @xmath155 and changed current density in steps @xmath156 at a fixed applied field .
as before , to protect the syf dynamics from collapsing into collinear static state , we assumed small thermal fluctuations corresponding to @xmath157 . at each step
we simulated the dynamics of coupled cfl s spins and calculated psd .
as in ref .  , we assumed that the input current is split between a load with resistance @xmath158 and pillar with resistance @xmath159 .
hence , voltage on the pillar has been calculated as @xmath160 $ ] , where we assumed @xmath161 , and @xmath162 is the cross - section of the pillar ( ellipsoid with the major and minor axes equal to @xmath163 and @xmath164  nm , respectively ) .
then , at a given @xmath55 we calculated voltage in the frequency domain , @xmath165 , using fast fourier transformation over the period @xmath166 following the equilibration time of @xmath167
. the power spectral density has been defined as @xmath168 , where @xmath169 .
power spectral density calculated for the spin valve with syf free layer at @xmath157 and @xmath113 ( a ) and @xmath170 ( b ) .
( c ) and ( d ) show steady time evolution of spins @xmath22-components in a time window of @xmath171 after the equilibration at @xmath113 and @xmath172 .
( e ) and ( f ) show steady time evolution of the time - dependent part of spin valve resistance in a time window of @xmath173 after the equilibration at @xmath113 and @xmath174 ( both @xmath50 and @xmath51 precess out of layer plane ) and @xmath175 ( @xmath51 performs out - of - plane precessions while @xmath50 precesses in the layer s plane ) , respectively .
panels ( g ) and ( h ) depict trajectories of @xmath50 ( red solid line ) and @xmath51 ( black dashed line ) corresponding to resistance oscillations ( e ) and ( f ) , respectivelly . ]
figures  [ fig : psd](a ) and ( b ) show psd calculated at @xmath113 and @xmath170 , respectively .
the former case corresponds to that part of the diagram in fig .
[ fig : sf_dynam1](a ) , which includes opp modes , while in the latter case we observed ipp only .
let us analyze first the situation in fig .
[ fig : psd](a ) .
when current passes through the corresponding threshold value , both spins start precessing in the layers plane , similarly as shown in fig .
[ fig : sf_dynam1](b ) .
apart from the main peak in the psd at @xmath176 , two additional minor peaks close to @xmath177 are visible .
we attribute them to the oscillations of precessional amplitudes of both spins . with increasing amplitude of the current density ,
the precessional angles of both spins increase and their precessional frequencies slightly decrease .
moreover , with increasing current the frequencies of the minor peaks become closer and closer , until they finally coincide . at this point
the psd becomes widely distributed along the whole range of observed frequencies , which is an evidence of noisy variation of the resistance .
an example of spin dynamics in this region is shown in figs .
[ fig : psd](c ) and ( d ) which have been taken in a time window as large as @xmath171 after the equilibration period for @xmath113 and @xmath178 [ within the broad feature of psd in fig .
[ fig : psd](a ) ] .
firstly , the figures show that @xmath51 starts to perform out - of - plane precessions as a result of the competition between stt and rkky coupling .
secondly , one can note thermally activated random transitions of @xmath50 between opp and ipp modes .
these random transitions modify opp precessions of @xmath51 as well .
simultaneous dynamics of both spins causes chaotic variation of spin valve resistance and broadens the psd .
the quasi - chaotic feature of the spin dynamics in this range of current densities can be seen also on the spin trajectories , which cover almost the whole sphere ( not shown ) .
further increase in current density leads to stabilization of the opp mode of @xmath50 .
hence spin valve resistance becomes more periodic [ see fig .  [
fig : psd](e ) ] and psd reveals a narrow peak again . since both spins perform rather complicated dynamics including ipp but dominated by opp regime [ see fig .  [
fig : psd](g ) ] , we observe a blue - shift in psd with current , which is connected with a decrease in the precessional angles .
however , at a certain value of @xmath55 we notice an abrupt drop in the peak s frequency . at this current density
the stt acting on the left interface of layer @xmath3 starts to dominate the dynamics of @xmath50 and enables only small angle ipps along the @xmath75 direction , which modifies the trajectory of @xmath51 .
@xmath51 still remains in the opp regime [ see fig .  [
fig : psd](h ) ] and hence the blue - shift with current appears .
the fact that ipp of @xmath50 still influence the dynamics of the whole syf is also shown in fig .  [
fig : psd](f ) , which presents the dynamic part of the spin valve resistance at @xmath179 and @xmath113 . as a result of ipps of @xmath50 , amplitude of the resistance varies periodically .
in addition , comparison of figs .
[ fig : psd](e ) and ( f ) shows that the simultaneous opps of both spins lead to stronger variation of the resistance than in the case when the layers are in the ipp state .
contrary , at @xmath180 one observes only ipp modes of both spins similar to those shown in fig .
[ fig : sf_dynam1](e ) .
the in - plane precessional angle increases with current density and hence the peak frequency in psd decreases and becomes broader . in real systems , however , one might expect the peaks narrower than those obtained in the macrospin simulations , as observed in standard spin valves with a simple free layer  @xcite .
first , we derive some approximate expressions for critical current density needed to induce dynamics of cfl , derived from linearized llg equation . in metallic structures , the out - of - plane torque components are generally much smaller than the in - plane ones , and therefore will be omitted in the analytical considerations of this section ( @xmath181 ) .
the coupled llg equations in spherical coordinates can be then written as @xmath182 where @xmath183 is a 4-dimensional column vector which describes spin orientation in both layers constituting the cfl , and @xmath184 stands for the torque components , @xmath185 and @xmath186 , with @xmath187 and @xmath188 denoting unit vectors in local spherical coordinates associated with @xmath8 . in turn , the @xmath189 matrix @xmath190 takes the form @xmath191 static points of the cfl dynamics have to satisfy @xmath192 and @xmath193 for both @xmath194 and @xmath195 .
these conditions are obeyed in all collinear configurations , i.e. @xmath196 .
additional four trivial static points can be found in the out - of - plane configurations with @xmath197 and @xmath198 .
following ref .
we linearize eq .
( [ eq : llg_spherical ] ) by expanding @xmath199 into a series around the static points , which leads to @xmath200 where @xmath201 is a jacobian matrix of @xmath202 components .
the matrix product @xmath203 defines here the dynamic matrix @xmath204 .
this matrix allows one to study stability of cfl s spins in their static points .
if @xmath205 is negative , the static point is stable , otherwise it is unstable .
hence , the condition for critical current is  @xcite @xmath206 .    to obtain threshold current for dynamics onset of individual spins in the cfl , we assume first that one of the spins is fixed in its initial position and investigate stability of the second spin . the dynamic matrix @xmath207 becomes then reduced to a @xmath208 matrix .
considering initial position of saf with @xmath108 ( i.e. @xmath209 and @xmath210 ) , marked as @xmath117 , and polarizer @xmath211 , the stability condition leads to the following critical currents @xmath212 and @xmath213 for @xmath50 and @xmath51 , respectively : @xmath214\ , ,    \label{eq : critcurr1_du}\ ] ] with @xmath215 , and @xmath216\ , ,    \label{eq : critcurr2_du}\ ] ] with @xmath217 . in both above expressions @xmath218 , @xmath219 , and @xmath220
are taken in the considered static point , while the demagnetization field for the @xmath9-th layer is given by @xmath221 analogically , one can derive similar formulas for critical currents in the opposite @xmath222 magnetic configuration of the cfl .
now we relax the assumption that one of the spins is fixed , and consider both spins of the cfl as free .
then , we calculate the trace of the whole @xmath189 matrix , which leads to the following expression for critical current destabilizing the whole cfl : @xmath223\ , ,      \end{split}\ ] ] where @xmath224 .
since the spins of cfl are antiparallel in the considered static point , @xmath225 is independent of external magnetic field .
the above equation describes the critical current at which the cfl is destabilized as a rigid structure ( unaffected by external magnetic field along the @xmath22-axis ) .
numerical calculations presented below show that critical current is usually smaller than that given by eq.(12 ) . apparently , as shown by numerical simulations , there is a phase shift in initial precessions of @xmath50 and @xmath51 .
such a phase shift slightly perturbs initial antiparallel configuration and might reduce the critical current for the dynamics onset .    similar formula also holds for the opposite configuration ( @xmath119 ) , where the critical current is given by @xmath226\ , ,      \end{split}\ ] ] with @xmath227 .
now we discuss the theoretical results on critical currents in the context of those following from numerical simulations .
let us consider first the critical currents for individual spins of the syaf free layer , assuming that the second spin remains stable in its initial position , eqs .
( [ eq : critcurr1_du ] ) and  ( [ eq : critcurr2_du ] ) .
the corresponding results obtained from the formula derived above are presented in table  [ tab : crit ] , where we have omitted a weak dependence on @xmath21 .
for the studied structure with syaf free layer , @xmath212 is negative while @xmath213 is positive . from our analysis follows , that the current density at which dynamics appears in the simulations ( fig .
[ fig : saf_dynam1](a ) ) is higher than that given by @xmath212 . however , we checked numerically that the critical value @xmath212 agrees with the critical current obtained from simulations when assuming @xmath50 as free and fixing @xmath51 along @xmath96 .
.[tab : crit ] critical current densities in the units of @xmath228 , calculated according to equations  ( [ eq : critcurr1_du ] ) , ( [ eq : critcurr2_du ] ) , and  ( [ eq : iccfl ] ) for both syaf and syf free layers . [ cols="<,^,^,^,^ " , ]     following the above discussion of the cfl dynamics , one can understand the shift of the threshold current as follows .
initially , when the current density exceeds @xmath212 , @xmath50 becomes destabilized .
then , @xmath51 responses to the initial dynamics of @xmath50 with similar coherent precession .
however , @xmath51 should still be stable in its initial position at this current density and common precessions of the two coupled spin moments damps the initial dynamics .
accordingly , syaf ends up in the closest static state ( @xmath117 ) .
however , as the current density increases , the initial precessions of @xmath50 become more pronounced , which in turn means that the initial antiparallel configuration becomes distorted and @xmath51 becomes destabilized .
this results in coupled dynamics of both spins and finally leads to switching of the syaf structure .
on the other hand , we have also calculated the critical current for the whole syaf structure according to eq .
( [ eq : iccfl ] ) , and for the given structure we got @xmath229 shown in fig .
[ fig : saf_dynam1](a ) by the dashed vertical line ( see also table i ) .
equation  ( [ eq : iccfl ] ) describes stability of the whole cfl , and since the interlayer coupling is strong , @xmath229 corresponds to the current density at which both spins become destabilized simultaneously preserving their antiparallel orientation . as can be seen in fig .
[ fig : saf_dynam1](a ) , this is the maximum threshold current density for current - induced dynamics . because the rigid structure consisting of two antiparallel spins is not influenced by an external homogeneous magnetic field , there is no dependence of @xmath229 on @xmath21 .
nevertheless , from our numerical simulation follows that the threshold current for the syaf dynamics , @xmath230 , obeys the condition @xmath231 , provided that @xmath232 or ( as in our case ) @xmath213 has different sign .
when the syaf is in the @xmath119 configuration , the spin accumulation and spin current are different from those in the @xmath117 configuration ( at the same voltage ) .
this in turn leads to different spin torques , which is the reason why the critical currents destabilizing @xmath119 state are different from those for @xmath117 , as shown in table  [ tab : crit ] . from the critical currents one can expect
relatively symmetric hysteresis with applied current in structures with syaf .
in contrast , @xmath233 for the syf is negative , similarly as @xmath229 , but it is significantly larger , which indicates lack of hysteresis . to compare switching of the syaf and syf free layers from the @xmath117 to @xmath119 configurations with the opposite one ( @xmath119 to @xmath117 ) , we have simulated dynamics of the corresponding cfls assuming @xmath234 and varying current density .
the simulations have been performed in the quasistatic regime , i.e. , for each value of current density the spin dynamics was first equilibrated for @xmath110 and then averaged values of spin components and pillar resistance were calculated from the data taken for the next @xmath109 of dynamics . in order to prevent the system from collapsing into a static state with zero torque
, we have included a thermal stochastic field corresponding to @xmath157 [ see eq .
( [ eq : dth ] ) ] .
starting from @xmath155 and going first towards negative currents we have constructed the current dependence of the averaged resistance and related @xmath22-components of both spins , as shown in fig .
[ fig : hyster ] .
hysteresis loops of the resistance for the studied pillars with syaf ( a ) and syf ( d ) free layers .
panel ( b ) and ( c ) depict spin dynamics of @xmath50 and @xmath51 in syaf , respectivelly , corresponding to resistance loop ( a ) .
panels ( e ) and ( f ) show dynamics of @xmath50 and @xmath51 in syf , respectivelly , corresponding to resistance loop ( d ) .
the initial point of each hysteresis loop is marked with a dot .
the arrows indicate direction of the current change . figures ( g ) and ( h ) correspond to the upper parts of ( a ) and ( d ) , in which however the effects due to magnetostatic field of the reference layer to the cfl spins have been omitted . ]    for both syaf [ figs .
[ fig : hyster](a c ) ] and syf [ figs .
[ fig : hyster](d  f ) ] free layers , one can see relatively symmetric hysteresis with the current density . in both cases
direct switching from @xmath117 to @xmath119 state occurs at a current density comparable to @xmath229 .
in contrast , in the case of syf free layer , the second transition ( @xmath119 to @xmath117 ) appears at a current density which is very different from that predicted by the linearized llg model .
moreover , in both cases switching from @xmath119 to @xmath117 state does not appear directly , but through some precessional states . more precisely , as the positive current density increases , both spins start precessing in the layers plane prior to switching .
the in - plane precessions are connected with a significant drop in the resistance and with a reduction of the @xmath235-components .
the range of ipp regime is particularly large in the case of syf . from the analysis of spins trajectories
one may conclude that the angle of ipps increases with increasing current density , and after exceeding a certain threshold angle cfl switches to the @xmath117 configuration .
the other factor giving rise to the the difference in switching from @xmath119 to @xmath117 and from @xmath117 to @xmath119 follows from the fact that the magnetostatic interaction of the cfl s layers with the polarizer is different in the @xmath117 and @xmath119 states . to prove this we have constructed analogical hysteresis loops for syaf and syf free layers disregarding magnetostatic interaction with the @xmath1 layer ; see figs .
[ fig : hyster](g ) and ( h ) . for both syaf and syf free layers we observe now large decrease in @xmath127 for both switchings .
this implies that both switchings are realized _ via _ in - plane precessions , in contrast to the case when @xmath1 influences the cfl dynamics via the corresponding magnetostatic field .
while the hysteresis loop for syaf remains symmetric , the one for syf becomes highly asymmetric .
the asymmetry of syf loop is due to a significant asymmetry of stt in @xmath119 and @xmath117 states , which was previously shaded by the magnetostatic coupling with the layer @xmath1 .
we have studied current - induced dynamics of syaf and syf composite free layers . by means of numerical simulations we identified variety of dynamical regimes . the most significant difference between dynamics of syaf and
syf free layers concerns the evidence of self - sustained dynamics of both cfl spins . while in the case of syaf only coupled in - plane precessions in a narrow window of external parameters ( @xmath21 and @xmath55 )
are observed , syf free layer reveals more complex and richer dynamics , with the possibility of coupled out - of - plane precessions which might be interesting from the application point of view .
furthermore , as shown by numerical simulations , both syaf and syf are switchable back and forth without the need of external magnetic field . for syaf element two possible ways of switching have been identified .
since they lead to different switching times , their identification might be crucial for optimization of switching in real devices with syaf free layers .
however , one has to note that the diagrams shown in figs .
[ fig : saf_dynam1 ] and  [ fig : sf_dynam1 ] may be changed when magnetization in cfl becomes non - homogeneous .
a disadvantage of the studied structures is their relatively low efficiency of switching , i.e. high amplitude of critical current and long switching time . in order to show more sophisticated ways of tuning the cfl devices
, we have analyzed critical currents derived from the linearized llg equation .
the formula  ( [ eq : iccfl ] ) has been identified as the maximum value of critical current at which dynamics of the cfl structure should be observed .
this formula reveals some basic dependence of critical current on spin valve parameters , and therefore might be useful as an initial tool for its tuning .
however , in some cases non - linear effects in cfl dynamics might completely change the process of cfl switching , as shown by the presented numerical simulations . but
the effects of non - linear dynamics go beyond the simple approach of linearized llg equation , and their study requires more sophisticated nontrivial methods and/or numerical simulations .
this work was supported by polish ministry of science and higher education as a research project in years 2010 - 2011 , and partly by eu through the marie curie training network spinswitch ( mrtn - ct-2006 - 035327 ) .
the authors thank m. gmitra for helpful discussions .
one of us ( pb ) also thanks l. lpez daz , e. jaromirska , u. ebels , and d. gusakova for valuable suggestions .
torque acting on the left interface of @xmath3 is calculated from @xmath99 and @xmath98 components of @xmath236 , where @xmath237 is spin current vector in @xmath2 layer ( written in global frame ; shown in fig .
[ fig : scheme ] ) , and @xmath238 , where @xmath239 is the matrix of rotation by angle @xmath12 along the axis @xmath240 in the counterclockwise direction when looking towards origin of the coordinate system .
hence @xmath241 components can be written as      where @xmath243 are spherical coordinates of @xmath50 in the global frame .
similarly , we define torques amplitudes on the left interface of @xmath5 from the components of transformed spin current vector @xmath244 . in this case , however , @xmath245 is not written in the global frame , but in the local coordination system coordinate system connected with @xmath50 . to rotate local coordinate system of @xmath50 to local coordinate system of @xmath51 we need to know spherical angles @xmath246 and @xmath247 of vector @xmath51 in the local coordinate system of @xmath50 .
this might be done by transforming first @xmath51 vector to local coordinate system of @xmath50 as @xmath248 and calculate its angles @xmath246 and @xmath247 .
then we can calculate components of @xmath249 similarly as for the left interface      equation in @xmath4 , which is adjacent non - magnetic interface from the right - hand side of @xmath3 , are written in local coordinate system of @xmath50 . to apply the definition of @xmath251 and @xmath252 we need to rotate the local coordinate system so , that its @xmath98-axis will lie in the layer given by vectors @xmath50 and @xmath51 .
this might be done by single rotation of local coordinate system around its @xmath22-axis by angle @xmath253 , @xmath254 , where


binary motion is now well known to generate a spiral pattern in the circumstellar envelope of a mass losing star @xcite .
however , it is not generally recognized that the binary - induced pattern , vertically extended from the orbital plane , exhibits a ring - like pattern with an inclined viewing angle . at an inclination close to edge - on
, in particular , the binary - induced spiral - shell appears as eccentric arcs in the plane of the sky ( see e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
even at a low inclination ( close to face - on ) , the pitch angle of the spiral pattern decreases with radius from the center of binary orbital motion @xcite , such that the outer part of the spiral approximates a concentric ring structure .
this misidentification may occur , especially , for the outer ring - like patterns of ( pre-)planetary nebulae ( ppn / pn ) whose inner parts are contaminated by the bright , fast - moving lobes ( e.g. , ngc 7027 @xcite , egg nebula @xcite ) . because of the incompleteness of the observed ring - like patterns , whether they arise from spherical shells due to the pulsation of single stars or spiral - shell patterns due to binary orbital motion is inconclusive . in this paper
, we focus on the interpretation that the incomplete ring - like patterns observed in the past may be caused by the motion of ( unknown ) binary stars .    linking a binary scenario with the ring - like patterns observed on the circumstellar envelopes of asymptotic giant branch ( agb ) stars and their remnant material in the outer parts of ppn / pn has an interesting implication for the late phases of stellar evolution .
since the bipolar ( often multi - polar ) structures , including jets and tori , observed in the ppn / pn ( e.g. , * ? ? ?
* ; * ? ? ?
* ) likely originate from the central binary systems @xcite , a fundamental issue remains concerning the origin of the ring - like patterns in the outer parts of the ppn / pn .
such a ring - like pattern is commonly assumed as the consequence of the periodic mass loss of a _ single _ star in the previous agb phase , due to stellar pulsation ( e.g. , * ? ? ?
* ; * ? ? ?
* ) or thermal pulsation ( e.g. , * ? ? ?
* ) on the spherically symmetric wind .
for example , the archetypal agb star , irc+10216 , reveals the circumstellar pattern in the optical and infrared images @xcite , which is characterized by an interval corresponding to @xmath12501700 years .
interestingly , @xcite noted that this pattern interval of irc+10216 does not agree with either the stellar pulsation ( 649 days , * ? ? ? * ) or the thermal pulsation ( 600033000 years , * ? ? ?
similar discrepancy has also been found in several ppn / pn @xcite .
a binary - induced spiral - shell model , in which the pattern interval time is defined by the binary orbital period without a strong limitation on the scale , may possibly serve as an alternative explanation for the ring - like ( but not perfectly symmetric ) patterns formed in the agb phase .    in the binary model
, the observed spiral - shell pattern of a circumstellar envelope can constrain the binary mass and orbital properties . for the carbon - rich agb star , afgl 3068 , displaying a well - defined spiral circumstellar pattern in dust scattered light , @xcite estimated the orbital period simply from the arm interval and wind expansion speed , ignoring the inclination . for another carbon star , r sculptoris , exhibiting a spiral pattern in addition to a thin thermal pulse shell , @xcite characterized the thermal pulsation with the help of their modeled binary properties for the spiral pattern features . similarly , the circumstellar patterns in close binary systems may be used in determining the relative importance and observational consequence of the accretion disks surrounding the companion star , as theoretically investigated ( e.g. , * ? ? ?
* ; * ? ? ?
furthermore , @xcite found that the bow shock around a disk can lead to complicated time - dependent post - shock structures for orbital separations less than about 20 au .
as such , the features of a spiral - shell pattern provide indirect evidence for the presence of a binary companion to the agb star . in cases where the companion is obscured by the circumstellar envelope
, the detection of spiral features may be the only diagnostic probe for , otherwise hidden binary systems .    in a previous paper @xcite
, we took account of the effect of inclination angle , for the first time , in modeling an observed circumstellar pattern within the binary framework .
it was shown that the elongated spiral shape in the plane of the sky , seen in dust scattered light , can be used to constrain key binary quantities ( i.e. , inclination of the orbital plane , orbital period , companion mass , and mass ratio ) .
however , this analysis left a degeneracy in such model parameters since the apparent shape in the plane of the sky reflects both the binary motion and projection effect due to the inclination of the orbital plane . in this paper , we show that the degeneracy can be lifted by using gas kinematics from high resolution observations , which provide three - dimensional information .
to illustrate the modeling method , we apply the model to the object , cit 6 , by revisiting its incomplete ring - like pattern observed in molecular line emission with the spiral - shell due to the motion of ( hypothesized ) central binary stars .    in section [ sec : cit ] , we introduce the known properties of cit 6 and the existing molecular line observation with the very large array ( vla ) , which are used to constrain the model parameters . in section [ sec : geo ] , we formulate a simple analytic model based on the oblate spheroid geometry of an archimedes spiral - shell with a binary orbital inclination @xmath2 .
this simple analytic model can be easily applied to other sources in order to ascertain whether a spiral - shell pattern exists . in section [ sec : deg ] , a parameter space analysis is carried out following the method of @xcite based on the elongated shape of the pattern in the plane of the sky . in section [ sec : hyd ] , the results of hydrodynamic radiative transfer simulations are presented , from which we use molecular line kinematics to resolve the degeneracy in the parameter space analysis .
we find that the simulated spiral - shell model can reproduce the observational data in the channel maps and position - velocity ( p - v ) diagrams , providing a better fit than obtained by using a spherically symmetric shell model . in section [ sec : dis ]
, we discuss the advantage and limitation of the current model .
finally in section [ sec : sum ] , we summarize our findings and suggest future applications .
cit 6 ( also known as rw lmi , irc+30219 , iras 10131 + 3049 , afgl 1403 ) is an extreme carbon star and believed to be in the transition from the agb to the post - agb phase .
@xcite found the presence of a nascent bipolar nebula , providing evidence that the evolutionary phase of cit 6 lies just pass the tip of agb .
its effective temperature of 2800k @xcite and the bolometric luminosity of @xmath3 @xcite also correspond to the values at the tip of agb phase of the evolutionary track for a star with the initial mass of 23 , according to the single stellar evolution package of @xcite assuming a solar metallicity .
the distance to cit 6 is estimated to be @xmath4pc based on the pulsation period ( @xmath5 yr ) using a period - luminosity relation ( * ? ? ? * and references therein ) .
the importance of cit 6 stems from the existence of partial rings of molecular line emission based on vla interferometric observational data obtained by ( * ? ? ?
* see also the red - scaled image in figure[fig : chm ] ) at high angular resolution ( @xmath6 ) and high spectral resolution ( @xmath7 ) .
@xcite modeled this line emission pattern with four spherical shells , using the shell radii , expansion velocities , and central velocities as free parameters for the individual shells .
their model fitting results in different central velocities between the modeled shells , which may imply the departure from spherical symmetry of the circumstellar pattern about the central agb star .
@xcite remarked on the existence of the asymmetric features such as the spiral shape ( relatively well traced in the velocity channel of @xmath8 ) , stronger emission observed in the redshifted side , and lack of emission in the west and northwest near the central channels .
@xcite showed that the optical depth effect can not reproduce the asymmetric line profile of cit 6 for a physical structure assumed to be spherically symmetric .
their three - dimensional morphokinematic model instead suggested several incomplete shells for the structure of cit 6 .
here we suggest an alternative interpretation for the asymmetric shape of the circumstellar pattern as an indication for a central binary star system .
some evidence for binarity is indeed found in the hubble space telescope ( hst ) images with two compact cores separated by @xmath9 at the position angle ( pa , measured from north to east ) of 10  @xcite .
the northern lobe varies in its flux and polarization with the known pulsation period of the carbon star ( @xmath5 yr ; @xcite ) , while the bluer southern lobe is stable and has a polarization direction @xmath10 different from the northern red component @xcite .
these characteristics of the southern lobe , differing from those of the carbon star and the northern lobe , suggest that the blue southern lobe is illuminated by a companion star in the vicinity @xcite .
the existence of a companion star is supported by the consistence of the polarization direction of h@xmath11 emission @xcite to the blue component . from the observed flux at 4500
, @xcite inferred that the companion is likely later than spectral type g0 ( @xmath12 ) if it is a main - sequence star .
the orbital motion of a mass losing star in a binary system creates a spiral - shell pattern in its circumstellar envelope ( e.g. , * ? ? ?
we consider the overall morphology of the spiral - shell pattern as an oblate spheroid , following @xcite . in the coordinates having the vertical axis @xmath13 aligned with the orbital axis
, the oblate spheroid can be described by ( see figure[fig : geo ] ) @xmath14 where the axial ratio @xmath15 is defined by the ratio of pattern propagation speeds @xmath16 in the orbital plane and along the polar axis . here , the radial propagation speed of the spiral pattern in the orbital plane is affected by the orbital motion of the mass losing star with the orbital speed of @xmath17 to be @xmath18 @xcite , while the pattern propagation speed along the orbital axis is the average wind speed @xmath19 corresponding to the intrinsic wind speed .
we define new coordinates @xmath20 rotated by an inclination @xmath2 , in which the observer is at infinity in the direction of the @xmath21-axis .
a rotational transformation about the @xmath22-axis by the angle @xmath2 is defined by @xmath23 , @xmath24 , and @xmath25 .
substitution of these into equation ( [ eqn : obl ] ) yields @xmath26 where @xmath27 and @xmath28 refer to the tangential coordinate @xmath29 and the azimuthal angle in the @xmath22@xmath30 plane , respectively , defining @xmath31 and @xmath32 .
a factor @xmath33 is defined as @xmath34 , thus @xmath35 alternative forms of equation ( [ eqn : inc ] ) are @xmath36 or @xmath37 because the oblate spheroid geometry originates from the latitudinal variation of the wind velocity , these coordinates @xmath27 and @xmath21 can be linearly transformed to the tangential and line - of - sight velocities of the wind .
thus , the tangential component is @xmath38 representing the pattern propagation speed in the inclined plane within the @xmath39-channel .
the second term in the right - hand - side of equation ( [ eqn : v_t ] ) shows the inclination effect , and the last term is related to the channel velocity representing the height of the relevant material from the midplane . from equation ( [ eqn : v_t ] ) , the line - of - sight velocity can be written as @xmath40 at a given @xmath41 .
the maximum line - of - sight velocity is @xmath42 at the orbital center ( i.e. , at the center of mass in the circular orbit case ) .
a circumstellar spiral pattern due to the orbital motion of the agb star ( orbital radius @xmath43 ; orbital velocity @xmath17 ) for a given @xmath21 is defined by a differential equation in the polar coordinates @xmath44 as , @xmath45 where the tangential velocity @xmath41 represents the pattern propagation velocity in the layer within the @xmath39-channel . the velocity @xmath41 is given by equation ( [ eqn : v_t ] ) as a function of @xmath28 with given wind velocity , orbital velocity , and inclination angle . as the velocity of an agb wind is nearly constant beyond the wind acceleration region ( typically less than several times the stellar radius ; @xcite ) , the integrated form of equation ( [ eqn : dif ] ) well approximates to @xmath46 for the pattern located at a large radius @xmath27 . here
, @xmath47 is the distance to the source .
in order to obtain approximate constraints on the binary properties , we explore the possible parameter space for the system of interest , following the method described in @xcite using an overall oblate spheroidal morphology of the pattern .
four observed quantities are required , corresponding to the projected axial ratio @xmath48 , the arm pattern spacing @xmath49 along the major axis , the projected binary separation @xmath50 , and the angular position of the stars relative to the major axis of the elliptical pattern shape @xmath51 . here
, @xmath52 refers to a position angle measured from the north to the east direction and hereafter @xmath28 is defined as the angle relative to the line of nodes of the binary orbit ( @xmath53 ; the intersection of the orbital plane with the plane of the sky ) measured in the counterclockwise direction .
refer to @xcite for the details of this method using an overall oblate spheroid shape .
we first identify the line of nodes of the binary orbit as the major axis of the elongated spiral pattern of the hc@xmath0n(43 ) line emission in figure[fig : chm ] .
because of the insufficient sensitivity of observation , it is challenging to precisely determine the angular direction of the major axis only by ellipse fitting .
thus , we also use the fact that the p - v diagram is most symmetric along the line of nodes and most asymmetric in the perpendicular direction ( middle and right panels of figure[fig : pvd ] , respectively ; see also section [ sec : inc ] ) . from the resulting ellipse
fit with the major axis at @xmath54 , we measure the axis ratio @xmath48 of 1.15 and the pattern spacing @xmath49 of @xmath55 along the major axis .
we assume the two compact components from the optical images are members of a binary ( or their surrounding gaseous envelopes ) following @xcite , which provides a projected separation of @xmath56 and a position angle of the binary at @xmath57 .
the line of nodes nearly coincides with the alignment of the two stars , implying that @xmath50 is approximately equivalent to the actual binary separation .
in addition to these four observed quantities , the distance of 400pc is adopted ( see section [ sec : cit ] ) , thus constraining the binary properties in a parameter space spanning the average wind velocity @xmath19 versus the orbital velocity @xmath17 of the mass losing star ( figure[fig : deg ] ) .
the gray - colored area in figure[fig : deg ] shows the possible sets of binary parameters .
the first constraint for the gray area is the mass of the evolved carbon star ( @xmath58 , black solid line ) greater than the mean mass of white dwarfs @xcite .
in particular , the stellar luminosity of @xmath3 for cit 6 suggests the core mass of @xmath59 with 2% accuracy based on @xcite , providing the lower limit on the current mass , @xmath58 , of the agb star .
the other constraint stems from the maximum line - of - sight velocity of the pattern @xmath60 corresponding to the half width of the hc@xmath0n(43 ) line , measured as 18  with an uncertainty of @xmath61 considering the spectral resolution . in previous papers @xcite , the lower limit of the maximum line - of - sight velocity
was loosely constrained by the condition @xmath62 , which is satisfied at any inclination .
we further calculate the maximum line - of - sight velocity taking account of the inclination effect . from equation ( [ eqn : vrm ] ) , the measured line - of - sight velocity range @xmath63 is marked in figure[fig : deg ] as a set of yellow dotted lines .
the plausibility of a particular model parameter set can be easily checked by the simple analytic model in section [ sec : geo ] . among the possible models in figure[fig :
deg ] , the parameter set marked by a plus symbol is applied to the analytic model .
the model consists of an estimate of the agb star mass of @xmath64 , companion mass @xmath65 , binary separation @xmath66au , inclination @xmath67 , and average wind velocity @xmath68 with the distance @xmath69pc .
the result is overlaid in each velocity channel map of figure[fig : chm ] by a gray - colored solid line , showing a reasonable match with the outer edge of the observed pattern ( red - scaled ) , particularly for the middle - row panels displaying the near - central velocity channels .
furthermore , this inclined spiral can explain the wiggle feature of the observed pattern in the polar coordinates , clearly seen in figure[fig : pol ] . without considering the inclination effect ( i.e. , assuming face - on ; @xmath70 ) , the pattern propagation speed @xmath41 in each @xmath39-velocity channel is constant in equation ( [ eqn : v_t ] ) , indicating that the pattern should appear as a straight line in the polar coordinates with the slope determined by equation ( [ eqn : pol ] ) . for comparison ,
a concentric spherical shell pattern forms horizontal lines in the @xmath28-@xmath27 coordinates in figure[fig : pol ] .
therefore , the wiggle feature , which is also found in other sources ( e.g. , r sculptoris , see figure 2a in @xcite ) , is a natural outcome of the inclined spiral model , and is distinguishable from the face - on view of the spiral model for a binary system and the spherical shell model for a single star .
the region delineated by the gray - colored area in the parameter space ( figure[fig : deg ] ) , based on the elongated spiral shape projected on the plane of the sky , can be sharpened by taking advantage of three - dimensional information provided by the molecular line kinematics .
we explore the parameter space in the vicinity of the @xmath71  line with hydrodynamic simulations followed by radiative transfer calculations for the hc@xmath0n molecular line .
the parameters for the model which best matches the aspects of the vla observation @xcite , described in sections [ sec : inc ] and [ sec : res ] , are summarized in table[tab : par ] ( and by a plus sign in figure[fig : deg ] ) .
we note that these parameters are not independent of each other due to their constraints imposed by the parameter space analysis .
the basic equations of hydrodynamics are integrated using the flash3 code @xcite based on a piecewise parabolic method @xcite .
we construct a centrally concentrated , static grid utilizing the block - structured adaptive mesh refinement paramesh package and adjusting the refinement level based on the distance from the stars .
the maximum refinement level of eight with the grids of @xmath72 per block is used for the simulation domain of 8400au@xmath738400au@xmath734200au ( @xmath74 at 400pc , corresponding to the image size of the vla observation in consideration ) in three - dimensional cartesian coordinates , assuming the mirror symmetry about the orbital @xmath22@xmath30 plane coordinates defined in section [ sec : geo ] and figure[fig : geo ] .
the prime notation is omitted from now on . ] .
the spatial gridding corresponds to a resolution of 1au ( @xmath75 ) in the central region , assigning a sufficient number of grid cells at the wind generating radius @xmath76 within which physical conditions are reset every simulational timestep , in order to assure the ejection of gas from the mass losing star in a spherical shape .
the gas temperature over the simulation domain is scaled by the temperature @xmath77 at the wind generating radius @xmath76 from the mass losing star .
we set @xmath78 , and @xmath77 is calculated by the temperature distribution @xmath79 of an adiabatic gaseous medium with a @xmath80 density distribution .
the observed effective temperature @xmath81 is set at the assumed stellar photospheric radius @xmath82 or 2.59au .
an adiabatic equation of state with the adiabatic index @xmath83 of 1.4 is used .
the gravitational potentials of the stars are treated as plummer potentials ( see , e.g. , * ? ? ?
* ) using the pointmass implementation of the flash3 code with modifications to include the gravitational softening radii and the motion of the two stars . the gravitational softening radius of the mass losing star does not change the result as it is smaller than the wind generating radius .
the gravitational softening radius of the companion is set to be comparable to the hoyle - lyttleton accretion radius .
an additional simulation with the companion s softening radius reduced by a factor of 2 is compared with the model presented in this paper , in order to verify its insignificant contribution to the overall shape of the arm pattern .
the accretion disk of the companion is expected to have a negligible size if it is gauged by the impact parameter ( 0.004 au with the large separation of our binary system ) in @xcite .
resolving the accretion disk for a wide binary system is computationally very expensive and beyond the scope of this paper .
the flash3 results are used as inputs to the radiative transfer calculation for hc@xmath0n lines with the sparx code .
sparx , standing for simulation package for astrophysical radiative transfer , is a software package for calculating molecular excitation and radiative transfer of ( dust ) continuum and molecular lines at millimeter / submillimeter wavelengths .
it adopts the accelerated monte carlo ( amc ) algorithm developed by @xcite , and finds full non - local thermodynamical equilibrium ( non - lte ) solutions of molecular level population and radiation fields iteratively and consistently . instead of assuming that the molecular excitation is simply in equilibrium corresponding to the local temperature ,
non - lte demands only statistical equilibrium of the level population through collisional and radiative excitation and de - excitation processes , which can be solved with the detailed balance equation .
additional ( de-)excitation mechanisms , such as infrared radiative pumping , have not been incorporated in the package .    as a post - processing tool , the sparx requires inputs of molecular gas density , temperature , velocity as well as the linewidth distribution , which can be provided by the results of a hydrodynamic simulation .
knowledge of the molecular abundance distribution and parameters such as dipole moments , line frequencies , and collisional rates are also required by the sparx .
the final sparx output includes synthetic images which can be used for generating spectra and p - v diagrams , as input for a telescope filter for simulated observations , or for other scientific analysis .
details of the package such as the code implementation and its verification with benchmark cases will be presented in liu et al .
( 2013 , in preparation ) .    for the sparx calculation
, the hc@xmath0n molecular parameters are obtained from the leiden atomic and molecular database ( lamda ) @xcite .
the hc@xmath0n abundance of @xmath84 is adopted throughout the envelope @xcite , while a spherical central region of 2   was carved out with hc@xmath0n abundance set to zero in order to mimic the central hole seen in the observed emission of hc@xmath0n ( also found in its parent molecule cn ; @xcite ) . for computational efficiency , we use the spatial and spectral resolutions of @xmath85 and 0.4  in the sparx calculation , which are sufficient to simulate the vla observation of @xcite summarized in section[sec : cit ] . to further compare the sparx output with the vla observation
, we perform a convolution with a gaussian beam of @xmath86 and the rebinning to 1  over the synthetic images .
the inferred binary parameters are sensitive to the orbital inclination as illustrated in figure[fig : deg ] .
hence , determination of the inclination can lift the degeneracy that remain in the parameter space analysis based only on the pattern shape in the plane of the sky . with the binary parameters in table[tab : par ] , figure[fig : inc ] compares the images for the non - zero velocity channels and the p - v diagrams along the major ( i.e. , @xmath22-axis ; line of nodes ) and minor ( i.e. , @xmath30-axis ; perpendicular to the line of nodes ) axes for five different inclinations of 0  ( face - on ) , 30 , 50 , 70 , and 90   ( edge - on ) .
the channel images ( top row ) illustrate the change of pattern shape from a spiral in the face - on view to ellipses in the edge - on view .
in particular , the spiral pattern is split along @xmath87-axis in the image for @xmath88 and attached to the adjacent component in the image for @xmath89 , making a new ring .
the p - v diagrams of figure[fig : inc ] ( middle and bottom rows ) reveal notable features .
the middle row show that , regardless of inclination , the p - v diagram along the line of nodes ( i.e. , @xmath22-axis ) appears relatively symmetric about the central velocity because of the mirror symmetry of the overall oblate spheroid geometry about the orbital plane .
in contrast , the p - v diagram along the minor axis highly depends on the inclination angle ( bottom row ) .
one important implication is that the pa dependence of the p - v diagram can be used to determine the line of nodes of the binary orbit ( as in section [ sec : deg ] ) .
moreover , the comparison between the p - v diagrams along major and minor axes helps to infer the inclination angle of the system from the observed data cube .
for example , the pa dependence of the p - v diagram is minimized in the face - on case .
the p - v diagram along the minor axis ( bottom row ) provides further information on the orbital inclination . at a mid - range inclination ( see e.g. @xmath90 ) , it reveals an asymmetry in intensity between the redshifted and blueshifted components , the degree of asymmetry of which differs with position , and thus the overall shape of p - v diagram is tilted .
hence , the degree of tilt and spectral asymmetry is dependent on the inclination . on the other hand , the maximum line - of - sight velocity ( i.e. , @xmath39 at @xmath91 ) increases with inclination as predicted in equation ( [ eqn : vrm ] ) .
this increase of maximum velocity makes the outline of p - v diagram at the high end velocity at each @xmath30 position ( @xmath92 ) shallower in the face - on case ( @xmath93 ) and becoming steeper as the inclination increases .
all these features of the p - v diagrams can be used to constrain the inclination angle of the binary system .    in the framework of the binary model ,
the companion creates a gravitationally - induced density wake in a flatter spiral shape , overlapped with the spiral - shell pattern due to the orbital motion of the agb star @xcite .
the overlapped region appears as knots in the plane of the sky at the pa depending on the inclination ( see figure 10 of * ? ? ?
* ) because the knots represent the companion s wake relatively thin in the vertical direction from the orbital plane @xcite . in figure[fig : inc ] , the top panels for a non - zero velocity channel show such knots at high inclinations ( @xmath94 in this model ) at the pa varying with the inclination ( marked by arrows ) . in the p - v diagram along the minor axis ( bottom row )
, knots appearing with low inclinations ( @xmath95 ) also shift the positions depending on the inclination . on the other hand , in the p - v diagram along the line of nodes of orbits ( middle row )
all knotty structures are present at the zero velocity , confirming the existence of the companion s wake on the orbital plane .    in short , the shape of p - v diagram , the degree of spectral asymmetry , the position of knots both in individual channels and in individual p - v diagrams provide important clues in identifying the best model .
it is also worth noting that the velocity dispersion of individual arm component is dominated by the velocity variation due to the binary orbital motion ( @xmath96 ) . for the binary system in table[tab : par ] with an inclination angle of 60 , the above features reasonably match the vla observation .
figure[fig : chm ] shows the modeled hc@xmath0n(43 ) line emission with the parameters in table [ tab : par ] in green color , overlaid with the observational data displayed in red color .
a good correlation exhibits between the observed and modeled data in the individual velocity channels , albeit the outermost arm in the observation is barely detected at the @xmath97 level , where the noise level @xmath98mjy / beam refers the root - mean - square of the line - free channels in the vla observation .
the observed channel images ( red ) range between 1 and 5mjy / beam . to present the pattern shape beyond the sensitivity limit of the observation
, the modeled channel images ( green ) are displayed from 0.1 to 2mjy / beam .
all arc patches in the observed line emission are well located on the modeled spiral pattern . at the central channels ( middle panels of figure[fig : chm ] )
the model reveals elongated spiral patterns , while at higher velocity channels ( top and bottom panels ) ring - like patterns are reproduced . here , we emphasize that the ring - like features appear in the spiral - shell models with considerably large inclination angles ( see * ? ? ?
it suggests that ring - like patterns observed in the agb circumstellar envelopes and in the outer parts of ppn / pn do not necessarily imply an underlying spherically symmetric nature .
the inclination angle of the presented model is 60 .
the binary model in figure[fig : chm ] shows the breakup of rings and their unequal spacing between adjacent arms .
this appears in the intermediate range of inclination angles , reflecting the transition from a spiral shape to a ring shape , as explained in section [ sec : inc ] . in the @xmath8 channel , for example
, the ring - like pattern is incomplete revealed by the gap in the east at @xmath99 ( along the east side of minor axis ; see also the first panel of figure[fig : pvd ] ) .
the ring radius of the northern part of the gap is larger than the radius of the southern part , in which the model and observation match sensibly .
the companion s gravitational wake appears as knots in the modeled channel images located in specific pas ( marked by arrows in figure[fig : chm ] ) as determined by the inclination . in the zero - velocity channel ( middle panel of figure [ fig : chm ] ) representing the midplane , the knots are arranged along the line of nodes of the orbits ( i.e. , @xmath100 and 190 ) .
the knots shift to the east in the redshifted velocity channels and to the west in the blueshifted channels .
the direction of shift of the knots with the channel velocity is determined by the orientation of the companion s orbit .
the model knots eventually converge on the east and west ends of the ring pattern at the maximum velocity channels , which results in the thickening and brightening of the pattern at such locations .
the pas of knots in the modeled channel images are consistent with bright emission spots in the corresponding vla images .
for example , the observed image ( in figure[fig : chm ] ) exhibits thickening and brightening features in the ring pattern on the east and west ends at the @xmath101 and @xmath102 channels , respectively .
another example is found at the southwest in the @xmath103 and @xmath8 channels , where the detected emission in the outermost ring ( red ) corresponds to the position of bright knotty structure in the model ( green ) . at the locations of these knots , the arm ( in red and green images )
is found considerably inside the analytic archimedes spiral model without considering the companion s gravitational wake ( gray solid line ) .
this decrease of arm radii is explained by ( * ? ? ?
* figure 9 and text ) , which showed that a shock front ( appearing as knots in the channel images ) at the location of the companion s wake reduces the speed of the downstream flow and , therefore , the pattern propagation speed .
we have shown that a numerical hydrodynamic radiative transfer model for cit 6 demonstrates that the observed incomplete ring patches can be naturally understood within the framework of a single spiral - shell pattern induced by a binary motion .
in contrast , the previous spherical shell model @xcite and incomplete shell model @xcite required the use of individual shell parameters without providing physical understanding for their variations .
moreover , the binary model reveals that the knotty structures result from the overlapping between a spiral - shell shaped pattern due to the orbital motion of the mass losing star and a more flattened spiral pattern onto the orbital plane due to the companion s motion . in this picture , the knots are aligned in a specific pa , which systematically changes with the velocity channel .
this implies that such regularly - aligned knots are a natural consequence of the model .
the total flux integrated over the full simulation domain is related to the mass loss rate and the hc@xmath0n abundance for a given temperature distribution .
we have adopted the mass loss rate of @xmath104 and relative molecular abundance of @xmath84   from @xcite . unlike @xcite , who assumed the local thermal equilibrium with a single temperature , the temperature distribution in our model is determined by an adiabatic gas equation of state . the resulting total flux is smaller than the vla observation by a factor of 2 .
this implies that a complete modeling of these parameters requires studies involving a further fine - tuning of the input parameters as well as the introduction of molecular chemistry .
such studies will also necessitate incorporating observations of various molecular species , transitions , and angular resolutions .
however , these parameters merely scale the flux level , but do not change the shape of the spiral pattern which is the focus of this paper .
the fact that the observed redshifted emission component is stronger than the blueshifted component is , however , not reproduced in our spiral model , albeit a small magnitude of spectral asymmetry is found depending on the inclination .
this implies that our model is incomplete and that the asymmetry along the line - of - sight does not originate solely from the binary motion
. the observed spectral asymmetry may agree with a presence of a bow shock due to the interaction between the circumstellar and interstellar media if the space motion away from us is significant ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
however , the vla data does not show a shift of the local standard of rest ( lsr ) velocity from the central velocity of the spiral - shell pattern .
other possible effects which may enhance the asymmetry include physical and chemical inhomogenities in the surrounding matter or in the dust formation phase at the onset of the stellar wind ( e.g. , * ? ? ?
* ; * ? ? ?
* ) if propagated to large scales .
the missing emission in the west side of the central channels is also beyond the scope of this paper .
this observed feature can not be explained simply by a spiral - shell model nor by a spherical shell model .
this spatial asymmetry is also found in a different transition of hc@xmath0n and in its parent molecule cn @xcite .
we speculate that such lack of emission could possibly be caused by a chemical effect reducing the molecular abundance , rather than being affected by the gas density .    compared to the other model sets in the parameter space determined by the projected pattern shape ( figure[fig : deg ] ) , for instance at @xmath105
, the presented model shows a better match with the observation for the pattern features explained in section [ sec : inc][sec : res ] in both channel maps and p - v diagrams .
however , the difference of our fiducial model from the other comparison model sets are not decisive for the given uncertainties in the model parameters .
in particular , observational uncertainties due to the limited sensitivity in determining the elongated shape and the arm spacing , as well as the determination of the source distance , stellar luminosity , photospheric temperature , and binary separation enter .
follow - up observations with a higher sensitivity for a chemically stable molecule such as co would address some of these uncertainties .
we have developed a method to investigate the incomplete ring - like feature of a spiral - shell pattern in an outflowing circumstellar envelope , due to the orbital motion of central binary stars with an inclination with respect to the orbital plane , providing an alternative explanation to models invoking periodic mass loss events of a single star .
as a first step in this method , a parameter space study is performed as in @xcite and improved upon by taking account of the inclination angle for the constraint on the wind expansion velocity .
the results of such an analysis significantly narrow down the parameter space that can provide a fit to the observed pattern shape in the plane of the sky .
this reduction in the parameter space is important for the second step where a modest number of hydrodynamic and radiative transfer investigations can be carried out .
the pattern shape , degree of asymmetry , position of knots in individual velocity channels and p - v diagrams are used via detailed modeling to constrain binary properties .
we have shown that gas kinematics is important in constraining the binary properties , which would otherwise remain degenerate under the constraints from the parameter space analysis based only on the shape of the circumstellar spiral - shell pattern in the plane of the sky .
we have also demonstrated that a simple analytic analysis for the elongated spiral model based on an archimedes spiral - shell pattern viewed as a function of inclination angle for comparison to individual velocity channels is useful for determining the applicability of the model to observational data .
finally , observed circumstellar patterns often reveal wiggles with varying slope in the polar coordinates that can be adequately described by the elongated spiral model , not by a face - on archimedes spiral nor by spherical shells .
we have used the agb star , cit 6 , as the first example to illustrate a method which can be used to help to constrain the binary properties given the existence of high resolution and sensitivity vla data for the hc@xmath0n(43 ) molecular line .
this first success of applying the binary - induced spiral - shell model to an observed incomplete ring - like pattern opens the possibility of connecting the ring - like patterns commonly found in the agb circumstellar envelopes and in the outer parts of ppn / pn and pointing to the conceivable presence of central binary systems .
we acknowledge a detailed comments from an anonymous referee , and thank dr .
mark claussen for providing the calibrated fits image of the vla data .
this research is supported by the theoretical institute for advanced research in astrophysics ( tiara ) in the academia sinica institute of astronomy and astrophysics ( asiaa ) .
the computations presented here have been performed through the asiaa / tiara computing resource .
the sparx code development is supported by the coordinated / coherent hydrodynamic and astrophysical research , modeling and synthesis ( charms ) project in asiaa .
the flash code is developed by the doe - supported asc / alliance center for astrophysical thermonuclear flashes at the university of chicago .
balick , b. , & frank , a.  2002 , , 40 , 439 binney , j. , & tremaine , s.  2008 , galactic dynamics ( princeton , nj : princeton univ . press ) bond , h.  1996 , hst press release stsci - pr96 - 05 colella , p. , & woodward , p.  r.  1984 , journal of computational physics , 54 , 174 chau , w. , zhang , y. , nakashima , j .-
i . , deguchi , s. , & kwok , s.  2012 , , 760 , 66 claussen , m.  j. , sjouwerman , l.  o. , rupen , m.  p. , et al .
2011 , , 739 , l5 cohen , m.  1979 , , 186 , 837 cohen , m. , & hitchon , k.  1996 , , 111 , 962 decin , l. , royer , p. , cox , n.  l.  j. , et al .
2011 , , 534 , a1 de marco , o.  2009 , , 121 , 316 dinh - v .- trung , & lim , j.  2009 , , 701 , 292 edgar , r.  g. , nordhaus , j. , blackman , e.  g. , & frank , a.  2008 , , 675 , l101 fryxell , b , .
olson , k. , ricker , p. , et al .
2000 , , 131 , 273 he , j.  h.  2007 , , 467 , 1081 hfner , s.  2007 , why galaxies care about agb stars : their importance as actors and probes , 378 , 145 hogerheijde , m.  r. , & van der tak , f.  f.  s.  2000 , , 362 , 697 huggins , p.  j.  2007 , , 663 , 342 huarte - espinosa , m. , carroll - nellenback , j. , nordhaus , j. , frank , a. , & blackman , e.  g.  2013 , , 1448 hurley , j.  r. , pols , o.  r. , & tout , c.  a.  2000 , , 315 , 543 kim , h.  2011 , , 739 , 102 kim , h. , & taam , r.  e.  2012a , , 744 , 136 kim , h. , & taam , r.  e.  2012b , , 759 , 59 kim , h. , & taam , r.  e.  2012c , , 759 , l22 leo , i.  c. , de laverny , p. , mkarnia , d. , de medeiros , j.  r. , & vandame , b.  2006 , , 455 , 187 ladjal , d. , barlow , m.  j. , groenewegen , m.  a.  t. , et al .
2010 , , 518 , l141 le bertre , t.  1992 , , 94 , 377 lindqvist , m. , schier , f.  l. , lucas , r. , & olofsson , h.  2000 , , 361 , 1036 loup , c. , forveille , t. , omont , a. , & paul , j.  f.  1993 , , 99 , 291 madej , j. , naleyty , m. , & althaus , l.  g.  2004 , , 419 , l5 maercker , m. , mohamed , s. , vlemmings , w.  h.  t. , et al .
2012 , , 490 , 232 mastrodemos , n. , & morris , m.  1999 , , 523 , 357 mattsson , l. , hfner , s. , & herwig , f.  2007 , , 470 , 339 mauron , n. , & huggins , p.  j.  2006 , , 452 , 257 mayer , a. , jorissen , a. , kerschbaum , f. , et al .
2011 , , 531 , l4 monnier , j.  d. , tuthill , p.  g. , & danchi , w.  c.  2000 , , 545 , 957 morris , m.  1987 , , 99 , 1115 nordhaus , j. , & blackman , e.  g.  2006 , , 370 , 2004 paczyski , b.  1971 , , 21 , 271 sahai , r. , morris , m.  r. , & villar , g.  g.  2011 , , 141 , 134 sahai , r. , trauger , j.  t. , watson , a.  m. , et al .
1998 , , 493 , 301 schmidt , g.  d. , hines , d.  c. , & swift , s.  2002 , , 576 , 429 schier , f.  l. , van der tak , f.  f.  s. , van dishoeck , e.  f. , & black , j.  h.  2005 , , 432 , 369 soker , n.  1994 , , 270 , 774 soker , n.  2000 , , 312 , 217 su , k.  y.  l.  2004 , asymmetrical planetary nebulae iii : winds , structure and the thunderbird , 313 , 247 theuns , t. , & jorissen , a.  1993 , , 265 , 946 trammell , s.  r. , dinerstein , h.  l. , & goodrich , r.  w.  1994 , , 108 , 984 trammell , s.  r. , & goodrich , r.  w.  1997 , aas , 191 , 47.15 wareing , c.  j. , zijlstra , a.  a. , speck , a.  k. , et al .
2006 , , 372 , l63 wilkin , f.  p.  1996
, , 459 , l31 willson , l.  a. , & hill , s.  j.  1979 , , 228 , 854 woitke , p. , & niccolini , g.  2005 , , 433 , 1101 wood , p.  r.  1979 , , 227 , 220 zhang , y. , kwok , s. , & dinh - v - trung 2009 , , 691 , 1660