In order to gain a better understanding of the influence of cosmic rays (CRs) and magnetic fields in the disk-halo interface of edge-on spiral galaxies, we investigate the radio continuum halo, the magnetic field, and the transport processes of the CRs of the edge-on spiral galaxy NGC 4217 using CHANG-ES radio data at two frequencies, 6 GHz (C-band) and 1.5 GHz (L-band), and supplemental LOFAR data at 150 MHz and X-ray Chandra data. NGC 4217 shows a large-scale X-shaped magnetic field structure, covering a major part of the galaxy with a mean total magnetic field strength in the disk of 9 micro Gauss (via equipartition). Using rotation measure synthesis (RM-synthesis) at C-band, we found that the direction of the disk magnetic field is pointing inward. A helical outflow structure is furthermore present in the northwestern part of the galaxy, which is extended nearly 7 kpc into the halo. More polarized emission is observed on the approaching side of the galaxy. With a simplified galaxy disk model, we are able to explain that finding and predict that roughly 75% of edge-on spiral galaxies will show higher polarized intensity on the approaching side. Many loop and shell structures are found throughout the galaxy in total intensity at C-band. A superbubble-like structure is prominent in total and polarized intensity, as well as in Halpha and optical dust filaments, being a possible result of concentrated star formation in the disk. The flux density contribution of the disk in comparison to the halo decreases toward lower frequencies. Total intensity profiles at the three radio frequencies were fit with two-component exponential functions. The frequency dependence of the resulting scale heights between C-band and L-band suggests advection to be the main CR transport process. The 1D CR transport modeling (SPINNAKER) shows that advection appears to be more important than diffusion.

NGC 4631 is an interacting galaxy which exhibits one of the largest gaseous halos observed among edge-on galaxies. We aim to examine the synchrotron and polarization properties of its disk and halo emission with new radio continuum data. Radio continuum observations of NGC 4631 were performed with the Karl G. Jansky Very Large Array at C-band (5.99 GHz) in the C & D array configurations, and at L-band (1.57 GHz) in the B, C, & D array configurations. The Rotation Measure Synthesis algorithm was utilized to derive the polarization properties. We detected linearly polarized emission at C-band and L-band. The magnetic field in the halo is characterized by strong vertical components above and below the central region of the galaxy. The magnetic field in the disk is only clearly seen in the eastern side of NGC 4631, where it is parallel to the plane of the major axis of the galaxy. We detected for the first time a large-scale, smooth Faraday depth pattern in a halo of an external spiral galaxy, which implies the existence of a regular (coherent) magnetic field. A quasi-periodic pattern in Faraday depth with field reversals was found in the northern halo of the galaxy. The field reversals in the northern halo of NGC 4631, together with the observed polarization angles, indicate giant magnetic ropes (GMRs) with alternating directions. To our knowledge, this is the first time such reversals are observed in an external galaxy.

NGC 4631 is an interacting galaxy that exhibits one of the largest, gaseous halos observed among edge-on galaxies. We aim to examine the synchrotron and cosmic-ray propagation properties of its disk and halo emission with new radio continuum data. Radio continuum observations of NGC 4631 were performed with the Karl G. Jansky Very Large Array at C-band (5.99 GHz) in the C and D array configurations, and at L-band (1.57 GHz) in the B, C, and D array configurations. Complementary observations of NGC 4631 with the Effelsberg telescope were performed at 1.42 and 4.85 GHz. The interferometric total intensity data were combined with the single-dish Effelsberg data in order to recover the missing large-scale total power emission. The thermal and nonthermal components of the total radio emission were separated by estimating the thermal contribution through the extinction-corrected H$\alpha$ emission. The H$\alpha$ radiation was corrected for extinction using a linear combination of the observed H$\alpha$ and 24 $\mu$m data. NGC 4631 has a global thermal fraction at 5.99 (1.57) GHz of 14$\pm$3% (5.4$\pm$1.1%). The mean scale heights of the total emission in the radio halo (thick disk) at 5.99 (1.57) GHz are $1.79\pm0.54$ kpc ($1.75\pm0.27$ kpc) and have about the same values for the synchrotron emission. The total magnetic field of NGC 4631 has a mean strength of $\rm{\langle B_{eq}\rangle} \simeq 9 \rm{\mu G}$ in the disk, and a mean strength of $\rm{\langle B_{eq}\rangle}~\simeq 7~\rm{\mu G}$ in the halo. We also studied a double-lobed background radio galaxy southwest of NGC 4631, which is an FR~II radio galaxy according to the distribution of spectral index across the lobes. From the halo scale heights we estimated that the radio halo is escape-dominated with convective cosmic ray propagation, and conclude that there is a galactic wind in the halo of NGC 4631.

Rotation measure (RM) synthesis maps of NGC 4631 show remarkable sign reversals with distance from the minor axis in the northern halo of the galaxy on kpc scales. To explain this new phenomenon, we solve the dynamo equations under the assumption of scale invariance and search for rotating logarithmic spiral solutions. Solutions for velocity fields representing accretion onto the disk, outflow from the disk, and rotation-only in the disk are found that produce RM with reversing signs viewed edge-on. Model RM maps are created for a variety of input parameters using a Faraday screen and are scaled to the same amplitude as the observational maps. Residual images are then made and compared to find a best fit. Solutions for rotation-only, i.e. relative to a pattern uniform rotation, did in general, not fit the observations of NGC 4631 well. However, outflow models did provide a reasonable fit to the magnetic field. The best results for the specific region that was modelled in the northern halo are found with accretion. Since there is abundant evidence for both winds and accretion in NGC 4631, this modelling technique has the potential to distinguish between the dominant flows in galaxies.

