Nontrivial Polydispersity Exponents in Aggregation Models

We consider the scaling solutions of Smoluchowski's equation of irreversible aggregation, for a non gelling collision kernel. The scaling mass distribution f(s) diverges as s^{-tau} when s->0. tau is non trivial and could, until now, only be computed by numerical simulations. We develop here new general methods to obtain exact bounds and good approximations of $\tau$. For the specific kernel KdD(x,y)=(x^{1/D}+y^{1/D})^d, describing a mean-field model of particles moving in d dimensions and aggregating with conservation of ``mass'' s=R^D (R is the particle radius), perturbative and nonperturbative expansions are derived. For a general kernel, we find exact inequalities for tau and develop a variational approximation which is used to carry out the first systematic study of tau(d,D) for KdD. The agreement is excellent both with the expansions we derived and with existing numerical values. Finally, we discuss a possible application to 2d decaying turbulence.Comment: 16 pages (multicol.sty), 6 eps figures (uses epsfig), Minor corrections. Notations improved, as published in Phys. Rev. E 55, 546

Kinetic Anomalies in Addition-Aggregation Processes

We investigate irreversible aggregation in which monomer-monomer, monomer-cluster, and cluster-cluster reactions occur with constant but distinct rates K_{MM}, K_{MC}, and K_{CC}, respectively. The dynamics crucially depends on the ratio gamma=K_{CC}/K_{MC} and secondarily on epsilon=K_{MM}/K_{MC}. For epsilon=0 and gamma<2, there is conventional scaling in the long-time limit, with a single mass scale that grows linearly in time. For gamma >= 2, there is unusual behavior in which the concentration of clusters of mass k, c_k decays as a stretched exponential in time within a boundary layer k<k* propto t^{1-2/gamma} (k* propto ln t for gamma=2), while c_k propto t^{-2} in the bulk region k>k*. When epsilon>0, analogous behaviors emerge for gamma<2 and gamma >= 2.Comment: 6 pages, 2 column revtex4 format, for submission to J. Phys.

Communities of Local Optima as Funnels in Fitness Landscapes

We conduct an analysis of local optima networks extracted from fitness landscapes of the Kauffman NK model under iterated local search. Applying the Markov Cluster Algorithm for community detection to the local optima networks, we find that the landscapes consist of multiple clusters. This result complements recent findings in the literature that landscapes often decompose into multiple funnels, which increases their difficulty for iterated local search. Our results suggest that the number of clusters as well as the size of the cluster in which the global optimum is located are correlated to the search difficulty of landscapes. We conclude that clusters found by community detection in local optima networks offer a new way to characterize the multi-funnel structure of fitness landscapes

Effects of in-chain and off-chain substitutions on spin fluctuations in the spin-Peierls compound CuGeO_3

The effect of in-chain and off-chain substitutions on 1D spin fluctuations in the spin-Peierls compound CuGeO_3 has been studied using Raman scattering in order to understand the interplay between defect induced states, enhanced spin-spin correlations and the ground state of low dimensional systems. In-chain and off-chain substitutions quench the spin-Peierls state and induce 3D antiferromagnetic order at T\leq 5 K. Consequently a suppression of a 1D gap-induced mode as well as a constant intensity of a spinon continuum are observed at low temperatures. A 3D two-magnon density of states now gradually extends to higher temperatures T\leq 60K compared with pure CuGeO_3. This effect is more pronounced in the case of off-chain substitutions (Si) for which a N\'eel state occurs over a larger substitution range, starting at very low concentrations. Besides, additional low energy excitations are induced. These effects, i.e. the shift of a dimensional crossover to higher temperatures are due to an enhancement of the spin-spin correlations induced by a small amount of substitutions. The results are compared with recent Monte Carlo studies on substituted spin ladders, pointing to a similar instability of coupled, dimerized spin chains and spin ladders upon substitution.Comment: 14 pages, 6 eps figures, to be published in PR

Orbital-selective Mott transitions in the anisotropic two-band Hubbard model at finite temperatures

The anisotropic degenerate two-orbital Hubbard model is studied within dynamical mean-field theory at low temperatures. High-precision calculations on the basis of a refined quantum Monte Carlo (QMC) method reveal that two distinct orbital-selective Mott transitions occur for a bandwidth ratio of 2 even in the absence of spin-flip contributions to the Hund exchange. The second transition -- not seen in earlier studies using QMC, iterative perturbation theory, and exact diagonalization -- is clearly exposed in a low-frequency analysis of the self-energy and in local spectra.Comment: 4 pages, 5 figure

Metal--Insulator Transitions in the Falicov--Kimball Model with Disorder

The ground state phase diagrams of the Falicov--Kimball model with local disorder is derived within the dynamical mean--field theory and using the geometrically averaged (''typical'') local density of states. Correlated metal, Mott insulator and Anderson insulator phases are identified. The metal--insulator transitions are found to be continuous. The interaction and disorder compete with each other stabilizing the metallic phase against occurring one of the insulators. The Mott and Anderson insulators are found to be continuously connected.Comment: 6 pages, 7 figure

