We solve the equations of motion of a one-dimensional planar Heisenberg (or Vaks-Larkin) model consisting of a system of interacting macro-spins aligned along a ring. Each spin has unit length and is described by its angle with respect to the rotational axis. The orientation of the spins can vary in time due to random forcing and spin-spin interaction. We statistically describe the behaviour of the sum of all spins for different parameters. The term "domino model" in the title refers to the interaction among the spins. We compare the model results with geomagnetic field reversals and find strikingly similar behaviour. The aggregate of all spins keeps the same direction for a long time and, once in a while, begins flipping to change the orientation by almost 180 degrees (mimicking a geomagnetic reversal) or to move back to the original direction (mimicking an excursion). Most of the time the spins are aligned or anti-aligned and deviate only slightly with respect to the rotational axis (mimicking the secular variation of the geomagnetic pole with respect to the geographic pole). Reversals are fast compared to the times in between and they occur at random times, both in the model and in the case of the Earth's magnetic field.

Geomagnetism is characterized by intermittent polarity reversals and rapid fluctuations. We have recently proposed a coupled macro-spin model to describe these dynamics based on the idea that the whole dynamo mechanism is described by the coherent interactions of many small dynamo elements. In this paper, we further develop this idea and construct a minimal model for magnetic variations. This simple model naturally yields many of the observed features of geomagnetism: its time evolution, the power spectrum, the frequency distribution of stable polarity periods, etc. This model has coexistent two phases; i.e. the cluster phase which determines the global dipole magnetic moment and the expanded phase which gives random perpetual perturbations that yield intermittent polarity flip of the dipole moment. This model can also describe the synchronization of the spin oscillation. This corresponds to the case of sun and the model well describes the quasi-regular cycles of the solar magnetism. Furthermore, by analyzing the relevant terms of MHD equation based on our model, we have obtained a scaling relation for the magnetism for planets, satellites, sun, and stars. Comparing it with various observations, we can estimate the scale of the macro-spins.

For the initial fields of the density contrast and peculiar velocity, we theoretically calculate the differential and integral length scales, i.e., statistical measures that respectively characterize the small- and large-scale fluctuations of a random field. These length scales and the associated mass scales explain the length and mass scales observed for (1) halos of young galaxies at z > 5, (2) halos of galaxies at z = 0, and (3) the largest structures in the galaxy distribution at z = 0. We thereby discuss that such observed scales are fossil imprints of the characteristic scales of the initial fields.

Using the initial peculiar velocity field, we analytically study the hierarchical formation of gravitationally bound objects. The field is smoothed over a scale that corresponds to the mass of a given class of objects. Through the Zel'dovich approximation, the smoothed field determines how the objects cluster together to form a new class of more massive objects. The standard cosmological parameters lead to the evolution of primordial clouds with 10^6 M(sun) -> galaxies with 10^12 M(sun) -> clusters of galaxies with 10^15 M(sun) -> superclusters of galaxies with 10^16 M(sun). The epochs obtained for the formation of these classes of objects are consistent with observations.

There is observed a trend that a lower mass galaxy forms stars at a later epoch. This downsizing of star-forming galaxies has been attributed to hydrodynamical or radiative feedback processes that regulate star formation. However, here we explain the downsizing by gravitational processes alone, in the bottom-up scenario where galaxies evolve from subgalactic-scale objects. Within a region of the initial density field that is to evolve into a lower mass galaxy, subgalactic-scale fluctuation is of a smaller amplitude. The formation of subgalactic-scale objects, i.e., gravitational collapse of the subgalactic-scale fluctuation, and the subsequent onset of star formation accordingly occur at a later epoch for a lower mass galaxy. As a function of galaxy mass, we calculate the peak epoch of formation of subgalactic-scale objects. The peak epoch is consistent with the peak epoch of star formation derived from observations. [abridged]

Why does the Zel'dovich approximation (ZA) work well to describe gravitational collapse in the universe? This problem is examined by focusing on its dependence on the dimensionality of the collapse. The ZA is known to be exact for a one-dimensional collapse. We show that the ZA becomes progressively more accurate in the order of three-, two-, and one-dimensional collapses. Furthermore, using models for spheroidal collapse, we show that the ZA remains accurate in all collapses, which become progressively lower dimensional with the passage of time. That is, the ZA is accurate because the essence of the gravitational collapse is incorporated in the ZA.

Galaxy interactions are studied in terms of the starburst-Seyfert connection. The starburst requires a high rate of gas supply. Since the efficiency for supplying the gas is high in a galaxy interaction, although the companion is not necessarily discernible, Seyfert galaxies with circumnuclear starbursts are expected to be interacting. Since the large amounts of circumnuclear gas and dust obscure the broad-line region, they are expected to be observed as Seyfert 2. The active galactic nucleus itself does not require a high rate of gas supply. Seyfert galaxies without circumnuclear starbursts are not necessarily expected to be interacting even at the highest luminosities. They are not necessarily expected to evolve from Seyfert galaxies with circumnuclear starbursts. We derive these and other theoretical expectations and confirm them with statistics on observational data of magnitude-limited samples of Seyfert galaxies.

According to N-body numerical simulations, the radial pairwise velocities of galaxies have negative skewness in the quasi-nonlinear regime. To understand its origin, we calculate the probability distribution function of the radial pairwise velocity using the Zel'dovich approximation, i.e., an analytical approximation for gravitational clustering. The calculated probability distribution function is in good agreement with the result of N-body simulations. Thus the negative skewness originates in relative motions of galaxies in the clustering process that the infall dominates over the expansion.

