Coagulation of submicron-sized dust grains into porous aggregates is the initial step of dust evolution in protoplanetary disks. Recently, it has been pointed out that negative charging of dust in the weakly ionized disks could significantly slow down the coagulation process. In this paper, we apply the growth criteria obtained in Paper I to finding out a location ("frozen" zone) where the charging stalls dust growth at the fractal growth stage. For low-turbulence disks, we find that the frozen zone can cover the major part of the disks at a few to 100 AU from the central star. The maximum mass of the aggregates is approximately 10^{-7} g at 1 AU and as small as a few monomer masses at 100 AU. Strong turbulence can significantly reduce the size of the frozen zone, but such turbulence will cause the fragmentation of macroscopic aggregates at later stages. We examine a possibility that complete freezeout of dust evolution in low-turbulence disks could be prevented by global transport of dust in the disks. Our simple estimation shows that global dust transport can lead to the supply of macroscopic aggregates and the removal of frozen aggregates on a timescale of 10^6 yr. This overturns the usual understanding that tiny dust particles get depleted on much shorter timescales unless collisional fragmentation is effective. The frozen zone together with global dust transport might explain "slow" (\sim 10^6 yr) dust evolution suggested by infrared observation of T Tauri stars and by radioactive dating of chondrites.

Collisional growth of submicron-sized dust grains into macroscopic aggregates is the first step of planet formation in protoplanetary disks. These grains are expected to carry nonzero negative charges in the weakly ionized disks, but its effect on their collisional growth has not been fully understood so far. In this paper, we investigate how the charging affects the evolution of the dust size distribution properly taking into account the charging mechanism in a weakly ionized gas as well as porosity evolution through low-energy collisions. To clarify the role of the size distribution, we divide our analysis into two steps. First, we analyze the collisional growth of charged aggregates assuming a monodisperse (i.e., narrow) size distribution. We show that the monodisperse growth stalls due to the electrostatic repulsion when a certain condition is met, as is already expected in the previous work. Second, we numerically simulate dust coagulation using Smoluchowski's method to see how the outcome changes when the size distribution is allowed to freely evolve. We find that, under certain conditions, the dust undergoes bimodal growth where only a limited number of aggregates continue to grow carrying the major part of the dust mass in the system. This occurs because remaining small aggregates efficiently sweep up free electrons to prevent the larger aggregates from being strongly charged. We obtain a set of simple criteria that allows us to predict how the size distribution evolves for a given condition. In Paper II (arXiv:1009.3101), we apply these criteria to dust growth in protoplanetary disks.

Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. We point out an important implication of this discovery for dust growth in protoplanetary disks.

In a general metric theory of gravitation in four dimensions, six polarizations of a gravitational wave are allowed: two scalar and two vector modes, in addition to two tensor modes in general relativity. Such additional polarization modes appear due to additional degrees of freedom in modified theories of gravitation or theories with extra dimensions. Thus, observations of gravitational waves can be utilized to constrain the extended models of gravitation. In this paper, we investigate detectability of additional polarization modes of gravitational waves, particularly focusing on a stochastic gravitational-wave background, with laser-interferometric detectors on the Earth. We found that multiple detectors can separate the mixture of polarization modes in detector outputs, and that they have almost the same sensitivity to each polarization mode of stochastic gravitational-wave background.

We investigate the nonlinear evolution of black hole ringdown in the framework of higher-order metric perturbation theory. By solving the initial-value problem of a simplified nonlinear field model analytically as well as numerically, we find that (i) second-order quasinormal modes (QNMs) are indeed excited at frequencies different from those of first-order QNMs, as predicted recently. We also find serendipitously that (ii) late-time evolution is dominated by a new type of power-law tail. This ``second-order power-law tail'' decays more slowly than any late-time tails known in the first-order (i.e., linear) perturbation theory, and is generated at the wavefront of the first-order perturbation by an essentially nonlinear mechanism. These nonlinear components should be particularly significant for binary black hole coalescences, and could open a new precision science in gravitational wave studies.

Quasinormal ringing of acoustic black holes in Laval nozzles is discussed. The equation for sounds in a transonic flow is written into a Schr\"{o}dinger-type equation with a potential barrier, and the quasinormal frequencies are calculated semianalytically. From the results of numerical simulations, it is shown that the quasinormal modes are actually excited when the transonic flow is formed or slightly perturbed, as well as in the real black hole case. In an actual experiment, however, the purely-outgoing boundary condition will not be satisfied at late times due to the wave reflection at the end of the apparatus, and a late-time ringing will be expressed as a superposition of "boxed" quasinormal modes. It is shown that the late-time ringing damps more slowly than the ordinary quasinormal ringing, while its central frequency is not greatly different from that of the ordinary one. Using this fact, an efficient way for experimentally detecting the quasinormal ringing of an acoustic black hole is discussed.