Quantum critical point in a periodic Anderson model

We investigate the symmetric Periodic Anderson Model (PAM) on a three-dimensional cubic lattice with nearest-neighbor hopping and hybridization matrix elements. Using Gutzwiller's variational method and the Hubbard-III approximation (which corresponds to the exact solution of an appropriate Falicov-Kimball model in infinite dimensions) we demonstrate the existence of a quantum critical point at zero temperature. Below a critical value $V_c$ of the hybridization (or above a critical interaction $U_c$) the system is an {\em insulator} in Gutzwiller's and a {\em semi-metal} in Hubbard's approach, whereas above $V_c$ (below $U_c$) it behaves like a metal in both approximations. These predictions are compared with the density of states of the $d$- and $f$-bands calculated from Quantum Monte Carlo and NRG calculations. Our conclusion is that the half-filled symmetric PAM contains a {\em metal-semimetal transition}, not a metal-insulator transition as has been suggested previously.Comment: ReVteX, 10 pages, 2 EPS figures. Minor corrections made in the text and in the figure captions from the first version. More references added. Accepted for publication in Physical Review

Asymptotics of self-similar solutions to coagulation equations with product kernel

We consider mass-conserving self-similar solutions for Smoluchowski's coagulation equation with kernel $K(\xi,\eta)= (\xi \eta)^{\lambda}$ with $\lambda \in (0,1/2)$. It is known that such self-similar solutions $g(x)$ satisfy that $x^{-1+2\lambda} g(x)$ is bounded above and below as $x \to 0$. In this paper we describe in detail via formal asymptotics the qualitative behavior of a suitably rescaled function $h(x)=h_{\lambda} x^{-1+2\lambda} g(x)$ in the limit $\lambda \to 0$. It turns out that $h \sim 1+ C x^{\lambda/2} \cos(\sqrt{\lambda} \log x)$ as $x \to 0$. As $x$ becomes larger $h$ develops peaks of height $1/\lambda$ that are separated by large regions where $h$ is small. Finally, $h$ converges to zero exponentially fast as $x \to \infty$. Our analysis is based on different approximations of a nonlocal operator, that reduces the original equation in certain regimes to a system of ODE

A Survey of Numerical Solutions to the Coagulation Equation

We present the results of a systematic survey of numerical solutions to the coagulation equation for a rate coefficient of the form A_ij \propto (i^mu j^nu + i^nu j^mu) and monodisperse initial conditions. The results confirm that there are three classes of rate coefficients with qualitatively different solutions. For nu \leq 1 and lambda = mu + nu \leq 1, the numerical solution evolves in an orderly fashion and tends toward a self-similar solution at large time t. The properties of the numerical solution in the scaling limit agree with the analytic predictions of van Dongen and Ernst. In particular, for the subset with mu > 0 and lambda < 1, we disagree with Krivitsky and find that the scaling function approaches the analytically predicted power-law behavior at small mass, but in a damped oscillatory fashion that was not known previously. For nu \leq 1 and lambda > 1, the numerical solution tends toward a self-similar solution as t approaches a finite time t_0. The mass spectrum n_k develops at t_0 a power-law tail n_k \propto k^{-tau} at large mass that violates mass conservation, and runaway growth/gelation is expected to start at t_crit = t_0 in the limit the initial number of particles n_0 -> \infty. The exponent tau is in general less than the analytic prediction (lambda + 3)/2, and t_0 = K/[(lambda - 1) n_0 A_11] with K = 1--2 if lambda > 1.1. For nu > 1, the behaviors of the numerical solution are similar to those found in a previous paper by us. They strongly suggest that there are no self-consistent solutions at any time and that runaway growth is instantaneous in the limit n_0 -> \infty. They also indicate that the time t_crit for the onset of runaway growth decreases slowly toward zero with increasing n_0.Comment: 41 pages, including 14 figures; accepted for publication in J. Phys.

Comparison of Variational Approaches for the Exactly Solvable 1/r-Hubbard Chain

We study Hartree-Fock, Gutzwiller, Baeriswyl, and combined Gutzwiller-Baeriswyl wave functions for the exactly solvable one-dimensional $1/r$-Hubbard model. We find that none of these variational wave functions is able to correctly reproduce the physics of the metal-to-insulator transition which occurs in the model for half-filled bands when the interaction strength equals the bandwidth. The many-particle problem to calculate the variational ground state energy for the Baeriswyl and combined Gutzwiller-Baeriswyl wave function is exactly solved for the~$1/r$-Hubbard model. The latter wave function becomes exact both for small and large interaction strength, but it incorrectly predicts the metal-to-insulator transition to happen at infinitely strong interactions. We conclude that neither Hartree-Fock nor Jastrow-type wave functions yield reliable predictions on zero temperature phase transitions in low-dimensional, i.e., charge-spin separated systems.Comment: 23 pages + 3 figures available on request; LaTeX under REVTeX 3.