trapped ultra - cold atoms have become model systems of choice for simulating physical effects from condensed matter  @xcite to cosmology  @xcite . the recently achieved ability to create synthetic vector potentials  @xcite acting on neutral atoms
has increased the versatility of the atomic - physics simulation toolkit even further .
it is now possible to simulate magnetic fields by inducing spatially varying @xmath0 ( i.e. , scalar ) gauge potentials  @xcite , and pseudo - spin splittings can be created in spinor gases  @xcite using spatially constant vector potentials having a ( possibly non - abelian ) matrix structure .
these advances have stimulated a host of theoretical works studying , e.g. , the effect of uniform @xmath1 gauge potentials on the behavior of quantum particles subject to uniform ordinary magnetic fields  @xcite , or proposing the use of staggered effective spin - dependent magnetic fields in optical lattices  @xcite to simulate a new class of materials called topological insulators  @xcite that exhibit the quantum spin hall effect  @xcite . very recently , the non - quantized intrinsic spin hall effect  @xcite has been realized experimentally in a quantum gas  @xcite , and the authors of this paper outline the way forward to reaching conditions where the quantum spin hall effect could be observed .
furthermore , newly demonstrated methods to simulate strong - enough magnetic fields to probe ultra - cold atom gases in the ordinary quantum - hall regime  @xcite are expected to be adaptable for the purpose of generating spin - dependent quantizing magnetic fields  @xcite , which opens up another avenue towards the exploration of quantum - spin - hall physics .
part of the motivation for our present theoretical work arises from these rapid developments of experimental capabilities .
the ordinary quantum hall ( qh ) effect  @xcite occurs because particles confined to move in two spatial dimensions and subject to a strong perpendicular magnetic field develop incompressibilities at integer , and certain fractional , values of the landau - level filling factor  @xcite .
while the landau quantization of single - particle energies is the origin of the integer qh effect , incompressibility at fractional filling factors is caused by the discrete spectrum of interaction energies for two particles occupying states from the same landau level  @xcite . in the conceptually simplest realization of the quantum _
spin _ hall ( qsh ) effect  @xcite , particles exhibit an integer qh effect due to a spin - dependent perpendicular magnetic field that points in opposite directions for the two opposite - spin components .
it seems then quite straightforward to conjecture  @xcite that a fractional version of the qsh effect should exist that mirrors features of the ordinary fractional qh effect in multi - component systems  @xcite . following this line of thought ,
some previous discussions of a putative fractional qsh physics  @xcite have been based on an _ ad hoc _ adaptation of trial wave functions first proposed in ref .
however , unless only particles with the same spin interact , such an approach is fraught with difficulty  @xcite .
here we revisit the question of how a fractional qsh effect can arise in an interacting ( pseudo-)spin-1/2 system that experiences a spin - dependent quantizing magnetic field .
in particular , we elucidate the effect of interactions between particles having opposite spin .
we find that such inter - species interactions significantly alter the expected qsh physics , but they also open up new opportunities for tailoring the properties of quantum many - particle states .
our study is complementary to recent investigations of fractional qsh phases  @xcite that arise in materials with exotic topological band structures  @xcite or strained graphene  @xcite .
the results obtained here are relevant for electronic systems as well as for ultra - cold bosonic or fermionic atoms .
particular attention is paid to trapped bosons .
the article is organized as follows . in sec .
[ sec : hamiltonian ] , we introduce the basic model description of an interacting system of ( pseudo-)spin-1/2 particles that are subject to a spin - dependent magnetic field .
the single - particle states are given in the representation of spin - dependent guiding - center and landau - level quantum numbers .
the eigenvalue problem of two interacting particles is solved  for both cases of equal and opposite - spin particles  in the subsequent sec .  [
sec : twopart ] .
when the two interacting particles have opposite spin , important differences arise with respect to the classic results obtained  @xcite for spinless ( or same - spin ) particles .
we explore the ramifications of this fact by numerical exact - diagonalization studies with up to 6 bosons for which results are presented in sec .
[ sec : numfew ] .
our conclusions are summarized in sec .
[ sec : conc ] .
we consider a gas of particles ( e.g. , atoms ) that carry a ( pseudo-)spin-1/2 degree of freedom and are confined to move in the @xmath2 plane .
a spin - dependent vector potential @xmath3 is presumed to be generated ( e.g. , by optical means in an atom gas  @xcite ) . here
@xmath4 denotes the diagonal pauli matrix , and the @xmath5 are cartesian unit vectors in real space .
the vector potential ( [ eq : acgauge ] ) is abelian and gives rise to a spin - dependent magnetic field perpendicular to the @xmath2 plane : @xmath6 .
this situation of opposite - spin particles being subjected to oppositely directed magnetic fields corresponds directly to setups considered for a semiconductor heterostructure  @xcite and in neutral - atom systems  @xcite .
a somewhat related study in the context of cold bosonic gases was given in ref .
@xcite , only that there the two spin components also experience a large zeeman - like energy shift and , therefore , this work focused only on the dynamics of a single component .    without loss of generality
, we will assume @xmath7 from now on . also note that , with unit conventions chosen in this article , the ` magnetic - field ' magnitude @xmath8 is related to a fundamental ( ` magnetic ' ) length scale @xmath9 .
the most general form of the many - body hamiltonian that describes our system of interest is @xmath10 , where @xmath11 ^ 2 + { \mathcal v}({{\ensuremath{\bm{\mathrm{r } } } } } ) -\mu \ , \hat 1 \right\ } \hat\psi({{\ensuremath{\bm{\mathrm{r } } } } } ) \quad , \\ \label{eq : model_int } { \mathcal h}_{\mathrm{int } } & = & \int d^2 r \,\
, \sum_{\sigma , \sigma^\prime = \pm}\ , \frac{g_{\sigma\sigma^\prime}}{2 } \ , \hat\psi_{\sigma}^\dagger({{\ensuremath{\bm{\mathrm{r } } } } } ) \hat\psi_{\sigma^\prime}^\dagger({{\ensuremath{\bm{\mathrm{r } } } } } ) \hat\psi_{\sigma^\prime}({{\ensuremath{\bm{\mathrm{r}}}}})\hat\psi_{\sigma}({{\ensuremath{\bm{\mathrm{r } } } } } ) \quad .\end{aligned}\ ] ] here @xmath12 denotes the @xmath13 unit matrix , and @xmath14 is a two - spinor of ( bosonic or fermionic ) annihilation operators for a particle located at position @xmath15 and having a definite pseudo - spin-1/2 projection along the @xmath16 axis .
furthermore , the relation @xmath17 is implicit in the formalism .
as it is often useful , we also give an equivalent expression for the interaction part of the hamiltonian , @xmath18 with @xmath19 , @xmath20 , and @xmath21 . depending on whether the asymmetric interaction couplings @xmath22 , @xmath23 are positive or negative , different physical regimes may be reached  @xcite.)([eq : ints_2 ] )
can be related to that which is often adopted in the atom - gas literature  @xcite by setting @xmath24 , @xmath25 , and @xmath26 .
see also ref .
note , however , the different parameterization used in ref .
@xcite where @xmath27 are interaction constants associated with the atomic spin-1 degree of freedom from which the pseudo - spin-1/2 components are derived .
our notation is related to theirs via @xmath28 , @xmath29 , and @xmath30 . ]    in the case where @xmath31 , the system reduces to two independent two - dimensional ( electron or atom ) gases that are each subject to a perpendicular magnetic field .
known phenomena associated with the fractional qh effect  @xcite will then be exhibited by the individual systems . however , as seen from our study presented in secs .
[ sec : twopart ] and [ sec : numfew ] below , the behavior of the system with @xmath32 departs from the previously considered  @xcite two - component fractional - qh physics because of the very different type of constraints that is placed on the orbital motion of particles subject to oppositely directed magnetic fields .
we now consider single - particle states associated with spin component @xmath33 .
the kinetic momentum is @xmath34 , and straightforward calculation establishes the commutation relations @xmath35=\sigma\ , i \,\frac{\hbar^2}{l_{\mathcal b}^2 } \ , \delta_{\sigma\sigma'}\ , \epsilon_{\alpha\alpha ' } \quad , \ ] ] where @xmath36 .
operators for the guiding - center locations can then be defined in the usual manner  @xcite , @xmath37/\hbar$ ] , and their components satisfy the commutation relations @xmath38=-\sigma\ , i \ , l_{\mathcal b}^2 \
, \delta_{\sigma\sigma'}\,\epsilon_{\alpha\alpha ' } \quad .\ ] ] moreover , we find @xmath39=0 $ ] .
following the familiar approach  @xcite , we define the harmonic - oscillator landau - level ladder operator for states with spin @xmath33 via @xmath40 similarly , the ladder operator operating within a landau level for spin component @xmath33 is @xmath41    we can express the kinetic energy and the @xmath42 component of angular momentum in terms of the ladder operators [ @xmath43 : @xmath44 ^ 2}{2 m } & = & \hbar\omega_{\mathrm{c } } \left ( a_\sigma^\dagger a_\sigma + \frac{1}{2 } \right ) \quad , \\
\label{eq : angmom } l_z & \equiv & x\ , p_y - y\ , p_x = \hbar\ , \sigma \left ( b_\sigma^\dagger b_\sigma - a_\sigma^\dagger a_\sigma \right ) \quad .\end{aligned}\ ] ] landau - level eigenstates are generated via @xmath45 { \left | \mathrm{vac } \right \rangle } \,\
, , \ ] ] where @xmath46 and @xmath47 is the state that is annihilated by all ladder operators @xmath48 and @xmath49 . in the following ,
we will focus on the case where all particles are in the lowest landau level , i.e. , when @xmath50 .
before presenting a formal analysis of the interacting two - particle system subject to a strong spin - dependent magnetic field in the following subsection , we provide a heuristic argument for how the cases where the two particles feel the same and opposite magnetic fields differ .
consider two particles , located at @xmath51 and @xmath52 , respectively , that interact via a generic potential @xmath53 .
the corresponding first - quantized two - particle hamiltonian reads @xmath54
^ 2}{2 m }   + \frac{\left [ { { \ensuremath{\bm{\mathrm{p}}}}}_2 - { { \ensuremath{\bm{\mathrm{\cal{a}}}}}}^{(\sigma_2)}({{\ensuremath{\bm{\mathrm{r}}}}}_{2})\right]^2}{2 m } + v({{\ensuremath{\bm{\mathrm{r}}}}}_1 - { { \ensuremath{\bm{\mathrm{r}}}}}_2 ) \quad , \ ] ] with the spin - dependent vector potentials from eq .
( [ eq : acgauge ] ) .
it is straightforward to show that the sum of kinetic - energy contributions for each particle can be re - arranged in terms of the linear combinations @xmath55 yielding the expression @xmath56 ^
2}{2 m } + \frac{\left[{{\ensuremath{\bm{\mathrm{p}}}}}_{\sigma_1 \sigma_2 } - 2 { { \ensuremath{\bm{\mathrm{\cal{a}}}}}}^{(+ ) } ( { { \ensuremath{\bm{\mathrm{r}}}}}_{\sigma_1\sigma_2})\right]^2}{2 m}+v({{\ensuremath{\bm{\mathrm{r}}}}}_1 - { { \ensuremath{\bm{\mathrm{r}}}}}_2 ) .
\nonumber \\\end{aligned}\ ] ]    the two - particle problem for particles with the same spin reduces to two independent single - particle problems in the center - of - mass ( com ) and relative - coordinate degrees of freedom ( @xmath57 and @xmath58 , respectively ) because @xmath59 . in that case , only the relative - coordinate degree of freedom feels the interaction potential @xmath60 , and it can be minimized by placing two particles away from each other . in the specific case of @xmath61 ,
the zero - energy states of the two - particle system are of the form @xmath62 , where @xmath63 is a commonly used complex notation for the position of particle @xmath64  @xcite .
the non - negative integers @xmath65 and @xmath66 correspond to the quantized values of com angular momentum and relative angular momentum , respectively  @xcite .
quite a different situation arises for opposite - spin particles . the kinetic energy of the two - particle system decouples in the coordinates @xmath67 and @xmath68 , motivating the proposal of trial wave functions  @xcite @xmath69 .
however , @xmath70 still couples the two - particle coordinates @xmath67 and @xmath68 and , as a result , the proposed wave function is energetically not favorable for interacting particles  @xcite .
rigorous examination of the interacting two - particle system in the opposite - spin configuration ( see below ) shows that energy eigenstates are not eigenstates of com angular momentum or relative angular momentum and , furthermore , have an unusual distribution .      to gain a deeper understanding of the effect of two - particle interactions , we follow the basic approach employed by previous studies of the fractional qh effect  @xcite and find the interaction potential in the representation of lowest - landau - level states .
the starting point of such an analysis is the fourier decomposition of a spin - dependent interaction potential given by @xmath71 because its matrix elements can then be directly related to the corresponding matrix elements of the exponential in the integrand of ( [ eq : fourierpot ] ) . to find the latter for states in the lowest landau level
, we use the relation @xmath72 and the straightforwardly obtained expressions @xmath73 \ , , \\   & & \frac{l_{\mathcal b}^2}{\hbar}\left({{\ensuremath{\bm{\mathrm{q } } } } } \times \hat{{{\ensuremath{\bm{\mathrm{z}}}}}}\right)\cdot \left ( \sigma_1\ , { { \ensuremath{\bm{\mathrm{\pi}}}}}_{1}^{(\sigma_1 ) } - \sigma_2\ , { { \ensuremath{\bm{\mathrm{\pi}}}}}_{2}^{(\sigma_2)}\right ) = \frac{l_{\mathcal b}}{\sqrt{2 } } \left [ q_x \left ( a_{\sigma_1 } - a_{\sigma_2 } + a_{\sigma_1}^\dagger - a_{\sigma_2}^\dagger\right ) \right .
\nonumber \\ & & \hspace{4.5 cm } \left .
+ \sigma_1 i \ ,
q_y \left ( a_{\sigma_1 }   - \sigma_1\sigma_2\ , a_{\sigma_2 } - a_{\sigma_1}^\dagger + \sigma_1 \sigma_2 \ ,   a_{\sigma_2}^\dagger\right ) \right ] .
\end{aligned}\ ] ] introducing spin - resolved ladder operators for com and relative angular momentum , @xmath74 and analogous ladder operators for com energy and relative - motion energy , @xmath75 we find using @xmath76 @xmath77 \times\left\ { \begin{array}{cl } { \ensuremath{\mathrm{e}}}^{i \bar q_\sigma l_{\mathcal b}\ , b^{\mathrm{(r)}\dagger}_{\sigma\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{i q_\sigma l_{\mathcal b}\ , b^{\mathrm{(r)}}_{\sigma\sigma } } & \mbox{for}\ , \sigma_1 = \sigma_2 \equiv \sigma \\ { \ensuremath{\mathrm{e}}}^{i q_x l_{\mathcal b}\ , b^{\mathrm{(r)}\dagger}_{\sigma , -\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{i q_x l_{\mathcal b } \ , b^{\mathrm{(r)}}_{\sigma,-\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{\sigma q_y l_{\mathcal b}\ , b^{\mathrm{(c)}\dagger}_{\sigma,-\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{-\sigma q_y l_{\mathcal b}\ , b^{\mathrm{(c)}}_{\sigma,-\sigma } }   & \mbox{for } \ , \sigma_1= -\sigma_2 \equiv \sigma \end{array } \right . , \\[0.1 cm ] { \ensuremath{\mathrm{e}}}^{i\frac{l_{\mathcal b}^2}{\hbar}\left({{\ensuremath{\bm{\mathrm{q } } } } } \times \hat{{{\ensuremath{\bm{\mathrm{z}}}}}}\right)\cdot \left ( \sigma_1\ , { { \ensuremath{\bm{\mathrm{\pi}}}}}_{1}^{(\sigma_1 ) } - \sigma_2\ , { { \ensuremath{\bm{\mathrm{\pi}}}}}_{2}^{(\sigma_2)}\right ) } = { \ensuremath{\mathrm{e}}}^{-\frac{q^2\ , l_{\mathcal b}^2}{2 } }   \nonumber\\[0.2 cm ] \times \left\ { \begin{array}{cl } { \ensuremath{\mathrm{e}}}^{i \bar q_\sigma l_{\mathcal b}\ , a^{\mathrm{(r)}\dagger}_{\sigma\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{i q_\sigma l_{\mathcal b}\ , a^{\mathrm{(r)}}_{\sigma\sigma } } & \mbox{for}\ , \sigma_1=\sigma_2 \equiv \sigma \\ { \ensuremath{\mathrm{e}}}^{i q_x l_{\mathcal b}\ , a^{\mathrm{(r)}\dagger}_{\sigma,-\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{i q_x l_{\mathcal b } \ , a^{\mathrm{(r)}}_{\sigma,-\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{\sigma q_y l_{\mathcal b}\ , a^{\mathrm{(c)}\dagger}_{\sigma,-\sigma}}\,\ , { \ensuremath{\mathrm{e}}}^{-\sigma q_y l_{\mathcal b}\ , a^{\mathrm{(c)}}_{\sigma,-\sigma } }   & \mbox{for } \ ,
\sigma_1= -\sigma_2 \equiv \sigma \end{array } \right . .
\label{eq : spinexp2}\end{aligned}\ ] ] inspection of eqs .
( [ eq : spinexp1])([eq : spinexp2 ] ) reveals a very important formal difference between cases when the interacting particles have equal or opposite spin . for a pair of particles with the same spin
, the interaction only couples to the relative - motion landau - level degrees of freedom .
in contrast , for particles with opposite spin , the interaction involves both the relative - motion _ and the com _ degrees of freedom .
this latter situation is unlike any other encountered previously in the context of fractional qh physics .    with the expressions ( [ eq : spinexp1])([eq : spinexp2 ] ) , we are now able to express the interaction potential for a pair of particles having spin @xmath78 and @xmath79 , respectively , in the basis of com - angular - momentum and relative - angular - momentum eigenstates from the lowest landau level given by @xmath80^{m_\mathrm{c}}}{\sqrt{m_\mathrm{c } ! } } \ , \frac{\left [ b^{\mathrm{(r)}\dagger}_{\sigma_1\sigma_2 } \right]^{m_\mathrm{r}}}{\sqrt{m_\mathrm{r } ! } } \ , { \left | 0 \right \rangle}_1\otimes { \left | 0 \right \rangle}_2 \quad .\ ] ] for clarity , the cases where the interacting particles have equal or opposite spin will be discussed separately .
the fact that only the relative - angular - momentum operators enter the expression ( [ eq : spinexp1 ] ) for @xmath81 implies that the interaction matrix is diagonal in com space .
straightforward calculation yields @xmath82 where @xmath83 in terms of the generalized laguerre polynomial @xmath84 . using the result ( [ eq : spinlike ] ) and the relation ( [ eq : fourierpot ] ) for a contact interaction where @xmath85 yields the well - known expression  @xcite @xmath86 for the interaction matrix elements .
the remarkable result ( [ eq : usualpp ] ) underpins the basic description of fractional - qh physics  @xcite .
it implies that the two - particle eigenstates are also eigenstates of com and relative angular momentum .
furthermore , the energy spectrum of two particles with spin @xmath33 from the lowest landau level is two - valued : states with @xmath87 and arbitrary @xmath65 have energy @xmath88 , and all other states have energy @xmath89 .      using ( [ eq : spinexp1 ] ) for the case @xmath90 , we find @xmath91 the contact - interaction matrix element for opposite - spin particles is then calculated as @xmath92 with @xmath93}{2\sqrt{m!\ , m ' ! } } \,\ , \gamma\left ( \frac{m+m'+1}{2}\right)\ ] ] in terms of the euler gamma function @xmath94 .
thus we find that the interaction matrix for two particles from the lowest landau level with opposite spin is nondiagonal in the com - angular - momentum and relative - angular - momentum spaces .
this is markedly different from the case of same - spin particles .
panel a : eigenvalues @xmath95 of the opposite - spin two - particle interaction matrix [ cf.eq .
( [ eq : oppspin ] ) ] in units of @xmath96 , sorted by magnitude .
data are shown for various values of the angular - momentum cutoff @xmath97 ( blue ) , 20 ( red ) , 30 ( green ) , and @xmath98 .
the straight line is a plot of @xmath99 with @xmath100 .
panel b : energy spectrum obtained for a system of two particles with opposite spin by exact diagonalization .
a finite system size is imposed by limiting the number of modes available in angular - momentum space for each particle to @xmath101 .
note the @xmath101-dependence of the obtained values . the data for @xmath102 are also shown as the magenta data points in panel a and exhibit excellent agreement with the power - law - type distribution predicted from the solution in com and relative angular - momentum space .
green stars show the energy calculated for two - particle versions of trial states  @xcite @xmath69 with @xmath103 and @xmath104 , @xmath105 , @xmath106 .
higher @xmath66 is observed to correlate with lower energy , but there are many states even lower in energy than the trial state with largest @xmath66 that is compatible with the finite systems size .
panel c : comparison of two - particle densities of states for same - spin case ( blue arrows indicating delta functions ) and for opposite - spin case ( red curve).,title="fig : " ]   panel a : eigenvalues @xmath95 of the opposite - spin two - particle interaction matrix [ cf.eq .
( [ eq : oppspin ] ) ] in units of @xmath96 , sorted by magnitude .
data are shown for various values of the angular - momentum cutoff @xmath97 ( blue ) , 20 ( red ) , 30 ( green ) , and @xmath98 .
the straight line is a plot of @xmath99 with @xmath100 .
panel b : energy spectrum obtained for a system of two particles with opposite spin by exact diagonalization .
a finite system size is imposed by limiting the number of modes available in angular - momentum space for each particle to @xmath101 .
note the @xmath101-dependence of the obtained values . the data for @xmath102 are also shown as the magenta data points in panel a and exhibit excellent agreement with the power - law - type distribution predicted from the solution in com and relative angular - momentum space .
green stars show the energy calculated for two - particle versions of trial states  @xcite @xmath69 with @xmath103 and @xmath104 , @xmath105 , @xmath106 .
higher @xmath66 is observed to correlate with lower energy , but there are many states even lower in energy than the trial state with largest @xmath66 that is compatible with the finite systems size .
panel c : comparison of two - particle densities of states for same - spin case ( blue arrows indicating delta functions ) and for opposite - spin case ( red curve).,title="fig : " ] + [ 0.2 cm ]   panel a : eigenvalues @xmath95 of the opposite - spin two - particle interaction matrix [ cf.eq .
( [ eq : oppspin ] ) ] in units of @xmath96 , sorted by magnitude .
data are shown for various values of the angular - momentum cutoff @xmath97 ( blue ) , 20 ( red ) , 30 ( green ) , and @xmath98 .
the straight line is a plot of @xmath99 with @xmath100 .
panel b : energy spectrum obtained for a system of two particles with opposite spin by exact diagonalization .
a finite system size is imposed by limiting the number of modes available in angular - momentum space for each particle to @xmath101 .
note the @xmath101-dependence of the obtained values . the data for @xmath102 are also shown as the magenta data points in panel a and exhibit excellent agreement with the power - law - type distribution predicted from the solution in com and relative angular - momentum space .
green stars show the energy calculated for two - particle versions of trial states  @xcite @xmath69 with @xmath103 and @xmath104 , @xmath105 , @xmath106 .
higher @xmath66 is observed to correlate with lower energy , but there are many states even lower in energy than the trial state with largest @xmath66 that is compatible with the finite systems size .
panel c : comparison of two - particle densities of states for same - spin case ( blue arrows indicating delta functions ) and for opposite - spin case ( red curve).,title="fig : " ]    straightforward diagonalization of the matrix ( [ eq : oppspin ] ) yields the two - particle eigenenergies @xmath107 when both particles have opposite spin .
figure  [ 2particlestheo]a shows a logarithmic plot of the @xmath107 , ordered by decreasing magnitude , for different values @xmath108 of the cut - off value for com and relative angular momentum .
we observe an exponential dependence of the sorted eigenvalues as a function of the scaled index @xmath109 , which translates into a power - law density of states @xmath110 here @xmath111 is the area corresponding to the cut - off in com and relative angular momentum , and @xmath112 has been determined numerically .
the numerical data deviate from eq .
( [ eq : dos ] ) close to the maximum energy @xmath113 , where the density of states reaches zero , and for small energy where it becomes cutoff dependent .
none of the individual eigenvalues is strictly independent of the cutoff , which indicates that there are no compact eigenstates .
figure  [ 2particlestheo]c illustrates the different density - of - states behavior for interacting two - particle systems for the two cases of particles having the same and opposite spin , respectively .
it is also useful to look at the distribution of eigenvalues over total angular momentum .
we do this with a different numerical scheme using exact diagonalization of the two - particle hilbert space on a disk , as it preserves the @xmath42 component of angular momentum as a good quantum number .
( details are given in the following section . ) the spectrum for @xmath114 is shown in fig .
[ 2particlestheo]b .
note the dependence of the eigenvalues on the systems size ( ie , the cut - off in angular momentum of available landau - level states ) .
any systematic difference between the results given in figs .
[ 2particlestheo]a and [ 2particlestheo]b is probably at least in part due to the fact that the representation using the com and relative angular - momentum basis assumes an infinite number of single - particle angular - momentum modes to be available to the particles .
nevertheless , when the energy eigenvalues obtained for the finite system are plotted alongside the results for the analytical model ( see magenta data points in fig .  [ 2particlestheo]a ) , both are seen to exhibit the same exponential behavior . for comparision , the energies calculated for proposed trial states  @xcite are also shown in fig .
[ 2particlestheo]b as green stars .
while the interacting two - particle problem has lent itself to analytical study , the behavior of systems with three or more interacting particles either requires approximate , e.g. , variational , approaches or must be done numerically . as the complications encountered already for the case of two interacting particles with opposite spin stymie progress for the variational option
, we follow the numerical route here .
we start by representing the schrdinger field operator for a particle at position @xmath15 with spin @xmath33 projected onto the lowest spin - related landau level , @xmath115 where @xmath116 creates a particle in component @xmath33 with angular momentum @xmath117 in the state @xmath118 . substituting this into eq .
( [ eq : model_int ] ) , we get @xmath119 where @xmath120^\ast \left[\phi^{(\sigma')}_{0 , m_2}({{\ensuremath{\bm{\mathrm{r } } } } } ) \right]^\ast \phi^{(\sigma')}_{0,m_3}({{\ensuremath{\bm{\mathrm{r}}}}})\ , \phi^{(\sigma)}_{0 , m_4}({{\ensuremath{\bm{\mathrm{r } } } } } ) \quad .\ ] ] for same - spin particles , i.e. , @xmath121 , we obtain @xmath122 in contrast , for the matrix element involving opposite - spin particles ( @xmath123 ) , we find @xmath124 the way indices are distributed in the arguments of the @xmath125-functions in eqs .
( [ eq : sameint ] ) and ( [ eq : oppositeint ] ) implies that the system s total angular momentum @xmath126 [ cf .
( [ eq : angmom ] ) for the definition of @xmath127 is a conserved quantity in the presence of interactions .
cold - atom systems are usually studied while trapped by an external potential of tunable strength .
to model this situation , we introduce the second - quantized form of a parabolic potential in the representation of lowest - landau - level states , @xmath128 where @xmath129 in terms of the harmonic - trap frequency @xmath130 .
switching on the trap will lift degeneracies of few - particle states and serve to identify the most compact ground states of our systems of interest .
spectrum for various four - particle systems ( i.e. , @xmath131 ) .
energies are given in units of the intra - species haldane - pseudopotential energy scale @xmath132 .
a : single - component system with @xmath133 , @xmath134 .
the four - particle laughlin state is the zero - energy state with the smallest total angular momentum @xmath135 .
b : system with @xmath136 and @xmath137 , @xmath138 ( no interspecies interaction ) .
c : same situation as for b but with a finite trapping potential ( @xmath139 ) switched on in addition , revealing the energy degeneracies in b. the lowest - energy state is a superposition of two - particle laughlin states in each component .
d : same situation as for b but with finite interspecies interaction @xmath140 in addition.,width=453,height=340 ]    we calculate the few - particle energy spectra and associated eigenstates for @xmath141 in the fock basis of lowest - landau - level states for the two spin components .
we focus here on the case of bosonic particles to be directly applicable to currently studied ultra - cold atom systems , but our general conclusions apply to systems of fermionic particles as well .
results obtained for systems with @xmath131 are shown in fig .
[ fig : n2n2spec ] .
panel a shows the situation where only particles from a single component are present , which is analogous to the previously considered case of spinless bosons @xcite .
the zero - energy state at lowest total angular momentum has @xmath142 and corresponds to the filling - factor-@xmath143 laughlin state  @xcite .
zero - energy eigenstates at higher magnitudes of total angular momentum correspond to edge excitations of the laughlin state  @xcite .
the @xmath144 state has an energy of @xmath145 , where @xmath146 .
when particles occupy states in both components , the situation becomes complex . without interaction between the different spin species , states of each component
will be the ones that are obtained by diagonalising the interacting hamiltonian within that component .
the entire system is then essentially an independent superposition of eigenstates for the individual spin species .
however , in contrast to ordinary multi - component qh states discussed , e.g. , in refs .
@xcite , the total angular momenta for states from different components have opposite sign .
therefore , e.g. , the combination of laughlin states in each component with the same number of particles has zero total angular momentum .
this case is illustrated in fig .
[ fig : n2n2spec]b . however , the superpositions of edge excitations with same magnitude of excess angular momentum for the opposite - spin laughlin states will also be zero - energy , zero - angular - momentum eigenstates . to reveal the associated degeneracies of the spectrum shown in fig .
[ fig : n2n2spec]b , we obtained the energy eigenvalues in the presence of a parabolic confinement .
see fig .
[ fig : n2n2spec]c .
notice the band of low - lying energy levels separated by a gap from higher - energy states .
the lowest - energy @xmath144 state is the superposition of the two - particle laughlin states for the two spin species .
the other states in the low - energy band correspond to edge excitations of this configuration .
figure  [ fig : n2n2spec]d illustrates the dramatic effect of interactions between opposite - spin particles .
the spectrum seen there has to be compared with that given in panel  b where only particles with the same spin interact .
note the disappearance of energy gaps and accumulation of states at low energy , reflecting the characteristic features of the opposite - spin two - particle interaction spectrum shown in fig .
[ 2particlestheo]b . clearly , the system is not incompressible anymore , and no qh - related physics can be expected to occur .
low - lying energy levels for a system with @xmath147 in the sector of total angular momentum @xmath144 . in the calculation ,
lowest - landau - level states with @xmath148 have been included .
different panels correspond to different interspecies - interaction strengths .
a : no inter - species interactions ( @xmath138 ) .
a finite trapping potential lifts the energy degeneracies seen at @xmath149 and singles out a unique lowest - energy state .
for small @xmath150 , the latter turns out to be the superposition of laughlin states for each individual component . at @xmath151
it becomes an incompressible state with a single laughlin quasi - particle in each component .
finally , at @xmath152 both components are bose - condensed in the lowest landau level .
modest interspecies - interaction strengths ( @xmath153 in panel b and @xmath154 in panel c ) cause avoided crossings but preserve the incompressible nature of the states seen in panel a. a stronger interspecies interaction ( @xmath155 in panel d ) washes out that picture completely.,width=453,height=340 ]    in figure  [ fig : energy3 ] , the interplay between interactions and confinement is elucidated .
panels  a  d show the evolution of low - lying few - particle eigenstates as the confinement strength is varied for situations with different magnitude of interaction strength between opposite - spin particles .
panel  a corresponds to the case with @xmath138 . due to the occurrence of level crossings ,
the character of the lowest - energy ( ground ) state is found to be different for regimes associated with weak , intermediate , and strong confinement .
analogous behavior has been discussed previously for ordinary ( spinless ) few - boson fractional qh systems  @xcite . in our case depicted in fig .
[ fig : energy3]a , the ground state in the weak - confinement regime corresponds to a superposition of three - particle laughlin states for filling factor @xmath143 in the individual pseudospin components .
after the first level crossing , each component turns out to be in the laughlin - quasiparticle state  @xcite and , after another level crossing , each spin component has its three particles occupying the lowest state defined by the parabolic confinement potential .
( our description of the ground states found in the three different regimes is supported by the analysis of real - space density and angular - momentum distribution functions .
see the following subsection for details . )
switching on interactions between opposite - spin particles turns crossings into anti - crossings
. figures .
[ fig : energy3]b and [ fig : energy3]c depict situations where interactions between same - spin particles are still dominant .
the existence of anticrossings enables smooth transitions between the different ground states that would not be possible in the case of simple crossings as seen , e.g. , in panel  a. independent tuning of interactions between opposite - spin particles can therefore be used to enable engineering of quantum many - particle states in ways not anticipated in previous work  @xcite .
strong interactions between opposite - spin particles are again seen to fundamentally alter the character of the system s ground and excited states . in fig .
[ fig : energy3]d , the strengths of interactions between same - spin and opposite - spin particles are equal .
the variation of few - particle states as a function of confinement strength is seen to be almost uniform , again pointing to the loss of distinctiveness for few - particle states in the presence of inter - species interactions .
furthermore , energy differences between low - lying states are much reduced as compared to the situation depicted in panel a of the same figure , which is a reflection of the unusual distribution of energy eigenvalues found for the interacting opposite - spin two - particle system .
the one - particle density profiles in coordinate space and in angular - momentum space are useful quantities to enable greater understanding of the properties of specific many - body quantum states  @xcite . in the basis of lowest - landau - level states from the two spin components , the single - particle density matrix of a many - particle state @xmath156 has matrix elements @xmath157 in terms of this quantity
, we can define the angular - momentum distribution for each spin component , @xmath158 and also the spin - resolved single - particle density profile in real space , @xmath159^\ast \phi^{(\sigma)}_{0 , m'}({{\ensuremath{\bm{\mathrm{r } } } } } ) \quad .\ ] ] in the following , we focus on the properties of the lowest - energy ( ground ) state in the different regimes associated with small , intermediate , and strong confinement strength for the systems whose energy spectra are shown in fig .
[ fig : energy3 ] .
cross - sectional pseudo - spin @xmath160 density profiles of the few - particle ground state associated with the lowest - lying energy level shown in the corresponding panels a  d of fig .
[ fig : energy3 ] , aggregated as a function of the confinement - potential strength @xmath150 . in panel
a ( only particles with same spin interact ) , sharp transitions occur between the fqh ( laughlin ) state in the regime of small @xmath150 , a laughlin - quasiparticle - type state for intermediate @xmath150 , and the gaussian bose - einstein - condensed state at high @xmath150 . for moderate interaction strength between opposite - spin components ( repulsive in panel b , attractive in panel c ) ,
transitions become smooth crossovers associated with anticrossings in fig .
[ fig : energy3 ] .
stronger interactions strengths between the spin components significantly change the character of the few - particle state at small @xmath150 ( panel d).,title="fig : " ]   cross - sectional pseudo - spin @xmath160 density profiles of the few - particle ground state associated with the lowest - lying energy level shown in the corresponding panels a  d of fig .
[ fig : energy3 ] , aggregated as a function of the confinement - potential strength @xmath150 . in panel
a ( only particles with same spin interact ) , sharp transitions occur between the fqh ( laughlin ) state in the regime of small @xmath150 , a laughlin - quasiparticle - type state for intermediate @xmath150 , and the gaussian bose - einstein - condensed state at high @xmath150 .
for moderate interaction strength between opposite - spin components ( repulsive in panel b , attractive in panel c ) , transitions become smooth crossovers associated with anticrossings in fig .  [
fig : energy3 ] .
stronger interactions strengths between the spin components significantly change the character of the few - particle state at small @xmath150 ( panel d).,title="fig : " ] +   cross - sectional pseudo - spin @xmath160 density profiles of the few - particle ground state associated with the lowest - lying energy level shown in the corresponding panels a  d of fig .
[ fig : energy3 ] , aggregated as a function of the confinement - potential strength @xmath150 . in panel
a ( only particles with same spin interact ) , sharp transitions occur between the fqh ( laughlin ) state in the regime of small @xmath150 , a laughlin - quasiparticle - type state for intermediate @xmath150 , and the gaussian bose - einstein - condensed state at high @xmath150 . for moderate interaction strength between opposite - spin components ( repulsive in panel b , attractive in panel c ) , transitions become smooth crossovers associated with anticrossings in fig .  [
fig : energy3 ] .
stronger interactions strengths between the spin components significantly change the character of the few - particle state at small @xmath150 ( panel d).,title="fig : " ]   cross - sectional pseudo - spin @xmath160 density profiles of the few - particle ground state associated with the lowest - lying energy level shown in the corresponding panels a  d of fig .
[ fig : energy3 ] , aggregated as a function of the confinement - potential strength @xmath150 . in panel
a ( only particles with same spin interact ) , sharp transitions occur between the fqh ( laughlin ) state in the regime of small @xmath150 , a laughlin - quasiparticle - type state for intermediate @xmath150 , and the gaussian bose - einstein - condensed state at high @xmath150 . for moderate interaction strength between opposite - spin components ( repulsive in panel b , attractive in panel c ) , transitions become smooth crossovers associated with anticrossings in fig .
[ fig : energy3 ] .
stronger interactions strengths between the spin components significantly change the character of the few - particle state at small @xmath150 ( panel d).,title="fig : " ]    figure  [ fig : densities3 ] shows the real - space profile of @xmath161 across a diameter of the disk - shaped three - particle systems associated with the ground - state levels shown in fig .  [
fig : energy3 ] . in the absence of interactions between opposite - spin particles ,
each component realizes correlated few - particle states of the type that have been found in previous work  @xcite .
the few - particle filling - factor-1/2 fqh state is the ground state for a weak confinement potential .
increasing the trapping - potential strength favors more compact correlated states , hence , at a critical value of @xmath150 , a transition occurs to a three - particle version of the laughlin - quasiparticle excited state . at
even higher @xmath150 , the system transitions to the gaussian bose - einstein - condensate state .
the sharpness of the transitions reflects the existence of level crossings in fig .
[ fig : energy3]a .
practically , simple variation of @xmath150 would not lead to any such transitions because there is no mechanism for the system to switch between different many - particle states . to make adiabatic passage between different many - particle states possible
, some symmetry of the system needs to broken , and previous work has proposed scenarios for achieving this in the single - component case  @xcite .
one - particle angular - momentum distribution for pseudo - spin @xmath160 particles for the ground states of systems whose energy spectra are shown in fig .  [
fig : energy3 ] .
compare also with the real - space density profiles shown in fig .
[ fig : densities3 ] .
note that the single - particle angular momentum cut - off at @xmath162 defines the sample size for vanishing @xmath150 in situations where opposite - spin particles interact ( panels b  d).,title="fig : " ]   one - particle angular - momentum distribution for pseudo - spin @xmath160 particles for the ground states of systems whose energy spectra are shown in fig .
[ fig : energy3 ] . compare also with the real - space density profiles shown in fig .
[ fig : densities3 ] .
note that the single - particle angular momentum cut - off at @xmath162 defines the sample size for vanishing @xmath150 in situations where opposite - spin particles interact ( panels b  d).,title="fig : " ] +   one - particle angular - momentum distribution for pseudo - spin @xmath160 particles for the ground states of systems whose energy spectra are shown in fig .  [ fig : energy3 ] .
compare also with the real - space density profiles shown in fig .
[ fig : densities3 ] .
note that the single - particle angular momentum cut - off at @xmath162 defines the sample size for vanishing @xmath150 in situations where opposite - spin particles interact ( panels b  d).,title="fig : " ]   one - particle angular - momentum distribution for pseudo - spin @xmath160 particles for the ground states of systems whose energy spectra are shown in fig .
[ fig : energy3 ] .
compare also with the real - space density profiles shown in fig .
[ fig : densities3 ] .
note that the single - particle angular momentum cut - off at @xmath162 defines the sample size for vanishing @xmath150 in situations where opposite - spin particles interact ( panels b  d).,title="fig : " ]    for our system of interest , an additional possibility arises from the ability to tune the interaction strength between the two spin components . as seen in panels b and c of fig .
[ fig : energy3 ] , a moderate value of @xmath163 turns the crossings occurring in panel a into anti - crossings , thus , different many - particle states are now adiabatically connected .
concomitantly , there is a continuous evolution of the spin - resolved one - particle density profile as a function of the confinement strength seen in figs .
[ fig : densities3]b and [ fig : densities3]c . when interactions between same - spin and opposite - spin particles have the same magnitude , the density profile changes significantly ( see fig .  [
fig : densities3]d ) , which indicates that the character of many - particle ground states is very different from a fractional - qsh state .
investigation of the one - particle angular - momentum - state distribution for the few - particle ground states discussed so far further solidifies our conclusions .
see figure  [ fig : momentum3 ] . in the absence of interactions between opposite - spin particles , the characteristic distributions for few - particle versions of the laughlin and laughlin - quasiparticle states emerge at low and intermediate values of @xmath150 . in the limit of strong trapping potential ,
the system condenses into the @xmath164 state .
switching on moderate repulsive ( attractive ) interaction strength between opposite - spin particles smoothens the transitions and also shifts the critical values of @xmath150 to larger ( smaller ) values .
when interactions among same - spin and between opposite - spin particles have equal magnitude , the one - particle momentum distribution of the ground state differs markedly from that associated with a fractional - qh state .
it appears that strong inter - component interactions favor a state with increased occupation of high - angular - momentum states , spreading out the particles more evenly across the accessible sample size and leading to an accumulation at the system s boundary . in the limit of vanishingly small trapping - potential strength @xmath150 ,
the latter is defined by the cut - off for single - particle angular momentum applied in our calculations .
we have considered the interplay of landau quantization and spin - dependent interactions in systems where particles with same spin feel the same strong magnetic field whereas particles with opposite spin are subject to magnetic fields with the same magnitude but opposite direction .
it has been expected  @xcite that such systems exhibit the fractional qsh effect , but we find that interactions between particles with opposite spin weaken or destroy features associated with fractional - qsh physics .
similar behavior has been seen in numerical studies of lattice realizations of fractional - qsh systems  @xcite .
we have elucidated how behavior that is very different from ordinary two - component fractional - qh systems is rooted in the drastically different spectral properties of two - particle interactions for particles feeling the same versus opposite magnetic - field directions .
thus any feasible route towards realizing the fractional qsh effect using a spin - dependent uniform magnetic field  @xcite should strive to eliminate interactions between the opposite - spin components . if the opposite - spin interaction strength is weak , adiabatic passage between different correlated many - particle states is facilitated by adjusting the strength of a trapping potential .
our conclusions are supported by numerically obtained real - space - density profiles and angular - momentum - state occupation distributions for few - particle systems .
the latter could also be utilized as blueprints for classifying images of correlated ultra - cold atom states .
part of the motivation for this project came about from stimulating conversations that one of us ( uz ) had with j.  j.  heremans and r.  winkler at the 2011 gordon godfrey workshop on spins and strong correlations ( sydney , australia , 24  28 october 2011 ) .
we would also like to thank m.  fleischhauer and a.  h.  macdonald for useful discussions .
of was supported by the marsden fund council from government funding ( contract no .
mau1205 ) , administered by the royal society of new zealand .
10 url # 1#1urlprefix[2][]#2 bloch i , dalibard j and nascimbene s 2012 _ nat .
phys . _ * 8 * 267276


