453 research outputs found
Tumbling of a rigid rod in a shear flow
The tumbling of a rigid rod in a shear flow is analyzed in the high viscosity
limit. Following Burgers, the Master Equation is derived for the probability
distribution of the orientation of the rod. The equation contains one
dimensionless number, the Weissenberg number, which is the ratio of the shear
rate and the orientational diffusion constant. The equation is solved for the
stationary state distribution for arbitrary Weissenberg numbers, in particular
for the limit of high Weissenberg numbers. The stationary state gives an
interesting flow pattern for the orientation of the rod, showing the interplay
between flow due to the driving shear force and diffusion due to the random
thermal forces of the fluid. The average tumbling time and tumbling frequency
are calculated as a function of the Weissenberg number. A simple cross-over
function is proposed which covers the whole regime from small to large
Weissenberg numbers.Comment: 22 pages, 9 figure
Scaling in the vicinity of the four-state Potts fixed point
We study a self-dual generalization of the Baxter-Wu model, employing results
obtained by transfer matrix calculations of the magnetic scaling dimension and
the free energy. While the pure critical Baxter-Wu model displays the critical
behavior of the four-state Potts fixed point in two dimensions, in the sense
that logarithmic corrections are absent, the introduction of different
couplings in the up- and down triangles moves the model away from this fixed
point, so that logarithmic corrections appear. Real couplings move the model
into the first-order range, away from the behavior displayed by the
nearest-neighbor, four-state Potts model. We also use complex couplings, which
bring the model in the opposite direction characterized by the same type of
logarithmic corrections as present in the four-state Potts model. Our
finite-size analysis confirms in detail the existing renormalization theory
describing the immediate vicinity of the four-state Potts fixed point.Comment: 19 pages, 7 figure
Extended surface disorder in the quantum Ising chain
We consider random extended surface perturbations in the transverse field
Ising model decaying as a power of the distance from the surface towards a pure
bulk system. The decay may be linked either to the evolution of the couplings
or to their probabilities. Using scaling arguments, we develop a
relevance-irrelevance criterion for such perturbations. We study the
probability distribution of the surface magnetization, its average and typical
critical behaviour for marginal and relevant perturbations. According to
analytical results, the surface magnetization follows a log-normal distribution
and both the average and typical critical behaviours are characterized by
power-law singularities with continuously varying exponents in the marginal
case and essential singularities in the relevant case. For enhanced average
local couplings, the transition becomes first order with a nonvanishing
critical surface magnetization. This occurs above a positive threshold value of
the perturbation amplitude in the marginal case.Comment: 15 pages, 10 figures, Plain TeX. J. Phys. A (accepted
Effective Field Theory of the Zero-Temperature Triangular-Lattice Antiferromagnet: A Monte Carlo Study
Using a Monte Carlo coarse-graining technique introduced by Binder et al., we
have explicitly constructed the continuum field theory for the zero-temperature
triangular Ising antiferromagnet. We verify the conjecture that this is a
gaussian theory of the height variable in the interface representation of the
spin model. We also measure the height-height correlation function and deduce
the stiffness constant. In addition, we investigate the nature of defect-defect
interactions at finite temperatures, and find that the two-dimensional Coulomb
gas scenario applies at low temperatures.Comment: 26 pages, 9 figure
Monte Carlo Renormalization of the 3-D Ising model: Analyticity and Convergence
We review the assumptions on which the Monte Carlo renormalization technique
is based, in particular the analyticity of the block spin transformations. On
this basis, we select an optimized Kadanoff blocking rule in combination with
the simulation of a d=3 Ising model with reduced corrections to scaling. This
is achieved by including interactions with second and third neighbors. As a
consequence of the improved analyticity properties, this Monte Carlo
renormalization method yields a fast convergence and a high accuracy. The
results for the critical exponents are y_H=2.481(1) and y_T=1.585(3).Comment: RevTeX, 4 PostScript file
d=2 transverse-field Ising model under the screw-boundary condition: An optimization of the screw pitch
A length-N spin chain with the \sqrt{N}(=v)-th neighbor interaction is
identical to a two-dimensional (d=2) model under the screw-boundary (SB)
condition. The SB condition provides a flexible scheme to construct a d\ge2
cluster from an arbitrary number of spins; the numerical diagonalization
combined with the SB condition admits a potential applicability to a class of
systems intractable with the quantum Monte Carlo method due to the
negative-sign problem. However, the simulation results suffer from
characteristic finite-size corrections inherent in SB. In order to suppress
these corrections, we adjust the screw pitch v(N) so as to minimize the
excitation gap for each N. This idea is adapted to the transverse-field Ising
model on the triangular lattice with N\le32 spins. As a demonstration, the
correlation-length critical exponent is analyzed in some detail
Gap of the Linear Spin-1 Heisenberg Antiferromagnet: A Monte Carlo Calculation
We have performed Monte Carlo calculations of the energies of several low-lying energy states of one-dimensional, spin-1 Heisenberg antiferromagnets with linear sizes up to n=32. Our results support Haldane’s prediction that a gap exists in the excitation spectrum for n→∞. .A
Monte Carlo Calculation of Free Energy, Critical Point, and Surface Critical Behavior of Three-Dimensional Heisenberg Ferromagnets
A transfer-matrix Monte Carlo technique is developed to compute the free energy of three-dimensional, classical Heisenberg ferromagnets. From the free energy of systems with periodic and antiperiodic boundary conditions, helicity moduli are calculated. From these the critical couplings for simple-cubic (sc) and face-centered-cubic lattices are estimated by use of finite-size scaling. For the simple-cubic lattice, the critical dimension of the surface magnetization is estimated with standard Monte Carlo methods, yielding a result in excellent agreement with the ε-expansion work of Diehl and Nüsser
The Magnetization of the 3D Ising Model
We present highly accurate Monte Carlo results for simple cubic Ising
lattices containing up to spins. These results were obtained by means
of the Cluster Processor, a newly built special-purpose computer for the Wolff
cluster simulation of the 3D Ising model. We find that the magnetization
is perfectly described by , where
, in a wide temperature range .
If there exist corrections to scaling with higher powers of , they are very
small. The magnetization exponent is determined as (6). An
analysis of the magnetization distribution near criticality yields a new
determination of the critical point: ,
with a standard deviation of .Comment: 7 pages, 5 Postscript figure
Geometric properties of two-dimensional O(n) loop configurations
We study the fractal geometry of O() loop configurations in two dimensions
by means of scaling and a Monte Carlo method, and compare the results with
predictions based on the Coulomb gas technique. The Monte Carlo algorithm is
applicable to models with noninteger and uses local updates. Although these
updates typically lead to nonlocal modifications of loop connectivities, the
number of operations required per update is only of order one. The Monte Carlo
algorithm is applied to the O() model for several values of , including
noninteger ones. We thus determine scaling exponents that describe the fractal
nature of O() loops at criticality. The results of the numerical analysis
agree with the theoretical predictions.Comment: 18 pages, 6 figure
- …