Geometric Phase of a Transmon in a Dissipative Quantum Circuit

Superconducting circuits reveal themselves as promising physical devices with multiple uses. Within those uses, the fundamental concept of the geometric phase accumulated by the state of a system shows up recurrently, as, for example, in the construction of geometric gates. Given this framework, we study the geometric phases acquired by a paradigmatic setup: a transmon coupled to a superconductor resonating cavity. We do so both for the case in which the evolution is unitary and when it is subjected to dissipative effects. These models offer a comprehensive quantum description of an anharmonic system interacting with a single mode of the electromagnetic field within a perfect or dissipative cavity, respectively. In the dissipative model, the non-unitary effects arise from dephasing, relaxation, and decay of the transmon coupled to its environment. Our approach enables a comparison of the geometric phases obtained in these models, leading to a thorough understanding of the corrections introduced by the presence of the environment.

Superconducting circuits reveal themselves as promising physical devices with multiple uses.Within those uses, the fundamental concept of the geometric phase accumulated by the state of a system shows up recurrently, as, for example, in the construction of geometric gates.Given this framework, we study the geometric phases acquired by a paradigmatic setup: a transmon coupled to a superconductor resonating cavity.We do so both for the case in which the evolution is unitary and when it is subjected to dissipative effects.These models offer a comprehensive quantum description of an anharmonic system interacting with a single mode of the electromagnetic field within a perfect or dissipative cavity, respectively.In the dissipative model, the non-unitary effects arise from dephasing, relaxation, and decay of the transmon coupled to its environment.Our approach enables a comparison of the geometric phases obtained in these models, leading to a thorough understanding of the corrections introduced by the presence of the environment.

I. INTRODUCTION
Significant advancements in coherent superconducting circuits have enabled the development of a diverse range of qubit designs, encompassing the classic flux [1][2][3][4], phase [5][6][7][8], and charge [9][10][11] qubits, as well as more contemporary designs like transmon [12] and fluxonium [13] circuits.The mathematical representation and quantum dynamics of these qubits, when transversely coupled to resonators, are described by circuit quantum electrodynamics (cQED) [14][15][16].Although these qubits are designed to behave as a two-level system within a superconducting circuit, they inherently possess multiple additional energy levels that can influence their interactions with other components.Focusing on the transmon qubit, one of its main advantages is its longer coherence timescale when compared to other circuits, which is essential for performing quantum computations and quantum error correction operations [17][18][19].This increased coherence can be attributed to the non-harmonic energy level structure, which helps suppress certain types of decohering effects [20].
Within the landscape of physical systems provided by cQED setups, the objects known as geometric phases (GPs) have lately played an important role, mainly for implementing measurements and other quantum operations and therefore allowing GPs to be harnessed for various applications, such as geometric gates.
The idea that the phase acquired by the state of a quantum system can be decomposed into a dynamical and a geometrical component originated with Berry's theoretical work, where it was constrained to the context of adiabatic, cyclic, unitary evolution [21].Subsequently, the concept of GP has been extended to nonadiabatic cyclic, non-cyclic, and even non-unitary evolutions [22][23][24][25][26][27][28][29].These generalizations consistently reduce to less comprehensive results as additional conditions are met.The GP has also been elucidated as a consequence of quantum kinematics, interpreted in terms of a parallel transport condition dependent solely on the geometry of the Hilbert space, from which its name is derived [30].
As evident in the extensive literature, GPs have evolved into not only a fruitful avenue for exploring fundamental aspects of quantum systems but also a topic of technological interest.For instance, due to its resilience to fluctuations in a coupled bath, the GP has been proposed as a significant resource for constructing phase gates [31][32][33] for quantum information processing.Superconducting circuits have been extensively investigated as the physical system allowing this aim [34][35][36][37][38].In a more fundamental approach, Berry phase has also been theoretically studied [39] and measured [40,41] in several circuit architectures.The vacuum GP accumulated by a superconducting artificial atom that was interacting with a single mode of a microwave cavity was measured [42] as well, while the corrections on the GP introduced by transitions to higher exited levels of the transmon were examined in [43].
When dealing with non-unitary dynamics leading to mixed states, the GP needs further generalization from its pure-state definition.A well defined proposal that applies under these conditions was presented in [28].Thereafter, this definition has been applied to measure the corrections induced on the GP in non-unitary evolutions [44] and to explain the noise effects observed in the GP of a superconducting qubit [41,45].Particularly, the GP of a two-level system under the influence of an external environment has been studied in a wide variety of scenarios [46,47].Even though the GP is not an immediate reflection of the dynamics and can therefore remain robust to the effect of the environment, it differs, in the general case, from that accumulated by the associated closed system, as the evolution is now affected by non-unitary effects such as decoherence and dissipation.Under suitable conditions, the non-unitary GP can be measured through interferometric (atomic interference) [48,49], spin echo [41], and NRM [44,50] experiments.
In a previous study [51], we thoroughly examined the GP accumulated by a two-level system (TLS) that was interacting with a single mode of the quantized electromagnetic field within a dissipative cavity, a physical setup known as a dissipative Jaynes-Cummings (JC) model.Addressing the scenario frequently encountered in semiconductor cavity quantum electrodynamics (QED) [52], the interaction between the atom-mode system and its environment was characterized by the flow of photons through the cavity mirrors and the continuous, incoherent pumping spontaneously exciting the TLS.
In the present paper, we extend the work on the JC model to encompass the scenario of a nonlinear transmon coupled to a transmission line or resonator.Additionally, both the transmon and resonator are coupled to two semi-infinite waveguides serving as the surrounding environments.We will investigate the dynamics of the composite system, both in the qubit sector and in the two-excitation sector.Restricting to the one-excitation sector will allow for direct comparison of the results previously obtained, which implies the comparison of two different architectures in which atom-cavity dynamics emerge.It is worth highlighting a major difference between both studies even at this early stage, which is the non-monotonic behavior encountered in the GP under certain environmental conditions.Thereafter, to further explore the richer nature of the transmon atom, we also examine the GP and its environmentally induced corrections when the two-excitation levels are involved in the dynamics.
In the next section, we will introduce the Hamiltonian describing the non-harmonic transmon-field system under investigation and the coupling of the composite system to the environment.In this section, we will also present some insights about the dynamic evolution of the system and the definition of the geometric phase.In Section III, we will describe the one-excitation subspace dynamics and the correspondence with previous results.The two-excitation space and the role of charging energy and non-harmonicity is discussed in Section IV.Section V summarizes our main conclusions.