understanding of fully developed turbulence is a complex , rich , and challenging problem . among the most important features of such behavior are cascades of energy and intermittency .
the former one brings energy from large scales , responsible for a creation of turbulence , toward smaller scales , in which viscosity begins to play a major role and dissipation effects dominate , whereas the latter one , being itself an irregular alternation of phases of some dynamics , in fact means that very rare configurations of a system give the most significant contributions to statistical distributions . for turbulence
this fact is revealed in anomalous scaling , which implies singular behavior of various statistical quantities as functions of the integral turbulence scales  @xcite .    another very interesting problem is a turbulent advection of an impurity field .
both natural experiments and numerical simulations suggest that deviations from the classical kolmogorov theory are even more strongly pronounced for passively advected fields than for the velocity field itself  @xcite .
a turbulent environment in such models may be introduced by some `` synthetic '' velocity field with prescribed statistics or by stochastic navier - stokes equation  @xcite .
models of the former type are more tractable from a mathematical point of view , whereas the latter models bear a closer resemblance to the real world .
a very powerful theoretical framework , which has worked very well in statistical physics in order to investigate different models of critical phenomena , characterized by scaling behavior , is the field theoretic renormalization group ( rg ) and the operator product expansion ( ope ) ; it has become an indispensable tool when dealing with systems where no clear separation of scales occurs , see the monographs  @xcite . in the rg+ope scenario ,
an anomalous scaling emerges as a consequence of the existence of some composite fields ( `` composite operators '' in the quantum - field terminology ) in the model with negative scaling dimensions  @xcite .
the main advantage of this approach with respect to the turbulence is that it is based on a microscopic model and allows one to construct a systematic perturbation expansion for the anomalous exponents , similar to the famous @xmath4-expansion in theory of critical phenomena .
the procedure has been successful in verifying kolmogorov scaling and provides an efficient theoretical tool for calculation of universal quantities  @xcite . in a number of papers
the rg+ope approach has been applied to the case of passive advection by kraichnan s ensemble ( velocity field is taken isotropic , gaussian , not correlated in time , having a power - like correlation function , and fluid assumed to be incompressible ) , see  @xcite , and by its numerous generalizations : large - scale anisotropy , helicity , compressibility , finite correlation time , non - gaussianity , and a more general form of the nonlinearity  @xcite .
a similar approach can be generalized as well to the case of non - gaussian velocity field , governed by the stochastic navier - stokes equation   both for the scaling behavior of velocity field itself , and for an advection of passive fields by this ensemble  @xcite .    up to now a majority of works on fully developed turbulence
have been concerned with an incompressible fluid .
nevertheless , several results for compressible fluids have also been obtained  @xcite .
all of them hints at large influence of compressibility both to velocity field itself and passively advected quantities .
in particular , the transition from a turbulent to a certain purely chaotic state may occur at large degrees of compressibility  @xcite .
in other papers corrections in the mach number to the incompressible scaling regime were studied  @xcite .
thus , obtained corrections become arbitrarily large and destroy the incompressible scaling regime for fixed mach numbers and large distances , what can be explained by an existence of crossover to another , yet unknown regime .
the works  @xcite were devoted to a compressible fluid itself .
the results are rather controversial , particularly the model , considered in  @xcite , appears to be in fact unrenormalizable . from a general point of view further investigations of compressibility are therefore called for .    in this paper
we present an application of the field theoretic renormalization group onto the scaling regimes of a compressible fluid , whose behavior is governed by a proper generalization of the stochastic navier stokes equation  @xcite .
a resulting stationary scaling regime in this approach is associated with an infrared ( ir ) attractive fixed point of the corresponding multiplicatively renormalizable field theoretic model . in a similar fashion as has been done for an incompressible case  @xcite , the double expansion scheme is employed . here the formal expansion parameters are @xmath2 , which describes the scaling behavior of a random force , and @xmath5 , i.e. , a deviation from the space dimension @xmath0 .    following a procedure of previous works  @xcite , first the stochastic navier - stokes equation is discussed . after establishing an existence of necessary fixed points of navier - stokes equation ,
the advection of scalar fields is explored . since rg functions of the parameters , entering into navier - stokes equation , do not depend on the parameters , connected with advection - diffusion equation ,
this is a possible and probably the easiest way .
the paper is organized as follows .    in sec .
[ sec : model ] a detailed description of the stochastic navier - stokes equation for a compressible fluid is given .
[ sec : qft ] is devoted to field theoretic formulation of the model and the corresponding diagrammatic technique .
in particular , possible types of divergent green functions at @xmath1 and @xmath0 and the necessity of introducing a new coupling constant at @xmath0 are discussed . in sec .  [
sec : renorm ] the renormalizability of the model is established and one - loop explicit expressions for the renormalization constants are derived . much attention is paid here to technical aspects of performed calculations . in sec .
[ sec : rc ] the rg functions ( anomalous dimensions and @xmath6 functions ) are derived and examined .
ir asymptotic behavior , obtained as a result of solution of rg equations , is discussed .
it is shown that , depending on two exponents @xmath2 and @xmath4 , the rg equations possess an ir attractive fixed point , which implies existence of a scaling regime in the inertial range .
the corresponding scaling dimensions of all fields and parameters of the model are presented .    in sec .
[ sec : ops1 ] advection of a passive scalar ( density ) field by compressible velocity field , which obeys navier - stokes equation , is analyzed .
field theoretic formulation of the full model is presented .
it is shown that the full model is multiplicatively renormalizable ; the existence of a scaling regime in the ir range is established .
the renormalization of composite operators is carried out .
an inertial - range behavior of various correlation functions is studied by means of the ope .
it is shown that leading terms of the inertial - range behavior are determined by the contributions of the operators built solely of the scalar fields itself . as a result ,
the ir behavior of the pair correlation functions of the composite operators is power - like with negative critical dimensions   the situation , called anomalous scaling , occurs .    sec .
[ sec : conc ] is reserved for conclusion .
the navier - stokes equation for a viscid compressible fluid can be written in the following form  @xcite : @xmath7 v_{k } + \mu_0 \partial_{i}\partial_{k } v_{k } - \partial_{i } p + \eta_{i } , \nonumber \\ { } \label{ns}\end{aligned}\ ] ] where the differential operator on the right hand side @xmath8 is the lagrangian ( convective ) derivative , @xmath9 is a fluid density field , @xmath10 is a velocity field , @xmath11 , @xmath12 , @xmath13 is the laplace operator , @xmath14 is a pressure field , and @xmath15 is a density of an external force per unit volume .
the fields @xmath10 , @xmath16 , @xmath9 and @xmath14 depend on @xmath17 with @xmath18 , where @xmath19 is a dimensionality of space .
the constants @xmath20 and @xmath21 are two independent molecular viscosity coefficients  @xcite ; in   we have explicitly separated the transverse and longitudinal components of the viscous term .
summations over repeated vector indices are always implied in this work .    to get a closed system of equations , the model  ( [ ns ] )
must be augmented by two additional equations , namely a continuity equation and an equation of state between deviations @xmath22 and @xmath23 from the equilibrium values .
they explicitly read @xmath24 and @xmath25    in order to obtain the renormalizable field , a theoretic model expression  ( [ ns ] ) has to be divided by @xmath9 , and fluctuations in viscous terms have to be neglected  @xcite .
further , by using expressions   and   the problem can be recast in the form of two coupled equations : @xmath26 v_{k}\ ! + \
! \mu_0 \partial_{i}\partial_{k } v_{k } -\ !
\partial_{i } \phi\ ! + \ !
f_{i } , \label{anu } \\
\nabla_{t } \phi & = & -c_{0}^{2 } \partial_{i}v_{i}. \label{anu1}\end{aligned}\ ] ] here a new scalar field @xmath27 is related to the density fluctuations via the relation @xmath28 .
a parameter @xmath29 is an adiabatic speed of sound , @xmath30 denotes the mean value of @xmath9 , and @xmath31 is a density of the external force per unit mass .    in the stochastic formulation of the problem the turbulence
is modeled by an external force   it is assumed to be a random variable , which mimics the input of energy into the system from the outer large scale @xmath32 .
its precise form is believed to be unimportant and is usually considered to be a random gaussian variable with zero mean and prescribed correlation function  @xcite . for the use of the standard rg technique this correlator must exhibit a power law asymptotic behavior at large wave numbers  @xcite . in the case of compressible fluid
it should be naturally augmented with a non - transverse component , hence the simplest way is to choose it in the form  @xcite @xmath33 where the argument is given by @xmath34 here , @xmath35 and @xmath36 are the transverse and longitudinal projectors , @xmath37 , amplitude @xmath38 is a free parameter , an amplitude @xmath39 is a coupling constant ( formal expansion parameter in the ordinary perturbation theory ) ; the relation @xmath40 sets in the typical ultraviolet ( uv ) momentum scale @xmath41 , which is a reciprocal of the dissipation length scale .
a parameter @xmath42 provides an infrared regularization ; its precise form is unessential and the sharp cut - off is the simplest choice for calculation purposes .
an exponent @xmath2 provides the uv regularization and therefore plays a role of a formally small expansion parameter  @xcite .
the most realistic ( physical ) value is obtained in the limit @xmath43 , when the function in  ( [ power ] ) can be interpreted as power - like representation of the dirac function @xmath44 : physically it corresponds to the idealized picture of the energy input from infinitely large scales .
the galilean invariance for the model ( [ ns ] ) is ensured when the function   is delta - correlated in time  @xcite .
according to the general theorem  @xcite , the stochastic problem   is equivalent to the field theoretic model with a doubled set of fields @xmath45 and de dominicis - janssen action functional , written in a compact form as    @xmath46v_k    + u_0 \nu_0 \partial_i \partial_k v_k - \partial_i \phi    \biggl\ } \nonumber \\   &    + \phi'[-\partial_t \phi + v_j\partial_j \phi + v_0 \nu_0 \partial^2 \phi - c_0 ^ 2 ( \partial_i v_i ) ] ,    \label{eq : action } \end{aligned}\ ] ]    where @xmath47 is the correlation function  ( [ power ] ) . here
a condensed notation has been employed , in which integrals over the spatial variable @xmath48 and the time variable @xmath49 , as well as summation over repeated indices , are implicitly assumed .
moreover , we have introduced a new dimensionless parameter @xmath50 and introduced a new term @xmath51 with another positive dimensionless parameter @xmath52 , which is needed to ensure multiplicative renormalizability .
the action  ( [ eq : action ] ) is amenable to the standard methods of the quantum field theory , such as feynman diagrammatic technique and renormalization group procedure .    in a standard approach , if quantum field methods are applied to the stochastic differential equations , the space dimension @xmath19 plays a passive role and an actual perturbative parameter is @xmath2 ; for more details see the monographs  @xcite .
our approach closely follows the analysis of the incompressible navier - stokes equation near space dimension @xmath53  @xcite . in this case
an additional divergence appear in the green function @xmath54 .
it can be absorbed by a suitable local counterterm @xmath55 , and a regular expansion in both @xmath2 and @xmath56 must be constructed . up to now the present model  ( [ eq : action ] ) has been investigated in fixed space dimension @xmath1 , for which the action  ( [ eq : action ] ) contains all terms that can be generated during the renormalization procedure  @xcite . however , from the dimensional analysis ( see below ) it follows that in @xmath0 an additional divergence appears in a similar fashion in the green function @xmath54 .
therefore , to keep the model renormalizable in @xmath0 , the kernel function @xmath57 in  ( [ ff ] ) has to be generalized into the function @xmath58 , where @xmath59 a new term on the right hand side with an additional coupling constant @xmath60 absorbs divergent contributions from @xmath54 .
in contrast to the two - dimensional incompressible case  @xcite no momentum dependence is needed here .
the field theoretic formulation  ( [ eq : action ] ) means that various correlation and response functions of the original stochastic problem are represented by functional averages over the full set of fields with functional weight @xmath61 , and in this sense they can be viewed as green functions of the field theoretic model @xcite .
the perturbation theory of the model can be expressed according to the usual feynman diagrammatic expansion  @xcite .
bare propagators are read off from the inverse matrix of the gaussian ( free ) part of the action functional , while a nonlinear part of the differential equation defines the interaction vertices . a graphical representation of the propagator functions
is depicted in fig .
[ fig : prop_vertex ] , and of the vertices   in fig .  [ fig:2 ] .
from   it follows that the feynman diagrammatic technique for this model contains two vertices , @xmath62 and @xmath63 .
the propagator functions in the frequency - momentum representation read @xmath64 where the vector indices of the fields and the projectors are omitted . for convenience
following abbreviations have been used @xmath65 and @xmath66+g_{20}\nu_0 ^ 3 .
\label{energies}\end{aligned}\ ] ]    in the limit @xmath67 propagator functions @xmath68 and @xmath69 become purely transverse , and all mixed propagators except @xmath70 vanish .
moreover , the scalar fields @xmath71 and @xmath72 decouple from the velocity fields @xmath73 and @xmath74   it is impossible to construct a diagram with only velocity fields @xmath73 and @xmath75 as external lines , containing internal lines with fields @xmath71 or @xmath72 . thus , as it should be from physical point of view , the well - known feynman rules for the incompressible fluid  @xcite are obtained .    [ cols="^ " , ]      canonical dimensions of the new fields and parameters of the full model  ( [ fact ] ) are given in table  [ table1 ] , where we have introduced a new dimensionless parameter @xmath76 with @xmath20 from  ( [ ns ] ) .
the expression  ( [ index ] ) for the formal index of uv divergence remains valid , but the summation should now run over the full set of six fields @xmath77 . the rules ( 1 )   ( 6 ) from sec .
[ sec : canon ] should be generalized and augmented as follows :    1 .
all the 1-irreducible green functions without the response field @xmath78 vanish identically and require no counterterms .
2 .   using integration by parts the derivative at the vertex @xmath79 can be moved onto the field @xmath78 , therefore expression   should be modified as @xmath80 since the field @xmath78 can enter the counterterms only in the form of spatial derivatives , the counterterm @xmath81 to the 1-irreducible green function @xmath82 with @xmath83 , @xmath84 is forbidden .
also this requires that the counterterms to the 1-irreducible function @xmath85 with @xmath86 , @xmath87 necessarily reduce to the form @xmath88 .
galilean symmetry wherein allows them enter the counterterms only in the form of invariant combination @xmath89 .
hence , they are also forbidden .
3 .   as consequence of the linearity of the original stochastic equation  ( [ density1 ] ) with respect to the field @xmath90 one obtains that for any 1-irreducible function @xmath91 , where @xmath92 is the total number of bare propagators @xmath93 entering the diagram .
this fact is very important for the renormalizability of the model : without the restriction @xmath94 , we the infinity number of superficially divergent 1-irreducible functions @xmath95 would proliferate , and hence the lack of renormalizability would follow .    from these rules
we finally conclude that superficial divergences can be present only in 1-irreducible green function @xmath82 with the only counterterm @xmath96 .
it is naturally reproduced as multiplicative renormalization of the diffusion coefficient , @xmath97 .
no renormalization of the fields @xmath78 and @xmath90 is needed : @xmath98 .
the renormalized analog of the action functional  ( [ fact ] ) has the form    @xmath99    where @xmath100 is the action  ( [ raction ] ) ,    @xmath101    @xmath102 here stands for the covariance of stochastic force  .      the one - loop approximation for the 1-irreducible response function @xmath82 can be formally written as @xmath103{g2_1.eps } } ,   \label{dyson}\end{aligned}\ ] ] where , as earlier in expressions
, @xmath104 represents an external momentum ; single solid line denotes the bare propagator @xmath105 from ( [ lines ] ) , double solid line with a slash denotes the bare propagator @xmath106 from ( [ lines3 ] ) , the slashed end corresponds to the field @xmath78 .
the dots with three attached fields @xmath78 , @xmath90 and @xmath107 denote the vertex  ( [ vertex1 ] ) .
the constant @xmath108 is to be found from the requirement of uv finiteness of 1-irreducible green function @xmath82 . like for the original navier - stokes model ,
the divergent part of the considered feynman diagram is independent on @xmath109 and therefore can be calculated directly at @xmath110 ; see the discussion in sec .
[ sec : c ] .
the analytical expression for the diagram @xmath111 , entering into expression  , is @xmath112,\end{aligned}\ ] ] where @xmath113 and @xmath114 are two vertices of the type  ; scalar coefficients @xmath115 and @xmath116 of the propagator @xmath117 are defined in  .    since in the leading - order approximation the renormalization constant @xmath118 in the bare term of  ( [ dyson ] )
is taken only in the first order in coupling constants @xmath119 and @xmath120 , during the actual calculation all other renormalization constants in the diagram @xmath111 , entering , for example , into functions @xmath115 and @xmath116 , should be replaced with unities .
integration over the frequency gives @xmath121 } ; \nonumber \\ \int\frac{\drm\omega}{2\pi } \frac{b(k ) } { -i\omega+w_0\nu ( { \bf p}+{\bf k})^{2 } } & = \frac{1}{2\nu^2uk^2[uk^2+w_0({\bf p+k})^2]}. \label{tmp}\end{aligned}\ ] ] divide expression   into two parts : @xmath122 } ; \nonumber \\ \check{i}_2&=-\int\frac{\drm^d{\bf k}}{(2\pi)^{d}}p_a(p+k)_c \times \frac{q_{ac}({\bf
k})}{2\nu^2uk^2[uk^2+w_0({\bf p+k})^2]}.\end{aligned}\ ] ] we are interesting in the term , proportional to @xmath123 , therefore in calculation of @xmath124 one may immediately put @xmath125 and , using expressions   and   and notation @xmath126 , gets @xmath127 in calculation of @xmath128 the second expression of   has to be expanded up to the linear term in @xmath104 ; in the end we get @xmath129 finally , collecting @xmath124 and @xmath128 yields    @xmath130.\end{aligned}\ ] ]    then the renormalization constant @xmath118 [ see  ] should be chosen as @xmath131\frac{g_1}{y } \nonumber \\ - \frac{1}{2dw } \left [ \frac{d-1}{1+w } + \frac{u - w}{u(u+w)^{2 } } \right]\frac{g_2}{\varepsilon } , \label{zk}\end{aligned}\ ] ] while the corresponding anomalous dimension is @xmath132g_1 \nonumber \\
+ \frac{1}{2dw } \left [ \frac{d-1}{1+w } + \frac{u - w}{u(u+w)^{2 } } \right]g_2 , \label{gk}\end{aligned}\ ] ] with the possible corrections coming from higher orders terms .
the function @xmath133 for the new dimensionless parameter @xmath134 has the form @xmath135 see sec .
[ sec : rgrg ] .
now the coordinates @xmath136 of fixed points fp@xmath137 and fp@xmath138 [ see   and  ] are substituted into the equation @xmath139 at @xmath0 . to this end we rewrite expression for @xmath140 at @xmath141 : @xmath142 . \label{u1}\end{aligned}\ ] ] from equation   it is clear that the only positive solution for both fp@xmath137 and fp@xmath138 is @xmath143    functions  ( [ betagw ] ) do not depend on @xmath134 , therefore new eigenvalue of the matrix  ( [ omega ] ) , corresponding to this parameter , coincides with the diagonal element @xmath144 at point @xmath145 : @xmath146}>0 \quad\text{for { fp$ { \textrm{iii}}$}}. \label{wir}\end{aligned}\ ] ]    from the inequalities   it follows that fixed points with the coordinates   and   and @xmath147 are ir attractive in the full space of couplings @xmath148 and governs the ir asymptotic behavior of the full - scale model  ( [ fact ] ) .
the critical dimensions of the passive fields @xmath90 and @xmath78 are obtained from the table  [ table1 ] and the expression  ( [ krit ] ) for @xmath149 .
for fixed point fp@xmath138 they are the same as in the case @xmath1 , namely @xmath150 for another point fp@xmath137 they are @xmath151      in the following , the central role will be played by composite fields ( `` composite operators '' ) built solely of the basic fields @xmath90 : @xmath152 in general , a local composite operator is a polynomial constructed from the primary fields @xmath153 and their finite - order derivatives at a single space - time point @xmath17 . due to coincidence of the field arguments , new uv divergences arise in the green functions with such objects  @xcite .    the total canonical dimension of arbitrary 1-irreducible green function @xmath154 that includes one composite operator @xmath155 and arbitrary number of primary fields @xmath156 ( the formal index of uv divergence ) is given by the relation @xmath157 where @xmath158 are the numbers of the fields entering into @xmath159 , @xmath160 are their total canonical dimensions , @xmath161 is the canonical dimension of the operator , and the summation over all types of the fields is implied .    during renormalization
procedure operators can mix with each other , @xmath162 and @xmath163 is the renormalization matrix .
but in the simplest case of operators   matrix @xmath163 is diagonal , i.e. , @xmath164 . in particular
this means that the critical dimension of the operator is given by the expression  ( [ krit ] ) .
superficial uv divergences , whose removal requires counterterms , can be present only in those functions @xmath159 for which the index of divergence @xmath165 is a non - negative integer .
for the operators of the form   one has @xmath166 . due to the linearity of our model in @xmath90 , the number of the fields @xmath90 in any 1-irreducible function with the operator
@xmath167 can not exceed their number in the operator itself . from the analysis of expression  ( [ d - comp ] ) it therefore follows that the superficial divergence can only be present in the 1-irreducible function with @xmath168 and @xmath169 for all other types of the fields @xmath156 .
for this function @xmath170 and the corresponding counterterm has the form @xmath171 ; hence , operators   are multiplicatively renormalizable , @xmath172 .    introduce @xmath173 : the @xmath174 term of the expansion in @xmath175 of the generating functional of the 1-irreducible green functions with one composite operator @xmath167 and any number of fields @xmath90 :    @xmath176    renormalization constants @xmath177 are determined by the requirement that the 1-irreducible functions   be uv finite in renormalized theory .
the one - loop approximation for the 1-irreducible function @xmath178 can be formally written as @xmath179{comp_density.eps } } } \label{gamma2}\end{aligned}\ ] ] the first term is the tree ( loop - less ) approximation , double solid lines with a slash denotes the propagators @xmath180 , single solid line corresponds to the propagator @xmath181 , @xmath182 is the symmetry coefficient of the given graph , and the dot with two attached lines in the top of the diagram denotes the operator vertex , i.e. , the variational derivative @xmath183 the contribution of a specific diagram into the functional   for any composite operator @xmath155 is represented in the form @xmath184 where @xmath185 is the vertex factor  , @xmath186 is the diagram itself , and the product @xmath187 corresponds to external tails .
the divergence of the considering graph is logarithmic , hence one should set all the external frequencies and momenta equal to zero .
therefore the analytical expression of the diagram is given by @xmath188,\end{aligned}\ ] ] where @xmath189 and @xmath190 are two vertices  ; scalar coefficients @xmath115 and @xmath116 of the propagator @xmath117 are defined in   with the replacement of original bare parameters to their renormalized counterparts . as @xmath191 , only second term in   gives non - vanishing contribution .
integration over the frequency gives @xmath192 contracting tensor indices , using  , and collecting all the factors one obtains that expression   takes form @xmath193 \biggl\ } , \label{gamma22}\end{aligned}\ ] ] in notation @xmath194 and up to a finite part and higher - order corrections .
the renormalization constants @xmath195 are found from the requirement that the renormalized analog @xmath196 of the function ( [ gamma1 ] ) be uv finite in terms of renormalized parameters .
in contrast to expressions   , in this case renormalization constants @xmath177 belong not to some parameters , but to the green functions itself .
this means that using loop expansion   one finds in fact not renormalization constants @xmath177 , but inversed one @xmath197 .
taking into account minus sign in the exponent , from   it follows that in the ms scheme renormalization constants take form @xmath198 the corresponding anomalous dimensions are @xmath199 with higher - order corrections in @xmath119 and @xmath120 .
the critical dimensions of the operators @xmath174 from the expression  ( [ krit ] ) are readily derived @xmath200= n\delta_{\theta } + \gamma_{n}^{*}. \label{krop}\ ] ] substituting of fixed - point values fp@xmath137 and fp@xmath138 into  ( [ gaman ] ) finally gives that for fixed point fp@xmath137 critical dimensions are @xmath201&=-n+\frac{n\varepsilon}{4}-\frac{n(n-1)}{3}\varepsilon ; \label{krop1}\end{aligned}\ ] ] for fixed point fp@xmath138 they differs from the case @xmath1 and are equal to @xmath202&=-n+\frac{ny}{6}-\frac{2n(n-1)}{3}\frac{\alpha y(y-\varepsilon)}{y(2+\alpha)-3\varepsilon}.\end{aligned}\ ] ] both expressions   and   suppose higher - order corrections in @xmath2 and @xmath4 . for both cases fp@xmath137 and fp@xmath138 the dimensions are negative , i.e. , `` dangerous '' in the sense of operator product expansion  @xcite , and decrease as @xmath203 grows .
the measurable quantities and , therefore , the objects of interest are equal - time pair correlation functions of two ( uv finite ) renormalized local composite operators @xmath204 .
from the dimensionality considerations ( see table  [ table1 ] ) it follows that @xmath205 where @xmath206 and @xmath161 are frequency and total canonical dimensions of the correlation function , @xmath207 , and @xmath208 is a function of dimensionless variables .
if the correlation function  ( [ parr0 ] ) is multiplicatively renormalizable , in the ir region it fulfills the differential equation  , which describes the ir scaling behavior .
that is , the behavior of the function @xmath208 for @xmath209 is determined by the ir attractive fixed points fp@xmath137 and fp@xmath138 of the rg equation .
the solution of this equation leads to the following asymptotic expression : @xmath210 .
\label{parr}\ ] ] here @xmath211 is the critical dimension of the correlation function , given by simple sum of the dimensions of the operators ; @xmath212 is unknown scaling function with completely ( both canonically and critically ) dimensionless arguments , and @xmath213 is invariant speed of sound .    for the composite operator @xmath214 ,
expression  ( [ parr ] ) yields @xmath215 , \label{parr2}\ ] ] where critical dimensions @xmath216 for two scaling regimes are given by expressions  ( [ krop1 ] ) and  ( [ krop2 ] ) .    representation   holds for @xmath217 and any fixed value of @xmath218 .
the inertial - convective range @xmath219 corresponds to the additional condition @xmath220 .
the behavior of the functions @xmath212 at @xmath221 can be studied by means of the operator product expansion ; see  @xcite . according to the ope
, the equal - time product @xmath222 of two renormalized operators for @xmath223 and @xmath224 has the representation    @xmath225 f(t,{\bf x } ) , \label{ope}\ ] ]    where @xmath226 are numerical coefficient functions analytical in @xmath218 and @xmath213 and @xmath155 are all possible renormalized local composite operators allowed by the symmetry .
the correlation function  ( [ parr ] ) is obtained by averaging  ( [ ope ] ) with the weight @xmath227 , where @xmath228 is the renormalized action functional  ( [ fact ] ) .
mean values @xmath229 appear on the right hand side .
their asymptotic behavior at small @xmath230 is found from the corresponding rg equations and has the form    @xmath231 , \label{mean}\ ] ]    with another set of scaling functions @xmath232 .
since the diagrams of the perturbation theory have finite limits both for @xmath233 and @xmath234 , we may assume that the functions @xmath235 are restricted for all values of @xmath236 and can be estimated by some constants .
moreover , for invariant variable @xmath213 ir asymptotic behavior together with requirement of its dimensionless gives @xmath237 where @xmath236 is renormalized speed fo sound .
this means that @xmath238 .
taking into account  , for non - local scaling regime fp@xmath138 one obtains that for @xmath239 ( i.e. , including the most realistic case @xmath43 ) the argument @xmath240 becomes small for fixed @xmath236 and @xmath241 , and the function @xmath232 can be replaced by its finite limit value @xmath242 . for local scaling regime fp@xmath137 from   it follows that at @xmath243 function @xmath232 can be replaced by its finite limit value @xmath244 . from these two remarks
we finally conclude that in the ir range for both local and non - local scaling regimes up to a different constants we can write @xmath245 combining rg representation  ( [ parr2 ] ) with information gained from ope  ( [ ope ] ) and expression  ( [ mean2 ] ) gives the desired asymptotic behavior of the scaling functions : @xmath246 \simeq \sum_{f } a_{f } [ mr , c(r)]\times ( mr)^{\delta_{f } } , \label{fin}\ ] ] where the summation runs over all galilean invariant scalar operators ( including operators with derivatives , etc . ) , with the coefficient functions @xmath247 analytical in their arguments .
the leading contribution into the sum   is given by the operator with the lowest ( minimal ) critical dimension ; the others can be considered as corrections .
the anomalous scaling ( i.e. , singular behavior as @xmath221 ) results from the contributions of the operators with negative critical dimensions . from   and
it is easily seen that for both scaling regimes all the operators @xmath174 have negative dimensions , and the spectrum of their dimensions is not restricted from below .
fortunately , due to the linearity of initial stochastic equation   in the field @xmath90 , the number of such fields on the right hand side of expression   can not exceed their number on the left hand side .
this means that for a given correlation function only a finite number of those operators can contribute to the ope . for correlation functions  ( [ parr2 ] )
these operators are those for which @xmath248 .
the leading term of the behavior as @xmath221 is given by the operator with the maximum possible @xmath249 and without any derivatives , so the final expression has the form @xmath250    the fact that the leading term in the ope is given by the operator from the same family with the summed exponent together with inequality @xmath251 can be interpreted as the statement that the correlations of the scalar field in the model   show multifractal behavior ; see  @xcite .
in this paper , which is a natural extension of  @xcite , the stochastic navier stokes equation for a _
compressible _ fluid has been studied using the field theoretic approach .
in contrast to previous work , the model is considered in special space dimension @xmath0 . in this case
additional uv divergence is observed in 1-irreducible green function @xmath252 .
this fact significantly affects rg analysis , regarding both technical aspects and results .
in particular , renormalization group technique with double expansion scheme is needed . in the one loop approximation , depending on the exponent @xmath2 and deviation from the dimension of @xmath48 space @xmath5 , the model possesses two stable nontrivial fixed points in the ir region ( i.e. , two possible nontrivial scaling regimes )   local one , called in text fp@xmath137 , and nonlocal , called fp@xmath138 .
this shows that the simple analysis around @xmath1 , which indicates existence of only one nontrivial fixed point , corresponding to non - local scaling regime , can be regard as incomplete in this case .
the crossover between local and non - local regimes occurs along the line @xmath253 , which is in accordance with  @xcite .
the new ( local ) regime , which takes place at @xmath0 , should by continuity move into @xmath1 at @xmath243 .
nevertheless , at physical values @xmath254 and @xmath255 new local regime fp@xmath137 is unstable , and non - local regime fp@xmath138 is realized , therefore the main conclusions of the previous papers  @xcite about existence and ir attraction of this fixed point remain true
.    the model of passive scalar advection of density field by this velocity ensemble has been also analyzed .
the full stochastic problem can be formulated as field theoretic model , which appear multiplicatively renormalizable .
the new parameter @xmath256 does not affect to rg functions of the navier - stokes equation itself , so the critical behavior of this model is also described by two fixed points , local one and non - local one .
the inertial range ( @xmath219 ) behavior of correlation functions has been studied by means of the ope technique .
existence of anomalous scaling , i.e. , singular power - like dependence on the integral scale @xmath32 , has been established .
the corresponding anomalous exponents have been identified with critical dimensions of certain composite operators and calculated in the leading one - loop approximation .    as in the case @xmath1 ,
results of the present model have some similarities and some differences with simplified models , which describe compressible turbulent velocity field as some gaussian statistical ensemble  @xcite .
in contrast to  @xcite ( gaussian ensemble with finite correlation time ) , the present model is manifestly galilean covariant , and this fact holds in all orders of the perturbation theory . as in the rapid - change ( zero correlation time ) models ,
critical dimensions are insensitive to the specific choice of the random force .
this is a consequence of the fact that propagator @xmath257 does not enter into the relevant green functions .    in principle
, the anomalous exponents could depend on some dimensionless parameters like @xmath258 , @xmath259 , or @xmath260 .
after the rg procedure , these parameters are replaced with the corresponding invariant variables .
existence of the unique ir attractive fixed point for each pair @xmath261 shows that in the ir range these invariant variables tend to their fixed - point values @xmath262 , @xmath263 , @xmath264 irrespective of the initial values @xmath258 , @xmath259 , and @xmath260 .
this is a consequence of the one - loop approximation ; in higher - loop level overlaps may occur in the picture of domains of ir stability .
in contrast to the case @xmath1 , at @xmath0 anomalous dimensions of composite operators   and   do not grow with @xmath38 without bound .
this is a consequence of eliminating of poles in @xmath4 near @xmath0 , which leads to a significant improvement of situation near @xmath1 .
the situation is similar to consideration of navier - stokes equation for incompressible fluid near @xmath53 , where renormalization procedure at @xmath53 is a reason for for the kolmogorov constant s better agreement with experimental results  @xcite .    as in the case
@xmath1 , at @xmath0 the coordinates of both nontrivial fixed points are valid for all @xmath265 , have finite limits at @xmath266 , but ceases to exist for purely potential forcing , i.e. , for @xmath267 . this fact means that under the one - loop approximation both fixed points fp@xmath137 and fp@xmath138 appears to disappear or lose its stability . corresponding scaling regimes consequently undergoes some qualitative changeover , which is possibly accompanied by phase transition to a purely chaotic state observed in  @xcite for a simplified model .    to clarify this problem
it is necessary to go beyond the one - loop approximation and to discuss the existence , stability and the dependence on @xmath38 of fixed points at least at the two - loop level , which seems to be a difficult technical task .
also it would be very interesting to investigate a scalar admixture in case of tracer field , or passively advected vector fields .
another very important task is to develop compressible navier - stokes equation near @xmath53 .
such analysis may show another types of ir behavior or another dependence of different parameters like @xmath38 , viscosity ratios , etc .
this work is left for the future and is partly in progress .
the authors are indebted to m.  yu .
nalimov , l.  ts .
adzhemyan , m.  hnati , j.  honkonen , and v.  kultty for discussions .
the work was supported by vega grant no .
1/0222/13 of the ministry of education , science , research and sport of the slovak republic , and by the russian foundation for basic research within the project 16 - 32 - 00086 .
n.  m.  g. acknowledges the support from the saint petersburg committee of science and high school .
n.  m.  g. and m.  m.  k. were also supported by the dmitry zimin s `` dynasty '' foundation .
k. gawdzki and a.  kupiainen , phys .
lett . * 75 * , 3834 ( 1995 ) ; d.  bernard , k.  gawdzki , and a.  kupiainen , phys .
e * 54 * , 2564 ( 1996 ) ; m.  chertkov , g.  falkovich , i.  kolokolov , and v.  lebedev , phys .
e * 52 * , 4924 ( 1995 ) ; m.  chertkov and g.  falkovich , phys . rev
. lett . * 76 * , 2706 ( 1996 ) ; a.  pumir , europhys . lett . * 34 * , 25 ( 1996 ) ; * 37 * , 529 ( 1997 ) ; phys . rev .
e * 57 * , 2914 ( 1998 ) .
l.  ts .
adzhemyan , n.  v. antonov , v.  a.  barinov , yu .  s. kabrits , and a.  n. vasilev , phys .
e * 63 * , 025303(r ) ( 2001 ) ; e * 64 * , 019901(e ) ( 2001 ) ; e * 64 * , 056306 ( 2001 ) .
n.  v.  antonov , phys .
e * 60 * , 6691 ( 1999 ) ; l.  ts .
adzhemyan , n.  v.  antonov , and j.  honkonen , phys .
e * 66 * , 036313 ( 2002 ) .
d.  ronis , phys .
a * 36 * , 3322 ( 1987 ) ; l.  ts .
adzhemyan , j.  honkonen , m.  v.  kompaniets , a.  n.  vasilev , phys .
e * 68 * , 055302(r ) ( 2003 ) ; l.  ts .
adzhemyan , m.  hnatich and j.  honkonen , eur .
j b * 73 * , 275 ( 2010 ) .            c.  de  dominicis and p.  c.  martin , phys .
a * 19 * , 419 ( 1979 ) ; p.  l.  sulem , j.  d.  fournier , and u.  frisch , lecture notes in physics , * 104 * , 321 ( 1979 ) ; j.  d.  fournier and u.  frisch , phys .
a * 28 * , 1000 ( 1983 ) ; l.  ts .
adzhemyan , a.  n.  vasilev , and yu .  m.  pismak , theor .
* 57 * , 1131 ( 1983 ) .
adzhemyan , a.  n.  vasilev , and m.  hnatich , theor .
phys . * 74 * , 115 ( 1988 ) ; l.  ts .
adzhemyan , n.  v.  antonov , and t.  l.  kim , theor .
math . phys . * 100 * , 1086 ( 1994 ) ; n.  v.  antonov , s.  v.  borisenok , and v.  i.  girina , theor . math
. phys . * 106 * , 75 ( 1996 ) ; n.  v.  antonov and a.  n.  vasilev , theor .
phys . * 110 * , 97 ( 1997 ) .