Dynamics of selfavoiding tethered membranes. I. Model A dynamics (Rouse model)

The dynamical scaling properties of selfavoiding polymerized membranes with internal dimension D are studied using model A dynamics. It is shown that the theory is renormalizable to all orders in perturbation theory and that the dynamical scaling exponent z is given by . This result applies especially to membranes (D=2) but also to polymers (D=1). Received: 5 September 1997 / Accepted: 17 November 1997     

Quantum rotors with regular frustration and the quantum Lifshitz point

We have discussed the zero-temperature quantum phase transition in n-component quantum rotor Hamiltonian in the presence of regular frustration in the interaction. The phase diagram consists of ferromagnetic, helical and quantum paramagnetic phase, where the ferro-para and the helical-para phase boundary meets at a multicritical point called a (d,m) quantum Lifshitz point where (d,m) indicates that the m of the d spatial dimensions incorporate frustration. We have studied the Hamiltonian in the vicinity of the quantum Lifshitz point in the spherical limit and also studied the renormalisation group flow behaviour using standard momentum space renormalisation technique (for finite n). In the spherical limit ()one finds that the helical phase does not exist in the presence of any nonvanishing quantum fluctuation for m =d though the quantum Lifshitz point exists for all d > 1+m/2, and the upper critical dimensionality is given by d _u = 3 +m/2. The scaling behaviour in the neighbourhood of a quantum Lifshitz point in d dimensions is consistent with the behaviour near the classical Lifshitz point in (d+z) dimensions. The dynamical exponent of the quantum Hamiltonian z is unity in the case of anisotropic Lifshitz point (d>m) whereas z=2 in the case of isotropic Lifshitz point (d=m). We have evaluated all the exponents using the renormalisation flow equations along-with the scaling relations near the quantum Lifshitz point. We have also obtained the exponents in the spherical limit (). It has also been shown that the exponents in the spherical model are all related to those of the corresponding Gaussian model by Fisher renormalisation. Received: 23 December 1997 / Received in final form: 6 January 1998 / Accepted: 7 January 1998     

Dynamics of polymeric manifolds in melts: the Hartree approximation

The Martin-Siggia-Rose functional technique and the selfconsistent Hartree approximation is applied to the dynamics of a D-dimensional manifold in a melt of similar manifolds. The generalized Rouse equation is derived and its static and dynamic properties are studied. The static upper critical dimension, d _uc =2D/(2-D), discriminates between Gaussian (or screened) and non-Gaussian regimes, whereas its dynamical counterpart, , discriminates between Rouse- and renormalized-Rouse behavior. The Rouse modes correlation function in a stretched exponential form and the dynamical exponents are calculated explicitly. The special case of linear chains D=1 shows agreement with Monte-Carlo simulations. Received: 22 May 1998 / Received in final form: 31 August 1998 / Accepted: 8 September 1998     

Single-cluster Monte Carlo dynamics for the Ising model

We present an extensive study of a new Monte Carlo acceleration algorithm introduced by Wolff for the Ising model. It differs from the Swendsen-Wang algorithm by growing and flipping single clusters at a random seed. In general, it is more efficient than Swendsen-Wang dynamics ford>2, giving zero critical slowing down in the upper critical dimension. Monte Carlo simulations give dynamical critical exponentsz _w=0.33±0.05 and 0.44+0.10 ind=2 and 3, respectively, and numbers consistent withz _w=0 ind=4 and mean-field theory. We present scaling arguments which indicate that the Wolff mechanism for decorrelation differs substantially from Swendsen-Wang despite the apparent similarities of the two methods.     

Numerical diagonalization analysis of the ground-state superfluid-localization transition in two dimensions

Ground state of the two-dimensional hard-core-boson system in the presence of the quenched random chemical potential is investigated by means of the exact-diagonalization method for the system sizes up to L=5. The criticality and the DC conductivity at the superfluid-localization transition have been controversial so far. We estimate, with the finite-size scaling analysis, the correlation-length and the dynamical critical exponents as and z=2, respectively. The AC conductivity is computed with the Gagliano-Balseiro formula, with which the resolvent (dynamical response function) is expressed in terms of the continued-fraction form consisting of Lanczos tri-diagonal elements. Thereby, we estimate the universal DC conductivity as . Received 19 August 1998     

A geometric generalization of field theory to manifolds of arbitrary dimension

We introduce a generalization of the O(N) field theory to N-colored membranes of arbitrary inner dimension D. The O(N) model is obtained for , while leads to self-avoiding tethered membranes (as the O(N) model reduces to self-avoiding polymers). The model is studied perturbatively by a 1-loop renormalization group analysis, and exactly as .Freedom to choose the expansion point D, leads to precise estimates of critical exponents of the O(N) model. Insights gained from this generalization include a conjecture on the nature of droplets dominating the 3d-Ising model at criticality; and the fixed point governing the random bond Ising model. Received: 15 October 1998 / Accepted: 4 November 1998     

Topology invariance in percolation thresholds

An universal invariant for site and bond percolation thresholds ( and respectively) is proposed. The invariant writes where and are positive constants, and d the space dimension. It is independent of the coordination number, thus exhibiting a topology invariance at any d. The formula is checked against a large class of percolation problems, including percolation in non-Bravais lattices and in aperiodic lattices as well as rigid percolation. The invariant is satisfied within a relative error of for all the twenty lattices of our sample at d=2, d=3, plus all hypercubes up to d=6. Received: 7 July 1997 / Accepted: 5 November 1997     

Sensitivity to initial conditions in the Bak-Sneppen model of biological evolution

We consider biological evolution as described within the Bak and Sneppen 1993 model. We exhibit, at the self-organized critical state, a power-law sensitivity to the initial conditions, calculate the associated exponent, and relate it to the recently introduced nonextensive thermostatistics. The scenario which here emerges without tuning strongly reminds of that of the tuned onset of chaos in say logistic-like one-dimensional maps. We also calculate the dynamical exponent z. Received: 5 November 1997 / Received in final form: 11 November 1997 / Accepted: 19 November 1997     

Critical dynamics of the three-dimensional Ising model: A Monte Carlo study

Extensive Monte Carlo simulations have been performed to analyze the dynamical behavior of the three-dimensional Ising model with local dynamics. We have studied the equilibrium correlation functions and the power spectral densities of odd and even observables. The exponential relaxation times have been calculated in the asymptotic one-exponential time region. We find that the critical exponentz=2.09 ±0.02 characterizes the algebraic divergence with lattice size for all observables. The influence of scaling corrections has been analyzed. We have determined integrated relaxation times as well. Their dynamical exponentz _int agrees withz for correlations of the magnetization and its absolute value, but it is different for energy correlations. We have applied a scaling method to analyze the behavior of the correlation functions. This method verifies excellent scaling behavior and yields a dynamical exponentz _scal which perfectly agrees withz.     

Circular-like maps: sensitivity to the initial conditions, multifractality and nonextensivity

Dissipative one-dimensional maps may exhibit special points (e.g., chaos threshold) at which the Lyapunov exponent vanishes. Consistently, the sensitivity to the initial conditions has a power-law time dependence, instead of the usual exponential one. The associated exponent can be identified with 1/(1-q), where q characterizes the nonextensivity of a generalized entropic form currently used to extend standard, Boltzmann-Gibbs statistical mechanics in order to cover a variety of anomalous situations. It has been recently proposed (Lyra and Tsallis, Phys. Rev. Lett. 80, 53 (1998)) for such maps the scaling law , where and are the extreme values appearing in the multifractal function. We generalize herein the usual circular map by considering inflexions of arbitrary power z, and verify that the scaling law holds for a large range of z. Since, for this family of maps, the Hausdorff dimension d_f equals unity for all z in contrast with q which does depend on z, it becomes clear that d_f plays no major role in the sensitivity to the initial conditions. Received 5 February 1999     

Uniform susceptibility of classical antiferromagnets in one and two dimensions in a magnetic field

We simulated the field-dependent magnetization m(H,T) and the uniform susceptibility of classical Heisenberg antiferromagnets in the chain and square-lattice geometry using Monte Carlo methods. The results confirm the singular behavior of at small T,H: and , where D=3 is the number of spin components, J ₀=zJ, and z is the number of nearest neighbors. A good agreement is achieved in a wide range of temperatures T and magnetic fields H with the first-order 1/D expansion results (D.A. Garanin, J. Stat. Phys. 83, 907 (1996)). Received 20 March 2000     

Interfacial reaction kinetics

We study irreversible A-B reaction kinetics at a fixed interface separating two immiscible bulk phases, A and B. Coupled equations are derived for the hierarchy of many-body correlation functions. Postulating physically motivated bounds, closed equations result without the need for ad hoc decoupling approximations. We consider general dynamical exponent z, where is the rms diffusion distance after time t. At short times the number of reactions per unit area, , is 2nd order in the far-field reactant densities . For spatial dimensions dabove a critical value , simple mean field (MF) kinetics pertain, where Q_b is the local reactivity. For low dimensions , this MF regime is followed by 2nd order diffusion controlled (DC) kinetics, , provided . Logarithmic corrections arise in marginal cases. At long times, a cross-over to 1st order DC kinetics occurs: . A density depletion hole grows on the more dilute A side. In the symmetric case (), when the long time decay of the interfacial reactant density, , is determined by fluctuations in the initial reactant distribution, giving . Correspondingly, A-rich and B-rich regions develop at the interface analogously to the segregation effects established by other authors for the bulk reaction . For fluctuations are unimportant: local mean field theory applies at the interface (joint density distribution approximating the product of A and B densities) and . We apply our results to simple molecules (Fickian diffusion, z=2) and to several models of short-time polymer diffusion (z>2). Received 8 June 1998 and Received in final form 10 September 1999     

Effect of hydrodynamic flow on kinetics of nematic-isotropic transition in liquid crystals

We investigate kinetics of nematic-isotropic transition by solving the hydrodynamic equations for the nematic tensor order parameter and the fluid velocity in two space dimension (x-y plane). Numerical results indicate that nematic directors tend to align parallel to the x-y plane when hydrodynamic flow is incorporated. Late stage growth exponents, for the correlation length and for the number of topological defects, are not significantly altered by hydrodynamic flow. However, in contrast to the case without flow, the relation holds well, which may indicate the validity of dynamical scaling for the case with hydrodynamic flow. Received: 8 September 1997 / Received in final form: 23 October 1997 /Accepted: 3 November 1997     

Nuclear fusion in charge-asymmetric muonic molecules

Nuclear fusion reactions in hydrogen-lithium muonic molecules, (where h=p,d,t) are considered and fusion rates from rotational states J=0 of the molecules are presented. Results obtained depend on the isotopic composition of the molecules and range between and . The upper limit for fusion rates from rotational states J=0 of hydrogen-helium muonic molecules, and , equal , is also found. Received: 4 December 1997 / Revised: 30 April 1998 / Accepted: 7 May 1998     

Noise-assisted mound coarsening in epitaxial growth

Two types of mechanisms are proposed for mound coarsening during unstable epitaxial growth: stochastic, due to deposition noise, and deterministic, due to mass currents driven by surface energy differences. Both yield the relation H=(RWL)² between the typical mound height W, mound size L, and the film thickness H. An analysis of simulations and experimental data shows that the parameter R saturates to a value which discriminates sharply between stochastic () and deterministic () coarsening. We derive a scaling relation between the coarsening exponent 1/z and the mound-height exponent which, for a saturated mound slope, yields . Received: 11 November 1997 / Revised in final form: 28 November 1997 / Accepted: 28 November 1997     

Scaling of Loop-Erased Walks in 2 to 4 Dimensions

We simulate loop-erased random walks on simple (hyper-)cubic lattices of dimensions 2, 3 and 4. These simulations were mainly motivated to test recent two loop renormalization group predictions for logarithmic corrections in d=4, simulations in lower dimensions were done for completeness and in order to test the algorithm. In d=2, we verify with high precision the prediction D=5/4, where the number of steps n after erasure scales with the number N of steps before erasure as n∼N ^D/2. In d=3 we again find a power law, but with an exponent different from the one found in the most precise previous simulations: D=1.6236±0.0004. Finally, we see clear deviations from the naive scaling n∼N in d=4. While they agree only qualitatively with the leading logarithmic corrections predicted by several authors, their agreement with the two-loop prediction is nearly perfect.     

Born effective charge reversal and metallic threshold state at a band insulator-Mott insulator transition

We study the quantum phase transition between a band (“ionic”) insulator and a Mott-Hubbard insulator, realized at a critical value in a bipartite Hubbard model with two inequivalent sites, whose on-site energies differ by an offset . The study is carried out both in D=1 and D=2 (square and honeycomb lattices), using exact Lanczos diagonalization, finite-size scaling, and Berry's phase calculations of the polarization. The Born effective charge jump from positive infinity to negative infinity previously discovered in D=1 by Resta and Sorella is confirmed to be directly connected with the transition from the band insulator to the Mott insulating state, in agreement with recent work of Ortiz et al. In addition, symmetry is analysed, and the transition is found to be associated with a reversal of inversion symmetry in the ground state, of magnetic origin. We also study the D=1 excitation spectrum by Lanczos diagonalization and finite-size scaling. Not only the spin gap closes at the transition, consistent with the magnetic nature of the Mott state, but also the charge gap closes, so that the intermediate state between the two insulators appears to be metallic. This finding, rationalized within Hartree-Fock as due to a sign change of the effective on-site energy offset for the minority spin electrons, underlines the profound difference between the two insulators. The band-to-Mott insulator transition is also studied and found in the same model in D=2. There too we find an associated, although weaker, polarization anomaly, with some differences between square and honeycomb lattices. The honeycomb lattice, which does not possess an inversion symmetry, is used to demonstrate the possibility of an inverted piezoelectric effect in this kind of ionic Mott insulator. Received 21 May 1999     

Instanton Calculus for the Self-Avoiding Manifold Model

We compute the normalisation factor for the large order asymptotics of perturbation theory for the self-avoiding manifold (SAM) model describing flexible tethered (D-dimensional) membranes in d-dimensional space, and the ε-expansion for this problem. For that purpose, we develop the methods inspired from instanton calculus, that we introduced in a previous publication (Nucl. Phys. B 534 (1998) 555), and we compute the functional determinant of the fluctuations around the instanton configuration. This determinant has UV divergences and we show that the renormalized action used to make perturbation theory finite also renders the contribution of the instanton UV-finite. To compute this determinant, we develop a systematic large-d expansion. For the renormalized theory, we point out problems in the interplay between the limits ε→ 0 and d→∞, as well as IR divergences when ε=0. We show that many cancellations between IR divergences occur, and argue that the remaining IR-singular term is associated to amenable non-analytic contributions in the large-d limit when ε=0. The consistency with the standard instanton-calculus results for the self-avoiding walk is checked for D=1.