The CHANG-ES (Continuum Halos in Nearby Galaxies) survey of 35 nearby edge-on galaxies is revealing new and sometimes unexpected and startling results in their radio continuum emission. The observations were in wide bandwidths centered at 1.6 and 6.0 GHz. Unique to this survey is full polarization data showing magnetic field structures in unprecedented detail, resolution and sensitivity for such a large sample. A wide range of new results are reported here, some never before seen in any galaxy. We see circular polarization and variability in active galactic nuclei (AGNs), in-disk discrete features, disk-halo structures sometimes only seen in polarization, and broad-scale halos with reversing magnetic fields, among others. This paper summarizes some of the CHANG-ES results seen thus far. Released images can be found at https://www.queensu.ca/changes.

We propose a new systematic method of studying correlations between parameters that describe an astronomical (or any) physical system. We recall that behind Dimensionless scaling laws in complex, self-interacting physical objects lies a rigorous theorem of Dimensional analysis, known widely as the Buckingham theorem. Once a {\it catalogue} of properties and forces that define an object or physical system is established, the theorem allows one to select a complete set of Dimensionless quantities or {\it Numbers} on which structure must depend. The internal structure takes the form of a functionally defined manifold in the space of these Numbers. Simple and familiar examples are discussed by way of introduction. Correlations in properties of astronomical objects can be sought either through the constancy of these Numbers or between pairs of the Numbers. In either case, within errors, the functional dependences take on an absolute numerical character. As our principal application, we study a well defined sample of galaxies in order to reveal the implied Tully Fisher and Baryonic Tully Fisher relations. We find that $L\,\propto\,v_{rot}^4$ for the former and $M_b\,\propto\,v_{rot}^3$ for the latter, suggesting that these relations may have different causal origins.

This paper extends the results of an earlier paper on scale invariant galactic spiral magnetic modes to time dependent, scale invariant, spirals. The examples given are all exact in that they may be described in terms of hypergeometric functions. We restrict the discussion to an infinitely conducting medium in order to avoid earlier approximation, which limited the solutions to cones lying within about twenty degrees to the plane. The magnetic disc spirals, `X type" poloidal fields, and the recent discovery of rotation measure screens in edge-on galactic halos were all recovered in such solutions.

We detect 5 galaxies in the CHANG-ES (Continuum Halos in Nearby Galaxies -- an EVLA Survey) sample that show circular polarization (CP) at L-band in our high resolution data sets. Two of the galaxies (NGC~4388 and NGC~4845) show strong Stokes $V/I\,\equiv\,m_C\,\sim\,2$\%, two (NGC~660 and NGC~3628) have values of $m_C\sim \,0.3$\%, and NGC~3079 is a marginal detection at $m_C\sim \,0.2$\%. The two strongest $m_C$ galaxies also have the most luminous X-ray cores and the strongest internal absorption in X-rays. We have expanded on our previous Faraday conversion interpretation and analysis and provide analytical expressions for the expected $V$ signal for a general case in which the cosmic ray electron energy spectral index can take on any value. We provide examples as to how such expressions could be used to estimate magnetic field strengths and the lower energy cutoff for CR electrons. Four out of our detections are {\it resolved}, showing unique structures, including a {\it jet} in NGC~4388 and a CP `conversion disk' in NGC~4845. The conversion disk is inclined to the galactic disk but is perpendicular to a possible outflow direction. Such CP structures have never before been seen in any galaxy to our knowledge. None of the galaxy cores show linear polarization at L-band. Thus CP may provide a unique probe of physical conditions deep into radio AGNs.

We give a selection of exact dynamos in axial symmetry on a galactic scale. These include some steady examples, at least one of which is wholly analytic in terms of simple functions and has been discussed elsewhere. Most solutions are found in terms of special functions, such as associated Lagrange or hypergeometric functions. They may be considered exact in the sense that they are known to any desired accuracy in principle. The new aspect developed here is to present scale invariant solutions with zero resistivity that are self-similar in time. The time dependence is either a power law or an exponential factor, but since the geometry of the solution is self-similar in time we do not need to fix a time to study it. Several examples are discussed. Our results demonstrate (without the need to invoke any other mechanisms) X-shaped magnetic fields and (axially symmetric) magnetic spiral arms (both of which are well observed and documented) and predict reversing rotation measures in galaxy halos (now observed in the CHANG-ES sample) as well as the fact that planar magnetic spirals are lifted into the galactic halo.

We study the global consequences in the halos of spiral galaxies of the steady, axially symmetric, mean field dynamo. We use the classical theory but add the possibility of using the velocity field components as parameters in addition to the helicity and diffusivity. The analysis is based on the simplest version of the theory and uses scale-invariant solutions. The velocity field (subject to restrictions) is a scale invariant field in a `pattern' frame, in place of a full dynamical theory. The `pattern frame' of reference may either be the systemic frame or some rigidly rotating spiral pattern frame. One type of solution for the magnetic field yields off-axis, spirally wound, magnetic field lines. These predict sign changes in the Faraday screen rotation measure in every quadrant of the halo of an edge-on galaxy. Such rotation measure oscillations have been observed in the CHANG-ES survey.

