Monte Carlo investigation of phase transitions in the frustrated Heisenberg model on a triangular lattice

The critical properties and phase transitions of the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice have been investigated using the Monte Carlo method with a replica algorithm. The critical temperature has been determined and the character of the phase transitions has been analyzed using the method of fourth-order Binder cumulants. A second-order phase transition has been found in the three-dimensional frustrated Heisenberg model on a triangular lattice. The static magnetic and chiral critical exponents of the heat capacity α, the susceptibility γ and γ _k , the magnetization β and β _k , the correlation length ν and ν _k , as well as the Fisher exponents η and η _k , have been calculated in terms of the finite-size scaling theory. It has been demonstrated that the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice forms a new universality class of the critical behavior.     

Computer simulation of the frustrated antiferromagnetic Heisenberg model on a layered triangle lattice

Critical properties of the 3D frustrated Heisenberg model on a triangle latticeare investigated using a replica Monte-Carlo method that considers the interaction between next nearest neighbors. Static magnetic and chiral critical indices for heat capacity α, susceptibility γ, γ _k , magnetization β, β _k , and correlation radius ν are calculated using the theory of finite-size scaling.     

A frustrated heisenberg antiferromagnet on a triangular lattice with next-nearest neighbor interactions

Using the Monte Carlo method, we study the critical properties of the three-dimensional frustrated Heisenberg model on a triangular lattice with allowance for next-nearest neighbor interactions. Using the theory of finite-size scaling, we calculate the static magnetic and chiral critical exponents of heat capacity α, susceptibility γ, γ _k , magnetization β, β _k , and correlation length ν.     

Towards Cohomology of Renormalization: Bigrading the Combinatorial Hopf Algebra of Rooted Trees

The renormalization of quantum field theory twists the antipode of a noncocommutative Hopf algebra of rooted trees, decorated by an infinite set of primitive divergences. The Hopf algebra of undecorated rooted trees, ℋ _R , generated by a single primitive divergence, solves a universal problem in Hochschild cohomology. It has two nontrivial closed Hopf subalgebras: the cocommutative subalgebra ℋ_ladder of pure ladder diagrams and the Connes–Moscovici noncocommutative subalgebra ℋ_CM of noncommutative geometry. These three Hopf algebras admit a bigrading by n, the number of nodes, and an index k that specifies the degree of primitivity. In each case, we use iterations of the relevant coproduct to compute the dimensions of subspaces with modest values of n and k and infer a simple generating procedure for the remainder. The results for ℋ_ladder are familiar from the theory of partitions, while those for ℋ_CM involve novel transforms of partitions. Most beautiful is the bigrading of ℋ _R , the largest of the three. Thanks to Sloane's superseeker, we discovered that it saturates all possible inequalities. We prove this by using the universal Hochschild-closed one-cocycle B ₊, which plugs one set of divergences into another, and by generalizing the concept of natural growth beyond that entailed by the Connes–Moscovici case. We emphasize the yet greater challenge of handling the infinite set of decorations of realistic quantum field theory. Received: 31 January 2000 / Accepted: 7 July 2000     

A Dynamical Uncertainty Principle in von Neumann Algebras by Operator Monotone Functions

Suppose that A ₁,…,A _N are observables (selfadjoint matrices) and ρ is a state (density matrix). In this case the standard uncertainty principle, proved by Robertson, gives a bound for the quantum generalized variance, namely for det {Cov  _ρ (A _j ,A _k )}, using the commutators [A _j ,A _k ]; this bound is trivial when N is odd. Recently a different inequality of Robertson-type has been proved by the authors with the help of the theory of operator monotone functions. In this case the bound makes use of the commutators [ρ,A _j ] and is non-trivial for any N. In the present paper we generalize this new result to the von Neumann algebra case. Nevertheless the proof appears to simplify all the existing ones.     

Realizability of Point Processes

