Critical dynamics of cluster algorithms in the dilute Ising model

Autocorrelation times for thermodynamic quantities atT _C are calculated from Monte Carlo simulations of the site-diluted simple cubic Ising model, using the Swendsen-Wang and Wolff cluster algorithms. Our results show that for these algorithms the autocorrelation timesdecrease when reducing the concentration of magnetic sites from 100% down to 40%. This is of crucial importance when estimating static properties of the model, since the variances of these estimators increase with autocorrelation time. The dynamical critical exponents are calculated for both algorithms, observing pronounced finite-size effects in the energy autocorrelation data for the algorithm of Wolff. We conclude that, when applied to the dilute Ising model, cluster algorithms become even more effective than local algorithms, for whichincreasing autocorrelation times are expected.     

Inequalities among static and dynamic critical-point exponents

The mode-mode coupling theoretical expression for the decay rate of order-parameter fluctuations in simple fluids near the critical point is evaluated with the aid of the recently extended Ornstein-Zernike theory. This gives that ψ?ν for the critical-point exponent ψ which characterize the behavior of the thermal conductivity. Also, the results for sound-wave absorption and dispersion near the gas-liquid critical point need not depend on a single reduced frequency variable since scaling may be violated.     

Some relations between dimension and Lyapounov exponents

We consider differentiable maps and compact invariant sets. We introduce dimensional quantities related to the ergodic invariant measures, and prove some simple relations.     

Anomalous dimension exponents on fractal structures for the Ising and three-state Potts model

We provide an overall picture of the magnetic critical behavior of the Ising and three-state Potts models on fractal structures. The results brought out from Monte Carlo simulations for several Hausdorff dimensions between 1 and 3 show that this behavior can be understood in the framework of weak universality. Moreover, the maxima of the susceptibility follow power laws in a very reliable way, which allows us to calculate the ratio of the exponents γ/ν and the anomalous dimension exponent η in a reliable way. At last, the evolution of these exponents with the Hausdorff dimension is discussed.     

Critical exponents for a three-dimensional impure Ising model in the five-loop approximation

The renormalization-group functions governing the critical behavior of the three-dimensional weakly-disordered Ising model are calculated in the five-loop approximation. The random fixed point location and critical exponents for impure Ising systems are estimated by means of the Padé-Borel-Leroy resummation of the renormalization-group expansions derived. The asymptotic critical exponents are found to be γ=1.325 ± 0.003, η=0.025 ± 0.01, ν= 0.671 ± 0.005, α=?0.0125 ± 0.008, β=0.344 ± 0.006, while for the correction-to-scaling exponent, a less accurate estimate ω=0.32 ± 0.06 is obtained.     

Critical exponents of the Ising model on low-dimensional fractal media

The critical behavior of the Ising model on fractal substrates with noninteger Hausdorff dimension d_H<2 and infinite ramification order is studied by means of the short-time critical dynamic scaling approach. Our determinations of the critical temperatures and critical exponents β, γ, and ν are compared to the predictions of the Wilson-Fisher expansion, the Wallace-Zia expansion, the transfer matrix method, and more recent Monte Carlo simulations using finite-size scaling analysis. We also determined the effective dimension (d_ef), which plays the role of the Euclidean dimension in the formulation of the dynamic scaling and in the hyperscaling relationship d_ef=2β/ν+γ/ν. Furthermore, we obtained the dynamic exponent z of the nonequilibrium correlation length and the exponent θ that governs the initial increase of the magnetization. Our results are consistent with the convergence of the lower-critical dimension towards d=1 for fractal substrates and suggest that the Hausdorff dimension may be different from the effective dimension.     

Exact relations for correlation functions in the general Ising model

The exact relations for the correlation functions in the Ising model are obtained in the simplest form by the Green functions method for the arbitrary spin S and for the arbitrary coordination lattice number P.     

A model for distinguishing between static and dynamic exposure of personal thermoluminescence dosimeters

