# Application of Perturbation Theory to a Master Equation

Authors: B. M. Villegas-Martínez; F. Soto-Eguibar; H. M. Moya-Cessa
Journal: Advances in Mathematical Physics
 (2016)

---

## Abstract

We develop a matrix perturbation method for the Lindblad master equation. The first- and second-order corrections are obtained and the method is generalized for higher orders. The perturbation method developed is applied to the problem of a lossy cavity filled with a Kerr medium; the second-order corrections are estimated and compared with the known exact analytic solution. The comparison is done by calculating theQ-function, the average number of photons, and the distance between density matrices.


---

## Body

## 1. Introduction

Open systems, that is, systems that interact with an environment, represent an important problem in many branches of physics such as cosmology [1–3], quantum information [4], quantum optics [5], and condensed matter [6]. The irreversible loss of information and the dissipative process generated by a reservoir are a sign that the environment plays an important role in the dynamics of physical systems [5, 7, 8]. Open system dynamics are often determined by an equation of motion for the density matrix, the master equation. Usually these equations require certain assumptions to provide the dynamic description [9–12]. The Lindblad master equation represents one of the simplest cases that describe losses in an adequate form. Unfortunately, even this master equation is difficult to treat and to solve in an exact form [13–15]; hence, it is often required to apply a perturbative treatment. Many perturbative methods have been developed to solve particular problems modeled by the Lindblad master equation, such as a two-level nonlinear quantum system, a single-mode field in a lossy cavity, two-level atom coupling to a Bose-mode environment, and a single atom coupling to a mode of a lossy cavity [16–18]. In fact, it has been shown that even though decoherence takes place, the reconstruction of quasi-probability distribution functions may be achieved in atom-field [19, 20] or laser-trapped ion interactions [21].In this work, we show that it is possible to implement a matrix perturbation method on the Lindblad master equation that allows us to determine in a simple and effective form themth-order correction. The present paper is divided into two parts. In the first one, we develop a perturbative method for the Lindblad master equation by using superoperator techniques in terms of matrices; we obtain the first- and second-order corrections; these results lead to the generalization of the method to further order corrections. In the second part, the validity of the method is verified through a particular problem concerning a lossy cavity filled with a Kerr medium that has exact analytical solution. We obtain an approximate solution using the method proposed in this contribution and compare it with the exact result. The verification of both results is based on a comparative evaluation of the system via Q-function, the average photon number, and the measure of the distance between its density matrices.


## 2. Lindblad Master Equation

The Lindblad master equation, which describes the interaction between a given system and its environment at zero-temperature, is given by [5, 6, 22](1)
d
ρ
d
t
=
S
^
+
λ
L
^
ρ
.In this equation, ρ is the density matrix; the superoperator S
^ is (we have set ħ
=
1)(2)
S
^
ρ
=
-
i
H
,
ρ
,with H being the interaction (time-independent) Hamiltonian; the superoperator L
^, which describes the interaction of the system with the environment, is given by(3)
L
^
ρ
=
γ
2
a
ρ
a
†
-
a
†
a
ρ
-
ρ
a
†
a
,with a and a
† being the usual creation and annihilation operators and γ is the rate at which the system loses energy; and, finally, λ is a perturbation parameter (effectively, the perturbation parameter is λ
γ as we consider small γ’s).The formal solution to the master equation is(4)
ρ
t
=
e
S
^
+
λ
L
^
t
ρ
0
,where ρ
(
0
) is the density matrix of the initial state of the system.


### 2.1. First-Order Correction

