We present solutions of force-free pulsar magnetospheres coupled with a uniform external magnetic field aligned with the dipole magnetic moment of the pulsar. The inclusion of the uniform magnetic field has the following consequences: The equatorial current sheet is truncated to a finite length, the fraction of field lines that are open increases, and the open field lines are confined within a cylindrical surface instead of becoming radial. A strong external magnetic field antiparallel to the dipole allows for solutions where the pulsar magnetic field is fully enclosed within an ellipsoidal surface. Configurations of fully enclosed or confined magnetospheres may appear in a double neutron star (DNS) where one of the components is a slowly rotating, strongly magnetised pulsar and the other a weakly magnetised millisecond pulsar. Depending on the geometry, twisted field lines from the millisecond pulsar could form an asymmetric wind. More dramatic consequences could appear closer to merger: there, the neutron star with the weaker magnetic field may undergo a stage where it alternates between an open and a fully enclosed magnetosphere releasing up $10^{38}$erg. Under favourable combinations of magnetic fields DNSs can spend up to 10$^3$ yr in the coupled phase, implying a Galactic population of $0.02$ systems. Potential observational signatures of this interaction are discussed including the possibility of powering recurring Fast Radio Bursts (FRBs). We also note that the magnetic interaction cannot have any significant effect on the spin evolution of the two pulsars.

Consideration of the entropy production in the creation of the CMB leads to a simple model of the evolution of the universe during this period which suggests a connection between the small observed acceleration term and the early inflation of a closed universe. From this we find an unexpected relationship between the Omega's of cosmology and calculate the total volume of the universe.

Lema\^{i}tre-Tolman-Bondi (LTB) solutions have traditionally been confined to systems with no pressure in which the gravity is due to massive dust, but the solutions are little changed in form if, as in cosmology, the pressure is uniform in space at each comoving time. This allows the equations of cosmology to be deduced in a manner that more closely resembles classical mechanics. It also gives some inhomogeneous solutions with growing condensations and black holes. We give criteria by which the sizes of different closed models of the universe can be compared and discuss conditions for self-closure of inhomogeneous cosmologies with a $\Lambda$-term.

We show that, provided the principal axes of the second velocity moment tensor of a stellar population are generally unequal and are oriented perpendicular to a set of orthogonal surfaces at each point, then those surfaces must be confocal quadric surfaces and the potential must be separable or Stackel. This is true under the mild assumption that the even part of the distribution function is invariant under time reversal $v_i \rightarrow -v_i$ of each velocity component. In particular, if the second velocity moment tensor is everywhere exactly aligned in spherical polar coordinates, then the potential must be of Stackel form (excepting degenerate cases where two or more of the semiaxes of ellipsoid are everywhere the same). The theorem also has consequences for alignment in cylindrical polar coordinates, which is used in the popular Jeans Anisotropic Models (JAM). We analyse data on the radial velocities and proper motions of a sample of $\sim 7400$ stars in the stellar halo of the Milky Way. We provide the distributions of the tilt angles or misalignments from the spherical polar coordinate systems. We show that in this sample the misalignment is always small (usually within $3^\circ$) for radii between 7 and 12 kpc. The velocity anisotropy is very radially biased ($\beta \approx 0.7$), and almost invariant across the volume in our study. Finally, we construct a triaxial stellar halo in a triaxial NFW dark matter halo using a made-to-measure method. Despite the triaxiality of the potential, the velocity ellipsoid of the stellar halo is nearly spherically aligned within $\sim6^\circ$ for large regions of space, particularly outside the scale radius of the stellar halo. We conclude that the second velocity moment ellipsoid can be close to spherically aligned for a much wider class of potentials than the strong constraints that arise from exact alignment might suggest.