The orbital decay and tidal disruption of a star cluster in a galaxy is studied in an analytical manner. Owing to dynamical friction, the star cluster spirals in toward the center of the galaxy. Simultaneously, the galactic tidal field strips stars from the outskirts of the star cluster. Under an assumption that the star cluster undergoes a self-similar evolution, we obtain the condition and timescale for the star cluster to reach the galaxy center before its disruption. The result is used to discuss the fate of so-called intermediate-mass black holes with >10^3 M(sun) found recently in young star clusters of starburst galaxies and also the mass function of globular clusters in galaxies.

Mass segregation in a star cluster is studied in an analytical manner. We consider a two-component cluster, which consists of two types of stars with different masses. Plummer's model is used for the initial condition. We trace the overall behaviors of the probability distribution functions of the two components and obtain the timescale of mass segregation as a simple function of the cluster parameters. The result is used to discuss the origin of a black hole with mass of > 1000 M(sun) found in the starburst galaxy M82.

Following the discovery of a black hole (BH) with a mass of 10^3-10^6 M(sun) in a starburst galaxy M82, we study formation of such a BH via successive merging of stellar-mass BHs within a star cluster. The merging has a runaway characteristic. This is because massive BHs sink into the cluster core and have a high number density, and because the merging probability is higher for more massive BHs. We use the Smoluchowski equation to study analytically the evolution of the BH mass distribution. Under favorable conditions, which are expected for some star clusters in starburst galaxies, the timescale of the runaway merging is at most of order 10^7 yr. This is short enough to account for the presence of a BH heavier than 10^3 M(sun) in an ongoing starburst region.

Far-infrared flux densities are newly extracted from the IRAS database for the RSA and CfA complete samples of Seyfert galaxies. These data are used to classify the Seyfert galaxies into those where the far-infrared continuum emission is dominated by the active galactic nucleus (AGN), circumnuclear starburst, or host galaxy. While AGN-dominant objects consist of comparable numbers of Seyfert 1 and 2 galaxies, starburst- and host-dominant objects consist preferentially of Seyfert 2 galaxies. Thus, in addition to the dusty torus, the circumnuclear starburst region and host galaxy are important in hiding the broad-line region. Morphologically, starburst-dominant Seyfert galaxies are of later types and more strongly interacting than AGN-dominant Seyfert galaxies. In a later-type galaxy, the AGN central engine has a lower Eddington luminosity, and the gaseous content is higher. The gas is efficiently supplied to the starburst via a galaxy-galaxy interaction. Morphologies of host-dominant Seyfert galaxies are of various types. Since starbursts in Seyfert galaxies are older than those in classical starburst galaxies, we propose an evolution from starburst to starburst-dominant Seyfert to host-dominant Seyfert for a late-type galaxy. An evolution from AGN-dominant Seyfert to host-dominant Seyfert is proposed for an early-type galaxy. These sequences have durations of a few x 10^8 yr and occur repeatedly within a galaxy during its evolution from a late type to an early type.

Old starburst galaxies are deficient in O stars, and hence do not exhibit strong line emission in the optical regime. However, there remain many B stars, which are expected to heat dust grains and generate strong continuum emission in the far-infrared. The IRAS data for a statistically complete sample of nearby galaxies reveal for the first time that such objects are as numerous as 30-40% of the local galaxy population.

By using the flux ratio [FeII]8617/[OI]6300, we demonstrate that most of the interstellar dust grains survive in shocks associated with supernova remnants and Herbig-Haro objects. The [FeII]/[OI] flux ratio is sensitive to the gas-phase Fe/O abundance ratio, but is insensitive to the ionization state, temperature, and density of the gas. We calculate the [FeII]/[OI] flux ratio in shocks, and compare the results with the observational data. When only 20% of iron is in the gas phase, the models reproduce most successfully the observations. This finding is in conflict with the current consensus that shocks destroy almost all the grains and 100% of metals are in the gas phase. We comment on previous works on grain destruction, and discuss why grains are not destroyed in shocks.

In the paper "Excitation Mechanism of Near-Infrared [Fe II] Emission in Seyfert and Starburst Galaxies" by Hideaki Mouri, Kimiaki Kawara, and Yoshiaki Taniguchi (ApJ, 528, 186 [2000]), the two panels of Figure 6 were printed into different sizes, as the result of an error in the printing process. The correct version of Figure 6 appears here.

Dusty tori around active galactic nuclei (AGNs) play an important role in the classification of Seyfert galaxies. Therefore, physical properties of dusty tori are of great interest. We briefly introduce three statistical studies investigating properties of dusty tori; 1) physical sizes of dusty tori based on water-vapor maser emission, 2) ionization condition of the inner wall of tori based on high-ionization emission lines, and 3) viewing angle toward dusty tori based on mid-infrared color.

We propose a new diagnostic of the ``dusty torus'' model for Seyfert nuclei. Dust grains in the torus are heated by the nuclear continuum, and reradiate mostly in the mid-infrared wavelengths. From the torus geometry, it is predicted that the emission at lambda < 10 micron has strong dependence on the viewing angle. Since the dependence is predicted to be insignificant at lambda > 10 micron, we study the flux ratio between 3.5 micron (L band) and 25 micron; R(L,25) = log [(nu_3.5 um S_nu_3.5 um)/(nu_25 um S_nu_25 um)]. In three different samples (optically selected, X-ray selected, and infrared selected samples) of Seyfert galaxies, the observed values of R(L,25) between type 1 Seyferts (S1s) and type 2 Seyferts (S2s) are found to be clearly separated; R(L,25) > -0.6 for S1s while R(L,25) < -0.6 for S2s. This implies universality of their torus properties. With this result and the other observational characteristics, we investigate the most plausible torus model among those presented in Pier & Krolik (1992, 1993).