Spiral Galaxies, no dark matter

Spiral Galaxies

Galactic rotation curves without dark matter.

A new paradigm for the universe.

https://msp.warwick.ac.uk/~cpr/paradigm

ISBN: 9781973129868

I keep dipping into this book.

Each time with new understanding. Each time with a new perspective.

It is a wonderful work, with compelling arguments.

Chapter 2, Sciama’s principle finishes with:

Sciama’s initiative, to base a dynamical theory on Mach’s principle as formulated in Sciama’s principle, has never been followed up and this approach to dynamics remains dormant. One of the aims of this book is to reawaken this approach.

One of the aims of blume is to help understanding of such theories.

In particular, help my own understanding with computer simulations.

This project will explore some of the mathematics of the book.

It will require astropy.

Solar system, galactic and local galaxy visualsation and simulation.

Simulation of gravitational waves generated by black hole mergers using the physics of the book.

Simulations of gravitational wave from a new galactic arrival.

2021 update

The Geometry of the Universe

Colin Rourke

https://www.worldscientific.com/worldscibooks/10.1142/12195

ISBN: 978-981-123-386-9 (hardcover)

A new book, with a bold claim. A stunning story, a revival of some old ideas and new brilliance too.

Spiral

Spiral galaxies? Why spirals?

The problem noted by observational scientist Vera Rubin was that stars in the outer reaches of spiral galaxies are moving far too fast.

This observation is what lead to the dark matter assumption: there must be some matter we can’t see that is dragging the stars along.

Colin Rourke, says nonesense, there must be a giant rotating mass at the centre, dragging things along.

Einstein’s general relativity withstands double pulsar’s Scrutiny

https://physics.aps.org/articles/v14/173

Strong-Field Gravity Tests with the Double Pulsar

13. Kramer et al.

Phys. Rev. X 11, 041050 (2021)

Published December 13, 2021

Binary pulsars provide strong tests of general relativity. The idea is that there are two, rapidly rotating, neutron stars in a close orbit around each other.

The Hulse-Taylor binary system has long been the best observed system and the strongest test.

A new system PSR J0737-3039A/B has now taken that crown as the best observed system. Lets call it Jumb0 737, from the J designate, or Jumbo for short.

The rotations of the Jumbo are incredibly stable, at the level of an atomic clock, allowing very precise measurements of the time delays of the pulses of energy as the system makes each rotation.

In addition, Jumbo is nearer than the Hulse-Taylor system and there is much less uncertainty in our estimate of its distance from earth.

It all gets complicated, since as although the system is in our own galaxy, the intervening spacetime modulates the signal we receive.

The new paper explains the precise measurements that can now be made and how variations on general relativity, can, in some cases, be tested.

In particular, the paper now mentions that rotation of the faster rotating body now has to be taken into account in the equations of state of the entire system. The other body is rotating around one hundred times slower so the effect can largely be ignored.

The aproach that is taken is to assume the Lense-Thiring effect.

This is discussed in some detail in chapter 3 of the book, the biggest blunder…. Particular attention should be paid to the discussion around page 60.

There are some quite technical arguments. The Lense Thiring effect is intuitively compelling. General relativity tells us that if space-time is a vacuum, then the Kerr metric is the only solution that fits Einstein’s equations.

This is problematic for the book’s arguments as the Kerr metric falls off as 1/(r**3) whereas the book argues it should fall off as 1/r.

Rourke makes the assumption that linear motion has no inertial effect and notes that you can change angular momentum by adding a linear motion, whereas angular velocity cannot be changed in this way.

Is space time a vacuum?

The impasse is resolved if we assume that space is not in fact a vacuum.

For starters, there is an awful lot of microwave radiation buzzing around: the cosmic microwave background.

The nature of this background depends critically on our assumptions about the geometry of the universe: big bang, or static universe?

What would a static universe look like given the effects of special relativity, simplifying a universe to the set of galaxies above a certain size?

This model will instead show how by assuming this 1/r relation for the effect of a masses rotation on its surroundings, produces galactic rotation curves very much in line with observations.

What about Jumbo?

So the cool news is that we now get estimates of the angular velocity of each body and it’s moment of inertia.

The latter adds some uncertainty to model fit as there is uncertainty of the exact distribution of the matter within each neutron star, which is important to know, as the model being used assumed it is angular momentum that drives the rotational frame dragging.

As noted above, Rourke argues that it is angular velocity, rather than angular momentum that matters in the calculation, in short the matter distribution within the black hole is not required for his model.

Further, that the Lense-Thiring effect drops off as 1/r, not the 1/r**3 that is presumably being used by the new paper. Unless the orbits are highly oblique, this difference is not going to be detectable with a single system.

This should mean that a fit using Rourke’s model reduces the uncertainty in all the parameters that the fitting process estimates, since the uncertainty in the matter distribution no longer comes into play, just the uncertainty in angular velocity.

Another project would be to fit the model to the latest Hulse-Taylor data and see what changes.

I am also curious how Rourke’s model affects the long term evolution of these systems. My hunch (actually I think I read something along these lines in the book) is that the feedback from the rotation keeps these binary systems stable and that it is highly improbable that they will in-spiral and coalesce.

Which raises the interesting question of what is the source of the waves that our gravitational wave detectors are seeing?

Which reminds me, I need to work on the grb module.

What about the spiral module?

The idea is to create a visualisation of the formation of spiral galaxies with a omega m / r model.

It would be good to also be able to model binary systems while we are at it.

class gotu.spiral.Cosmo(cosmo=None)

Mimic an astropy Cosmology object

It’s a Sciama-DeSitter Universe.

Static, but doesn’t look static.

Old and big in a state of equilibrium.