To get the first-order correction to the nonperturbed solution to the master equation, we expand the exponential in (4) in Taylor series and keep only first-order terms in λ:(5)
ρ
t
≈
∑
n
=
0
∞
t
n
n
!
S
^
n
+
λ
∑
n
=
1
∞
t
n
n
!
∑
m
=
0
n
-
1
S
^
m
L
^
S
^
n
-
m
-
1
ρ
0
.We simplify the above expression using the matrix method [23]; we define a triangular 2
×
2 array of superoperators, where the diagonal elements are given by the nonperturbed system and the superior triangle contains the perturbation:(6)
M
ρ
=
S
^
L
^
0
S
^
ρ
.Equation (5) can then be written as(7)
ρ
t
≈
e
S
^
t
+
λ
e
M
t
1,2
ρ
0
,where (
e
M
t
)
1,2 stands for the element (
1,2
) of the matrix e
M
t. In what follows, we will denote the matrix elements of a matrix by a couple of integers subindexes. In this expression, we have separated the first-order approximation into two parts; one involves only the system and the other tells us how the environment affects the system. Thus, we can split the density matrix into one part concerning the nonperturbed system and a small contribution in terms of λ referent to first-order perturbation; so we can write (7) as(8)
ρ
t
≈
ρ
0
t
+
λ
ρ
P
1,2
,where 2
×
2 perturbed density matrix is given by(9)
ρ
P
=
ρ
1,1
ρ
1,2
ρ
2,1
ρ
2,2
.Deriving (7) and (8) with respect to time and equating λ terms, we obtain the differential equation(10)
d
d
t
ρ
P
1,2
=
M
e
M
t
ρ
0
1,2or the equivalent systems of differential equations(11)
d
d
t
ρ
1,1
ρ
1,2
ρ
2,1
ρ
2,2
=
M
e
M
t
ρ
0
0
0
ρ
0
.This system of equations can be solved and it is easy to see that the differential equation(12)
d
d
t
ρ
P
=
M
ρ
Pis also satisfied. By doing a matrix multiplication in the above expression, we can prove that the first-order solution is related to the second column of perturbed density matrix; that is,(13)
d
d
t
ρ
1,2
ρ
2,2
=
M
ρ
1,2
ρ
2,2
.Doing the transformation ϕ
1,2
=
e
-
S
^
t
ρ
1,2, we have(14)
d
d
t
ϕ
1,2
=
e
-
S
^
t
L
^
e
S
^
t
ρ
0
,and so(15)
ρ
1,2
=
e
S
^
t
∫
0
t
e
-
S
^
t
1
L
^
e
S
^
t
1
d
t
1
ρ
0
.


### 2.2. Second-Order Correction

The second-order correction to the nonperturbed solution to the master equation may be obtained if we take into account the terms inλ
2 from the Taylor series expansion (5); we have then(16)
ρ
t
≈
∑
n
=
0
∞
t
n
n
!
S
^
n
ρ
0
+
λ
∑
n
=
1
∞
t
n
n
!
∑
m
=
0
n
-
1
S
^
m
L
^
S
^
n
-
m
-
1
ρ
0
+
λ
2
∑
n
=
2
∞
t
n
n
!
∑
m
=
1
n
-
1
∑
j
=
0
n
-
m
-
1
S
^
n
-
m
-
j
-
1
L
^
S
^
j
L
^
S
^
m
-
1
ρ
0
.In this case, we increase the dimension of our superoperators matrix as follows:(17)
M
ρ
=
S
^
L
^
0
0
S
^
L
^
0
0
S
^
ρ
.In the previous section, it is demonstrated that the element 1,2 of M
n gave us all terms of the first-order correction. The second order will be a similar situation: all the information will be in the element 1,3 of the new M raised to the power n; indeed, we get(18)
ρ
t
≈
e
S
^
t
+
λ
e
M
t
1,2
+
λ
2
e
M
t
1,3
ρ
0
.Following the same steps that are in the first-order case, we can write the density matrix as(19)
ρ
t
≈
ρ
0
t
+
λ
ρ
P
1,2
+
λ
2
ρ
P
1,3
,where the solution to ρ
1,3 will be associated with the third column of the perturbed density matrix(20)
d
d
t
ρ
1,3
ρ
2,3
ρ
3,3
=
M
ρ
1,3
ρ
2,3
ρ
3,3
.Solving the system of equations through the transformations ϕ
1,3
=
e
-
S
^
t
ρ
1,3 and ϕ
2,3
=
e
-
S
^
t
ρ
2,3, we get the second-order correction(21)
ρ
1,3
=
e
S
^
t
∫
0
t
∫
0
t
1
e
-
S
^
t
1
L
^
e
S
^
t
1
e
-
S
^
t
2
L
^
e
S
^
t
2
d
t
2
d
t
1
ρ
0
.


### 2.3. Higher Orders

