Relating the curvature of De Sitter Universe to Open Quantum Lamb Shift Spectroscopy

In this paper, we explore the connection between the curvature of the background De Sitter space-time with the spectroscopic study of entanglement of two atoms. Our set up is in the context of an Open Quantum System (OQS), where the two atoms, each having two energy levels and represented by Pauli spin tensor operators projected along any arbitrary direction. The system mimics the role of a pair of freely falling Unruh De-Witt detectors, which are allowed to non-adiabatically interact with a conformally coupled massless probe scalar field which has the role of background thermal bath. The effective dynamics of this combined system takes into account of the non-adiabatic interaction, which is commonly known as the Resonant Casimir Polder Interaction (RCPI) with the thermal bath. Our analysis reveals that the RCPI of two stable entangled atoms in the quantum vacuum states in OQS depends on the de Sitter space-time curvature relevant to the temperature of the thermal bath felt by the static observer. We also find that, in OQS, RCPI produces a new significant contribution appearing in the effective Hamiltonian of the total system and thermal bath under consideration. We find that the Lamb Shift is characterized by a decreasing inverse square power-law behavior, $L^{-2}$, when inter atomic Euclidean distance, $L$, is much larger than a characteristic length scale, $k$, which is the inverse surface gravity of the background De Sitter space. If the background space-time would have been Minkowskian this shift decreases as, $L^{-1}$, and is independent of temperature. Thus, we establish a connection between the curvature of the De Sitter space-time with the Lamb Shift spectroscopy.


Contents
1 Introduction 1 2 Geometry of de Sitter space 5 2.1 Two atomic Wightman functions for probe external field in de Sitter space 6 3 Open quantum system (OQS) of two entangled atoms 10 3.1 Two atomic OQS model    55 1 Introduction Quantum theory of open system has been popular in recent times mainly because it provides a platform to study time evolution of system interacting weakly with its environment.
See for examples ref. [1] for the development of this subject. This development in Open Quantum SYstems (OQS) plays a crucial role for studying physical quantum phenomena in nature because they are mostly not isolated from its environment. It is important to develop a theoretical framework for treating these non-adiabatic interactions in order to obtain an accurate understanding of quantum mechanical systems. Theoretical tools and techniques developed in the context of OQS have been very successful in the context of quantum optics, quantum measurement theory, quantum statistical mechanics, quantum information theory, quantum thermodynamics, quantum cosmology, quantum biological systems. In a most generalised prescription, the time evolution of OQS is described by the non-adiabatic interactions between the physical system and its thermal environment. Consequently, the dynamical behaviour of the OQS cannot be accurately described using unitary time evolution operators alone after integrating out the bath degrees of freedom from the environment. The time evolution of OQS can be explicitly determined by solving the effective master equations of motion, also known as Gorini-Kossakowski-Sudarshan-Lindblad (GSKL) master equations, from which one can understand the non-unitary time evolution of the reduced density matrix of the system. In this context, the time evolution actually describes the time dependent behaviour of the system at different stages over time and also the dynamical behaviour of the observables that are associated with the OQS. The theory of OQS treats the system with dependent degrees of freedom as a subsystem in a much larger thermal bath [2]. Due to the complicated structure of the environmental degrees of freedom finding exact analytical solution of Gorini-Kossakowski-Sudarshan-Lindblad (GSKL) master equations are extremely difficult for practical purposes. Due to this difficulty, a variety of approaches have been developed in the present context. In this connection, a common objective is to derive the reduced time dependent description of the OQS wherein the dynamical behaviour of the quantum system are considered explicitly and the corresponding bath dynamics are described implicitly to know about the underlying physics of OQS under consideration. The seemingly obvious approach to deal with such systems, which is simulating evolution of both system and environment, would be naive as its complexity evolves exponentially. The OQS can no longer be defined by a pure state and to study the modified time evolutionary dynamics in presence of weak interaction between the system and thermal bath, density matrix formalism is required. We can study smaller subsystems to get across this problem by incorporating probabilistic description, where the quantum state of the subsystem is described by density operator. In an OQS, the effects of dissipation and decoherence is introduced by the environment degrees of freedom. The induced effects of decoherence and dissipation owing to the system is introduced by an operator more famously known as the Lindbladian [3]. While dealing with OQS, following crucial assumptions play significant role to describe the underlying physics:

Assumption I:
The combination of the system and the thermal environment is treated as a large closed system. Therefore we can assume that the time evolution is given by a unitary transformation generated by a global Hamiltonian.

Assumption II:
The interaction between system and environment is considered as Markovian, which describes the state of the quantum mechanical system in the next instant and dependent only on the current moment, not in the past. In short, the interaction between the system and thermal bath describes a phenomena without memory. This approximation is justified when the OQS under consideration has enough time for the system to achieve equilibrium state before being quantum mechanically perturbed again by nonadiabatic interactions with its thermal environment.

Assumption III:
The interaction between the system and thermal bath is considered to be weak in nature, which implies that the only change over time we see originates in the open system. This helps us to treat the time evolution of the system in a perturbative manner.

