Numerical evidence for Goldstone-mode effects in the three-dimensional O(4)-model

We investigate the three-dimensional O(4)-model on 24^3--96^3 lattices as a function of the magnetic field H. Below the critical point, in the broken phase, we confirm explicitly the H^{1/2} dependence of the magnetization and the corresponding H^{-1/2} divergence of the longitudinal susceptibility, which are predicted from the existence of Goldstone modes. At the critical point the magnetization follows the expected finite-size-scaling behavior, with critical exponent delta = 4.87(1).Comment: LATTICE99 (spin models), added preprint number and corrected acknowledgement

Dewetting hydrodynamics in 1+1 dimensions

A model for the phase transition between partial wetting and dewetting of a substrate has been formulated that explicitly incorporates the hydrodynamic flow during the dewetting process in 1 + 1 dimensions. The model simulates a fluid layer of finite thickness on a substrate in coexistence with a dry part of the substrate and a gas phase above the substrate. Under nonequilibrium ''dewetting'' conditions, the front between the dry part and the wet part of the surface moves towards the wet part inducing hydrodynamic flow inside the wet layer. In more general terms, the model handles two immiscible fluids with a freely movable interface in an inhomogeneous external force field. Handling the interface by a new variant of the phase-field model, we obtain an efficient code with well-defined interfacial properties. In particular, the (free) energy can be chosen at will. We demonstrate that our model works well in the viscosity range of creeping flow and we give qualitative results for the higher Reynolds numbers. Connections to experimental realizations are discussed

Influence of external flows on crystal growth: numerical investigation

We use a combined phase-field/lattice-Boltzmann scheme [D. Medvedev, K. Kassner, Phys. Rev. E {\bf 72}, 056703 (2005)] to simulate non-facetted crystal growth from an undercooled melt in external flows. Selected growth parameters are determined numerically. For growth patterns at moderate to high undercooling and relatively large anisotropy, the values of the tip radius and selection parameter plotted as a function of the Peclet number fall approximately on single curves. Hence, it may be argued that a parallel flow changes the selected tip radius and growth velocity solely by modifying (increasing) the Peclet number. This has interesting implications for the availability of current selection theories as predictors of growth characteristics under flow. At smaller anisotropy, a modification of the morphology diagram in the plane undercooling versus anisotropy is observed. The transition line from dendrites to doublons is shifted in favour of dendritic patterns, which become faster than doublons as the flow speed is increased, thus rendering the basin of attraction of dendritic structures larger. For small anisotropy and Prandtl number, we find oscillations of the tip velocity in the presence of flow. On increasing the fluid viscosity or decreasing the flow velocity, we observe a reduction in the amplitude of these oscillations.Comment: 10 pages, 7 figures, accepted for Physical Review E; size of some images had to be substantially reduced in comparison to original, resulting in low qualit

Phase Field Modeling of Fast Crack Propagation

We present a continuum theory which predicts the steady state propagation of cracks. The theory overcomes the usual problem of a finite time cusp singularity of the Grinfeld instability by the inclusion of elastodynamic effects which restore selection of the steady state tip radius and velocity. We developed a phase field model for elastically induced phase transitions; in the limit of small or vanishing elastic coefficients in the new phase, fracture can be studied. The simulations confirm analytical predictions for fast crack propagation.Comment: 5 pages, 11 figure

Percolation and Critical Behaviour in SU(2) Gauge Theory

The paramagnetic-ferromagnetic transition in the Ising model can be described as percolation of suitably defined clusters. We have tried to extend such picture to the confinement-deconfinement transition of SU(2) pure gauge theory, which is in the same universality class of the Ising model. The cluster definition is derived by approximating SU(2) by means of Ising-like effective theories. The geometrical transition of such clusters turns out to describe successfully the thermal counterpart for two different lattice regularizations of (3+1)-d SU(2).Comment: Lattice 2000 (Finite Temperature), 4 pages, 4 figures, 2 table

A polarizable interatomic force field for TiO2_2 parameterized using density functional theory