II. TRANSMON WITH ATOMIC-KERR INTERACTION
Due to their substantial dimensions, stemming from the necessity of maintaining low charging energy (via large capacitance), transmon qubits inherently lend themselves to capacitive coupling with microwave resonators.This coupling is reflected in the transmon Hamiltonian Ĥ = 4E c (n−n g ) 2 −E J cos φ by the substitution of the classical voltage source V g with the resonator, a quantized gate voltage n g → −n r , representing the charge bias of the transmon due to the resonator.The Hamiltonian of the combined system is therefore [12] where n = Q/2e is the charge number operator, and φ = (2π/Φ 0 ) Φ (mod 2π) is the phase operator, defined by the charge and phase operators, respectively, of the quantum circuit.The charging energy is E c = e 2 /2C Σ , with C Σ = C J + C S being the sum of the junction's capacitance C J and the shunt capacitance C S .The operator nr can be written as nr = m nm , with nm = (C g /C m ) Qm /2e being the contribution to the charge bias to the mth resonator mode.In this expression, C g is the capacitance of the gate and C m is the corresponding associated resonator mode capacitance.It is usual to consider C g ≪ C Σ , C m .When assuming that the transmon frequency is much closer to one of the resonator modes than all the other modes, it is possible to truncate the sum over m to a single term.In this single-mode approximation for the resonator or cavity, the Hamiltonian reduces to a single oscillator of frequency denoted by ω r coupled to a transmon.
Expressed in terms of creation and annihilation operators, in the single-mode approximation, the Hamiltonian for the transmon-resonator cavity reduces to where the term n2 r has been absorbed in the charging energy term of the resonator mode and therefore leads to a renormalization of the resonator frequency, which we omit for simplicity.The frequency of the mode of interest is ω r , and the coupling constant between the artificial atom and the resonator is given by the relation g = ω r Cg /CΣ ( E J /2Ec) 1 /4 ( πZr /R K ) 1 /2 , where Z r is the characteristic impedance of the resonator mode and R K = h/2e 2 is the resistance quantum.The above Hamiltonian can be simplified further in the experimentally relevant situation where the coupling constant is much smaller than the system frequencies, |g| ≪ ω r , ω q .After rotatingwave approximation, the Hamiltonian reads This is the Hamiltonian we shall consider to study the interaction between the artificial atom (in the transmon regime E J /E c ≫ 1) and the electromagnetic field mode of the resonator.In this Hamiltonian, the term proportional to E c is the so-called Kerr-like interaction term, or non-harmonic term.
A. Coupling to the Environment Hitherto, our focus has been on quantum systems completely isolated from the surrounding environment.Nevertheless, a comprehensive portrayal of quantum circuits necessitates consideration of the manner in which these systems engage with their environment, encompassing both measurement apparatus and control circuitry.Indeed, the environment assumes a dual function in quantum technology; portraying quantum systems as entirely isolated is not only impossible due to inevitable coupling with undesirable environmental degrees of freedom but also renders a perfectly isolated system impractical for manipulation.Such a system would lack utility since we would be devoid of the means to control or observe it.Given these considerations, in this section, we investigate our quantum system coupled to external semi-infinite transmission lines that represent the measurement and control mechanisms while constituting the primary environment leading to photon loss and spontaneous decay of the transmon as their main environmental effect, respectively.We also consider the existence of a flux line for tunability of the transmon that leads to dephasing due to flux noise.
To study the open system comprising the transmonresonator and transmission lines (reservoir), we will employ the conventional formalism provided by Lindblad master equations.In this context, we can model the transmission lines as a set of harmonic oscillators, similar to the approach taken in quantum Brownian motion models, which are paradigmatic examples of open quantum systems.Assuming, therefore, that the transmon and the resonator are coupled to independent baths (this system is illustrated in Figure 1), the master equation for the composite system can be expressed as follows: where ρ is the density matrix of the composite system (transmon-resonator), and ĤK is the Hamiltonian of Equation (3).While this equation may suggest that dissipative processes independently impact the transmon and the resonator, the entanglement introduced by ĤK implies that events such as the loss of a resonator photon can result in qubit relaxation.In this master equation, the coefficient κ is the photon decay rate; γ represents the relaxation rate of the artificial atom, which is related to the qubit-environment coupling strength evaluated at the qubit frequency; and γ φ is the pure dephasing rate that superconducting quantum circuits can also suffer, caused, for example, by fluctuations of parameters controlling their transition frequency and by dispersive coupling to other degrees of freedom in their environment.The symbol D[ Ô]ρ in Equation ( 4) represents the dissipator where {. ; .}, is the anticommutator.We will work on the basis B = {|m n⟩}, where m refers to the mth energy eigenstate of the transmon and n is the nth Fock state of the mode field in the resonator.
Transmon capacitively coupled to the resonator and both of them capacitively coupled to transmission lines.The coefficient κ is associated with the coupling between the resonator and the readout transmission line and represents the photon decay rate.In addition, γ is related to the qubit-environment coupling strength introduced by the control mechanism and represents the relaxation rate.Finally, γφ is the pure dephasing rate emerging due to the flux noise taking place in the flux line that allows for tunability of the transmon frequency.

