# Dispersion properties of transverse waves in electrically polarized BECs

###### Abstract

Further development of the method of quantum hydrodynamics in application for Bose-Einstein condensates (BECs) is presented. To consider evolution of polarization direction along with particles movement we have developed corresponding set of quantum hydrodynamic equations. It includes equations of the polarization evolution and the polarization current evolution along with the continuity equation and the Euler equation (the momentum balance equation). Dispersion properties of the transverse waves including the electromagnetic waves propagating through the BECs are considered. To this end we consider full set of the Maxwell equations for description of electromagnetic field dynamics. This approximation gives us possibility to consider the electromagnetic waves along with the matter waves. We find a splitting of the electromagnetic wave dispersion on two branches. As a result we have four solutions, two for the electromagnetic waves and two for the matter waves, the last two are the concentration-polarization waves appearing as a generalization of the Bogoliubov mode. We also obtain that if the matter wave propagate perpendicular to external electric field when dipolar contribution does not disappear (as it follows from our generalization of the Bogoliubov spectrum). In this case exist a small dipolar frequency shift due to transverse electric field of perturbation.

###### pacs:

03.75.Kk 03.75.Hh^{†}

^{†}preprint: APS/123-QED

## I Introduction

At studying of electrically polarized Bose-Einstein condensates (BECs) the generalization of the Gross-Pitaevskii (GP) equation Goral PRA 00 ()-Fischer PRA 06R () and some its further generalization Wang NJP 08 ()-Lima PRA 12 () have been used, see also review papers Lahaye RPP 09 ()-Baranov CR 12 (). This nonlinear Schrodinger equation defines the complex scalar wave function , which describes evolution of the particle concentration . When direction of polarization of each particle d is not changing during the considering processes and all of them have same direction (the aligned dipoles) we can write formula for the density of polarization . Thus, the polarization density changes due to particles movement, and consequently, at changing of concentration. The waves of polarization can easily propagate even in the ferroelectric materials as well as the spin waves in ferromagnetic materials. All the more we can expect existence of such waves in the polarized BECs. To consider evolution of polarization direction along with particles movement we have developed the set of quantum hydrodynamic equations, which includes equations of the polarization evolution and the polarization current evolution along with the continuity equation and the Euler equation (the momentum balance equation) Andreev arxiv 12 02 ()-Andreev PRB 11 (). Earlier we developed similar approach for the electric dipolar ultracold fermions Andreev RPJ 13 ().

There is a field in dipolar BECs studying, where evolution of dipole directions is under consideration. This is research of spinor BECs, which was suggested by T. Ohmi and K. Machida Machida JPSJ 98 () and T. L. Ho Ho PRL 98 () in 1998 and was recently reviewed in Refs. Kawaguchi Ph Rep 12 (), Stamper-Kurn RMP 13 (). We should mention that this field is dedicated to magnetic dipole BECs, but the dipole-dipole (spin-spin) interaction is neglected there.

It was marked in Ref. Andreev arxiv 12 02 () that the full and correct form of Hamiltonian of the electric dipole interaction, which also corresponds to the Maxwell equations, is

(1) |

Justification of this formula in accordance with classical electrodynamics Landau 2 () is presented in the Appendix A.