To distinguish between static and dynamic (normal) exposure of personal TL dosimeters, a model of radiation deposition and of TL light transport in the TLD dosimeter is proposed. The RADOS dosimeter badge using MCP-N (LiF:Mg,Cu,P) TL detectors with standard filters replaced by special Pb and Cu filters with a pattern of holes or inserts was modelled. The photon radiation transport in the dosimeter and energy deposition in the TL detector were simulated by the Penelope Monte Carlo transport code. The model of TL light transport within the TL pellet takes into account the distribution of energy deposition in the TL detector, light self-absorption in the detector and reflection of TL light of the heating planchet. The shape and hole pattern of the filters were optimized with respect to best distinction between static and dynamic exposures. The results of calculations were verified experimentally by exposing RADOS badges with modified filters to beams of low energy X-rays directed at various angles.     

Rigorous bounds of the Lyapunov exponents of the one-dimensional random Ising model

We find analytic upper and lower bounds of the Lyapunov exponents of the product of random matrices related to the one-dimensional disordered Ising model, using a deterministic map which transforms the original system into a new one with smaller average couplings and magnetic fields. The iteration of the map gives bounds which estimate the Lyapunov exponents with increasing accuracy. We prove, in fact, that both the upper and the lower bounds converge to the Lyapunov exponents in the limit of infinite iterations of the map. A formal expression of the Lyapunov exponents is thus obtained in terms of the limit of a sequence. Our results allow us to introduce a new numerical procedure for the computation of the Lyapunov exponents which has a precision higher than Monte Carlo simulations.     

Asymptotic behavior of Mayer cluster sums for the one-dimensional Ising model

The properties of the high-field polynomialsL _n (u) for the one-dimensional spin 1/2 Ising model are investigated. [The polynomialsL _n (u) are essentially lattice gas analogues of the Mayer cluster integralsb _n (T) for a continuum gas.] It is shown thatu ^?1 L _n (u) can be expressed in terms of a shifted Jacobi polynomial of degreen?1. From this result it follows thatu ^?1 L _n (u); n=1, 2,... is a set of orthogonal polynomials in the interval (0, 1) with a weight functionω(u)=u, andu ^?1 L _n (u) hasn?1 simple zerosu _n (v); v=1, 2,...,n?1 which all lie in the interval 0<u<1. Next the detailed behavior ofL _n (u) asn→∞ is studied. In particular, various asymptotic expansions forL _n (u) are derived which areuniformly valid in the intervalsu<0, 0<u<1, andu>1. These expansions are then used to analyze the asymptotic properties of the zeros {u _n (v); v=1, 2,...,n?1}. It is found that $$\begin{array}{*{20}c} {u_n (v) \sim \tfrac{1}{4}({{j_{1,v} } \mathord{\left/ {\vphantom {{j_{1,v} } n}} \right. \kern-\nulldelimiterspace} n})^2 [1 - ({{j_{1,v}^2 } \mathord{\left/ {\vphantom {{j_{1,v}^2 } {12}}} \right. \kern-\nulldelimiterspace} {12}})n^{ - 1} + ({{j_{1,v}^2 } \mathord{\left/ {\vphantom {{j_{1,v}^2 } {700)( - 3 + 2j_{1,v}^2 )n^{ - 4} }}} \right. \kern-\nulldelimiterspace} {700)( - 3 + 2j_{1,v}^2 )n^{ - 4} }}} \\ { + ({{j_{1,v}^2 } \mathord{\left/ {\vphantom {{j_{1,v}^2 } {20160)(40 + 4j_{1,v}^2 - j_{1,v}^4 }}} \right. \kern-\nulldelimiterspace} {20160)(40 + 4j_{1,v}^2 - j_{1,v}^4 }})n^{ - 6} + \cdot \cdot \cdot ]} \\ {u_n (n - v) \sim 1 - ({{j_{0,v}^2 } \mathord{\left/ {\vphantom {{j_{0,v}^2 } 4}} \right. \kern-\nulldelimiterspace} 4})n^{ - 2} + ({{j_{0,v}^2 } \mathord{\left/ {\vphantom {{j_{0,v}^2 } {48)( - 2 + j_{0,v}^2 )n^{ - 4} }}} \right. \kern-\nulldelimiterspace} {48)( - 2 + j_{0,v}^2 )n^{ - 4} }}} \\ { + ({{j_{0,v}^2 } \mathord{\left/ {\vphantom {{j_{0,v}^2 } {2880)(2 + 9j_{0,v}^2 - 2j_{0,v}^4 )n^{ - 6} + \cdot \cdot \cdot }}} \right. \kern-\nulldelimiterspace} {2880)(2 + 9j_{0,v}^2 - 2j_{0,v}^4 )n^{ - 6} + \cdot \cdot \cdot }}} \\ \end{array} $$ asn→∞v fixed, wherej _k,v denotes thevth zero of the Bessel functionJ _k(^z)     

