Casimir Effect in Stochastic Semi-Classical Gravity

M.Hossein. B.Khoshnevis Sadaf Ebadi M. Reza Lahooti Reza Pirmoradian
Abstract

This article aims to examine the Casimir effect in the framework of stochastic semi-classical gravity. We commence with the semi-classical Einstein-Langevin equation, which introduces a first-order correction to the semi-classical gravity theory. Subsequently, we analyze the alteration in the Casimir force caused by this type of correction. Our results demonstrate that the corrections exhibit significant sensitivity to both the distance between the two parallel plates and the order parameter of the weak field. This finding underscores the nuanced interplay between these factors and the overall behavior of the system, providing valuable insights into the underlying dynamics.

1 Introduction

In 1948 a force of attraction between two parallel, neutral and conducting plates in a vacuum was predicted by Casimir which was named Casimir force [1, 2]. The Casimir effect itself arises from quantum fluctuations of the electromagnetic field between two uncharged, conducting plates placed very close to each other in a vacuum. These fluctuations cause an attractive force between the plates, which can be calculated using quantum field theory. The Casimir effect stands out as one of the most remarkable manifestations of non-trivial properties of the vacuum state in quantum field theory [3, 4] and this effect can be interpreted as the polarization of the vacuum by the boundary conditions or geometry. A straightforward computation of the Casimir pressure yields infinite sums, which therefore requires regularization [5, 6, 7, 8]. Several regulators may be applied, each yielding the same result; Green’s function [3, 9], stress tensor methods [10, 12, 13, 11], Krein space method [14, 15, 16, 17]. The renormalization of the Casimir effect is a fascinating topic in quantum field theory and condensed matter physics. Formally, in the context of the Casimir effect, renormalization involves the following steps; Regularization: This step involves introducing a cutoff parameter to make the integrals finite. For example, one might impose a maximum frequency or momentum to avoid divergences. Renormalization: After regularization, the physical quantities are redefined to absorb the divergences. This is done by introducing counterterms into the Hamiltonian or Lagrangian of the system. Physical Predictions: The renormalized quantities are then used to make physical predictions that are finite and meaningful. Renormalization of the Casimir effect has several important applications and implications such as nanotechnology, where understanding the Casimir effect is crucial for designing nanoscale devices, as the forces can significantly affect the behavior of components at small scales. The Casimir effect can be generalized to other types of fields and materials, providing insights into phenomena such as van der Waals forces and surface interactions [18].
The relationship between the Casimir effect and weak gravity is an intriguing area of study that bridges quantum field theory and general relativity. In the regime of weak gravity, where gravitational fields are small perturbations to flat spacetime, the modifications to vacuum fluctuations become significant [19]. By studying the Casimir effect in curved spacetime, one can investigate how general relativity impacts quantum electrodynamics, this may provide insights into the quantum aspects of gravity. Theoretical predictions suggest that gravitational waves, ripples in spacetime caused by massive accelerating objects, could induce a Casimir-like force between closely spaced mirrors [20]. This potential gravitational Casimir effect could provide a novel experimental avenue for testing quantum gravity theories and searching for gravitons, the hypothetical quantum particles that mediate gravity [21]. Such studies highlight the interplay between quantum mechanics and general relativity and open new pathways for exploring fundamental physics (see also [22]). Basically, this phenomenon could, in principle, occur for any quantized field so that if gravity has a quantum nature, then gravitational waves should also generate Casimir-like forces. In the regime of the semi-classical theory of gravity which describes a classical gravitational field coupled to a quantum matter field [23], vacuum polarization caused by the gravitational field is also referred to as the Casimir effect. The interaction between gravity and quantum fields was first explored in the context of cosmological back-reaction problems [24, 25]. In the context of quantum field theory, zero-point field fluctuations contribute to the vacuum energy. Despite being fundamental, this energy remains unobservable in laboratory settings. Nevertheless, the vacuum energy density, more precisely the vacuum expectation value of the stress-energy tensor, is incorporated into Einstein’s equation as it pertains to the vacuum state of the scalar field, thereby influencing various cosmological phenomena [26, 27] as

T^μνvacsubscriptdelimited-⟨⟩subscript^𝑇𝜇𝜈𝑣𝑎𝑐\displaystyle\langle\hat{T}_{\mu\nu}\rangle_{vac}⟨ over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ start_POSTSUBSCRIPT italic_v italic_a italic_c end_POSTSUBSCRIPT =ρvacgμνabsentsubscriptdelimited-⟨⟩𝜌𝑣𝑎𝑐subscript𝑔𝜇𝜈\displaystyle=-\langle\rho\rangle_{vac}g_{\mu\nu}= - ⟨ italic_ρ ⟩ start_POSTSUBSCRIPT italic_v italic_a italic_c end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT
Rμν12gμνR+Λgμνsubscript𝑅𝜇𝜈12subscript𝑔𝜇𝜈𝑅Λsubscript𝑔𝜇𝜈\displaystyle R_{\mu\nu}-\dfrac{1}{2}g_{\mu\nu}R+\Lambda g_{\mu\nu}italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_R + roman_Λ italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT =8πGTμνvac,absent8𝜋𝐺subscriptdelimited-⟨⟩subscript𝑇𝜇𝜈𝑣𝑎𝑐\displaystyle=8\pi G\langle T_{\mu\nu}\rangle_{vac},= 8 italic_π italic_G ⟨ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ start_POSTSUBSCRIPT italic_v italic_a italic_c end_POSTSUBSCRIPT , (1)

Rμνsubscript𝑅𝜇𝜈R_{\mu\nu}italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT, R𝑅Ritalic_R, ΛΛ\Lambdaroman_Λ and G𝐺Gitalic_G are, respectively, Ricci tensor, Ricci scalar, cosmological constant and Newton’s gravitational constant. The classical metric’s equation of motion is the semi-classical Einstein’s equation which describes the back reaction of the matter field on the spacetime; it is a generalized form of the Einstein’s equation where the energy source is supposed to be the expectation value in some quantum state of the matter stress energy tensor operator. The lowest order quantum stress energy fluctuations of the matter field as a source of classical stochastic fluctuations of the gravitational field is described by the semi-classical Einstein-Langevin equation [28].

On the other hand, stochastic gravity deals with the effects of quantum fluctuations in spacetime, particularly focusing on the gravitational back-reaction from quantum fields, which can result in phenomena such as the generation of gravitational waves or the understanding of black hole thermodynamics. Stochastic gravity starts with quantum field theory in curved spacetime [29, 30, 31, 32, 33]. In fact the physical observable that measures the stress-energy fluctuations may be considered as an advantage of the derivation of Einstein-Langevin equation and of the stochastic gravity theory [34]. In the context of stochastic gravity, Casimir force might be interpreted as part of the broader discussion of quantum field theory in curved spacetime, where the random nature of quantum fluctuations interacts with gravitational fields. The interplay between quantum forces, such as the Casimir force, and gravity in stochastic gravity frameworks provides important insights into the behavior of the quantum vacuum and its potential influence on large-scale cosmic structures. The aim of this paper is to investigate the Casimir force in the linear form in the stochastic gravity framework.

The paper is structured as follows: Sections 2 and 3 review the mathematical foundations underpinning this study, focusing on semi-classical gravity and the Casimir force. Section 4 applies these results to explore the Casimir force within the weak field limit of gravity. Finally, the paper wraps up with a discussion of the conclusions derived from the investigation.

2 Semi-classical Gravity: Mathematical Framework

Semi-classical gravity is a theoretical approach that unifies classical general relativity with quantum field theory. In this approach, gravity is treated classically through Einstein’s field equations, while other fields, such as the electromagnetic field, are quantized. The Casimir effect, which arises from quantum fluctuations of the electromagnetic field between two uncharged, closely spaced conducting plates, can be studied within this context. By examining the Casimir effect in semi-classical gravity, one can explore how quantum field fluctuations influence and are influenced by the curvature of spacetime. This interplay offers insights into the quantum aspects of gravity and the effects of quantum vacuum energy on gravitational systems. In fact Semi-classical gravity describes the interaction of gravitational field, which is assumed to be a classical field, with quantum matter fields [35]. The theory at least for the linear matter field, is mathematically well-defined and fairly well-understood [32]. The classical action for a massive real scalar ϕ(x)italic-ϕ𝑥\phi(x)italic_ϕ ( italic_x ) is given by

