Oscillatory instability in a driven granular gas

We discovered an oscillatory instability in a system of inelastically colliding hard spheres, driven by two opposite "thermal" walls at zero gravity. The instability, predicted by a linear stability analysis of the equations of granular hydrodynamics, occurs when the inelasticity of particle collisions exceeds a critical value. Molecular dynamic simulations support the theory and show a stripe-shaped cluster moving back and forth in the middle of the box away from the driving walls. The oscillations are irregular but have a single dominating frequency that is close to the frequency at the instability onset, predicted from hydrodynamics.Comment: 7 pages, 4 figures, to appear in Europhysics Letter

Geometrical optics of large deviations of fractional Brownian motion

It has been shown recently that the optimal fluctuation method -- essentially geometrical optics -- provides a valuable insight into large deviations of Brownian motion. Here we extend the geometrical optics formalism to two-sided, $-\infty<t<\infty$, fractional Brownian motion (fBM) on the line, which is "pushed" to a large deviation regime by imposed constraints. We test the formalism on three examples where exact solutions are available: the two- and three-point probability distributions of the fBm and the distribution of the area under the fBm on a specified time interval. Then we apply the formalsim to several previously unsolved problems by evaluating large-deviation tails of the following distributions: (i) of the first-passage time, (ii) of the maximum of, and (iii) of the area under, fractional Brownian bridge and fractional Brownian excursion, and (iv) of the first-passage area distribution of the fBm. An intrinsic part of a geometrical optics calculation is determination of the optimal path -- the most likely realization of the process which dominates the probability distribution of the conditioned process. Due to the non-Markovian nature of the fBm, the optimal paths of a fBm, subject to constraints on a finite interval $0<t\leq T$, involve both the past $-\infty<t<0$ and the future $T<t<\infty$.Comment: 8 pages, 6 figure

Wavelength limits on isobaricity of perturbations in a thermally unstable radiatively cooling medium

Nonlinear evolution of one-dimensional planar perturbations in an optically thin radiatively cooling medium in the long-wavelength limit is studied numerically. The accepted cooling function generates in thermal equilibrium a bistable equation of state $P(\rho)$. The unperturbed state is taken close to the upper (low-density) unstable state with infinite compressibility ($dP/d\rho= 0$). The evolution is shown to proceed in three different stages. At first stage, pressure and density set in the equilibrium equation of state, and velocity profile steepens gradually as in case of pressure-free flows. At second stage, those regions of the flow where anomalous pressure (i.e. with negative compressibility) holds, create velocity profile more sharp than in pressure-free case, which in turn results in formation of a very narrow (short-wavelength) region where gas separates the equilibrium equation of state and pressure equilibrium sets in rapidly. On this stage, variation in pressure between narrow dense region and extended environment does not exceed more than 0.01 of the unperturbed value. On third stage, gas in the short-wavelength region reaches the second (high-density) stable state, and pressure balance establishes through the flow with pressure equal to the one in the unperturbed state. In external (long-wavelength) regions, gas forms slow isobaric inflow toward the short-wavelength layer. The duration of these stages decreases when the ratio of the acoustic time to the radiative cooling time increases. Limits in which nonlinear evolution of thermally unstable long-wavelength perturbations develops in isobaric regime are obtained.Comment: 21 pages with 7 figures, Revtex, accepted in Physics of Plasma

Logarithmically Slow Expansion of Hot Bubbles in Gases

We report logarithmically slow expansion of hot bubbles in gases in the process of cooling. A model problem first solved, when the temperature has compact support. Then temperature profile decaying exponentially at large distances is considered. The periphery of the bubble is shown to remain essentially static ("glassy") in the process of cooling until it is taken over by a logarithmically slowly expanding "core". An analytical solution to the problem is obtained by matched asymptotic expansion. This problem gives an example of how logarithmic corrections enter dynamic scaling.Comment: 4 pages, 1 figur

Reducing multiphoton ionization in a linearly polarized microwave field by local control

We present a control procedure to reduce the stochastic ionization of hydrogen atom in a strong microwave field by adding to the original Hamiltonian a comparatively small control term which might consist of an additional set of microwave fields. This modification restores select invariant tori in the dynamics and prevents ionization. We demonstrate the procedure on the one-dimensional model of microwave ionization.Comment: 8 page

Time-resolved extinction rates of stochastic populations

Extinction of a long-lived isolated stochastic population can be described as an exponentially slow decay of quasi-stationary probability distribution of the population size. We address extinction of a population in a two-population system in the case when the population turnover -- renewal and removal -- is much slower than all other processes. In this case there is a time scale separation in the system which enables one to introduce a short-time quasi-stationary extinction rate W_1 and a long-time quasi-stationary extinction rate W_2, and develop a time-dependent theory of the transition between the two rates. It is shown that W_1 and W_2 coincide with the extinction rates when the population turnover is absent, and present but very slow, respectively. The exponentially large disparity between the two rates reflects fragility of the extinction rate in the population dynamics without turnover.Comment: 8 pages, 4 figure

Navier-Stokes hydrodynamics of thermal collapse in a freely cooling granular gas

We employ Navier-Stokes granular hydrodynamics to investigate the long-time behavior of clustering instability in a freely cooling dilute granular gas in two dimensions. We find that, in circular containers, the homogeneous cooling state (HCS) of the gas loses its stability via a sub-critical pitchfork bifurcation. There are no time-independent solutions for the gas density in the supercritical region, and we present analytical and numerical evidence that the gas develops thermal collapse unarrested by heat diffusion. To get more insight, we switch to a simpler geometry of a narrow-sector-shaped container. Here the HCS loses its stability via a transcritical bifurcation. For some initial conditions a time-independent inhomogeneous density profile sets in, qualitatively similar to that previously found in a narrow-channel geometry. For other initial conditions, however, the dilute gas develops thermal collapse unarrested by heat diffusion. We determine the dynamic scalings of the flow close to collapse analytically and verify them in hydrodynamic simulations. The results of this work imply that, in dimension higher than one, Navier-Stokes hydrodynamics of a dilute granular gas is prone to finite-time density blowups. This provides a natural explanation to the formation of densely packed clusters of particles in a variety of initially dilute granular flows.Comment: 18 pages, 19 figure

Emergence of stability in a stochastically driven pendulum: beyond the Kapitsa effect

We consider a prototypical nonlinear system which can be stabilized by multiplicative noise: an underdamped non-linear pendulum with a stochastically vibrating pivot. A numerical solution of the pertinent Fokker-Planck equation shows that the upper equilibrium point of the pendulum can become stable even when the noise is white, and the "Kapitsa pendulum" effect is not at work. The stabilization occurs in a strong-noise regime where WKB approximation does not hold.Comment: 4 pages, 7 figure