B. Dynamic Evolution of the Open System
To address the dynamics of the transmon-field composite system, we numerically solve the master Equation ( 4), constrained to the subspace with two or less excitations, so that the base B reduces to {|00⟩, |10⟩, |01⟩, |20⟩, |11⟩, |02⟩}.At a given instant t, the state of the system is described by a 6 × 6 density matrix which can be decomposed into blocks as in Equation (6).Elements belonging to different blocks satisfy differential equations that are decoupled so, by taking an initial condition with vanishing ρ ij elements for the offdiagonal blocks, the state of the system remains block diagonal along the whole evolution.
In this way, when the system is prepared in a oneexcitation state |Ψ(0)⟩ ∈ {|10⟩, |01⟩}, the evolution remains restricted to the subspace form by the 2 × 2 block and ρ 00 .In that case, the dynamics renders independent on the value of the capacitance energy E c generating the anharmonicity.On the other hand, by setting the initial state as a pure state within the 3 × 3 block, all the matrix elements in the block-diagonal subspace are involved in the evolution.

C. Geometric Phase in the Open System
As already noted in the introductory Section I, a generalized definition of a GP that is suitable to be computed for a mixed state under non-unitary evolution was proposed in [27].It reads where ψ k (t) are the instantaneous eigenvectors of the density matrix, and ϵ k (t) are the corresponding eigenvalues.This formula provides a well defined GP that, although defined for non-degenerate but otherwise general mixed states, when computed over pure states under unitary evolution, reduces to the expression defined over the most general unitary evolution of a pure state |Ψ(t)⟩ [23,53].It is also manifestly gauge invariant and therefore depends solely on the path traced by the state in the ray space.When dealing with pure initial states, ϵ k (0) = 1 for the specific k labeling the initial state and vanishes for all k ′ ̸ = k.Therefore, Equation ( 7) reduces to a simpler form where |ψ + (t)⟩ is the eigenvector of ρ(t) that coincides, at t = 0, with the initial state.This is the eigenvector such that ϵ + (0) = 1.Equation ( 9) has the exact same functional form of the GP defined over unitary evolutions, for which the only difference is that the pure state involved is the eigenstate |ψ + (t)⟩ of the density matrix and not the state of the system itself, which is now a mixed state.
We will restrict our analysis to pure initial states, so that the usual analysis applied to pure-state GPs can be immediately extrapolated by observing the behavior of the density matrix eigenstate |ψ + (t)⟩.

