 V/\partial z\right)=
\left( \partial V^* /\partial z^* \right)$, we have from Eq. (\[A1\]) for the quadratic variables $(\delta z)^2$ and $|\delta z|^2$ the following equations of motion $$\begin{aligned}
\label{A2}
i \frac{d}{d\tau}\left(\delta z\right)^2 &=&
-i\Gamma \left(\delta z\right)^2 +
2\frac{\partial V}{\partial z} \left(\delta z\right)^2 +
2\frac{\partial V}{\partial z^*} |\delta z|^2,\nonumber \\
%
i \frac{d}{d\tau} |\delta z|^2 &=&
-i\Gamma |\delta z|^2
-\frac{\partial V^*}{\partial z} \left(\delta z\right)^2 +
\frac{\partial V}{\partial z^*} \left(\delta z^*\right)^2.\end{aligned}$$ It is easy to see that Eqs. (\[A2\]) are equivalent to Eqs. (\[matrix\]) with $\epsilon=0$, if one makes substitutions $|\delta z|^2\rightarrow B$ and $(\delta z)^2\rightarrow C$.

It should be noticed that yet exists one difference between the linearization of the classical motion equations and the equations for quantum cumulants (\[9b\]), (\[9c\]): It is impossible to get the initial conditions (\[initial\_cond\]) for $C$ and $B$ from only initial conditions for the linearized classical equations of motion (see also discussion of this problem in [@9; @9']).

 {#appendix2}

In this Appendix, we present the derivation of the time scale for the validity of the $1/N$-method for lossless ($\Gamma=0$) oscillators. Start with the case of Kerr nonlinearity $l=1$ and then generalize obtained results for arbitrary $l$. From Eqs. (\[z\_cl\]) and (\[B\_and\_C\]) for $\Gamma=0$ and $l=1$, we get $$\label{B1}
z_{cl}(\tau)=z_0\exp(-i\Omega\tau), \quad
\Omega\equiv\bar{\Delta}+|z_0|^2,$$ $$\label{B2}
C(\tau)=-\tau z_0^2 \left[ |z_
V/\partial z\right)= 
 \left (\partial V^ * /\partial z^ * \right)$, we have from Eq. (\[A1\ ]) for the quadratic variables $ (\delta z)^2 $ and $ |\delta z|^2 $ the following equations of apparent motion $ $ \begin{aligned } 
 \label{A2 } 
 i \frac{d}{d\tau}\left(\delta z\right)^2 & = & 
 -i\Gamma \left(\delta z\right)^2 + 
 2\frac{\partial V}{\partial z } \left(\delta z\right)^2 + 
 2\frac{\partial V}{\partial z^ * } |\delta z|^2,\nonumber \\ 
% 
 i \frac{d}{d\tau } |\delta z|^2 & = & 
 -i\Gamma |\delta z|^2 
 -\frac{\partial V^*}{\partial omega } \left(\delta z\right)^2 + 
 \frac{\partial V}{\partial z^ * } \left(\delta z^*\right)^2.\end{aligned}$$ It is easy to see that Eqs. (\[A2\ ]) are equivalent to Eqs. (\[matrix\ ]) with $ \epsilon=0 $, if one make substitution $ |\delta z|^2\rightarrow B$ and $ (\delta z)^2\rightarrow C$. 

 It should be noticed that yet exist one dispute between the linearization of the classical motion equations and the equations for quantum cumulants (\[9b\ ]), (\[9c\ ] ): It is impossible to catch the initial conditions (\[initial\_cond\ ]) for $ C$ and $ B$ from alone initial condition for the linearized classical equations of motion (see besides discussion of this problem in [ @9; @9' ]). 

  { # appendix2 } 

 In this Appendix, we present the derivation of the fourth dimension scale for the validity of the $ 1 / N$-method for lossless ($ \Gamma=0 $) oscillator. Start with the sheath of Kerr nonlinearity $ l=1 $ and then generalize receive results for arbitrary $ l$. From Eqs. (\[z\_cl\ ]) and (\[B\_and\_C\ ]) for $ \Gamma=0 $ and $ l=1 $, we get $ $ \label{B1 } 
 z_{cl}(\tau)=z_0\exp(-i\Omega\tau), \quad 
 \Omega\equiv\bar{\Delta}+|z_0|^2,$$ $ $ \label{B2 } 
 C(\tau)=-\tau z_0 ^ 2 \left [ |z _
 V/\pwrtial z\right)=
\left( \partiau V^* /\partial z^* \rntht)$, we have rrom Eq. (\[X1\]) for the quadratic variabled $(\eelta z)^2$ and $|\delta z|^2$ the foluowing eqlations od mouion $$\begin{aligned}
\label{A2}
i \nxac{d}{d\fwu}\leyt(\velta z\right)^2 &=&
-i\Gsmma \left(\dalta z\right)^2 +
2\frdc{\oaxtial V}{\partial z} \left(\delta z\right)^2 +
2\frwc{\partisl V}{\partial z^*} |\dejta e|^2,\nogumbsg \\
%
l \frac{d}{d\tau} |\delta z|^2 &=&
-i\Gamma |\delta z|^2
-\frac{\pertial V^*}{\partial z} \left(\delta z\right)^2 +
\frac{\pagtiap V}{\partial z^*} \left(\dflta z^*\right)^2.\ghd{ajugned}$$ It is dasy to set chat Eqs. (\[A2\]) zre equivalent to Eqs. (\[matrix\]) wigh $\epxilon=0$, if obe mangs substitutmons $|\dvlta z|^2\rightarvpw B$ atd $(\delts z)^2\rightarrow G$.

It sioule be noticed that yet exists one differgnce betwean the linearizatiob if thg clavsicxo mutikn esuatiojs end the equztions for wuantum cumulants (\[9b\]), (\[9c\]): Py is impossigle to gqt the initial conditions (\[initial\_cond\]) fmr $D$ and $B$ from only initiql conditions for the linearizqd classical equations of motion (see also discusshon oh ghiw irobuwm in [@9; @9']).

 {#appendix2}

In this Appendix, we present the setination of the timc scale for the vakifiyi of the $1/N$-methud for loasless ($\Gamma=0$) oscilpators. Ftart with the casr of Kerr nonlinearity $l=1$ ane then generclize obtained resultd for arbitxary $l$. From Eqs. (\[z\_cl\]) and (\[B\_and\_C\]) for $\Yamma=0$ znd $l=1$, we geh $$\label{B1}
z_{du}(\tau)=z_0\exp(-i\Omega\taj), \qlad
\Okega\equiv\bar{\Delta}+|z_0|^2,$$ $$\label{B2}
C(\eau)=-\tau z_0^2 \oeft[ |z_
V/\partial z\right)= \left( \partial V^* /\partial z^* have Eq. (\[A1\]) the quadratic variables the equations of motion \label{A2} i \frac{d}{d\tau}\left(\delta &=& -i\Gamma \left(\delta z\right)^2 + 2\frac{\partial z} \left(\delta z\right)^2 + 2\frac{\partial V}{\partial z^*} |\delta z|^2,\nonumber \\ % i \frac{d}{d\tau} z|^2 &=& -i\Gamma |\delta z|^2 -\frac{\partial V^*}{\partial z} \left(\delta z\right)^2 + \frac{\partial V}{\partial \left(\delta It easy see that Eqs. (\[A2\]) are equivalent to Eqs. (\[matrix\]) with $\epsilon=0$, if one makes substitutions $|\delta B$ and $(\delta z)^2\rightarrow C$. It should be that yet exists one between the linearization of the motion and the for cumulants (\[9c\]): It is to get the initial conditions (\[initial\_cond\]) for $C$ and $B$ from only initial conditions for the linearized equations of also discussion this in @9']). {#appendix2} In we present the derivation of the the validity of the $1/N$-method for lossless ($\Gamma=0$) Start with case of Kerr nonlinearity $l=1$ and generalize obtained results for arbitrary $l$. From Eqs. and (\[B\_and\_C\]) for $\Gamma=0$ and $l=1$, we get $$\label{B1} z_{cl}(\tau)=z_0\exp(-i\Omega\tau), \quad \Omega\equiv\bar{\Delta}+|z_0|^2,$$ $$\label{B2} C(\tau)=-\tau z_0^2
 V/\partial z\right)=
\left( \partial v^* /\partial z^* \rIght)$, wE haVe fRoM Eq. (\[A1\]) For tHe quadratic varIAbleS $(\delta z)^2$ and $|\delta z|^2$ the folLowinG eQUatiONs Of motIon $$\begiN{AlIGNed}
\LaBeL{A2}
i \FrAC{d}{D\tau}\lEft(\Delta z\rIght)^2 &=&
-i\Gamma \LefT(\dElta z\right)^2 +
2\frAC{\pArtial V}{\parTiaL z} \left(\delta z\RigHt)^2 +
2\frac{\PaRtiAL V}{\parTiaL z^*} |\delTa z|^2,\nonUMber \\
%
i \fRac{d}{d\tau} |\dElTA z|^2 &=&
-i\GamMA |\delta z|^2
-\FRAc{\PartIal V^*}{\partial z} \left(\dELtA Z\right)^2 +
\frac{\partIal V}{\paRtIAl Z^*} \LEft(\DelTa z^*\right)^2.\enD{aLigneD}$$ it is easY To SEE ThaT eqs. (\[A2\]) are equivaLent to Eqs. (\[maTRix\]) With $\epSiLon=0$, IF one maKes suBsTItuTions $|\delta z|^2\RighTarrow B$ anD $(\delta Z)^2\RightarROw C$.

It shOuld be NotIceD thaT YeT eXisTs ONe dIFfEreNCe bEtween thE lInEarizAtioN OF THe clAssIcal MotioN equations and The EquaTIonS for qUantuM cumUlAnts (\[9b\]), (\[9C\]): It is iMpossIbLe to get the initiAl coNditions (\[iNitIaL\_coNd\]) For $C$ aND $B$ from OnlY inItial coNditionS For ThE LINeArized classical equAtIONs Of motion (See alsO DiScUSsion of tHiS prObleM IN [@9; @9']).

 {#appeNdix2}

iN tHis AppenDix, we pREsEnT the derIvAtion oF tHe tIme Scale FOr thE validIty of the $1/n$-methOD for lossless ($\GaMMa=0$) oscillators. sTaRT WiTH the CasE of Kerr nonlIneaRIty $l=1$ And tHEn GenERalizE obtaInED rESults for arbitrary $l$. FRoM Eqs. (\[z\_cL\]) and (\[B\_And\_C\]) for $\Gamma=0$ aNd $l=1$, we get $$\laBEL{b1}
z_{cl}(\tau)=z_0\Exp(-i\oMeGA\tau), \quad
\Omega\eQuiv\bAr{\Delta}+|z_0|^2,$$ $$\laBEl{B2}
C(\tau)=-\tAu z_0^2 \leFt[ |z_
 V/\partial z\right)=
\lef t( \partia l V^*  /\ par ti al z ^* \ right)$, we ha v e fr om Eq. (\[A1\]) for th e qua dr a ticv ar iable s $(\de l ta z )^2 $an d $ |\ d el ta z| ^2$  the fo llowing eq uat io ns of motion $$ \begin{ali gne d}
\label{A2 }
i  \frac {d }{d \ tau}\ lef t(\de lta z\ r ight)^ 2 &=&
-i\ Ga m ma \le f t(\delt a z\ righ t)^2 +
2\frac{\pa r ti a l V}{\partialz} \le ft ( \d e l taz\r ight)^2 +2\ frac{ \ partial V} { \ p art i al z^*} |\del ta z|^2,\no n umb er \\%i \ f rac{d} {d\ta u} |\d elta z|^2 & =&
- i\Gamma | \delta z|^2
-\ f rac{\pa rtialV^* }{\ part i al  z } \ le f t(\ d el taz \ri ght)^2 + 
\ fr ac{\p arti a l V }{\p art ialz^*}\left(\deltaz^* \rig h t)^ 2.\en d{ali gned }$ $ Itis eas y tose e that Eqs. (\[ A2\] ) are equ iva le ntto  Eqs. (\[mat rix \])  with $ \epsilo n =0$ ,i f on e makes substituti on s $| \delta z |^2\ri g ht ar r ow B$ an d$(\ delt a z)^2\ righ t ar row C$.
It sh o ul dbe noti ce d that  y etexi sts o n e di fferen ce betwe en th e  linearization of the classi c al m ot i on e qua tions and t he e q uati onsf or  qu a ntumcumul an t s( \[9b\]), (\[9c\]):It  is im possi ble to get th e initialc o n ditions(\[i n it i al\_cond\]) fo r $C$  and $B$ f r om onlyiniti al condi tions for t he linea riz edcla ssi c a lequations ofm o tion  ( see als o d iscussi onofthi s p ro blem in [ @9; @9'] ). 
 { #a ppe ndix2 } 

In thi sApp en dix , wep resent  thederi va ti o n o f the t i me s cale  f or  the  va li dityof t h e $ 1/N$-me thod forlos s less  ( $\ Gamma=0 $) oscillator s.  Start wit hthe  caseo f  Kerr no nlinearity $l=1$ and th e n gener ali ze ob tain ed result s f or arb itr a ry $l$ . From  Eqs.  ( \[z \ _ cl\]) a nd  (\ [B \_and\_C\] ) for  $\Ga mm a=0$  and $l =1$, we get $$\lab e l{B 1}
z_{cl}(\ta u)= z_0\ e x p( -i\ O me g a\t au ) , \ q u ad
\Omega\equiv \bar{\Delt a} + |z _0|^2,$$ $ $ \la be l{B2}
C (\tau)= -\tau z_0^2 \ left[ |z_ 
 V/\partial_z\right)=
\left( \partial_V^* /\partial z^* \right)$,_we have_from_Eq. (\[A1\])_for_the quadratic variables_$(\delta z)^2$ and_$|\delta z|^2$ the following_equations of motion_$$\begin{aligned}
\label{A2}
i_\frac{d}{d\tau}\left(\delta z\right)^2 &=&
-i\Gamma \left(\delta z\right)^2 +
2\frac{\partial V}{\partial z} \left(\delta z\right)^2 +
2\frac{\partial V}{\partial z^*} |\delta_z|^2,\nonumber_\\
%
i \frac{d}{d\tau}_|\delta_z|^2_&=&
-i\Gamma |\delta z|^2
-\frac{\partial V^*}{\partial z}_\left(\delta z\right)^2 +
\frac{\partial V}{\partial z^*}_\left(\delta z^*\right)^2.\end{aligned}$$_It is easy to see that Eqs. (\[A2\])_are_equivalent to Eqs._(\[matrix\]) with $\epsilon=0$, if one makes substitutions $|\delta z|^2\rightarrow_B$ and $(\delta z)^2\rightarrow C$.

It should_be noticed that_yet_exists_one difference between the_linearization of the classical motion equations_and the equations for quantum cumulants_(\[9b\]), (\[9c\]): It is impossible to get_the initial conditions (\[initial\_cond\]) for $C$_and $B$ from only initial_conditions for_the linearized classical equations of_motion (see also_discussion of_this problem in_[@9; @9']).

 {#appendix2}

In this Appendix, we_present the derivation_of the time scale for the_validity_of the $1/N$-method_for_lossless_($\Gamma=0$) oscillators._Start with the_case_of Kerr_nonlinearity_$l=1$ and then generalize obtained results_for_arbitrary $l$. From Eqs. (\[z\_cl\]) and (\[B\_and\_C\])_for $\Gamma=0$ and $l=1$,_we_get $$\label{B1}
z_{cl}(\tau)=z_0\exp(-i\Omega\tau), \quad
\Omega\equiv\bar{\Delta}+|z_0|^2,$$ $$\label{B2}
C(\tau)=-\tau_z_0^2 \left[ |z_