The primary argument levelled against Milgrom's MOND is that it has no theoretical support, even though considerable effort has been expended in attempting to provide it. Against that criticism, MOND irrefutably enjoys an expanding portfolio of success and so is almost certainly tapping into something fundamental. But what? Over roughly the same period that MOND has been a topic of controversy, Baryshev, Sylos Labini and others have been claiming, with equal controversy in earlier years, that, on medium scales at least, material in the universe is distributed in a quasi-fractal $D\approx 2$ fashion. There is a link: if the idea of a quasi-fractal $D\approx 2$ universe on medium scales is taken seriously, then there is an associated characteristic mass surface density, $\Sigma_F$ say, and an associated characteristic acceleration scale $a_F = 4\pi G \,\Sigma_F$. The whole success of MOND is predicated upon the idea of a critical acceleration scale, $a_0$. It is an obvious step to make the association $a_0 \sim a_F$ and then to consider the MOND critical acceleration boundary simply as a marker for a characteristic mass surface density boundary separating 'galaxy' from an environment characterized by $\Sigma_F$. This provides a route to the synthesis of conservative MOND from first principles. The radial acceleration relation (RAR) for conservative MOND when applied to the SPARC sample is the unity line. There is no mass discrepancy.

The claim that the large scale structure of the Universe is hierarchical has a very long history going back at least to Charlier's papers of the early 20th century. In recent years, the debate has centered largely on the works of Sylos Labini, Joyce, Pietronero and others, who have made the quantative claim that the large scale structure of the Universe is quasi-fractal with fractal dimension $D\approx 2$. There is now a concensus that this is the case on medium scales, with the main debate revolving around what happens on the scales of the largest available modern surveys. This paper, which is a realization of a worldview which is deeply rooted in the ideas of Leibniz and Mach shows, as a very special case of a general formalism, that such a fractal $D\approx 2$ world is necessarily a world of dynamical equilibrium. We use the cosmology to write down a simple galaxy model as an arbitrary spherically symmetric bounded perturbation of the $D\approx 2$ equilibrium environment. A necessary mass continuity condition on the perturbation boundary between galactic interior and the exterior environment immediately defines that boundary as a critical acceleration boundary with $a_0 \sim 4\pi G \Sigma_F$ where $a_0$ is the critical acceleration parameter and $\Sigma_F$ is the characteristic mass surface density of the $D \approx 2$ fractal environment. The MOND acceleration condition $V_0^2/R_0 = a_0$ for circular orbits on the critical boundary follows in a straightforward manner.

The SPARC sample consists of 175 nearby galaxies with modern surface photometry at $3.6\,\mu m$ and high quality rotation curves. The sample has been constructed to span very wide ranges in surface brightness, luminosity, rotation velocity and Hubble type, thereby forming a representative sample on galaxies in the nearby Universe. To date, the SPARC sample is the largest collection of galaxies with both high-quality rotation curves and NIR surface photometry. The neo-MOND model used here to analyse the SPARC sample recognizably conforms to the general pattern of the classical MOND algorithm, with the primary difference that, whereas the classical MOND model is purely phenomonological, the neo-MOND model is a special case of a general theory motivated by the ideas of Leibniz and Mach (not discussed here). The consequent main results can be broadly summarized as follows: (1) neo-MOND provides the basis for the derivation of a complete theory of the baryonic Tully-Fisher relation; (2) the details of the derivation provide a means of setting absolute distance scales for disc galaxies independently of standard candles and the photometric method; (3) subsequent determinations of dynamical mass (computed directly from neo-MOND fits to SPARC rotation curves) track photometric mass (estimated from SPARC surface photometry) across the whole SPARC sample in a statistically perfect way. To summarize, if the input to neo-MOND is disk dynamical data, then the output is whole-disk mass data together with absolute distance scales.

We show how the hypothesis that galaxy discs conform to self-similar dynamics leads to an objective identification, in the data, of an annular region of the optical disc which is such that, corresponding to the classical Tully-Fisher scaling law defined on the exterior annular boundary, there is a similar scaling law defined on the interior annular boundary. This result is confirmed at the level of statistical certainty over several large ORC samples. Furthermore, the same analysis provides insight into the uncertainties associated with the best way of defining Vmax, the rotation velocity used for the Tully-Fisher scaling law and Rmax, the galaxy radius at which Vmax is measured. Finally, as a direct consequence, we are led to a refined Tully-Fisher law which is largely insensitive to the means by which a galaxy's rotation velocity is defined.