There are various situations in which it is natural to ask whether a given collection of k functions, ρ _j (r ₁,…,r _j ), j=1,…,k, defined on a set X, are the first k correlation functions of a point process on X. Here we describe some necessary and sufficient conditions on the ρ _j ’s for this to be true. Our primary examples are X=ℝ ^d , X=ℤ ^d , and X an arbitrary finite set. In particular, we extend a result by Ambartzumian and Sukiasian showing realizability at sufficiently small densities ρ ₁(r). Typically if any realizing process exists there will be many (even an uncountable number); in this case we prove, when X is a finite set, the existence of a realizing Gibbs measure with k body potentials which maximizes the entropy among all realizing measures. We also investigate in detail a simple example in which a uniform density ρ and translation invariant ρ ₂ are specified on ℤ; there is a gap between our best upper bound on possible values of ρ and the largest ρ for which realizability can be established.     

Sum-of-Squares Results for Polynomials Related to the Bessis–Moussa–Villani Conjecture

We show that the polynomial S _m,k (A,B), that is the sum of all words in noncommuting variables A and B having length m and exactly k letters equal to B, is not equal to a sum of commutators and Hermitian squares in the algebra R〈X,Y〉, where X ²=A and Y ²=B, for all even values of m and k with 6≤k≤m−10, and also for (m,k)=(12,6). This leaves only the case (m,k)=(16,8) open. This topic is of interest in connection with the Lieb–Seiringer formulation of the Bessis–Moussa–Villani conjecture, which asks whether Tr (S _m,k (A,B))≥0 holds for all positive semidefinite matrices A and B. These results eliminate the possibility of using “descent + sum-of-squares” to prove the BMV conjecture.     

Two energy scales and breakdown of renormalizability in the Kondo problem

It is shown that a second energy scale T ₀ ≫ T _k arises in the Kondo problem. Perturbation theory is valid only in the region T > T ₀. For this reason, the transition from weak to strong coupling occurs at temperatures much higher than the Kondo temperature T _k. Pis’ma Zh. éksp. Teor. Fiz. 70, No. 2, 106–111 (25 July 1999)     

Modeling 1/f noise using moving averages

A model of 1/f noise is considered, based on moving averages of ordern. The coefficientsα _k defining the model are calculated numerically using Seidel iteration which turns out to converge rapidly. The convergence is independent ofn which seems to be caused by the fact that the nonlinear problem solved is defined by a self-similar matrix. The coefficientsα _k appear to approach, with indefinitely growingn, valuesα _k =1/√k and thus the model has a kind of Fourier invariance. Physical interpretation of the invariance is suggested as well as of coefficientsα _k describing long-range correlations. Fractal sets of dimensiond=2.5 are proposed to play certain role in explaining the latter. This work was partly supported by the Scientific Grant Agency of the Ministry of Education of Slovak Republic (VEGA) under the Grant No. 1/3143/96 and by the Slovak Literary Fund.     

D-dimensional massless particle with extended gauge invariance

We propose the model ofD-dimensional massless particle whose Lagrangian is given by theN-th extrinsic curvature of world-line. The system hasN+1 gauge degrees of freedom constitutingW-like algebra; the classical trajectories of the model are space-like curves which obey the conditionsk _N+a=k_N−a, k_2N =0,a=1, ...,N−1,N≤[(D−2)/2], while the firstN curvaturesk _i remain arbitrary. We show that the model admits consistent formulation on the anti-DeSitter space. The solutions of the system are the massless irreducible representations of Poincaré group withN nonzero helicities, which are equal to each other. Presented at the 9th Colloquium “Quantum Groups and Integrable Systems”, Prague, 22–24 June 2000.     

Vortex Condensates¶for the SU(3) Chern–Simons Theory

We investigate SU(3)-periodic vortices in the self-dual Chern–Simons theory proposed by Dunne in [13, 15]. At the first admissible non-zero energy level E= 2 π, and for each (broken and unbroken) vacuum state φ₍₀₎ of the system, we find a family of periodic vortices asymptotically gauge equivalent to φ₍₀₎, as the Chern–Simons coupling parameter k→ 0. At higher energy levels, we show the existence of multiple gauge distinct periodic vortices with at least one of them asymptotically gauge equivalent to the (broken) principal embedding vacuum, when k→ 0. Received: 23 October 1999 / Accepted: 14 March 2000     

Bottomonium γ(5<Emphasis Type="Italic">S</Emphasis>) decays into <Emphasis Type="Italic">BB</Emphasis> and <Emphasis Type="Italic">BB</Emphasis>π