A variation of Newton's method of mapping Kepler's orbits into orbits in the simple harmonic oscillator is shown to give Henon's Isochrone. The statistical mechanics of a micro-canonical ensemble of isochrone oscillators shows that the temperature reaches a maximum as a function of the energy and then declines to zero at the escape energy. In that declining region adding heat (energy) decreases the temperature, as occurs in star clusters. We then define the internal temperature of an ensemble of binary stars all at the same (negative) energy and show that they too get cooler when heated and hotter when cooled. When the internal temperature of a binary is less that the temperature of the stars it interacts with, then on average heat will flow into it, making it less bound and of still lower temperature. Conversely hard binaries have higher internal temperatures than the local stars, so they lose energy and become hotter and yet more strongly bound, a process invoked by Michel Henon in his explanation of star-cluster evolution. Finally we give an isochronal variation of Newton's exactly soluble N-body problem.

Large contributions to the near closure of the Universe and to the acceleration of its expansion are due to the gravitation of components of the stress-energy tensor other than its mass density. To familiarise astronomers with the gravitation of these components we conduct thought experiments on gravity, analogous to the real experiments that our forebears conducted on electricity. By analogy to the forces due to electric currents we investigate the gravitational forces due to the flows of momentum, angular momentum, and energy along a cylinder. Under tension the gravity of the cylinder decreases but the 'closure' of the 3-space around it increases. When the cylinder carries a torque the flow of angular momentum along it leads to a novel helical interpretation of Levi-Civita's external metric and a novel relativistic effect. Energy currents give gravomagnetic effects in which parallel currents repel and antiparallel currents attract, though such effects must be added to those of static gravity. The gravity of beams of light give striking illustrations of these effects and a re-derivation of light bending via the gravity of the light itself. Faraday's experiments lead us to discuss lines of force of both gravomagnetic and gravity fields. A serious conundrum arrises if Landau and Lifshitz's definition of Gravitational force is adopted.

We study gravitational waves to first and second order in amplitude in vacuum asymptotically flat spacetimes. The Einstein equations are solved to first order and these solutions are superposed to form a time-symmetric ingoing and then outgoing pulse regular everywhere. The waves are assumed to have odd-parity and a non-vanishing angular momentum which keeps them away from the axis at all times. The averaged energy of the waves is evaluated. The relevant Einstein equation is then solved to second order in the amplitude. The influence of the angular momentum of the waves on the rotation of local inertial frames with respect to the frames at great distances is analyzed. The rotation of the frames occurs even in the region around the origin where spacetime is almost flat. The rotation is without time delay as it follows from the constraint equation. The results are illustrated graphically for various values of the "harmonic index m" corresponding to azimuthal rotation and the "harmonic index l" describing the latitudinal rotation of the waves. The apparent motions of the fixed stars on the celestial sphere as seen through rotating waves from the local inertial frame are calculated and displayed.

A new formula is given for the fast linear gravitational dragging of the inertial frame within a rapidly accelerated spherical shell of deep potential. The shell is charged and is electrically accelerated by an electric field whose sources are included in the solution.

Allan Sandage was an observational astronomer who was happiest at a telescope. On Hubble's sudden death Allan Sandage inherited the programmes using the world's largest optical telescope at Palomar to determine the distances and number counts of galaxies. Over many years he greatly revised the distance scale and, on re-working Hubble's analysis, discovered the error that had led Hubble to doubt the interpretation of the galaxies' redshifts as an expansion of the universe. Sandage showed that there was a consistent age of Creation for the stars, the elements, and the Cosmos. Through work with Baade and Schwarzschild he discovered the key to the interpretation of the colour-magnitude diagrams of star clusters in terms of stellar evolution. With others he founded Galactic Archaeology, interpreting the motions and elemental abundances of the oldest stars in terms of a model for the Galaxy's formation. He published several fine atlasses and catalogues of galaxies and a definitive history of the Mount Wilson Observatory.

Exact relativistic force free fields with cylindrical symmetry are explored. Such fields are generated in the interstellar gas via their connection to pulsar magnetospheres both inside and outside their light cylinders. The possibility of much enhanced interstellar fields wound on cylinders of Solar system dimensions is discussed but these are most likely unstable.

