Singularity of the specific heat of two-dimensional random Ising models

The singularity of the specific heat is studied for the dilution (J>J'>0) type and Gaussian type random Ising models using the Pfaffian method numerically. The type of singularity at the paramagnetic-ferromagnetic phase boundary is studied using the standard regression method using data up to system size. It is shown that the logarithmic type singularity is more reliable than the double-logarithmic type and cusp type singularities. The critical temperatures are estimated accurately for both the dilution type and Gaussian type random Ising models. A phase diagram relating strength of the randomness and temperature is also presented. Received: 26 February 1998 / Revised: 15 May 1998 / Accepted: 25 June 1998     

Phase transition in the Blume-Capel model with second neighbour interaction

A two dimensional antiferromagnetic spin-1 Ising model with negative next- nearest neighbour interaction (J ₂ <0) and under an external magnetic field is investigated by two methods: The mean-field theory and Finite-Size-Scaling based on transfer matrix (TMFSS) calculations. The ground state diagrams exhibit several new phases including frustrated ones. At finite temperature we obtain by these two methods quite rich phase diagrams, with several multicritical points. While Mean field approximation yields phase diagrams which are sometimes even qualitatively incorrect, accurate results are obtained from transfer matrix finite size scaling calculations. For a certain range of interaction parameters, the model is shown to violate the ordinary universality hypothesis. Received: 3 November 1997 / Revised: 31 March 1998 / Accepted: 7 April 1998     

On the incommensurate phase in modulated Heisenberg chains

Using the density matrix renormalization group method (DMRG) we calculate the magnetization of frustrated S=1/2 Heisenberg chains for various modulation patterns of the nearest neighbour coupling: commensurate, incommensurate with sinusoidal modulation and incommensurate with solitonic modulation. We focus on the order of the phase transition from the commensurate dimerized phase (D) to the incommensurate phase (I). It is shown that the order of the phase transition depends sensitively on the model. For the solitonic model in particular, a k-dependent elastic energy modifies the order of the transition. Furthermore, we calculate gaps in the incommensurate phase in adiabatic approximation. Received: 9 March 1998 / Accepted: 17 April 1998     

Low temperature thermodynamics of inverse square spin models in one dimension

We present a field-theoretic renormalization group calculation in two loop order for classical O(N)-models with an inverse square interaction in the vicinity of their lower critical dimensionality one. The magnetic susceptibility at low temperatures is shown to diverge like with a=(N-2)/(N-1) and . From a comparison with the exactly solvable Haldane-Shastry model we find that the same temperature dependence applies also to ferromagnetic quantum spin chains. Received: 20 February 1998 / Revised: 27 April 1998 / Accepted: 30 April 1998     

Monte-Carlo study of correlations in quantum spin chains at non-zero temperature

Antiferromagnetic Heisenberg spin chains with various spin values (S=1/2,1,3/2,2,5/2) are studied numerically with the quantum Monte-Carlo method. Effective spin S chains are realized by ferromagnetically coupling n=2S antiferromagnetic spin chains with S=1/2. The temperature dependence of the uniform susceptibility, the staggered susceptibility, and the static structure factor peak intensity are computed down to very low temperatures, . The correlation length at each temperature is deduced from numerical measurements of the instantaneous spin-spin correlation function. At high temperatures, very good agreement with exact results for the classical spin chain is obtained independent of the value of S. For the S=2 chain which has a gap , the correlation length and the uniform susceptibility in the temperature range are well predicted by the semi-classical theory of Damle and Sachdev. Received: 23 December 1997 / Revised and Accepted: 11 March 1998     

On evaluation of transverse spin fluctuations in quantum magnets

A numerical method is described for evaluating transverse spin correlations in the random phase approximation. Quantum spin-fluctuation corrections to sublattice magnetization are evaluated for the antiferromagnetic ground state of the half-filled Hubbard model in two and three dimensions in the whole U/t range. Extension to the case of defects in the AF is also discussed for spin vacancies and low-U impurities. In the limit, the vacancy-induced enhancement in the spin fluctuation correction is obtained for the spin-vacancy problem in two dimensions, for vacancy concentration up to the percolation threshold. For low-U impurities, the overall spin fluctuation correction is found to be strongly suppressed, although surprisingly spin fluctuations are locally enhanced at the low-U sites. Received 27 April 1998 and Received in final form 13 August 1998     

Quantum rotors with regular frustration and the quantum Lifshitz point

We have discussed the zero-temperature quantum phase transition in n-component quantum rotor Hamiltonian in the presence of regular frustration in the interaction. The phase diagram consists of ferromagnetic, helical and quantum paramagnetic phase, where the ferro-para and the helical-para phase boundary meets at a multicritical point called a (d,m) quantum Lifshitz point where (d,m) indicates that the m of the d spatial dimensions incorporate frustration. We have studied the Hamiltonian in the vicinity of the quantum Lifshitz point in the spherical limit and also studied the renormalisation group flow behaviour using standard momentum space renormalisation technique (for finite n). In the spherical limit ()one finds that the helical phase does not exist in the presence of any nonvanishing quantum fluctuation for m =d though the quantum Lifshitz point exists for all d > 1+m/2, and the upper critical dimensionality is given by d _u = 3 +m/2. The scaling behaviour in the neighbourhood of a quantum Lifshitz point in d dimensions is consistent with the behaviour near the classical Lifshitz point in (d+z) dimensions. The dynamical exponent of the quantum Hamiltonian z is unity in the case of anisotropic Lifshitz point (d>m) whereas z=2 in the case of isotropic Lifshitz point (d=m). We have evaluated all the exponents using the renormalisation flow equations along-with the scaling relations near the quantum Lifshitz point. We have also obtained the exponents in the spherical limit (). It has also been shown that the exponents in the spherical model are all related to those of the corresponding Gaussian model by Fisher renormalisation. Received: 23 December 1997 / Received in final form: 6 January 1998 / Accepted: 7 January 1998     

Suppression of antiferromagnetic correlations by quenched dipole-type impurities

The effects of quenched dipole moments on a two-dimensional Heisenberg antiferromagnet are found exactly, by applying the renormalization group to the appropriate classical non-linear sigma model. Such dipole moments represent random fields with power law correlations. At low temperatures, they also represent the long range effects of quenched random strong ferromagnetic bonds on the antiferromagnetic correlation length, , of a two-dimensional Heisenberg antiferromagnet. It is found that the antiferromagnetic long range order is destroyed for any non-zero concentration, x, of the dipolar defects, even at zero temperature. Below a line , where T is the temperature, is independent of T, and decreases exponentially with x. At higher temperatures, it decays exponentially with , with an effective stiffness constant , which decreases with increasing x/T. The latter behavior is the same as for annealed dipole moments, and we use our quenched results to interpolate between the two types of averaging for the problem of ferromagnetic bonds in an antiferromagnet. The results are used to estimate the three-dimensional Néel temperature of a lamellar system with weakly coupled planes, which decays linearly with x at small concentrations, and drops precipitously at a critical concentration. These predictions are shown to reproduce successfully several of the prominent features of experiments on slightly doped copper oxides. Received 22 October 1998     

Interaction versus dimerization in one-dimensional Fermi systems

In order to study the effect of interaction and lattice distortion on quantum coherence in one-dimensional Fermi systems, we calculate the ground state energy and the phase sensitivity of a ring of interacting spinless fermions on a dimerized lattice. Our numerical DMRG studies, in which we keep up to 1000 states for systems of about 100 sites, are supplemented by analytical considerations using bosonization techniques. We find a delocalized phase for an attractive interaction, which differs from that obtained for random lattice distortions. The extension of this delocalized phase depends strongly on the dimerization induced modification of the interaction. Taking into account the harmonic lattice energy, we find a dimerized ground state for a repulsive interaction only. The dimerization is suppressed at half filling, when the correlation gap becomes large. Received: 11 February 1998 / Revised: 1st April 1998 / Accepted: 30 April 1998     

A model for anomalous directed percolation

We introduce a model for the spreading of epidemics by long-range infections and investigate the critical behaviour at the spreading transition. The model generalizes directed bond percolation and is characterized by a probability distribution for long-range infections which decays in d spatial dimensions as . Extensive numerical simulations are performed in order to determine the density exponent and the correlation length exponents and for various values of . We observe that these exponents vary continuously with , in agreement with recent field-theoretic predictions. We also study a model for pairwise annihilation of particles with algebraically distributed long-range interactions. Received: 4 September 1998 / Accepted: 22 September 1998     

A geometric generalization of field theory to manifolds of arbitrary dimension

We introduce a generalization of the O(N) field theory to N-colored membranes of arbitrary inner dimension D. The O(N) model is obtained for , while leads to self-avoiding tethered membranes (as the O(N) model reduces to self-avoiding polymers). The model is studied perturbatively by a 1-loop renormalization group analysis, and exactly as .Freedom to choose the expansion point D, leads to precise estimates of critical exponents of the O(N) model. Insights gained from this generalization include a conjecture on the nature of droplets dominating the 3d-Ising model at criticality; and the fixed point governing the random bond Ising model. Received: 15 October 1998 / Accepted: 4 November 1998     

Quadrupolar order in the Heisenberg ferromagnet with the single-ion cubic anisotropy

The ground state properties of S =2 ferromagnets with isotropic Heisenberg exchange (J) and single-ion cubic anisotropy (D) are studied. The perturbation theory for is used to find an effective Hamiltonian up to the fourth order for 1, 2 and 3 dimensions. It is shown that in opposition to the MFA prediction there is the quadrupolar long range order at T = 0 in the non-magnetic state of the system without a quadrupolar type of interaction. The effect is a consequence of the quantum nature of the model. Received: 19 February 1998 / Accepted: 17 March 1998     

A DMRG study of the -symmetric Heisenberg chain

The spin one-half Heisenberg chain with U q [ SU (2)] symmetry is studied via density-matrix renormalization. Ground-state energy and q-symmetric correlation functions are calculated for the non-Hermitian case with integer r. This gives bulk and surface exponents for (para)fermionic correlations in the related Ising and Potts models. The case of real q corresponding to a diffusion problem is treated analytically. Received: 18 February 1998 / Accepted: 17 March 1998     

Phase diagram of randomly polymerized membrane

Using a replica formalism, a generalization of a recent mean field model corresponding to the observed wrinkling transition in randomly polymerized membranes is presented. In this model we study the effects of global fluctuations of the surface normals to the flat membrane, which can be introduced by a random local field. In absence of these global fluctuations, we show that, the model exhibits both continuous and discontinuous transitions between flat and wrinkled phases, contrary to what has been predicted by Bensimon et al. and Attal et al. Phase diagrams both in replica symmetry and in breaking of replica symmetry in sense of Almeida and Thouless are given. We have also investigated the effects of global fluctuations on the replica symmetry phase diagram. We show that, the wrinkled phase is favored and the flat phase is unstable. For large global fluctuations, the transition between wrinkled and flat phases becomes first order. Received: 3 December 1997 / Revised: 31 March 1998 / Accepted: 3 August 1998     

High-temperature series analysis of the p-state Potts glass model on d-dimensional hypercubic lattices

We analyze recently extended high-temperature series expansions for the “Edwards-Anderson” spin-glass susceptibility of the p-state Potts glass model on d-dimensional hypercubic lattices for the case of a symmetric bimodal distribution of ferro- and antiferromagnetic nearest-neighbor couplings . In these star-graph expansions up to order 22 in the inverse temperature , the number of Potts states p and the dimension d are kept as free parameters which can take any value. By applying several series analysis techniques to the new series expansions, this enabled us to determine the critical coupling K_c and the critical exponent of the spin-glass susceptibility in a large region of the two-dimensional (p,d)-parameter space. We discuss the thus obtained information with emphasis on the lower and upper critical dimensions of the model and present a careful comparison with previous estimates for special values of p and d. Received: 25 May 1998 / Revised and Accepted: 11 August 1998     

Theory of site-disordered magnets

In realistic spinglasses, such as , and , magnetic atoms are located at random positions. Their couplings are determined by their relative positions. For such systems a field theory is formulated. In certain limits it reduces to the Hopfield model, the Sherrington-Kirkpatrick model, and the Viana-Bray model. The model has a percolation transition, while for RKKY couplings the “concentration scaling” occurs. Within the Gaussian approximation the Ginzburg-Landau expansion is considered in the clusterglass phase, that is to say, for not too small concentrations. Near special points, the prefactor of the cubic term, or the one of the replica-symmetry-breaking quartic term, may go through zero. Around such points new spin glass phases are found. Received: 27 April 1998 / Received in final form: 27 July 1998 / Accepted: 13 August 1998     

Numerical diagonalization analysis of the ground-state superfluid-localization transition in two dimensions

Ground state of the two-dimensional hard-core-boson system in the presence of the quenched random chemical potential is investigated by means of the exact-diagonalization method for the system sizes up to L=5. The criticality and the DC conductivity at the superfluid-localization transition have been controversial so far. We estimate, with the finite-size scaling analysis, the correlation-length and the dynamical critical exponents as and z=2, respectively. The AC conductivity is computed with the Gagliano-Balseiro formula, with which the resolvent (dynamical response function) is expressed in terms of the continued-fraction form consisting of Lanczos tri-diagonal elements. Thereby, we estimate the universal DC conductivity as . Received 19 August 1998     

Free energy of semiflexible polymers and structure of interfaces

The free energy of semiflexible polymers is calculated as a functional of the compositional scalar order parameter and the orientational order parameter of second-rank tensor S_ij on the basis of a microscopic model of wormlike chains with variable segment lengths. We use a density functional theory and a gradient expansion to evaluate the entropic part of the free energy, which is given in a power series of .The interaction term of the free energy is derived with a random phase approximation. For the rigid rod limit, the nematic-isotropic transition point is given by , N and w being the degree of polymerization and the anisotropic interaction parameter, respectively, and the degree of ordering at the transition point is 0.33448. We also find that the contour length of polymer chains becomes larger in a nematic phase than in an isotropic phase. Interface profiles are obtained numerically for some typical cases. In the neighborhood of isotropic-isotropic interfaces, polymer chains tend to align parallel to the interface on the polymer-rich side and perpendicular on the poor side. When an isotropic region and a nematic region coexist, orientational order parallel to the interface is preferred in the nematic region. Received: 28 May 1998 / Revised: 12 August 1998 / Accepted: 8 September 1998