Nature of vibrational eigenmodes in topologically disordered solids

We use a local projectional analysis method to investigate the effect of topological disorder on the vibrational dynamics in a model glass simulated by molecular dynamics. Evidence is presented that the vibrational eigenmodes in the glass are generically related to the corresponding eigenmodes of its crystalline counterpart via disorder-induced level-repelling and hybridization effects. It is argued that the effect of topological disorder in the glass on the dynamical matrix can be simulated by introducing positional disorder in a crystalline counterpart.Comment: 7 pages, 6 figures, PRB, to be publishe

Persistence in higher dimensions : a finite size scaling study

We show that the persistence probability $P(t,L)$, in a coarsening system of linear size $L$ at a time $t$, has the finite size scaling form $P(t,L)\sim L^{-z\theta}f(\frac{t}{L^{z}})$ where $\theta$ is the persistence exponent and $z$ is the coarsening exponent. The scaling function $f(x)\sim x^{-\theta}$ for $x \ll 1$ and is constant for large $x$. The scaling form implies a fractal distribution of persistent sites with power-law spatial correlations. We study the scaling numerically for Glauber-Ising model at dimension $d = 1$ to 4 and extend the study to the diffusion problem. Our finite size scaling ansatz is satisfied in all these cases providing a good estimate of the exponent $\theta$.Comment: 4 pages in RevTeX with 6 figures. To appear in Phys. Rev.

Persistence in a Stationary Time-series

We study the persistence in a class of continuous stochastic processes that are stationary only under integer shifts of time. We show that under certain conditions, the persistence of such a continuous process reduces to the persistence of a corresponding discrete sequence obtained from the measurement of the process only at integer times. We then construct a specific sequence for which the persistence can be computed even though the sequence is non-Markovian. We show that this may be considered as a limiting case of persistence in the diffusion process on a hierarchical lattice.Comment: 8 pages revte

Persistence of a Continuous Stochastic Process with Discrete-Time Sampling: Non-Markov Processes

We consider the problem of `discrete-time persistence', which deals with the zero-crossings of a continuous stochastic process, X(T), measured at discrete times, T = n(\Delta T). For a Gaussian Stationary Process the persistence (no crossing) probability decays as exp(-\theta_D T) = [\rho(a)]^n for large n, where a = \exp[-(\Delta T)/2], and the discrete persistence exponent, \theta_D, is given by \theta_D = \ln(\rho)/2\ln(a). Using the `Independent Interval Approximation', we show how \theta_D varies with (\Delta T) for small (\Delta T) and conclude that experimental measurements of persistence for smooth processes, such as diffusion, are less sensitive to the effects of discrete sampling than measurements of a randomly accelerated particle or random walker. We extend the matrix method developed by us previously [Phys. Rev. E 64, 015151(R) (2001)] to determine \rho(a) for a two-dimensional random walk and the one-dimensional random acceleration problem. We also consider `alternating persistence', which corresponds to a < 0, and calculate \rho(a) for this case.Comment: 14 pages plus 8 figure

Economics of controlled urea fertilizer in cropping systems

Non-Peer ReviewedNitrogen fertilizers are chemical compounds given to plants to promote growth. However, their application can impact the environment through ammonia and greenhouse gas emissions, and through leaching which can lead to deterioration of water quality. To understand the nitrogen dynamics in the soil, a set of experiments were designed by Agriculture and Agri-Food Canada (AAFC) research scientists with funding received from the Environmental Technology Assessment for Agriculture (ETAA) and GAPS programs of AAFC, and Agrium, a fertilizer company. Research trials were conducted to assess the impact of urea with comparison being made between ESN (environmentally sensitive nitrogen) and uncoated urea. Trials were conducted in five provinces and eight research sites. Harrow was the only eastern Canada site. Experiments also included several fertilizer application rates, and at least two tillage systems – conventional and reduced tillage. Zone tillage was also included at Harrow, Ontario. The general conclusion from the foregoing set of data and analyses is that ESN application is a better economic choice for certain crops in some regions, but not in all regions or for all crops. Even a single crop was not found to generate a positive producer surplus in all regions. Reasons for these differences need further investigation

