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).     

Infrared asymptotic freedom of a hierarchical φ 3 6 lattice theory

For the weakly coupled lattice ₃ ⁶ theory in a hierarchical model approximation a nonperturbative renormalization group analysis in the spirit of Gawedzki and Kupiainen is performed to study the flow of the effective actions. We deduce a domain of attraction to the tricritical (Gaussian) fixed point. The two relevant coupling constants of the problem are controlled by analytic continuation to complex domains, tracing their images under the renormalization group iterations.     

Nucleon octet in a relativistic logarithmic potential

A simple independent-quark-model based on the Dirac equation with logarithmic potential is used to calculate several properties of octet baryons such as magnetic moment, the axial vector coupling constantg _A (n) for neutronβ-decay and the charge radius of the proton. In view of the simplicity of the model, the results obtained are quite good.     

Return to equilibrium in theXY model

We prove that the locally perturbedXY model returns to equilibrium under the unperturbed evolution but the unperturbed model does not necessarily approach equilibrium under the perturbed evolution. In fact this latter property is false for perturbation by a local magnetization. The failure is directly attributable to the formation of bound states. If the perturbation is quadratic these problems are reduced to spectral analysis of the one-particle Hamiltonian. We demonstrate that the perturbed Hamiltonian has a finite set of eigenvalues of finite multiplicity together with some absolutely continuous spectrum. Eigenvalues can occur in the continuum if, and only if, the perturbation dislocates the system. Singular continuous spectrum cannot occur.     

Stationary velocity distribution in an external field: A one-dimensional model

The velocity distribution of a charged hard rod coupled to an external field and moving in a neutral equilibrium hard rod gas is studied on the basis of Boltzmann's equation. The exact stationary solution is found. Above a threshold value the field becomes effective in the high-velocity region slowing down the decay of the velocity distribution. The drift velocity and the mean kinetic energy are discussed as functions of the field.     

高温等离子体中固体弹丸消融过程的研究

本文求得了弹丸消融过程的自洽解析解。结果表明:Parks的弹丸消融理论需要修正;弹丸消融速率由消融物云和表面蒸发层共同控制;入射电子通过消融物云时慢化效应并不显著;消融物离开弹丸表面时一般呈部分电离态;聚变堆条件下,弹丸表面可能会出现激波。     

Asymptotic behavior of fluctuations for the 1D Ising model in zero-temperature limit

Fluctuation of the average spin for one-dimensional Ising spins with nearest neighbor interactions are studied. The distribution function for the average spin is calculated for a finite volume, finite temperature, and finite magnetic field. As the volume increases and the temperature diminishes at zero magnetic field, there are two limits in which the probability distribution shows quite different behaviors: in the thermodynamic limit as the volume goes to infinity for finite temperature, small deviations of the fluctuations are described by a Gaussian distribution, and in the limit as the temperature vanishes for a finite volume, the ground states are realized with probability one. The crossover between these limits is analyzed via a ratio of the correlation length to the volume. The helix-coil transition in a polypeptide is discussed as an application.     

Fast vectorized algorithm for the Monte Carlo simulation of the random field Ising model

An algorithm for the simulation of the 3-dimensional random field Ising model with a binary distribution of the random fields is presented. It uses multi-spin coding and simulates 64 physically different systems simultaneously. On one processor of a Cray YMP it reaches a speed of 184 million spin updates per second. For smaller field strength we present a version of the algorithm that can perform 242 million spin updates per second on the same machine.     

Droplet dynamics for asymmetric Ising model

Nucleation from a metastable state is studied for an anisotropic Ising model at very low temperatures. It turns out that the critical nucleus as well as configurations on a typical path to it differ from the Wulff shape of an equilibrium droplet.     

Ising model in a quasiperiodic transverse field,percolation, and contact processes in quasiperiodic environments

Quantum Ising models in a transverse field are related to continuous-time percolation processes whose oriented percolation versions are contact processes. We study such models in the presence of quasiperiodic disorder and prove localization in the ground state, no percolation, and extinction, respectively, for sufficiently large disorder.