We seek the conditions for a {\it steady} mean field galactic dynamo. The parameter set is reduced to those appearing in the $\alpha^2$ and $\alpha/\omega$ dynamo, namely velocity amplitudes, and the ratio of sub-scale helicity to diffusivity. The parameters can be allowed to vary on conical spirals. We analyze the mean field dynamo equations in terms of scale invariant logarithmic spiral modes and special exact solutions. Compatible scale invariant gravitational spiral arms are introduced and illustrated in an appendix, but the detailed dynamical interaction with the magnetic field is left for another work. As a result of planar magnetic spirals `lifting' into the halo, multiple sign changes in average rotation measures forming a regular pattern on each side of the galactic minor axis, are predicted. Such changes have recently been detected in the CHANG-ES survey.

Tidal disruption events (TDEs) occur when a star or sub-stellar object passes close enough to a galaxy's supermassive black hole to be disrupted by tidal forces. NGC 4845 (d=17 Mpc) was host to a TDE, IGR J12580+0134, detected in November 2010. Its proximity offers us a unique close-up of the TDE and its aftermath. We discuss new Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (JVLA) observations, which show that the radio flux from the active nucleus created by the TDE has decayed in a manner consistent with predictions from a jet-circumnuclear medium interaction model. This model explains the source's broadband spectral evolution, which shows a spectral peak that has moved from the submm (at the end of 2010) to GHz radio frequencies (in 2011-2013) to <1 GHz in 2015. The milliarcsecond-scale core is circularly polarized at 1.5 GHz but not at 5 GHz, consistent with the model. The VLBA images show a complex structure at 1.5 GHz that includes an east west extension ~40 milliarcsec (3 pc) long as well as a resolved component 52 milliarcsec (4.1 pc) northwest of the flat-spectrum core, which is all that can be seen at 5 GHz. If ejected in 2010, the NW component must have had v=0.96 c over five years. However, this is unlikely, as our model suggests strong deceleration to speeds < 0.5c within months and a much smaller, sub-parsec size. In this interpretation, the northwest component could have either a non-nuclear origin or be from an earlier event.

We report on C-band (5 - 7 GHz) observations of the galaxy, NGC~2992, from the CHANG-ES sample. This galaxy displays an embedded nuclear double-lobed radio morphology within its spiral disk, as revealed in linearly polarized emission but {\it not} in total intensity emission. The radio lobes are kpc-sized, similar to what has been observed in the past for other Seyfert galaxies, and show ordered magnetic fields. NGC~2992 has shown previous evidence for AGN-related activity, but not the linearly polarized radio features that we present here. We draw attention to this galaxy as the first clear example (and prototype) of bipolar radio outflow that is revealed in linearly polarized emission only. Such polarization observations, which are unobscured by dust, provide a new tool for uncovering hidden weak AGN activity which may otherwise be masked by brighter unpolarized emission within which it is embedded. The radio lobes observed in NGC~2992 are interacting with the surrounding interstellar medium and offer new opportunities to investigate the interactions between nuclear outflows and the ISM in nearby galaxies. We also compare the radio emission with a new CHANDRA X-ray image of this galaxy. A new CHANG-ES image of NGC~3079 is also briefly shown as another example as to how much more obvious radio lobes appear in linear polarization as opposed to total intensity.

We present an analytic model of a magnetized galactic halo surrounding a Mestel gravitating disc. The magnetic field is taken to be in energy equipartition with the pressure dominant rotating halo gas ({\it not} with the cosmic rays), and the whole system is in a steady state. A more flexible `anisotropic equipartition' model is also explored. A definite pressure law is required to maintain the equilibrium, but the halo density is constant. The velocity/magnetic system is scale-free. The objective is to find the rotational velocity lag in such a halo. The magnetic field is not force-free so that angular momentum may be transported from the halo to the intergalactic medium. We find that the `X'-shaped structure observed for halo magnetic fields can be obtained together with a simple analytic formula for the rate of decline of the velocity with height $z$. The formula also predicts the change in lag with radius, $r$.

