О допустимости некорректного выбора объема квантования

Несмотря на, казалось бы, исчерпывающее описание процессов параметрического рассеяния света в литературе ( и цит. лит.), представляет интерес вновь возвратиться к анализу этого явления и на этом примере проследить особенности выбора объема квантования. В настоящей работе поставленная задача решается двумя способами. В первом случае рассматривается некорректный, но более удобный с  практической точки зрения выбор объема квантования. Во втором случае используется более точное решение, справедливое для любого объема квантования. При сравнении, как ни странно, оказывается, что результаты получаются почти одинаковые. Это позволяет надеяться на то, что и в других задачах условия выбора объема квантования могут быть более мягкими.

В процессе параметрического рассеяния света фотоны накачки распадаются на пары сигнальных ($s$) и холостых ($i$) фотонов с частотами

Кроме того, существует жесткая связь между поляризациями сигнального и холостого фотонов.

Исчезновение фотона накачки описывается оператором уничтожения $\partial p,$ а рождение вместо него сигнального и холостого фотонов --- операторами рождения $\widehat a_s^{+}$ и $\widehat b_i^{+}.$ Это поясняет структуру трехмодового гамильтониана взаимодействия процесса

noindent параметрического рассеяния ( и цит. лит.)

где $\chi^{(2)}$ --- квадратичная нелинейность, а эрмитово сопряженный оператор H.c. описывает обратный параметрическому рассеянию процесс --- рождение фотона накачки при одновременном исчезновении сигнального и холостого фотонов, который также возможен.

Поскольку эффективность параметрического рассеяния мала, истощением накачки в первом приближении можно пренебречь, считая ее амплитуду постоянной. С учетом того, что для реального наблюдения параметрического рассеяния фотонов в накачке должно быть много, ее можно описывать классически, заменив оператор $\widehat a_p$ в гамильтониане на постоянную комплексную амплитуду $\varepsilon_p$:

где

отрицательно--частотные операторы поля, эрмитово положительно операторам поля $\widehat{E}^{(+)}_{s,i}(\vec{r},t).$ Частоты $\omega_s=\frac{ck_s}{n_s}$, $\omega_i=\frac{ck_i}{n_i}$ меньше частоты накачки $\omega_p$; $k_{s,i}=|\vec{k}_{s,i}|;$ $n_s$ и $n_i$ --- показатели преломления нелинейной среды; $V$ --- объем нелинейной среды; $\varepsilon_p(\vec{r},t)$ --- классическая комплексная амплитуда накачки, которую далее полагаем плоской и монохроматической, распространяющейся по направлению оси $z$: $\varepsilon_p(\vec{r},t)=\varepsilon_pe^{i(k_pz-\omega_pt)},$ $k_p=|\vec{k}_p|.$

Объем квантования зададим равным объему нелинейной среды $V$, полагая, что она имеет форму прямоугольного параллелепипеда, хотя формально объем квантования должен быть существенно больше объема среды, так как на границах объема квантования амплитуды полей должны либо полностью спадать, либо не нарушать периодической структуры. Оси системы координат направим параллельно осям параллелепипеда, а ее центр расположим в центре параллелепипеда.

Динамика системы в представлении Гейзенберга описывается простыми соотношениями (см., например, ):

где $C_{k_sk_i}=\frac{\chi^{(2)}\varepsilon_p\sqrt{\omega_s\omega_i}}{i\cdot4\pi}.$

Это выражение справедливо, когда сигнальное ($s$) и холостое ($i$) излучения имеют одинаковые поляризации (вырождены по поляризации). Первое слагаемое в нем соответствует вырождению ${\vec{k}_s=\vec{k}_i}.$ Это слагаемое равно нулю в случае взаимноортогональных поляризаций сигнального и холостого излучений, который и будет рассматриваться далее.

Интегрирование по $\textbf{r}$ дает

Интервалы между соседними модами $(k_{xs}+k_{xi})_\textrm{min}=\frac{2\pi}{L_x}$ и аналогично для остальных компонент. Пусть $\widetilde{{\bf k}}_{xs}+\widetilde{{\bf k}}_{xi}=0,$ тогда $\mathop{\rm sinc}\nolimits\left[(k_{xs}+k_{xi})\frac{L_x}{2}\right]=1.$ При $\widetilde{{\bf k}}_{xs}\ne k_{xs}$ и (или) $\widetilde{{\bf k}}_{xi}\ne k_{xi}$ имеем $(k_{xs}+k_{xi})\frac{L_x}{2}=\pi m,$ $m={\pm}1,{\pm}2,\dots,$ и $\mathop{\rm sinc}\nolimits\left[(k_{xs}+k_{xi})\frac{L_x}{2}\right]=0.$

Итак,