We employ the Union compilation of Type Ia supernovae with a maximum likelihood analysis to search for a dark energy dipole. To approach this problem, we present a simple, computationally efficient, and largely model independent method. We opted to weight each SN by its observed error estimate, so poorly measured SNe that deviate substantially from the Hubble law do not produce spurious results. We find, with very low significance, a dipole in the cosmic acceleration directed roughly towards the cosmic microwave background dipole, but this is almost certainly coincidental.

Following Prendergast we study the relativistically expanding electromagnetic fields generated by an axisymmetric explosion of magnetic energy in a small volume. The magnetic field expands uniformly either within a cone or in all directions and it is therefore accompanied by an electric field. In the highly conducting plasma the charges move to annul the electric field in the frame of the moving plasma. The solutions presented are analytical and semi-analytical. We find that the time-scale for the winding up of the initial magnetic field is crucial, as short time-scales lead to strong radiant fields. Assuming a magnetic field of $10^{13}Gauss$ emerging from a magnetosphere of $10^{9}cm$ we end with a jet when confined by a pressure environment that falls more slowly than $r^{-4}$. The jet carries energy of $10^{51}erg$, which is mostly due to differential rotation at the base.

A useful crude approximation for Abelian functions is developed and applied to orbits. The bound orbits in the power-law potentials A*r^{-alpha} take the simple form (l/r)^k = 1 + e cos(m*phi), where k = 2 - alpha > 0 and 'l' and 'e' are generalisations of the semi-latus-rectum and the eccentricity. 'm' is given as a function of 'eccentricity'. For nearly circular orbits 'm' is sqrt{k}, while the above orbit becomes exact at the energy of escape where 'e' is one and 'm' is 'k'. Orbits in the logarithmic potential that gives rise to a constant circular velocity are derived via the limit of small alpha. For such orbits, r^2 vibrates almost harmonically whatever the 'eccentricity'. Unbound orbits in power-law potentials are given in an appendix. The transformation of orbits in one potential to give orbits in a different potential is used to determine orbits in potentials that are positive powers of r. These transformations are extended to form a group which associates orbits in sets of six potentials, e.g. there are corresponding orbits in the potentials proportional to r, r^{-2/3}, r^{-3}, r^{-6}, r^{4/3} and r^{-4}. A degeneracy reduces this to three, which are r^{-1}, r^2 and r^{-4} for the Keplerian case. A generalisation of this group includes the isochrone with the Kepler set.