We report a classical interatomic force field for TiO$_2$, which has been parameterized using density functional theory forces, energies, and stresses in the rutile crystal structure. The reliability of this new classical potential is tested by evaluating the structural properties, equation of state, phonon properties, thermal expansion, and some thermodynamic quantities such as entropy, free energy, and specific heat under constant volume. The good agreement of our results with {\em ab initio} calculations and with experimental data, indicates that our force-field describes the atomic interactions of TiO$_2$ in the rutile structure very well. The force field can also describe the structures of the brookite and anatase crystals with good accuracy.Comment: Accepted for publication in Phys. Rev. B; Changes from v1 include multiple minor revisions and a re-write of the description of the force field in Section II

Length Scales and Power Laws in the Two-Dimensional Forest-Fire Model

We re-examine a two-dimensional forest-fire model via Monte-Carlo simulations and show the existence of two length scales with different critical exponents associated with clusters and with the usual two-point correlation function of trees. We check resp. improve previously obtained values for other critical exponents and perform a first investigation of the critical behaviour of the slowest relaxational mode. We also investigate the possibility of describing the critical point in terms of a distribution of the global density. We find that some qualitative features such as a temporal oscillation and a power law of the cluster-size distribution can nicely be obtained from such a model that discards the spatial structure.Comment: 20 pages plain TeX, 7 figures included using psfig.sty, PostScript for the complete paper also available at http://www.physik.fu-berlin.de/~ag-peschel/papers/forest2d.ps.gz , extra software at http://www.physik.fu-berlin.de/~ag-peschel/software/forest2d.html ; main change: inclusion of further data in the determination of nu_T in Section 2.1 + some small changes; final version to appear in Physica

Comments on Sweeny and Gliozzi dynamics for simulations of Potts models in the Fortuin-Kasteleyn representation

We compare the correlation times of the Sweeny and Gliozzi dynamics for two-dimensional Ising and three-state Potts models, and the three-dimensional Ising model for the simulations in the percolation prepresentation. The results are also compared with Swendsen-Wang and Wolff cluster dynamics. It is found that Sweeny and Gliozzi dynamics have essentially the same dynamical critical behavior. Contrary to Gliozzi's claim (cond-mat/0201285), the Gliozzi dynamics has critical slowing down comparable to that of other cluster methods. For the two-dimensional Ising model, both Sweeny and Gliozzi dynamics give good fits to logarithmic size dependences; for two-dimensional three-state Potts model, their dynamical critical exponent z is 0.49(1); the three-dimensional Ising model has z = 0.37(2).Comment: RevTeX, 4 pages, 5 figure

Simulations of a single membrane between two walls using a Monte Carlo method

Quantitative theory of interbilayer interactions is essential to interpret x-ray scattering data and to elucidate these interactions for biologically relevant systems. For this purpose Monte Carlo simulations have been performed to obtain pressure P and positional fluctuations sigma. A new method, called Fourier Monte-Carlo (FMC), that is based on a Fourier representation of the displacement field, is developed and its superiority over the standard method is demonstrated. The FMC method is applied to simulating a single membrane between two hard walls, which models a stack of lipid bilayer membranes with non-harmonic interactions. Finite size scaling is demonstrated and used to obtain accurate values for P and sigma in the limit of a large continuous membrane. The results are compared with perturbation theory approximations, and numerical differences are found in the non-harmonic case. Therefore, the FMC method, rather than the approximations, should be used for establishing the connection between model potentials and observable quantities, as well as for pure modeling purposes.Comment: 10 pages, 10 figure

Correlated percolation and the correlated resistor network

We present some exact results on percolation properties of the Ising model, when the range of the percolating bonds is larger than nearest-neighbors. We show that for a percolation range to next-nearest neighbors the percolation threshold Tp is still equal to the Ising critical temperature Tc, and present the phase diagram for this type of percolation. In addition, we present Monte Carlo calculations of the finite size behavior of the correlated resistor network defined on the Ising model. The thermal exponent t of the conductivity that follows from it is found to be t = 0.2000 +- 0.0007. We observe no corrections to scaling in its finite size behavior.Comment: 16 pages, REVTeX, 6 figures include