In this paper, a quantitative characterization for the evolutionary sequence of stellar self-gravitating system is investigated, focusing on the pre-collapse stage of the long-term dynamical evolution. In particular, we consider the quasi-equilibrium behaviors of the N-body systems in the setup of the so-called Antonov problem, i.e., self-gravitating N-body system confined in an adiabatic wall and try to seek a possible connection with thermostatistics of self-gravitating systems. For this purpose, a series of long-term N-body simulations with various initial conditions are performed. We found that a quasi-equilibrium sequence away from the thermal equilibrium can be characterized by the one-parameter family of the stellar models. Especially, the stellar polytropic distribution satisfying the effective equation of state $P\propto\rho^{1+1/n}$ provides an excellent approximation to the evolutionary sequence of the N-body system. Based on the numerical results, we discuss a link between the quasi-equilibrium state and the generalized thermostatistics by means of the non-extensive entropy.

After introducing the fundamental properties of self-gravitating systems, we present an application of Tsallis' generalized entropy to the analysis of their thermodynamic nature. By extremizing the Tsallis entropy, we obtain an equation of state known as the stellar polytrope. For a self-gravitating stellar system confined within a perfectly reflecting wall, we discuss the thermodynamic instability caused by its negative specific heat. The role of the extremum as a quasi-equilibrium is also demonstrated from the results of N-body simulations.

With particular attention to the recently postulated introduction of a non-extensive generalization of Boltzmann-Gibbs statistics, we study the long-term stellar dynamical evolution of self-gravitating systems on timescales much longer than the two-body relaxation time. In a self-gravitating N-body system confined in an adiabatic wall, we show that the quasi-equilibrium sequence arising from the Tsallis entropy, so-called stellar polytropes, plays an important role in characterizing the transient states away from the Boltzmann-Gibbs equilibrium state.

We revisit the issues on the thermodynamic property of stellar self-gravitating system arising from Tsallis' non-extensive entropy. Previous papers (Taruya & Sakagami, Physica A 307 (2002) 185 (cond-mat/0107494); ibid. (2002) in press (cond-mat/0204315)) have revealed that the extremum-state of Tsallis entropy characterized by the so-called stellar polytrope has consistent thermodynamic structure, which predicts the thermodynamic instability due to the negative specific heat. However, their analyses heavily relies on the old Tsallis formalism using standard linear mean values. In this paper, extending our previous study, we focus on the equilibrium structure based on the standard framework by means of the normalized q-expectation values. It then turns out that the new extremum-state of Tsallis entropy essentially remains unchanged from the previous result, i.e., the stellar quasi-equilibrium distribution can be described by the stellar polytrope. While the thermodynamic stability for a system confined in an adiabatic wall completely agrees with the previous study and thereby the stability/instability criterion remains unchanged, the stability analysis reveals a new equilibrium property for the system surrounded by a thermal bath. In any case, the stability/instability criteria are consistently explained from the presence of negative specific heat and within the formalism, the stellar polytrope is characterized as a plausible non-extensive meta-equilibrium state.

We investigate the thermodynamic properties of stellar self-gravitating system arising from the Tsallis generalized entropy. In particular, physical interpretation of the thermodynamic instability, as has been revealed by previous paper(Taruya & Sakagami, cond-mat/0107494, Physica A 307, 185 (2002)), is discussed in detail based on the non-extensive thermostatistics. Examining the Clausius relation in a quasi-static experiment, we obtain the standard result of thermodynamic relation that the physical temperature of the equilibrium non-extensive system is identified with the inverse of the Lagrange multiplier, $T_{phys}=1/\beta$. Using this relation, the specific heat of total system is computed, and confirm the common feature of self-gravitating system that the presence of negative specific heat leads to the thermodynamic instability. In addition to the gravothermal instability discovered previously, the specific heat shows the curious divergent behavior at the polytrope index $n>3$, suggesting another type of thermodynamic instability. Evaluating the second variation of free energy, we find that the marginal stability condition indicated from the specific heat can be exactly recovered from the second variation of free energy. Thus, the stellar polytropic system is consistently characterized by the non-extensive thermostatistics as a plausible thermal equilibrium state. We also clarify the non-trivial scaling behavior appeared in specific heat and address the origin of non-extensive nature in stellar polytrope.

We propose an experimental model using the Laval nozzle of a sonic analogue of the Hawking radiation. We observe a power spectrum of the out-going wave emitted from the vicinity a sonic horizon in place of a created particle number. Our treatment is based on classical theory, and it will make experiments to be easier. This experimental possibility is a great advantage of our model.

We present a first physical application of Tsallis' generalized entropy to the thermodynamics of self-gravitating systems. The stellar system confined in a spherical cavity of radius $r_e$ exhibits an instability, so-called gravothermal catastrophe, which has been originally investigated by Antonov (1962) and Lynden-Bell & Wood (1968) on the basis of the maximum entropy principle for the phase-space distribution function. In contrast to previous analyses using the Boltzmann-Gibbs entropy, we apply the Tsallis-type generalized entropy to seek the equilibrium criteria. Then the distribution function of Vlassov-Poisson system can be reduced to the stellar polytrope system. Evaluating the second variation of Tsallis entropy and solving the zero eigenvalue problem explicitly, we find that the gravothermal instability appears in cases with polytrope index $n>5$. The critical point characterizing the onset of instability are obtained, which exactly matches with the results derived from the standard turning-point analysis. The results give an important suggestion that the Tsallis entropy is indeed applicable and viable to the long-range nature of the self-gravitating system.