Решение системы $\eqref{beleq11}$ $\eqref{beleq12}$имеет вид

Подставляя эти решения в выражения для отрицательно операторов поля, получим операторы поля в кристалле в виде суперпозиции мод, удовлетворяющих условию

Рассмотрим вырожденный по частоте режим: $\widetilde{\omega}_s=\widetilde{\omega}_i=\omega_p,$ тогда $\theta_s=\theta_i=\theta,$ следовательно, $\widetilde{k}_s=\widetilde{k}_i=k,$ и

где $v$ --- скорость света в среде.

Полагая $z=vt+z_0=\frac{vt}{2},$ имеем

и аналогично для $\widehat{b}_{k_\bot}$ и соответствующих эрмитово операторов.

Заметим, что $e^{i\vec{k}\vec{r}}=e^{i\bigl(k+k^2_\bot/(2k)\bigr)L_z},$ поскольку

Сравним теперь полученный результат с решением той же задачи также в приближении заданной накачки, но в представлении Гейзенберга. При этом будем полагать, что параметрическое рассеяние описывается квазилинейными уравнениями параболического типа (см., например, ):

где $\beta$ --- комплексный коэффициент нелинейной связи, зависящий от амплитуды, фазы и профиля накачки, а также от вида и геометрии взаимодействия, $\widehat{E}^{(+)}$ и $\widehat{E}^{(-)}$ --- операторы положительно- и отрицательно частей поля в представлении Гейзенберга.

Решение предыдущего уравнения можно записать для фурье $\widehat{E}^{(+)}_j(\mbox{$\varkappa\hspace{-6.7pt}\varkappa$},z)$ оператора $\widetilde{E}^{(+)}_j(\vec{r}_\bot,z)$ в виде

где $\mbox{$\varkappa\hspace{-6.7pt}\varkappa$}$ --- переменная фурье-преобразования, $\alpha=2\frac{\gamma}{v},$ $\theta=\frac{k^2_\bot}{k},$ $\beta=\frac{\chi^{(2)}\omega\varepsilon_p}{i\cdot4\pi v},$ $\varphi$ --- аргумент $\beta$, т. е. фаза накачки.

С точностью до фазового сомножителя $e^{i(k^2_\bot/2k)L_z},$ связанного с дополнительным фазовым набегом наклонных пучков, решение, полученное в первом случае (21), совпадает с решением, полученном во втором (23). При описании прямого детектирования параметрически рассеянного света вблизи нелинейного кристалла наличие этого сомножителя никак не повлияет на результат, поскольку квадратичный детектор фотонов реагирует на не зависящий от фазы квадрат модуля комплексной амплитуды поля или оператор числа фотонов в квантовом описании. Однако при удалении детекторов от параметрического генератора ситуация изменится, ведь фазовый сомножитель $e^{i(k^2_\bot/2k)L_z}$ фактически эквивалентен дополнительной линзе в параксиальном приближении, т. е. ход лучей будет иным. Более того, существенным будет искажение и в случае приготовления сжатых состояний света посредством параметрического рассеяния, поскольку этот процесс очень чувствителен в фазе взаимодействующих мод (см., например, ), а также в интерференционных схемах. Поэтому рассмотренное приближение следует применять с осторожностью, зная его некорректность в отношении этого фазового сомножителя. Вообще же говоря, коррекция фазы не представляет сложности введением дополнительных линз. Например, повышение эффективности квантового сжатия в расходящихся пучках легко осуществляется собирающей линзой с фокусным расстоянием $f',$ значительно превышающим линейные размеры кристалла $L_x,L_y,L_z.$ Оптическая ось линзы при этом должна совпадать с осью $z$ (см., например, ).

Таким образом, задача многомодового параметрического рассеяния света решена в двух вариантах: минимального (некорректного) объема квантования в представлении Шрёдингера и более корректного, не зависящего выбора определенного объема квантования, в представлении Гейзенберга; получены почти одинаковые результаты. Это позволяет надеяться на то, что, хотя бы в первом приближении, можно решать квантовые задачи несколько проще, чем это требуется при абсолютно правильно выбранном объеме квантования.

Итак, на основании исследования особенностей выбора объема квантования в процессе многомодового параметрического рассеяния света можно сделать вывод о допустимости в первом приближении использовать ради простоты не совсем корректный объем, меньший традиционно применяемого, но существенно упрощающего вычисления. Хотя это и выходит за рамки настоящей работы, отметим, что такое упрощение возможно окажется полезным при выборе объема квантования и в случае существенно неоднородных сред с заведомо малыми размерами: фотонные кристаллы, гетероструктуры, одиночные квантовые точки, квантовые ямы и др. (см., например, и цит. лит.). Именно эти среды составляют основу приоритетных исследований в настоящее время.