Sm[g,ϕ]=12d4xg[gμνμϕνϕ+(m2+ξR)ϕ2]subscript𝑆𝑚𝑔italic-ϕ12superscript𝑑4𝑥𝑔delimited-[]superscript𝑔𝜇𝜈subscript𝜇italic-ϕsubscript𝜈italic-ϕsuperscript𝑚2𝜉𝑅superscriptitalic-ϕ2S_{m}[g,\phi]=-\dfrac{1}{2}\int d^{4}x\sqrt{-g}[g^{\mu\nu}\nabla_{\mu}\phi% \nabla_{\nu}\phi+(m^{2}+\xi R)\phi^{2}]italic_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_g , italic_ϕ ] = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x square-root start_ARG - italic_g end_ARG [ italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ϕ ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ϕ + ( italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ξ italic_R ) italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] (2)

where ξ𝜉\xiitalic_ξ is called conformal coupling. The Klein-Gordon equation can be written for a real scalar field as

(m2ξR)ϕ^=0,superscript𝑚2𝜉𝑅^italic-ϕ0(\square-m^{2}-\xi R)\hat{\phi}=0,( □ - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ξ italic_R ) over^ start_ARG italic_ϕ end_ARG = 0 , (3)

note that =μμsubscript𝜇superscript𝜇\square=\nabla_{\mu}\nabla^{\mu}□ = ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∇ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT. The stress-energy tensor is given by

Tμν2detgδSmδgμν,subscript𝑇𝜇𝜈2𝑑𝑒𝑡𝑔𝛿subscript𝑆𝑚𝛿superscript𝑔𝜇𝜈T_{\mu\nu}\equiv\dfrac{-2}{\sqrt{-detg}}\dfrac{\delta S_{m}}{\delta g^{\mu\nu}},italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ≡ divide start_ARG - 2 end_ARG start_ARG square-root start_ARG - italic_d italic_e italic_t italic_g end_ARG end_ARG divide start_ARG italic_δ italic_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_δ italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT end_ARG , (4)

which leads to

Tμν[g,ϕ]=μϕνϕ12gμν(ρϕρϕ+m2ϕ2)+ξ(gμνμν+Gμν)ϕ2,subscript𝑇𝜇𝜈𝑔italic-ϕsubscript𝜇italic-ϕsubscript𝜈italic-ϕ12subscript𝑔𝜇𝜈superscript𝜌italic-ϕsubscript𝜌italic-ϕsuperscript𝑚2superscriptitalic-ϕ2𝜉subscript𝑔𝜇𝜈subscript𝜇subscript𝜈subscript𝐺𝜇𝜈superscriptitalic-ϕ2T_{\mu\nu}[g,\phi]=\nabla_{\mu}\phi\nabla_{\nu}\phi-\dfrac{1}{2}g_{\mu\nu}(% \nabla^{\rho}\phi\nabla_{\rho}\phi+m^{2}\phi^{2})+\xi(g_{\mu\nu}\square-\nabla% _{\mu}\nabla_{\nu}+G_{\mu\nu})\phi^{2},italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g , italic_ϕ ] = ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ϕ ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ϕ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( ∇ start_POSTSUPERSCRIPT italic_ρ end_POSTSUPERSCRIPT italic_ϕ ∇ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT italic_ϕ + italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_ξ ( italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT □ - ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT + italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (5)

Gμνsubscript𝐺𝜇𝜈G_{\mu\nu}italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT is the Einstein tensor. In the semi-classical theory the gravitational field equation will not be the classical equation instead, the source is supposed to be the quantum expectation value T^uv[g]delimited-⟨⟩subscript^𝑇𝑢𝑣delimited-[]𝑔\langle\hat{T}_{uv}\rangle[g]⟨ over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_u italic_v end_POSTSUBSCRIPT ⟩ [ italic_g ]. The semi-classical Einstein equation for the metric guvsubscript𝑔𝑢𝑣g_{uv}italic_g start_POSTSUBSCRIPT italic_u italic_v end_POSTSUBSCRIPT is given by

18πG(Gμν[g]+Λgμν)2(αAμν+βBμν)[g]=T^μνR[g],18𝜋𝐺subscript𝐺𝜇𝜈delimited-[]𝑔Λsubscript𝑔𝜇𝜈2𝛼subscript𝐴𝜇𝜈𝛽subscript𝐵𝜇𝜈delimited-[]𝑔delimited-⟨⟩subscriptsuperscript^𝑇𝑅𝜇𝜈delimited-[]𝑔\dfrac{1}{8\pi G}(G_{\mu\nu}[g]+\Lambda g_{\mu\nu})-2(\alpha A_{\mu\nu}+\beta B% _{\mu\nu})[g]=\langle\hat{T}^{R}_{\mu\nu}\rangle[g],divide start_ARG 1 end_ARG start_ARG 8 italic_π italic_G end_ARG ( italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g ] + roman_Λ italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) - 2 ( italic_α italic_A start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + italic_β italic_B start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) [ italic_g ] = ⟨ over^ start_ARG italic_T end_ARG start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ [ italic_g ] , (6)

T^μνR[g]delimited-⟨⟩subscriptsuperscript^𝑇𝑅𝜇𝜈delimited-[]𝑔\langle\hat{T}^{R}_{\mu\nu}\rangle[g]⟨ over^ start_ARG italic_T end_ARG start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ [ italic_g ] is the expectation value of T^μνRsubscriptsuperscript^𝑇𝑅𝜇𝜈\hat{T}^{R}_{\mu\nu}over^ start_ARG italic_T end_ARG start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT defined by

T^μνR[g](x)=T^μν[g](x)+FμνC[g](x)I^,subscriptsuperscript^𝑇𝑅𝜇𝜈delimited-[]𝑔𝑥subscript^𝑇𝜇𝜈delimited-[]𝑔𝑥subscriptsuperscript𝐹𝐶𝜇𝜈delimited-[]𝑔𝑥^𝐼\hat{T}^{R}_{\mu\nu}[g](x)=\hat{T}_{\mu\nu}[g](x)+F^{C}_{\mu\nu}[g](x)\hat{I},over^ start_ARG italic_T end_ARG start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g ] ( italic_x ) = over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g ] ( italic_x ) + italic_F start_POSTSUPERSCRIPT italic_C end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g ] ( italic_x ) over^ start_ARG italic_I end_ARG , (7)

where I^^𝐼\hat{I}over^ start_ARG italic_I end_ARG is the identity operator and FμνC[g]subscriptsuperscript𝐹𝐶𝜇𝜈delimited-[]𝑔F^{C}_{\mu\nu}[g]italic_F start_POSTSUPERSCRIPT italic_C end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT [ italic_g ] are some symmetric tensor counter terms that depend on the regulator and local function of the metric. Left hand side of the equation (6) may be derived from the following gravitational action

