Archive for March, 2017

Black holes are a theoretical concept invented by Oppenheimer-Snyder who showed that pure Einstein equations lead to a collapse of a star to a pure singularity (Oppenheimer_1939).  In 2014 Laura Mersini-Houghton shows that with additional considerations a black hole cannot form because the collapsed star must have a finite radius (1408.1449).  Her papers are published in Physics Letters B so this is considered seriously.  I think this makes sense.  The idea of black holes is based only on a very simple model whose conclusions were just extrapolated wildly.  Just like the Big Bang singularity these probably are figments of our imagination.  Singularities in spacetime have never been observed.  It’s fairly crazy to extrapolate about such things.  Of course experimenters have ‘found’ black holes.  Experimenters make their living finding things that theorists tell them to find.

Read Full Post »

This is a very simple simulation that is not meant to reproduce an exact fit to CMB anisotropy but to get some feel for whether CMB anisotropy occurs only with gravitational redshift as the only source of anisotropy which it does not try to do exactly.  A static Einstein universe is an excellent model now that we know that the original reasons for abandoning it were not quite right (Eddington’s instability result).  The six parameter Lambda-CDM model fits the anisotropy well which is the model to overthrow.  So is it possible to do this?  I believe that a static Einstein universe is an infinitely superior model (as it gets rid of the cosmological constant problem since quantum field theory on static Einstein universe is actually close to the measured value of around 10^{-52} cm and the quantum gravity problem is not as difficult here since for example the Wheeler-deWitt equation can be solved explicitly).  So here is a little simulation that shows that there is some interesting anisotropy in the simplest type of model for gravitational redshift effect.  Distribute a uniform set of points over a round three-sphere and assign them masses by an exponential distribution.  Choose a point on the three-sphere as a reference pont and draw a two-sphere around it and histogram the geodesics meeting the 2-sphere.  Assume the background radiation has 2.7 K uniform background so that the anisotropy is due only to the gravitational redshift of photons that scatter on the fixed objects (scattered uniformly in the three-sphere).  Smoothen the three-sphere histogram and extract the spherical harmonics coefficients for them and plot the power spectrum.  We ask:  do we see an anisotropy that gives any sense of the anisotropy.  The answer to this question is yes (in a relatively quick implementation).  This model can be refined to fit the actual CMB anisotropy I am confident.  If I don’t succeed, someone else should do this because the Big Bang cosmology is insanely ambitious in claiming anything about the ‘origins of the universe’.  A static Einstein universe would be infinitely better science with what we actually can know about the universe without any ability to do any experiments.  Of course I am not doing any experiments either but I am doing real work:  find a tight fit to the CMB anisotropy in a SIMPLE model.  This cannot be that hard.

The code is in octave:  The first example of anisotropy for the model is this little graph:


pkg load statistics
pkg load gsl
pkg load image

X = mvnrnd( [0,0,0,0],diag([1,1,1,1]),n);
for i=1:n,

# stereographic projection + normalize

function z=pwrrnd(k,n)

masses = exprnd(1,n)

function d=direction(x1,x2)

function d=distance(x1,x2)

function f=gforce( x1,x2, m1,m2)

for a=1:n,
force_on_x = 0

x1 = X(a,:)
m1 = masses(a)
for b=1:n,
x2 = X(b,:)
if d>0 and d<pi,
force_on_x = force_on_x + gforce(x1,x2,m1,m2)
x = x+0.5*force_on_x*timestep^2/m1
x = x/norm(x,2)

z = zeros(n,3);
for i=1:n,
zp=[ y1/(1-y4), y2/(1-y4),y3/(1-y4) ];
z(i,:) = zp/norm(zp,2);

function t=fitangle(x)
while t>2*pi,
while t<0,

function b=bin_s2point( z )
theta = atan( z(2)/z(1));
phi = atan( sqrt(z(1)^2+z(2)^2)/z(3));
phi = fitangle(phi);
thetan = floor(theta*100/(2*pi));
phin = floor(phi*100/(2*pi));
# create a spherical grid in longitudes/latitudes
g = zeros(q0,q0);