Two-and three-body decays of γ(5S) into BB, BB*, B*B*, B _s B _s , B _s B _s ^*, and BB*π, B*B*π are evaluated using the theory developed earlier for dipion-bottomonium transitions. The theory contains only two parameters—vertex masses M _br and M _ω—known from the dipion spectra and width. Predicted values of Γ_tot(5S) and six partial widths Γ _k (5S), k = BB, BB*, ... are in agreement with the experiment. The decay widths Γ_5S (πBB*) and Γ_5S (πB*B*) are also calculated and found to be on the order of 10 keV. The text was submitted by the authors in English.     

Trajectory attractor of a reaction-diffusion system with a series of zero diffusion coefficients

We study a reaction-diffusion system of N equations with k nonzero and N − k zero diffusion coefficients. More exactly, the first k equations of the system contain the terms a _i Δu _i − f _j (u, v), i = 1, …, k, with the diffusion coefficient a _i > 0. The right-hand sides of the other N − k equations contain only nonlinear interaction functions −h _j (u, v), j = k + 1, …, N, with zero diffusion. Here u = (u ₁, …, u _k ) and v = (υ _k+1, …, υ _N ) are unknown concentration vectors. Under appropriate assumptions on the interaction functions f(·) and h(·), we construct the trajectory attractor of this reaction-diffusion system. We also find the trajectory attractors , δ = (δ ₁, …, δ _k ), for the analogous reaction-diffusion systems having the terms δ _j Δυ _j − h _j (u, v), j = k + 1, …, N, with small diffusion coefficients δ _j ⩾ 0 in the last N − k equations. We prove that the trajectory attractors converge to (in an appropriate topology) as δ → 0+. Dedicated to the memory of Vladimir Borovikov Partially supported by the Russian Foundation for Basic Research (projects nos. 08-01-00784 and 07-01-00500).     

Neighborhood preferences in random matching problems

We consider a class of random matching problems where the distance between two points has a probability law which, for a small distance l, goes like l^r. In the framework of the cavity method, in the limit of an infinite number of points, we derive equations for p_k, the probability for some given point to be matched to its kth nearest neighbor in the optimal configuration. These equations are solved in two limiting cases: r = 0 -- where we recover p _k = 1/2^k, as numerically conjectured by Houdayer et al. and recently rigorously proved by Aldous -- and r→ + ∞. For 0 < r < + ∞, we are not able to solve the equations analytically, but we compute the leading behavior of p_k for large k. Received 14 February 2001     

High-energy strong interactions: from ‘hard’ to ‘soft’

