 stationary configuration. Let us take a constant background and examine small perturbations of the form $$a = a_0 + a_1({\bf r},t)~\textrm{and}~{\bf U} = {\bf U}_1({\bf r},t)\label{eq:35}$$ so we have $$\begin{aligned}
2 a_0  \frac{\partial a_1}{\partial t} + a_0^2 \nabla \cdot {\bf U}_1 &=& 0,\label{eq:36} \\
  \frac{\partial {\bf U_1}}{\partial t} + \nabla \left ( V_1 -  \frac{1}{2 a_0} \nabla^2 a_1 \right ) &=& 0. \label{eq:37}\end{aligned}$$ Now if we take the time derivative of Eq. (\[eq:36\]) and the divergence of Eq. (\[eq:37\]), we can combine the equations to yield $$2 a_0 \frac{\partial^2 a_1}{\partial t^2} - a_0^2 \nabla^2 V_1 + \frac{a_0}{2} \nabla^4 a_1 = 0
  \label{eq:38}$$ and $$\frac{\partial^2 a_1}{\partial t^2} - \frac{\alpha}{2} a_1 |a_0|^\alpha + \frac{1}{4} \nabla^4 a_1 = 0.
  \label{eq:39}$$ If we expand the perturbations in Fourier components we obtain the following dispersion relation $$\omega^2 =  \frac{k^4}{4} - \frac{\alpha}{2} |a_0|^\alpha 
  \label{eq:40}$$ where the first term is the standard result for the deBroglie wavelength of a particle and the second term is due to the self-gravity of the perturbation.

We can be a bit more sophisticated now and assume that small perturbations lie near a static solution so $$a = a_0({\bf r}) + a_1({\bf r},t)~\textrm{and}~{\bf U} = {\bf U}_1({\bf r},t)\label{eq:41}$$ thus we have $$\begin{aligned}
  2 a_0  \frac{\partial a_1}{\partial t} + \nabla \cdot \
stationary configuration. Let us take a constant setting and test small perturbations of the form $ $ a = a_0 + a_1({\bf r},t)~\textrm{and}~{\bf uranium } = { \bf U}_1({\bf r},t)\label{eq:35}$$ so we have $ $ \begin{aligned } 
 2 a_0   \frac{\partial a_1}{\partial t } + a_0 ^ 2 \nabla \cdot { \bf U}_1 & = & 0,\label{eq:36 } \\ 
   \frac{\partial { \bf U_1}}{\partial t } + \nabla \left (V_1 -   \frac{1}{2 a_0 } \nabla^2 a_1 \right) & = & 0. \label{eq:37}\end{aligned}$$ immediately if we assume the time derived function of Eq.   (\[eq:36\ ]) and the divergence of Eq.   (\[eq:37\ ]), we can combine the equation to yield $ $ 2 a_0 \frac{\partial^2 a_1}{\partial t^2 } - a_0 ^ 2 \nabla^2 V_1 + \frac{a_0}{2 } \nabla^4 a_1 = 0 
   \label{eq:38}$$ and $ $ \frac{\partial^2 a_1}{\partial t^2 } - \frac{\alpha}{2 } a_1 |a_0|^\alpha + \frac{1}{4 } \nabla^4 a_1 = 0. 
   \label{eq:39}$$ If we expand the perturbation in Fourier components we receive the following dispersion sexual intercourse $ $ \omega^2 =   \frac{k^4}{4 } - \frac{\alpha}{2 } |a_0|^\alpha 
   \label{eq:40}$$ where the first term is the standard result for the deBroglie wavelength of a atom and the second term is due to the self - graveness of the perturbation. 

 We can be a act more sophisticated immediately and assume that small perturbations dwell near a static solution so $ $ a = a_0({\bf r }) + a_1({\bf r},t)~\textrm{and}~{\bf U } = { \bf U}_1({\bf r},t)\label{eq:41}$$ thus we have $ $ \begin{aligned } 
   2 a_0   \frac{\partial a_1}{\partial t } + \nabla \cdot \
 stwtionary configuration. Ltt us take a consjabt bacngrouns and exxmine small perturbations of tye foem $$a = a_0 + a_1({\bf r},t)~\textrm{avd}~{\bf U} = {\bv U}_1({\bf r},t)\oabeo{wq:35}$$ so we hebe $$\begik{clignsf}