refpoint = [1,0,0,0]
for i = 1:n,
cds = bin_s2point(z(i,:));
sdir = sign(direction(refpoint,X(i,:));
dist = distance(refpoint,X(i,:));
dg = zeros(q0,q0);
dg(cds(1)+1,cds(2)+1) = 1;
dg = imsmooth( dg, “Gaussian”, ceil(m1),ceil(m1));
g = g+dg;

#g = (g-min(min(g)));


# value at a spherical harmonic

function v=sphCoeff( l,m, g)
v = 0;
for j=1:q0,
for k=1:q0,
w1 = cos(2*pi*j/100);
w2 = exp(1i*2*pi*k/100);
fnval = gsl_sf_legendre_sphPlm( l, m, w1 )*w2;
v = v + fnval*g(j,k);

function v=sumSphCoeff( l, g)
for m=0:l
v = v+sphCoeff(l,m,g);

valsByL = zeros(1,100);

for k=1:100,

plot( 1:100, valsByL)



Read Full Post »

I took Differential Geometry at Princeton in my sophomore year in 1992 and have become incredibly rusty so it is time to do some computations in General Relativity to get back in shape.  An elementary problem: find the light-cone geodesics for the metric

g = dt^2 - e^{-2Ht} dx^2

which satisfy g_{ab} \dot{x}^a \dot{x}^b=0.  The Christoffel symbols are computed from

\Gamma_{abc} = \frac{1}{2}( g_{ca,b} + g_{cb,a} - g_{ab,c})

There are too many of these to compute exhaustively so we note first that

g_{ab,c} = 0 unless a=b=2, c=1

Therefore the only nonzero \Gamma_{abc} are those with two 2 and one 1 in the indices.  For example,

\Gamma_{212} = \frac{1}{2} ( g_{21,2} + g_{22,1} - g_{12,2}) = \frac{1}{2} g_{22,1} = -H e^{-2Ht}


$\Gamma_{122}=\Gamma_{212}=\Gamma_{221} = -H e^{-2Ht}$

Then we have to raise the first index which has the solution

\Gamma^1_{ab} = g^{11} \Gamma_{1ab} + g^{21} \Gamma_{2ab} = \Gamma_{1ab} = -He^{-2Ht}

\Gamma^2_{ab} = g^{12}\Gamma_{2ab} + g^{22}\Gamma_{2ab} = e^{+2Ht} \Gamma_{2ab} = - e^{+2Ht} H e^{-2Ht} = -H

Now the geodesic equation is

\frac{ d^2{x}^c}{ds^2} + \Gamma^c_{ab} \frac{dx^a}{ds} \frac{dx^b}{ds} = 0

For coordinates t,x we have

\ddot{t} + (-He^{-2Ht}\dot{x}^2 = 0

\ddot{x} - 2H\dot{t}\dot{x} = 0

I don’t know how to solve this system of equations in general.  But the light cone condition g_{ab}\dot{x}^a\dot{x}^b =0 in this case is

\dot{t}^2 = e^{2Ht} \dot{x}^2

which we plug into the first equation \ddot{t} - H\dot{t}^2=0.  This is \dot{z} = Hz^2 solved by integrating \dot{z}{z^2} = H, so z=-\frac{1}{Hs} and t = -\frac{1}{H} \log(s).

The second equation now simplifies when we plug in \dot{t} = -\frac{1}{Hs}.

\ddot{x} +2/s \dot{x} = 0

Again \dot{z} = -(2/s)z and \log(z) = -\log(s^2) +C and x = \int^{s} \exp(-r^2) dr.

Read Full Post »

This is explained in Burbridge-Hoyle 1998 (Burbidge_1998_ApJ_509_L1).  So issues regarding Helium abundances do not require Big Bang nucleosynthesis.  This is useful to reconsider the Einstein static model as the cosmological model.  Not only quantum field theory but the constructive Euclidean quantum field theory can be done in a static universe without problems (Jaffee-Ritter etc.)  One should not need expansion, dark energy or dark matter with a fixed geometry of spacetime.  More useful would be a cleaner theory of quantization in a fixed geometry where the cosmological constant problem disappears immediately.

Read Full Post »

The two original explanations of the observed distance-redshift relation are:

  • Doppler effect (universe is expanding)
  • Tired light theory (light traveling from distance sources getting absorbed and re-emitted by dust)

My new explanation — several years old — the universe is static with fixed curvature.  My preferred model is a membrane in a four-sphere but the explanation is equally good in the original Einstein static three-sphere.  This is that spherical waves do not have an inverse relation between frequency and wavelength.  It is easiest to see this by considering a circle of radius R versus a three sphere of the same radius.  The frequencies of standing waves on the circle have eigenvalues of Laplacian $C n^2$ while on a three-sphere they are Cn(n+2). For light traveling on a sphere of a large distance, the resulting wavelength will seem to grow linearly with distance while there has been no actual physical feature of the universe that has changed.  Thus my explanation of the redshift is that it is a MATHEMATICAL ARTIFACT of expecting an inverse wavelength-redshift relationship while in fact it is just an effect of curvature.

One could then bring up the evidence for acceleration on measurements of supernovae.  Well this is actually a real redshift but this is actually graviational redshift.  Heavier objects have a higher redshift due to time dilation according to Einstein, and in fact there is work addressing this issue: gravitationalredshift.  So once you correct for this there is no acceleration of any expansion.

DARK ENERGY does not exist.  It is simply the curvature of a static spherical eternal universe. And the Big Bang singularity does not exist either.  By the way, the cosmic background radiation was predicted by steady state theorists as the temperature more accurately before Gamow which is well known:

Read Full Post »

The classical velocity orbit formula is

v = \sqrt{ GM/r}

where v is the orbital velocity and M is the central mass.  It is NOT v = Cr^{-3/2}.  It is v= C r^{-1/2}.  So the issue is what is a reasonable model of mass distribution of a galaxy in terms of radius r.  The most reasonable hypothesis is that it is M(r) = C r^2 assuming it’s two dimensional.  Then the expected velocity rotation curve is

v= C r^{1/2}

So all the graphs that people present showing anomalous rotational velocities are severely misleading because the observed galaxy rotations are actually extremely close to the above relation.  A discrepancy between the observed velocities and the above relation can be easily explained by GRAVITATIONAL REDSHIFT.  Gravitational redshift formula is approximately (see Gravitational redshift)

z_{approx} = \frac{1}{2} GM/{c^2 r}

So the approximate velocity due to gravitational redshift is just this multiplied by the speed of light.  I will illustrate by a simple procedure in R how to get a v = C\sqrt{r} from actual data so that the relation has statistical significance with R^2 > 0.8-0.9 on actual data which can be found here:Galaxy rotation dataGalaxy rotation dataGalaxy rotation data.  So take for example the data for DDO 64

r v verr
1 1.5 6.3 4.9
2 4.5 7.6 3.0
3 7.5 9.9 2.6
4 10.5 13.8 1.4
5 13.5 9.1 2.5
6 16.5 16.6 3.2
7 19.5 14.2 2.5
8 25.5 13.3 1.2
9 28.5 22.4 2.4
10 31.5 25.2 5.0
11 34.5 32.0 4.0
12 40.5 30.8 1.7
13 43.5 41.6 2.9
14 46.5 44.6 2.3
15 49.5 36.7 4.0
16 52.5 51.2 3.9
17 55.5 44.9 1.8
18 58.5 53.6 13.0
19 61.5 48.1 9.8
20 64.5 59.8 9.0

The columns are radius, velocity, velocity error.  We conside a power law relationship so in R summarize linear regression of log-velocity versus log-radius:

> summary(lm(log(v)~log(r)))

lm(formula = log(v) ~ log(r))

     Min       1Q   Median       3Q      Max 
-0.59982 -0.09126  0.01336  0.19051  0.54301 

            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  1.02706    0.21837   4.703 0.000177 ***
log(r)       0.66710    0.06551  10.183 6.75e-09 ***
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.2843 on 18 degrees of freedom
Multiple R-squared:  0.8521,    Adjusted R-squared:  0.8439 
F-statistic: 103.7 on 1 and 18 DF,  p-value: 6.752e-09

So we see that the relationship that fits (with very significant p-value) is a power law v = C r^{0.667}. This is a bit higher than the simple classical model suggests, i.e. v= C r^{0.5}. So let’s assume that this deviation is due to gravitational redshift, which by the same model will be based on the mass-by-radius formula for a two-disc, M(r) = Cr^2 where the constant C is different in varoous occurrences above since we are interested in the power of r in this exercise. Well, a little experimentation gives us:

> summary(lm(log(v-0.0059*r^2)~log(r)))

lm(formula = log(v - 0.0059 * r^2) ~ log(r))

     Min       1Q   Median       3Q      Max 
-0.64957 -0.08524 -0.00853  0.19738  0.35968 

            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  1.27580    0.20315   6.280 6.38e-06 ***
log(r)       0.50058    0.06095   8.214 1.68e-07 ***
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.2645 on 18 degrees of freedom
Multiple R-squared:  0.7894,    Adjusted R-squared:  0.7777 
F-statistic: 67.46 on 1 and 18 DF,  p-value: 1.681e-07

So if we consider the action of gravitational redshift to lower the velocities we recover approximately v = C r^{1/2} quite easily. It is not hard to check that this is true about the other examples as well. The significance of this exercise is that in fact there is no missing mass at all if we make the assumption that the mass is approxmately distributed as M(r)= Cr^2 which can be obtained from the luminous matter without the need for any mysterious dark matter. We do not need to modify gravity laws (which is what many people had done to explain the anomalous rotation curves). We don’t need to introduce dark matter either. The problem is that the ‘expected Keplerian velocities’ are a red herring since they assume that objects distributed inside the galaxy can treat the center with the same mass. This assumption is silly. A point at a given radius from the center should experience the center of mass as a function of radius obviously. So this tells us that standard Newtonian gravity with Einstein’s gravitational redshift as correction WITHOUT dark matter explains the galaxy rotations without any problems. There is no evidence of dark matter from galaxy rotation velocities.

Another example to ensure that we are not making an unjustified extrapolation.  Data for UGC 1281:


Without any gravitational redshift adjustment, we get a law v \sim C r^{0.7416}:

> summary(lm(log(v-0.00*r^2)~log(r)))

lm(formula = log(v - 0 * r^2) ~ log(r))

     Min       1Q   Median       3Q      Max 
-0.79011 -0.22565  0.00927  0.20012  0.73392 

            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   0.2390     0.3605   0.663    0.514    
log(r)        0.7416     0.1040   7.129 2.28e-07 ***
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.3665 on 24 degrees of freedom
Multiple R-squared:  0.6793,    Adjusted R-squared:  0.6659 
F-statistic: 50.83 on 1 and 24 DF,  p-value: 2.279e-07

With a loss of goodness-of-fit but without losing significance of the power law fit, we can obtain v_{adj} \sim C r^{1/2}:

> summary(lm(log(v-0.0035*r^2)~log(r)))

lm(formula = log(v - 0.0035 * r^2) ~ log(r))

     Min       1Q   Median       3Q      Max 
-1.01665 -0.23224  0.02452  0.20987  0.70148 

            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   0.7161     0.4291   1.669 0.108115    
log(r)        0.5051     0.1238   4.080 0.000431 ***
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.4363 on 24 degrees of freedom
Multiple R-squared:  0.4095,    Adjusted R-squared:  0.3849 
F-statistic: 16.65 on 1 and 24 DF,  p-value: 0.0004306

We are treating gravitational redshift with mass distribution M(r) = Cr^3 for simplicity. Our central point here is that galaxy rotation curves in the classical model with mass distribution M(r)=Cr^2 should follow v \sim Cr^{1/2} rather than v \sim Cr^{-1/2} so the anomaly in these rotation curves is small perturbation in the power law from expected \beta=1/2. Gravitational redshift can account for the discrepancy. There should be no need to modify gravitational laws and no need to hypothesize dark matter to explain the rotation curves.


Read Full Post »

The problem of quantum gravity is impossibly hard because it is too general.  Quantum gravity problem is not trivial but not nearly as hard if we consider a static spacetime such as Einstein static model \mathbf{R}\times S^3(a).  Quantum field theory on this specific spacetime is essentially a solved problem.  In fact, even Euclidean quantum field theory following Osterwalder-Schrader quantization and Wightman axioms has been worked out in this case ( jaffee-ritter-qft1jaffee-ritter-qft2 osterwalder-schrader-2osterwalder-schrader-1 ).  In presentations of quantum field theory on curved spacetimes the Einstein static spacetime appears often as a relatively simple example (see Fewster lecturenote-3908 (1)).

The trouble with quantum gravity is that it is only an inaccessibly hard problem because it is the wrong question on arbitrary Lorentzian spacetimes.  The problem is a bad cosmology: there is no expansion or any early universe.  The problem of quantum gravity is that it is attempting to solve an irrelevant and overly general problem of an expanding spacetime while in fact the problem can be solved when the spacetime is static and globally hyperbolic with a compact Cauchy surface.  It is in this setting where quickly we discover that there is no dark energy necessary in the sense that vacuum energy (which is well-defined for static spacetimes not suffering ambiguities of more general Lorentzian situations) is approximately equal to measured cosmological constant.  The dark energy is just the curvature of a static universe.  If one proceeds by a Wheeler-deWitt equation for quantum gravity, then this can also be solved explicitly via Huen functions etc.  Thus specialization of the quantum gravity problem to a specific static spacetime resolves the most pressing issues of quantum gravity.  It seems quite dangerous for thousands of physicists getting lost in speculative scenarios which will most likely have the same fate as the thousands of papers on solutions of Einstein equations without physical relevance.  The quantum gravity problem is a GLOBAL problem in the sense that it must necessarily be affected by the geometry of the entire cosmos.  It will be very unlikely to have a local solution that has merit.  The hardness of the problem seems to be due to trying to force local considerations lead to gravity.

Read Full Post »

Older Posts »