III. ONE-EXCITATION STATE: DISSIPATIVE JAYNES-CUMMINGS MODEL
It will be useful to start by exploring some aspects of the better-known subspace accessed when preparing the system in a state with only one excitation.In this case, as stated in Section II B, the state of the system remains constrained to the 2 × 2 block spanned by {|10⟩, |01⟩}, with the only exception of showing population exchange to ρ 00 .
On general grounds, the populations ρ 11 and ρ 22 of energy levels with one excitation oscillate while decaying, whereas the population of the vacuum state ρ 00 grows.Depending on the parameters, the asymptotic state can be a pure |00⟩ state or a mixed state with non-vanishing but suppressed ρ 11 and ρ 22 populations.The only nonzero coherences ρ 12/21 (t) increase in absolute value up to a maximum value and then vanish asymptotically.It is useful to notice already at this point that the relation between the parameters defining the unitary evolution and the parameters defining the environmental effects results in a dynamic that exhibits specific characteristics.This is noticeable in Figure 2, where the absolute value of the only non-vanishing coherences reaches a minimum |ρ 12 | ∼ 0 at t ∼ 7.5 τ , where it is smaller than the asymptotic value.This fact will be shown in different ways when observing the evolution of the eigenstate |ψ + (t)⟩ and the GP.In this archetypal example, we use ω q = 2π × 6 GHz, g = 2π × 166.85 MHz, and κ = 2π × 30 MHz, which constitute up-to-date typical values [54][55][56][57], while we consider a non-dispersive detuning ∆ = ω q − ω r originating from a resonator frequency ω r = 2π × 5.97 GHz.In what follows, we will keep the typical values for the artificial atom frequency ω q and the atom-mode coupling g while inspecting the dynamics arising in different conditions defined by different parameter values.In this sense, the atom-mode detuning ∆ will be modified from closer-to-resonance values ∆ ∼ O(10) MHz to values within the dispersive regime ∆ ∼ O(1) GHz.The artificial atom decay rate will also be increased to γ = 2π × 30 MHz.On the other hand, along the entire work, time is measured in units of τ = 2π/Ω, with the JC-Rabi frequency Ω = ∆ 2 + 4g 2 .
The eigenvector |ψ + (t)⟩ of the density matrix that is involved in the expression for the geometrical phase belongs, in this case, to the 2 × 2 block and can thus be observed on the Bloch sphere.In order to inspect the dependence of the dynamics with the detuning ∆, Figure 3 shows the evolution of |ψ + (t)⟩ on the Bloch sphere for a system that is prepared in a state |Ψ(0)⟩ = |01⟩ with one field excitation for three different relations of the detuning ∆ to the frequency ω q associated with the artificial atom.The remaining conditions are taken to be equal in all three cases.Panel (a) exhibits the case with ∆ = 0.0017ω q , panel (b) displays the case where ∆ = 0.017ω q , and panel (c) shows the case with ∆ = 0.17ω q .In all plots, the time is given by the color according to the color bar on the very right of the figure.(a-c) show the cases with ∆ = 0.0017 ωq, ∆ = 0.017 ωq, and ∆ = 0.17 ωq, respectively.The environment and remaining features of the system are the same in all three panels.The environment is characterized by a photon loss rate κ = 0.005 ωq and negligible atom decay and dephasing rates, and the atom-field coupling considered satisfies the relation g = 0.028 ωq.
Figure 3 shows that the path described by |ψ + (t)⟩ on the Bloch sphere is in all three cases a spiral that starts in the south pole of the sphere.The axis along which the spiral winds and moves differs in all three cases.Under the conditions in panel (a) of Figure 3, in which the system is closer to resonance, the spiral axis is almost the x-axis and the curve moves little along it.Thus, the state traces a path that slightly deviates from vertical rings.When increasing the detuning, the axis of the spiral tilts and the turns separate from each other, as visible in panel (b).By further increasing the detuning, the axis of the spiral gets closer and closer to the z-axis and the initial turns get again closer while spreading for longer times.
In some cases, as displayed in panels (a) and (b) of Figure 3, this behavior implies the exploration of different hemispheres of the sphere, in which the winding starts on one side of a certain (different in each case) great circle and crosses to the other side at some point.When this happens, the GP accumulated changes sign.In order to see this, Figure 4 shows both the GP accumulated as a function of time (a), and the corresponding path traced by the |ψ + (t)⟩ eigenstate on the Bloch sphere (b), for the characteristic example depicted in Figures 2 and 3b.
As expected, the accumulated GP has a sign while the state is winding the sphere on one hemisphere, and it changes sign as soon as the winding takes place on the other hemisphere.It can be seen from Figure 2 that this occurs at the specific time instant in which the coherences reach their minimum value.Therefore, we can classify the dynamics of the system initialized with one excitation as two kinds.One kind of evolution is that in which the coherences reach a minimum lower than the asymptotic value, the winding of the eigenstate |ψ + (t)⟩ changes hemisphere, and the GP changes sign, whereas the other is that in which none of these things happen.
To better explore under which physical circumstances the GP accumulation is non-monotone, we study the dependence in three characteristics of the system and environment.These are: the initial state of the system, which we take to be either |10⟩ or |01⟩; the main source of decoherence and dissipative effects, which we consider to be either the photon loss κ or the atom spontaneous decay γ; and the detuning ∆.The results of this examination are displayed in Figures 5 and 6, which show the GP accumulated in time for different combinations of the parameters characterizing the system.In both figures, solid lines represent the GP accumulated by the dissipative system, and the unitary results are introduced as dotted lines for reference.
Figure 5 corresponds to the previously explored case in which the photon loss process is the main source of dissi-pation.Three different ratios of the detuning to the frequency associated with the artificial atom ∆ = 0.0017 ω q , ∆ = 0.017 ω q , and ∆ = 0.034 ω q are displayed.Panel (a) shows the GP accumulated by a system prepared in the first excited level of the transmon atom |Ψ(0)⟩ = |10⟩, and, in panel (b), the initial state has a single field excitation |Ψ(0)⟩ = |01⟩.The two kinds of evolution, giving rise to monotonic or non-monotonic GPs, are clearly observed in Figure 5.In panel (a), the GP accumulated by an initial |10⟩ state is softer than the unitary result due to the environmental effects on the dynamics.When the state reaches the steady state and therefore stops moving on the ray space, the GP settles.On the other hand, in panel (b), the GP accumulated by an initial |01⟩ is non-monotonic, with the change of direction found sooner for smaller ∆/ω q ratios.Therefore, the dynamics leading to non-monotonic GPs are found when the main environmental effects are those affecting the initial excitation of the system.It is worth noticing that the results in Figure 5 are in full agreement with the results obtained in [51], in which the explored situation was that of a system prepared in an initial |10⟩ state and afterwards evolving in an imperfect cavity, corresponding to the case displayed in panel (a).
Likewise, Figure 6, where the GP accumulated in time is shown for the same three ratios of the detuning to the frequency of the transmon ∆ = 0.0017 ω q , ∆ = 0.017 ω q , and ∆ = 0.034 ω q , but the main environmental effect is the atom spontaneous decay.Consistent with the statement that non-monotonic GPs are found when the main environmental effects are those affecting the initial excitation of the system, in Figure 6 a change in the GP direction is only observed in panel (a), in which the initial state is the first excited level of the transmon atom and vacuum field |Ψ(0)⟩ = |10⟩.As was also found in Figure 5, the strongest the atom decay rate in relation to the frequencies associated with unitary evolution, the sooner the GP accumulation changes sign.
Comparing to previous results in [51], the absence of non-monotonic behavior found there can be explained as a combination of both the initial state in which the atomphoton system was prepared and the main environmental phenomena affecting it when the physical architecture is semiconductor cavities.In that case, the effect considered in Figure 6 was absent.
The evolution in Figure 6 can be re-observed, giving emphasis to the time instants by displaying it in the same manner as Figure 4.The analog plots compose Figure 7, which thus shows the GP accumulated as a function of time (a) and the corresponding path traced by the |ψ + (t)⟩ eigenstate on the Bloch sphere (b) for a system prepared in the first excited level of the transmon atom and vacuum field, with detuning-to-frequency rate ∆ = 0.017 ω q , and atom decay rate γ = 0.005 ω q .
Once again, the change in the direction of the GP coincides in time with the moment in which the path traced by |ψ + (t)⟩ crosses from one side to the opposite of a great circle.
In order to re-state the description in terms of the behavior of the coherence, we go back to the initially considered case in which the initial state of the system has a single field excitation |Ψ(t)⟩ = |01⟩.The analysis performed indicates there are two main situations in which the GP accumulated by this state will be monotonic: (a) if the relation between the environmental effects is such that the photon loss results are negligible in comparison with the atom decay rate, and (b) if, even though the main source of decoherence were the photon loss, the unitary parameters are strong enough to prevent the hemisphere crossing until the steady state is achieved.
Figure 8 shows the density matrix elements evolution for both these situations in panels (a) and (b), respectively.
It is immediately noticed that, in none of those cases, the absolute value of the coherences goes below its asymptotic value.In panel (a), where the photon loss is negligible in comparison with the atom decay rate, the excited populations ρ 11 and ρ 22 decay oscillating as the ground state gets populated.Different from what was observed in Figure 2, the absolute value of the non-zero coherences |ρ 12 | never reached a minimum below its asymptotic value.A similar behavior is observed in panel (b), where the increase in the detuning ratio ∆/ω q results in a relatively less strong environment that requires observation over longer timescales.As in panel (a), in this case the coherences absolute value is never below the asymptotic value.