Over the past two decades or so, it has become increasingly appparent that, out to quite large distances, galaxies are distributed in a quasi-fractal fashion with fractal dimension $D \approx 2$. Whether or not this behaviour continues onto indefinitely large scales is a matter of live debate and is a question that can only be settled if, at some scale, there is an unambiguous transition to homogeneity. This point has not yet been reached and may never be. This paper has been written on the basis of the tentative hypothesis that quasi-fractal $D\approx 2$ behaviour is a persistent and fundamental feature of galaxy distribution on all scales and addresses the question of the origins of this putative fractality. Given this tentative hypothesis then, except for the device of putting the fractal behaviour into the initial conditions - which is to by-pass the question, there is no obvious explanation within the framework of conventional cosmology. So we adopt the position that fractality is per-se a signature that different thinking is required. We find that a beautiful solution flows deductively from elementary observations about the world in which we live when we take seriously a view of space and time that can be traced, via Mach, Berkeley and Leibnitz to Aristotle. In summary, we find that a globally inertial space and time can be irreducibly associated with a fractal, D=2, distribution of material if we are prepared to review our understanding of what is meant by the notion of metric on astrophysical scales.

We test the periodicity of quasar redshifts in the 2dF and SDSS surveys. In the overall surveys redshift peaks are already apparent in the brighter quasars. But by analyzing sample areas in detail it is shown that the redshifts fit very closely the long standing Karlssson formula and strongly suggest the existence of preferred values in the distribution of quasar redshifts. We introduce a powerful new test for groups of quasars of differing redshifts which not only demonstrates the periodicity of the redshifts, but also their physical association with a parent galaxy. Further such analyses of the large area surveys should produce more information on the properties of the periodicity.

In previous work we have shown how a worldview that has its origins in the ideas of Aristotle, Leibniz and Mach leads to a quasi-classical (that is, one-clock) metric theory of gravitation (astro-ph/0107397) which, for example, when applied to model low surface brightness spirals (astro-ph/0306228), produces results that have, hitherto, only been matched by Milgrom's MOND algorithm. In this paper we show how the natural generalization of this worldview into a properly relativistic two-clock theory, applied to model a spherically symmetric gravitational source, produces results that cannot be distinguished from the canonical picture for all the standard local tests and which, when interpreted as a radiation model, produces no dipole radiation. Furthermore, although black-holes within this picture have an event horizon at the usual Schwarzschild radius, they do not have an essential singularity at the origin - the solutions are perfectly regular there.

In previous work (astro-ph/0306228) we have established that a particular quasi-classical theoretical viewpoint about the nature of gravitation in galaxy discs provides extremely high concordance (comparable with MOND) with the observed dynamics in optical discs. In this paper, we consider the possibility that the optical disc and the flat radio disc represent distinct dynamical regions which are inherently described by distinct phases of the same over-arching theory of which the quasi-classical theory is but one phase. This is a defining moment for the theory under discussion since the only natural generalization available is that from the existing quasi-classical (one-clock) model into a quasi-relativistic (two-clock) model. The quasi-relativistic disc turns out to be a stationary hyperbolic disc within which information propagates along characteristics and within which the rotation curve is necessarily exactly flat. The transition between the two components of the disc is managed by physically determined jump conditions which are represented by a set of non-linear algebraic equations. These equations predict the existence of a forbidden region within the parameter space which coincides exactly with that part of the parameter space in which there is an already observed highly significant underpopulation of galaxy discs (astro-ph/0107300).

