]> raman_wigner_diff Joel Corney Fibre soliton (sagnac configuation) squeezing calculation in Wigner representation. Raman noise modelled by including a heap (10) of cross-propagating phonon modes. This is a really cumbersome thing to do at the moment in xmds. It would be a lot easier if one could simulate integro-differential equations in xmds (with nonlocal terms), or if arrays of fields could be defined. Also there is something funny going on with the initialization of the cross-propagating modes - it changes with the longitudinal (forward-propagating) coordinate, if the cross-propagating mode has some initial dynamics. x yes yes yes 10000 10 20 40 /***********/ /***********/ /*************/ /*************/ main t &Nt; (-&tmax;,&tmax;) 1 1 1 1 /***MAIN***/ main complex phi dphi psi dpsi no &tmax; /***CROSS***/ /*Initialise the phonon variables */ cross complex b0 b1 b2 b3 b4 b5 b6 b7 b8 b9 db0 db1 db2 db3 db4 db5 db6 db7 db8 db9 b_lo0 b_lo1 b_lo2 b_lo3 b_lo4 b_lo5 b_lo6 b_lo7 b_lo8 b_lo9 db_lo0 db_lo1 db_lo2 db_lo3 db_lo4 db_lo5 db_lo6 db_lo7 db_lo8 db_lo9 /***HEAVISIDE***/ /*Step function*/ Heaviside double step_function &tmax; /***PHONONS***/ Phonons /*Constant operator for phonon free evolution*/ complex expikt0 expikt1 expikt2 expikt3 expikt4 expikt5 expikt6 expikt7 expikt8 expikt9 /******************/ /******************/ SIIP 12.57 500 25 100 25 100 yes L 4 main cross Heaviside Phonons /*****************/ /*Integrate phonon variables*/ /*****************/ cross t /**************/

/**************/ raman_wigner_diff.xsil no 32 main cross I I_phi I_psi bm0 bm1 bd0 bd1 b_lom0 b_lom1 b_lod0 b_lod1 no 0 N_sym,dN_sym complex p = sqrt(alpha)*phi-sqrt(1.0-alpha)*psi; /* Beam splitter*/ complex dp = sqrt(alpha)*dphi-sqrt(1.0-alpha)*dpsi; /* Beam splitter*/ N_sym = ~p*p; /*Symmetrically ordered*/ dN_sym = ~dp*dp+2*real(~p*dp); /*Symmetrically ordered*/ no N NN N0 double N_s = N_sym; double dN_s = dN_sym; double c1 = &Nt;/n; /*&Nt;Int[delta(0)/n dxi], from same-point commutator */ double c2 = 1/n; /*Int[delta(xi)/n dxi], from different-point commutator*/ N0 = N_s; /*mean number in output port*/ N= dN_s- c1/2; /*Normally ordered fluctuations*/ NN = dN_s*dN_s - c1*dN_s + (c1*c1-c2*c1)/4; /*Normally ordered fluctuations squared*/ yes 256 main cross Phonons I_phik I_psik db1k db9k b1k b9k exp1 exp9 yes 0 main cross norm_phi norm_psi omega_phi omega_psi