Linear perturbations of Friedmann-Robertson-Walker universes with any curvature and cosmological constant are studied in a general gauge without decomposition into harmonics. Desirable gauges are selected as those which embody best Mach's principle: in these gauges local inertial frames can be determined instantaneously via the perturbed Einstein field equations from the distributions of energy and momentum in the universe. The inertial frames are identified by their `accelerations and rotations' with respect to the cosmological frames associated with the `Machian gauges'. In closed spherical universes, integral gauge conditions are imposed to eliminate motions generated by the conformal Killing vectors. The meaning of Traschen's integral constraint vectors is thus elucidated. For all three types of FRW universes the Machian gauges admit much less residual freedom than the synchronous or generalized harmonic gauge. Mach's principle is best exhibited in the Machian gauges in closed spherical universes. Independent of any Machian motivation the general perturbation equations and discussion of gauges are useful for cosmological perturbation theory.

We present a geometrodynamical method for determining distances to orbital streams of HI gas in the Galaxy. The method makes use of our offset from the Galactic centre and assumes that the gas comprising the stream nearly follows a planar orbit about the Galactic centre. We apply this technique to the Magellanic Stream and determine the distances to all points along it; a consistency check shows that the angular momentum is approximately constant. Applying this technique to the Large Magellanic Cloud itself gives an independent distance which agrees within its accuracy of around 10%. Relaxing the demand for exact conservation of energy and angular momentum at all points along the stream allows for an increase in orbital period between the lagging end and the front end led by the Magellanic Clouds. Similar methods are applicable to other long streams of high-velocity clouds, provided they also nearly follow planar orbits; these would allow otherwise unknown distances to be determined.

Variational principles for magnetohydrodynamics were introduced by previous authors both in Lagrangian and Eulerian form. In this paper we introduce simpler Eulerian variational principles from which all the relevant equations of barotropic magnetohydrodynamics can be derived. The variational principle is given in terms of six independent functions for non-stationary barotropic flows and three independent functions for stationary barotropic flows. This is less then the seven variables which appear in the standard equations of barotropic magnetohydrodynamics which are the magnetic field $\vec B$ the velocity field $\vec v$ and the density $\rho$. The equations obtained for non-stationary barotropic magnetohydrodynamics resemble the equations of Frenkel, Levich & Stilman \cite{FLS}. The connection between the Hamiltonian formalism introduced in \cite{FLS} and the present Lagrangian formalism (with Eulerian variables) will be discussed. Finally the relations between barotropic magnetohydrodynamics topological constants and the functions of the present formalism will be elucidated.

We consider the possibility that the Galactic neutral hydrogen stream Complex A and the stellar Orphan stream are related, and use this hypothesis to determine possible distances to Complex A and the Orphan stream, and line-of-sight velocities for the latter. The method presented uses our current knowledge of the projected positions of the streams, as well as line-of-sight velocities for Complex A, and we show that a solution exists in which the two streams share the same orbit. If Complex A and the Orphan stream are on this orbit, our calculations suggest the Orphan stream to be at an average distance of 20 kpc, with heliocentric radial velocities of approximately -110 km/s. Complex A would be ahead of the Orphan stream in the same wrap of the orbit, with an average distance of 10 kpc, which is consistent with the distance constraints determined through interstellar absorption line techniques.

The magnetic fields and energy flows in an astronomical jet described by our earlier model are calculated in detail. Though the field distribution varies with the external pressure function p(z), it depends only weakly on the other boundary conditions. Individual fieldlines were plotted; the lines become nearly vertical at the bottom and are twisted at the top. An animation of a fieldline's motion was made, which shows the line being wound up by the accretion disc's differential rotation and rising as a result of this. The distribution of Poynting flux within the jet indicates that much of the energy flows up the jet from the inside of the accretion disc but a substantial fraction flows back down to the outside.

We use Sloan Digital Sky Survey Data Release 5 photometry and spectroscopy to study a tidal stream that extends over 50 degrees in the North Galactic Cap. From the analysis of the path of the stream and the colors and magnitudes of its stars, the stream is about 20 kpc away at its nearest detection (the celestial equator). We detect a distance gradient -- the stream is farther away from us at higher declination. The contents of the stream are made up from a predominantly old and metal-poor population that is similar to the globular clusters M13 and M92. The integrated absolute magnitude of the stream stars is estimated to be M_r = -7.5. There istentative evidence for a velocity signature, with the stream moving at -40 km/s at low declinations and +100 km/s at high declinations. The stream lies on the same great circle as Complex A, a roughly linear association of HI high velocity clouds stretching over 30 degrees on the sky, and as Ursa Major II, a recently discovered dwarf spheroidal galaxy. Lying close to the same great circle are a number of anomalous, young and metal-poor globular clusters, including Palomar 1 and Ruprecht 106.

Systems whose potential energies consists of pieces that scale as r^-2 together with pieces that scale as r^2, show no violent relaxation to Virial equilibrium but may pulsate at considerable amplitude for ever. Despite this pulsation these systems form lattices when the non-pulsational `energy' is low, and these disintegrate as that energy is increased. The `specific heats' show the expected halving as the `solid' is gradually replaced by the `fluid' of independent particles. The forms of the lattices are described here for N ~ 20 and they become hexagonal close packed for large N. In the larger N limit, a shell structure is formed. Their large N behaviour is analogous to a gamma=5/3 polytropic fluid with a quasi-gravity such that every element of fluid attracts every other in proportion to their separation. For such a fluid, we study the `rotating pulsating equilibria' and their relaxation back to uniform but pulsating rotation. We also compare the rotating pulsating fluid to its discrete counter part, and study the rate at which the rotating crystal redistributes angular momentum and mixes as a function of extra heat content.