Exact relations between multifractal exponents at the Anderson transition

Two exact relations between mutlifractal exponents are shown to hold at the critical point of the Anderson localization transition. The first relation implies a symmetry of the multifractal spectrum linking the exponents with indices q<1/2 to those with q>1/2. The second relation connects the wave-function multifractality to that of Wigner delay times in a system with a lead attached.     

The infinite cluster method in the two-dimensional Ising model

By studying infinite clusters in the two dimensional ferromagnetic Ising model some new results on the problem of existence of non-translation invariant equilibrium states are obtained. Furthermore a new proof of a theorem by Abraham and Reed is given.This work was supported in part by C.N.R., G.N.F.M.     

A mean-field equation of motion for the dynamic Ising model

A mean-field type of approximation is used to derive two differential equations, one approximately representing the average behavior of the Ising model with Glauber (spin-flip) stochastic dynamics, and the other doing the same for Kawasaki (spin-exchange) dynamics. The proposed new equations are compared with the Cahn-Allen and Cahn-Hilliard equations representing the same systems and with information about the exact behavior of the microscopic models.     

Empirical relations between noncommuting observables

A relation between noncommuting 1-0 quantum observables (i.e., projections) is introduced, being the state vector of the system. This relation extends the empirical implication between commuting projections. An operational interpretation of the new relation is given, which can be expressed also in counterfactual terms. It is shown that a relation proposed some years ago by Hardegree, namely the Sasaki arrow , can be interpreted in terms of the relation ; furthermore, this new relation turns out to be successful also in cases in which the Sasaki arrow fails.     

Monte Carlo simulation of cluster kinetics in three-dimensional Ising model

Simulations in lattices of size 100³ at the critical point of the three-dimensional Glauber kinetic Ising model indicate that the 1935 Becker-Doring equation no longer works there: The growth rates decay in time. These conclusions confirm those in two dimensions.     

A-B site-bond correlated percolation for antiferromagnetic Ising model: The Bethe cluster approximation

We study the site-bond percolation problem for clusters of holes and particles with antiferromagnetic order by means of the Bethe cluster approximation. We find that the droplets (i.e.P _B =1?e ^?|K|/2) diverge at the antiferromagnetic critical pointH=0,T=T _c; however forH≠0 they diverge along a percolation line which is different from the Antiferromagnetic Phase Boundary except atT=0.     

Yang-Baxter and other relations for free-fermion and Ising models

Eight-vertex, free fermion, and Ising models are formulated using a convention that emphasizes the algebra of the local transition operators that arise in the quantum inverse method. Equivalent classes of models, are investigated, with particular emphasis on the role of the star-triangle relations. Using these results, a natural and symmetrical parametrization is introduced and Yang-Baxter relations are constructed in an elementary way. The paper concludes with a consideration of duality, which links the present work to a recent paper of Baxter on the free fermion model.