The generalization of the method for higher-order corrections can be obtained directly from the results of the first-order and second-order corrections. So, following the same steps that take us to expression (21), we define the semi-infinite superoperators array:(22)
M
=
S
^
L
^
⋯
0
0
0
S
^
L
^
⋯
0
⋮
0
S
^
⋱
⋮
⋮
0
0
⋱
L
^
0
⋮
0
0
S
^
,and the mth-order correction can be expressed as(23)
ρ
1
,
m
+
1
=
e
S
^
t
∫
0
t
∫
0
t
1
⋯
∫
0
t
m
-
1
e
-
S
^
t
1
L
^
e
S
^
t
1
e
-
S
^
t
2
L
^
e
S
^
t
2
⋯
-
e
-
S
^
t
m
L
^
e
S
^
t
m
d
t
m
⋯
d
t
2
d
t
1
ρ
0
.


## 2.1. First-Order Correction

To get the first-order correction to the nonperturbed solution to the master equation, we expand the exponential in (4) in Taylor series and keep only first-order terms in λ:(5)
ρ
t
≈
∑
n
=
0
∞
t
n
n
!
S
^
n
+
λ
∑
n
=
1
∞
t
n
n
!
∑
m
=
0
n
-
1
S
^
m
L
^
S
^
n
-
m
-
1
ρ
0
.We simplify the above expression using the matrix method [23]; we define a triangular 2
×
2 array of superoperators, where the diagonal elements are given by the nonperturbed system and the superior triangle contains the perturbation:(6)
M
ρ
=
S
^
L
^
0
S
^
ρ
.Equation (5) can then be written as(7)
ρ
t
≈
e
S
^
t
+
λ
e
M
t
1,2
ρ
0
,where (
e
M
t
)
1,2 stands for the element (
1,2
) of the matrix e
M
t. In what follows, we will denote the matrix elements of a matrix by a couple of integers subindexes. In this expression, we have separated the first-order approximation into two parts; one involves only the system and the other tells us how the environment affects the system. Thus, we can split the density matrix into one part concerning the nonperturbed system and a small contribution in terms of λ referent to first-order perturbation; so we can write (7) as(8)
ρ
t
≈
ρ
0
t
+
λ
ρ
P
1,2
,where 2
×
2 perturbed density matrix is given by(9)
ρ
P
=
ρ
1,1
ρ
1,2
ρ
2,1
ρ
2,2
.Deriving (7) and (8) with respect to time and equating λ terms, we obtain the differential equation(10)
d
d
t
ρ
P
1,2
=
M
e
M
t
ρ
0
1,2or the equivalent systems of differential equations(11)
d
d
t
ρ
1,1
ρ
1,2
ρ
2,1
ρ
2,2
=
M
e
M
t
ρ
0
0
0
ρ
0
.This system of equations can be solved and it is easy to see that the differential equation(12)
d
d
t
ρ
P
=
M
ρ
Pis also satisfied. By doing a matrix multiplication in the above expression, we can prove that the first-order solution is related to the second column of perturbed density matrix; that is,(13)
d
d
t
ρ
1,2
ρ
2,2
=
M
ρ
1,2
ρ
2,2
.Doing the transformation ϕ
1,2
=
e
-
S
^
t
ρ
1,2, we have(14)
d
d
t
ϕ
1,2
=
e
-
S
^
t
L
^
e
S
^
t
ρ
0
,and so(15)
ρ
1,2
=
e
S
^
t
∫
0
t
e
-
S
^
t
1
L
^
e
S
^
t
1
d
t
1
ρ
0
.


## 2.2. Second-Order Correction