We have observed the Virgo Cluster spiral galaxy, NGC~4845, at 1.6 and 6 GHz using the Karl G. Jansky Very Large Array, as part of the `Continuum Halos in Nearby Galaxies -- an EVLA Survey' (CHANG-ES). The source consists of a bright unresolved core with a surrounding weak central disk (1.8 kpc diameter). The core is variable over the 6 month time scale of the CHANG-ES data and has increased by a factor of $\approx$ 6 since 1995. The wide bandwidths of CHANG-ES have allowed us to determine the spectral evolution of this core which peaks {\it between} 1.6 and 6 GHz (it is a GigaHertz-peaked spectrum source).We show that the spectral turnover is dominated by synchrotron self-absorption and that the spectral evolution can be explained by adiabatic expansion (outflow), likely in the form of a jet or cone. The CHANG-ES observations serendipitously overlap in time with the hard X-ray light curve obtained by Nikolajuk \& Walter (2013) which they interpret as due to a tidal disruption event (TDE) of a super-Jupiter mass object around a $10^5\, M_\odot$ black hole. We outline a standard jet model, provide an explanation for the observed circular polarization, and quantitatively suggest a link between the peak radio and peak X-ray emission via inverse Compton upscattering of the photons emitted by the relativistic electrons. We predict that it should be possible to resolve a young radio jet via VLBI as a result of this nearby TDE.

We present the first maps of NGC 3044 and NGC 4157 at $\lambda\,450 \mu$m and $\lambda\,850 \mu$m from the JCMT as well as the first maps at 617 MHz from the GMRT. High latitude emission has been detected in both the radio continuum and sub-mm for NGC 3044 and in the radio continuum for NGC 4157, including several new features. The dust spectrum at long wavelengths required fitting with a two-temperature model for both galaxies, implying the presence of cold dust. Dust masses are $M_d\,=\,1.6\,\times\,10^8 M_\odot$ and $M_d\,=\,2.1\,\times\,10^7 M_\odot$ for NGC 3044 and NGC 4157, respectively, and are dominated by the cold component. There is a clear correlation between the 617 MHz and $\lambda\,850 \mu$m emission in the two galaxies. In the case of NGC 3044, this implies a relation between the non-thermal synchrotron emission and cold dust. The 617 MHz component represents an integration of massive star formation over the past $10^{7-8}$ yr and the $\lambda \,850 \mu$m emission represents heating from the diffuse interstellar radiation field (ISRF). The 617 MHz -- $\lambda\,850 \mu$m correlation improves when a smoothing kernel is applied to the $\lambda\,850 \mu$m data to account for differences between the CR electron diffusion scale and the mean free path of an ISRF photon to dust. The best-fit relation is $L_{617_{\rm MHz}}\,\propto\,{L_{850 \mu{\rm m}}}^{2.1\,\pm\,0.2}$ for NGC 3044. If variations in the cold dust emissivity are dominated by variations in dust density, and the synchrotron emission depends on magnetic field strength (a function of gas density) as well as CR electron generation (a function of massive star formation rate and therefore density via the Schmidt law) then the expected correlation for NGC 3044 is $L_{617_{\rm MHz}}\,\propto\,{L_{850 \mu{\rm m}}}^{2.2}$, in agreement with the observed correlation.

We study transient spiral structures in an isothermal, thin, galactic disc. We find no instability that can grow everywhere from infinitesimal disturbances, but spiral structure does grow in the disc due to an arbitrarily strong, asymmetric, central source. An initially finite spiral undergoes transient swing amplification as it is gradually wound-up by differential rotation. An independent sequence in negative time describes a leading spiral swinging to a trailing spiral. The dynamical coupling is established between the swinging potential and the arm particles, by ensuring that this potential constrains a locally rotating distribution function centred on the arms. This swing amplification propagates in radius at the constant rotational speed of the disc, and leaves behind an exponential density decline in space and time.

We begin by recalling the isothermal, collisionless, disc-halo. The disc component is the Mestel disc. Subsequently we introduce spiral arms to such an isothermal disc-halo system that are co-moving in the mean with an axi-symmetric background. These correspond to a similar disturbance in the halo, which is comprised of spiral structures on cones. The arms are necessarily transient due to the differential winding in the disc and their gradual destruction is described. Although the spiral potentials are weak compared to the axi-symmetric potential the arms are not propagating waves on the background, but rather co-move with it. They have an effect disproportionate to their relative magnitude on the gas distribution in the disc. The gas accumulates on the outside leading edge of the 'stellar' arm and an arm-interarm modulation of up to 100% is possible. Compatible isothermal, scale-free, distribution functions are found either exactly or approximately for all of the collisionless components of the disc-halo system. Repeated episodes of winding arms can produce an exponential disc.

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile the most comprehensive data set yet obtained for the study of radio halo properties. This is the first survey of radio halos to include all polarization products. In this first paper, we outline the scientific motivation of the survey, the specific science goals, and the expected improvements in noise levels and spatial coverage from the survey. Our goals include investigating the physical conditions and origin of halos, characterizing cosmic ray transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal radio emission with high-energy gamma-ray models. The sample size allows us to search for correlations between radio halos and other properties, including environment, star formation rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction steps and present the first results of the survey for the galaxy, NGC 4631.

We present the first results from the CHANG-ES survey, a new survey of 35 edge-on galaxies to search for both in-disk as well as extra-planar radio continuum emission. The motivation and science case for the survey are presented in a companion paper (Paper I). In this paper (Paper II), we outline the observations and data reduction steps required for wide-band calibration and mapping of EVLA data, including polarization, based on C-array test observations of NGC 4631. With modest on-source observing times (30 minutes at 1.5 GHz and 75 minutes at 6 GHz for the test data) we have achieved best rms noise levels of 22 and 3.5 $\mu$Jy beam$^{-1}$ at 1.5 GHz and 6 GHz, respectively. New disk-halo features have been detected, among them two at 1.5 GHz that appear as loops in projection. We present the first 1.5 GHz spectral index map of NGC 4631 to be formed from a single wide-band observation in a single array configuration. This map represents tangent slopes to the intensities within the band centered at 1.5 GHz, rather than fits across widely separated frequencies as has been done in the past and is also the highest spatial resolution spectral index map yet presented for this galaxy. The average spectral index in the disk is $\bar\alpha_{1.5 GHz}\,=\,-0.84\,\pm\,0.05$ indicating that the emission is largely non-thermal, but a small global thermal contribution is sufficient to explain a positive curvature term in the spectral index over the band. Two specific star forming regions have spectral indices that are consistent with thermal emission. Polarization results (uncorrected for internal Faraday rotation) are consistent with previous observations and also reveal some new features. On broad scales, we find strong support for the notion that magnetic fields constrain the X-ray emitting hot gas.

In this paper we suggest the existence in the central regions of spiral galaxies of collisionless, scale-free, rigidly rotating, self-gravitating discs with spiral symmetry. Such discs must be truncated at a finite radius, and they must be stabilized and rendered self-similar by a suitable halo. The halo and the rotating disc share the self-similar class and must form together to arrive at the suggested state. We make comparisons with the well-known rigidly rotating, Kalnajs discs; one of which is axi-symmetric and finite while the other is infinite and decomposed into spiral modes. We find the self-consistent, self-similar, distribution functions in one and two dimensions in a rigidly rotating, collisionless system. In the case of two dimensions we deduce the self-consistency condition for discrete spiral arms. We give an estimate of the disturbance created in the halo by the presence of the disc, and argue that the halo itself should be close to self-similarity. A very weak cusp in the halo may be necessary. The necessary spatial coincidence of the halo results in a kind of disc-halo `conspiracy'. Finally the disc equations are formulated in `spiral' coordinates, and the passage to an approximately discrete `line spiral' is given as an example. Although in two dimensions the collisionless particles enter and leave the arms in non-linear epicycles, they move approximately parallel to the arms in the line spiral limit. The spiral pattern is however in rigid rotation. Aperiodic spiral arms are suggested wherein discontinuities may be coarse-grained to appear as collisionless shocks.

In this paper, we study the distribution functions that arise naturally during self-similar radial infall of collisionless matter. Such matter may be thought of either as stars or as dark matter particles. If a rigorous steady state is assumed, then the system is infinite and is described by a universal distribution function given the self-similar index. The steady logarithmic potential case is exceptional and yields the familiar Gaussian for an infinite system with an inverse-square density profile. We show subsequently that for time-dependent radial self-similar infall, the logarithmic case is accurately described by the Fridmann and Polyachenko distribution function. The system in this case is finite but growing. We are able to embed a central mass in the universal steady distribution only by iteration, except in the case of massless particles. The iteration yields logarithmic corrections to the massless particle case and requires a `renormalization' of the central mass. A central spherical mass may be accurately embedded in the Fridmann and Polyachenko growing distribution however. Some speculation is given concerning the importance of radial collisionless infall in actual galaxy formation.