IV. TWO-EXCITATION SPACE: ROLE OF CHARGING ENERGY AND NON-HARMONICITY
When turning to higher excited initial states, the evolution of the system takes place in the full six-dimensional Hilbert space described in Section II B. As stated there, if the system is prepared in a state with a defined number of excitations, the originally vanishing coherences remain zero at all following times.
On general grounds, the excited populations decay as lower-energy states get populated and, in the same way as described in Section III, the system evolves to an asymptotic state that can be either the pure ground state |00⟩ or a mixed state with non-vanishing but suppressed populations on exited states.The non-zero coherences increase in absolute value up to a maximum value, to vanish asymptotically afterwards.
In order to make the most simple generalization pos- sible of the one-excitation case, we consider in this section a system that is prepared in a two-excitation state |Ψ(0)⟩ = |11⟩, with excitations of a different nature.Figure 9 shows an specific example depicting the above described behavior.It can be observed in panel (a) of Figure 9 that the population ρ 44 (t) associated with the |11⟩ state decreases and shows non-harmonic oscillations.Along with this decrement, an immediate increment of the ρ 33 (t) and ρ 55 (t) populations, associated with the remaining twoexcitation states |20⟩ and |02⟩, takes place.These are suppressed within a few Rabi periods.Due to the effect of the environment, the ρ 11 (t) and ρ 22 (t) elements associated with one-excitation states also show an incre- ment and afterwards decrement while oscillate harmonically.The timescale of the increment-decrement of these elements is longer than the timescale shown by the twoexcitation populations.Finally, the ρ 00 (t) population associated with the ground state monotonically increases along the whole evolution up to a steady value ∼1.The non-zero coherences are now four (and their corresponding complex conjugates).Panel (b) of Figure 9 shows the evolutions of these elements.In all cases, the absolute value of these elements show an initial increment and asymptotic vanishing.Noticeably, the timescale as-sociated with the ρ 12 (t) coherence belonging to the oneexcitation subspace is larger than the timescale in which the coherences ρ 34 , ρ 35 , and ρ 45 , belonging to the twoexcitation block of ρ(t), evolve.
In order to reproduce the analysis performed for evolutions constrained to the one and zero-excitation subspace, we turn now to the observation of the GP accumulated by the state of the system in time.With this purpose, Figure 10 shows the GP as a function of time for systems characterized by different ratios of the detuning ∆ to the frequency ω q associated with the transmon first transition, but otherwise equal.On each panel, the main environmental effect is of a different kind: in panel (a) the main environmental effect is the spontaneous decay of the transmon, whereas in panel (b) the main environmental effect is the photon loss.
In panel (a) of Figure 10, where the initial |11⟩ state evolves subjected to the spontaneous decay of the artificial atom, the GP accumulated along the evolution differs from the unitary GP, in that the non-unitary increment is slower and softer up to a point in which it is completely stopped by the meeting of a steady state that does not move on the ray space any more.The bigger the atom-decay rate relative to the unitary parameters, the strongest the effect of the environment on the dynamics, which is reflected in the same way by the GP.
In addition, panel (b) shows the evolution of the same system for the case in which the main environmental effect is the photon loss.Under these circumstances, the non-monotonic behavior of the GP already observed in Section III is recovered.If the detuning is small enough, the environment has the effect of changing the sign in the GP accumulated, which afterwards tends to an asymptotic state as the system reaches a steady state that does not move in the ray space.However, the differences observed when comparing to Figure 5 are not only quantitative but also qualitative: the GP acquired for ∆ = 0.0017 ω q not only changes direction but also does so more than once.Instead of observing a single minimum step, there are two.
With the anharmonicity between levels being one of the main differences between these state subspaces, we further inspect the effect of the ratio E c /ω q with focus on this difference in the behaviors observed in Figures 5  and 10.