Sg[g]=d4xg[116πG(R2Λ)+αCμνρσCμνρσ+βR2],subscript𝑆𝑔delimited-[]𝑔superscript𝑑4𝑥𝑔delimited-[]116𝜋𝐺𝑅2Λ𝛼subscript𝐶𝜇𝜈𝜌𝜎superscript𝐶𝜇𝜈𝜌𝜎𝛽superscript𝑅2S_{g}[g]=\int d^{4}x\sqrt{-g}[\dfrac{1}{16\pi G}(R-2\Lambda)+\alpha C_{\mu\nu% \rho\sigma}C^{\mu\nu\rho\sigma}+\beta R^{2}],italic_S start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT [ italic_g ] = ∫ italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x square-root start_ARG - italic_g end_ARG [ divide start_ARG 1 end_ARG start_ARG 16 italic_π italic_G end_ARG ( italic_R - 2 roman_Λ ) + italic_α italic_C start_POSTSUBSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUPERSCRIPT + italic_β italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , (8)

where Cμνρσsubscript𝐶𝜇𝜈𝜌𝜎C_{\mu\nu\rho\sigma}italic_C start_POSTSUBSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUBSCRIPT is the Weyl tensor. Aμνsubscript𝐴𝜇𝜈A_{\mu\nu}italic_A start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT and Bμνsubscript𝐵𝜇𝜈B_{\mu\nu}italic_B start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT are tensors that are derived from derivatives of the functional with respect to the metric of the quadratic terms in the curvature in (8) which are explicitly given by

Aμνsubscript𝐴𝜇𝜈\displaystyle A_{\mu\nu}italic_A start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT =12gμνCτϵρσCτϵρσ2RμρστRνρστ+4RμρRνρ23RRμν2Rμν+23μνR+13gμνR,absent12subscript𝑔𝜇𝜈subscript𝐶𝜏italic-ϵ𝜌𝜎superscript𝐶𝜏italic-ϵ𝜌𝜎2subscript𝑅𝜇𝜌𝜎𝜏superscriptsubscript𝑅𝜈𝜌𝜎𝜏4subscript𝑅𝜇𝜌subscriptsuperscript𝑅𝜌𝜈23𝑅subscript𝑅𝜇𝜈2subscript𝑅𝜇𝜈23subscript𝜇subscript𝜈𝑅13subscript𝑔𝜇𝜈𝑅\displaystyle=\dfrac{1}{2}g_{\mu\nu}C_{\tau\epsilon\rho\sigma}C^{\tau\epsilon% \rho\sigma}-2R_{\mu\rho\sigma\tau}R_{\nu}^{\rho\sigma\tau}+4R_{\mu\rho}R^{\rho% }_{\nu}-\dfrac{2}{3}RR_{\mu\nu}-2\square R_{\mu\nu}+\dfrac{2}{3}\nabla_{\mu}% \nabla_{\nu}R+\dfrac{1}{3}g_{\mu\nu}\square R,= divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_τ italic_ϵ italic_ρ italic_σ end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_τ italic_ϵ italic_ρ italic_σ end_POSTSUPERSCRIPT - 2 italic_R start_POSTSUBSCRIPT italic_μ italic_ρ italic_σ italic_τ end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ρ italic_σ italic_τ end_POSTSUPERSCRIPT + 4 italic_R start_POSTSUBSCRIPT italic_μ italic_ρ end_POSTSUBSCRIPT italic_R start_POSTSUPERSCRIPT italic_ρ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT - divide start_ARG 2 end_ARG start_ARG 3 end_ARG italic_R italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - 2 □ italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + divide start_ARG 2 end_ARG start_ARG 3 end_ARG ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_R + divide start_ARG 1 end_ARG start_ARG 3 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT □ italic_R ,
Bμνsubscript𝐵𝜇𝜈\displaystyle B_{\mu\nu}italic_B start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT =12gμνR22RRμν+2μνR2gμνR,absent12subscript𝑔𝜇𝜈superscript𝑅22𝑅subscript𝑅𝜇𝜈2subscript𝜇subscript𝜈𝑅2subscript𝑔𝜇𝜈𝑅\displaystyle=\dfrac{1}{2}g_{\mu\nu}R^{2}-2RR_{\mu\nu}+2\nabla_{\mu}\nabla_{% \nu}R-2g_{\mu\nu}\square R,= divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_R italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + 2 ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_R - 2 italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT □ italic_R , (9)

Rμνρσsubscript𝑅𝜇𝜈𝜌𝜎R_{\mu\nu\rho\sigma}italic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUBSCRIPT and Rμνsubscript𝑅𝜇𝜈R_{\mu\nu}italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT are Riemann and Ricci tensors respectively.

Stochastic gravity is an advanced framework that blends quantum field theory with general relativity to investigate the impact of quantum fluctuations on gravitational phenomena. Within this framework, zero-point field fluctuations contribute to the vacuum energy, which remains unobservable in laboratory settings. By rewriting the equation (1) as the form of the Einstein-Langevin equation which is central to stochastic gravity given by

Gμν+Λgμν=8πGTμν+8πGd4xNμν(x,x)subscript𝐺𝜇𝜈Λsubscript𝑔𝜇𝜈8𝜋𝐺delimited-⟨⟩subscript𝑇𝜇𝜈8𝜋𝐺superscript𝑑4superscript𝑥subscript𝑁𝜇𝜈𝑥superscript𝑥G_{\mu\nu}+\Lambda g_{\mu\nu}=8\pi G\langle T_{\mu\nu}\rangle+\sqrt{8\pi G}% \int d^{4}x^{\prime}\,N_{\mu\nu}(x,x^{\prime})italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + roman_Λ italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = 8 italic_π italic_G ⟨ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ + square-root start_ARG 8 italic_π italic_G end_ARG ∫ italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT )

where Tμνdelimited-⟨⟩subscript𝑇𝜇𝜈\langle T_{\mu\nu}\rangle⟨ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ is the expectation value of the stress-energy tensor. The term Nμν(x,x)subscript𝑁𝜇𝜈𝑥superscript𝑥N_{\mu\nu}(x,x^{\prime})italic_N start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) is the noise kernel, defined as:

Nμν(x,x)=δTμν(x)δTμν(x)subscript𝑁𝜇𝜈𝑥superscript𝑥delimited-⟨⟩𝛿subscript𝑇𝜇𝜈𝑥𝛿subscript𝑇𝜇𝜈superscript𝑥N_{\mu\nu}(x,x^{\prime})=\langle\delta T_{\mu\nu}(x)\delta T_{\mu\nu}(x^{% \prime})\rangleitalic_N start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = ⟨ italic_δ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x ) italic_δ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ (10)

which captures the vacuum fluctuations in the stress-energy tensor. Metric perturbations around Minkowski spacetime can be represented as hμν=gμνημνsubscript𝜇𝜈subscript𝑔𝜇𝜈subscript𝜂𝜇𝜈h_{\mu\nu}=g_{\mu\nu}-\eta_{\mu\nu}italic_h start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - italic_η start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT. The fluctuation-dissipation relation,

