#### 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 the -function, the average number of photons, and the distance between density matrices.

#### 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 the th-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 -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]In this equation, is the density matrix; the superoperator is (we have set )with being the interaction (time-independent) Hamiltonian; the superoperator , which describes the interaction of the system with the environment, is given bywith and 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 iswhere 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 :We simplify the above expression using the matrix method [23]; we define a triangular array of superoperators, where the diagonal elements are given by the nonperturbed system and the superior triangle contains the perturbation:Equation (5) can then be written aswhere stands for the element of the matrix . 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) aswhere perturbed density matrix is given byDeriving (7) and (8) with respect to time and equating terms, we obtain the differential equationor the equivalent systems of differential equationsThis system of equations can be solved and it is easy to see that the differential equationis 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,Doing the transformation , we haveand so

##### 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 from the Taylor series expansion (5); we have thenIn this case, we increase the dimension of our superoperators matrix as follows:In the previous section, it is demonstrated that the element of 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 of the new raised to the power ; indeed, we getFollowing the same steps that are in the first-order case, we can write the density matrix aswhere the solution to will be associated with the third column of the perturbed density matrixSolving the system of equations through the transformations and , we get the second-order correction

##### 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:and the th-order correction can be expressed as

#### 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]with , , and the parameter being the ratio between the cavity decay and the Kerr medium constant. The superoperators , , , and are defined asand they satisfy the commutation relationsAs initial state condition of the system, one can assume a coherent state; that is, ; so, with the help of the commutation relations, one getswhereand withThe 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 for th-order correction and taking into account that the superoperator will be defined by the sum of superoperators and ; taking this into consideration, we get the second-order correctionThe 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 matrixwhere

#### 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 Husimi -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 -function, we findwithFor the approximated density matrix, we getwithThe Husimi function, as it evolves in time, is plotted in Figure 1 for . 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 . 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 .

**(a)**

**(b)**

**(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 obtainand, with the second-order approximated solution, we get These results are presented in Figure 3. The parameters chosen are and , 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.

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 bywhere is a parameter that evaluates the closeness of and . Both matrices will be similar if or completely different if .

Figure 4 illustrates the numerical evaluation of parameter upon , for 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 -function and the average number of photons.

#### 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 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.

#### Competing Interests

The authors declare that they have no competing interests.

#### Acknowledgments

B. M. Villegas-Martínez acknowledges CONACYT for support.