Aims. In this paper we continue our study of density cusps that may contain central black holes. Methods. We recall our attempts to use distribution functions with a memory of self-similar relaxation, but mostly they apply only in restricted regions of the global system. We are forced to consider related distribution functions that are steady but not self-similar. Results. One remarkably simple distribution function that has a filled loss cone describes a bulge that transits from a near black hole domain to an outer 'zero flux' regime where$\rho\propto r^{-7/4}$. The transition passes from an initial inverse square profile through a region having a 1/r density profile. The structure is likely to be developed at an early stage in the growth of a galaxy. A central black hole is shown to grow exponentially in this background with an e-folding time of a few million years. Conclusions. We derive our results from first principles, using only the angular momentum integral in spherical symmetry. The initial relaxation probably requires bar instabilities and clump-clump interactions.

Aims: In this paper we study density cusps that may contain central black holes. The actual co-eval self-similar growth would not distinguish between the central object and the surroundings. Methods: To study the environment of a growing black hole we seek descriptions of steady `cusps' that may contain a black hole and that retain at least a memory of self-similarity. We refer to the environment in brief as the `bulge' and on smaller scales, the `halo'. Results: We find simple descriptions of the simulations of collisionless matter by comparing predicted densities, velocity dispersions and distribution functions with the simulations. In some cases central point masses may be included by iteration. We emphasize that the co-eval self-similar growth allows an explanation of the black hole bulge mass correlation between approximately similar collisionless systems. Conclusions: We have derived our results from first principles assuming adiabatic self-similarity and either self-similar virialisation or normal steady virialisation. We conclude that distribution functions that retain a memory of self-similar evolution provide an understanding of collisionless systems. The implied energy relaxation of the collisionless matter is due to the time dependence. Phase mixing relaxation may be enhanced by clump-clump interactions.