An interpretation of Mach's Principle led us to consider if it was possible to have a globally inertial universe that was irreducibly associated with a non-trivial global matter distribution, Roscoe (GRG,2002,34,5,577-602, astro-ph/0107397). This question received a positive answer, subject to the condition that the global matter distribution is necessarily fractal, D=2. The purpose of the present paper is to show how general gravitational processes arise in this universe. We illustrate the theory by using it to model an idealized spiral galaxy. One particular subclass of solutions, corresponding to logarithmic spirals, has already been extensively tested in Roscoe A&A,1999,343,788-800 (astro-ph/0107305), and shown to resolve dynamical data over large samples of ORCs with a very high degree of statistical precision. However, this latter analysis led directly to the discovery of a major new phenomenology in spiral discs - that of discrete dynamical classes - comprehensively confirmed in Roscoe A&A,2002,385,431-453 (astro-ph/0107300) over four large independent samples of ORCs. In this paper, we analyse the theory to show how the discrete dynamical classes phenomenology has a ready explanation in terms of an algebraic consistency condition which must necessarily be satisfied. Of equal significance, we apply the theory with complete success to the detailed modelling of a sample of eight Low Surface Brightness spirals (LSBs) which, hitherto, have been succesfully modelled only by Milgrom's MOND algorithm. We are able to conclude that the essence of the MOND algorithm must be contained within the presented theory.

This first of two papers is a speculation on the nature of Mach's Principle, derived from the formulation of two questions: (A) Can a globally inertial space and time be associated with a non-trivial global matter distribution? (B) If so, what are the general properties of such a global distribution? The answers are yes subject to the condition that mass is distributed fractally with D=2. Such a prediction is supported, for example, in Joyce, Montuori, Labini astro-ph/9901290. The second paper then shows how gravitational phenomena arise as consequences of perturbations of this equilibrium D=2 distribution, and applies this theory to spiral galaxy dynamics. Specifically, the second paper derives a gravitational model for an idealized spiral, defined to have perfect cylindrical symmetry. One of the main results of this analysis is that, *in the disc*, V_rot = A/R^alpha where (A, alpha) are constants which vary between spirals. See astro-ph/0107305, which concludes that a power-law model is powerfully confirmed for ORCs as a statistical generality over large samples, and astro-ph/0107300 for new phenomonology discovered via this analysis. Thus, unlike MOND, we have only been concerned with model predictions made for idealized matter distributions, and tested over very large numbers of ORCs. But, like those of MOND, the predictions are very powerfully supported on the data. For this reason, we hypothesize that perspectives of MOND are buried within our own formalism. This remains a topic for future work.

We consider the four large optical rotation curve samples of Mathewson, Ford & Buchhorn (900 ORCs), Mathewson & Ford (1200 ORCs), Dale, Giovanelli, Haynes et al (494 ORCs) and Courteau (300 ORCs). We show how it is possible to partition any given ORC into two distinct dynamical regions - an interior region and an exterior region. Each of the samples considered is partitioned in this way, so that we are able to consider the dynamics on the exterior parts of the ORCs for each of the samples. This technique was originaly introduced by Roscoe astro-ph/0107305. We then find that the dynamics in the exterior part of spiral discs fall into one of four sharply defined classes. Taking account of the various correlations which exist between luminosity properties and dynamical properties of discs, our final phenomonological conclusion is that spiral discs in the samples considered appear constrained to occupy one of four distinct planes, F(M,S,a) = k0, k1, k2, k3, where M = absolute magnitude, S = surface brightness, "a" is a parameter computed directly from each ORC and k0, k1, k2, k3 are well defined numerical constants. These results are statistically supported at the level of virtual certainty. We suggest three possible mechanisms for this phenomonology. The overall analysis implies the existence of a second generation of Tully-Fisher methods in which the standard TF relation is augmented with a second relation defining where on an ORC the linewidth should be measured.

The 900 ORCs of Mathewson, Ford & Buchhorn which were folded by Persic & Salucci (astro-ph/9502091) are analysed. A technique for removing that part of an ORC most dominated by the central region it introduced (see astro-ph/0107300 for a far more comprehensive discussion of this). The analysis then shows that a power-law model of the form V = A R^b, where the parameters (A,b) differ between ORCs, provides an extremely good description of the generic behaviour of ORCs when these central regions are removed. It is shown that a very powerful correlation exists between the free parameters (A,b), and that parameter A can be very accurately modelled in terms of the luminosity properties of galaxies. The net conclusion is that the power-law model is, at the very least, an extremely good first-order model for the behaviour of disc dynamics in the outer regions of ORCs.