Рассмотрена система, состоящая из двух идентичных искусственных атомов (кубитов), нерезонансно взаимодействующих посредством вырожденных двухфотонных переходов с модой теплового квантового поля идеального микроволнового резонатора при наличии керровской нелинейности. Для рассматриваемой модели получено точное решение квантового уравнения Лиувилля для полной матрицы плотности системы «два кубита + мода поля резонатора». Для решения квантового уравнения эволюции использовано представление «одетых» состояний, то есть собственных функций гамильтониана. Найден полный набор «одетых» состояний рассматриваемой модели. С его помощью первоначально найдено решение уравнения эволюции для перепутанных начальных состояний кубитов и фоковских состояний поля резонатора, то есть состояний с определенным числом фотонов в резонаторной моде. Указанное решение использовано для построения точного решения квантового уравнения Лиувилля в случае теплового состояния поля резонатора. Усреднением полной матрицы плотности по переменным поля резонатора найдена редуцированная матрица плотности двух кубитов. Двухкубитная матрица плотности использована для вычисления параметра перепутывания кубитов в аналитическом виде для двух типов начальных перепутанных состояний кубитов белловского типа. В качестве количественного критерия перепутывания кубитов выбран параметр Переса-Хородецких, или отрицательность. Проведено численное моделирование временной зависимости параметра перепутывания кубитов для различных параметров модели и начальных состояний кубитов. Наиболее интересным представляется результат, заключающийся в том, что для некоторых параметров модели учет керровской нелинейности приводит к существенной стабилизации начального перепутывания кубитов, а также к исчезновению эффекта мгновенной смерти перепутывания.
Исследуется динамика трехкубитной многофотонной модели Тависа-Каммингса. Модель включает изолированный кубит и два кубита, захваченных в независимых одномодовых идеальных резонаторах. Захваченные кубиты резонансно взаимодействуют с тепловым или вакуумным полем своего резонатора посредством многофотонных переходов, а их начальные состояния соответствуют истинно перепутанным состояниям Гринбергера-Хорна-Цайлингера и Вернера. Для модели найдено точное решение в виде оператора эволюции, полученного из нестационарного уравнения Шредингера. С его помощью построено решение уравнения Лиувилля для матрицы плотности и вычислены меры перепутывания - отрицательность и верность. Результаты численного моделирования демонстрируют существенное различие в зависимости от типа поля. Для вакуумных полей увеличение многофотонного порядка
m не стабилизирует перепутывание, а лишь сокращает период колебаний отрицательности и верности. В случае же тепловых полей рост
m повышает максимальную верность и эффективно стабилизирует перепутывание кубитов. Таким образом, многофотонный порядок служит эффективным механизмом управления перепутанностью кубитов именно в условиях тепловых полей в резонаторах.