The second-order correction to the nonperturbed solution to the master equation may be obtained if we take into account the terms inλ
2 from the Taylor series expansion (5); we have then(16)
ρ
t
≈
∑
n
=
0
∞
t
n
n
!
S
^
n
ρ
0
+
λ
∑
n
=
1
∞
t
n
n
!
∑
m
=
0
n
-
1
S
^
m
L
^
S
^
n
-
m
-
1
ρ
0
+
λ
2
∑
n
=
2
∞
t
n
n
!
∑
m
=
1
n
-
1
∑
j
=
0
n
-
m
-
1
S
^
n
-
m
-
j
-
1
L
^
S
^
j
L
^
S
^
m
-
1
ρ
0
.In this case, we increase the dimension of our superoperators matrix as follows:(17)
M
ρ
=
S
^
L
^
0
0
S
^
L
^
0
0
S
^
ρ
.In the previous section, it is demonstrated that the element 1,2 of M
n gave us all terms of the first-order correction. The second order will be a similar situation: all the information will be in the element 1,3 of the new M raised to the power n; indeed, we get(18)
ρ
t
≈
e
S
^
t
+
λ
e
M
t
1,2
+
λ
2
e
M
t
1,3
ρ
0
.Following the same steps that are in the first-order case, we can write the density matrix as(19)
ρ
t
≈
ρ
0
t
+
λ
ρ
P
1,2
+
λ
2
ρ
P
1,3
,where the solution to ρ
1,3 will be associated with the third column of the perturbed density matrix(20)
d
d
t
ρ
1,3
ρ
2,3
ρ
3,3
=
M
ρ
1,3
ρ
2,3
ρ
3,3
.Solving the system of equations through the transformations ϕ
1,3
=
e
-
S
^
t
ρ
1,3 and ϕ
2,3
=
e
-
S
^
t
ρ
2,3, we get the second-order correction(21)
ρ
1,3
=
e
S
^
t
∫
0
t
∫
0
t
1
e
-
S
^
t
1
L
^
e
S
^
t
1
e
-
S
^
t
2
L
^
e
S
^
t
2
d
t
2
d
t
1
ρ
0
.


## 2.3. Higher Orders

The generalization of the method for higher-order corrections can be obtained directly from the results of the first-order and second-order corrections. So, following the same steps that take us to expression (21), we define the semi-infinite superoperators array:(22)
M
=
S
^
L
^
⋯
0
0
0
S
^
L
^
⋯
0
⋮
0
S
^
⋱
⋮
⋮
0
0
⋱
L
^
0
⋮
0
0
S
^
,and the mth-order correction can be expressed as(23)
ρ
1
,
m
+
1
=
e
S
^
t
∫
0
t
∫
0
t
1
⋯
∫
0
t
m
-
1
e
-
S
^
t
1
L
^
e
S
^
t
1
e
-
S
^
t
2
L
^
e
S
^
t
2
⋯
-
e
-
S
^
t
m
L
^
e
S
^
t
m
d
t
m
⋯
d
t
2
d
t
1
ρ
0
.


## 3. Lossy Cavity Filled with a Kerr Medium

To demonstrate the accuracy and capability of the method, we obtain the perturbative solution to the master equation of a Kerr medium filling an optical cavity with losses. The exact analytic solution for the master equation in this case is [24](24)
ρ
e
x
a
τ
=
e
S
^
τ
e
L
^
τ
exp
⁡
1
-
e
-
2
τ
λ
+
i
R
^
2
λ
+
i
R
^
J
^
ρ
0
,with τ
=
χ
t, λ
=
γ
/
χ, and the parameter λ being the ratio between the cavity decay and the Kerr medium constant. The superoperators S
^, J
^, L
^, and R
^ are defined as(25)
S
^
ρ
=
-
i
a
†
2
a
2
,
ρ
,
J
^
ρ
=
2
λ
a
ρ
a
†
,
L
^
ρ
=
-
λ
a
†
a
ρ
+
ρ
a
†
a
,
R
^
ρ
=
a
†
a
ρ
-
ρ
a
†
a
,and they satisfy the commutation relations(26)
S
^
,
J
^
ρ
=
2
i
R
^
J
^
ρ
,
L
^
,
J
^
ρ
=
2
λ
J
^
ρ
,
R
^
,
J
^
ρ
=
0
.As initial state condition of the system, one can assume a coherent state; that is, ρ
(
0
)
=
α
α; so, with the help of the commutation relations, one gets(27)
ρ
e
x
a
τ
=
e
-
α
2
∑
n
=
0
∞
∑
m
=
0
∞
C
n
,
m
n
m
,where(28)
C
n
,
m
=
α
n
α
∗
m
n
!
m
!
e
-
i
n
n
-
1
-
m
m
-
1
τ
D
n
,
mand with(29)
D
n
,
m
=
e
-
λ
n
+
m
τ
exp
⁡
λ
α
2
1
-
e
-
2
τ
λ
+
i
n
-
m
λ
+
i
n
-
m
.The exact density matrix shows how the initial coherent state structure is lost due to dissipation of energy generated by the cavity walls and the quadratic terms associated with the nonlinear medium.


### 3.1. Second-Order Correction

The approximate solution for the cavity problem is found using the expression formth-order correction and taking into account that the superoperator L
^
ρ will be defined by the sum of superoperators J
^
ρ and L
^
ρ; taking this into consideration, we get the second-order correction(30)
ρ
t
≈
e
S
^
τ
ρ
0
+
λ
e
S
^
τ
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
+
L
^
τ
ρ
0
+
λ
2
e
S
^
τ
L
^
2
τ
2
2
!
ρ
0
+
λ
2
e
S
^
τ
1
2
!
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
2
+
2
2
i
R
^
τ
e
-
2
i
R
^
τ
-
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
ρ
0
+
λ
2
e
S
^
τ
L
^
τ
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
ρ
0
.The difference between (24) and (30) is that for the approximated solution λ has been considered as a perturbation parameter. Considering as initial state a coherent state, we arrived at the approximate density matrix(31)
ρ
a
p
r
o
τ
=
e
-
α
2
∑
p
=
0
∞
∑
q
=
0
∞
F
p
,
q
p
q
,where(32)
F
p
,
q
=
α
p
α
∗
q
p
!
q
!
e
-
i
p
p
-
1
-
q
q
-
1
τ
1
+
H
p
,
q
,
H
p
,
q
=
λ
2
τ
2
2
p
+
q
-
2
α
2
e
-
i
p
-
q
τ
s
i
n
c
p
-
q
τ
2
-
λ
p
+
q
τ
+
2
λ
α
2
τ
e
-
i
p
-
q
τ
s
i
n
c
p
-
q
τ
-
λ
2
α
2
e
-
2
i
p
-
q
τ
1
-
e
2
i
p
-
q
τ
+
2
i
p
-
q
τ
p
-
q
2
.


## 3.1. Second-Order Correction

The approximate solution for the cavity problem is found using the expression formth-order correction and taking into account that the superoperator L
^
ρ will be defined by the sum of superoperators J
^
ρ and L
^
ρ; taking this into consideration, we get the second-order correction(30)
ρ
t
≈
e
S
^
τ
ρ
0
+
λ
e
S
^
τ
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
+
L
^
τ
ρ
0
+
λ
2
e
S
^
τ
L
^
2
τ
2
2
!
ρ
0
+
λ
2
e
S
^
τ
1
2
!
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
2
+
2
2
i
R
^
τ
e
-
2
i
R
^
τ
-
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
ρ
0
+
λ
2
e
S
^
τ
L
^
τ
1
-
e
-
2
i
R
^
τ
2
i
R
^
J
^
ρ
0
.The difference between (24) and (30) is that for the approximated solution λ has been considered as a perturbation parameter. Considering as initial state a coherent state, we arrived at the approximate density matrix(31)
ρ
a
p
r
o
τ
=
e
-
α
2
∑
p
=
0
∞
∑
q
=
0
∞
F
p
,
q
p
q
,where(32)
F
p
,
q
=
α
p
α
∗
q
p
!
q
!
e
-
i
p
p
-
1
-
q
q
-
1
τ
1
+
H
p
,
q
,
H
p
,
q
=
λ
2
τ
2
2
p
+
q
-
2
α
2
e
-
i
p
-
q
τ
s
i
n
c
p
-
q
τ
2
-
λ
p
+
q
τ
+
2
λ
α
2
τ
e
-
i
p
-
q
τ
s
i
n
c
p
-
q
τ
-
λ
2
α
2
e
-
2
i
p
-
q
τ
1
-
e
2
i
p
-
q
τ
+
2
i
p
-
q
τ
p
-
q
2
.


## 4. Comparison between Exact and Approximate Solution

A simple and direct form to visualize the evolution of a cavity-Kerr system in phase space is calculating a quasi-probability function. The HusimiQ-function is the simplest of all quasi-probability functions and it is defined as the expectation value of the density matrix in a coherent base [25, 26]. Thus, if we use the exact density matrix to evaluate the Q-function, we find(33)
Q
e
x
a
β
=
1
π
e
-
α
2
+
β
2
∑
n
=
0
∞
∑
l
=
0
∞
α
β
∗
n
α
∗
β
l
n
!
l
!
e
-
λ
n
+
l
τ
D
n
,
l
,with(34)
D
n
,
l
=
e
-
i
n
n
-
1
-
m
m
-
1
τ
exp
⁡
λ
α
2
1
-
e
-
2
τ
λ
+
i
n
-
m
λ
+
i
n
-
m
.For the approximated density matrix, we get(35)
Q
a
p
r
o
β
=
1
π
e
-
α
2
+
β
2
∑
p
=
0
∞
∑
l
=
0
∞
α
β
∗
p
α
∗
β
l
p
!
l
!
e
-
i
p
p
-
1
-
l
l
-
1
τ
1
+
G
p
,
l
,with(36)
G
p
,
l
=
λ
2
τ
2
2
p
+
l
-
2
α
2
e
-
i
p
-
l
τ
s
i
n
c
p
-
l
τ
2
-
λ
p
+
l
τ
+
2
λ
α
2
τ
e
-
i
p
-
l
τ
s
i
n
c
p
-
l
τ
+
4
τ
λ
2
α
2
e
-
i
p
-
l
τ
e
-
i
p
-
l
τ
-
s
i
n
c
p
-
l
τ
2
i
p
-
l
.The Husimi function, as it evolves in time, is plotted in Figure 1 for λ
=
0.05. In Figure 1(a), we have the numerical results for the exact solution, and, in Figure 1(b), we have the results for the second-order correction solution. It seems that in short times the result is the same and for larger times some differences appear; in Figure 2, we show the evolution of the Husimi function for λ
=
0.13. Apparently, if we increase the value of lambda, the range of time in which the exact and the approximated solutions are similar becomes smaller. The reason for this behavior is easy to understand, when it is realized that the real perturbation parameter is λ
τ and not just l
a
m
b
d
a.Figure 1
Evolution of theQ-function with α
=
2 and λ
=
0.05. (a) Exact solution and (b) approximated solution.
(a)
(b)Figure 2
Evolution of theQ-function with α
=
2 and λ
=
0.13. (a) Exact solution and (b) approximated solution.
(a)
(b)As a second way to test the accuracy of the perturbation approximation, we proceed to calculate the mean photon number, which is a relevant physical quantity of the Kerr lossy cavity. Using the exact solution, we obtain(37)
N
e
x
a
=
α
2
e
-
2
λ
τ
,and, with the second-order approximated solution, we get (38)
N
a
p
r
o
=
α
2
1
-
2
λ
τ
+
-
2
λ
τ
2
2
!
.These results are presented in Figure 3. The parameters chosen are α
=
2 and λ
=
0.05, 0.07, 0.09, 0.11, and 0.13. The solid lines represent the exact solution, whereas the dotted lines show the results of the perturbative solution.Figure 3
The mean photon numberN versus time τ, with α
=
2 and λ
= 0.05, 0.07, 0.09, 0.11, and 0.13. The solid lines are for the exact solution, while the dotted lines are for the perturbative solution.The figure shows that approximate solutions with different values ofλ are indeed a good approximation for the problem of Kerr lossy cavity for short times.Finally, as another measure of proximity for the solutions, we evaluate the distance between the exact density matrix and the approximated density matrix [27]. The geometrical measure of distance between two density matrixes is given by(39)
ρ
1
-
ρ
2
2
=
1
2
T
r
ρ
1
-
ρ
2
2
=
1
-
F
,where F
=
T
r
ρ
1
ρ
2
+
(
1
/
2
)
2
-
T
r
ρ
1
2
-
T
r
ρ
2
2 is a parameter that evaluates the closeness of ρ
1 and ρ
2. Both matrices will be similar if F
=
1 or completely different if F
=
0.Figure4 illustrates the numerical evaluation of parameter F upon τ, for α
=
2 and different values of λ. This plot shows that density matrices are very similar for short times; these results coincide and are in agreement with the temporal behavior of the Q-function and the average number of photons.Figure 4
Plot ofF
=
T
r
ρ
1
ρ
2
+
(
1
/
2
)
2
-
T
r
ρ
1
2
-
T
r
ρ
2
2 versus τ for different values of λ.


## 5. Conclusions

In summary, we can conclude that in the examined case, a lossy cavity filled with a Kerr medium, the matrix perturbative method gives good results. When time grows the results start to differ, but that is not surprising since the real measure of the perturbation is not given by onlyλ, and it is given by the product λ
τ. The same logical behavior is observed in the case of the Husimi function and the parameter F for different values of lambda, where we found good agreement for short times and some differences when the time grows. The second-order result is enough to precisely reproduce the exact solution for this specific system; when studying more complicated systems, the contributions of orders superior to the first one could be relevant.



---

*Source: 9265039-2016-06-30.xml*