The vectorial velocity is given as a function of the position of a particle in orbit when a Newtonian central force is supplemented by an inverse cubic force as in Newton's theorem on revolving orbits. Such expressions are useful in fitting orbits to radial velocities of orbital streams. The Hamilton-Laplace-Runge-Lenz eccentricity vector is generalised to give a constant of the motion for these systems and an approximate constant for orbits in general central potentials. A related vector is found for Hooke's centred ellipse.

A broad swathe of astrophysical phenomena, ranging from tubular planetary nebulae through Herbig-Haro objects, radio-galaxy and quasar emissions to gamma-ray bursts and perhaps high-energy cosmic rays, may be driven by magnetically-dominated jets emanating from accretion disks. We give a self-contained account of the analytic theory of non-relativistic magnetically dominated jets wound up by a swirling disk and making a magnetic cavity in a background medium of any prescribed pressure, p(z). We solve the time-dependent problem for any specified distribution of magnetic flux P(R,0) emerging from the disk at z=0, with any specified disk angular velocity Omega_d(R). The physics required to do this involves only the freezing of the lines of force to the conducting medium and the principle of minimum energy.

We introduce a three independent functions variational formalism for stationary and non-stationary barotropic flows. This is less than the four variables which appear in the standard equations of fluid dynamics which are the velocity field $\vec v$ and the density $\rho$. It will be shown how in terms of our new variable the Euler and continuity equations can be integrated in the stationary case.

We examine an inconsistency in theories of violent-relaxation by Lynden-Bell and Nakamura. The inconsistency arises from the non-transitive nature of these theories: a system that under goes a violent-relaxation, relaxes and then upon an addition of energy, undergoes violent-relaxation once again would settle in an equilibrium state that is different from the one that is predicted had the system would go directly from the initial to the final states. We conclude that a proper description of the violent-relaxation process cannot be achieved by equilibrium statistical mechanics approach, but instead a dynamical theory for the coarse-grained phase-space density is needed.

Explicit expressions are found for the axisymmetric metric perturbations of the closed, flat and open FRW universes caused by toroidal motions of the cosmic fluid. The perturbations are decomposed in vector spherical harmonics on 2-spheres, but the radial dependence is left general. Solutions for general odd-parity $l$-pole perturbations are given for either angular velocities or angular momenta prescribed. In particular, in case of closed universes the solutions require a special treatment of the Legendre equation.

If perturbations beyond the horizon have the velocities prescribed everywhere then the dragging of inertial frames near the origin is suppressed by an exponential factor. However if perturbations are prescribed in terms of their angular momenta there is no such suppression. We resolve this paradox and in doing so give new explicit results on the dragging of inertial frames in closed, flat and open universe with and without a cosmological constant.

Starting from the problem of two fixed centres we find a simple derivation of its third integral in terms of the scalar product of the angular momenta about the two fixed centres. This is then generalised to find the general form of the potential in which an exact third integral exists.

We show that accretion disks with magnetic fields in them ought to make jets provided that their electrical conductivity prevents slippage and there is an ambient pressure in their surroundings. We study equilibria of highly wound magnetic structures. General Energy theorems demonstrate that they form tall magnetic towers whose height grows with every turn at a velocity related to the circular velocity in the accretion disk. The pinch effect amplifies the magnetic pressures toward the axis of the towers whose stability is briefly considered.