δTμν(x)δTμν(x)=12(δTμνδgμν(x))δ4(xx),delimited-⟨⟩𝛿subscript𝑇𝜇𝜈𝑥𝛿subscript𝑇𝜇𝜈superscript𝑥12𝛿delimited-⟨⟩subscript𝑇𝜇𝜈𝛿superscript𝑔𝜇𝜈𝑥superscript𝛿4𝑥superscript𝑥\langle\delta T_{\mu\nu}(x)\delta T_{\mu\nu}(x^{\prime})\rangle=\frac{1}{2}% \left(\frac{\delta\langle T_{\mu\nu}\rangle}{\delta g^{\mu\nu}(x)}\right)% \delta^{4}(x-x^{\prime}),⟨ italic_δ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x ) italic_δ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ( italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( divide start_ARG italic_δ ⟨ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ⟩ end_ARG start_ARG italic_δ italic_g start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT ( italic_x ) end_ARG ) italic_δ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x - italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (11)

provides a connection between the fluctuations in the stress-energy tensor and the dissipative dynamics of metric perturbations. These are the key formulas that illustrate how stochastic gravity integrates quantum effects into classical gravitational theory, providing a robust framework for exploring quantum fluctuations’ influence on gravity [36, 37]. Stochastic semi-classical gravity aims to extend the semi-classical theory by accounting for these fluctuations in a self-consistent way. We use this approach to study the Casimir effect.

3 Casimir Effect in Minkowski Space-time

The Casimir effect has been studied for different fields in different background geometries [6, 7, 30, 38, 2]. For two large, perfectly conducting plates with surface area A𝐴Aitalic_A, separated by a distance L𝐿Litalic_L, where A>>Lmuch-greater-than𝐴𝐿\sqrt{A}>>Lsquare-root start_ARG italic_A end_ARG > > italic_L, allowing us to neglect edge contributions. Plates are in xy𝑥𝑦x-yitalic_x - italic_y plane at Z=0𝑍0Z=0italic_Z = 0 and Z=L𝑍𝐿Z=Litalic_Z = italic_L. The attractive force per unit area, i.e. the pressure between two infinitely large, neutral parallel planes made of an ideal metal at zero temperature with Dirichlet boundary condition for a massless minimally coupled field in the Minkowski space-time, where the Casimir energy per unit area is C=π2720L3subscript𝐶superscript𝜋2720superscript𝐿3\mathcal{E}_{C}=-\dfrac{\pi^{2}}{720L^{3}}caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT = - divide start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 720 italic_L start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG, is given by [39]

P(L)Mink=π2240L4.𝑃subscript𝐿𝑀𝑖𝑛𝑘superscript𝜋2240superscript𝐿4P(L)_{Mink}=\dfrac{-\pi^{2}}{240L^{4}}.italic_P ( italic_L ) start_POSTSUBSCRIPT italic_M italic_i italic_n italic_k end_POSTSUBSCRIPT = divide start_ARG - italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 240 italic_L start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG . (12)

For a massless conformally coupled scalar field in the flat space-time, the stress-energy-momentum tensor in Minkowski space-time in the massless conformal coupling is given by

Tμν=23μϕνϕ+16ημναϕαϕ13ϕμνϕ.subscript𝑇𝜇𝜈23subscript𝜇italic-ϕsubscript𝜈italic-ϕ16subscript𝜂𝜇𝜈subscript𝛼italic-ϕsuperscript𝛼italic-ϕ13italic-ϕsubscript𝜇subscript𝜈italic-ϕT_{\mu\nu}=\dfrac{2}{3}\partial_{\mu}\phi\partial_{\nu}\phi+\dfrac{1}{6}\eta_{% \mu\nu}\partial_{\alpha}\phi\partial^{\alpha}\phi-\dfrac{1}{3}\phi\partial_{% \mu}\partial_{\nu}\phi.italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = divide start_ARG 2 end_ARG start_ARG 3 end_ARG ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ϕ ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ϕ + divide start_ARG 1 end_ARG start_ARG 6 end_ARG italic_η start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT italic_ϕ ∂ start_POSTSUPERSCRIPT italic_α end_POSTSUPERSCRIPT italic_ϕ - divide start_ARG 1 end_ARG start_ARG 3 end_ARG italic_ϕ ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ϕ . (13)

The vacuum energy density of the scalar field in the presence of boundaries is the mean value of the 00-component of the energy-momentum tensor in the vacuum state. In the conformally coupled case it is given by

0|T^00|0=560|0ϕ^0ϕ^|0160|iϕ^iϕ^|0130|ϕ^02ϕ^|0quantum-operator-product0subscript^𝑇00056quantum-operator-product0subscript0^italic-ϕsubscript0^italic-ϕ016quantum-operator-product0subscript𝑖^italic-ϕsubscript𝑖^italic-ϕ013quantum-operator-product0^italic-ϕsubscriptsuperscript20^italic-ϕ0\langle 0|\hat{T}_{00}|0\rangle=\dfrac{5}{6}\langle 0|\partial_{0}\hat{\phi}% \partial_{0}\hat{\phi}|0\rangle-\dfrac{1}{6}\langle 0|\partial_{i}\hat{\phi}% \partial_{i}\hat{\phi}|0\rangle-\dfrac{1}{3}\langle 0|\hat{\phi}\partial^{2}_{% 0}\hat{\phi}|0\rangle⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ = divide start_ARG 5 end_ARG start_ARG 6 end_ARG ⟨ 0 | ∂ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG ∂ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG | 0 ⟩ - divide start_ARG 1 end_ARG start_ARG 6 end_ARG ⟨ 0 | ∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG ∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG | 0 ⟩ - divide start_ARG 1 end_ARG start_ARG 3 end_ARG ⟨ 0 | over^ start_ARG italic_ϕ end_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG | 0 ⟩ (14)

it can be easily shown that

i(ϕ^iϕ^)=iϕ^iϕ^subscript𝑖^italic-ϕsubscript𝑖^italic-ϕsubscript𝑖^italic-ϕsubscript𝑖^italic-ϕ\partial_{i}(\hat{\phi}\partial_{i}\hat{\phi})=-\partial_{i}\hat{\phi}\partial% _{i}\hat{\phi}∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( over^ start_ARG italic_ϕ end_ARG ∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG ) = - ∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG ∂ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over^ start_ARG italic_ϕ end_ARG (15)

thus, the expectation value of the energy-momentum tensor (00-component) leads to

0|T^00|0=[560x0y+16ixiy]0|ϕ^(x)ϕ^(y)|0.quantum-operator-product0subscript^𝑇000delimited-[]56subscriptsuperscript𝑥0subscriptsuperscript𝑦016subscriptsuperscript𝑥𝑖subscriptsuperscript𝑦𝑖quantum-operator-product0^italic-ϕ𝑥^italic-ϕ𝑦0\langle 0|\hat{T}_{00}|0\rangle=\bigg{[}\dfrac{5}{6}\partial^{x}_{0}\partial^{% y}_{0}+\dfrac{1}{6}\partial^{x}_{i}\partial^{y}_{i}\bigg{]}\langle 0|\hat{\phi% }(x)\hat{\phi}(y)|0\rangle.⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ = [ divide start_ARG 5 end_ARG start_ARG 6 end_ARG ∂ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∂ start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 6 end_ARG ∂ start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∂ start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] ⟨ 0 | over^ start_ARG italic_ϕ end_ARG ( italic_x ) over^ start_ARG italic_ϕ end_ARG ( italic_y ) | 0 ⟩ . (16)

To compute the above relation we use the Feynman propagator is defined by

iGF(x,x)𝑖subscript𝐺𝐹𝑥superscript𝑥\displaystyle iG_{F}(x,x^{{}^{\prime}})italic_i italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) =0|T(ϕ^(x)ϕ^(x)|0\displaystyle=\langle 0|T(\hat{\phi}(x)\hat{\phi}(x^{{}^{\prime}})|0\rangle= ⟨ 0 | italic_T ( over^ start_ARG italic_ϕ end_ARG ( italic_x ) over^ start_ARG italic_ϕ end_ARG ( italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) | 0 ⟩
=θ(tt)G+(x,x)+θ(tt)G(x,x)absent𝜃𝑡superscript𝑡superscript𝐺𝑥superscript𝑥𝜃superscript𝑡𝑡superscript𝐺𝑥superscript𝑥\displaystyle=\theta(t-t^{{}^{\prime}})G^{+}(x,x^{{}^{\prime}})+\theta(t^{{}^{% \prime}}-t)G^{-}(x,x^{{}^{\prime}})= italic_θ ( italic_t - italic_t start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) italic_G start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) + italic_θ ( italic_t start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT - italic_t ) italic_G start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) (17)

where the Wightman functions G±superscript𝐺plus-or-minusG^{\pm}italic_G start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are defined as follows

G+(x,x)superscript𝐺𝑥superscript𝑥\displaystyle G^{+}(x,x^{{}^{\prime}})italic_G start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) =0|ϕ^(x)ϕ^(x)|0,absentquantum-operator-product0^italic-ϕ𝑥^italic-ϕsuperscript𝑥0\displaystyle=\langle 0|\hat{\phi}(x)\hat{\phi}(x^{{}^{\prime}})|0\rangle,= ⟨ 0 | over^ start_ARG italic_ϕ end_ARG ( italic_x ) over^ start_ARG italic_ϕ end_ARG ( italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) | 0 ⟩ ,
G(x,x)superscript𝐺𝑥superscript𝑥\displaystyle G^{-}(x,x^{{}^{\prime}})italic_G start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) =0|ϕ^(x)ϕ^(x)|0.absentquantum-operator-product0^italic-ϕsuperscript𝑥^italic-ϕ𝑥0\displaystyle=\langle 0|\hat{\phi}(x^{{}^{\prime}})\hat{\phi}(x)|0\rangle.= ⟨ 0 | over^ start_ARG italic_ϕ end_ARG ( italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) over^ start_ARG italic_ϕ end_ARG ( italic_x ) | 0 ⟩ . (18)

θ(t)𝜃𝑡\theta(t)italic_θ ( italic_t ) is the Heaviside step function, it is then straightforward to verify that