Using well-known identity (see for instance Refs. O'Dell PRL 04 (), Bao JCP 10 ())

(2) |

we can see that the Hamiltonian (1) differs from usually used one

(3) |

Moreover, the Hamiltonian (1) is the original form of the potential energy of electric dipole interaction. Necessity to consider the Hamiltonian of electric dipole interaction in the form (1) caused by the fact that it must accord to the Maxwell equations. This connection exists since Maxwell equations describe the electric field and its connection with the sources. In our case source of the field is the density of electric polarization. Action of the electric field on the density of polarization comes in to equations of motion via the force field. As the result it describes interaction of the polarization (the electric dipole moments).

It has been shown in Refs. Santos PRL 03 ()-Giovanazzi EPJD 04 () that using of the non-linear Schrodinger equation, which generalized the GP equation for the electrically polarized BECs, and using of the formula (3) for description of the dipole-dipole interaction lead to the appearance of the contribution of the equilibrium polarization in the dispersion dependence of the Bogoliubov mode. This contribution leads to the anisotropy of the dispersion dependence. The anisotropy of the sound velocity in dipolar Bose gases has recently been confirmed experimentally (see Ref. Bismut PRL 12 ()).

In Refs. Andreev arxiv Pol (), Andreev RPJ 12 () we considered the evolution of electric dipole directions and got dispersion dependencies for two waves. That means what account of polarization evolution leads to appearing of new wave in the BECs. This is a polarization wave, which is an analog of spin waves in systems of spinning particles.

At the first step of studying of the polarization direction evolution we considered propagation of waves parallel to direction of electric field, so their dispersion reveals no angle dependence Andreev arxiv Pol (), Andreev RPJ 12 (). As an intermediate step in studying of the polarization evolution in BECs we considered fully polarized BECs. Usual theory of fully polarized BECs is based on formula (3) (see for instance Lahaye RPP 09 (), Baranov P R 08 (), Baranov CR 12 ()). It has the following consequences. Dipolar part of is proportional to , so it can change sign, where is the angle between the direction of wave propagation k and the direction of external field. If the polarization part is negative at positive and small enough short-range part of the spectrum, one finds the roton instability Lahaye RPP 09 (), Baranov P R 08 (), Baranov CR 12 (). Our consideration of fully polarized electric dipolar BECs is based on formula (1). Consequently the dipolar part of is proportional . So it is positive for all angles and shows no instability of three dimensional BECs with the repulsive short range interaction Andreev 1401.4842 (). Moreover full electric dipole interaction (1) leads to stabilization of the Bogoliubov spectrum with attractive short range interaction Andreev 1401.4842 (), since positive dipolar part can exceed negative contribution of the short range interaction. We obtained that consideration of the full Hamiltonian (1) is essential for fully polarized electric BECs, the full Hamiltonian of spin-spin interaction was considered for the fully polarized magnetic BECs as well. In both cases we have one wave solution only, as it was presented in Refs. Andreev 1401.4842 (). However the dipolar part of dispersion dependence changes due to account of the delta-function part of the full potential. As a result we find that the dipolar part in the electric dipolar BECs becomes positive for all angles (all directions of wave propagation). We have different picture for the magnetic dipolar BECs, where we have another coefficient in front of the delta-function (see formula (83.9) in Ref. LL 4 ()). Consequently the dipolar part of spectrum is also different. It is negative for all angles. Hence, there is no roton instability for the electrically polarized BECs, but the roton instability exists in the magnetically polarized BECs at all angles for small enough constant of the short-range interaction.

Consideration of the full electric dipolar BEC theory including the dipoles direction evolution was presented in Ref. Andreev arxiv 12 02 (), where we applied our model to calculation of dispersion of waves propagating at different angles . Thus Ref. Andreev arxiv 12 02 () contains generalization of results obtained in Ref. Andreev arxiv Pol (). Waves considered in Refs. Andreev arxiv 12 02 (), Andreev arxiv Pol (), Andreev RPJ 12 () were longitudinal, that means perturbations of the electric field is parallel to the direction of wave propagation . Many features of the quantum hydrodynamics of dipolar BECs and meaning of considered approximations were described in Ref. Andreev arxiv 12 02 () as well.

In this paper we suggest more general model in comparison with our recent papers. We consider full set of the Maxwell equations instead of the pair quasi-electrostatic equations describing the electric field considered earlier. This approximation gives us possibility to consider electromagnetic waves along with the matter waves. Below we will show that in the electrically polarized BECs we have a splitting of the electromagnetic waves on two branches. So, we have four solution, two for the electromagnetic waves and two for the matter waves, the last two are waves of concentration-polarization. It is important and interestingly to admit that in this approximation we find anisotropy in the dispersion dependencies for all four waves.

This paper is organized as follows. In Sec. II we present the quantum hydrodynamics model for description of the electrically polarized BECs. In Sec. III we present dispersion equation and describe the method to get it. In Sec. IV we present analysis of the dispersion equation. In Sec. V brief discussion of obtained results is presented.

## Ii Basic equations

In our previous paper Andreev arxiv 12 02 (), Andreev 2013 non-int GP () we have presented the derivation of the QHD equation for the electrically polarized BECs. These equations were directly derived from the microscopic many-particle Schrodinger equation by means the QHD method. Method of quantum hydrodynamics was suggested in 1999 for derivation of hydrodynamic equations for quantum plasmas MaksimovTMP 1999 (). In 2008 this method was applied to ultracold quantum gases Andreev PRA08 (). Later it was used for system of particles baring electric dipole moment Andreev arxiv 12 02 (), Andreev RPJ 12 (), Andreev PRB 11 (). Let us now present the set of the quantum hydrodynamics equations.

The first equation of the QHD equations system is the continuity equation

(4) |

showing conservation of the particle number, and giving time evolution of the particle concentration.

The momentum balance equation (Euler equation) for the polarized BECs has the form

(5) |

where

(6) |

and

(7) |

are the numerical coefficients. In equation (5) we define a parameter as (7). This definition differs from the one in the work Andreev PRA08 (). Here we put multiplier within the definition of to make equations looks better. Terms proportional to appear as a result of usage of the quantum kinematics. They are usually called the quantum Bohm potential. The first two terms at the right-hand side of the equation (5) are first terms of expansion of the quantum stress tensor. They occur due to account of the short-range interaction potential . The interaction potential determinates the macroscopic interaction constants and . Term proportional contains higher spatial derivation in compare with the Gross-Pitaevskii term . This is why it is called a non-local interaction. Contributions of the third order by the radius of the short range interaction (terms proportional ) in non-linear properties of quantum gases were considered in Refs. Andreev MPL B 12 (), Zezyulin EPJ D 13 () for non-polarized quantum gases. Influence of on dispersion of waves in dipolar fermions was described in Ref. Andreev RPJ 13 (). A non-local interaction for BECs being used as an analogue gravity system in Ref. Sarkar ariv 13 (). The last term of the equation (5) describes force field that affects the dipole moment in a unit of volume as the effect of the external electrical field and the field produced by other dipoles. It is written using the self-consistent field approximation Andreev arxiv 12 02 (), Andreev PRB 11 (), MaksimovTMP 1999 (). is the tensor of the kinetic pressure, which depends on particle thermal velocities and makes no contribution into the BEC dynamics at near zero temperatures. In equation (5) and below E is the sum of external and internal electric fields, where internal electric field is created by dipoles. P is the density of electric dipole moment (the polarization of medium).

The first order by the interaction radius constant for dilute gases may be presented in the form

(8) |

where is the scattering length, so we see , where is the Gross-Pitaevskii interaction constant L.P.Pitaevskii RMP 99 (); Andreev PRA08 ().

We have also the field equations

(9) |

and

(10) |

These equations show relation between the polarization P and the electric field E caused by the polarization P. Equation (9) allows to consider propagation of waves parallel to the equilibrium polarization only. For consideration of waves propagating in arbitrary direction we have to consider the couple of equations (9) and (10).

Having a short-range interaction in microscopic many-particle Schrodinger equation we obtain corresponding force field in the Euler equation as result of application of the QHD method for derivation of equations for collective evolution description Andreev PRA08 (). Using the fact that interaction is short-range we find the force field as a series. The force field also appears as the divergence of a second rank tensor, which is the quantum stress tensor . Considering the first two non-zero terms in the series, existing for spherically symmetric potential of the short-range interaction for particles in the BEC state, we get the Gross-Pitaevskii term and a non-local term containing the third spatial derivative of the square of the particle concentration , with . As a consequence we see that there is no corresponding non-linear Schrodinger equation, since there is no Cauchy integral of the Euler equation. Nevertheless a formal representation of the hydrodynamic equations in form of a non-linear Schrodinger equation was made in Ref. Andreev PRA08 ().

One may try to include non-local terms in the GP equation or the corresponding Lagrangian including spatial derivatives of the particle concentration. For instance one may consider instead of . However, a microscopic derivation of QHD equations shows that the Euler equation has more complicate form. We found instead of corresponding to a non-linear Schrodinger equation.

If particles do not contain the dipole moment, the continuity equation and the momentum balance equation form a closed set of equations. If the dipole moment is taken into account in the momentum balance equation when a new physical quantity emerges. This is the polarization vector field . It makes the set of equations to become incomplete.

If we have a system of fully polarized particles pictured on Fig. (1) when polarization changes due to motion of particles in space. Consequently we have , where is the electric dipole moment of particles, l is the direction of all dipoles. A non-linear Schrodinger equation can be derived for system of fully polarized dipoles from the first principal derived QHD equations Andreev PRA08 (). In an integral form it can be written as

(11) |

with is the macroscopic wave function, is the chemical potential, is the potential of an external field, is the constant of short-range interaction, is the dipole electric moment of a single particle, is the angle between and direction of external electric field, is the mass of a particle and is the Planck constant divided by . This is a generalized GP equation, where we included the Dirac delta function part of the full potential of electric dipole interaction, as we did earlier in Appendix B of Ref. Andreev arxiv 12 02 (). Let us admit that the macroscopic wave function appears to be defined in hydrodynamic variables as follows with is the potential of the velocity field .

The Delta function term in equation (11) can be combined with the short range interaction term presented by the term with the cubic nonlinearity. So we might obtain shifted constant of the short range interaction . But this is unwise move. We can go other way allowing to simplify equation (11) without extra approximations. We can explicitly introduce the electric field created by dipoles (see Appendix B formula (54)) and rewrite the non-linear Schrodinger equation in a non-integral form Andreev 2013 non-int GP ()

(12) |

We do not included the chemical potential in equation (12). Equation (12) contains an extra function , which obey the quasi electrostatic Maxwell equations (9) and (10). So we have a set of coupled equation for medium motion (12) and field (9) and (10). The non-integral form (12) of the generalized GP equation (11) was introduced in Ref. Andreev 2013 non-int GP (). See also discussion of this equation in Ref. Andreev 1401.4842 (). Let us admit that in this paper we discuss theory of dipolar BECs in model of point-like particles. Finite size particles were considered in Ref. Andreev 1401.4842 (), where a radius of molecules explicitly comes in the model.

The short range interaction constant in equation (12) does not depend on electric dipole moment . Hence the corresponding scattering length does not depend on either.

The dipole dependent scattering length for the short range part of interaction in dipolar BECs is under discussion in Ref. Bortolotti PRL 06 (). Authors consider a potential, which is close to the real one, with simplified treatment of the short range interaction, which considered as infinite repulsion at a cutoff radius (see formula (2) of Ref. Bortolotti PRL 06 ()). Considering GP equation they pay attention to an appropriate pseudopotential, its part proportional to the delta function is considered to be dependent on the dipole moment. So they have a dipole dependent scattering length. Due to numerical analysis it is hard to get an appropriate picture of dependence of the scattering length on the dipole moment and the cutoff radius. Nevertheless consideration of a close to realistic potential is related to a theory of finite size dipolar particles presented in Andreev 1401.4842 (), where dependence on particle finite size and dipole moment is presented explicitly and analytically. No pseudopotential was used in Ref. Andreev 1401.4842 () applying the quantum hydrodynamic method. Application of the dipole dependent scattering length along with the no cutoff of the dipole-dipole potential for finite particles means that contains a contribution aimed to cancel an endowment of dipole-dipole interaction inside of the sphere with the cutoff radius. contains dependence on cutoff radius as well.

Refs. Bortolotti PRL 06 () and Kanjilal PRA 07 () presented appearing at consideration of dipolar BECs in terms of the pseudopotential. This dependence is mentioned in reviews on dipolar BECs (see for instance formula (3) of Ref. Baranov P R 08 () and text after it, and formula (4.6) of Ref. Lahaye RPP 09 () and its discussion below). Unfortunately dipole dependence of the interaction constant is not included in discussions of spectrum of bulk collective Bogoliubov excitations. These discussions (see Ref. Baranov P R 08 () text after formula (6), Baranov CR 12 () text after formula (11), Lahaye RPP 09 () subsection (5.1)) of dipole dependence of the Bogoliubov spectrum happens as there is no dipole dependent scattering length , but usual scattering length of dipoleless particles.

Getting back to our analysis we separated the short range interaction and the dipole-dipole interaction with the beginning of our paper. So we have the scattering length , which does not depend on . All dependencies on the electric dipole moment are included in the dipole-dipole interaction.

Generalized integral GP equation, but with the shorted potential of dipole-dipole interaction (3), still finds a lot of applications. Frequencies of collective excitations in trapped dipolar BECs were calculated in the Thomas-Fermi regime in Ref. van Bijnen PRA 10 (). One-dimensional solitons on a weak two-dimensional square and triangular optical lattice potentials placed perpendicular to the polarization direction were considered in Ref. Adhikari PL A 12 (). Order and chaos in the dipolar BECs are investigated in Ref. Koberle NJP 09 () considering behavior of bifurcations. A dipolar boson-fermion mixture and its stability was studied in Ref. Adhikari PRA 13 (). When stability of dipolar Bose and Fermi gases was investigated in Ref. Adhikari JP B 13 (). Solitons in dipolar boson-boson mixtures were described in Ref. Adhikari JP B 14 (). Influence of dipole-dipole interaction on the formation of vortices in a rotating dipolar BECs atoms in quasi two-dimensional geometry was considered in Ref. Kumar JP B 12 (). Solitons in dipolar BECs trapped in the two dimensional plane were studied in Ref. Kumar JP B 13 (). Modification of the properties of super-fluid vortices by dipole-dipole interactions within the context of a two-dimensional atomic Bose gas of fully polarized dipoles was considered in Ref. Mulkerin PRL (). A discussion of linear properties of dipolar BEC based on integral GP equation with shorted dipole-dipole potential (3) is also presented in Ref. Mulkerin PRL (). Judging upon formula (2) in Ref. Qiuzi Li PRB 11 () we assume that potential (3) was used for consideration of collective modes in fermionic dipolar liquid.

We can also find explicit examples of using of identity (2) when authors work with the shorted potential of dipole-dipole interaction (3) see text before formula (418) in review Kawaguchi Ph Rep 12 (), see formula (36) in Ref. Lazur PRA 10 (), formula (8) in Ref. Eberlein PRA 05 (), formula (8) and text before this formula in Ref. O'Dell PRL 04 (), and formula (4) in Ref. Cai PRA 10 ().

Other approaches to dipolar BECs were developed in last years Lima PRA 12 (), Lima PRA 10 (), Babadi PRA 12 (). Some of them we briefly describe below. Various beyond mean-field effects on Bose gases at zero temperature were considered within the realm of the Bogoliubov-de Gennes theory Lima PRA 12 (). For the homogeneous system, the condensate depletion, the ground-state energy, the equation of state, and the speed of sound are discussed in detail. Authors analyze the influence of quantum fluctuations on properties of quantum gases. During calculations the condensate density being replaced by the total number density, which includes distribution of particles on states with non minimal energy. The fluctuation Hamiltonian density is included. The corrected sound velocity due to beyond mean-field equation of state, which is a generalization of the Lee-Huang-Yang quantum corrected equation of state, is presented by formula (34). Corresponding hydrodynamic equations are presented by formulas (32) and (33) of Ref. Lima PRA 12 (). Applying the many-particle Schrodinger equation for ultracold dipolar Fermi gases in anisotropic traps, which are polarized in the z direction, the couple of quantum hydrodynamic equations, i.e. the continuity and Euler equations, was obtained via the equation of motion for the one-body density matrix Lima PRA 10 (). Collective motion of polarized dipolar Fermi gases in the hydrodynamic regime was considered in Ref. Lima PRA 10R () applying a variational time-dependent Hartree-Fock approach.

Collective excitations of quasi-two-dimensional trapped dipolar fermions in framework of the collisional Boltzmann-Vlasov equation obtained by reducing the evolution equations of the nonequilibrium Green functions to quantum kinetic equations were considered in Ref. Babadi PRA 12 (). This analysis is performed at finite temperatures. Transition from collisionless to hydrodynamic regime was considered.

In this paper we are interested in dipole direction evolution instead of consideration of fully polarized BECs. So let us make a step towards the main topic. The variation of dipole direction can be different in different areas of space. Consequently we have different magnitude of polarization in different points of space. This variations can cause a turning of the polarization vector P relative to the direction of the external electric field (see Fig. (2 a)) or it can cause change of module of P without change of direction (see Fig. (2 b)).

If particles possessing electric dipole moment are in different cells of an optical lattice in the Mott insulator regime when they ability to get close to each other is limited. However there is no such limitation in the superfluid phase or when particles are in the absence of a lattice. So particles can get close to each other. In this case magnitude of the dipole-dipole interaction increases and direction of dipoles can change due to the dipole-dipole interaction. It gives possibility for polarization wave propagation, which are waves related to evolution of dipole direction Fig. (3).

All of it engage us to include dipole direction evolution in our model. Consequently we derive next equation of the chain of QHD equations, which is the equation for time evolution of the polarization (density of the electric dipoles) P. Explicit definitions of the QHD variables allowing to derive the set of QHD equations are presented in Appendix B.

Next equation we need for studying of the collective excitation dispersion is the equation of polarization evolution

(13) |

where is the current of polarization.

The equation (13) does not contain information about the effect of the interaction on the polarization evolution. So we should consider equation for evolution, which can be derived by means of the quantum hydrodynamic method Andreev PRA08 (), Andreev PRB 11 (). Using the self-consistent field approximation of the dipole-dipole interaction we obtain the equation for the polarization current evolution

(14) |

where presents the contribution of the thermal movement of the polarized particles into the dynamics of . As we deal with the BECs below, the contribution of may be neglected. The last term of the formula (14) includes both an external electrical field and the self-consistent field that particle dipoles create. This term contain numerical constant . The second group of terms in the left-hand side of equation (14) has kinematic nature. It is an analog of in the Euler equation (5). The fourth and fifth group of terms are the contribution of the quantum Bohm potential in the polarization evolution. Their appearance is related to the de-Broglie wave nature of quantum particles. The first term in the right-hand side of equation (14) describes the short range interaction. We restrict ourself with the first order by the interaction radius approximation for the short range interaction in the equation for evolution. That corresponds the Gross-Pitaevskii approach to the short range interaction .

The quasi-electrostatic Maxwell equations (9) and (10) emerge in the considered non-relativistic limit and allow to study longitudinal waves only. When we use term longitudinal we assume that direction of the electric field perturbation is parallel to the direction of the wave propagation. This case has been considered in Refs. Andreev arxiv 12 02 ()-Andreev RPJ 12 (). It is well-known that equations (9) and (10) are the part of the set of the Maxwell equations. Therefore, if we want to consider the transverse waves we should use the full set of the Maxwell equations, which has well-known form

(15) |

(16) |

(17) |

and

(18) |

where is the speed of light, and B is the vector of magnetic field. From equation (17) we see . Hence we have contribution of along with . When can be considered by means quasi electrostatic equations (9) and (10) only. The fundamental fact that the time evolution of polarization leads to generation of the magnetic field is also included in described model by means of equation (18).

## Iii Collective excitations in the electric polarized BEC: Dispersion equation

We can analyze the linear dynamics of collective excitations in the polarized BEC using the QHD equations (4), (5), (13), (14) and the Maxwell equations (15)-(18). Let us assume the system is placed in an external electrical field . The values of concentration and polarization for the system in an equilibrium state are constant and uniform and equilibrium velocity field , magnetic field and tensor equal zero.

We consider the small perturbation of equilibrium state

(19) |

Substituting these relations into system of equations (4), (5), (13), (14) and (15)-(18) and neglecting nonlinear terms, we obtain a system of the linear homogeneous equations in partial derivatives with constant coefficients. Passing to the following representation for small perturbations

yields the homogeneous system of algebraic equations. The electric field strength is assumed to have a nonzero value. Expressing all the quantities entering the system of equations in terms of the electric field, we come to the equation

(20) |

where

(21) |

with . is the dielectric permittivity tensor of electric dipolar BECs.

We suppose that the external electric field parallel to the z-axis. Thereby, the equilibrium polarization also parallel to the z-axis . In this case tensor (21) becomes more simple

(22) |

where

(23) |

Dispersion equation emerges as condition of existence of nonzero perturbations of electric field. That means the determinant of matrix (22) equals to zero: . Using formula we find an explicit form of the dispersion equation

(24) |

## Iv Dispersion dependencies

### iv.1 Matter waves

Equation (24) simplifies in the low frequency limit ()

(25) |

If is not close to zero, when the first term in the large brackets is mach more than the second one. Consequently equation (25) simplifies to

(26) |

at we have

(27) |

Equation (27) was obtained in Refs. Andreev arxiv 12 02 (), Andreev arxiv Pol (), Andreev RPJ 12 (). Equation (26) was also obtained in Ref. Andreev arxiv 12 02 (), its properties were considered under condition of maximal polarization contribution, that corresponds . So properties of equation (27) were studied in full details in Ref. Andreev arxiv 12 02 ().

Equation (26) has following solution

(28) |

In the case of wave propagation parallel to the equilibrium polarization solution we get solution obtained in Refs. Andreev arxiv 12 02 (), Andreev arxiv Pol (), Andreev RPJ 12 (), and at great length considered in Ref. Andreev arxiv 12 02 (). The contribution of equilibrium polarization becomes smaller at increasing of angle from to and vanishes at . However equation (26) does not describe area around . This area is described by equation (32) and below.

Analysis of formula (28) was presented in Ref. Andreev arxiv 12 02 () for . This analysis includes contribution of the short range interaction up to the third order by the interaction radius terms proportional . So let us discuss some properties of formula (28) in the first order by the interaction radius of the short-range interaction (neglecting terms proportional ). We also assume and use . Applying all this approximations we obtain two separate formulas from formula (28):

(29) |

and

(30) |

Assuming we can rewrite dipolar part of spectrum as . Formula (30) is a generalization of the Bogoliubov spectrum for dipolar particles. Extra solution (29) appears as a consequence of polarization evolution. Solutions (29) and (30) are presented on Fig. (4).

Solution (30) is rather different from the well-known result containing signature of the roton instability (see for instance formula (5.1) in Ref. Lahaye RPP 09 (), formula next to formula (6) in Ref. Baranov P R 08 (), formula (11) in Ref. Baranov CR 12 (), and formula (416) in Ref. Kawaguchi Ph Rep 12 ()). Experimental results in Ref. Bismut PRL 12 () were obtained for the magnetic dipolar BECs. Consequently our analysis does not relate to this experiment. However we should admit that our previous results for the fully polarized BECs based on the full potential of magnetic dipole interaction gives an angle dependence of the spectrum corresponding to one obtained experimentally in Ref. Bismut PRL 12 (). Hence our result (30) does not contradict to presence of the roton instability in spectrum of magnetic dipolar BECs. The Bogoliubov spectrum for the dipolar BECs presented in Refs. Lahaye RPP 09 (), Baranov P R 08 (), Baranov CR 12 (), Kawaguchi Ph Rep 12 () has the following form

(31) |

is the dipolar coupling constant in the SI units, where is the vacuum permittivity, or in the CGS units. It has been applied to both the magnetic and the electric dipolar BECs. In our recent papers we shown that generalization of the Bogoliubov spectrum are different for the electric and the magnetic dipolar BECs even for the fully polarized BECs Andreev 1401.4842 (), Andreev 2013 non-int GP (). These generalizations differ from formula (31). These differences appear due to the consideration of the full potentials of dipole-dipole interactions, which contain the Dirac delta function term (see formulas (1) and (2) for the electric dipoles). Coefficients before the delta functions are different for the electric and magnetic dipole interactions. The delta function term in the electric dipole interaction gives an additional repulsion, when the delta function term in the magnetic dipole interaction gives an extra attraction (see also discussions in Ref. Andreev 1401.4842 ()). Nevertheless formula (30) corresponds to the result obtained for the fully polarized electric dipolar BECs Andreev 1401.4842 (), Andreev 2013 non-int GP (), where the full potential of dipole-dipole interaction (1) was applied.

At equation (25) appears as

(32) |

We did not have such limit case as equation (32) when we considered longitudinal waves Andreev arxiv 12 02 (), Andreev arxiv Pol (), Andreev RPJ 12 (). Appearance of equation (32) is related to consideration of the full set of Maxwell equations (15)-(18) and presence of the perturbation of the electric field perpendicular to the direction of wave propagation. As a signature of presence of the transverse electric field we have the light speed in equation (32) and its solutions below.

Explicit form of equation (32) is

(33) |

Neglecting contributions of the quantum Bohm potential (terms proportional ) and the third order by the interaction radius (terms proportional ) we find the following solution

(34) |

Let us admit that solution (34) is obtained for . If formula (34) simplifies to and , the first of them corresponds to the Bogoliubov spectrum. Getting the first order amendment on , in the case of small polarization, we obtain

(35) |

and

(36) |

In opposite limit, when , we find (modified Bogoliubov mode) and

(37) |

The last solution has module of frequency much smaller than frequencies of the Bogoliubov mode. Solution (37) shows an aperiodic damping and no wave solution since . On the other hand it might cause instability of dipolar BEC cloud to collapse.

### iv.2 Electromagnetic waves

In this subsection we do not use assumption of small frequencies. On the contrary we focus our attention on properties of high frequency excitations. If formula (24) gives one electromagnetic wave solution only: . In this case we do not have influence of the medium on properties of the electromagnetic waves. Let us admit that the model under consideration does not account resonance interaction of the electromagnetic waves with medium when frequencies of the waves cause electron transitions within atoms.

Under the condition or equation (24) assumes the form

(38) |

From this equation we find usual the dispersion dependence of the light , and equation describes the two matter waves, which appear instead of the Bogoliubov mode in the electrically polarized BECs, see formulas (28)-(30) of this paper and Refs. Andreev arxiv 12 02 (), Andreev arxiv Pol (), Andreev RPJ 12 ().

In the opposite case, when equation (24) simplifies to

(39) |

In this case we have two independent equation. One of them describes changeless dispersion of the light. The second equation

(40) |

is the equation of the third degree of . So, it contains three wave branches. We can admit that it contains both modified by the BEC dispersion of the light and the two matter wave. This equation has no low frequencies limitation as (26), thus we have here high frequencies ”tail” in the dispersion of the matter waves.

Let us consider influence of the BEC on the light as a small correction. Therefore we put in of equation (40), and accounting that , , , in the result we find

(41) |

where . Expanding formula (41) in a series by we get

(42) |

with . We can see that the polarized BEC caused splitting of the light on two waves at light propagation perpendicular to the direction of the equilibrium polarization. Magnitude of splitting is equal to

(43) |

From this formula we find so frequency of the second wave is increased in compare with the : . However this shift is very small for cm, CGS units, g.

## V Conclusion

We have considered transverse waves in the electrically polarized BEC. To do it we have used whole set of the Maxwell equation for description of electromagnetic field caused by the dynamic of electric dipoles of the medium. In comparison to the case of longitudinal waves considered in our previous papers we have found an advantage, which is a possibility of consideration of the light propagating through the electrically polarized BECs. Thus, we have considered dispersion of the light propagating through the medium and we have found the splitting of the light on two waves in the case when the light propagate at the angle to the direction of the equilibrium polarization. We have calculated the magnitude of frequencies splitting for the case of the light propagation perpendicular to the direction of equilibrium polarization. There is no splitting in the case of the light propagation along the direction of the equilibrium polarization, there is also no contribution of the medium in the light dispersion dependence in this case.

As in the case of longitudinal waves we have got two matter waves, which exist instead of the Bogoliubov mode existing in the unpolarized BECs and fully polarized BECs. We have found analytical solution for the plane matter waves propagating in the three dimensional BECs and shown the influence of the anisotropy on the value of the frequencies.

In the result, we have developed generalization of the quantum hydrodynamic equations for the electrically polarized BECs including transverse wave propagation. Using it we have shown that in the electrically polarized BEC situated in the external electric field there are four wave, two of them are high frequencies waves, and them associated with the light. Two more wave are the anisotropic matter waves.

## Vi Appendix A: potential of electric dipole interaction

The potential of electric field created by a single electric dipole appears as

(46) |

(see for instance Landau 2 ()). We can take the derivative in formula (46) and obtain more explicit form of . However it is more illuminating to keep in form (46). The potential energy of an electric dipole in the electric field of another dipole is , where is the electric field created by the electric dipole in the position of the dipole . Using (46) we find

(47) |

In tensor notations it can be presented as . Let us have a look on the differential operator in formula (47). It is a tensor operator emerges as . The explicit form of this operator is presented by formula (2). We may now apply formula (2) to get explicit form of electric field created by dipoles. As the result we find difference with usually used formula (3).

Real atoms and molecules have a finite size, which can be presented by hard-core short range potential. It leads to cutting off of potential of the dipole-dipole interaction at small interparticle distances

(48) |

where is the potential of the long range dipole-dipole interaction, and is the potential of the short range interaction presenting a field surrounding a neutral particle. Potential is short range means that at . See Ref. Andreev 1401.4842 () for the Bogoliubov spectrum of the finite size dipolar BECs.

For non-dipolar BECs we have

(49) |

So we have the following constant of the short range interaction for the cut off potential .

## Vii Appendix B: definition of quantum hydrodynamic variables

The concentration of particles is defined as the quantum average of the concentration operator in the coordinate representation with the sum over all particles in the system:

(50) |

where is a vector in dimensional configurational space containing set of coordinates of all particles in the system under consideration, is the volume element in dimensional space, is the wave function describing exact evolution of particle quantum system. obeys the many-particle Schrodinger equation. The operator of particle concentration contains the Dirac delta functions, which provide projection of dimensional quantum dynamics in three dimensional physical space, where real particles move.

The particle flux or the momentum density appears as

(51) |

where is the momentum operator of the th particle, and stands for the complex conjugation. The particle flux allows to derive the Euler equation applying .

Polarization appears in the Euler equation as follows

(52) |

where is the operator of electric dipole moment of a neutral particle.

The polarization current or polarization flux is

(53) |

The internal electric field caused by dipoles has the following explicit form

(54) |

To derive QHD equations we need the Schrodinger equation , which governs evolution of the many-particle wave function. We also need explicit form of the Hamiltonian of the system under consideration. Which within the quasi-static approximation is

(55) |

The first term in the Hamiltonian is the operator of the kinetic energy. The second term represents the interaction between the dipole moment and the external electric field. The subsequent terms represent the short-range and the long range dipole-dipole interactions between neutral particles baring electric dipole moment. The Green function for the dipole-dipole interaction reads as .

At derivation of the QHD equations we use quasi static approximation for the dipole-dipole interaction. It means that at each moment of time we consider dipole-dipole interaction between motionless dipoles. This interaction changes state of the translational motion of molecules and direction of their dipoles. At next moment of time particles are in new positions and having new direction of dipoles. In this new state particles interact as motionless, etc. Quasi static interaction is an interaction via the electric field created by motionless dipoles. Slow motion of dipoles creates small magnetic field and electromagnetic radiation. Which can be accounted by consideration of the full set of Maxwell equations, as we do it in the end of section II.

More details on derivation of quantum collective observable evolution by quantum hydrodynamic method can be found in the following Refs. Andreev PRB 11 (), MaksimovTMP 1999 () (focused on dipole evolution), Andreev PRA08 () (focused on BECs evolution).

## References

- (1) K. Goral, K. Rzazewski, and T. Pfau, Phys. Rev. A 61, 051601(R) (2000).
- (2) L. Santos, G.V. Shlyapnikov, P. Zoller, and M. Lewenstein, Phys. Rev. Lett. 85, 1791 (2000).
- (3) S. Yi and L. You, Phys. Rev. A, 61, 041604(R) (2000).
- (4) C. Ticknor, R. M. Wilson, and J. L. Bohn, Phys. Rev. Lett. 106, 065301 (2011).
- (5) Uwe R. Fischer, Phys. Rev. A 73, 031602(R) (2006).
- (6) Daw-Wei Wang, New J. Phys. 10, 053005 (2008).
- (7) Chao-Chun Huang, Daw-Wei Wang, and Wen-Chin Wu, Phys. Rev. A 81, 043629 (2010).
- (8) R. M. Wilson, S. T. Rittenhouse and J. L. Bohn, New J. Phys. 14, 043081 (2012).
- (9) A. R. P. Lima and A. Pelster, Phys. Rev. A 86, 063609 (2012).
- (10) T. Lahaye, C. Menotti, L. Santos, M. Lewenstein and T. Pfau, Rep. Prog. Phys. 72, 126401 (2009).
- (11) M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen(De), and U. Sen, Advances in Physics, 56, 243 (2007).
- (12) Lincoln D. Carr, David DeMille, Roman V. Krems and Jun Ye, New J. Phys. 11, 055049 (2009).
- (13) K.-K. Ni, S. Ospelkaus, D. J. Nesbitt, J. Ye and D. S. Jin, Phys. Chem. Chem. Phys. 11, 9626 (2009).
- (14) M. A. Baranov, Physics Reports 464, 71 (2008).
- (15) M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, Chem. Rev., 112, 5012 (2012).
- (16) P. A. Andreev and L. S. Kuz’menkov, Eur. Phys. J. D 67, 216 (2013).
- (17) P. A. Andreev and L. S. Kuzmenkov, arXiv:1106.0822.
- (18) P. A. Andreev, Russian Physics Journal 54, 1360 (2012).
- (19) P. A. Andreev, L. S. Kuzmenkov and M. I. Trukhanova, Phys. Rev. B 84, 245401 (2011).
- (20) P. A. Andreev, Russian Physics Journal 55, 1190 (2013).
- (21) T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998).
- (22) T. L. Ho, Phys. Rev. Lett. 81, 742 (1998).
- (23) Yuki Kawaguchi, Masahito Ueda, Physics Reports 520, 253 (2012).
- (24) Dan M. Stamper-Kurn, Masahito Ueda, Rev. Mod. Phys. 85, 1191 (2013).
- (25) L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Butterworth-Heinemann, 1975).
- (26) D. H. J. O’Dell, S. Giovanazzi, and C. Eberlein, Phys. Rev. Lett. 92, 250401 (2004).
- (27) W. Bao, Y. Cai, and H. Wang, J. Comput. Phys. 229, 7874 (2010).
- (28) L. Santos, G. V. Shlyapnikov, and M. Lewenstein, Phys. Rev. Lett. 90, 250403 (2003).
- (29) D. H. J. O’Dell, S. Giovanazzi, and G. Kurizki, Phys. Rev. Lett. 90, 110402 (2003).
- (30) S. Giovanazzi and D. H. J. O’Dell, Eur. Phys. J. D 31, 439 (2004).
- (31) G. Bismut, B. Laburthe-Tolra, E. Marechal, P. Pedri, O. Gorceix, and L. Vernac, Phys. Rev. Lett. 109, 155302 (2012).
- (32) P. A. Andreev, L. S. Kuz’menkov, arXiv:1401.4842.
- (33) V. B. Berestetskii, E. M. Lifshitz, L. P. Pitaevskii, Quantum Electrodynamics. Vol. 4 (2nd ed.), Butterworth-Heinemann (1982).
- (34) P. A. Andreev, Mod. Phys. Lett. B 27, 1350096 (2013).
- (35) L. S. Kuz’menkov, S. G. Maksimov, Theor. Math. Phys., 118, 287 (1999).
- (36) P. A. Andreev, L. S. Kuzmenkov, Phys. Rev. A 78, 053624 (2008).
- (37) P. A. Andreev, L. S. Kuzmenkov, Mod. Phys. Lett. B 26, 1250152 (2012).
- (38) K. V. Zezyulin, P. A. Andreev, and L. S. Kuz’menkov, Eur. Phys. J. D 67, 140 (2013).
- (39) Supratik Sarkar, A. Bhattacharyay, arXiv:1307.8231.
- (40) F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999).
- (41) D. C. E. Bortolotti, S. Ronen, J. L. Bohn, and D. Blume, Phys. Rev. Lett. 97, 160402 (2006).
- (42) K. Kanjilal, John L. Bohn, and D. Blume, Phys. Rev. A 75, 052705 (2007).
- (43) R. M. W. van Bijnen, N. G. Parker, S. J. J. M. F. Kokkelmans, A. M. Martin, and D. H. J. O’Dell, Phys. Rev. A 82, 033612 (2010).
- (44) S. K. Adhikari, P. Muruganandam, Phys. Lett. A 376, 2200 (2012).
- (45) Patrick Koberle, Holger Cartarius, Tomaz Fabcic, Jorg Main and Gunter Wunner, New J. Phys. 11, 023017 (2009).
- (46) S. K. Adhikari, Phys. Rev. A 88, 043603 (2013).
- (47) S. K. Adhikari, J. Phys. B 46, 115301 (2013).
- (48) S. K. Adhikari and Luis E. Young-S., J. Phys. B 47, 015302 (2014).
- (49) R. Kishor Kumar, P. Muruganandam, J. Phys. B: At. Mol. Opt. Phys. 45, 215301 (2012).
- (50) R. Kishor Kumar, P. Muruganandam, B. A. Malomed, J. Phys. B: At. Mol. Opt. Phys. 46, 175302 (2013).
- (51) B. C. Mulkerin, R. M. W. van Bijnen, D. H. J. O’Dell, A. M. Martin, N. G. Parker, arXiv:1310.3640.
- (52) Qiuzi Li, E. H. Hwang, and S. Das Sarma, Phys. Rev. B 82, 235126 (2011).
- (53) V. Yu. Lazur, S. I. Myhalyna, and O. K. Reity, Phys. Rev. A 81, 062707 (2010).
- (54) C. Eberlein, S. Giovanazzi, D. H. J. O’Dell, Phys. Rev. A 71, 033618 (2005).
- (55) Y. Cai, M. Rosenkranz, Z. Lei, and W. Bao, Phys. Rev. A 82, 043623 (2010).
- (56) A. R. P. Lima, A. Pelster, Phys. Rev. A 81, 063629 (2010).
- (57) M. Babadi and E. Demler, Phys. Rev. A 86, 063638 (2012).
- (58) A. R. P. Lima, A. Pelster, Phys. Rev. A 81, 021606(R) (2010).