An electromagnetic field of simple algebraic structure is simply derived. It turns out to be the G=0 limit of the charged rotating Kerr-Newman metrics. These all have gyromagnetic ratio 2, the same as the Dirac electron. The charge and current distributions giving this high gyromagnetic ratio have charges of both signs rotating at close to the velocity of light. It is conjectured that something similar may occur in the quantum electrodynamic charge distribution surrounding the point electron.

It is argued that jet modelers have given insufficient study to the natural magneto-static configurations of field wound up in the presence of a confining general pressure. Such fields form towers whose height grows with each twist at a velocity comparable to the circular velocity of the accretion disk that turns them. A discussion of the generation of such towers is preceded by a brief history of the idea that quasars, active galaxies, and galactic nuclei contain giant black holes with accretion disks.

Optically thin cooling gas at most temperatures above 30K will make condensations by pressure pushing material into cool dense regions. This works without gravity. Cooling condensations will flatten and become planar/similarity solutions. Most star formation may start from cooling condensations - with gravity only important in the later stages. The idea that some of the dark matter could be pristine white dwarfs that condensed slowly on to planetary sized seeds without firing nuclear reactions is found lacking. However, recent observations indicate fifty times more halo white dwarfs than have been previously acknowledged; enough to make the halo fraction observed as MACHOS. A cosmological census shows that only 1% of the mass of the Universe is of known constitution.

We produce the simple formula for the radius R of a cold body of mass M, (4 pi rho/3)^{1/3} R = I.M^{1/3} / [1+ (M/Mp)^{2/3}], where I is one for white dwarfs with non-relativistic electrons, but depends on the ratio of the total mass to Mch, the Chandrasekhar mass. Mp is the mass of the planet of maximum radius and rho is the density of the material at low pressure. Mp and rho are derived in terms of fundamental constants and Mp/Mch is basically given by the three halves power of the fine structure constant. The aim of the discussion is to emphasise physical principles and to connect the masses derived to the fundamental constants, so mathematical niceties are sacrificed and replaced by interpolations between simple exact limits.

In the electromagnetic fields t^2 B = b(r/t), E = r x B/t trajectories of non-relativistic charged particles conserve (r-vt)^2. The transformation r'=r/t t'=1/t maps such trajectories into orbits in the constant magnetic field all of which conserve kinetic energy. The transformation may also be used to transform any fields obeying Curl E = B/dt, div B = 0 into others and relates the particle trajectories in them.

Although magnetic fields in interstellar matter were postulated almost fifty years ago, magnetohydrodynamic theory was then much hampered by our inability to see what the magnetic field configurations were like and, after a decade of innovative development, cynics, not without some justification, began to claim that anything can happen when magnetism and an imaginative theorist get together. Thus cosmic lightning in particular received a bad press. More recently great advances in observational techniques that we shall hear of from Title, Beck, Moran and Mirabel have enabled us to see not only the sun's magnetic field with unprecedented clarity but the fields in galaxies, quasars and microquasars are now measured and not merely figments of fertile imaginations.

Application of Newton's ideas from "Principia" gives many new results in mechanics. Here we explore the question ``What form of extra force will maintain the magnitude of a vector constant of the motion while changing its direction?''

Evolution with mass segregation and the evolution of the rotation of cores are both discussed for self-similar core collapse. Evolution with angular velocity proportional to the square root of the density is predicted. On the Dynamical Main Sequence of globular clusters the energy emission from binaries balances the energy expended in expanding the halo. Newton's exactly solved N-body problem is then given, along with recent generalisations, all of which have no violent relaxation, but a new type of statistical equilibrium is discussed. Finally, we set the creation of streams in the Galaxy's halo in the historical context of their discovery.

The stability of self-similar power-law disks to non-axially symmetrical disturbances of zero frequency (omega = 0) is discussed. It is established that marginally unstable modes either have (omega = 0) or consist of continua with modes of all frequencies becoming unstable together. A careful study of all modes close to omega = 0 yields a remarkable conclusion.

