the discovery of persistence has recently generated considerable interest in understanding the statistics of first passage problems in spatially extended nonequilibrium systems , both theoretically and experimentally .
the definition of persistence is as follows .
let @xmath11 be a stochastic variable fluctuating in space and time according to some dynamics .
the persistence probability is simply the probability @xmath12 that at a fixed point in space , the quantity @xmath13 $ ] does not change up to time @xmath7 . in many systems of physical interest a power law decay , @xmath14 , is observed , where @xmath15 is the persistence exponent and is , in general , nontrivial .
the nontriviality of @xmath15 emerges as a consequence of the coupling of the field @xmath11 to its neighbours , since such coupling implies that the stochastic process at a fixed point in space and time is non - markovian .
persistence phenomena have been widely studied in recent years @xcite .
theoretical and computational studies include spin systems in one @xcite and higher @xcite dimensions , diffusion fields @xcite , fluctuating interfaces @xcite , phase - ordering dynamics @xcite , and reaction - diffusion systems @xcite .
experimental studies include the coarsening dynamics of breath figures @xcite , soap froths @xcite , and twisted nematic liquid crystals @xcite .
persistence in nonequilibrium critical phenomena has also been studied in the context of the global order parameter , @xmath16 , ( e.g.  the total magnetization of a ferromagnet ) , regarded as a stochastic process @xcite .    in the present work we consider spatially extended systems with a nonequilibrium field @xmath11 , which takes discrete values , at each lattice site @xmath17 .
the field then evolves in time @xmath7 through interactions with its neighbours .
the persistence probability at time @xmath7 is defined as the fraction of sites in which the stochastic field @xmath11 did not change its value in the time interval @xmath18 $ ] .
the physical interpretation of @xmath11 could be for , for example , the coarsening spin field in the ising model after being quenched to low temperature from an initial high temperature , the sign of a diffusion field starting from a random initial configuration or the sign of the height , relative to the mean height , of a fluctuating interface . as the stochastic field evolves in time
, such systems develop regions of persistent and nonpersistent sites .
since the number of persistent sites decays with time according to @xmath6 , the persistent clusters shrink in size and number with a corresponding growth in the size of the nonpersistent regions .
recently , manoj and ray ( mr ) @xcite have studied such spatially extended systems in one dimension ( 1d ) within the context of the @xmath19 reaction diffusion model , which is equivalent to the 1d ising model .
they found that the length scale which characterises the interval sizes between persistent clusters _ apparently _ has a different time - dependence from the length scale characterising the walker separations ( or spin correlations , in the ising representation ) , and futhermore seems to depend on the initial walker density .
if true , this result would be surprising .
naively , one would expect the initial walker density to be irrelevant to the asymptotic dynamics , and the coarsening length scale , set by the mean distance between walkers , to be the relevant length scale for all spatial correlations .
we shall show that the former expectation is correct , but the latter , in general , is not .    in this paper
we generalise , and reinterpret , this study in the context of the @xmath2-state potts model , which has @xmath2 distinct but equivalent ordered phases .
the @xmath19 model corresponds to the @xmath20 potts model , i.e.  the ising model , with the walkers identified as the domain walls between up and down ising spins . in the @xmath2-state potts model
the random walkers represent domain walls between potts states . at each time step , every walker hops randomly left or right .
persistent potts sites are those at which the potts state has never changed , i.e.   sites which have never been crossed by a walker .
if two walkers occupy the same site they either _ coalesce _ , with probability @xmath21 , or _
, with probability @xmath22 , these numbers being the probabilities that the states on the farther sides of the walkers are the same ( annihilation ) or different ( coalescence ) .
the persistence probability ( the fraction of sites that have never been jumped over by a walker ) decays as a power of time , with a @xmath2-dependent exponent , @xmath23 .    among the questions which emerge naturally in such a study
are ( i ) what is the dominant length scale in the problem , as far as the persistent structures are concerned ?
( ii ) what is the nature of the spatial distribution of the persistent sites ? ( iii ) what is the average size of a persistent cluster ?
an important point to make at the outset concerns the different length scales associated with the walkers and with the persistent sites .
the walker density decays as @xmath24 for any @xmath2 , so the mean distance between walkers grows as @xmath25 .
the density of persistent sites decays as @xmath26 , so the mean distance between these sites grows as @xmath27 .
this simple observation immediately suggests that the spatial structure of the persistent sites for @xmath8 , where @xmath28 is the larger length scale , will be very different from when @xmath29 , where @xmath30 is the larger length .
this is precisely what we find : the characteristic length scale , @xmath31 , controlling the distribution of the _ intervals _ between clusters of persistent sites , is given by @xmath32 .
for the ising case studied earlier @xcite , one has @xcite @xmath33 , which is ` close ' to @xmath34 .
we believe the proximity of these two exponents is responsible for the apparently non - universal behavior in the ising system reported in mr
.    derrida et al .
@xcite have obtained an exact expression for @xmath15 for 1d potts models with arbitrary @xmath2 : @xmath35 ^ 2\ .
\label{theta(q)}\ ] ] the value of @xmath2 corresponding to @xmath36 is @xmath37 = 2.70528\ldots$ ] , so @xmath38 for all integer @xmath39 .
note that the probabilistic algorithm for implementing the potts model through the annihilation or coalescence of random walkers ( domain walls ) allows @xmath2 to be a real number , @xmath40 , while an equivalent ising spin representation ( see section iv ) of the potts persistence problem even allows @xmath41 ! by this means means we can explore a range of @xmath15 above and below 1/2 .
our main result is that the scale length controlling the distribution of interval sizes between persistent clusters is given by @xmath10 when @xmath29 , but by @xmath9 when @xmath8 .
we find no evidence for any dependence of the asymptotic scaling distribution on the initial walker density ( other than through nonuniversal amplitudes ) .
this paper is organised in the following manner . in section
ii we outline a scaling phenomenology within which the above questions can be addressed . in section iii we give exact results for the @xmath42 potts model . finally , in section iv , we present extensive numerical simulations for the @xmath2-state potts model which confirm our predictions based on general scaling arguments and the exact large-@xmath2 results .
section v concludes with a discussion and summary of the results .
the basic scaling phenomenology was introduced by manoj and ray , although we will adopt a slightly different notation .
we are interested in the distribution of the sizes of the _ intervals _ between persistent sites .
we define an interval size , @xmath1 , as the number of _ non - persistent _ sites between two consecutive persistent sites .
if the interval size is zero , the sites belong to the same _ cluster _ of persistent sites , and we henceforth consider only intervals of non - zero size . we define @xmath0 to be the number of intervals ( per site ) of size @xmath1 at time @xmath7
a natural dynamical scaling assumption is @xmath43 where @xmath44 is the ` characteristic ' length scale at time @xmath7 .
( a dynamical exponent @xmath45 is conventionally defined via @xmath46 , rather than @xmath47 .
here we are following the notation of mr ) .
the rationale for the scaling form ( [ scaling ] ) is as follows .
the number of non - persistent sites ( per site ) is given by @xmath48 but since the number of persistent sites ( per site ) , @xmath12 , decays to zero as @xmath6 , and @xmath49 , it follows that @xmath50 for @xmath51 . converting the sum ( [ np ] ) to an integral , with lower limit zero ( valid as @xmath52 provided the integral converges )
, we see that the prefactor @xmath53 in ( [ scaling ] ) is precisely what is needed to satisfy the required condition @xmath54 for @xmath51 .
consider next the number , @xmath55 , of persistent clusters per site .
since the number of clusters is equal to the number of intervals , @xmath56 is given by @xmath57 converting the sum to an integral ( with lower limit zero ) , using ( [ scaling ] ) and assuming the integral converges , gives @xmath58 .
thus the mean distance between persistent clusters , i.e.  the mean interval size , increases as @xmath47 .
is this reasonable ?
let us recall that there are two length scales in the system , @xmath59 , the mean distance between walkers , and @xmath60 , the mean distance between persistent sites . to make further progress we make the following two assumptions ,
both of which are confirmed by our numerical studies , by exact results for @xmath61 , and by heuristic arguments to be expounded below : ( i ) the mean cluster size tends to a constant for @xmath51 , and ( ii ) the length scale @xmath62 that controls the interval size distribution is the _ larger _ of @xmath30 and @xmath28 , i.e.   @xmath63 .    from assumption ( i )
we deduce that @xmath64 , i.e.  @xmath65 , which is consistent with assumption ( ii ) provided @xmath8 .
what if @xmath29 ?
then we still require @xmath64 , but now @xmath66 , so the result @xmath58 , derived from ( [ cl ] ) breaks down . going back to ( [ cl ] ) , we infer that this breakdown requires that the conversion of the sum to an integral be invalid  the sum must have its dominant contribution from @xmath1 of order unity , rather than @xmath1 of order @xmath47 .
this in turn requires that the scaling function @xmath67 in ( [ scaling ] ) have a singular small-@xmath17 limit of the form @xcite @xmath68 with @xmath69 .
using this form in ( [ scaling ] ) , with @xmath66 , and inserting the result into ( [ cl ] ) , gives @xmath70 . comparing this with @xmath64 fixes the value of @xmath71 : @xmath72    we first present a heuristic argument for the result @xmath73 . consider first the case @xmath8 .
in this regime , the distance between clusters is much greater that the distance between walkers , i.e.  there are , on average , many ( @xmath74 ) random walkers between each consecutive pair of persistent clusters , as illustrated in figure [ fig0 ] .
however , walkers separated by distances large compared to @xmath10 are essentially uncorrelated , because any correlations are mediated by the random walkers , and correlations between walkers decay on the length scale @xmath10 .
therefore one expects the intervals between clusters , measured on the scale @xmath9 , to be independent random variables , and the interval size distribution to scale with this length .
this argument also suggests that the locations of the clusters as a function of position on the lattice are described by a poisson process , i.e.  the interval size distribution is an exponential , @xmath75 , with @xmath76 , for @xmath8 .
these expectations are borne out by our numerical studies .    for @xmath29 ,
the opposite is true : there are many clusters between each pair of walkers ( see the @xmath77 snapshot in figure [ fig0 ] ) .
the largest length scale is set by the walker spacing , @xmath59 , and the interval size distribution will scale with this length . on smaller scales , there will be residual structure in the cluster distribution left from an earlier epoch when the number of walkers was larger .
indeed , on these smaller scales the persistent sites form a fractal set @xcite .
consider two sites separated by a distance @xmath78 .
the probability , @xmath79 , that _ both _ are persistent has the scaling form , @xmath80 where @xmath81 for @xmath82 , since for @xmath83 the sites are uncorrelated .
the residual structure remaining on scales @xmath84 , however , means that if the first site is persistent , which occurs with probability of order @xmath26 , the probability that the second one is also persistent depends only on @xmath78 .
this implies @xmath85 for @xmath86 .
the number of persistent sites within a distance @xmath87 of a given persistent site can therefore be estimated as @xmath88 , where @xmath89 is the fractal dimension of the persistent set .
clearly this result only makes sense for @xmath29 , since @xmath90 can not be negative .
this suggests @xmath91 for @xmath92 , corresponding to point - like objects , i.e.  isolated finite clusters . for @xmath61 , for example , where @xmath93 , we find that there is typically one cluster of persistent sites , with a fixed mean size , per scale length @xmath94 .
we now present a heuristic argument in support of our claim that the mean cluster size , @xmath95 , approaches a constant as @xmath51 for all values of @xmath15 .
the initial steps of the argument follow the approach of mr .
the total number of clusters per site can be written as @xmath96 where @xmath97 is the fraction of walker sites which have never been visited by a walker , i.e.  the persistence probability for walker sites .
such sites form ` spacers ' between each pair of persistent potts states in a cluster , there being exactly one more persistent potts site than persistent walker site in each cluster .
the mean cluster size is @xmath98^{-1}\ . \label{lav}\ ] ] we expect , on universality grounds , that the exponents describing the decay of @xmath12 and @xmath97 should be the same ( i.e.  @xmath15 ) , but the ( nonuniversal ) amplitudes will be different , i.e.  @xmath99 , while @xmath100 , with @xmath101 . inserting these forms in ( [ lav ] ) gives the limiting value of the mean cluster size as @xmath102 .
this number is nonuniversal and is determined by the initial distribution of cluster sizes .
many of the general features of the dynamics for @xmath104 are exemplified by the @xmath105 limit , for which the walker dynamics reduces to @xmath106 , i.e.  the walkers always aggregate on contact . at @xmath107
all sites are persistent .
the random walkers initially present divide the sites into clusters of persistent sites .
clearly no cluster can increase in size .
consider a cluster of initial size @xmath108 .
we first calculate the probability density , @xmath109 , that the cluster survives to time @xmath7 and has size @xmath110 .
the key point is that , for @xmath106 dynamics , we need only consider the two walkers on either side of the initial cluster , since subsequent coalescence processes do not affect the random walk dynamics of these two walkers .
the cluster and the two walkers can therefore be treated in isolation .    for simplicity
we begin by treating the continuum limit in which the walkers are regarded as continuous - time random walkers on a continuous space .
this should be correct in the limit @xmath111 , an assertion we can subsequently test . let the ends of the initial cluster ( and the two walkers ) be located at @xmath112 and @xmath113 .
each walker obeys an equation of the form @xmath114 , where @xmath115 is a gaussian white noise with mean zero and correlator @xmath116 .
we first write down the probability distribution @xmath117 of the rightmost excursion , @xmath118 , up to time @xmath7 , of the left walker .
an elementary calculation gives @xmath119 for the cluster to survive , we require @xmath120 , so in the limit @xmath121 the exponential factor in ( [ p_r ] ) can be dropped .
the probability distribution , @xmath122 , of the leftmost excursion of the right walker ( from its initial position ) is given by a similar expression .
clearly @xmath123 is the residual size of the cluster at time @xmath7 , where @xmath124 means the cluster has disappeared .
the probability density @xmath109 that the cluster has survived and has size @xmath110 is therefore given , for @xmath121 , by @xmath125 the probability , @xmath126 , that the cluster survives till time @xmath7 is given by @xmath127    immediately we see that the mean interval between surviving clusters grows as @xmath7 ( @xmath128 , since @xmath129 ) , which is much larger than the mean interval between walkers , which grows as @xmath10 .
the mean length of surviving clusters , of given initial size @xmath108 , in the @xmath51 limit is @xmath130 if the initial clusters have a distribution of sizes , a straightforward calculation gives @xmath131 , where @xmath132 indicates an average over this distribution , while the fraction of clusters which survive is @xmath133 .
the above calculations demonstrate that the mean cluster size approaches a constant at late times , and that the scale length for the sizes of the intervals between clusters grows as @xmath7 , and not as the naive scaling length @xmath10 associated with the underlying domain wall coarsening .
in fact we can calculate the interval size distribution , @xmath0 , explicitly for this model .
first recall that the surviving clusters are separated , at late times , by many ( @xmath134 ) walkers .
the fate of a given cluster depends only on the nearest walker on either side .
the motion of these walkers is uncorrelated with that of the walkers bordering other clusters , because the correlation length for walkers grows only as @xmath10 .
therefore we can assume , in the scaling limit @xmath135 , @xmath51 , with @xmath136 fixed , that clusters survive _ independently _ with probability @xmath137 ( where we have specialized to initial clusters of fixed size @xmath108 ) .
the probability distribution of the interval size @xmath1 between neighboring clusters is therefore exponential , with mean @xmath138 , while @xmath139 , which has the scaling form ( [ scaling ] ) with @xmath140 .
this exponential form is in excellent agreement with simulation data presented in section iv .    since the initial clusters used in the @xmath61 simulations are quite small ( 2 or 4 ) , the continuum limit is not expected to be quantitatively correct . to conclude this section , therefore , we quote the asymptotic mean cluster size for clusters of arbitrary initial size @xmath108 .
the details are given in the appendix : the result is @xmath141 , which generalizes the continuum result @xmath142 .
let us recall the zero - temperature coarsening dynamics of the 1d @xmath2-state potts model , starting from a random initial configuration . to maximize the speed of the program , we adopt two - sublattice parallel updating .
the zero - temperature dynamics proceeds as follows : in each time step every spin on one of the sublattices changes its colour to that of one of its two neighbours with equal probability , the two sublattices being updated alternately .
the dynamics of such systems can equivalently be formulated in terms of the motions of the domain walls as a reaction - diffusion model .
the domain walls can be considered as random walkers performing independent random walks .
whenever two walkers meet , they annihilate ( @xmath19 ) with probability @xmath22 or aggregate ( @xmath143 ) , to become a single walker , with probability @xmath21 .
the persistence is the probability that a fixed point in space has not been traversed by any random walker . for @xmath20 ,
the particles only annihilate and hence this is equivalent to the glauber model , whereas in the @xmath61 limit the walkers only aggregate .
this algorithm is , however , restricted to modelling @xmath2-state potts models with @xmath40 , since @xmath144 generates negative probabilities , while the general result ( [ theta(q ) ] ) allows any real @xmath145 . in order to simulate values of @xmath2 in the interval @xmath146 we can map the potts model ( as far as persistence properties are concerned ) on to an ising representation with a fraction @xmath147 of the sites initially pointing up , and initially persistent , while the down spins are nonpersistent . studying the persistence of the minority ( majority ) spins
then corresponds to the case @xmath148 ( @xmath144 ) @xcite .
the domain walls between the up and down spins behave as random walkers , annihilating each other on contact with probability one . as far as the persistence is concerned , the dynamics of the potts and ising models are completely equivalent @xcite . for @xmath148 , results obtained from both types of simulation are entirely consistent .
the numerical simulations are performed on a 1d lattice of size @xmath149 with periodic boundary conditions .
the ising spins , or potts variables , occupy the odd sites of the lattice , while the random walkers are restricted to even - numbered sites . therefore the effective size of the lattice , in terms of ising / potts spins , is @xmath150 .
each walker jumps to an adjacent even - numbered site with equal probability , turning any persistent site it hops over into a nonpersistent site .
the positions of all walkers are updated simultaneously at each monte - carlo step .
the initial configurations are chosen to eliminate the possibility of any ` crossover ' of random walkers when they jump .
this can be done by placing them only on sites whose positions on the lattice are of the form @xmath151 , where @xmath1 is an integer . the walkers then occupy subsets of the sites @xmath151 and @xmath152 alternately .    for the potts simulations ( @xmath148 ) , the lattice is initially completely persistent .
the random walkers are periodically laid down , in clusters of uniform initial size @xmath108 .
the walkers then execute random walks , aggregating or annihilating according to the prescribed probabilities , @xmath21 and @xmath22 respectively .
the simulations are performed for initial cluster sizes of 2 and 4 .
the values of @xmath2 chosen are @xmath153 and 4 , for which the corresponding values of @xmath15 are @xmath154 , @xmath155 , and @xmath156 .
for @xmath157 the ising representation is employed , in which persistent sites ( up spins ) are periodically laid down on the lattice , with a fraction @xmath147 of all sites persistent , and the domain walls execute random walks annihilating according to @xmath19 .
the value @xmath77 was chosen , for which @xmath158 .
the simulations were performed for initial cluster sizes 4 and 8 .
the ising simulations were also run for @xmath159 and @xmath160 .
the results for scaling functions are consistent with the potts data , but the statistics are poorer due to the smaller number of clusters generated .
for this reason , the data presented in figs .
3 - 9 were generated from the potts runs .
all the simulations ( potts and ising ) were run for @xmath161 monte - carlo steps , and the results averaged over 30 independent runs .
the simulations investigate the limiting values of the average sizes of the persistent clusters , and the distribution , @xmath0 , of intervals between persistent clusters with emphasis on its dynamic scaling form .
while the first of these lends itself readily to direct measurement , the dynamic scaling form for @xmath0 manifests itself most clearly in the data through a study of the ( complement of the ) cumulative distribution .
we define @xmath162 inserting the scaling form ( [ scaling ] ) and converting the sum to an integral in the scaling limit @xmath135 , @xmath51 with @xmath163 fixed but arbitrary , gives @xmath164    we now discuss the numerical results for the @xmath61 model . in figure [ fig1 ]
we plot the mean cluster size against time for initial cluster sizes of 2 and 4 .
the numerical results are summarized in table 1 .
the numerics clearly support the exact results given in section iii [ and derived in the appendix as equation([l_av ] ) ] , namely @xmath165 where @xmath108 is the initial cluster size . in order to investigate the dynamic scaling form ( [ newscaling ] ) for @xmath61 , we plot , in figure [ fig2 ] , @xmath166 against @xmath136 for @xmath167 500 , 1000 , 2000 , 4000 , 8000 , for initial cluster size 2 , where we have used @xmath168 for @xmath61 .
excellent data collapse occurs in agreement with the dynamic scaling form ( [ newscaling ] ) .    [ cols="^,^,^",options="header " , ]     .
average cluster size at late times , @xmath169 , for various values of @xmath2 and various initial cluster sizes @xmath108 , where p or i indicate whether potts or ising representations were employed . in all cases
the data were averaged over times between 1000 and 9000 monte carlo steps .
in section iii we argued that the scaling function @xmath67 in ( [ scaling ] ) should be a simple exponential , which implies that @xmath170 in ( [ newscaling ] ) is also a simple exponential . to test this prediction , we plot @xmath171 against @xmath136 in figure [ fig3 ] .
the data lie on the expected straight line except at the latest times , where the deviations from the line are scatter due to statistical noise .
very similar results are obtained for initial cluster size 4 , confirming that the scaling function @xmath67 is independent of the initial walker density .    in figure
[ fig4 ] the mean cluster size is plotted against time for various finite values of @xmath2 and initial cluster size @xmath108 .
the numerical results indicate that the average cluster size tends to a constant ( see table i ) in all cases , as expected on the basis of the heuristic argument in section ii . for @xmath172 , @xmath173 and the dominant length scale
is still determined by the mean distance between persistent clusters .
we therefore anticipate a scaling form of the kind given by ( [ newscaling ] ) with @xmath174 .
the scaling plot presented in figure [ fig5 ] shows good data collapse with this value of @xmath45 , while figure [ fig6 ] indicates that @xmath170 , and therefore @xmath67 , is again a simple exponential .
the data collapse is not as good as for @xmath61 , which we attribute to a poorer separation of the length scales @xmath28 and @xmath30 at the time scales available .
this problem becomes more acute for @xmath175 ( see the discussion below ) .    for @xmath175 , @xmath176
, so we plot ( figure [ fig7 ] ) @xmath177 against @xmath163 with @xmath178 .
this time the scaling collapse is definitely not as good as for @xmath61 and @xmath172 , reflecting an apparent deviation from the dynamic scaling form ( [ newscaling ] ) .
we believe that this deviation reflects the proximity of the length scales @xmath179 and @xmath59 , which are not sufficiently well separated on the timescales achievable in the simulation , and would disappear for asymptotically large times . for example ,
@xmath180 at @xmath181 , while the corresponding ratio is about 14 for @xmath172 and 90 for @xmath61 .
a similar effect may explain the apparent deviations from universality ( e.g.  an apparent dependence of @xmath45 on the initial walker density ) in the @xmath20 simulations of mr .
in this case @xmath33 , and @xmath182 at @xmath183 , the largest time reached by the mr simulations . so in this case
also there is not a strong separation of length scales at the times achieved in the simulations ( a point recently emphasized by manoj and ray @xcite ) .
figure [ fig8 ] gives the log - linear plot for @xmath175 .
again , the data are consistent with an exponential scaling function , except at small scaling variable where there is a perceptible deviation from a straight line . in this region , however , the data also do not scale at all well either ( see inset in figure [ fig8 ] ) , which may be another manifestation of the poor separation of length scales : the scaling limit requires @xmath135 , @xmath51 with @xmath163 fixed . in particular the condition @xmath184 , which should be satisfied for good scaling ,
is violated at the small values of @xmath163 where the upturn in the data occurs .
we will return to this point in the discussion .    to investigate the spatial distribution of the sites in the regime @xmath29 we study @xmath77 , for which @xmath185 . in figure [ fig4 ]
we show that the mean cluster size is a constant for initial cluster sizes of 4 and 8 . in the regime
@xmath29 the density of the persistent sites decays more rapidly than that of the walkers .
the dominant length scale is therefore no longer given by the persistence length , @xmath60 , but is given instead by the mean separation , @xmath59 , of the the walkers .
hence we expect asymptotic dynamic scaling of the form ( [ scaling ] ) and ( [ newscaling ] ) with @xmath66 . in figure [ fig9 ]
we plot @xmath186 against @xmath187 for an initial cluster size @xmath188 .
the numerical results give excellent data collapse , supporting the proposed dynamic scaling form .
apart from a change of scale , very similar results are obtained for @xmath189 , supporting the universality of the scaling function , i.e.  the independence of @xmath170 from the initial walker density .
we choose @xmath77 rather than @xmath20 in order that the length scales @xmath28 and @xmath30 be reasonably well separated at late times : @xmath190 at @xmath181 .
figure [ fig9 ] reveals a pronounced upturn at small @xmath1 , which can be seen more clearly in the log - linear plot of figure [ fig10 ] and its inset , suggesting a divergence for @xmath191 .
this is , in fact , expected from the analysis of section ii , which predicts , for @xmath29 , @xmath192 for @xmath86 , with @xmath193 , i.e.  @xmath71 is in the range @xmath69 .
the scaling function @xmath170 is given by @xmath194 for @xmath86 , i.e.  @xmath195 , where @xmath196 is the fractal dimension of the persistent set on scales smaller than @xmath30 [ equation ( [ df ] ) ] . according to this prediction , the product @xmath197 should approach a constant at small @xmath17 .
this product is shown in figure [ fig11 ] .
the small-@xmath1 divergence has clearly been removed , the function approaching a constant at small @xmath17 as predicted ( the erratic behaviour at very small @xmath17 is due to statistical noise ) .    for @xmath77
all of the numerical simulations have been performed using the ising representation , with periodic initial conditions in which clusters of @xmath108 ` up spins ' are placed at uniform intervals in a ` down spin ' background , and occupy a fraction @xmath147 of the sites . for @xmath148 , the potts representation has mainly been used , but we expect all universal properties , such as exponents and scaling functions , to be the same for the two representations . this expectation is confirmed by the results for @xmath5 and the scaling function @xmath170 . in order to obtain a _ precise _ correspondence between the ising and the potts representations , however , it is necessary to run the ising simulations with random initial conditions , rather than periodic initial conditions which we have mostly employed , and to scale the axes appropriately with @xmath2 to account for the different interval sizes in the two cases .
this is because the reaction - diffusion representation of the potts model , in which walkers annihilate or coalesce with probabilities which depend on @xmath2 , is only an exact representation if the potts states occur in a completely random sequence .
this is not true for the periodic initial conditions employed in the ising representation , whose correlations spoil the exact correspondence between ising and potts simulations . our expectation ,
then , is that the ising representation with random and periodic initial conditions should give the same _ scaling functions _ as the reaction - diffusion implementation of the potts model , but non - universal amplitudes will be different for the periodic initial conditions .    in order to test the exact equivalence between random ising and potts simulations we plot , in figure [ fig12 ] , the results for the @xmath172 model using the potts ( reaction - diffusion ) method .
the lattice is initially all persistent and the initial cluster size is @xmath198 .
this corresponds to a sequence of pairs of potts states in which the members of a pair are in same state , but there are no correlations between pairs ( except that neighboring pairs must be in different states ) .
the random initial conditions for the ising simulation are generated from the potts initial state by setting all pairs in potts state ` 1 ' ( say ) to be up spins and all other pairs to be down spins . in the potts simulations we keep track of all persistent sites , while in the ising simulation we track the persistence of the spins initially ` up ' .
the interval sizes ( between persistent clusters ) are then naturally larger by a factor @xmath2 in the ising simulations relative to the potts simulations .
figure [ fig12 ] shows both scaling functions in log - linear form , for @xmath199 , with this factor of @xmath2 scaled out ( i.e. @xmath200 is plotted against @xmath201 for these data ) .
the data overlap almost perfectly .
the ising data for a periodic initial condition are consistent with the same scaling function ( exponential ) , but with a different amplitude .
in this paper we have investigated the nature of persistent structures in the coarsening dynamics of 1d potts models .
a central concept has been the existence of two characteristic length scales , the mean separation , @xmath59 , of domain walls ( or ` walkers ' ) and the ` persistence length ' , @xmath60 , which measures the mean separation of persistent sites or clusters .
the focus of our attention has been the distribution function , @xmath0 , for the number of intervals ( between clusters of persistent sites ) of length @xmath1 .
this distribution has the scaling form ( [ scaling ] ) , with characteristic length scale , @xmath47 , which is the _ larger _ of @xmath30 and @xmath28 , i.e.  @xmath4 . within the 1d potts model
both regimes , @xmath29 and @xmath8 , can be accessed by varying the number of potts states , @xmath2 . from the general result , ( [ theta(q ) ] ) , we see that @xmath5 is a monotonically increasing function with , @xmath202 marking the boundary between the two regimes .
the regime @xmath8 is conceptually simpler .
the @xmath61 limit can been solved exactly ( section iii ) , and results for the long - time limit of the mean cluster size obtained .
the fact that the locations of the surviving clusters are statistically independent in the scaling regime leads to the result that the scaling function @xmath67 in ( [ scaling ] ) is a simple exponential .
we have argued that the same result should hold for all @xmath8 , since in this regime there are , on average , many walkers between each pair of persistent clusters .
the _ relevant _ walkers for this argument ( those which are turning persistent sites in neighboring clusters nonpersistent ) are uncorrelated since their separations , of order @xmath28 , are much larger than the typical separation , @xmath30 , of neighboring walkers . on this basis
we expect the asymptotic scaling function @xmath67 to be exponential for all @xmath8 ( i.e.  @xmath203 ) .
the data for @xmath172 ( figure [ fig6 ] ) and @xmath175 ( figure [ fig8 ] ) are consistent with an exponential form for @xmath204 , although the scaling is not perfect and there is a small upturn in the scaling function at small scaling variable for @xmath175 .
we attribute these features to an imperfect separation of length scales on the time scales achieved in the simulations , and conjecture that the true scaling function is exponential for all @xmath104 .
the case @xmath29 is more tricky . in this case
the persistent clusters outnumber the walkers .
the scaling function @xmath67 is clearly _ not _ a simple exponential , though it seems ( figure [ fig10 ] ) to have an exponential tail .
there is a small - argument singularity of the form @xmath192 , with @xmath193 .
this is related to the fractal dimension , @xmath196 , of the persistent sites : @xmath205 .
note that the borderline , @xmath206 , between the two regimes occurs at @xmath207 .
the existence of a small-@xmath17 singularity for @xmath29 raises the possibility of an alternative scenario for @xmath104 , in which the @xmath208 singularity , with @xmath209 , persists for @xmath8 , where it becomes an integrable singularity .
the small-@xmath17 singularity in @xmath170 would then take the form of a cusp : @xmath210 .
we have not been able to rule out this scenario from the data , but think it unlikely for the reasons given elsewhere in this paper .
further insight could be obtained if it were possible to perform an expansion around the soluble @xmath61 limit to first order in @xmath147 , but we leave this as a challenge for the future .
we conclude by discussing briefly the possibility of the existence of these two qualitatively different regimes in spatial dimension @xmath211 .
first we generalize the result ( [ df ] ) , relating @xmath90 and @xmath15 , to any dimension . starting from ( [ both ] ) , the result @xmath212 follows generally , and the number of persistent sites within a distance @xmath87 of a given site is estimated as @xmath213 , where @xmath214 .
generalizing still further , if the coarsening exponent is @xmath215 , rather than @xmath34 , the result is @xmath216 ( the @xmath217 version of this result is given in @xcite , based on the same reasoning ) .
clearly this result requires @xmath218 , since @xmath90 can not be negative .
if this inequality is violated , as in the 1d potts model with @xmath104 , the persistent sites no longer have a fractal structure but become point - like objects , with mean density @xmath219 , where @xmath220 is the coarsening length scale .
cases where @xmath221 are easy to find , e.g.  in the coarsening of the 2d ising model @xcite , the 2d diffusion equation @xcite , or the time - dependent ginzburg - landau equation in 2d . these all exhibit fractal persistent structures with the expected fractal dimension @xcite
. it would be interesting to look for examples , in addition to the 1d potts model , where one can have @xmath222 .
this work was supported by epsrc ( uk ) .    for a recent review
see s. n. majumdar , curr .
india * 77 * , 370 ( 1999 ) .
b. derrida , a. j. bray , and c. godrche , j. phys .  a * 27 * , l357 ( 1994 ) .
b. derrida , v. hakim , and v. pasquier , phys .
lett .   * 75 * , 751 ( 1995 ) ; j. stat .
85 , 763 ( 1996 ) .
d. stauffer , j. phys .
a * 27 * , 5029 ( 1994 ) .
s. n. majumdar , c. sire , a. j. bray , and s. j. cornell , phys .
lett .  * 77 * , 2867 ( 1996 ) ; b. derrida , v. hakim , and r. zeitak , ibid .
* 77 * , 2871 ( 1996 ) .
j. krug , h. kallabis , s. n. majumdar , s. j. cornell , a. j. bray , and c. sire , phys .
e * 56 * , 2702 ( 1997 ) .
a. j. bray , b. derrida , and c. godrche , europhys .
* 27 * , 175 ( 1994 ) ; p. l. krapivsky , e. ben - naim , and s. redner , phys .
e * 50 * , 2474 ( 1994 ) ; s. n. majumdar and c. sire , phys .
* 77 * , 1420 ( 1996 ) .
j. cardy , j. phys .
a * 28 * , l19 ( 1995 ) ; e. ben - naim , phys .
e * 53 * , 1566 ( 1996 ) ; m. howard , j. phys .
a * 29 * , 3437 ( 1996 ) ; s. j. cornell and a. j. bray , phys .
e * 54 * , 1153 ( 1996 ) .
m. marcos - martin , d. beysens , j - p .
bouchaud , c. godrche , and i. yekutieli , physica ( amsterdam ) 214d , 396 ( 1995 ) .
w. y. tam , r. zeitak , k. y. szeto , and j. stavans , phys .
lett .  * 78 * , 1588 ( 1997 ) .
b. yurke , a. n. pargellis , s. n. majumdar , and c. sire , phys .  rev .
e * 56 * , r40 ( 1997 ) .
s. n. majumdar , a. j. bray , s. j. cornell , and c. sire , phys .
lett .  * 77 * , 3704 ( 1996 ) ; k. oerding , s. j. cornell , and a. j. bray , phys .
e * 56 * , r25 ( 1997 ) .
g. manoj and p. ray , preprint , cond - mat/9901130 ( 1999 ) ; preprint , cond - mat/0003203 ( 2000 ) .
d. h. zanette , phys .
e * 55 * , 2462 ( 1997 ) .
g. manoj and p. ray , j. phys .
a * 33 * , l109 ( 2000 ) .
s. jain and h. flynn , preprint , cond - mat/0004148 ( 2000 ) .
e. ben - naim , l. frachebourg , and p. l. krapivsky , phys .
e * 53 * , 3078 ( 1996 ) .
s. j. cornell and a. j. bray , unpublished .
see http://theory.ph.man.ac.uk/@xmath223bray/pictures.html for some illustrations of the persistent structures which arise .
in this appendix we compute the probability that a cluster of initial size @xmath108 survives to time @xmath7 , and the mean size of surviving clusters , for the process @xmath106 .    first note that initially there is a random walker at each end of the domain . at each time
step the walkers move independently left or right with probability @xmath34 , so we can treat each walker independently .
consider , therefore , a single walker moving at discrete time steps on a discrete 1d lattice , starting , at time @xmath107 , at position @xmath78 .
let the ` origin ' be the point @xmath224 , and let @xmath225 be the probability that the walker has not yet reached the origin at time @xmath7 .
clearly , @xmath226,\ \ \ r \ge 2\ .
\label{recurrence}\end{aligned}\ ] ] we are interested only in the limit @xmath51 . in this limit , we know from standard random walk theory that every @xmath225 decays like @xmath24 , with an @xmath78-dependent amplitude . to leading order in @xmath24 , therefore , the @xmath7-dependence drops out of equations ( [ recurrence ] ) , which then become equations for the amplitudes . by inspection ,
the solution in this regime is @xmath227    now consider a walker starting immediately to the right of a cluster of @xmath108 persistent sites .
the probability that after @xmath7 steps the walker has jumped over exactly @xmath78 of these ( making them nonpersistent ) is @xmath228 , where the final result follows from ( [ leading ] ) , to leading order for large @xmath7 .
the same result holds for a walker starting immediately to the left of the cluster . the probability that @xmath110 sites remain persistent after time @xmath7 is @xmath229 times
the number of ways of partitioning the cluster of length @xmath108 into 3 sections , with the central section of length @xmath110 ( and zero lengths are allowed for the outer sections ) .
this number is @xmath230 .
so the probability of the cluster surviving and having length @xmath110 is @xmath231 a generalization of ( [ pl0 ] ) to the discrete system .


there have been a number of attempts to study hypernuclei containing not only nucleons but also hyperons with strangeness @xcite .
it has been of high interest to investigate how the presence of strangeness broadens and/or modifies the knowledge achieved in the field of ordinary nuclear physics .
it is a fundamental problem to understand properties of a quantum many - body system in terms of basic interactions between constituent particles . in hypernuclear physics a unified treatment of all the baryon - baryon interactions , especially nucleon - nucleon ( @xmath4 ) and hyperon - nucleon ( @xmath5 ) interactions ,
will be necessary if one wish to understand general properties of hypernuclei . however , our knowledge on the @xmath5 interaction is quite inadequate in contrast to that of the @xmath6 interaction due to experimental difficulties .
in such a situation many - body theoretical studies of hypernuclei could provide alternative information on the @xmath5 interaction .
it has been discussed that physics of hypernuclei is different from conventional nuclear physics in some aspects .
first , the channel couplings such as the @xmath7-@xmath8 coupling play a significant role in the structure of hypernuclei in contrast to the very limited role of the @xmath9-@xmath10 coupling appeared in nuclear many - body systems @xcite .
second , the anti - symmetric spin - orbit force , which is forbidden in ordinary nuclei , emerges in hypernuclei in addition to the symmetric spin - orbit force .
the spin - orbit splitting of @xmath11 single - particle levels in @xmath11 hypernuclei was considered experimentally to be much smaller than their nucleonic counterpart @xcite . however , some recent studies as in ref .
@xcite suggest that the splitting would not be so small as discussed earlier
. it will be highly desirable for the spin - orbit splitting to be established in experiment .
the magnitude of the spin - orbit splitting will put further constraints on the @xmath12 interaction .
there have been several possibilities of solving the structure of hypernuclei .
one possibility is to treat few - body or light hypernuclei by solving the faddeev equation @xcite or by applying the cluster model @xcite .
another possibility is microscopic studies of hypernuclei with a larger mass number .
while a few lower partial waves play a decisive role in few - body or light hypernuclei , higher partial waves could give some contributions in heavier hypernuclei .
these microscopic studies of hypernuclei could make us understand interesting physics of their structure itself , and also help us in determining the @xmath5 interaction .
microscopic structure calculation of heavier hypernuclei is usually performed by introducing a @xmath12 effective interaction . for this purpose
the @xmath13-matrix theory has been applied in two ways : in some works an effective @xmath5 interaction has been derived by simulating a nuclear matter @xmath13-matrix . in this approach the non - locality and energy dependence of the @xmath13-matrix
are neglected , and the effective @xmath5 interaction is represented in the two- or three - range gaussian form with the adjustable fermi momentum @xmath14 .
this type of the effective interaction has been referred to as the yng interaction @xcite .
the yng interaction has been successfully used in many of hypernuclear structure calculations @xcite .
it seems to us that the yng interaction has made an important contribution to the considerable progress in overall description of hypernuclei .
the other approach , in the framework of the @xmath13-matrix theory , is to calculate the @xmath13-matrix in finite hypernuclei by treating rigorously the pauli exclusion operator and making a self - consistent calculation of single - particle potential of hyperons @xcite .
in contrast to studies of using the yng interaction this approach possibly gives a better description of the state dependence of the effective @xmath5 interaction .
it has been established that the @xmath13-matrix can be used as a basic ingredient in constructing a microscopic many - body theory starting with bare interactions .    in spite of a marked success in the applications of the @xmath13-matrix theory
, we may say from a formal point of view that the @xmath13-matrix is energy - dependent and non - hermitian , and does not have the property of decoupling between a model space of low - lying two - particle states and its complement .
the @xmath13-matrix is not considered to be an effective interaction in a formal sense of the effective interaction theory . in order to derive an effective interaction
we should add some higher - order corrections such as folded diagrams @xcite .
it would be desirable if we could have a theory describing many - body systems in terms of an energy - independent and hermitian effective interaction with the property of decoupling .
two of the present authors , k.  s. and r.  o. , proposed a many - body theory , the unitary - model - operator approach ( umoa ) that was formulated on the basis of such an effective interaction .
the umoa was applied to finite nuclei , @xmath15o and @xmath16ca @xcite , and some improvements have been attained .
the method may be viewed as an alternative way of deriving an effective interaction in finite many - body systems .
a unified study of hypernuclei and ordinary nuclei would be helpful in understanding many - body systems of baryons in a microscopic way .
therefore , we extended the formulation of the umoa , in the previous work @xcite , to the description of hypernuclei and made a calculation of the properties in @xmath17o with a realistic @xmath12 interaction .
the purposes of this work are to present a general formulation of the extended umoa and apply it to @xmath17o by employing some realistic @xmath12 interactions .
this paper is organized as follows : in section 2 we present a general formulation of the umoa for deriving the @xmath18 effective interaction . in section 3 the approximation procedure
is given for making actual calculations .
in section 4 we apply the umoa to @xmath17o by employing the nijmegen soft - core ( nsc ) @xcite , jlich model -  ( j ) and model - b ( jb ) @xcite @xmath12 potentials . in section 5
we make some concluding remarks .
we present a formulation of the umoa for applying to a calculation of the properties in @xmath11 hypernuclei , especially , closed - shell nucleus plus one @xmath11 systems .
a main purpose is to give a method of calculating an effective interaction between @xmath11 and @xmath19 . in the derivation of the @xmath8 effective interaction
we should note that , because of the small mass difference of about 77mev between @xmath11 and @xmath20 , the coupling of the @xmath8 and @xmath7 channels plays an important role .
we first consider a two - body hamiltonian of a @xmath7-@xmath8 coupled system given by @xmath21 where @xmath22 @xmath23 and @xmath24 the terms @xmath25 and @xmath26 for @xmath27 and @xmath19 are the kinetic and single - particle potential energies of @xmath28 and @xmath19 , respectively , and @xmath29 is the bare two - body @xmath12 interaction including the @xmath8 and @xmath7 channels .
the terms @xmath30 and @xmath31 are the rest masses of @xmath32 and @xmath33 , respectively .
the hamiltonian @xmath34 describes a subsystem of interacting two particles @xmath8 or @xmath7 in the one - body potentials @xmath35 , @xmath36 and @xmath37 .    in order to derive a @xmath8 effective interaction
we introduce a unitarily transformed @xmath5 interaction given by @xmath38 where @xmath39 is the correlation operator defined in a space of two - particle states of @xmath8 and @xmath7 .
the @xmath39 is an anti - hermitian operator satisfying @xmath40 in the formulation of the umoa the correlation operator @xmath39 is determined by the equation of decoupling for @xmath41 between a certain model space consisting of low - momentum two - particle states and its complement .
this determination of the correlation operator @xmath39 is one of the characteristics of the present approach .    for the determination of a @xmath8 effective interaction
, we introduce a model space consisting of low - momentum @xmath8 states .
let @xmath42 be the projection operator onto the model space .
the complement of the @xmath43 space includes a space of high - momentum @xmath8 states referred to as the @xmath44 space and a space of all the @xmath7 states referred to as the @xmath45 space .
the condition of decoupling for the transformed interaction @xmath41 is given by @xmath46 hereafter , we refer to the above equation as the decoupling equation .
we assume that the one - body part @xmath47 in eq .
( [ eq : tvyn ] ) is decoupled as @xmath48 on the above assumption the decoupling equation ( [ eq : dec ] ) becomes @xmath49 in general , the decoupling equation ( [ eq : decr ] ) does not determine @xmath39 uniquely .
the usual restrictive conditions are @xmath50 these conditions are called the minimal effect requirements @xcite , which mean that the transformation exp@xmath51 does not induce unnecessary transformation within each of the @xmath52 and @xmath53 spaces .    in the general theory of constructing a hermitian effective interaction , the solution for @xmath54 in eq .
( [ eq : decr ] ) has been known @xcite and is given by @xmath55 with latexmath:[\ ] ] the @xmath59 in eq .
( [ eq : oyn ] ) is the biorthogonal state of @xmath60 satisfying latexmath:[\[\begin{aligned } \label{eq : bio } \langle \overline{\psi_{{\sst { \it \lambda n}}k } }    where @xmath60 is the @xmath52-space component of @xmath57 defined by latexmath:[\[\begin{aligned } \label{eq : psc }     that the solution @xmath39 in eq .
( [ eq : syn ] ) is determined dependently on the choice of a set of @xmath63 eigenstates @xmath64 as given in eq .
( [ eq : oyn ] ) .
we choose a set of @xmath65 so that they have the largest @xmath52-space overlaps among all the eigenstates in eq .
( [ eq : eige ] ) .
this is a usual choice in the derivation of the effective interaction .    with the solution @xmath39 in eq .
( [ eq : syn ] ) , the @xmath8 effective interaction in the @xmath52 space is given by @xmath66 it has been known that this @xmath8 effective interaction can be written explicitly @xcite as @xmath67 where @xmath68 ( @xmath69 ) and @xmath70 ( @xmath71 ) are defined through the eigenvalue equation @xmath72 and @xmath73 is given by @xmath74 the @xmath75 is a @xmath42-space operator that agrees with the effective interaction of non - hermitian type .
equation ( [ eq : evln ] ) is a formula for converting the non - hermitian effective interaction to the hermitian one .
the term @xmath76 in eq .
( [ eq : shyn ] ) is introduced as a self - consistent potential of @xmath32 calculated with the @xmath8 effective interaction @xmath77 as @xmath78 ( @xmath79 ) and @xmath80 are the single - particle states of @xmath32 and @xmath19 , respectively .    from eqs .
( [ eq : evln])-([eq : nhevln ] ) we see that the @xmath8 effective interaction @xmath81 can be calculated once the operator @xmath82 is given .
the operator @xmath82 is determined completely , as given in eq .
( [ eq : oyn ] ) , by solving the eigenvalue equation ( [ eq : eige ] ) and selecting @xmath63 eigenstates @xmath65 with the largest @xmath52-space overlaps .
it is noted here that eqs .
( [ eq : syn])-([eq : ul ] ) give a set of equations that determines the @xmath8 effective interaction @xmath77 self - consistently with the single - particle potential @xmath76 .
here we discuss a calculation procedure for obtaining the @xmath8 effective interaction acting in a small model space of low - lying @xmath8 states .
we use the harmonic oscillator ( h.o . ) wave functions as basis states .
it is desirable that the single - particle potentials of @xmath83 , @xmath33 and @xmath19 in intermediate states are taken into account in a self - consistent way .
however , as our main concern is to derive a @xmath8 effective interaction , the self - consistent treatment of @xmath35 would have predominant importance in the present calculation .
for this reason we make a self - consistent calculation only for @xmath35 , and as for @xmath37 we employ a fixed potential calculated previously for closed shell nuclei in the umoa @xcite .
the effect of the single - particle potential @xmath36 on the @xmath8 effective interaction has not necessarily been clarified , although some calculations with @xmath36 have been made for the @xmath13-matrices in hyperonic nuclear matter @xcite . in the actual calculations made so far the single - particle potential
@xmath36 has been considered to be rather shallow @xcite .
although we have not known clearly how much the effect of @xmath36 on the @xmath8 effective interaction , we here assume @xmath84 .    in solving the self - consistent equation for @xmath35 with @xmath81 we should be careful in the following two respects : ( a ) in the calculation of energies of intermediate states in propagators we should use @xmath35 that is given for @xmath83 states including high - momentum states . in order to do this
we need a @xmath8 effective interaction acting in a large model space , because in general the @xmath44 space contains high - momentum @xmath8 states and the @xmath85-space effective interaction @xmath86 is not always well behaved .
we thus make a calculation by introducing two model spaces , namely , large and small model spaces .
the @xmath8 effective interaction in each of the model spaces is calculated by following a two - step procedure . in the first - step procedure
, we define a large model space and calculate the @xmath8 effective interaction and single - particle potential of @xmath83 in high - lying states
. the @xmath8 effective interaction acting in this large model space will not be suitable for the effective interaction in a usual sense . by introducing a small model space and using the @xmath8 effective interaction determined in the first - step calculation
, we proceed to the second - step calculation where we calculate a @xmath8 effective interaction acting in this small model space of low - lying @xmath87 states .
( b ) it is important in actual calculations how to determine the @xmath88 and @xmath85 spaces because an effective interaction is determined dependently on the choice of the @xmath88 and @xmath85 spaces .
if a state in the @xmath88 space mixes strongly with @xmath85-space states in the eigenstate in eq .
( [ eq : eige ] ) , the matrix element of @xmath89 or @xmath54 defined in sec .
2 becomes large and as a result higher - order many - body correlations @xcite give rise to non - negligible contributions .
in general , the strong mixing takes place when some of the @xmath88- and @xmath85-space states are quasi - degenerate in energy . for this reason
, we should choose the @xmath88 and @xmath85 spaces which are well separated in energy .    with due regard to points ( a ) and ( b ) we perform a calculation of the effective interaction through the two - step procedure .
in the following subsections we present the two - step method in detail .
let us write a two - particle state of @xmath83 and @xmath19 as the product of the h.o .
single - particle wave functions in the usual notation as @xmath90 we define a model space denoted by @xmath91 and its complement @xmath92 with a boundary number @xmath93 as @xmath94 we also define a space of all the @xmath7 states denoted by @xmath95 as @xmath96 note that the model space @xmath91 does not contain the @xmath7 states .
if we make a graph of two - particle states , the @xmath97 , @xmath92 and @xmath98 spaces are represented in fig .  1 .
the @xmath93 means a number such that @xmath99 gives the maximum kinetic energy of the @xmath8 two - particle states in the @xmath97 space .
the @xmath100 denotes the unit of the h.o .
energy which we use commonly for @xmath32 , @xmath33 and @xmath19 .
as we require to introduce a sufficiently large model space in the first - step calculation , we should choose @xmath93 as large as possible .
another requirement is that the states in the @xmath91 space should be separated in energy from the states in the @xmath101 space .
for this reason , we adopt a model space @xmath97 of a triangular shape .
this choice of the @xmath97 space guarantees the separation in energy .
[ fig:1 ] = 13 cm    in fig .
1 the spaces @xmath102 for @xmath103 and @xmath33 , specified by the numbers @xmath104 and @xmath105 , stand for the excluded spaces due to the pauli principle for nucleons . the number @xmath106 is defined as @xmath107 with the quantum numbers @xmath108 of the highest occupied orbit in the closed - shell core nucleus .
the @xmath109 means a number such that @xmath110 gives the maximum kinetic energy of the @xmath12 two - particle states in the @xmath111 and @xmath112 spaces .
we then write the @xmath111 and @xmath112 spaces as @xmath113 and @xmath114 we also define the sum of the excluded spaces as @xmath115 furthermore , by using the projection operators , @xmath101 , @xmath116 , @xmath95 and @xmath117 , we introduce projection operators onto the spaces of allowed @xmath12 states as @xmath118 @xmath119 and we also define @xmath120 considering the pauli - exclusion principle for nucleons the terms @xmath121 and @xmath122 should be removed in treating the problem of two - body correlations .
however , if we remove all the matrix elements of these terms , this leads to overmuch counting of the pauli principle .
the matrix elements which are diagonal in nucleon states should be restored even though nucleons are in occupied states .
therefore , instead of @xmath123 , we should use a @xmath12 interaction @xmath124 defined as @xmath125 where @xmath126 denotes the set of the h.o .
quantum numbers @xmath127 of an occupied state of a nucleon , and the operator @xmath128 is defined as @xmath129    the two - body @xmath12 equation to be solved in the first - step calculation is now written as @xmath130 where we define @xmath131 as @xmath132 the @xmath133 is the single - particle potential of @xmath32 to be calculated with the @xmath8 effective interaction self - consistently .
the number @xmath134 stands for a set of quantum numbers to specify a two - body @xmath12 eigenstate .
the unitarily transformed @xmath12 interaction in the first - step calculation is given by @xmath135 the decoupling equation for determining the correlation operator @xmath136 becomes @xmath137 the @xmath8 effective interaction @xmath138 in the @xmath97 space is then given by @xmath139 in the first - step procedure we solve the decoupling equation ( [ eq : dece1 ] ) by introducing some assumptions . as we take a sufficiently large number as @xmath140
, we may assume that the one - body potential of @xmath83 for a state in the @xmath101 space is negligible in comparison with the sum of the energies , @xmath141 .
therefore , @xmath142 takes nonzero values only for states in the @xmath91 space , and we have @xmath143 the definitions of @xmath97 and @xmath92 in eq .
( [ eq : dtps ] ) allow us to divide a space of @xmath8 states written in relative and center - of - mass ( c.m . )
coordinates into the same subspaces @xmath97 and @xmath92 as @xmath144 where @xmath145 and @xmath146 are the h.o .
quantum numbers of relative and c.m .
motions , respectively .
the @xmath147 and @xmath148 denote , respectively , spin and angular momentum defined through @xmath149 .
the isospin @xmath150 is omitted because @xmath150 takes a constant @xmath151 in the @xmath8 system . a @xmath7 state in the @xmath152 space in fig .  1 is also written by the product of the h.o . relative and c.m .
states .
the @xmath5 interaction @xmath153 in eq .
( [ eq : vbyn ] ) is not diagonal in the c.m .
quantum numbers @xmath154 because the projection operators @xmath97 , @xmath155 and @xmath156 are not diagonal in @xmath154 .
thus , the @xmath8 effective interaction @xmath157 is not diagonal in @xmath154 .
it has been shown , however , in the calculation of ordinary nuclei that @xmath154 non - diagonal elements of the @xmath13-matrix are generally small @xcite .
therefore , we make an assumption that the @xmath154 non - diagonal matrix elements of @xmath157 can be neglected .
the two - body eigenvalue equation ( [ eq : eige1 ] ) is solved for the relative motion with a fixed set of @xmath158 written by @xmath159 where we rewrite @xmath160 as @xmath161 the @xmath162 and @xmath163 are the kinetic energies of the @xmath5 relative and c.m .
motions , respectively .
the state @xmath164 expresses a superposition of the @xmath8 and @xmath7 h.o .
states and gives a general form of the eigenstate when we consider the @xmath7-@xmath8 coupling and the @xmath12 interaction with the tensor and antisymmetric spin - orbit couplings .
the @xmath165 is an additional quantum number specifying a two - body @xmath5 eigenstate .    in general
, the one - body part @xmath166 does not satisfies the assumption of decoupling in eq .
( [ eq : assume ] ) .
however , it has been shown in the calculation in nuclei that if we take a larger number as @xmath167 , the effect of violating the assumption of decoupling becomes considerably smaller @xcite .
we have confirmed that this non - decoupling effect becomes negligible in calculating nuclear effective interactions .
we here assume this to be the case also in the determination of the @xmath8 effective interaction .
furthermore , we introduce another approximation that the sum of @xmath168 and @xmath169 be diagonal in the c.m .
quantum numbers @xmath170 and @xmath171 .
this assumption has also been made sure to be acceptable in the previous works @xcite .    in order to solve the eigenvalue equation ( [ eq : eige1r ] ) we need to calculate the matrix elements of @xmath172 in the relative and c.m .
coordinates representation .
however , the calculation of these matrix elements is difficult because the term @xmath172 can not be separated into a sum of operators written in relative and c.m . coordinates .
therefore , as have made in the previous work @xcite , we employ the angle - average approximation and have @xmath173 where , as was assumed previously , @xmath174 for @xmath175 .
the @xmath176 is the h.o .
transformation bracket @xcite for the @xmath12 two - body state which depends on the rest masses for @xmath177 and @xmath33 . here
the angle - average approximation means that , when a function @xmath178 is given with an angular momentum @xmath179 and its @xmath180-component @xmath181 , we approximate @xmath178 as @xmath182 . in the present case
the matrix element on the left - hand side in eq .
( [ eq : uyn ] ) is dependent on the @xmath180-components @xmath183 and @xmath184 of @xmath148 and @xmath171 , respectively .
the average over @xmath183 and @xmath184 leads to the expression given in eq .
( [ eq : uyn ] ) .
this approximation is essentially the same as has been made to derive an angle - averaged pauli operator by wong and sauer @xcite .
the use of the basis states @xmath185 makes it difficult to represent the operator @xmath155 , but it has been confirmed that the angle - average technique works well .
we thus make an approximation to the operator @xmath155 as @xmath186 where the summation should be made on the conditions @xmath187 here the coefficient @xmath188 is defined as @xmath189 where in the summation the following conditions should be satisfied @xmath190 furthermore , we need an expression of @xmath191 given in eq .
( [ eq : qynd ] ) in the angle - average approximation @xmath192 where the summation should be made on the same conditions as in eq .
( [ eq : sumc1 ] ) .
the coefficient @xmath193 is given by @xmath194 where @xmath195 and @xmath196 should satisfy the conditions given in eq .
( [ eq : sumc2 ] ) .    once the set of the eigenstates @xmath197 in eq .
( [ eq : eige1r ] ) is solved , the @xmath8 effective interaction @xmath198 acting in the @xmath97 space is given according to eqs .
( [ eq : evln])-([eq : nhevln ] ) by substituting @xmath199 , where @xmath200 is defined as latexmath:[\[\begin{aligned } \label{eq : oyn1 } \hspace{-5mm}\omega _ { \sst { \it yn}}^{\sst ( 1 ) } = \sum_{k_{1},ll'ss'j_{r},n_{c}l_{c } } \overline{q}_{\sst { \it yn}}^{\sst ( 1 ) }    \langle \overline{{\it yn};k_{1},ll'ss'j_{r},n_{c}l_{c } }    the @xmath8 effective interaction @xmath157 thus determined , we calculate a new single - particle potential @xmath202 as an averaged matrix element of @xmath203 .
the @xmath202 can be given explicitly in terms of the effective interaction @xmath157 as @xmath204 in the above treatment of the single - particle potential of @xmath32 in the first - step calculation , we do not consider the one - body spin - orbit splitting .
the @xmath8 effective interaction @xmath205 in eq .
( [ eq : tvln1 ] ) is dependent on the single - particle potential @xmath133 , and @xmath133 is determined with @xmath205 as in eq .
( [ eq : ul1 ] ) .
therefore , the equation for determining @xmath205 and @xmath133 should be solved simultaneously .
the model space @xmath91 in the first - step procedure is chosen as a _ low - momentum _
space , but the @xmath91 space is still large when we wish to use @xmath205 in the calculation of bound or low excited states of @xmath32 hypernuclei . with @xmath205 ,
we further proceed to the calculation of the @xmath8 effective interaction acting in a smaller model space and correspondingly the single - particle potential of @xmath32 in low - lying states .
we separate the @xmath91 space into the @xmath206 , @xmath207 , @xmath208 and @xmath209 spaces specified by the numbers @xmath140 , @xmath106 , and @xmath210 as shown in fig .  2 .
the number @xmath211 is introduced to specify the uppermost bound state of @xmath32 . here
these spaces are defined as @xmath212 the @xmath213 space is the excluded space due to the pauli principle for nucleons .
we are principally interested in the @xmath8 effective interaction that gives a description of the bound states of @xmath83 .
it would be desirable that such a @xmath8 effective interaction satisfies the condition of decoupling between two @xmath8 states consisting of the bound and unbound states of @xmath83 .
this requirement is equivalent to the situation that the @xmath8 effective interactions should be decoupled between two @xmath8 states in the @xmath214 and @xmath209 spaces .
therefore , we solve the decoupling equation for @xmath215 between the @xmath8 states in the @xmath214 and @xmath209 spaces .
the space @xmath208 is not a space to be excluded due to the pauli principle .
however , some states in the @xmath208 and @xmath214 spaces overlap each other in energy . therefore
, we treat separately the effect of the coupling between the @xmath208 and @xmath214 spaces .
actually the interaction between states in the @xmath208 and @xmath214 spaces induces the core polarization in the core nucleus , as discussed in refs .
@xcite and @xcite .
we take into account the effect of the core polarization in the usual perturbative calculation .
= 6 cm    in the second - step calculation the unitarily transformed interaction is given by @xmath216 where @xmath217 the term @xmath218 is the self - consistent potential of @xmath32 to be determined in the second - step procedure .
the decoupling equation for determining the correlation operator @xmath219 in eq .
( [ eq : tvyn2 ] ) becomes @xmath220 we note here that the assumption of decoupling @xmath221 is satisfied exactly because the two two - body states in the @xmath214 and @xmath209 spaces do not contain a common single - particle state of @xmath83 or @xmath19 .
the two - body eigenvalue equation to be solved in the second - step calculation is written as @xmath222 where @xmath223 denotes the total angular momentum , and @xmath224 a set of quantum numbers of a @xmath8 two - body eigenstate .
it should be noted that this eigenvalue equation can be solved exactly in the @xmath225 space by the diagonalization of the hamiltonian written with the basis states @xmath226 if the @xmath227 space is chosen suitably as a small space .
once a set of the eigenstates @xmath228 in eq .
( [ eq : eige2 ] ) is solved , the @xmath8 effective interaction @xmath229 in the @xmath214 space is given from eqs .
( [ eq : evln])-([eq : nhevln ] ) with @xmath230 defined by @xmath231 is calculated with the @xmath8 effective interactions , @xmath229 and @xmath157 , as @xmath232 we solve the set of equations for determining @xmath229 and @xmath233 iteratively until the result converges .
we remark here that we employ @xmath157 as the effective interaction in the @xmath213 space in the calculation of @xmath233 . since the effective interaction in the @xmath213 space is unchanged in the transformation made in the second - step procedure , we may write the @xmath8 effective interaction in the @xmath213 space as @xmath234 .
we , therefore , have used partly the @xmath87 effective interaction @xmath157 in the calculation of @xmath233 as in eq .
( [ eq : ul2 ] ) .
first we note that the one - body hamiltonian @xmath235 in eq .
( [ eq : hln2 ] ) has non - diagonal terms because we did not make the hartree - fock calculation .
these terms could bring about sizable corrections to the single - particle energies of @xmath236 .
we adopt the same h.o .
energy both for @xmath83 and @xmath19 , though it has been considered that the wave function of @xmath83 is basically different from that of @xmath237 .
therefore , we must evaluate the corrections arising from the non - diagonal terms of @xmath235 . by taking the diagonal part of @xmath235 as the unperturbed hamiltonian
we calculate these corrections perturbatively up to second order as shown in diagram ( a ) of fig .  3 .    next
, the interaction between the @xmath8 states in the @xmath214 and @xmath208 spaces induces the core polarization in the core nucleus , and is important in the calculation of the corrections .
we here assume that these corrections can be calculated approximately , using the interaction @xmath238 it will be desirable to take into account the effect of intermediate nucleon states with high momentum . in this perturbative calculation
we include only the nucleon states belonging to the @xmath208 space , that is , the nucleon h.o .
state should satisfy @xmath239 , where @xmath240 are the h.o .
quantum numbers of a @xmath83 state . the validity of this truncation of nucleon states will be checked in the numerical calculation in the next section by increasing the number @xmath241 .
the diagrams corresponding these corrections are given in diagram ( b ) and ( c ) of fig .
we performed calculations employing the jlich model -  ( j ) , model - b ( jb ) @xcite potentials in the momentum - space representation made by b.  holzenkamp and a.  reuber . ] and nijmegen soft - core ( nsc ) @xcite potentials for the _ yn _ interaction .
we used the h.o .
basis states with a common @xmath242mev for @xmath32 , @xmath33 and @xmath19 .
the dependence of the effective interaction and single - particle energies on the value @xmath100 has been already discussed in @xmath15o @xcite and we confirmed that the dependence was quite small around @xmath243mev if perturbation corrections were included . as for the nucleon single - particle potential @xmath244 , we used the fixed data @xcite calculated using the paris potential @xcite , and we have assumed that @xmath245 .
the boundary numbers @xmath246 , @xmath247 and @xmath109 are taken as @xmath248 , @xmath249 and @xmath250 for @xmath0o .
= 13.8 cm    we present , in fig .  4 , the @xmath93 dependence of the single - particle energies of @xmath32 including the perturbation corrections shown in fig .  3 .
the parameter @xmath93 specifies the @xmath91 space in which the @xmath8 effective interaction @xmath229 is determined rigorously without any approximation as discussed in sec .
3 . therefore , as @xmath140 becomes larger , the calculated result does more accurate .
if the result converges at a certain number @xmath251 , we may say that the method of the two - step calculation of the @xmath8 effective interaction would be trustworthy . as shown in fig .  4 ,
the results are convergent and stable for @xmath252 for each @xmath12 potential .
if we take @xmath93 larger than 6 , our model space is considered to be sufficiently large .
we note that the model space with @xmath253 corresponds to that used in the @xmath13-matrix calculations made by halderson @xcite and by hao and kuo @xcite .
.dependence of typical diagonal matrix elements @xmath254 on @xmath255 in units of mev , where @xmath256 and @xmath257 denote the @xmath32 and @xmath19 single - particle orbits , respectively , labeled as @xmath258 , @xmath259 and @xmath260 , and @xmath223 is the total angular momentum .
[ cols="^,^,^,^,^,^,^,^ " , ]     next , the partial - wave contributions to the first - order potential energy of @xmath32 in the 0@xmath1 state are given for the different treatments of the pauli - blocking effect .
we see from table  4 that the contributions of the @xmath261 and @xmath262-@xmath263 channels are very different dependently on the @xmath5 interactions employed . in the @xmath261
channel the largest pauli - blocking effect is brought about in the result for the j potential . on the other hand , in the @xmath262-@xmath263 channel we have an extremely large pauli - blocking effect for the nsc potential .
this result shows that the nsc potential has the strong @xmath7-@xmath8 coupling in the @xmath262-@xmath263 channel as discussed in ref .
the absolute values contributing to the first - order potential energy of @xmath32 in the channels with angular momenta @xmath264 are much smaller than those in the @xmath261 and @xmath262-@xmath263 channels .    in table  5
we display the single - particle energies of @xmath32 with the contributions of the perturbation corrections given in fig .  3 , where the pauli - blocking effect is treated accordingly to case ( 3 ) .
from this table we see that the j and nsc potentials yield almost the same values of the first - order potential energy pe and the perturbation corrections in spite of having different partial - wave contributions as shown in table  4 .
the single - particle energies obtained for the jb potential are most attractive among the @xmath265 interactions used .
it should be noted here that the values of the ke plus pe in the 0@xmath266 states for three kinds of @xmath5 interactions are positive ( unbound ) .
therefore , the perturbation corrections play a decisive role when we argue that the @xmath32 in the 0@xmath266 states is bound or not .
correction ( a ) in fig .
3 is induced by the presence of non - diagonal terms in the one - body hamiltonian @xmath267 .
this effect should be treated in the hartree - fock ( hf ) calculation . however , in the present calculation , we treat this effect perturbatively .
these non - diagonal terms work to change the single - particle wave function of @xmath32 .
therefore , this procedure is indispensable in this approach since we have taken a fixed value of @xmath100 commonly for @xmath32 and @xmath19 .
this correction  ( a ) contributes most attractively to the single - particle energies among all the corrections in fig .
the magnitude of every correction is considerably small in comparison with the first - order potential energies for all the @xmath5 interactions used .
all the perturbation corrections contribute , however , attractively to the single - particle energies of @xmath32 , and do not change the mutual spacings of the single - particle states significantly .
this trend is very different from the situation observed in ordinary nuclei such as @xmath15o and @xmath16ca , in which the correction terms contribute to about 40@xmath268 of the calculated spin - orbit splitting of the 0@xmath266 or 0@xmath269 states @xcite .
the role of many - body correlations , such as the core polarization , does not seem to be important in the spin - orbit splitting of @xmath32 .
we have presented a formulation of the umoa to apply it to a structure calculation of hypernuclei .
we have introduced a unitary transformation exp(@xmath270 ) , with the correlation operator @xmath270 , that describes two - body correlation of hyperon and nucleon .
an equation has been given for determining the correlation operator @xmath270 as the equation of decoupling so that the unitarily transformed interaction does not have non - zero matrix elements between the model space @xmath88 of low - lying @xmath8 states and its complement @xmath85 .
the effective interaction has been given by the projected interaction onto the @xmath88 space .
the effective interaction thus determined has properties of being hermitian , energy - independent and decoupled between the @xmath88 and @xmath85 spaces . in comparison with the usual @xmath13-matrix approach , the use of such an effective interaction would have some advantages in describing many - body systems .
the umoa has been applied to the calculation of the structure of @xmath0o .
the @xmath8 effective interaction has been calculated by employing three @xmath12 potentials , namely , the j , jb and nsc potentials . with the @xmath8 effective interaction ,
we have calculated the single - particle energies of @xmath32 in the 0@xmath1 , 0@xmath2 and 0@xmath3 states .
we have seen that how to treat the pauli - blocking effect is very important , especially when we solve the problem of the @xmath7-@xmath8 coupling . in the usual treatment one
excludes all the interactions acting in a space of @xmath271 and @xmath272 states in which a nucleon is in occupied states .
however , this treatment of the pauli - blocking effect would lead to overmuch counting .
we have introduced a new treatment in which , while we exclude the interactions acting among forbidden @xmath5 states as in the usual treatment , we restore the interactions that are diagonal in occupied nucleon states .
we have observed that this treatment of the pauli - blocking effect has given rise to a significant effect as compared with the usual treatment .
the present results have shown that @xmath32 in the 0@xmath3 state is unbound for the j and nsc potentials and , on the other hand , bound for the jb potential .
the spin - orbit splitting of the single - particle levels of @xmath83 in the 0@xmath266 states has been given by 1.41 , 1.57 and 1.66 for the j , jb and nsc potentials , respectively .
these values are very large compared with the results in ref .
it will be of high interest for the single - particle nature of @xmath32 in @xmath0o to be established experimentally .    in the application of the umoa
some problems have still remained .
the unitarily transformed hamiltonian @xcite contains originally three - or - more - body interactions .
the effects of these many - body interactions have been evaluated in the previous works of calculating nuclear properties , and we have concluded that the contributions of these many - body terms are much smaller than those of two - body interactions but sizable in some cases .
we therefore should estimate the effects of these many - body interactions also in hypernuclei .
another refinement to be made is that , in the present calculation , the single - particle energies of @xmath33 in intermediate states are neglected .
however , the single - particle spectrum of @xmath33 has a possibility of giving rise to some contributions to the @xmath8 effective interaction .
although the properties of the single - particle energies of @xmath33 have not always been made clear , we should investigate how much the spectrum of @xmath33 affects in the determination of the @xmath8 effective interaction .
this problem is also a remaining task in the umoa .
the authors would like to thank y.  yamamoto and k.  miyagawa for useful discussions and their help for using the nijmegen and jlich @xmath265 potential codes .
one of the authors ( s.  f. ) would like to express his thank to k.  takada , m.  kamimura and y.  r.  shimizu for their continuous encouragement and instructive discussions .
a.  gal , advances in nuclear physics , vol .  8 , eds . m.  baranger and e.  vogt ( plenum , new york , 1975 ) p.  1 .
b.  povh , ann .
nuc .  part .
28 ( 1978 ) 1 .
h.  band et al . , prog .
81 ( 1985 ) .
h.  band , t.  motoba and j.  @xmath273ofka , int .  j.  mod .
a  5 ( 1990 ) 4021 , and references therein . t.  motoba , y.  akaishi and k.  ikeda , prog .
117 ( 1994 ) .
b.  f.  gibson and e.  v.  hungerford iii , phys .
257 ( 1995 ) 349 .
b.  f.  gibson , i.  r.  afnan , j.  a.  carlson and d.  r.  lehman , prog .
117 ( 1994 ) 339 .
w.  brckner et al . , phys .
. 79  b ( 1976 ) 157 .
r.  h.  dalitz , d.  h.  davis , t.  motoba and d.  n.  tovee , nucl .
phys . a  625 ( 1997 ) 71 .
k.  miyagawa , and w.  glckle , phys .
rev . c  48 ( 1993 ) 2576 ; + k.  miyagawa , h.  kamada , w.  glckle and v.  stoks , phys .
rev . c  51 ( 1995 ) 2905 . t.  yamada , h.  ohkura , t.  motoba and k.  ikeda , prog .
117 ( 1994 ) 65 .
e.  hiyama , m.  kamimura , t.  motoba , t.  yamada and y.  yamamoto , prog .
97 ( 1997 ) 881 . h.  band and y.  yamamoto , prog .
theor .  phys .
81 ( 1985 ) 9 .
y.  yamamoto , t.  motoba , h.  himeno , k.  ikeda and s.  nagata , prog .
117 ( 1994 ) 361 .
h.  band , prog .
66 ( 1981 ) 1349 .
d.  halderson , phys .
rev . c  48 ( 1993 ) 581 .
j.  hao , t.  t.  s.  kuo , a.  reuber , k.  holinde , j.  speth and d.  j.  millener , phys .
71 ( 1993 ) 1498 ; + t.  t.  s.  kuo and j.  hao , prog .
117 ( 1994 ) 351 .
m.  hjorth - jensen , a.  polls , a.  ramos and h.  mther , nucl .
phys . a  605 ( 1996 ) 458 .
i.  vidaa , a.  polls , a.  ramos and m.  hjorth - jensen , nucl .
phys . a  644 ( 1998 ) 201 .
t.  t.  s.  kuo and e.  osnes , in : folded - diagram theory of the effective interaction in nuclei , atoms and molecules , lecture notes in physics , vol .
264 , ( springer , berlin , 1990 ) .
j.  hao and t.  t.  s.  kuo , phys .
264 ( 1996 ) 233 .
k.  suzuki , r.  okamoto and h.  kumagai , phys .
c  36 ( 1987 ) 804 ; + k.  suzuki and r.  okamoto , prog .  theor .
( 1994 ) 1045 .
h.  kumagai , k.  suzuki and r.  okamoto , prog .
97 ( 1997 ) 1023 .
s.  fujii , r.  okamoto and k.  suzuki , prog .
99 ( 1998 ) 151 . m.  m.  nagels , t.  a.  rijken and j.  j.  de  swart , phys .
rev . d  15 ( 1977 ) 2547 ; + p.
m.  m.  maessen , t.  a.  rijken and j.  j.  de  swart , phys .
c  40 ( 1989 ) 2226 .
b.  holzenkamp , k.  holinde and j.  speth , nucl .
phys . a  500 ( 1989 ) 485 ; + a.  reuber , k.  holinde and j.  speth , nucl .
phys . a  570 ( 1994 ) 543 .
m.  lacombe , b.  loiseau , j.  m.  richard , r.  vinh  mau , j.  ct , p.  pirs and r.  de  tourreil , phys .
c  21 ( 1980 ) 861 . b.  h.  brandow , in : effective interactions and operators in nuclei , lecture notes in physics , vol .  40 , ed . b.  r.  barrett ( springer , berlin , 1975 ) p.  1 .
i.  shavitt and l.  t.  redmon , j. chem .
73 ( 1980 ) 5711 .
p.  westhaus , int .
j. quantum chem . 20
( 1981 ) 1243 .
k.  suzuki , prog .
68 ( 1982 ) 246 .
t.  t.  s.  kuo , p.  j.  ellis , j.  hao , z.  li , k.  suzuki , r.  okamoto and h.  kumagai , nucl .
phys . a  560 ( 1993 ) 621 .
y.  yamamoto and h.  band , prog .
83 ( 1990 ) 254 .
schulze , a.  lejeune , j.  cugnon , m.  baldo and u.  lombardo , phys .
b  355 ( 1995 ) 21 ; phys .
c  57 ( 1998 ) 704 .
k.  suzuki and r.  okamoto , prog .
76 ( 1986 ) 127 .
e.  m.  krenciglowa , c.  l.  kung , t.  t.  s.  kuo and e.  osnes , ann .  of  phys . 101
( 1976 ) 154 .
b.  buck and a.  c.  merchant , nucl .
a  600 ( 1996 ) 387 . c.  w.  wong , nucl .
phys . a  91 ( 1967 ) 399 .
p.  u.  sauer , nucl .
phys . a  150 ( 1970 ) 467 .