MATLAB/Divergent series investigations/Casimir effect

Introduction edit

While the template labels this as Original research, it is amateur original research. The greatest theoretical physicists on this planet have thought deeply about the Casimir effect for many decades. Our goal is to make it understandable to more people. It is Research because the outcome is unknown. The method proposed here for calculating the Casimir effect may prove entirely wrong.

But, if this method is correct, you can be the first to do it. All you need is a simple computer program that can add perhaps 100,000 terms of a w:divergent series. The code follows that already posted on Wikipedia's MATLAB/Divergent series investigations. In that code, a trick was employed to convert the obviously infinite (divergent) series,

1 + 2 + 3 + 4 + ...

into a series that converges to a number very close to -1/12:

1.0 + 2.0 + 2.9 + 3.8 + bigger numbers +...+ smaller numbers +...+ negative numbers +...+ very small numbers

Only this time, we need the triple sum to include all the waves that can exist inside a rectangular box. Each wave has positive and well-known "ground state" energy of  . Somehow these energies need to add up to a negative number (...because an infinite number of waves are present ! ? ).

From the wikipedia article edit

The Wikipedia article w:Casimir effect shows how the vacuum energy of a box is infinite, and proportional to the following divergent series, S:

Algebraic steps

Following Wikipedia's Derivation of Casimir effect assuming zeta-regularization we consider a box of area A and thickness a<<A. The energy density per unit area, <E>/A, is given by the following expression in the limit  :

 

So as to not get our "knickers in a twist" with w:zeta function regularization and the w:Riemann zeta function , we set s=0 exactly

 

We may further refine the integral with the change of variable to   where,

 
  =  
 

Hence we seek

 

As was done to 1+2+3+4+..., we "tame" this series with a damping factor, and later take   and   along a path that yields a finite answer for the real part of the series.

matlab code edit

see matlab code for these calculations (casimirSumIntegral.m)


Tame the infinity edit

In a [MATLAB/Divergent series investigations|previous exercise] the w:divergent series 1+2+3+... was "tamed" by introducing a damping term into an integral approximation to the series. Then, when the tam

 

Setting  , we have,

 

which is easily converted into a single integral by a sort of "spherical coordinates" with r2 = x2 + n2, dxdn=rdrdθ, and x=rcosθ. After integrating over θ from 0 to π/2, we have:

 

The real part vanishes if   is imaginary. Denoting the real and imaginary parts of  , we have,for example

 

On the physical nature of these equations edit

Even though we don't have an answer, our equations are making sense, especially if we focus attention on the real parts of   and  :

 


 

Even though   for all  , the expression for   contains large terms. The expression for   also contains large terms. Moreover, in the limit  , the contributions to both   and   come from large values of   and  . Large values of   and   correspond to high frequency (short wavelength). Yet, the Casimir effect is dominated by the lowest frequencies (longest wavelengths) that the box can contain , as can be seen by the following manipulation of these expressions:

Highlighting the importance of low-order modes edit

  Since  , we have  

     

The two terms in this expression are identical, except that one contains an integral over n, while the other is a sum. We can emphasize this by defining an operator that acts on any function f:

 

In our case,

 

If it can be established that this operator tends to suppress high order modes, the following scenario explains the Casimir effect: Empty space is occupied by an infinite number of ground state energy levels, and the shorter wavelengths contain most of the energy. A small metal box also contains an infinite amount of energy, with the shorter wavelengths containing most of the energy. The energy levels in the box are quantized (discrete), but the spectra of the box and empty space differ only at the lower frequencies. Hence, while the higher energy modes contain most of the energy, the difference in energy density between as small box and empty space is dominated by the lowest energy modes.

Unfortunately, f(x) is large and unpredictable. I need to reverse the operations?