2 a_0  \hrac{\partial a_1}{\pattial t} + a_0^2 \ndbla \cdot {\bf U}_1 &=& 0,\uayel{eq:36} \\
  \frac{\partial {\bf U_1}}{\partial t} + \nwbla \legt ( V_1 -  \frac{1}{2 a_0} \nwbla^2 w_1 \rifht ) &=& 0. \label{eq:37}\end{aligned}$$ Now if we take tie time derivatove of Eq. (\[eq:36\]) and the divergfnce of Eq. (\[eq:37\]), we can colbine the ewuatyins to yield $$2 a_0 \frac{\pagcial^2 a_1}{\partizl t^2} - a_0^2 \nabla^2 V_1 + \frac{a_0}{2} \nabla^4 a_1 = 0
  \layel{eq:38}$$ and $$\ftcx{\pagjial^2 a_1}{\partiao t^2} - \srac{\alpha}{2} a_1 |a_0|^\alpha + \xrac{1}{4} \nanla^4 a_1 = 0.
  \label{cq:39}$$ If we expand the perturbatmons in Fourier compjnents we oytain the following duspersimn ralatkin $$\umefa^2 =  \frac{k^4}{4} - \fcac{\alpha}{2} |a_0|^\ampha 
  \labeo{eq:40}$$ where the first tqgk is the stahdard wefult for the deBroglie wavelength of a kartidle and the second term is due to the self-grwvity of ehe perturbation.

We can be a bit more sophisticateg now xnd afrymf that small perturbations lie near a static fklitpon so $$a = a_0({\bf r}) + c_1({\bf r},t)~\textrm{and}~{\bg K} = {\ff U}_1({\bf r},t)\labeu{eq:41}$$ thba se have $$\begin{alignfd}
  2 a_0  \frac{\partial a_1}{\kartisl t} + \nabla \cdot \
stationary configuration. Let us take a constant examine perturbations of form $$a = = U}_1({\bf r},t)\label{eq:35}$$ so have $$\begin{aligned} 2 \frac{\partial a_1}{\partial t} + a_0^2 \nabla {\bf U}_1 &=& 0,\label{eq:36} \\ \frac{\partial {\bf U_1}}{\partial t} + \nabla \left ( - \frac{1}{2 a_0} \nabla^2 a_1 \right ) &=& 0. \label{eq:37}\end{aligned}$$ Now if we the derivative Eq. and the divergence of Eq. (\[eq:37\]), we can combine the equations to yield $$2 a_0 \frac{\partial^2 t^2} - a_0^2 \nabla^2 V_1 + \frac{a_0}{2} \nabla^4 = 0 \label{eq:38}$$ and a_1}{\partial t^2} - \frac{\alpha}{2} a_1 + \nabla^4 a_1 0. If expand the perturbations Fourier components we obtain the following dispersion relation $$\omega^2 = \frac{k^4}{4} - \frac{\alpha}{2} |a_0|^\alpha \label{eq:40}$$ where the term is result for deBroglie of particle and the is due to the self-gravity of can be a bit more sophisticated now and that small lie near a static solution so = a_0({\bf r}) + a_1({\bf r},t)~\textrm{and}~{\bf U} = U}_1({\bf r},t)\label{eq:41}$$ thus we have $$\begin{aligned} 2 a_0 \frac{\partial a_1}{\partial t} + \nabla \cdot \
 stationary configuration. LeT us take a coNstanT baCkgRoUnd aNd exAmine small pertURbatIons of the form $$a = a_0 + a_1({\bf r},t)~\teXtrm{aNd}~{\BF U} = {\bf u}_1({\Bf R},t)\labEl{eq:35}$$ so wE HaVE $$\BegIn{AlIgnEd}
2 A_0  \FrAc{\parTiaL a_1}{\partiAl t} + a_0^2 \nabla \cDot {\Bf u}_1 &=& 0,\label{eq:36} \\
  \frac{\PArTial {\bf U_1}}{\parTiaL t} + \nabla \left ( V_1 -  \FraC{1}{2 a_0} \nablA^2 a_1 \RigHT ) &=& 0. \labeL{eq:37}\End{alIgned}$$ NOW if we tAke the timE dERivatiVE of Eq. (\[eq:36\]) AND tHe diVergence of Eq. (\[eq:37\]), we cAN cOMbine the equatiOns to yIeLD $$2 a_0 \FRAc{\pArtIal^2 a_1}{\partiaL t^2} - A_0^2 \nablA^2 v_1 + \frac{a_0}{2} \nABlA^4 A_1 = 0
  \LAbeL{Eq:38}$$ and $$\frac{\partIal^2 a_1}{\partial T^2} - \FraC{\alpha}{2} A_1 |a_0|^\AlpHA + \frac{1}{4} \nAbla^4 a_1 = 0.
  \LaBEl{eQ:39}$$ If we expand The pErturbatiOns in FOUrier coMPonents We obtaIn tHe fOlloWInG dIspErSIon RElAtiON $$\omEga^2 =  \frac{k^4}{4} - \FrAc{\Alpha}{2} |A_0|^\alpHA 
  \LABel{eQ:40}$$ whEre tHe firSt term is the stAndArd rESulT for tHe deBRoglIe WavelEngth oF a parTiCle and the second Term Is due to thE seLf-GraViTy of tHE pertuRbaTioN.