We discuss the qualitative features of the recent data on multiparticle production observed at the LHC. The tolerable agreement with Monte Carlos based on LO DGLAP evolution indicates that there is no qualitative difference between ‘hard’ and ‘soft’ interactions; and that a perturbative QCD approach may be extended into the soft domain. However, in order to describe the data, these Monte Carlos need an additional infrared cutoff k _min with a value k _min∼2–3 GeV which is not small, and which increases with collider energy. Here we explain the physical origin of the large k _min. Using an alternative model which matches the ‘soft’ high-energy hadron interactions smoothly on to perturbative QCD at small x, we demonstrate that this effective cutoff k _min is actually due to the strong absorption of low k _t partons. The model embodies the main features of the BFKL approach, including the diffusion in transverse momenta, ln k _t , and an intercept consistent with resummed next-to-leading log corrections. Moreover, the model uses a two-channel eikonal framework, and includes the contributions from the multi-Pomeron exchange diagrams, both non-enhanced and enhanced. The values of a small number of physically-motivated parameters are chosen to reproduce the available total, elastic and proton dissociation cross section (pre-LHC) data. Predictions are made for the LHC, and the relevance to ultra-high-energy cosmic rays is briefly discussed. The low x inclusive integrated gluon PDF, and the diffractive gluon PDF, are calculated in this framework, using the parameters which describe the high-energy pp and p[`(p)]p\bar{p} ‘soft’ data. Comparison with the PDFs obtained from the global parton analyses of deep inelastic and related hard scattering data and from diffractive deep inelastic data looks encouraging.     

On Noncommutative Multi-Solitons

 We find the moduli space of multi-solitons in noncommutative scalar field theories at large θ, in arbitrary dimension. The existence of a non-trivial moduli space at leading order in 1/θ is a consequence of a Bogomolnyi bound obeyed by the kinetic energy of the θ=∞ solitons. In two spatial dimensions, the parameter space for k solitons is a K?hler de-singularization of the symmetric product (ℝ²) ^k /S _k . We exploit the existence of this moduli space to construct solitons on quotient spaces of the plane: ℝ²/ℤ _k , cylinder, and T ² . However, we show that tori of area less than or equal to 2πθ do not admit stable solitons. In four dimensions the moduli space provides an explicit K?hler resolution of (ℝ⁴) ^k /S _k . In general spatial dimension 2d, we show it is isomorphic to the Hilbert scheme of k points in ℂ ^d , which for d>2 (and k>3) is not smooth and can have multiple branches. Received: 29 May 2001 / Accepted: 16 August 2002 Published online: 7 November 2002 Communicated by R.H. Dijkgraaf     

Estimation of Kauzmann temperature and isenthalpic temperature of metallic glasses

We propose expressions for the estimation of the isenthalpic temperature T ₀ (T ₀ = αT _m , α is a semi-empirical parameter and 0 ⩽ α < 1, T _m is the solidus temperature) and the Kauzmann temperature T _k (T _k = T _m exp(α−1)) for glass forming alloys. It is found that T _k estimated by T _k = T _m exp(α−1) is in agreement with that directly calculated from the heat capacity data, indicating that T _k = T _m exp(α − 1) can be used to estimate T _k of glass forming alloys. T ₀ estimated by T ₀ = αT _m , on the other hand, widely deviates from that of directly calculated from the heat capacity data. This suggests that the enthalpy difference of the under-cooled liquid and the crystal might be a nonlinear function of the temperature below T _k . Moreover, the Gibbs free energy difference ΔG is not sensitive to the deviation of α.     

Spectrum of Kelvin-wave turbulence in superfluids

We derive a type of kinetic equation for Kelvin waves on quantized vortex filaments with random large-scale curvature, that describes step-by-step (local) energy cascade over scales caused by 4-wave interactions. Resulting new energy spectrum E _LN(k) ∝ k ^−5/3 must replace in future theory (e.g., in finding the quantum turbulence decay rate) the previously used spectrum E _KS(k) ∝ k ^−7/5, which was recently shown to be inconsistent due to nonlocality of the 6-wave energy cascade.     

Static charged spheres with anisotropic pressure in general relativity

We report a generalization of our earlier formalism [Pramana, 54, 663 (1998)] to obtain exact solutions of Einstein-Maxwell’s equations for static spheres filled with a charged fluid having anisotropic pressure and of null conductivity. Defining new variables: w=(4π/3)(ρ+ε)r ², u=4πξr ², v _r=4πp _r r ², v _⊥=4πp _⊥ r ²[ρ, ξ(=−(1/2)F ₁₄ F ¹⁴), p _r, p _⊥ being respectively the energy densities of matter and electrostatic fields, radial and transverse fluid pressures whereas ε denotes the eigenvalue of the conformal Weyl tensor and interpreted as the energy density of the free gravitational field], we have recast Einstein’s field equations into a form easy to integrate. Since the system is underdetermined we make the following assumptions to solve the field equations (i) u=v _r=(a ²/2κ)r ⁿ⁺², v _⊥=k ₁ v _r, w=k ₂ v _r; a ², n(>0), k ₁, k ₂ being constants with κ=((k ₁+2)/3+k ₂) and (ii) w+u=(b ²/2)r ⁿ⁺², u=v _r, v _⊥−v _r=k, with b and k as constants. In both cases the field equations are integrated completely. The first solution is regular in the metric as well as physical variables for all values of n>0. Even though the second solution contains terms like k/r ² since Q(0)=0 it is argued that the pressure anisotropy, caused by the electric flux near the centre, can be made to vanish reducing it to the generalized Cooperstock-de la Cruz solution given in [14]. The interior solutions are shown to match with the exterior Reissner-Nordstrom solution over a fixed boundary. Dedicated to Prof. F A E Pirani.