The O*0 attributes give the current estimates of the share of the critical mass (also at z=0) in various buckets.

Ogamma for photons (primarial CMB) Om for non-realtivistic matter Ob for baryonic Odm for dark matter. Ode for dark energy.

Dark matter is not part of this cosmology.

Instead assume this turns up as black holes in the centre of galaxies.

check_critical_density()

Critical density formulae

The critical density only depends on the Hubble distance.

It is the mass per unit volume such that the total mass has Schwartzchild radius equal to the Hubble distance.

As such, it is determined by the Hubble constant.

is_flat()

Not sure what the answer is to this.

Perhaps, “asymptotically flat”?

class gotu.spiral.LLPegasi

Beautiful spiral with an 810 year period.

Dust around a binary system showing a clear spiral structure.

The gap between bands and their rate of expansion are consistent with the pair’s 810 year orbit based on their angular separation.

Distance is 1300 * u.pc

temperature 1800K

class gotu.spiral.SkyMap(gals=None, n=1566, fudge=None)

Yet another table viewer?

parsing csv and figuring out what it means.

or just give me lists of fields i can pick form for any purpose?

Or yet another universe viewer.

And so we descend into the world of coordinate systems.

Something astronomers understand very well.

See gotu.wits for more on this subject.

Maybe SkyMap allows systems to evolve over time, according to the paradigm.

async cmbsim()

Calculate some CMB related numbers

This is really just a test of planck_radiance_law_wavelength, checking that the energy gives the same Ogamm0 as is implied by the temperature.

create_sample(gals=None, n=None, t0=0)

Create a new sample of galaxies

This also initialises some parameters, based on the sample.

The goal is to set the stellar masses to match the baryonic matter as given by the cosmology Ob0. This is then scaled so that the sum of the central masses equals the mass of the universe.

gals: if given, a list of galaxies to include.

find_distant(results)

Find objects we no longer care about

async local_mode_sim()

Simulation of the gravity waves from the local cluster

sample_mass()

Return the sum of the masses of the sample

In solar masses.

set_mcent()

Set up central mass of all spheres based on stellar mass

tofu(u)

Work in natural units, c=G=Hubble Distance=1

Time in Hubble times.

waves()

What waves would be generated by the catalog?

class gotu.spiral.SpanishDancer

Spanish Dancer

class gotu.spiral.Spiral

Model a spiral galaxy

Or any rotating mass?

Want to convert this to use astropy units.

accretion()

Accretion rate

amplitude(R=None)

amplitude of the gravitational wave wave

Returns the amplitude of the wave if there was just one such object in the sphere radius R. Scale by the actual density to get useful numbers.

Suppose rho is the number of such objects per unit volume, the the full wave will be:

a = 4 * pi rho Sum M r^2 dr / r for 0 < r < R

Which simplifies to:

a = 4 * pi * rho * M Sum r dr for 0 < r < R

or:

a = 2 * pi * rho * S * R ** 2 [1]

Now:

N = rho * 4/3 * pi * R**3

So if N==1, rho = 4 * pi / (3* R**3), substituting in [1]

a = 3 S / 2 R

Where S be the Schwartzchild radius of Mcent.

R is the radius out to which bodies can contribute to the sum.

To get the full wave, scale by the number of objects in the region of radius R.

bondi()

The Bondi Sphere

Equate the root mean square velocity of hydrogen atoms in the medium, sqrt(3kT/m_h) with the escape velocity sqrt(2GM/B).

critical_radius()

Radius of Mcent based on critical density

Convert Mcent to the radius of the sphere it would occupy based on critical density.

density()

Return density of Mcent if it is a black hole

eddington(z=None)

Eddington Sphere

Easiest to calculate if we know the redshift.

galaxy()

Set parameters for a galaxy

rmin_check()

The length of the roots of the spirals

This can be used to set the B value.

Assume that the spiral roots end at radius r0

And assume the roots are moving with the inertial frame at that radius.

The rate of precession will match that of the inertial frame at that radius.

schwartzchild()

Schwartzchild radius in lightyears

spheres3()

Bondi, Eddington and Schwartzchild

We have mass in lightyears, the Schwartzchild radius.

To get the mass in kilograms, S = 2GM/c^2, M = Sc^2/2G

tb()

blue shift period

tofu(u)

Work in natural units, c=G=Hubble Distance=1

Time in Hubble times.

v(r)

Velocity at radius r

A = 0.0005 K = Mcent CC = -10

??

vinert(r, v)

Inertial part of the velocity

Part of velocity relative to inertial frame.

Notes

K is central mass. A = 0.0005

z_calc()

1.27e7 theta / M

M mass in suns

class gotu.spiral.Sun

gotu.spiral.cpr()

Started as Mathematica code from the new paradigm.

adapted to python over time.

See spiral class for more information over time.

gotu.spiral.from_heasarc(kwargs)

Create a Spiral from a Heasarc record

gotu.spiral.hubble_tension(cmb=<Quantity 67.66 km / (Mpc s)>, near=None)

Explore hubbe tension

gotu.spiral.mass_length()

Equivalence of mass and length based on Schwartzchild radius

This is for use with astropy constants.

for example: >>> u.solMass.compose(mass_length())

gotu.spiral.near_galaxies()

parse galaxy.txt from

https://heasarc.gsfc.nasa.gov/w3browse/all/neargalcat.html

gotu.spiral.planck_radiance_law_wavelength(wavelength, T=None)

Energy emitted in wavelength at temperature T

wavelength: a float, wavelength in meters T: temperature, a quantity with units K

returns astropy Quantity object with units

gotu.spiral.schwartzchild(m)

Mass in kg to schwartzchild radius in lightyears