Assumption IV:
In the present context, additionally we have assumed that system and thermal bath are completely uncorrelated at the initial time.
In this paper, we have studied the connection of the curvature of de Sitter spacetime with the two atomic spectroscopy using the framework of OQS. To construct the theoretical set up we use a two body entangled system, which is described by two identical atoms. This two entangled OQS is considered as our theoretical probe which mimic the role of a pair of freely falling Unruh de Witt detectors. In this connection, it is important to note that, the concept of quantum entanglement describes a physical phenomenon which deals with strong correlations between the two atomic quantum states of the entangled particles [4,5] in OQS. Since this non-local property of quantum mechanics seemed very puzzling, Einstein, Podolsky and Rosen (EPR) argued for the existence of hidden variables in the context of quantum theory [6]. Later Bell proposed a set of inequalities to test their existence which when untrue, support the non-locality of quantum theory [7]. Bell's inequalities have to be violated in inflationary era to get persistent long range quantum correlation in the early universe cosmology and for this, the concept of quantum entanglement is commonly used. The violation of the Bell-inequality in de Sitter space has been addressed with axionic Bell pair to a great extent in several works [8,9]. Bell-inequality violation in various cosmological scenario with quenched time dependent mass profiles in de Sitter space has also been extensively studied in [10,11].
In the present context, we use mainly the concept of two body quantum entanglement in OQS which currently serves as one of the leading candidates to study the long range quantum correlation and many other unknown physics of the de Sitter space. The idea being that two or more entangled atoms which mimic the role of particle detectors in OQS can be used to measure the spacetime curvature of de Sitter spacetime and the consequent thermal effects, some of which we will discuss in detail in the following sections. In refs. [12][13][14][15], the authors have explored many more aspects of quantum entanglement in the background of a fluctuating scalar field in de Sitter space.In OQS the interaction of single particle with conformally coupled massless scalar field in de Sitter space is identical to the interaction of the same with thermal bath environment in Minkowski spacetime, making it difficult to distinguish between the two frameworks. This problem has led to further investigations with two body entangled OQS in order to shed light on the quantum structure and nature of the spacetime of our universe.
In a more general prescription, it is important to note that, when an accelerated particle detector moves through an external field in its vacuum state, spontaneous excitation may occur. The particle behaves as in the presence of a thermal bath, giving rise to an Unruh Temperature, which is proportional to its acceleration [16,17]. The excitation rate, of such a particle detector vanishes when it is at rest, but is found to be non-zero when uniformly accelerated. Same result is achieved if we analyse the vacuum state of the external conformally coupled scalar field in the frame of the accelerating particle detector. This is known as thermaliszation and shows that these accelerating detectors act as an OQS in the context of two body physics. Very recently, it has been explicitly shown that scalar fields coupled conformally to a de Sitter spacetime can be treated as an out-of-equilibrium system with a fluctuating background that introduced the thermal excitations in de Sitter space [18]. Such a detector-system combined set up sees the above thermalisation phenomena as a manifestation of the decoherence and dissipation due to the interactions with the surrounding. Furthermore, the approach followed in this work is also used to investigate the Gibbons-Hawking effect [19], which can be treated as a consequence of the thermalisation effect of the vacuum fluctuation of the external field in the de Sitter spacetime. These two phenomena encode the thermal behaviour of de Sitter spacetime.
Another very important outcome of vacuum fluctuations of quantised fields is Casimir effect which is heavily influenced by the curvature of de Sitter spacetime [20][21][22]. Casimir effect and the associated Casimir-Polder interaction (CPI) has been verified experimentally in multitudes of systems at microscopic and macroscopic levels. It has been used to study properties of quantum entanglement, long range effect of field correlations, Unruh effect etc. The spacetime curvature can alter the Casimir-Polder interaction [23][24][25] between the atoms [26] in two body OQS. The CPI in the de Sitter spacetime [27] has been investigated in detail in refs [15,[28][29][30]. In this paper, we aim to distinguish the curvature of de Sitter spacetime and Minkowski flat spacetime using Casimir-Polder interaction as a theoretical probe. In order to distinguish a conformally coupled massless scalar field with a de Sitter background to that of a Minkowski spacetime interacting with a thermal bath in which the excited Bell-states of the pair of atoms interacts with the background field placed at thermal bath via the exchange of real photons [31][32][33][34][35]. Additionally, in ref. [36] the authors have explicitly shown the use of Resonance Casimir-Polder interaction (RCPI) to detect curvature of de Sitter spacetime. In this paper, we use this methodology to know about the curvature of de Sitter by analysing the Lamb shift spectra computed from the two entangled atomic OQS. In [37], the authors have provided a methodology to distinguish between thermal Minkowski spacetime and de Sitter spacetime using the entanglement power. However, here we will use RCPI for the same purpose. Entanglement between quantum field modes, have been discussed in [38] (for scalar fields) and [39] (for fermionic fields). Also Unruh-DeWitt detectors and entanglement harvesting have been studied in case of different spacetimes in these beautiful papers [40][41][42][43]. Some other useful references regarding quantum entanglement are [44][45][46][47][48]. An excellent review of entanglment in expanding universe is given in [49], where they also briefly discuss the information encoded in such entanglement and other prospects. Dynamics of such a two qubit system in Unruh spacetime has been studied in [50]. In [51], the authors have calculated som information theoretic measures for such a system.
The various discussions in this paper can be summarised as follows: In section 2, we present a brief overview of the dS spacetime. The background metric has a hyperbolic geometry with a cosmological constant Λ as a solution of Einstein's field equations. First, we present the static dS metric. Then compute the two body two point Wightman functions for the external conformally coupled massless scalar field. Once we have these Wightman functions, we calculate the Fourier transform of these which we use later to calculate Lamb Shift Hamiltonian and Lindbladian.
In section 3, we describe our model two atomic 1 system conformally coupled to a scalar field in a thermal bath in de Sitter spacetime. Here, we allow our spins to have any orientation in the spacetime. Then we investigate the reduced dynamics of a pair of Unruh-De-Witt detectors or atoms to obtain th Master Equation. Then using the Fourier and Hilbert transforms of the Wightman functions (calculated in section 2) we write down the matrix elemnts of Lamb shift Hamiltonian and Lindbladian. The details of calculation is given in appendix A and B.
Further, in section 4, we have derived the expression for the energy shift to explicitly study the role of Lamb Shift spectroscopy [52,53] in the context of the two atomic entangled OQS. To serve this purpose, first of all we actually construct all possible two atomic entangled states out of the individual ground and excited states of the two atoms. By using the tensor product we have constructed four possible combination of the two body entangled states, which are ground, excited, symmetric and antisymmetric states. Further, using these two body entangled quantum states we compute the expectation values of the most relevant part of the effective Hamiltonian, known as the Lamb Shift Hamiltonian. Using these energy shifts we distinguish between the geometry of de Sitter spacetime from that of Minkowski spacetime. Here we have actually expressed the spectroscopic energy shifts in terms of the Euclidean distance between the two atoms separated in the space time and the inverse of the surface gravity of De Sitter space. Further, we have expressed the temperature of the thermal bath in terms of the surface gravity of de Sitter space. This allows us to connect macroscopic description of the environment in terms of the microscopic quantum mechanical observables of the two entangled OQS under consideration as well the geometry of our universe. Here particularly the microscopic observables can be characterised through few sets of spectroscopic integrals (see appendix B) which is in principle divergent in nature in QFT. To collect only the finite contributions out of these integrals we apply the Bethe regularisation technique. At the end one can express the surface gravity in terms of the curvature of de Sitter space, where curvature is related to the positive cosmological constant. This helps us to measure the curvature of de Sitter space from microscopic quantum spectroscopy.
FInally, we end the article with discussion given in conclusion and presenting all the necessary relevant calculations in the Apendices that follow.

Geometry of de Sitter space
The de Sitter space describes a universe with constant positive curvature and having the same degree of symmetry as a Minkowski spacetime, which also fits in with the picture of our current universe quite well. Since we want to work in this spaceyime, we will have a brief overview of the de Sitter spacetime in this section. Here our prime objective is to compute the two body two point (Wightman) correlation function between two atom conformally coupled with a massless probe external scalar field in de Sitter space. We will use the results obtained in this section in further sections to calculate our desired results. To serve this purpose we start with the background metric which is represented by the surface of the following hyperboloid: which describes a solution of the Einstein's field equations with the following radius of the hyperboloid: Here, Λ is the cosmological constant with positive signature in De Sitter space and the corresponding embedded metric in the five dimensional Minkowski space is given by [54,55]: By applying the following appropriate parametrisation one can express the five dimensional Minkowski flat metric in terms of a static four dimensional De Sitter metric, as given by: which are consistent with the equation of the surface of the hyperboloid in five dimension as expressed in Eq (2.1).
With the above parametrisation the static de Sitter metric in four dimensions can be expressed as: which is actually characterised by (t, r, θ, φ) in spherical polar coordinate in four dimensions. In the present context, to compute the explicit contributions of the two body Wight function of the probe scalar field present in the external thermal bath the geometry of the re-parametrised four dimensional De Sitter space play the most crucial role. In the next subsection, we actually compute these two body Wightman function in detail.

Two atomic Wightman functions for probe external field in de Sitter space
To compute the expression for the each of the entries of the two body Wight function of the probe scalar field present in the external thermal bath we use the four dimensional static de Sitter geometry of our space-time. In this set of coordinate system in four dimension, the Klein-Gordon field equation for the massless conformally coupled external probe scalar field for the non-adiabatic environment can be expressed as: where L 2 is the Laplacian differential operator in the three dimensions characterised by the coordinate (χ, θ, φ) , which is explicitly defined as: where we introduce a new coordinate χ which is related to the radial coordinate r as: r = sin χ. Now, in the present context of discussion we are not exactly interested to solve the wave function for the probe scalar field. Since our prime objective is to compute the two body two point correlation function for the probe external scalar field we construct the two body Green's function.
From the geometrical structure of the four dimensional static De Sitter metric it is obvious that the coordinate singularity, (1 − r 2 /α 2 ) −2 → ∞, which appears at the point r = α, and this is identified to be the cosmological horizon in the present context. Here it is important to note that, in flat space there is no problem to define the vacuum state (i.e. Minkowski vacuum) of the quantum field in open quantum mechanical system. But for curved space the definition of the vacuum state is more complicated in OQS than the result obtained in the context of flat space. For the computation of Wightman function from the present two entangled atomic open quantum set up, in the curved space we choose specifically the de Sitter invariant isometric SO(1, 4) group to connect the idea of spectroscopic energy shift with the geometry of De Sitter space-time. In this computation the isometric vacuum state (i.e. Bunch Davies and α-vacua) is actually identified with the quantum mechanical state of open system described by the conformally coupled massless probe scalar field. Now, the corresponding two point correlation function, often known as the Wightman function for massless probe scalar field takes the following form [54]: where, represents an infinitesimal constant, which is appearing in the representation of Wightman function in the i prescription. Also, we define the distance between two static atoms localized at the coordinates (r, θ, φ) and (r, θ , φ), appearing in this computation as: Here L represents the Euclidean distance between the coordinates (r, θ, φ) and (r, θ , φ), which is defined as: where, the angular difference ∆θ is defined as, ∆θ = θ − θ . Here both the signatures of the ∆θ is allowed in the present context depending on the relative ordering of the angular coordinates θ and θ . Additionally, the parameter in the present context represents an infinitesimal constant.The corresponding two body Wightman function between two spacetime points for massless probe scalar field can be expressed as 2 : where the components of the two body Wightman function can be expressed as: Two atomic Wightman function representing autocorrelation : From the two atomic two body system we get four Wightman functions in the present context due to the interaction between the two atoms with the external probe conformally coupled massless scalar field. The diagonal entries of the two body Wightman function represents the auto correlation function of the atom 1 and atom 2 with the external probe scalar field present in the thermal bath respectively. Also it is import to note that, these diagonal entries are same as we have considered two identical in our computation. On the other hand, the off diagonal entries of the Wightman function represent the cross correlation between the atom 1 and atom 2 with the external probe scalar field. Since in the present context the Wightman function is symmetric and constructed due to the interaction of two identical atoms with the external probe scalar field, the contributions appearing from the off-diagonal entries are exactly same.
Two atomic Wightman function representing crosscorrelation : where we use the following identity: where being an infinitesimal constant which we have already mentioned earlier. Furthermore, we have used the fact that since is an infinitesimal constant then one can approximate: sinh(i ) = i sin( )∼ i (2.13) cosh(i ) = cos( )∼ 1 (2.14) As a result, we get the following simplified result: Additionally, we have used the following definitions: with τ being the proper-time in the co-moving frame of the pair of atoms and k represents the surface gravity in the present context. Now that we have our Wightman functions, for later use, here we calculate the Fourier transform of the two point field correlation functions in frequency (ω) space for external massless probe scalar field as follows: where, we define the spectral function f (ω, L/2) as: In this expression the Euclidean distance L between the coordinates (r, θ, φ) and (r, θ , φ) is already defined earlier.

Two atomic OQS model
In this section, we investigate a model of OQS with entanglement. For simplicity, we consider a pair of identical static entangled atoms at an Euclidean distance L apart in the de Sitter spacetime. The two sets of internal energy levels of each atom are represented by: These are the ground-states and the excited-states for two atoms respectively. To avoid any further confusion here we note that, by the word "atom" we actually represented a very simplest spin bath open quantum model where two spin is immersed in the thermal bath. Here such spin pairs are characterised by the Pauli spin matrices. However, for a more general situation where we are interested in very complicated two atomic models of Hamiltonian the present analysis holds good, but solving that two body OQS problem itself extremely complicated. So for our better understanding we restrict ourself to a pair of spins which we are treating as a pair of atom, using which we will study various underlying physics of OQS. In this discussion, such pair of atoms are conformally coupled to a massless scalar field in the De Sitter background. The scalar field acts as a thermal bath for the pair of entangled atoms. In our discussion, the above two atomic entangled system are represented by a pair of Unruh-De-Witt detectors. These two identical atoms interact weakly with a quantised conformally coupled massless probe scalar field in its quantum mechanical vacuum state. Consequently, the corresponding two energy levels of the two atoms are identified as: Here, ω represents the renormalized energy level for two atoms, given by: Here K αα (±ω 0 ) for α ∈ {1, 2} are Hilbert transformations of two-point Wightmann functions discussed in previous section. Also, ω 0 represents the natural frequency of the two identical atoms. In this context, the atoms are characterised by the label α ∈ [1, 2] and σ α i ∀i ∈ [1, 2, 3] are the Pauli matrices.
The Hilbert space of such a system is bipartite in nature, i. e.
where H System and H Bath are the corresponding Hilbert spaces of the system and bath. Also, H Total represents the Hilbert space corresponding to the combined configuration of the system and the bath. This entangled two atomic OQS can be represented by the following total Hamiltonian: where I System and I Bath are the identity operators defined for the system and bath. When we are accessing the system we don't see anything from the bath and converse is also true.
Here, the identity operators as appearing in the system and bath corresponds to the no access. More precisely , if one observer is sitting on the system, made by atom 1 and atom 2, then that observer will not feel any further effect from the thermal bath environment. On the other hand, once the observer is sitting at the reference frame of thermal bath, that observe will not see any further effect from the two entangled OQS. However, in the interaction term both the system and thermal bath explicitly contribute and due to the entanglement one cannot separate their contribution from this term in the Hamiltonian.
Here τ represents the proper time in the co-moving frame of two atoms. To make things clear, we would like to highlight the fact that instead of proper time one could use the static time cooardinate(t). In that case, one might have to deal with the inherent curvature of de-Sitter spacetime explicitly in every calculation. However, the use of proper time incorporates this curvature effect directly through it's relation with static time, as will be discussed later.
More details of each of the terms of the total Hamiltonian is as follows: 1. System: The system Hamiltonian of the two entangled atoms are described by the linear combination of the individual contributions coming for each atom 3 : where, the sum is taken over two individual contribution appearing from the Hamiltonians of the two atoms. In the above expression the normal vectors for two atoms are represented by n α i ∀α = 1, 2, & i = 1, 2, 3 and the corresponding projection of Pauli-matrices are characterised by the direction cosine of Euler angles α, β and γ respectively. 3 Here we define the following direction cosines for two entangled atoms:

Thermal bath:
The thermal bath Hamiltonian in the present context is characterised by a free rescaled field Φ(x) = a(τ )χ(x), where the original massless scalar field χ(x) is conformally coupled with De Sitter background with scale factor a(τ ). The Hamiltonian in the de Sitter background can be written as: where Π Φ (τ, r, θ, φ) is the canonically conjugate momentum of the scalar field Φ(x).

System-thermal bath interaction:
The interaction Hamiltonian is characterised by the following expression: where, µ represents the coupling between the pair of atoms and the external massless scalar field placed at the thermal bath. While deriving the reduced dynamics of the two atomic detectors, we consider the weak coupling limiting approximation between the pair of atoms and the external massless probe scalar field. We assume that the coupling parameter µ is sufficiently small so that perturbation theory is applicable in the open quantum mechanical system under consideration.

Time dynamics of two atomic OQS
To construct the effective Hamiltonian from the time dynamics of the present system, we first consider that there is no correlation between the pair of atoms with the the external probe free scalar field. Hence, the thermal density-matrix of the combined system and thermal environment can be expressed in the following form: where ρ System (0) and ρ Bath = |0 0| represent the initial density matrix for the pair of atoms (system) and environment (bath) degrees of freedom. Here, |0 characterise the vacuum state of the external free scalar field. Now it is important to note that, in the interaction picture of OQS, the time evolution of the total density matrix can be written in the following form: One can further write the most general structure of the time evolved version of the total density matrix of the combined system and thermal bath, as given by: (3.12) Now, we have already started with the assumption that the non-adiabatic interaction between the system and thermal bath is switched on at t = 0 (from initial structure of the total density matrix) and prior to that the interaction between the system and the thermal bath is extremely weak or absent for which there is no correlation exists in the OQS under consideration i.e. ρ correlation (0) = 0 (3.13) The idea of having no correlation at the initial stage is not very restrictive in nature, since for the computational purpose we can always find a time prior to which the system and environment did not interact. In the weak coupling limiting situation, the time evolved density matrix can be further simplified as: which is valid through a time scale on which the perturbation theory in OQS is valid. Furthermore, we also assume that the time scale corresponding to the correlation, τ Bath , which is sometimes identified to be the relaxation time scale of the thermal bath, assumed to be extremely weak. Here one can write: Time evolution of the bath for time scale τ >> τ Bath : ρ Bath (τ ) ≈ ρ Bath (0) (3.15) In this context, the time evolved total density matrix can be represented as: where U Total represents the time evolution operator of OQS under consideration. Since, we are only interested in the reduced dynamics of the density matrix of the entangled atoms only, therefore, we trace out the external field (thermal bath) degrees of freedom from the total combined open quantum two atomic system. Consequently, the reduced density matrix of such a two atomic system can be expressed as: Further, doing simplifications one can write the reduced density matrix of the system after performing the partial trace operation over the bath degrees of freedom in the following simplified form: The subscript T in the above equations refers to the term Total. where |k represents all possible orthonormal basis states of the thermal bath defined in the Hilbert space H Bath .
Here, additionally we define the operator M k as: which is defined on the Hilbert space corresponding to the system under consideration i.e. H System . Since in the present context, the time evolution operator U T is unitary i.e. U † T U T = I T , which further implies that: Now, we explicitly mention about the properties of the system density matrix: 1. Hermiticity: The system density matrix is hermitian in nature which can be tested as: The last step can be justified using the fact that the density matrix at time t = 0 is hermitian.

Positivity:
The system density matrix also satisfies positivity property. Now, the non-unitary time evolution of the reduced density matrix in the weak coupling limiting situation can expressed in terms of the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation in OQS, as given by: where H eff is the effective Hamiltonian of the two atomic system under consideration, which incorporates the effect of inter atomic interaction aka Resonant Casimir Polder Interaction (RCPI). Also, the last term in the above mentioned evolution is known as the Lindbladian 4 , which describes the dissipative contribution due to the influence of the thermal bath on the two entangled atomic system. In the following subsections, we discuss about the effective Hamiltonian and the Lindbadian with greater detail.

Effective Hamiltonian construction
As we are primarily interested in entanglement properties of atoms we trace out the bath degrees of freedom. When one does that in the path integral formalism, one obtains an effective Hamiltonian. For our setup of a two qubit system interacting with the scalar field this is given by: Lamb shift=Heisenberg spin chain (3.23) where the first term in the effective Hamiltonian represents the system Hamiltonian of the two atomic system which are interacting with each other and one can treat this term to be a self interaction at the level of Hamiltonian. The specific details of this issue we have already mentioned earlier. The second part necessarily captures the interaction between the two atoms as it has contribution from both atoms. In condensed matter theory (e.g. in Ising model), one usually considers nearest neighbour interactions to study the quantum mechanical behaviour and phase transition. Our second part of the above Hamiltonian resembles with that well known Ising model Hamiltonian and mimics the role of an spin-spin interaction Hamiltonian in Heisenberg spin chain model which arises from the interaction between the two atoms and the external field. The terminology Lamb shift was first used in the context of a hydrogen-like atom where in the lowest order approximation in the fine structure constant energy level shift was determined from an effective Hamiltonian. Following the same, we and also in other refs. [27,36,50] people have used the terminology Lamb shift in the present context to describe the energy level shift or spectroscopic shift of the effective Hamiltonian describing spin-spin self interaction. Hence for usual practice, we will name this specific part of Hamiltonian, in our case also, as Lamb Shift Hamiltonian. Here,n is the the normal unit vector for individual atom under consideration. Also the angles between normal vector and the Pauli matrices are characterised by the three Euler angles, α, β and γ. Therefore, the Lamb Shift Hamiltonian [27,36,50] can be simplified in terms of the Eulerian angles as: Now, to know the explicit contribution in the effective Hamiltonian in the present two entangled atomic set up we define the following set of Pauli operators: which is actually expressed in terms of the tensor product of 2 × 2 identity matrix σ 0 and the three 2 × 2 Pauli matrices σ i , ∀i = 1, 2, 3, which satisfy the following conditions: Pauli matrix algebra : Here also it is obvious from the mathematical structures of these Pauli operators obey the following algebra: where (i, j) = 1, 2, 3 and (α, β) = 1, 2. With these set of definition of the Pauli matrices the Pauli operators can be expressed in terms of the tensor product as: In this context, one can obtain the matrix elements of H αβ appearing in the Lamb shift Hamiltonian from the two point correlator which is given as: The Fourier-transform of the above two point propagator in frequency (ω 0 ) space can be written as: where ω 0 represents the difference between the energy levels of the ground and the excited states of the atoms respectively. In turn, the elements of the effective Hamiltonian matrix H αβ ij are given by the following Hilbert transform of the Wightman function (two point correlator) as given by: where P is the Principal value of the integral. Now, the elements of co-efficient matrix H (αβ) ij of the effective Hamiltonian can be explicitly represented by the following expression: where, the quantities A (αβ) and B (αβ) for the two atomic system are defined as: Here, µ is the coupling parameter which represents the interaction strength between the system and the external thermal bath (conformally coupled scalar field) degrees of freedom.
Determining the structure of the elements of the co-efficient matrix H (αβ) ij in terms of the two atomic two point correlation function (Wightman function) of external free conformally coupled massless scalar field in de Sitter background actually fixes the structure of the effective Hamiltonian in the present context.
We define the following uantities (equality sign holds because of symmetry of Hilbert transformations as seen in Appendix A.
From equation 3.34, the non-vanishing Hamiltonian matrix elemnts are as follows : and these will explicitly contribute to the final expression for the expectation values of the Lamb shift Hamiltonian.

Lindbladian construction
The fluctuation-dissipation into the system is introduced by the additional contribution in the time-evolution equation of the reduce density-matrix often known as the Lindbladian [3] or the Lindblad operator in OQS. The second significant term in the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation is actually characterised as the Lindbladian or Quantum Dissipator of the OQS, which is for the present two atomic model can be written as: where, ρ System (τ ) is the reduced density matrix of the two entangled atomic system where we have trace over the external bath scalar field degrees of freedom. The co-efficient matrix C αβ ij is known as the Gorini-Kossakowski-Sudarshan-Lindblad matrix, which is constructed under the weak coupling limiting approximation on the coupling parameter µ, as appearing in the interaction Hamiltonian 5 . In the context of any OQS, Lindbladian captures the effect of dissipation implicitly describing by the system operator (n α i .σ α i ), which one can treat as an influence of the thermal bath on the two entangled atomic system under consideration.
Following the same procedure performed in the previous section to compute the coefficient matrix element , H , as appearing in the expression for the Linbladian can be expressed as: In the case of closed quantum mechanical system all the entries of the Gorini-Kossakowski-Sudarshan-Lindblad matrix is zero. Consequently, the time evolution of the reduced density matrix is described by the quantum Liouville equation, which is the quantum mechanical analog of the classical Liouville equation.
where, the quantitiesÃ (αβ) andB (αβ) for the two atomic system are defined as: In the later half of the paper we will explicitly compute the entries of the Gorini-Kossakowski Sudarshan-Lindblad matrix, C In this paper, one of our prime objectives are to find out the expectation value of the of the Lamb Shift Hamiltonian of the OQS described by two atoms, which are entangled with each other. Fixing the co-efficient matrix H (αβ) ij in terms of the two atomic two point correlation function (Wightman function) of external free conformally coupled massless scalar field in De Sitter background helps us to compute the analytical expression for the energy shift explicitly. To compute this expression the main ingredient is the all possible quantum mechanical states which we have to construct in the present context from the ground and excited states of the two atoms. Using these atomic states we construct four possible quantum mechanical two atomic entangled states for the combined system (two atomic system+thermal bath), which are ground state, excited state, symmetric state and antisymmetric state respectively. We explicitly do this computation in the later part of the paper.

Lamb Shift spectroscopy from OQS of two entangled atoms
Lamb shifts are obtained by taking the expectation value of Lamb shift Hamiltonian in all the states of the two atomic system. Here, we first present all the states for the atomic system formed from the individual atomic states.

Bethe regularised Lamb Shift Spectra
Using these spectroscopic integrals we obtain the simplified result for Lamb shifts

Bethe regularised energy shift from Ground & Excited states:
Here, we evaluate the contributions coming from the ground state and excited state expectation values, as given by: Under a certain case for which cos 2 α 1 cos 2 α 2 + cos 2 β 1 cos 2 β 2 = 0 we have, This directly implies that in this particular orientation of the qubits guided by the values of the α's and β's satisfying the above equation the ground and excited states does not contribute in the Lamb shift spectroscopy of two entangled qubits in OQS.

Bethe regularised energy shift from Symmetric & Antisymmetric states:
Now, from the present analysis we observe that the RCPI between the two entangled atoms in the De Sitter background is being contributed only by the symmetric and antisymmetric part of the Hamiltonian H LS as it consists of the term f (ω, L/2), which contains a measure of the Euclidean distance L between the two entangled atoms. This actually contributes in the shift in the energy levels or more precisely the inter atomic interaction energy computed from the symmetric (|S ) and antisymmetric (|A ) quantum state constructed out of two entangled atoms in De Sitter space. On the other hand, we have already seen that there is no such term present in the shift in the energy levels between the ground state and excited state constructed solely from the quantum states |G and |E for two atoms. This is appearing due to the non inter atomic interaction appearing between the uncorrelated two atomic quantum states in the second order of perturbation theory of OQS under consideration in this work. The presence of the Euclidean distance dependent term in the expression of the symmetric and the antisymmetric part of the shift in the energy level defines the gradient of a potential between the two atoms in the curved de Sitter space and is the manifestation of the RCPI between them. Therefore, only the terms which contributes towards the RCPI between the two entangled atoms are given by the following expression: δE S (L) = − µ 2 P 8π 2 [ΩΓ 1 (α 1 , α 2 , β 1 , β 2 )∆ 3 + cos 2 (α 1 ) + cos 2 (α 2 ) + cos 2 (β 1 ) + cos 2 (β 2 ) ∆ 2 ] δE A (L) = − µ 2 P 8π 2 [ΩΓ 1 (α 1 , α 2 , β 1 , β 2 )∆ 3 + cos 2 (α 1 ) + cos 2 (α 2 ) + cos 2 (β 1 ) + cos 2 (β 2 ) ∆ 2 ] (4.11) where we have introduced two functions Γ 1 (α 1 , α 2 , β 1 , β 2 ) and Γ 2 (α 1 , α 2 , β 1 , β 2 ), which are defined as: Using the result obtained in Appendix B, we get the following simplified result for the inter atomic interaction energy shift: We get the following characteristic features: 1. If we increase the value of the Euclidean distance (L), then the oscillation in the Lamb Shift increase with respect to the surface gravity in the region L ≤ k.
2. We also found that the magnitude of the saturation value of the Lamb Shift decrease with increasing value of the Euclidean distance (L) in the region L > k.
Further, it is important to note that, the final results of the inter atomic energy shifts depend on the background De Sitter metric through the following relation: This directly implies that the parameter k is directly related to the positive cosmological constant of De Sitter space. Consequently, one can theoretically probe De Sitter space using a pair of entangled atoms in OQS in presence of RCPI. Here it is important to note that the result obtained for inter atomic energy level shift for two detectors (two entangled atoms) can be interpreted as the energy level shift obtained for a single detector immersed in a thermal bath with temperature: which is interpreted as the Unruh Temperature. A freely falling observer under a steady acceleration observes this temperature in de-Sitter space. In this case, the inter atomic interaction exhibits non thermal behaviour and carrying non thermal fluctuation. Now, to understand the detailed physical features of the obtained result for energy level shift from RCPI in de Sitter space we consider two limiting situations, as given by:

Case I:
When the inter atomic distance is much larger than characteristic length scale k i.e. L >> k. In this case the two entangled atomic system is placed near to the cosmological horizon. In this limit, the energy-level shift from RCPI can be simplified as: This result shows that the energy level shift is non trivially dependent on the parameter k because of that fact that in this case curvature of the De Sitter space is significant.
Here we observe that the RCPI in the limit L >> k falls as 1/L 2 , which shows that the first order correction to the energy maintains the inverse square law in De Sitter space.

Case II:
When the inter atomic distance is of the order of the characteristic length scale k i.e. L ∼ k. In this case the two entangled atomic system is placed exactly at the cosmological horizon. In this limit, the energy-level shift from RCPI can be simplified as: This result shows that the energy level shift is dependent on the parameter k because of that fact that in this case curvature of the de Sitter space is comparable to the atomic distance, which will give rise to the following constraint condition: Here we observe that the RCPI in the limit L ∼ k falls as 1/L, which shows that the first order correction to the energy maintains the inverse law in Minkowski space.

Case III:
When the inter atomic distance is much smaller than characteristic length scale k i.e. L << k. In this case the two entangled atomic system is placed far from the cosmological horizon. In this case it is possible to find a local inertial frame of reference where all physical principles coincides with that in Minkowski space. In this limit, the energy-level shift from RCPI can be simplified as: This result shows that the energy level shift is independent on the parameter k because of that fact that in this case curvature of the De Sitter space is negligibly small. This result exactly matches with the result obtained for Minkowski space. Here we observe that the RCPI in the limit L ∼ k falls as 1/L, which shows that the first order correction to the energy maintains the inverse law in De Sitter space and exact Minkowski space.
Additionally, in all of these physical limits we found the following features: • Feature I: We also observe that the RCPI contains the Eulerian angles within it in both the limiting results. This suggests that the shift in energy is described by the Eulerian angles which means that the RCPI is a function of the direction along which the spin of the atoms is directed along. In de Sitter space the RCPI is dependent on how the spin of the two atoms are oriented. The orientation of the spins of the two atoms which is determined by the Euler angles of rotation α i , β i and γ i ∀i = 1, 2, 3, quantifies the RCPI along arbitrary direction of spin projection.
• Feature II: The Euler angles of rotation determines the manifestation of the thermal environment that the atoms see in their comoving frame. Furthermore, we observe that the pre-factors in the energy shifts determines the parameter k associated with the temperature [13,14,19,[57][58][59] of the thermal bath: where, the Gibbons-Hawking temperature and Unruh temperature are defined through the following expressions: with the proper acceleration given by: which is defined in the co-moving frame of the two entangled atoms for the given OQS under consideration.
• Feature III: Now, it is important to note that in De Sitter space the curvature can be quantified though the Ricci scalar, which can be further expressed in terms of the cosmological constant as: As a result, the Gibbons-Hawking temperature and Unruh temperature can be expressed in terms of the Curvature of De Sitter space as: Consequently, the temperature of thermal bath can be expressed in terms of the curvature of de Sitter space as: In this case, the RCPI can be expressed in terms of the curvature of the de Sitter space as: In the limiting situation, L >> k = 12 1 − R DS r 12 2 /R DS , the RCPI can be expressed in terms of the curvature of the de Sitter space as: Here one can consider another limiting situation, where L ∼ k = 12 1 − R DS r 12 2 /R DS .
For this case, RCPI can be expressed in terms of the curvature of the de Sitter space as: • Feature IV: Now, if we take the limit R DS → 0 ( i.e. Λ → 0 or α → ∞), then we get: In this case te RCPI will be reduced to the result obtained in the limiting situation L << k, which is exactly same result as obtained for the Minkowski space.
• Feature V: Here additionally it is important to note that, only if the Unruh temperature vanishes, the corresponding proper acceleration of static atom vanishes i.e. a = 0 and this can be obtained when the atoms are localised at r = 0. As a result, kinematic contribution will not appear in the expression for RCPI. However, in this case, the RCPI is still can be expressed in terms of the curvature of the De Sitter space, as in this case temperature of the thermal bath is quantified by the non vanishing Gibbons-Hawking temperature i.e.
In this case, the RCPI can be expressed in terms of the curvature of the de Sitter space as: In the limiting situation, L >> k = 12/R DS , the RCPI can be expressed in terms of the curvature of the de Sitter space as: Here one can consider another limiting situation, where L ∼ k = 12/R DS . For this case, RCPI can be expressed in terms of the curvature of the De Sitter space as: • Feature VI: Now, we compare the obtained results for Lamb Shift in De Sitter space with the result corresponding to the Minkowski space. For this purpose we consider a specific situation where two static atoms are interacting with the environment, where it is represented by the massless scalar field in OQS. In this system, the two point field correlation can be expressed in terms of the Wightman function given by: where L is the Euclidean distance between two atoms in OQS. Using this Wightman function we can carry forward the similar calculation for Lamb Shift from the RCPI in Minkowsi space, which will finally give rise to the following expression: From the above mentioned result it is clearly observed that the Lamb shift obtained from RCPI in Minkowski space not containing any contribution from the temperature of the thermal bath, T = 1/2πk and only depends on Euler angles and the Euclidean distance L. Also we found that this result exactly matches with the result obtained for the case L << k with two inertial atoms in the static patch of de Sitter space. This result additionally implies that the inter-atomic interaction between two atoms behave differently in Minkowski space and in the static patch of de Sitter space. For two or more atoms one can construct the ground, excited, symmetric and antisymmetric entangled states, out of which symmetric and antisymmetric entangled quantum states will give rise to non-zero spectroscopic shifts in the two limiting situation, L k and L k. Using these limiting results one can able to probe RCPI and distinguish between the static de Sitter and Minkowski flat space-time in a very precise manner. But instead of using two atoms if we use only one single atom within the framework of OQS then one cannot implement the methodology quantum mechanical entanglement in the present scenario. For the single atomic non-entangled case one can only construct the ground state and excited state, but there will be no symmetric and antisymmetric states exist, which are .the necessary ingredient to probe RCPI using which one can able to distinguish the effects of Minkowski flat space-time and static patch of de Sitter space-time. For two or more atomic scenario the atoms are entangled with each other and interacting with a thermal bath which is modelled with scalar field embedded in the static patch of de Sitter space. In this scenario, RCPI can be used to study the spectroscopic shift from the spin-spin self interaction from the effective part of the Hamiltonian. On the other hand, when we describe the present scenario with a single atom which is interacting with a thermal bath modelled by a scalar field in static de Sitter background due to the absence of having any quantum entanglement RCPI will be replaced by a very simple self interacting quadratic term. Consequently, using the spectroscopic shift for single atom-bath system one cannot distinguish between the Minkowski space-time and static patch of the de Sitter space. Here two or more atomic OQS is the only feasible option using which one can precisely distinguish the behaviour of RCPI in Minkowski space and in the static patch of de Sitter space due to having quantum entanglement between the atoms which are interacting with the thermal bath. In short, the phenomena of quantum entanglement play very significant role in the present context of discussion as this can be only explained through two or more atomic systems to probe the curvature of the space-time through RCPI. in the spectroscopic shift computation. But as the quantum entanglement is not there in single atomic system, one cannot probe the curvature of the space-time in the present context.
• To understand the obtained results in a more better way and to compare with the previously obtained results in the same literature we discuss now two limiting situations, which are appended below: 1. Limit I: In this case we fix the position of the two detectors at the angular positions (α 1 = π/2, β 1 = π/2, γ 1 = 0) and (α 2 = π/2, β 2 = π/2, γ 1 = 0) respectively. In this limit, H A matches with the model studied in previous works, but H LS does not.
Our work Limit I Limit II In table 4.1, we have explicitly mentioned the comparison between our work and previous works.
In fig: 1, the behaviour of the spectroscopic Lamb shifts have been plotted with respect to the Euclidean distance between the two qubits. It is clearly observed that for smaller values of L, the shifts fluctuate with a very large amplitude. The rate of fluctuation increases with increasing L with a decreasing amplitude and for asymptotically large value of L the shifts diminish to zero. A naive interpretation from this plot can be that to measure a finite shift, smaller values of L is preferred to larger values.
In fig: 2, the spectroscopic Lamb shifts are plotted against the temperature of the thermal bath. It can be seen that for lower temperature of the bath a finite value of the spectral shift is obtained. However, after a certain characteristic temperature the spectral shifts start decaying and finally reach zero. This finite value of the spectral shifts for smaller values of the temperature of the thermal bath suggests that by this indirect mechanism one can expect a large value of the inverse curvature parameter k for the static patch of de Sitter space, which is related to the temperature of the thermal bath via the relation T = 1/2πk.
In fig: 3, we have plotted the characteristics of the spectroscopic Lamb shifts with respect to the curvature of the de Sitter background. It can be seen that for lower values of the curvature a finite value of the spectral shift is obtained. However, after a certain characteristic value of R DS the spectral shifts start diminishing in an oscillatory fashion and gradually go to zero. For larger values of R DS , the spectral shifts oscillates rapidly to give negligible value. Although from this indirect detection mechanism, it is not possible to determine the exact value of the de Sitter curvature but one can predict that its value is smaller than a certain characteristic value and can never be larger than that.

Conclusion
In this paper, we have established a direct connection between the curvature of De Sitter space and atomic spectroscopy for an OQS described by two entangled atoms. For this purpose, without doing any terrestrial and space observations designing a laboratory atomic experiment is very useful to establish the connection between space time geometry and quantum mechanics. To summarise, in this work, we have addressed the following issues to implement the above mentioned idea: • To begin with, we have started our discussion with an OQS characterised by two entangled atoms. We have considered the two body quantum entanglement as in this situation it is allowed to exchange energy through Lamb Shift in terms of the geometry of De Sitter space-time. In this theoretical construction these two atomic pair represents Unruh-De-Witt detectors, which are considered to be conformally coupled to a background scalar field in thermal bath.
• The non-adiabatic interaction between the detectors and the thermal bath in OQS is characterised by RCPI, which plays the key ingredient to determine the curvature of De Sitter space from the Lamb Shift spectroscopy.
• In order to study the full dynamics of the two entangled atoms in OQS for any arbitrary position of two detectors, we have used a generalised Hamiltonian described by Pauli operators which include contributions from the Euler rotation angles due to projection in any arbitrary direction. This direction of projection actually playing the role of direction of observation of the atoms in the atomic detectors. In this discussion, the Lamb Shift Hamiltonian includes a term that arises from the interaction between the atomic detectors with the background test scalar field. This is the most significant term which can be experimentally probed using atomic spectroscopy to detect to geometry of De Sitter space-time. In the time dynamics of the reduced density matrix • To compare between the geometrical features of De Sitter and Minkowski flat space we have also computed the Lamb Shift from the two body entangled OQS set up in the thermal state for the flat case. In the case of De Sitter space the Lamb Shift is described by inverse square power law dependence on the Euclidean distance (L), which is characterised by the length scale associated with the breaking of local inertial description of the two entangled atomic OQS. On the other hand, in the Minkowski flat case we do not get any temperature dependence and in this case the the spectroscopic Lamb Shift is described by L −1 behaviour. From this discussion it is evident that, even both thermal Minkowski and De Sitter space satisfy similar kind of properties and cannot explicitly discriminated by a single external probe field, but using non adiabatic RCPI it is possible to differentiate between these two geometrical spacetimes.
• From the obtained result for the Lamb Shift it is evident that, if the geometry of the space-time is curved, particularly if it is De Sitter space then in the context of two body entangled OQS the non-adiabatic inter-atomic RCPI is purely characterised by three important contributions which are appended bellow: 1. The amplitude of the Lamb Shift is mainly characterised by the L −2 factor, which indicates the inverse square power law decay in the limiting situation L >> k, where the inter atomic distance is larger than the characteristic length scale k.
This factor is appearing in the coefficient of the cut-off dependent contribution after applying Bethe rugularisation procedure. However, in the limiting situation where the Bethe cut-off frequency is smaller than the natural frequency of the two entangled atomic system i.e. ω c << ω 0 , such contribution will not contribute in the amplitude of Lamb Shift.
• On the other hand, in case of flat space-time the amplitude of the Lamb Shift is proportional to L −1 , and the angular modulation factor D lying within the same window as mentioned above. Most importantly, in the final expression for the Lamb Shift as appearing in the case of flat space no signature of the actual origin of quantum state i.e. whether it is non-thermal or thermal, can be observable.
• Finally, we have added a short discussion regarding the Gibbons Hawking temperature and Unruh temperature that these detectors would measure when accelerating through the De Sitter background space-time. We have included this discussion regarding the equilibrium temperature of the thermal bath as it can be directly expressed in terms of the curvature of the De Sitter space and consequently the energy shift can be expressed in terms of these fundamental quantities.
The future prospects of our work are appended point-wise.
• It is completely true fact that in quantum field theory in de Sitter space (e.g. inflationary perturbations), 2-point functions encode the power spectrum while 3-point functions encode aspects of non-gaussianities. But to explain this quantum entanglement is not the necessary ingredient. Now one can consider the three atomic entanglement where one needs to consider pair wise entanglement and third atom is not entanglement. Considering this possibility one can construct ground, excited, symmetric and antisymmetric entangled states using which one can study the perturbation in static de Sitter background and after path integrating out the bath degrees of freedom (scalar field in our consideration) one can explicitly compute two, point function and the associated power spectrum and three point functions and the associated non-gaussianities in the present context. If this type computation is performed in the present context, then that will be surely very helpful to probe the underlying physics of quantum mechanical entanglement at the level of correlation functions. It is expected that the result obtained from one un-entangled atom and the entangled two/three or more atom at the level finding the quantum correlation function will be different and distinguishable. Now if one can able to test these possibilities through various observational probes, then one can say that whether quantum mechanical entanglement is at all be tested in near future or not. By following the referees suggestion regarding this issue we are now thinking of doing the same computation in near future and we believe this analysis will be helpful to explore various unknowns in this context. We have added this point in the discussion of future direction of the present work in the conclusion section 6 .
• Last but not the least, in the present context one can also study the de Sitter entropy from the gravitational entanglement in the static patch of de Sitter space-time. In this connection, one can use generalizations of Ryu-Takayanagi holographic entanglement [60,61] by following the references [62,63] and related works in this area. Also one can use, the concept of the well known ER=EPR [64], and related ideas on geometryentanglement in the context of de Sitter space. Till now we have not been used the concept and computational techniques of the gravitational entanglement to compute the de Sitter entropy. We have used the tools and techniques of quantum field theory to study the result of de Sitter entropy in various different contexts. It would be really good and interesting if the calculation can be performed in the gravitational side as well to give complete understanding and overview of the full theory 7 .
• In future we will extend our discussion for multi entangled (even and odd number of atoms) OQS to understand the connection between curvature of De Sitter space and Lamb Shift.
• Also, our plan also is to study the non-unitary dissipative time dynamics of the system from which one can compute the expression for the equilibrium temperature of the bath, which help us to check the consistency about the expression for the temperature obtained from gravity sector.
• One can also study various quantum information theoretic measure from the present multi entangled open quantum theory set up to know more about many body quantum entanglement.
• Study of quantum fluctuations and related inflationary perturbations from the OQS is not well established yet in the context of early universe cosmology. For this reason it is good to study the cosmological consequences from the correlation functions and comparison with the various observables from OQS to know about many more unknown physical facts.
• Many body localisation and eigenstate thermalisation, study of tensor networks [65][66][67] and physics of quantum chaos [18,[68][69][70][71][72][73][74]  Last but not the least, we would like to acknowledge our debt to the people belonging to the various part of the world for their generous and steady support for research in natural sciences.

A Coefficients of the Lamb shift Hamiltonian and Lindbladian
In this appendix, we will calculate the co-efficient matrix elements of the Lamb shift Hamiltonian and Lindbladian. For this purpose, we will use the Fourier transform of the two point correlation functions as given in the equation 2.18. For clarity, we re write those equations here again.
where, we define the spectral function f (ω, L/2) as: Now, the elements of co-efficient matrix H (αβ) ij of the effective Hamiltonian can be explicitly represented by the following expression: K αβ (±ω 0 )∀(α, β = 1, 2) represents the Hilbert transform of the two point function in Fourier space, which we have defined in the equation 3.33. These are given as follows: From the symmetry of Hilbert transformed correlations, it is easy to note that From equation A.5 Further using equation (A.6) in equation (A.8) we get the following simplified expression for A 1 , B 1 , A 2 and B 2 : , (A.9) , (A.10) , (A.11) . (A.12) We will calculate the above spectroscopic integrals explicitly in the next section.
Similarly, the elements of the Gorini-Kossakowski-Sudarshan-Lindblad matrix, C (αβ) ij , as appearing in the expression for the Linbladian can be expressed as: where, the quantitiesÃ (αβ) andB (αβ) for the two atomic system are defined as: Again from symmetry of Fourier transformed correlations it can be seen that (A.20)

B Calculation of Bethe regularised spectroscopic integrals
In the following subsections we explicitly compute the Bethe regularised integrals, which are very useful to compute the expressions for the energy shift from ground, excited, symmetric and antisymmetric state respectively. These integrals are appended below: Integral I : Integral II : Integral III :

B.1 Spectroscopic Integral I
In this subsection we explicitly compute the finite contribution from the following integral: It is important to note that in the large frequency range, −∞ < ω < ∞, one can further expand the integrand by taking large ω approximation as: This implies that, after taking large ω approximation the integrand of ∆ 1 becomes independent of the parameter k, which representing the surface gravity. Now, further using this approximation the integral ∆ 1 can be further simplified as: where we have decomposed the integrals into two parts, indicated by U 1 (ω 0 ) and U 2 (ω 0 ) in the parenthesis symbol. Now, here we see that in the large frequency range, −∞ < ω < ∞, we get: Now, we see here that both U 1 (ω 0 ) and U 2 (ω 0 ) gives divergent contributions in the frequency range, −∞ < ω < 0 and 0 < ω < ∞. To get the finite contributions from these integrals we introduce a cut-off regulator ω c , by following Bethe regularisation technique. After introducing this cut-off we get: Consequently, we get the following expression for the integral ∆ 1 , as given by:

B.2 Spectroscopic Integral II
In this subsection we explicitly compute the finite contribution from the following integral: It is important to note that in the large frequency range, −∞ < ω < ∞, one can further expand the integrand by taking large ω approximation as: This implies that, after taking large ω approximation the integrand of ∆ 2 becomes independent of the parameter k, which representing the surface gravity. Now, further using this approximation the integral ∆ 2 can be further simplified as: where we have decomposed the integrals into two parts, indicated by W 1 (ω 0 ) and W 2 (ω 0 ) in the parenthesis symbol. Now, here we see that in the large frequency range, −∞ < ω < ∞, we get: Now, we see here that both W 1 (ω 0 ) and W 2 (ω 0 ) gives divergent contributions in the frequency range, −∞ < ω < 0 and 0 < ω < ∞. To get the finite contributions from these integrals we introduce a cut-off regulator ω c , by following Bethe regularisation technique. After introducing this cut-off we get: Consequently, we get the following expression for the integral ∆ 2 , as given by: Now, if we further use the approximation that, ω c << ω 0 i.e. the Bethe regularised cut-off is smaller than the natural frequency of the two entangled atomic system under consideration, then we get 8 :

B.3 Spectroscopic Integral III
In this subsection we explicitly compute the finite contribution from the following integral: In the limit, ω c << ω 0 we can approximate the Taylor series expansion of the following function as: where, we define the spectral function f (ω, L/2) as: It is important to note that in the large frequency range, −∞ < ω < ∞, one can further expand the integrand by taking large ω approximation as: (B.20) This implies that, after taking large ω approximation the integrand of ∆ 3 becomes not independent of the parameter k, which representing the surface gravity. Now, further using this approximation the integral ∆ 3 can be further simplified as: where we have decomposed the integrals into two parts, indicated by Q 1 (ω 0 , k) and Q 2 (ω 0 , k) in the parenthesis symbol. Now, here we see that in the large frequency range, −∞ < ω < ∞, we get: Consequently, we get the following expression for the integral ∆ 3 , as given by: