Scaling and renormalization group in replica-symmetry-breaking space: evidence for a simple analytical solution of the Sherrington-Kirkpatrick model at zero temperature

Using numerical self-consistent solutions of a sequence of finite replica symmetry breakings (RSB) and Wilson's renormalization group but with the number of RSB steps playing a role of decimation scales, we report evidence for a nontrivial T-->0 limit of the Parisi order function q(x) for the Sherrington-Kirkpatrick spin glass. Supported by scaling in RSB space, the fixed point order function is conjectured to be q*(a)=sqrt[pi]/2 a/xi erf(xi/a) on 0 a at T =0 and xi approximately 1.13+/-0.01. Xi plays the role of a correlation length in a-space. q*(a) may be viewed as the solution of an effective 1D field theory.     

Interfacial reactions: mixed order kinetics and segregation effects

We study A-B reaction kinetics at a fixed interface separating A and B bulks. Initially, the number of reactions R(t) approximately tn(infinity)(A)n(infinity)(B) is second order in the far-field densities n(infinity)(A), n(infinity)(B). First order kinetics, governed by diffusion from the dilute bulk, onset at long times: R(t) approximately x(t)n(infinity)(A), where x(t) approximately t(1/z) is the rms molecular displacement. Below a critical dimension, d0) leads to anomalous decay of interfacial densities. Numerical simulations for z = 2 support the theory.     

Maximal height scaling of kinetically growing surfaces

The scaling properties of the maximal height of a growing self-affine surface with a lateral extent L are considered. In the late-time regime its value measured relative to the evolving average height scales like the roughness: h*(L) approximately L alpha. For large values its distribution obeys logP(h*(L)) approximately (-)A(h*(L)/L(alpha))(a). In the early-time regime where the roughness grows as t(beta), we find h*(L) approximately t(beta)[lnL-(beta/alpha)lnt+C](1/b), where either b = a or b is the corresponding exponent of the velocity distribution. These properties are derived from scaling and extreme-value arguments. They are corroborated by numerical simulations and supported by exact results for surfaces in 1D with the asymptotic behavior of a Brownian path.     

Critical magnetic Prandtl number for small-scale dynamo

We report a series of numerical simulations showing that the critical magnetic Reynolds number Rm(c) for the nonhelical small-scale dynamo depends on the Reynolds number Re. Namely, the dynamo is shut down if the magnetic Prandtl number Pr(m)=Rm/Re is less than some critical value Pr(m,c)< approximately 1 even for Rm for which dynamo exists at Pr(m)> or =1. We argue that, in the limit of Re-->infinity, a finite Pr(m,c) may exist. The second possibility is that Pr(m,c)-->0 as Re--> infinity, while Rm(c) tends to a very large constant value inaccessible at current resolutions. If there is a finite Pr(m,c), the dynamo is sustainable only if magnetic fields can exist at scales smaller than the flow scale, i.e., it is always effectively a large-Pr(m) dynamo. If there is a finite Rm(c), our results provide a lower bound: Rm(c) greater, similar 220 for Pr(m)< or =1/8. This is larger than Rm in many planets and in all liquid-metal experiments.     

Low-temperature broken-symmetry phases of spiral antiferromagnets

We study Heisenberg antiferromagnets with nearest- (J1) and third- (J3) neighbor exchange on the square lattice. In the limit of spin S-->infinity, there is a zero temperature (T) Lifshitz point at J(3)=1/4J(1), with long-range spiral spin order at T=0 for J3>1/4J(1). We present classical Monte Carlo simulations and a theory for T>0 crossovers near the Lifshitz point: spin rotation symmetry is restored at any T>0, but there is a broken lattice reflection symmetry for 0< or =T相似文献   

Solitary wave solutions of nonlocal sine-Gordon equations

In this paper a nonlocal generalization of the sine-Gordon equation, u(tt)+sin u=( partial differential / partial differential x) integral (- infinity ) (+ infinity )G(x-x('))u(x(') )(x('),t)dx(') is considered. We present a brief review of the applications of such equations and show that involving such a nonlocality can change features of the model. In particular, some solutions of the sine-Gordon model (for example, traveling 2pi-kink solutions) may disappear in the nonlocal model; furthermore, some new classes of solutions such as traveling topological solitons with topological charge greater than 1 may arise. We show that the lack of Lorenz invariancy of the equation under consideration can lead to a phenomenon of discretization of kink velocities. We discussed this phenomenon in detail for the special class of kernels G(xi)= summation operator (j=1) (N)kappa(j)e(-eta(j)mid R:ximid R:), eta(j)>0, j=1,2, em leader,N. We show that, generally speaking, in this case the velocities of kinks (i) are determined unambiguously by a type of kink and value(s) of kernel parameter(s); (ii) are isolated i.e., if c(*) is the velocity of a kink then there are no other kink solutions of the same type with velocity c in (c(*)- varepsilon,c(*)+ varepsilon ) for a certain value of varepsilon. We also used this special class of kernels to construct approximations for analytical and numerical study of the problem in a more general case. Finally, we set forth results of the numerical investigation of the problem with the kernel that is the McDonald function G(xi) approximately K(0)(mid R:ximid R:/lambda) (lambda is a parameter) that have applications in the Josephson junction theory. (c) 1998 American Institute of Physics.     

Disappearance of the gas-liquid phase transition for highly charged colloids

We calculate the full phase diagram of spherical charged colloidal particles using Monte Carlo free energy calculations. The system is described using the primitive model, consisting of explicit hard-sphere colloids and point counterions in a uniform dielectric continuum. We show that the gas-liquid critical point becomes metastable with respect to a gas-solid phase separation at colloid charges Q > or =20 times the counterion charge. Approximate free energy calculations with only one and four particles in the fluid and solid phases, respectively, are used to determine the critical line for highly charged colloids up to Q=2000. We propose the scaling law T*(c) approximately Q(1/2) for this critical temperature.     

Renormalization of pinned elastic systems: how does it work beyond one loop?

We study the field theories for pinned elastic systems at equilibrium and at depinning. Their beta functions differ to two loops by novel "anomalous" terms. At equilibrium we find a roughness zeta = 0.208 298 04 epsilon + 0.006 858 epsilon(2) (random bond), zeta = epsilon/3 (random field). At depinning we prove two-loop renormalizability and that random field attracts shorter range disorder. We find zeta = epsilon/3(1 + 0.143 31 epsilon), epsilon = 4 - d, in violation of the conjecture zeta = epsilon/3, solving the discrepancy with simulations. For long range elasticity zeta = epsilon/3(1 + 0.397 35 epsilon), epsilon = 2 - d, much closer to the experimental value (approximately 0.5 both for liquid helium contact line depinning and slow crack fronts) than the standard prediction 1/3.     

Autocorrelation exponent of conserved spin systems in the scaling regime following a critical quench

We study the autocorrelation function of a conserved spin system following a quench at the critical temperature. Defining the correlation length L(t) approximately t(1/z), we find that for times t' and t satisfying L(t')infinity limit, we show that lambda(')(c)=d+2 and phi=z/2. We give a heuristic argument suggesting that this result is, in fact, valid for any dimension d and spin vector dimension n. We present numerical simulations for the conserved Ising model in d=1 and d=2, which are fully consistent with the present theory.     

Ballistic annihilation with continuous isotropic initial velocity distribution

Ballistic annihilation with continuous initial velocity distributions is investigated in the framework of the Boltzmann equation. The particle density and the rms velocity decay as c approximately t(-alpha) and velocity approximately t(-beta), with the exponents depending on the initial velocity distribution and the spatial dimension d. For instance, in one dimension for the uniform initial velocity distribution beta = 0.230 472ellipsis. In the opposite extreme d-->infinity, the dynamics is universal and beta-->(1-2(-1/2))d(-1). We also solve the Boltzmann equation for Maxwell particles and very hard particles in arbitrary spatial dimension. These solvable cases provide bounds for the decay exponents of the hard sphere gas.     

Qualitative modification of the high energy atomic photoionization cross section

It is demonstrated that the nonrelativistic high energy omega-->infinity behavior of the photoionization cross section of an nl atomic subshell, sigma(nl)(omega), for l>0 is independent of l and is given by sigma(nl)(omega) approximately 1/omega(9/2), rather than the previously generally accepted sigma(nl)(omega) approximately 1/omega(l+7/2). Furthermore, for l = 1, although the exponent does not change, the coefficient is significantly altered. This modification of sigma(nl)(omega) is due to the interchannel interaction between ns photoionization channels and l not equal0 channels in the atom. As a result, for the photoionization of l not equal0 electrons, the single-particle approximation is never correct in the omega-->infinity limit. This has important consequences for sum rule calculations.     

Screening in ionic systems: simulations for the Lebowitz length

Simulations of the Lebowitz length, xiL (T, rho), are reported for the restricted primitive model hard-core (diameter a) 1:1 electrolyte for densities rho approximately < 4rho(c) and T(c) approximately < T approximately < 40T(c). Finite-size effects are elucidated for the charge fluctuations in various subdomains that serve to evaluate xiL. On extrapolation to the bulk limit for T approximately > 10T(c) the exact low-density expansions are seen to fail badly when rho > 1/10 rho(c) (with rho(c)a3 approximately = 0.08). At higher densities xiL rises above the Debye length, xiD proportional to square root(T/rho), by 10%-30% (up to rho approximately =1.3rho(c)); the variation is portrayed fairly well by the generalized Debye-Hückel theory. On approaching criticality at fixed rho or fixed T, xiL (T, rho) remains finite with xiL(c) approximately = 0.30a approximately = 1.3xiD(c) but displays a weak entropylike singularity.     

Universality class of criticality in the restricted primitive model electrolyte

The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using loci in the temperature-density or (T,rho) plane of isothermal rho(2-k) vs pressure inflections, of Q identical with(2)/ maxima, and of canonical and C(V) criticality, yield estimates of (T(c),rho(c)) to +/-(0.04,3)%. Extrapolated exponents and Q ratio are (gamma,nu,Q(c)) = [1.24(3), 0.63(3); 0.624(2)], which support Ising (n = 1) behavior with (1.23(9), 0.630(3); 0.623(6)), but exclude classical, XY (n = 2), self-avoiding walk (n = 0), and n = 1 criticality with potentials varphi(r)>Phi/r(4.9) when r-->infinity.     

Peculiar scaling of self-avoiding walk contacts

The nearest neighbor contacts between the two halves of an N-site lattice self-avoiding walk offer an unusual example of scaling random geometry: for N-->infinity they are strictly finite in number but their radius of gyration R(c) is power law distributed proportional to R(-tau)(c), where tau>1 is a novel exponent characterizing universal behavior. A continuum of diverging length scales is associated with the R(c) distribution. A possibly superuniversal tau = 2 is also expected for the contacts of a self-avoiding or random walk with a confining wall.     

Photon-stimulated ion desorption from DCOO/Si(100) near the O K-edge

Photon-stimulated ion desorption from deuterated formic acid chemisorbed on Si(100) has been studied using pulsed synchrotron radiation in the energy region of the oxygen 1s electron excitation. The O 1s electrons of hydroxyl oxygen and carbonyl oxygen could be selectively excited in the O K-edge region because the chemical environments are different. It is found that the CDO⁺ yield is enhanced at the O 1s(C---O)→σ*(C---O) resonance and the CD⁺ yield at the O 1s(C=O)→σ*(C---O) resonance. The results indicate that ion desorption is related both to the antibonding character of excited molecular orbitals and the local character of core hole orbitals.     

Unconventional superconductivity in electron-doped layered Li0.48(THF)yHfNCL.

We report magnetic susceptibility measurements on a layered superconductor Li0.48(THF)0.3HfNCl having Tc approximately 26 K. The present study revealed that (a) the Fermi level density of states is small, N*(EF) approximately 0.25 states/(eV spin f.u.), (b) mass enhancement is negligible, gamma; approximately 1, (c) electron-phonon coupling is weak, lambda(ep)<1, (d) exchange enhancement is negligible, 1/(1+F(a)0) approximately 1, and (e) electronic density parameter is large, r(2D)s approximately 10.3 (i.e., low-carrier density). It is difficult to explain the origin of the high Tc in terms of the conventional phonon (BCS) mechanism of superconductivity.     

Connectivity of growing random networks

A solution for the time- and age-dependent connectivity distribution of a growing random network is presented. The network is built by adding sites that link to earlier sites with a probability A(k) which depends on the number of preexisting links k to that site. For homogeneous connection kernels, A(k) approximately k(gamma), different behaviors arise for gamma<1, gamma>1, and gamma = 1. For gamma<1, the number of sites with k links, N(k), varies as a stretched exponential. For gamma>1, a single site connects to nearly all other sites. In the borderline case A(k) approximately k, the power law N(k) approximately k(-nu) is found, where the exponent nu can be tuned to any value in the range 2相似文献   

Persistence in systems with algebraic interaction

Persistence in coarsening one-dimensional spin systems with a power-law interaction r(-1-sigma) is considered. Numerical studies indicate that for sufficiently large values of the interaction exponent sigma (sigma > or =1/2 in our simulations), persistence decays as an algebraic function of the length scale L, P(L) approximately L(-theta). The persistence exponent theta is found to be independent on the force exponent sigma and close to its value for the extremal (sigma-->infinity) model, theta =0.175 075 88. For smaller values of the force exponent (sigma < 1/2), finite size effects prevent the system from reaching the asymptotic regime. Scaling arguments suggest that in order to avoid significant boundary effects for small sigma, the system size should grow as [O(1/sigma)](1/sigma).     

Can nonlinear elasticity explain contact-line roughness at depinning?

We examine whether cubic nonlinearities, allowed by symmetry in the elastic energy of a contact line, may result in a different universality class at depinning. Standard linear elasticity predicts a roughness exponent zeta = 1/3 (one loop), zeta = 0.388 +/- 0.002 (numerics) while experiments give zeta approximately = 0.5. Within functional renormalization group methods we find that a nonlocal Kardar-Parisi-Zhang-type term is generated at depinning and grows under coarse graining. A fixed point with zeta approximately = 0.45 (one loop) is identified, showing that large enough cubic terms increase the roughness. This fixed point is unstable, revealing a rough strong-coupling phase. Experimental study of contact angles theta near pi/2, where cubic terms in the energy vanish, is suggested.     

Scaling properties of three-dimensional magnetohydrodynamic turbulence

The scaling properties of three-dimensional magnetohydrodynamic turbulence with finite magnetic helicity are obtained from direct numerical simulations using 512(3) modes. The results indicate that the turbulence does not follow the Iroshnikov-Kraichnan phenomenology. The scaling exponents of the structure functions can be described by a modified She-Leveque model zeta(p) = p/9+1-(1/3)(p/3), corresponding to basic Kolmogorov scaling and sheetlike dissipative structures. In particular, we find zeta(2) approximately 0.7, consistent with the energy spectrum E(k) approximately k(-5/3) as observed in the solar wind, and zeta(3) approximately 1, confirming a recent analytical result.