We can bE a bit moRE soPhISTIcAted now and assume thAt SMAlL perturbAtions LIe NeAR a static SoLutIon sO $$A = A_0({\bf r}) + a_1({\Bf r},t)~\TExTrm{and}~{\bf u} = {\bf U}_1({\bf R},T)\lAbEl{eq:41}$$ thuS wE have $$\bEgIn{aLigNed}
  2 a_0  \fRAc{\paRtial a_1}{\Partial t} + \Nabla \CDot \
 stationary configuration.  Let us ta ke acon sta nt  bac kgro und and examin e  sma ll perturbations of th e for m$ $a = a_ 0 + a _1({\bf r} , t )~\ te xt rm{ an d }~ {\bfU}= {\bfU}_1({\bfr}, t) \label{eq:35 } $$  so we hav e $ $\begin{alig ned }
2 a_ 0 \f r ac{\p art ial a _1}{\p a rtialt} + a_0^ 2\ nabla\ cdot {\ b f  U }_1&=& 0,\label{eq:3 6 }\ \
  \frac{\par tial { \b f  U _ 1 }}{ \pa rtial t} +  \ nabla \left ( V_ 1 -   \ f rac{1}{2 a_0}  \nabla^2 a _ 1 \ right)&=& 0. \la bel{e q: 3 7}\ end{aligned }$$Now if we  taket he time derivat ive of  Eq . ( \[eq : 36 \] ) a nd the di ver g enc e of Eq.  ( \[ eq:37 \]), w e cancom bine  theequations toyie ld $ $ 2 a _0 \f rac{\ part ia l^2 a _1}{\p artia lt^2} - a_0^2 \n abla ^2 V_1 +\fr ac {a_ 0} {2} \ n abla^4  a_ 1 =  0
  \l abel{eq : 38} $$ a n d$$\frac{\partial^2  a _ 1 }{ \partial  t^2}-  \ fr a c{\alpha }{ 2}a_1| a _0|^\ alph a  +  \frac{1 }{4} \ n ab la ^4 a_1=0.
  \ la bel {eq :39}$ $  Ifwe exp and thepertu r bations in Fou r ier component s  w e ob t ainthe  followingdisp e rsio n re l at ion $$\om ega^2  =  \ f rac{k^4}{4} - \frac {\ alpha} {2} | a_0|^\alpha   \label{e q : 4 0}$$ whe re t h ef irst term is t he st andard res u lt for t he de Brogliewavelengt h of a par tic leand  th e se cond term isd u e to  t he self -gr avity o f t heper tur ba tion.

We  can beabi tmo resophi s ticatedno w a nd  as sumet hat sm all p ertu rb at i ons  lie ne a ra stat ic  s olut ion  s o $$a  = a _ 0({ \bf r})  + a_1({\ bfr },t) ~\ te xtrm{an d}~{\bf U} ={\ bf U}_1({\ bf  r} ,t)\la b e l{eq:41} $$ thus we have $$\begi n {aligne d}  2 a _0 \frac{\pa rti al a_1 }{\ p artial  t} +\nabl a\cd o t  \
 stationary_configuration. Let_us take a constant_background and_examine_small perturbations_of_the form $$a_= a_0 +_a_1({\bf r},t)~\textrm{and}~{\bf U} =_{\bf U}_1({\bf r},t)\label{eq:35}$$_so_we have $$\begin{aligned}
2 a_0  \frac{\partial a_1}{\partial t} + a_0^2 \nabla \cdot {\bf U}_1_&=&_0,\label{eq:36} \\
__\frac{\partial_{\bf U_1}}{\partial t} + \nabla_\left ( V_1 - _\frac{1}{2 a_0}_\nabla^2 a_1 \right ) &=& 0. \label{eq:37}\end{aligned}$$ Now_if_we take the_time derivative of Eq. (\[eq:36\]) and the divergence of Eq. (\[eq:37\]),_we can combine the equations to_yield $$2 a_0_\frac{\partial^2_a_1}{\partial_t^2} - a_0^2 \nabla^2_V_1 + \frac{a_0}{2} \nabla^4 a_1 =_0
  \label{eq:38}$$ and $$\frac{\partial^2 a_1}{\partial_t^2} - \frac{\alpha}{2} a_1 |a_0|^\alpha + \frac{1}{4}_\nabla^4 a_1 = 0.
  \label{eq:39}$$_If we expand the perturbations_in Fourier_components we obtain the following_dispersion relation $$\omega^2_= _\frac{k^4}{4} - \frac{\alpha}{2}_|a_0|^\alpha 
  \label{eq:40}$$ where the_first term is_the standard result for the deBroglie_wavelength_of a particle_and_the_second term_is due to_the_self-gravity of_the_perturbation.

We can be a bit more_sophisticated_now and assume that small perturbations lie_near a static solution_so_$$a = a_0({\bf r})_+ a_1({\bf r},t)~\textrm{and}~{\bf U} =_{\bf U}_1({\bf r},t)\label{eq:41}$$ thus we have_$$\begin{aligned}
 _2 a_0_ \frac{\partial a_1}{\partial t} + \nabla \cdot \