Magnetization densities as replica parameters: The dilute ferromagnet

In this paper we compute exactly the ground state energy and entropy of the dilute ferromagnetic Ising model. The two thermodynamic quantities are also computed when a magnetic field with random locations is present. The result is reached in the replica approach frame by a class of replica order parameters introduced by Monasson (1998) [5]. The strategy is first illustrated considering the SK model, for which we will show the complete equivalence with the standard replica approach. Then, we apply to the diluted ferromagnetic Ising model with a random located magnetic field, which is mapped into a Potts model.     

Domain growth in the three-dimensional dilute Ising model

We investigate the kinetics of domain growth in the three-dimensional Ising model with quenched random site dilution, using Monte Carlo simulation technique. A crossover from the power law growth regime to a much slower growth observed in our simulation is interpreted through the roughening of the interfaces by the quenched impurities. The results are also compared with the corresponding results in two dimensions.     

Exact results for a dilute Potts model

Exact results are obtained for the annealed, dilute,q-component Potts model on the decorated square lattice. The phase diagram is found to consist of a high-temperature region, a low-temperature region, and a two-phase region in between which arises only forq>4: exact expressions for the phase boundary and the critical probability are derived. At the critical point the specific heat is generally finite and has a cusp; the slope of the cusp is finite forq=4 and infinite (vertical) forq=2 and 3.Work supported in part by NSF Grant No. DMR 78-18808.     

Monte Carlo simulation of three-dimensional dilute Ising model

A multispin coding program for site-diluted Ising models on large simple cubic lattices is described in detail. The spontaneous magnetization is computed as a function of temperature, and the critical temperature as a function of concentration is found to agree well with the data of Marro et al.⁽⁴⁾ and Landau⁽³⁾ for smaller systems.The first successful epsilon expansion seems to be by D. E. Khmelnitskii,ZhETF 68:1960 (1975), English translationSov. Phys. JETP 41:981 (1975); for numerical estimates see K. E. Newman and E. K. Riedel,Phys. Rev. H25:264 (1982), for experiments see R. J. Birgenau, R. A. Cowley, G. Shirane and H. Yoshizawa,J. Stat. Phys. 34:817 (1984).     

A self-consistent Ornstein-Zernike approximation for the site-diluted Ising model

We propose a theory for the site-diluted Ising model which is an extension to disordered systems of the self-consistent Ornstein-Zernike approximation of Hoye and Stell. By using the replica method in the context of liquid-state theory, we treat the concentration of impurities as an ordinary thermodynamic variable. This approach is not limited to the weak-disorder regime or to the vicinity of the percolation point. A preliminary analysis using series expansion shows that it can predict accurately the dependence of the critical temperature on dilution and can reproduce the nonuniversal behavior of the effective exponents. The theory also gives a reasonable estimate of the percolation threshold.     

Spontaneous magnetization of randomly dilute ferromagnets

We consider Ising ferromagnets on random subgraphs of the square lattice. These are obtained by independent random selections either of sites or of bonds. We assume that for each site (or, respectively, bond) the probability of being selected exceeds the critical percolation probability. Then, at sufficiently low temperatures and zero external field, spontaneous magnetization occurs. Some further related results are obtained.     

Nonequilibrium kinetics: Exact and approximate solutions

The relaxation of an internal state distribution in the presence of an excess of an inert gas is considered. The explicit time dependence of the nonequilibrium contributions to the transition rate coefficients is approximated using the Kapral-Hudson-Ross method. The resulting solution contains cross-correlation terms which do not appear when a single reaction is considered. It is shown that the first term of a perturbation expansion of an exact formal solution gives the Kapral-Hudson-Ross solution for short times, and the Chapman-Enskog solution at long times if there is a wide separation in time scales. The Kapral-Hudson-Ross, Chapman-Enskog, and exact solutions are compared for a two-state, hard-sphere model system.Based on a dissertation by S. Hudson, Massachusetts Institute of Technology, 1972.National Science Foundation predoctoral fellow, 1967–1971.     

Exact solutions for correlation functions of an Ising nanocluster

Using a four-spin nanocluster as an example, a technique for developing exact expressions for correlation functions of Ising-type nanosystems is demonstrated.     

Exact and approximate results for the Anderson impurity model

This paper provides exact and rigorous upper and lower bounds to the ground state energy of the Single Impurity Anderson Model in the limit of infiniteU. The upper bound approximates the ground state energy very well and the corresponding state may be used as a starting point for further investigations. The energy spectrum of the Single Impurity Anderson Model is calculated exactly for a special non-trivial model density of states describing the bare band electrons. All approaches introduced here are compared to this exact result and to other ground state energy calculations.     

Exact and approximate solutions for the electron nonsinglet structure function in QED