Random Walks in Logarithmic and Power-Law Potentials, Nonuniversal Persistence, and Vortex Dynamics in the Two-Dimensional XY Model

The Langevin equation for a particle (`random walker') moving in d-dimensional space under an attractive central force, and driven by a Gaussian white noise, is considered for the case of a power-law force, F(r) = - Ar^{-sigma}. The `persistence probability', P_0(t), that the particle has not visited the origin up to time t, is calculated. For sigma > 1, the force is asymptotically irrelevant (with respect to the noise), and the asymptotics of P_0(t) are those of a free random walker. For sigma < 1, the noise is (dangerously) irrelevant and the asymptotics of P_0(t) can be extracted from a weak noise limit within a path-integral formalism. For the case sigma=1, corresponding to a logarithmic potential, the noise is exactly marginal. In this case, P_0(t) decays as a power-law, P_0(t) \sim t^{-theta}, with an exponent theta that depends continuously on the ratio of the strength of the potential to the strength of the noise. This case, with d=2, is relevant to the annihilation dynamics of a vortex-antivortex pair in the two-dimensional XY model. Although the noise is multiplicative in the latter case, the relevant Langevin equation can be transformed to the standard form discussed in the first part of the paper. The mean annihilation time for a pair initially separated by r is given by t(r) \sim r^2 ln(r/a) where a is a microscopic cut-off (the vortex core size). Implications for the nonequilibrium critical dynamics of the system are discussed and compared to numerical simulation results.Comment: 10 pages, 1 figur

Large-Deviation Functions for Nonlinear Functionals of a Gaussian Stationary Markov Process

We introduce a general method, based on a mapping onto quantum mechanics, for investigating the large-T limit of the distribution P(r,T) of the nonlinear functional r[V] = (1/T)\int_0^T dT' V[X(T')], where V(X) is an arbitrary function of the stationary Gaussian Markov process X(T). For T tending to infinity at fixed r we find that P(r,T) behaves as exp[-theta(r) T], where theta(r) is a large deviation function. We present explicit results for a number of special cases, including the case V(X) = X \theta(X) which is related to the cooling and the heating degree days relevant to weather derivatives.Comment: 8 page

Survival-Time Distribution for Inelastic Collapse

In a recent publication [PRL {\bf 81}, 1142 (1998)] it was argued that a randomly forced particle which collides inelastically with a boundary can undergo inelastic collapse and come to rest in a finite time. Here we discuss the survival probability for the inelastic collapse transition. It is found that the collapse-time distribution behaves asymptotically as a power-law in time, and that the exponent governing this decay is non-universal. An approximate calculation of the collapse-time exponent confirms this behaviour and shows how inelastic collapse can be viewed as a generalised persistence phenomenon.Comment: 4 pages, RevTe

Fraction of uninfected walkers in the one-dimensional Potts model

The dynamics of the one-dimensional q-state Potts model, in the zero temperature limit, can be formulated through the motion of random walkers which either annihilate (A + A -> 0) or coalesce (A + A -> A) with a q-dependent probability. We consider all of the walkers in this model to be mutually infectious. Whenever two walkers meet, they experience mutual contamination. Walkers which avoid an encounter with another random walker up to time t remain uninfected. The fraction of uninfected walkers is investigated numerically and found to decay algebraically, U(t) \sim t^{-\phi(q)}, with a nontrivial exponent \phi(q). Our study is extended to include the coupled diffusion-limited reaction A+A -> B, B+B -> A in one dimension with equal initial densities of A and B particles. We find that the density of walkers decays in this model as \rho(t) \sim t^{-1/2}. The fraction of sites unvisited by either an A or a B particle is found to obey a power law, P(t) \sim t^{-\theta} with \theta \simeq 1.33. We discuss these exponents within the context of the q-state Potts model and present numerical evidence that the fraction of walkers which remain uninfected decays as U(t) \sim t^{-\phi}, where \phi \simeq 1.13 when infection occurs between like particles only, and \phi \simeq 1.93 when we also include cross-species contamination.Comment: Expanded introduction with more discussion of related wor