Effect of the Anharmonicity
In order to examine the effect of the anhamonicity E c on the GP accumulated by the non-unitary system, in panel (a) of Figure 11 we reproduce Figure 10b for a different value of the anharmonicity E c .In doing so, the behavior of both the unitary and non-unitary GP accumulated in time is qualitatively modified depending on the relation between parameters.
With regard to the unitary results, displayed in the figure in dotted lines, the GP accumulated when decreasing the anharmonicity remains qualitatively equal to the previous situation for the smallest ∆ = 0.0017ω q ratio, as depicted by the purple dotted line.However, when increasing the detuning, the unitary GP reaches a regime in which it changes sign periodically, leading to a step-like oscillation around a fixed value.
Regarding the GP accumulated in the presence of an environment, the non-monotonic behavior observed in Figure 10b completely disappears, and the GP performs steps that get softer and smaller as time goes by up to the stationary regime in which the state stops moving in ray space, so no further GP is accumulated.
To further inspect this effect, panel (b) of Figure 11 shows the GP accumulated in time for a system that is prepared in the same |Ψ(0)⟩ = |11⟩ state with fixed detuning ratio ∆ = 0.0017ω q and environment conditions described by κ = 0.005 ω q , but different values E c /ω q = 0.003, 0.035, and 0.067 of the anharmonicity E c .With increasing anharmonicity, the non-monotonic behavior is recovered, and more changes are observed in the GP sign for greater E c values.
As a whole, the GP accumulated by an initial state with one photon and the transmon in its first excited level shows the non-monotonic behavior only when the main effect of the environment is the photon loss, ob-served when both the anharmonicity splitting the transmon levels and the effect of the environment are big enough, leading to a system that is closer to a two-level system.Reducing the anharmonicity leads to a degeneracy in the artificial atom states that prevents the environmental non-monotonic behavior while introducing qualitative changes in the unitary GP.