The wave functions of Boson and Fermion gases are known even when the particles have harmonic interactions. Here we generalise these results by solving exactly the N-body Schrodinger equation for potentials V that can be any function of the sum of the squares of the distances of the particles from one another in 3 dimensions. For the harmonic case that function is linear in r^2. Explicit N-body solutions are given when U(r) = -2M \hbar^{-2} V(r) = \zeta r^{-1} - \zeta_2 r^{-2}. Here M is the sum of the masses and r^2 = 1/2 M^{-2} Sigma Sigma m_I m_J ({\bf x}_I - {\bf x}_J)^2. For general U(r) the solution is given in terms of the one or two body problem with potential U(r) in 3 dimensions. The degeneracies of the levels are derived for distinguishable particles, for Bosons of spin zero and for spin 1/2 Fermions. The latter involve significant combinatorial analysis which may have application to the shell model of atomic nuclei. For large N the Fermionic ground state gives the binding energy of a degenerate white dwarf star treated as a giant atom with an N-body wave function. The N-body forces involved in these extraordinary N-body problems are not the usual sums of two body interactions, but nor are forces between quarks or molecules. Bose-Einstein condensation of particles in 3 dimensions interacting via these strange potentials can be treated by this method.

When a charged insulating spherical shell is uniformly accelerated, an oppositely directed electric field is produced inside. Outside the field is the Born field of a uniformly accelerated charge, modified by a dipole. Radiation is produced. When the acceleration is annulled by the nearly uniform gravity field of an external shell with a 1 + beta cos theta surface distribution of mass, the differently viewed Born field is static and joins a static field outside the external shell; no radiation is produced. We discuss gravitational analogues of these phenomena. When a massive spherical shell is accelerated, an untouched test mass inside experiences a uniform gravity field and accelerates parallelly to the surrounding shell. In the strong gravity regime we illustrate these effects using exact conformastatic solutions of the Einstein-Maxwell equations with charged dust. We consider a massive charged shell on which the forces due to nearly uniform electrical and gravitational fields balance. Both fields are reduced inside by the ratio of the g_00 inside the shell to that away from it. The acceleration of a free test particle, relative to a static observer, is reduced correspondingly. We give physical explanations of these effects.

Starting from Antonov's discovery that there is no maximum to the entropy of a gravitating system of point particles at fixed energy in a spherical box if the density contrast between centre and edge exceeds 709, we review progress in the understanding of gravitational thermodynamics. We pinpoint the error in the proof that all systems have positive specific heat and say when it can occur. We discuss the development of the thermal runaway in both the gravothermal catastrophe and its inverse. The energy range over which microcanonical ensembles have negative heat capacity is replaced by a first order phase transition in the corresponding canonical ensembles. We conjecture that all first order phase transitions may be viewed as caused by negative heat capacities of units within them. We find such units in the theory of ionisation, chemical dissociation and in the Van der Waals gas so these concepts are applicable outside the realm of stars, star clusters and black holes.

We outline reasoning that led to the current theory of quasars and look at George Contopoulos's place in the long history of the N-body problem. Following Newton we find new exactly soluble N-body problems with multibody forces and give a strange eternally pulsating system that in its other degrees of freedom reaches statistical equilibrium.

We solve the N-body problems in which the total potential energy is any function of the mass-weighted root-mean-square radius of the system of N point masses. The fundamental breathing mode of such systems vibrates non-linearly for ever. If the potential is supplemented by any function that scales as the inverse square of the radius there is still no damping of the fundamental breathing mode. For such systems a remarkable new statistical equilibrium is found for the other coordinates and momenta, which persists even as the radius changes continually.