Two exact, valid up to infinite perturbative order, numerical solutions of the Lipatov equation for the nonsinglet electron structure function in the QED are presented. One of them is of the Monte Carlo type and another is based on the numerical inversion of the Mellin transform. They agree numerically to a very high precision (better than 0.05%). Within the leading logarithmic approximation, the exact solution is compared with the perturbative second and third order exponentiated solutions. It is shown that the perturbative second order solution inspired by the Yennie-Frautschi-Suura (exclusive) exponentiation is much closer to the exact solution than the other ones. New compact analytical formula for the third order exponentiated solution is given. It is shown to be in perfect numerical agreement with the infinite order solution of the Monte Carlo and Mellin type.     

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.     

Exact results for the diluted Ising model with competing interactions

We consider the random-bond Ising model with the exchange integrals J > 0, ?J and 0 with the respective probabilities p, q and r, where p + q + r = 1. We give the exact value of the averaged internal energy and an exact upper bound to the averaged specific heat at temperature T determined by $\text{k}{}_{\mathrm{B}}\text{T =}\text{2J}\text{In[}\text{p}\text{(1 ? p ? r)}\text{]}\text{,}$ where k_B is the Boltsmann constant. We show that all the averaged correlation functions of even spins are non-negative at this temperature.     

Effective spin-glass Hamiltonian for the anomalous dynamics of the HMF model

We discuss an effective spin-glass Hamiltonian which can be used to study the glassy-like dynamics observed in the metastable states of the Hamiltonian mean field (HMF) model. By means of the Replica formalism, we were able to find a self-consistent equation for the glassy order parameter which reproduces, in a restricted energy region below the phase transition, the microcanonical simulations for the polarization order parameter recently introduced in the HMF model.     

Hydrodynamic interactions of dilute polymer solutions in elongational flow

Dumbbell models are only crude representations of actual polymer molecules, but their simplicity allows for explicit calculations which in many instances have shed light on the connection between molecular properties and rheological behavior. On the other hand, hydrodynamic interactions are known to strongly influence the dynamic response of polymer solutions and this makes the representation of the hydrodynamic drag an important aspect in the calculations. In the present work, the effects of the state of flow are incorporated explicitly in the frictional properties of the FENE model of a dilute polymer solution. Results for the steady elongational viscosity and the mean square end-to-end distance arc presented.     

Exact and approximate results of non-extensive quantum statistics

We develop an analytical technique to derive explicit forms of thermodynamical quantities within the asymptotic approach to non-extensive quantum distribution functions. Using it, we find an expression for the number of particles in a boson system which we compare with other approximate scheme (i.e. factorization approach), and with the recently obtained exact result. To do this, we investigate the predictions on Bose-Einstein condensation and the blackbody radiation. We find that both approximation techniques give results similar to (up to ) the exact ones, making them a useful tool for computations. Because of the simplicity of the factorization approach formulae, it appears that this is the easiest way to handle with physical systems which might exhibit slight deviations from extensivity. Received 19 August 1999 and Received in final form 1 November 1999     

Nonconcavity of the magnetization in Ising ferromagnets

Some Ising ferromagnets having nonconcave magnetization are presented as counterexamples to the often assumed case of concave magnetization.     

Rigorous mean-field model for coherent-potential approximation: Anderson model with free random variables

A model of a randomly disordered system with site-diagonal random energy fluctuations is introduced. It is an extension of Wegner'sn-orbital model to arbitrary eigenvalue distribution in the electronic level space. The new feature is that the random energy values are not assumed to be independent at different sites, but free. Freeness of random variables is an analog of the concept of independence for noncommuting random operators. A possible realization is the ensemble of randomly rotated matrices at different lattice sites. The one- and two-particle Green functions of the proposed Hamiltonian are calculated exactly. The eigenstates are extended and the conductivity is novanishing everywhere inside the band. The long-range behavior and the zero-frequency limit of the two-particle Green function are universal with respect to the eigenvalue distribution in the electronic level space. The solutions solve the CPA equation for the one- and two-particle Green function of the corresponding Anderson model. Thus our (multisite) model is a rigorous mean-field model for the (single-site) CPA. We show how the Lloyd model is included in our model and treat various kinds of noises.     

Exact and approximate solutions of Riemann problems in non-linear elasticity

Eulerian shock-capturing schemes have advantages for modelling problems involving complex non-linear wave structures and large deformations in solid media. Various numerical methods now exist for solving hyperbolic conservation laws that have yet to be applied to non-linear elastic theory. In this paper one such class of solver is examined based upon characteristic tracing in conjunction with high-order monotonicity preserving weighted essentially non-oscillatory (MPWENO) reconstruction. Furthermore, a new iterative method for finding exact solutions of the Riemann problem in non-linear elasticity is presented. Access to exact solutions enables an assessment of the performance of the numerical techniques with focus on the resolution of the seven wave structure. The governing model represents a special case of a more general theory describing additional physics such as material plasticity. The numerical scheme therefore provides a firm basis for extension to simulate more complex physical phenomena. Comparison of exact and numerical solutions of one-dimensional initial values problems involving three-dimensional deformations is presented.     

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.