V. DISCUSSION
In this paper we examined the dynamics beyond the two-level approximation of a transmon.Particularly, we studied the open dynamics of a nonlinear transmon coupled to a one-mode resonator and a transmission line.We have shown that the density matrix can be decomposed into blocks that satisfy differential equations that decoupled under particular initial conditions.Therefore, we were able to separately analyze one-excitation and two-excitation subspaces.We further explored the geometric phase accumulated by the state in order to have a better insight into the richer nature of the transmon artificial atom.
In the case of a one-excitation state, we retrieved results of the dissipative Jaynes-Cummings model.However, as we contemplated complex Linblad equations (with three different noise channels), we can complete already existing results.For example, we found the existence of a non-monotonic behavior in the accumulated geometric phase due to a combination of both the initial state and the main leading noise ruling the dynamics.
In the case of the two-excitation state, we studied the dissipative dynamics of the system in a regime where E c cannot be neglected.The open dynamics are more complex, but we were able to present studies on the numerically accumulated geometric phase on a 3 × 3 subspace and study its behavior as a function of the anharmonicity rates.Again, the GP accumulated by an initial state with one photon and the transmon in its first excited level shows the non-monotonic behavior only when the main effect of the environment is the photon loss.This can be better understood with the help of the complete analysis on the 2 × 2 subspace (where we can interpret results with the help of the Bloch sphere).
It is important to remark that in the near-resonance case, the accumulated geometric phase seems to remain robust, as previous results stated.This research was funded by Agencia Nacional de Promoción Científica y Tecnológica (ANPCyT), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), and Universidad de 728 Buenos Aires (UBA).
Figure 2 shows the explicit situation in which the state is prepared in a state |Ψ(0)⟩ = |01⟩ with one field excitation.

Figure 2 .
Figure 2. Dynamics of a system prepared in an initial |Ψ(0)⟩ = |01⟩ and characterized by a detuning ∆ = 0.017 ωq and an atom-field coupling g = 0.028 ωq.The environment is described by a photon loss rate κ = 0.005 ωq and negligible atom relaxation γ and dephasing γφ rates.Panel (a) displays the density matrix elements evolution, and panel (b) shows the Bloch sphere representation of the 2 × 2 block.In panel (a), the main plot displays the evolution of populations ρ00, ρ11, and ρ22 with dotted purple, dashed orange, and solid yellow lines, whereas the inset displays the absolute value |ρ12|.

Figure 3 .
Figure 3. Trajectory displayed on the Bloch sphere by the eigenstate |ψ+(t)⟩ of the density matrix for the non-unitary evolution of a system prepared in an initial state |Ψ(0)⟩ = |01⟩ for different ratios of the detuning ∆ to the artificial atom frequency ωq.Panels(a-c) show the cases with ∆ = 0.0017 ωq, ∆ = 0.017 ωq, and ∆ = 0.17 ωq, respectively.The environment and remaining features of the system are the same in all three panels.The environment is characterized by a photon loss rate κ = 0.005 ωq and negligible atom decay and dephasing rates, and the atom-field coupling considered satisfies the relation g = 0.028 ωq.

Figure 4 .
Figure 4. (a) Geometric phase accumulated over time by an initial state |Ψ(0)⟩ = |01⟩ and (b) Bloch sphere depicting the path traced by the density matrix eigenvalue |Ψ+(t)⟩ for a system with detuning ∆ = 0.017 ωq and atom-field coupling g = 0.028 ωq.The environment is characterized by a photon loss rate κ = 0.005 ωq for the artificial atom and negligible atom decay γ and dephasing γφ rates.The color depicts, in both panels, the time instant as indicated by the color bar on the right.

Figure 5 .
Figure 5. Geometric phase accumulated by systems with different values of detuning ∆, but otherwise equal.In panel (a) the system is prepared in an initial state |Ψ(0)⟩ = |10⟩, whereas panel (b) shows the case with |Ψ(0)⟩ = |01⟩.The detuning values ∆/ωq = 0.0017, 0.017, and 0.34 are depicted by the light blue, red, and purple solid lines, respectively.The atom-field coupling satisfies g = 0.028 ωq, and the environment is characterized by photon loss rate κ = 0.005 ωq and negligible atom-decay γ and dephasing γφ rates.The unitary results are included as dotted lines for reference, following the same ∆/ωq-to-color code.
Once again, panel (a) shows the case in which the system is prepared in the first excited level of the transmon atom |Ψ(0)⟩ = |10⟩, whereas in panel (b) the initial state has a single field excitation |Ψ(0)⟩ = |01⟩.