Slowly rotating collapsing spherical shells have flat spaces inside and the inertial frames there rotate at omega_s(t) relative to infinity. As first shown by Lindblom & Brill the inertial axes within the shell rotate rigidly without time delays from one point to another. Although the rotation rate of the inertial axes is changing the axes are inertial so, relative to them, there is neither an omega dot_s times r (Euler) fictitious force nor any other. However, Euler and other fictitious forces arise in the frame which is at rest with respect to infinity. An observer at the center who looks in one direction (phi = 0, theta = pi/2 say) fixed to infinity will see that the sky appears to rotate and can compare its apparent rotation with those of the local inertial frame and of the shell itself. By contrast in the electromagnetic analogue there is a time delay in the propagation of the magnetic field inside a rotating collapsing charged shell in flat space. We demonstrate this time delay by devising a null experiment in which the Larmor precession of a charged oscillator would be exactly cancelled by the rotation of the inertial frame but for the delay. In the ``combined'' problem of a collapsing charged shell we show that, due to the coupling of electromagnetic and gravitational perturbations, the instantaneous rotation of inertial frames inside the shell can be caused by pure electric currents in a non-rotating shell.

We present a useful analytic approximation to the solution of the Lane-Emden equation for infinite polytropic index - the isothermal sphere. The optimized expression obtained for the density profile can be accurate to within 0.04% within 5 core-radii and to 0.1% within 10 core-radii.

A general method is presented for constructing distribution functions for flat systems whose surface density and Toomre's Q number profile is given. The purpose of these functions is to provide plausible galactic models and assess their critical stability with respect to global non axi-symmetric modes. The derivation may be carried out for an azimuthal velocity distribution (or a given specific energy distribution) which may either be observed or chosen to match a specified temperature profile. Distribution functions describing stable models with realistic velocity distributions for power law disks, the Isochrone and the Kuzmin disks are provided. Specially simple inversion formulae are also given for finding distribution functions for flat systems whose surface densities are known.

Complete sequences of new analytic solutions of Einstein's equations which describe thin super massive disks are constructed. These solutions are derived geometrically. The identification of points across two symmetrical cuts through a vacuum solution of Einstein's equations defines the gradient discontinuity from which the properties of the disk can be deduced. The subset of possible cuts which lead to physical solutions is presented. At large distances, all these disks become Newtonian, but in their central regions they exhibit relativistic features such as velocities close that of light, and large redshifts. Sections with zero extrinsic curvature yield cold disks. Curved sections may induce disks which are stable against radial instability. The general counter rotating flat disk with planar pressure tensor is found. Owing to gravomagnetic forces, there is no systematic method of constructing vacuum stationary fields for which the non-diagonal component of the metric is a free parameter. However, all static vacuum solutions may be extended to fully stationary fields via simple algebraic transformations. Such disks can generate a great variety of different metrics including Kerr's metric with any ratio of a to m. A simple inversion formula is given which yields all distribution functions compatible with the characteristics of the flow, providing formally a complete description of the stellar dynamics of flattened relativistic disks. It is illustrated for the Kerr disk.

A relative mechanics with no absolute space is shown to be equivalent to Newtonian mechanics applied in a universe of zero net angular momentum. Closed spaces in General Relativity have no angular momentum and shrivel to one point as the mass-energy contained tends to zero, so obeying Mach's principle on the origin of inertia.

We have visually examined twelve Palomar red Plates for galaxies at low Galactic latitude b, where the Supergalactic Plane (SGP) is crossed by the Galactic Plane (GP), at Galactic longitude l ~135 degrees. The catalogue consists of 2575 galaxy candidates, of which 462 have major axis diameters d >= 0.8 arc min (uncorrected for extinction). Galaxy candidates can be identified down to |b| ~ 0 degrees. One of our galaxy candidates (J24 = Dwingeloo 1) has recently been discovered independently in 21cm by Kraan-Korteweg et al. (1994) as a nearby galaxy. Comparisons with the structures seen in the IRAS and UGC catalogues are made. We compare the success rate of identifying galaxies using the IRAS Point Source Catalogue under different colour selection criteria. The criteria that require both the 60 micron and 100 micron fluxes to be of high quality, have the highest probability of selecting a galaxy (with d >= 0.6 arc min), but at the expense of selecting a smaller number of galaxies in total.

We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.