Cosmological simulations consistently predict specific properties of dark matter halos, but these have not yet led to a physical understanding that is generally accepted. This is especially true for the central regions of these structures. Recently two major themes have emerged. In one, the dark matter halo is primarily a result of the sequential accretion of primordial structure (ie `Nature'); while in the other, dynamical relaxation (ie `Nurture') dominates at least in the central regions. Some relaxation is however required in either mechanism. In this paper we accept the recently established scale-free sub-structure of halos as an essential part of both mechanisms. Consequently; a simple model for the central relaxation based on a self-similar cascade of tidal interactions, is contrasted with a model based on the accretion of adiabatically self-similar, primordial structure. We conclude that a weak form of this relaxation is present in the simulations, but that is normally described as the radial orbit instability.

In this paper we reconsider a series expansion for a dark matter distribution function in the spherically symmetric anisotropic limit. We show here that the expansion may be renormalized so that the series does converge in time to an estimate of the steady state distribution function in the central regions. Subsequently we use this distribution function to discuss the nature of the central equilibrium and, by invoking stationarity of Boltzmann's H function as a measure of (thermodynamic) relaxation, we calculate the adiabatic variation in the local logarithmic slope of the mass density. Similarly the pseudo (phase-space) density variation with radius is calculated. These are compared to empirical fitting functions. There is general agreement on the inner part of the logarithmic slope of the density and of the inner profile of the pseudo-density power law, but coincident continuity with the outer power-laws is not yet achieved. Finally some suggestions are made regarding the actual microphysics acting during the non-equilibrium approach to relaxation. In particular a cascade regime is identified.

Simulated dark matter profiles are often modelled as a `NFW' density profile rather than a single power law. Recently, attention has turned to the rather rigorous power-law behaviour exhibited by the `pseudo phase-space density' of the dark matter halo, which is defined dimensionally in terms of the local density and velocity dispersion of the dark matter particles. The non-power-law behaviour of the density profile is generally taken to exclude simple scale-free, in-fall models; however the power-law behaviour of the `pseudo-density' is a counter indication. We argue in this paper that both behaviours may be at least qualitatively understood in terms of a dynamically evolving self-similarity, rather than the form for self-similar infall that is fixed by cosmological initial conditions. The evolution is likely due to collective relaxation such as that provided by the radial-orbit instability on large scales. We deduce, from a distribution function given by first order coarse-graining, both the NFW-type density profile and the power-law pseudo-density profile. The results are not greatly sensitive to variation about 3 in the power of the velocity dispersion used in the definition of the phase space pseudo-density. We suggest that the power 2 may create the more physical quantity, whose deviations from a power-law are a diagnostic of incomplete relaxation.

The radial orbit instability drives dark matter halos toward a universal structure. This conclusion, first noted by Huss, Jain, and Steinmetz, is explored in detail through a series of numerical experiments involving the collapse of an isolated halo into the non-linear regime. The role played by the radial orbit instability in generating the density profile, shape, and orbit structure is carefully analyzed and, in all cases, the instability leads to universality independent of initial conditions. New insights into the underlying physics of the radial orbit instability are presented.

We compare the conclusions reached using the coarse-graining technique employed by Henriksen (2004) for a one degree of freedom (per particle) collisionless system, to those presented in a paper by Binney (2004) based on an exact one degree of freedom model. We find agreement in detail but in addition we show that the isolated 1D system is self-similar and therefore unrelaxed. Fine graining of this system recovers much less prominent wave-like structure than in a spherically symmetric isotropic 3D system. The rate of central flattening is also reduced in the 1D system. We take this to be evidence that relaxation of collisionless systems proceeds ultimately by way of short wavelength Landau damping. N-body systems, both real and simulated, can be trapped in an incompletely relaxed state because of a break in the cascade of energy to small scales. This may be due to the rapid dissipation of the small scale oscillations in an isolated system, to the existence of conserved quantities such as angular momentum, or to the failure in simulations to resolve sub-Jeans length scales. Such a partially relaxed state appears to be the NFW state, and is to be expected especially in young systems. The NFW core is shown to be isolated. In non-isolated systems continuing coarse-grained relaxation should be towards a density core in solid body rotation.

We study analytically the coarse and fine-grained distribution function established by the self-similar infall of collisionless matter. We find this function explicitly for isotropic and spherically symmetric systems in terms of cosmological initial conditions. The coarse-grained function is structureless and steady but the familiar phase space sheet sub-structure is recovered in the fine-grained limit. By breaking the self-similarity of the halo infall we are able to argue for a central density flattening. In addition there will be an edge steepening. The best fitting analytic density function is likely to be provided by a high order polytrope fit smoothly to an outer power law of index -3 for isolated systems. There may be a transition to a -4 power law in the outer regions of tidally truncated systems. We find a progressive central flattening that is expected to end either in the non-singular isothermal sphere, or in non-singular metastable polytropic cores. Therefore a collisionless system may pass through a family of polytropes of increasing order, finally approaching the limit of the non-singular isothermal sphere, if the `violent' collective relaxation is frequently re-excited by `merger' events. Our results suggest that no physics beyond that of the moderate collective relaxation often known as `violent relaxation' (due we think to wave-particle scattering) is necessary to explain the nature of dark matter density profiles.

The self-similar infall model (SSIM) is normally discussed in the context of radial orbits in spherical symmetry. However it is possible to retain the spherical symmetry while permitting the particles to move in Keplerian ellipses, each having the squared angular momentum peculiar to their 'shell'. The spherical 'shell', defined for example by the particles turning at a given radius, then moves according to the radial equation of motion of a 'shell' particle. The 'shell' itself has no physical existence except as an ensemble of particles, but it is convenient to sometimes refer to the shells since it is they that are followed by a shell code. In this note we find the distribution of squared angular momentum as a function of radius that yields the NFW density profile for the final dark matter halo. It transpires that this distribution is amply motivated dimensionally. An effective 'lambda' spin parameter is roughly constant over the shells. We also study the effects of angular momentum on the relaxation of a dark matter system using a three dimensional representation of the relaxed phase space.

In this paper, we present our conclusions from the numerical study of the collapse of a destabilized collisionless stellar system. We use both direct integration of the Vlasov-Poisson equations and an N-body tree code to obtain our results, which are mutually confirmed. We find that spherical and moderately nonspherical collapse configurations evolve to new equilibrium configurations in which the velocity distribution approaches a Gaussian form, at least in the central regions. The evolution to this state has long been an open question, and in this work we are able to clarify the process responsible and to support predictions made from statistical considerations (Lynden-Bell 1967; Nakamura 2000). The simulations of merging N-body systems show a transition to a Gaussian velocity distribution that is increasingly suppressed as the initial separation of centres is increased. Possible reasons for this are discussed.

This poster discusses a possible explanation for the relationship between the mass of the central supermassive black hole and the velocity dispersion in the bulge of the host galaxy. We suppose that the black hole and the dark matter halo are forming simultaneously as matter falls in, and a self-similar system then exists in which the mass and the velocities of the system evolve as power-law functions of time. This leads naturally to a relationship between the black hole mass and the velocities in the halo which, with a reasonable choice of cosmological parameters, is in good agreement with the observed relationship. We also confirm this relationship with more robust numerical results.

Powerful, highly collimated jets, surrounded by bipolar molecular outflows, are commonly observed near Young Stellar Objects (YSOs). In the usual theoretical picture of star formation, a jet is ejected from a magnetized accretion disk, with a molecular outflow being driven either by the jet or by a wider wind coming from the disk. Here, we propose an alternative global model for the flows surrounding YSOs. In addition to a central accretion-ejection engine driving the jet, the molecular outflow is powered by the infalling matter and follows a circulation pattern around the central object without necessarily being entrained by a jet. It is shown that the model produces a heated pressure-driven outflow with magneto-centrifugal acceleration and collimation. We report solutions for the three different parts of this self-similar model, i.e. the jet, the infalling envelope and the circulating matter that eventually forms the molecular outflow. This new picture of the accretion/outflow phase provides a possible explanation for several observed properties of YSO outflows. The most relevant ones are the presence of high mass molecular outflows around massive protostars, and a realistic fraction (typically 0.1) of the accretion flow that goes into the jet.

In this article we consider the growth of seed black holes immersed in dark matter halos. We first investigate the adiabatic growth in various initial distribution functions (isothermal, power law, and NFW) and find the resulting density, radial velocity, and anisotropy profiles. In addition we estimate the growth rate for a given black hole mass in the corresponding adiabatically modified dark matter distribution function. Only in the isothermal case is there a convincing black hole mass-age relation. By calculating the line of sight velocity dispersion for the various cases as a function of the black hole mass, we find the predicted adiabatic $M_{bh}-\sigma$ relation; this never approaches the recently observed power law. We conclude by abandoning adiabaticity, suggesting that the black hole grows proportionally to the dark matter halo itself on a dynamic time scale. This allows us to relate the observed $M_{bh}-\sigma$ relation to the cosmological power spectrum on galactic scales by using dimensional scaling arguments.

The density profiles established in the self-similar dynamical phase of dark matter haloes recollapse are all close to isothermal. This is steeper than the predictions of some n-body simulations for the central regions of the halo, which are in turn steeper than the profiles observed in the centre of some galaxies; particularly dwarfs and low surface brightness galaxies. The outer regions of galaxies both as observed and as simulated have density profiles steeper than the self-similar profile. Nevertheless there appears to be an intermediate region in most galaxies in which the inverse square behaviour is a good description. The outer deviations can be explained plausibly in terms of the transition from a self-gravitating extended halo to a Keplerian flow onto a dominant central mass (the isothermal distribution can not be complete), but the inner deviations are more problematic. Rather than attack this question directly, we use in this paper a novel coarse-graining technique combined with a shell code to establish both the distribution function associated with the self-similar density profile and the nature of the possible deviations in the central regions. In spherical symmetry we find that both in the case of purely radial orbits and in the case of orbits with non-zero angular momentum the self-similar density profile should flatten progressively near the centre of the system. The NFW limit of -1 seems possible. In a section aimed at demonstrating our technique for a spherically symmetric steady state, we argue that a Gaussian distribution function is the best approximation near the centre of the system.

The formation and collapse of a protostar involves the simultaneous infall and outflow of material in the presence of magnetic fields, self-gravity, and rotation. We use self-similar techniques to self-consistently model the anisotropic collapse and outflow by a set of angle-separated self-similar equations. The outflow is quite strong in our model, with the velocity increasing in proportion to radius, and material formally escaping to infinity in the finite time required for the central singularity to develop. Analytically tractable collapse models have been limited mainly to spherically symmetric collapse, with neither magnetic field nor rotation. Other analyses usually employ extensive numerical simulations, or either perturbative or quasistatic techniques. Our model is unique as an exact solution to the non-stationary equations of self-gravitating MHD, which features co-existing regions of infall and outflow. The velocity and magnetic topology of our model is quadrupolar, although dipolar solutions may also exist. We provide a qualitative model for the origin and subsequent evolution of such a state. However, a central singularity forms at late times, and we expect the late time behaviour to be dominated by the singularity rather than to depend on the details of its initial state. Our solution may, therefore, have the character of an attractor among a much more general class of self-similarity.

We study a self-similar circulation model for protostellar bipolar outflows. The model is axisymmetric and stationary, and now includes Poynting flux. Compared to an earlier version of the model, this addition produces faster and more collimated outflows. Moreover the luminosity needed for the radiative heating is smaller. The model has been applied to both low and high mass protostars. Monte Carlo methods have been used to explore systematically the parameter space. An inflow/outflow pattern including collimation of high speed material and an infalling toroidal disc arises naturally. The disc shape depends on the imposed heating, but it is naturally Keplerian given the central point mass. Outflows can have large opening angles, that increase when magnetic field weakens. Massive protostars produce faster but less collimated outflows than less massive protostars. The model is now at a stage where synthetic CO spectra reproduce very well the observational features. The results strengthen the idea that radiative heating and Poynting flux are ultimately the energy sources driving the outflow.

Navarro, Frenk, and White have suggested that the density profiles of simulated dark matter halos have a ``universal'' shape so that a given halo can be characterized by a single free parameter which fixes its mass. In this paper, we revisit the spherical infall model in the hope of recognizing in detail the existence and origin of any such universality. A system of particles is followed from linear perturbation, through first shell crossing, then through an accretion or infall phase, and finally to virialization. During the accretion phase, the system relaxes through a combination of phase mixing, phase space instability, and moderate violent relation. It is driven quickly, by the flow of mass through its surface, toward self-similar evolution. The self-similar solution plays its usual role of intermediate attractor and can be recognized from a virial-type theorem in scaled variables and from our numerical simulations. The transition to final equilibrium state once infall has ceased is relatively gentle, an observation which leads to an approximate form for the distribution function of the final system. The infall phase fixes the density profile in intermediate regions of the halo to be close to r^{-2}. We make contact with the standard hierarchical clustering scenario and explain how modifications of the self-similar infall model might lead to density profiles in agreement with those found in numerical simulations.

Recent observations suggest that mass ejection and mass accretion both decline significantly with time during early protostellar evolution (Bontemps et al. 1996). In the present paper, we propose that this rapid decay of accretion/ejection activity is a direct result of the non-singular density profiles characterizing pre-collapse clouds. The radial density profiles of pre-stellar cores are indeed found to be much flatter than rho(r) ~ r(-2) at radii less than a few thousand AU (Ward-Thompson et al. 1994). Here we show, through Lagrangian analytical calculations, that the supersonic gravitational collapse of pre-stellar cloud cores with centrally peaked, but flattened density profiles leads to a transitory phase of energetic accretion immediately following the formation of the central hydrostatic protostar. Physically, the flat inner region collapses first nearly homologously to form a finite mass stellar nucleus, and the remaining cloud core material then accretes supersonically onto a non-zero point mass. Enhanced accretion persists as long as the gravitational pull of this initial point mass, which does not exist in the Shu singular solution, remains significant. We suggest that this epoch of vigorous accretion coincides with Class 0 protostars, which would explain their unusually powerful jets compared to the more evolved Class I objects. We also use a simple two-component power-law model to fit the diagrams of outflow power versus envelope mass observed by Bontemps et al. (1996), and suggest that rho Ophiuchi and Taurus young stellar objects follow different accretion histories because of differing initial conditions.

We consider the evolution of an initially cold, spherically symmetric system of self-gravitating particles. At early times the system undergoes self-similar collapse of the type described by Fillmore & Goldreich and Bertschinger. This stage of phase mixing soon gives way to a period of violent relaxation driven by an instability in the similarity solution. The onset of violent relaxation is illustrated by numerical simulations and supported by non-linear analysis of the ``scaled'' collisionless Boltzmann and Poisson equations.

A new model for the high-energy emission from pulsars is developed by considering charged particle motion in the fields of a spinning, highly magnetised and conducting sphere in vacuum. A generally applicable approximation to the particle motion in strong fields is developed and applied to the numerical modelling, and the radiation emitted by curvature emission is summed to generate light curves. The model predicts many of the observed features of pulsar light curves. This paper outlines the basic properties of the model; a subsequent paper will discuss the statistical properties of a population of model pulsars and apply the model to the known gamma-ray pulsars.

In this article we give an account of a method of smoothing spatial inhomogeneous data sets by using wavelet reconstruction on a regular grid in an auxilliary space onto which the original data is mapped. In a previous paper by the present authors, we devised a method for inferring the velocity potential from the radial component of the cosmic velocity field assuming an ideal sampling. Unfortunately the sparseness of the real data as well as errors of measurement require us to first smooth the velocity field as observed on a 3-dimensional support (i.e. the galaxy positions) inhomogeneously distributed throughout the sampled volume. The wavelet formalism permits us to introduce a minimal smoothing procedure that is characterized by the variation in size of the smothing window function. Moreover the output smoothed radial velocity field can be shown to correspond to a well defined theoretical quantity as long as the spatial sampling support satisfies certain criteria. We argue also that one should be very cautious when comparing the velocity potential derived from such a smoothed radial component of the velocity field with related quantities derived from other studies (e.g : of the density field).

The motions of comets and neutron stars have been integrated over five billion years in the Galactic potential to determine a gamma-ray burst distribution, presuming that bursts are the result of interactions between these two families of objects. The comets originate in two distinct populations - one from ejection by stars in the Galactic disk, and the other from ejection by stars in globular clusters. No choice of the free parameters resulted in agreement with both the isotropy data and the $\log(N>F) - \log(F)$ data.

We study systematically stationary solutions to the coupled Vlasov and Poisson equations which have `self-similar' or scaling symmetry in phase space. In particular, we find analytically {\it all} spherically symmetric distribution functions where the mass density and gravitational potential are strict power laws in $r$, the distance from the symmetry point. We treat as special cases, systems built from purely radial orbits and systems that are isotropic in velocity space. We then discuss systems with arbitrary velocity space anisotropy finding a new and very general class of distribution functions. These distributions may prove useful in modelling galaxies. Distribution functions in cylindrical and planar geometries are also discussed. Finally, we study spatially spheroidal systems that again exhibit strict power-law behaviour for the density and potential and find results in agreement with results published recently.

We propose that the recently reported MACHO/EROS events correspond to gravitational amplification by dark clouds rather than compact objects. These clouds must be very dense with $M\sim 0.1 ~M_\odot$ and $R\la 10^{14}~{\rm cm}$. In all likelihood, the clouds will be members of a family of objects with different sizes and masses. We therefore expect events of longer duration than the ones reported by the MACHO and EROS groups but with light curves very different from the ones derived assuming point mass lenses. We suggest that one such event has already been observed in radio measurements of the quasar 1502+106. The abundances of free electrons, metals, complex molecules, and dust grains are constrained to be very small suggesting that the clouds are formed from a primordial mixture of hydrogen and helium. Cosmic rays and background UV radiation ionize a halo around the cloud. Radio waves from distant sources will be scattered by the electrons in this halo, an effect which may have already been observed in quasars such as 1502+106. We argue that dark clouds are a viable alternative to compact objects for baryonic dark matter in the halo.

We propose that the prevalence of bipolarity in Young Stellar Objects is due to the fine tuning that is required for spherical accretion of an ambient medium onto a central node.It is shown that there are two steady modes that are more likely than radial accretion, each of which is associated with a hyperbolic central point in the meridional stream lines, and consequently with either an equatorial inflow and an axial ejection or vice versa. In each case the stream lines pass through a thick accretion torus, which is better thought of as a standing pressure wave rather than as a relatively inert Keplerian structure.We base our arguments on a simple analytic example,which is topologi cally generic,wherein each bipolarmode is created by the rebound of accreting matter under the action of the thermal,magnetic,turbulent and centrifugal pressures created in the flow. In both bipolar modes the presence of non-zero angular momentum implies axial regions wherein the pressure is first reduced below the value at infinity and then becomes negative, where the solution fails because rotating material can not enter this region without suction.The model thus has empty stems where the activity of the central source must dominate.So the basic engine of the bipolar flow discussed here is simply the rebound of freely falling material from a thick pressure disc into an axial low pressure region.The low mass,high velocity outflow must be produced in this region by an additional mechanism. This is reminiscent of the cored apple structure observed recently in the very young bipolar source VLA 1623.