Figure 6 .
Figure 6.Geometric phase accumulated by systems with different values of detuning ∆, but otherwise equal.In panel (a) the system is prepared in an initial state |Ψ(0)⟩ = |10⟩, whereas panel (b) shows the case with |Ψ(0)⟩ = |01⟩.The detuning values ∆/ωq = 0.0017, 0.017, and 0.34 are depicted by the light blue, red, and purple solid lines, respectively.The atom-field coupling satisfies g = 0.028 ωq, and the environment is characterized by a decay rate γ = 0.005 ωq and negligible photon loss κ and dephasing γφ rates.The unitary results are included as dotted lines for reference, following the same ∆/ωq-to-color code.

Figure 7 .
Figure 7. (a) Geometric phase accumulated over time by an initial state |Ψ(0)⟩ = |10⟩ and (b) Bloch sphere depicting the path traced by the density matrix eigenvalue |ψ+(t)⟩ for a system with detuning ∆ = 0.017 ωq and spin-field coupling g = 0.028 ωq.The environment is characterized by a relaxation rate γ = 0.005 ωq for the artificial atom and negligible photon loss κ and dephasing γφ rates.The color depicts, in both panels, the time instant, as indicated by the color bar on the right.

Figure 8 .
Figure 8. Density matrix elements evolution in time for two different conditions under which the GP does not change monotony.In both cases, the main plot displays the evolution of populations ρ00, ρ11, and ρ22 with dotted purple, dashed orange, and solid yellow lines, and the insert displays the absolute value |ρ12| (solid light-blue line) and the vanishing of the remaining coherences |ρi0| (blue dotted line).Also in both panels, the system is prepared in an initial |Ψ(0)⟩ = |01⟩ with atom-field coupling g = 0.028 ωq.Panel (a) shows the case with ∆ = 0.0017 ωq, γ = 0.005 ωq, and negligible photon loss κ rate, whereas panel (b) addresses the case with ∆ = 0.017 ωq, κ = 0.005 ωq, and negligible atom decay γ rate.Dephasing is considered a subleading process on both plots.

Figure 9 .
Figure 9. Dynamics of the system, depicted by the time evolution of the density matrix elements.Panel (a) shows the evolution of the populations ρii(t), with line stiles associated with matrix elements as follows: the dotted purple line shows ρ00, the dashed red line shows ρ11, the solid orange line shows ρ22, the dot-dashed green line shows ρ33, and the double dotdashed light blue line shows ρ44.Panel (b) displays the evolution, in absolute value, of the non-zero coherences ρij , i ̸ = j.In this panel, the dotted purple line, dashed red line, solid orange line, dot-dashed green line, and double dot-dashed light blue line correspond to the matrix elements ρ12, ρ34, ρ35, and ρ45.The considered system is prepared in an initial state |Ψ(0)⟩ = |11⟩ in the first excited level of the artificial atom and one field excitation.It is further characterized by a detuning value ∆ = 0.0017 ωq, atom-field coupling rate of g = 0.028 ωq, and anharmonicity Ec = 0.035 ωq.The environment is characterized by photon loss κ = 0.005 ωq rate and negligible atom decay γ and dephasing γφ rates.The unitary results are included as dotted lines for reference.

Figure 10 .
Figure 10.Geometric phase accumulated by a system prepared in an initial state with one excitation of each kind |Ψ(0)⟩ = |11⟩ and characterized by different values of detuning ∆, but otherwise equal.The detuning values ∆/ωq = 0.0017, 0.013, 0.015, and 0.017 are depicted by the purple, red, light blue, and blue solid lines, respectively, the atomfield coupling satisfies g = 0.028 ωq, and the anharmonicity is Ec = 0.035 ωq.In panel (a), the environment is characterized by an atom decay rate γ = 0.005 ωq and negligible photon loss κ and dephasing γφ rates.On the other hand, in panel (b), the environment is characterized by photon loss rate κ = 0.005 ωq and negligible atom decay γ and dephasing γφ ratios.The unitary results are included as dotted lines for reference, following the same ∆/ωq-to-color code.

Figure 11 .
Figure 11.Geometric phase accumulated by a system prepared in an initial state with one excitation of each kind |Ψ(0)⟩ = |11⟩.Panel (a) reproduces (b), changing the values of the anhamonicity to Ec = 0.003 ωq.The considered detuning values are ∆/ωq = 0.0017, 0.013, 0.015, and 0.017, depicted by the purple, red, light blue, and blue solid lines, respectively, and the atom-field coupling satisfies g = 0.028 ωq.The environment is characterized by a photon loss rate κ = 0.005 ωq and negligible atom decay γ and dephasing γφ rates.On the other hand, panel (b) shows the GP accumulated for systems with different values of the anharmonicity Ec/ ωq = 0.003, 0.035, and 0.067, depicted by the purple, red, and light blue solid lines, respectively.The system and environment are further characterized by a detuning ∆ = 0.0017 ωq, an atom-field coupling g = 0.028 ωq, photon loss rate κ = 0.005 ωq, and negligible atom decay γ and dephasing γφ rates.On both panels, the unitary results are included as dotted lines for reference, following the same (a) ∆/ωq-to-color and (b) Ec/ωq-to-color codes.