(t2i2)GF(x,x)=iδ4(xx).subscriptsuperscript2𝑡subscriptsuperscript2𝑖subscript𝐺𝐹𝑥superscript𝑥𝑖superscript𝛿4𝑥superscript𝑥\displaystyle(\partial^{2}_{t}-\partial^{2}_{i})G_{F}(x,x^{{}^{\prime}})=i% \delta^{4}(x-x^{{}^{\prime}}).( ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) = italic_i italic_δ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x - italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) . (19)

The general form of all Green’s functions in which we use the mode decomposition of ϕitalic-ϕ\phiitalic_ϕ can be demonstrated as

𝒢(x,x)=dnk(2π)nexp[ik.(xx)ik0(tt)](k0)2|k|2m2.\displaystyle\mathcal{G}(x,x^{{}^{\prime}})=\int\dfrac{d^{n}k}{(2\pi)^{n}}% \dfrac{exp[i\vec{k}.(\vec{x}-\vec{x}^{{}^{\prime}})-ik^{0}(t-t^{{}^{\prime}})]% }{(k^{0})^{2}-|\vec{k}|^{2}-m^{2}}.caligraphic_G ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) = ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_k end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_e italic_x italic_p [ italic_i over→ start_ARG italic_k end_ARG . ( over→ start_ARG italic_x end_ARG - over→ start_ARG italic_x end_ARG start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) - italic_i italic_k start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) ] end_ARG start_ARG ( italic_k start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - | over→ start_ARG italic_k end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (20)

where for Feynman propagator, we would shift the poles at k0=±(|k|2+m2)12superscript𝑘0plus-or-minussuperscriptsuperscript𝑘2superscript𝑚212k^{0}=\pm(|k|^{2}+m^{2})^{1\over{2}}italic_k start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT = ± ( | italic_k | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT off the real axis by the replacement m2m2iϵsuperscript𝑚2superscript𝑚2𝑖italic-ϵm^{2}\rightarrow m^{2}-i\epsilonitalic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT → italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_i italic_ϵ to recover GF(x,x)subscript𝐺𝐹𝑥superscript𝑥G_{F}(x,x^{{}^{\prime}})italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ). As a consequence, in our case, it becomes

GF(x,x)=14π21(xx)2.subscript𝐺𝐹𝑥superscript𝑥14superscript𝜋21superscript𝑥superscript𝑥2G_{F}(x,x^{{}^{\prime}})=\dfrac{1}{4\pi^{2}}\dfrac{1}{(x-x^{{}^{\prime}})^{2}}.italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( italic_x , italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) = divide start_ARG 1 end_ARG start_ARG 4 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( italic_x - italic_x start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (21)

The Casimir energy density arises from the quantum fluctuations of the electromagnetic field in the vacuum, influenced by the boundary conditions imposed by conducting surfaces or materials. On an interval, such as between two parallel conducting plates, these boundary conditions restrict the possible modes of the field, leading to a modification of the zero-point energy. The resulting Casimir energy density, which can be attractive or repulsive depending on the configuration, manifests as a measurable force between the plates, often referred to as the Casimir effect. This energy density is inversely proportional to the fourth power of the distance between the plates, highlighting its significance at microscopic scales. To calculate the Casimir force from the energy density, integrate the energy density over the spatial interval to obtain the total Casimir energy. Then, derive the Casimir force by taking the negative gradient (derivative) of the total Casimir energy with respect to the distance between the plates. As a result, by using the Eq.(16), the Casimir energy density on the interval 0<z<L0𝑧𝐿0<z<L0 < italic_z < italic_L can be obtained using the local version of the bellow equation

ϵ=0|T^00|0L0|T^00|0=italic-ϵsuperscriptquantum-operator-product0subscript^𝑇000𝐿superscriptquantum-operator-product0subscript^𝑇000absent\displaystyle\epsilon=\langle 0|\hat{T}_{00}|0\rangle^{L}-\langle 0|\hat{T}_{0% 0}|0\rangle^{\infty}=italic_ϵ = ⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT - ⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT = 132π2n01(nL)4132superscript𝜋2subscript𝑛01superscript𝑛𝐿4\displaystyle-\dfrac{1}{32\pi^{2}}\sum_{n\neq 0}\dfrac{1}{(nL)^{4}}- divide start_ARG 1 end_ARG start_ARG 32 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_n ≠ 0 end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG ( italic_n italic_L ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG
=\displaystyle== 116π2L4n=11n4=116π2L4ζ(4).116superscript𝜋2superscript𝐿4subscriptsuperscript𝑛11superscript𝑛4116superscript𝜋2superscript𝐿4𝜁4\displaystyle-\dfrac{1}{16\pi^{2}L^{4}}\sum^{\infty}_{n=1}\dfrac{1}{n^{4}}=-% \dfrac{1}{16\pi^{2}L^{4}}\zeta(4).- divide start_ARG 1 end_ARG start_ARG 16 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_L start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG italic_n start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG = - divide start_ARG 1 end_ARG start_ARG 16 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_L start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG italic_ζ ( 4 ) . (22)

where, the zeta function is given by

ζ(4)=n=11n4=π490.𝜁4subscriptsuperscript𝑛11superscript𝑛4superscript𝜋490\zeta(4)=\sum^{\infty}_{n=1}\dfrac{1}{n^{4}}=\dfrac{\pi^{4}}{90}.italic_ζ ( 4 ) = ∑ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG italic_n start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG = divide start_ARG italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 90 end_ARG . (23)

Thus, one obtains

ϵ=0|T^00|0L0|T^00|0=π21440L4.italic-ϵsuperscriptquantum-operator-product0subscript^𝑇000𝐿superscriptquantum-operator-product0subscript^𝑇000superscript𝜋21440superscript𝐿4\epsilon=\langle 0|\hat{T}_{00}|0\rangle^{L}-\langle 0|\hat{T}_{00}|0\rangle^{% \infty}=-\dfrac{\pi^{2}}{1440L^{4}}.italic_ϵ = ⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT - ⟨ 0 | over^ start_ARG italic_T end_ARG start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT | 0 ⟩ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT = - divide start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1440 italic_L start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG . (24)
Refer to caption
Figure 1: Casimir energy density as function of distance L

Now, we consider a factor of 2 to take into account the two polarizations of the electromagnetic field. Hence, the Casimir energy per unit area becomes as follows

C=π2720L3,subscript𝐶superscript𝜋2720superscript𝐿3\displaystyle\mathcal{E}_{C}=-\dfrac{\pi^{2}}{720L^{3}},caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT = - divide start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 720 italic_L start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG , (25)

and the Casimir force per unit area is given by

P(L)=π2240L4.𝑃𝐿superscript𝜋2240superscript𝐿4P(L)=-\dfrac{\pi^{2}}{240L^{4}}.italic_P ( italic_L ) = - divide start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 240 italic_L start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG . (26)
Refer to caption
Figure 2: Casimir force per unit area as function of distance L
Refer to caption
Figure 3: Casimir energy per unit area as function of distance L

4 Gravity and Casimir Force: Weak Field Limit

As mentioned, the Casimir effect, a quantum phenomenon resulting from vacuum fluctuations of the electromagnetic field, can be significantly influenced by the curvature of spacetime. In the context of curved space, the boundaries or geometries that constrain the field, such as conducting plates, experience modified quantum fluctuations due to the presence of gravitational fields. This leads to changes in the Casimir force, which can be either enhanced or reduced depending on the curvature. Studying the Casimir effect in curved space offers a unique window into understanding how quantum field theory and general relativity interact, providing insights into the quantum nature of gravity and the behavior of quantum fields in non-flat spacetime geometries. This intersection of quantum mechanics and general relativity helps advance our knowledge of fundamental physics and could potentially inform theories of quantum gravity. In a general curved space-time, we have shown μTμν=0subscript𝜇superscript𝑇𝜇𝜈0\nabla_{\mu}T^{\mu\nu}=0∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT = 0, and in a space-time having special symmetries or isometries (ς𝜍\varsigmaitalic_ς generates an isometry ςgμν=μςν+νςμ=0subscript𝜍subscript𝑔𝜇𝜈subscript𝜇subscript𝜍𝜈subscript𝜈subscript𝜍𝜇0\mathcal{L}_{\varsigma}g_{\mu\nu}=\nabla_{\mu}\varsigma_{\nu}+\nabla_{\nu}% \varsigma_{\mu}=0caligraphic_L start_POSTSUBSCRIPT italic_ς end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = ∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT + ∇ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = 0) as implied by the existence of one or more killing vector fields ςμsuperscript𝜍𝜇\varsigma^{\mu}italic_ς start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT, we have μ(Tμνςν)=0subscript𝜇superscript𝑇𝜇𝜈subscript𝜍𝜈0\nabla_{\mu}(T^{\mu\nu}\varsigma_{\nu})=0∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ( italic_T start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT italic_ς start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) = 0 and 𝒫ς𝑑VxTν0ςνsubscript𝒫𝜍differential-dsubscript𝑉𝑥subscriptsuperscript𝑇0𝜈superscript𝜍𝜈\mathcal{P}_{\varsigma}\equiv\int dV_{x}T^{0}_{\nu}\varsigma^{\nu}caligraphic_P start_POSTSUBSCRIPT italic_ς end_POSTSUBSCRIPT ≡ ∫ italic_d italic_V start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ς start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT is constant [40], where dVx=dn1x|g|12𝑑subscript𝑉𝑥superscript𝑑𝑛1𝑥superscript𝑔12dV_{x}=d^{n-1}x|g|^{1\over{2}}italic_d italic_V start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = italic_d start_POSTSUPERSCRIPT italic_n - 1 end_POSTSUPERSCRIPT italic_x | italic_g | start_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT.
Stochastic gravity, on the other hand, extends the integration of quantum field theory and general relativity by incorporating quantum fluctuations into the gravitational field equations through the Einstein-Langevin equation. This advanced framework considers the noise kernel, representing stress-energy tensor fluctuations, to explore gravitational phenomena such as spacetime stability, structure formation, and the backreaction of quantum effects on black holes. Combining these two concepts, the Casimir effect demonstrates the observable consequences of quantum fluctuations, while stochastic gravity provides a theoretical foundation for understanding their impact on the curvature of spacetime. Here, let us consider a Casimir apparatus of parallel plates separated by a distance L𝐿Litalic_L, with transverse dimensions a>>Lmuch-greater-than𝑎𝐿a>>Litalic_a > > italic_L. Brown and Maclay [41] showed that, in the z-direction (the direction of gravity), the vacuum expectation value stress-energy of a conformally invariant scalar field is given by

Tμν=CLdiag(1,1,1,3)delimited-⟨⟩superscript𝑇𝜇𝜈subscript𝐶𝐿𝑑𝑖𝑎𝑔1113\displaystyle\langle T^{\mu\nu}\rangle=\dfrac{\mathcal{E}_{C}}{L}diag(1,-1,1,3)⟨ italic_T start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT ⟩ = divide start_ARG caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG italic_d italic_i italic_a italic_g ( 1 , - 1 , 1 , 3 ) (27)

The Cartesian coordinates associated with the Casimir apparatus are (ξ~,η~,χ~~𝜉~𝜂~𝜒\tilde{\xi},\tilde{\eta},\tilde{\chi}over~ start_ARG italic_ξ end_ARG , over~ start_ARG italic_η end_ARG , over~ start_ARG italic_χ end_ARG) where ξ~~𝜉\tilde{\xi}over~ start_ARG italic_ξ end_ARG is normal to the plate, while η~~𝜂\tilde{\eta}over~ start_ARG italic_η end_ARG and χ~~𝜒\tilde{\chi}over~ start_ARG italic_χ end_ARG are parallel to the plate. The apparatus is oriented at an angle α𝛼\alphaitalic_α with respect to the direction of gravity. Where

z=ξ~cosα+η~sinα,y=η~cosαξ~sinα,x=χ~.formulae-sequence𝑧~𝜉𝛼~𝜂𝛼formulae-sequence𝑦~𝜂𝛼~𝜉𝛼𝑥~𝜒z=\tilde{\xi}\cos\alpha+\tilde{\eta}\sin\alpha,\hskip 11.38092pty=\tilde{\eta}% \cos\alpha-\tilde{\xi}\sin\alpha,\hskip 11.38092ptx=\tilde{\chi}.italic_z = over~ start_ARG italic_ξ end_ARG roman_cos italic_α + over~ start_ARG italic_η end_ARG roman_sin italic_α , italic_y = over~ start_ARG italic_η end_ARG roman_cos italic_α - over~ start_ARG italic_ξ end_ARG roman_sin italic_α , italic_x = over~ start_ARG italic_χ end_ARG . (28)

Let the center of the apparatus be located at (ξ~=ξ~0,η~=0,χ~=0formulae-sequence~𝜉subscript~𝜉0formulae-sequence~𝜂0~𝜒0\tilde{\xi}=\tilde{\xi}_{0},\tilde{\eta}=0,\tilde{\chi}=0over~ start_ARG italic_ξ end_ARG = over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , over~ start_ARG italic_η end_ARG = 0 , over~ start_ARG italic_χ end_ARG = 0). Gauge transformation of hμνsubscript𝜇𝜈h_{\mu\nu}italic_h start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT for stochastic gravity obeys the following equation

hμνhμν+μςν+νςμ.subscript𝜇𝜈subscript𝜇𝜈subscript𝜇subscript𝜍𝜈subscript𝜈subscript𝜍𝜇h_{\mu\nu}\rightarrow h_{\mu\nu}+\partial_{\mu}\varsigma_{\nu}+\partial_{\nu}% \varsigma_{\mu}.italic_h start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT → italic_h start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT + ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT . (29)

Thus, the action changes as follows

ΔW=2(dx)ςνμTμνΔ𝑊2𝑑𝑥subscript𝜍𝜈subscript𝜇superscript𝑇𝜇𝜈\displaystyle\Delta W=-2\int(dx)\varsigma_{\nu}\partial_{\mu}T^{\mu\nu}roman_Δ italic_W = - 2 ∫ ( italic_d italic_x ) italic_ς start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT (30)

To study the effect of gravity on the Casimir effect, we use Fermi coordinates. A Fermi coordinates system is the general relativistic generalization of an inertial coordinate frame. The ”constant field metric” is simply the Fermi coordinate metric for gravitating body [42]

ds2=(1+2gz)dt2+d𝐫𝟐.\displaystyle ds^{2}=-(1+2gz)dt^{2}+d\bf{r}^{{}^{\prime}2}.italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - ( 1 + 2 italic_g italic_z ) italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_d bold_r start_POSTSUPERSCRIPT start_FLOATSUPERSCRIPT ′ end_FLOATSUPERSCRIPT bold_2 end_POSTSUPERSCRIPT . (31)

We compute the additional gravitational energy, in terms of the gauge field ςμsubscript𝜍𝜇\varsigma_{\mu}italic_ς start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT, which carries us from isotropic coordinates to Fermi coordinates hμνF=hμνI+μςν+νςμsubscriptsuperscript𝐹𝜇𝜈subscriptsuperscript𝐼𝜇𝜈subscript𝜇subscript𝜍𝜈subscript𝜈subscript𝜍𝜇h^{F}_{\mu\nu}=h^{I}_{\mu\nu}+\partial_{\mu}\varsigma_{\nu}+\partial_{\nu}% \varsigma_{\mu}italic_h start_POSTSUPERSCRIPT italic_F end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_h start_POSTSUPERSCRIPT italic_I end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + ∂ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT + ∂ start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_ς start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT. Therefore, we have

h00I=gz,hijI=gzδij,h00F=gz,hijF=0.formulae-sequencesubscriptsuperscript𝐼00𝑔𝑧formulae-sequencesubscriptsuperscript𝐼𝑖𝑗𝑔𝑧subscript𝛿𝑖𝑗formulae-sequencesubscriptsuperscript𝐹00𝑔𝑧subscriptsuperscript𝐹𝑖𝑗0h^{I}_{00}=-gz,\hskip 11.38092pth^{I}_{ij}=-gz\delta_{ij},\hskip 11.38092pth^{% F}_{00}=-gz,\hskip 11.38092pth^{F}_{ij}=0.italic_h start_POSTSUPERSCRIPT italic_I end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT = - italic_g italic_z , italic_h start_POSTSUPERSCRIPT italic_I end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = - italic_g italic_z italic_δ start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT , italic_h start_POSTSUPERSCRIPT italic_F end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT = - italic_g italic_z , italic_h start_POSTSUPERSCRIPT italic_F end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = 0 . (32)

With the proper definition of gauge field for the change in the energy obtained from ΔWΔ𝑊\Delta Wroman_Δ italic_W, finally we get

ΔEg=6CLa2a2𝑑η~a2a2𝑑χ~14gcosα(2ξ~0L)2CLξ~0L2ξ~0+L2𝑑ξ~a2a2𝑑χ~12gcosα(a)ξ~2CLξ~0L2ξ~0+L2𝑑ξ~a2a2𝑑η~12g(ξ~cosα+η~sinα)(a)=AgCξ~0cosα=AgCz0Δsubscript𝐸𝑔6subscript𝐶𝐿subscriptsuperscript𝑎2𝑎2differential-d~𝜂subscriptsuperscript𝑎2𝑎2differential-d~𝜒14𝑔𝛼2subscript~𝜉0𝐿2subscript𝐶𝐿subscriptsuperscriptsubscript~𝜉0𝐿2subscript~𝜉0𝐿2differential-d~𝜉subscriptsuperscript𝑎2𝑎2differential-d~𝜒12𝑔𝛼𝑎~𝜉2subscript𝐶𝐿subscriptsuperscriptsubscript~𝜉0𝐿2subscript~𝜉0𝐿2differential-d~𝜉subscriptsuperscript𝑎2𝑎2differential-d~𝜂12𝑔~𝜉𝛼~𝜂𝛼𝑎𝐴𝑔subscript𝐶subscript~𝜉0𝛼𝐴𝑔subscript𝐶subscript𝑧0\displaystyle\begin{split}\Delta E_{g}&=\dfrac{6\mathcal{E}_{C}}{L}\int^{a% \over{2}}_{-{a\over{2}}}d\tilde{\eta}\int^{a\over{2}}_{-{a\over{2}}}d\tilde{% \chi}\dfrac{1}{4}g\cos\alpha(-2\tilde{\xi}_{0}L)-\dfrac{2\mathcal{E}_{C}}{L}% \int^{\tilde{\xi}_{0}+{L\over{2}}}_{\tilde{\xi}_{0}-{L\over{2}}}d\tilde{\xi}% \int^{a\over{2}}_{-{a\over{2}}}d\tilde{\chi}\dfrac{1}{2}g\cos\alpha(-a)\tilde{% \xi}\\ &-\dfrac{2\mathcal{E}_{C}}{L}\int^{\tilde{\xi}_{0}+{L\over{2}}}_{\tilde{\xi}_{% 0}-{L\over{2}}}d\tilde{\xi}\int^{a\over{2}}_{-{a\over{2}}}d\tilde{\eta}\dfrac{% 1}{2}g(\tilde{\xi}\cos\alpha+\tilde{\eta}\sin\alpha)(-a)=-Ag\mathcal{E}_{C}% \tilde{\xi}_{0}\cos\alpha=-Ag\mathcal{E}_{C}z_{0}\end{split}start_ROW start_CELL roman_Δ italic_E start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_CELL start_CELL = divide start_ARG 6 caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ∫ start_POSTSUPERSCRIPT divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_η end_ARG ∫ start_POSTSUPERSCRIPT divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_χ end_ARG divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_g roman_cos italic_α ( - 2 over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_L ) - divide start_ARG 2 caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ∫ start_POSTSUPERSCRIPT over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG italic_L end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG italic_L end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_ξ end_ARG ∫ start_POSTSUPERSCRIPT divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_χ end_ARG divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g roman_cos italic_α ( - italic_a ) over~ start_ARG italic_ξ end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - divide start_ARG 2 caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ∫ start_POSTSUPERSCRIPT over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG italic_L end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG italic_L end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_ξ end_ARG ∫ start_POSTSUPERSCRIPT divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - divide start_ARG italic_a end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT italic_d over~ start_ARG italic_η end_ARG divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g ( over~ start_ARG italic_ξ end_ARG roman_cos italic_α + over~ start_ARG italic_η end_ARG roman_sin italic_α ) ( - italic_a ) = - italic_A italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT over~ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_α = - italic_A italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW

Also we can write

δΔEgAδz0=ΔFA=gC.𝛿Δsubscript𝐸𝑔𝐴𝛿subscript𝑧0Δ𝐹𝐴𝑔subscript𝐶-\dfrac{\delta\Delta E_{g}}{A\delta z_{0}}=\dfrac{\Delta F}{A}=g\mathcal{E}_{C}.- divide start_ARG italic_δ roman_Δ italic_E start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG italic_A italic_δ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = divide start_ARG roman_Δ italic_F end_ARG start_ARG italic_A end_ARG = italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT . (33)

When we add it to the isotropic force (FI=2gCsuperscript𝐹𝐼2𝑔subscript𝐶F^{I}=-2g\mathcal{E}_{C}italic_F start_POSTSUPERSCRIPT italic_I end_POSTSUPERSCRIPT = - 2 italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT), we obtain the Fermi force as follows [42, 43, 22, 44]

FI+ΔFA=2gC+gC.superscript𝐹𝐼Δ𝐹𝐴2𝑔subscript𝐶𝑔subscript𝐶\dfrac{F^{I}+\Delta F}{A}=-2g\mathcal{E}_{C}+g\mathcal{E}_{C}.divide start_ARG italic_F start_POSTSUPERSCRIPT italic_I end_POSTSUPERSCRIPT + roman_Δ italic_F end_ARG start_ARG italic_A end_ARG = - 2 italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT + italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT . (34)

At the end, we get

gC=FFA.𝑔subscript𝐶superscript𝐹𝐹𝐴-g\mathcal{E}_{C}=\dfrac{F^{F}}{A}.- italic_g caligraphic_E start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT = divide start_ARG italic_F start_POSTSUPERSCRIPT italic_F end_POSTSUPERSCRIPT end_ARG start_ARG italic_A end_ARG . (35)
Refer to caption
Figure 4: change in Force as a function of distance L for different area A
Refer to caption
Figure 5: change in Force as a function of area A for different distance L
Refer to caption
Figure 6: Fermi force and change in force per unit area as function of distance L

Concluding Remarks

By incorporating the Einstein-Langevin equation to account for stochastic fluctuations, we have derived first-order corrections to the semi-classical gravity framework. Our results demonstrate that these corrections are highly sensitive to both the distance between parallel plates and the order parameter of the weak field. This sensitivity underscores the intricate interplay between quantum fluctuations and gravitational influences, revealing that even subtle variations in these parameters can substantially alter the Casimir force. These findings not only enhance our understanding of the Casimir effect within the context of stochastic gravity but also highlight the broader implications of quantum field fluctuations on macroscopic gravitational phenomena. Further exploration of this interplay could potentially unlock new avenues in quantum gravity research and its applications in modern physics.

In conclusion, our investigation into the effects of weak field stochastic gravity on the Casimir force has yielded results which illustrated in the figures. In Figure 1, Equation 3 is utilized to plot the Casimir energy density as a function of the distance between the plates. Figure 2 illustrates the Casimir force per unit area, derived from Equation 26, which is inversely proportional to the fourth power of the distance between the plates. Furthermore, as indicated by Equation 25, the Casimir force is shown to be inversely proportional to the cube of the distance between the plates, as depicted in Figure 3. Figure 4 employs Equation 33 to plot the change in force for different areas A𝐴Aitalic_A, while Figure 5 shows the change in force as a function of area A𝐴Aitalic_A for varying distances. Finally, Figure 6 presents the change in force per unit area (light-green-sea dashed) and the Fermi force per unit area (violet-red dot-dashed) as functions of the distance L𝐿Litalic_L.

Acknowledgments

The authors would like to thank M. Amraji for his corporation and interest in this work.

References

  • [1] H.B.G. Casimir, On the attraction between two perfectly conducting plates, Proc. K. Ned. Akad. Wet.51,793 (1948)
  • [2] M. Bordag et al., Advances in Casimir Effect, Oxford University Press, Oxford (2009).
  • [3] G. Plunien, B. Muller and W. Greiner, Phys. Rep. 134,87 (1986).
  • [4] V.M. Mostepanenko and N.N. Trunov, The Casimir Effect and its Applications, Oxford University Press, Oxford, England (1997).
  • [5] S. Leseduarte and A. Romeo, Ann. Phys. 250, 165 (1978).
  • [6] E. Elizalde, S.D. Odintsov, A. Romeo, A.A. Bystenko and S. Zerbini, Zeta Regularization Techniques with Applications, World Scientific, Singapore (1994).
  • [7] E. Elizalde, Ten Physical Applications of Spectral Zeta Functions (springer, Berlin, 1995).
  • [8] S.W. Hawking, Commun. Math. Phys. 55, 133 (1997).
  • [9] K.A. Milton, L.L. DeRaad and J. Schwinger, Casimir self-stress on a perfectly conducting spherical shell, Ann. Phys. 115 (2), 388-403.
  • [10] D. Deutsch and P. Candelas, Boundary effects in quantum field theory, Phys. Rev. D 20, 3063 (1979).
  • [11] Belyad, M., Tanhayi, M.R. On Casimir Energy and Mutual Information in Non-relativistic Backgrounds. Int J Theor Phys 59, 1905–1916 (2020). https://doi.org/10.1007/s10773-020-04462-9
  • [12] M. Belyad and M. R. Tanhayi, Casimir Energy in Non-relativistic Backgrounds: Numerical Approach, International Journal of Mathematical Modelling and Computations, IJM2C-2305-1274, doi: 10.30495/ijm2c.2023.1985644.1274
  • [13] Tavakoli ,F., Pirmoradian, R. , Parsabod, I. Weak Gravitational Wave and Casimir Energy of a Scalar Field. Int J Theor Phys 55, 4100–4108 (2016). https://doi.org/10.1007/s10773-016-3038-3.
  • [14] H. Khosravi, M. Naseri, S. Rouhani and M. V. Takook, “Casimir Effect In Krein Space Quantization,” Phys. Lett. B 640, 48-51 (2006) doi:10.1016/j.physletb.2006.07.027 [arXiv:gr-qc/0604036 [gr-qc]].
  • [15] Pirmoradian .Reza , Tavakoli. Farhad, (2016), Vacuum energy through Krein quantization approach, Canadian Journal of Physics. 95. 119-124. 10.1139/cjp-2016-0013.
  • [16] Pejhan H, Tanhayi MR, Takook MV. Casimir effect for a scalar field via Krein quantization. Annals of Physics. 2014 Feb 1;341:195-204.
  • [17] B. Forghan, M. R. Tanhayi, Mod. Phys. Lett. A 30 , 1550126 (2015) Krein regularization of quantum chromodynamics.
  • [18] Ichinose S. Renormalization Group Approach to Casimir Effect and the Attractive and Repulsive Forces in Substance. arXiv preprint arXiv:1203.2708. 2012 Mar 13.
  • [19] Milton KA. The Casimir Effect: Physical Manifestations of Zero-Point Energy. World Scientific; 2001.
  • [20] Bordag M, Klimchitskaya GL, Mohideen U, Mostepanenko VM. Advances in the Casimir Effect. Oxford University Press; 2009.
  • [21] Decca RS, Lopez D, Fischbach E, Krause DE. Tests of New Physics from Precision Casimir Force Measurements. Phys Rev Lett. 2007;97(22):223602.
  • [22] Tanhayi, M.R., Pirmoradian, R. Weak Gravitational Field and Casimir Energy. Int J Theor Phys 55, 766–773 (2016). https://doi.org/10.1007/s10773-015-2714-z
  • [23] L.H. Ford, Spacetime in Semiclassical Gravity, arXiv:gr-qc/0504096v2 (2005)
  • [24] V.N. Lukash, A.A. Starobinsky. Isotropization of cosmological expansion due to particle creation effect. Zh. Eksp. Teor. Fiz., 66, 1515, 1974       (Sov. Phys. - JETP, 39, 742, 1974); L. Grishchuk, Ann. N.Y. Acad. Sci. 302, 439 (1976); B. L. Hu and L. Parker, Phys. Lett. A 63, 217(1977); B. L. Hu and L. Parker, Phys. Rev. D 17,933 (1978); J. B. Hartle, Phys. Rev. D 20, 1772 (1979).
  • [25] J. B. Hartle, Effective-Potential Approach to Graviton Production in the Early Universe, Phys. Rev. Lett. 39, 1373(1977); M. Fischetti, J. B. Hartle and B. L. Hu, Phys. Rev. D 20, 1757 (1979); J. B. Hartle and B. L. Hu, Phys. Rev. D 20, 1772 (1979).
  • [26] R.L. Jaffe, The Casimir Effect and the Quantum Vacuum,Phys. arXiv:hep-th/0503158, Rev. D72 (2005) 021301
  • [27] S. Weinberg, The cosmological constant problem, Rev. Mod. Phys. 61, 1 (1989)
  • [28] R. Martin and E. Verdaguer, On the semiclassical Einstein-Langevin equation, PHYS LETT B, 465(1-4), 1999,
  • [29] B.S. DeWitt, Phys. Rep. 10C, 295 (1975)
  • [30] N.D. Birrell and P.C.W. Davies, Quantum fields in curved space (Cambridge University Press, Cambridge, England, 1982)
  • [31] S.A. Fulling, Aspects of quantum field theory in curved spacetime (Cambridge University Press, Cambridge, England, 1989).
  • [32] R.M. Wald, Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics. The University of Chicago Press, Chicago (1994)
  • [33] A.A.Grib, S.G.Mamayev, and V.M.Mostepanenko, Vacuum Quantum Effects in Strong Fields, Friedmann Laboratory Publishing, St. Petersburg, 1994.
  • [34] B.L. Hu, Stochastic Gravity: A Primer with Applications, arxiv: gr-qc /0211090v1 (2002).
  • [35] B.L. Hu and E. Verdaguer, Recent Advances in Stochastic Gravity: Theory and Issues.Advances in the Interplay Between Quantum and Gravity Physics, Volume 60 of the series NATO Science Series pp 133-218 (2002).
  • [36] B.L. Hu, Stochastic Gravity, arXiv: gr-qc/9902064v1, Int.J.Theor.Phys.38:2987-3037 (1999).
  • [37] B.L. Hu, A Kinetic Theory Approach to Quantum Gravity, arXiv: gr-qc/0204069 ,Int.J.Theor.Phys.41:2091-2119(2002).
  • [38] K. Kirsten, Spectral Functions in Mathematics and Physics, Chapman and Hall/CRC, Boca Raton, FL (2002)
  • [39] Ph. Burda, Casimir Effect for a Massless Minimally Coupled Scalar Field Between Parallel Plates in de Sitter Spacetime. ISSN 0021-3640, JETP Letters, Vol. 93, No. 11, pp. 632-637. arXiv:1101.2624 [hep-th], (2011).
  • [40] L. Parker, D. Toms, Quantum Field Theory in Curved Spacetime: Quantized Fields and Gravity. (Cambridge University Press, Cambridge, 2009).
  • [41] L.S. Brown, G.J. Maclay, Vacuum stress between conducting plates: an image solution, Physical Review 184 (5), 1272 (1969).
  • [42] S.A. Fulling, K.A. Milton, P. Parashar, A. Romeo, K.V. Shajesh, J. Wagner, How does Casimir energy fall? Phys,Rev, D 76 (2), 025004. arxiv: hep-th/ 0702091 (2007).
  • [43] S.Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, New York,Wiley (1972).
  • [44] J. Schwinger, Particles, Sources, and Fields, Vol. 1, (Addison-Wesley, 1970).