V. Laliena

Hybrid Monte Carlo algorithm for the double exchange model

Abstract The Hybrid Monte Carlo algorithm is adapted to the simulation of a system of classical degrees of freedom coupled to non self-interacting lattices fermions. The diagonalization of the Hamiltonian matrix is avoided by introducing a path-integral formulation of the problem, in d+1 Euclidean space–time. A perfect action formulation allows to work on the continuum Euclidean time, without need for a Trotter–Suzuki extrapolation. To demonstrate the feasibility of the method we study the Double Exchange Model in three dimensions. The complexity of the algorithm grows only as the system volume, allowing to simulate in lattices as large as 16 3 on a personal computer. We conclude that the second order paramagnetic–ferromagnetic phase transition of Double Exchange Materials close to half-filling belongs to the Universality Class of the three-dimensional classical Heisenberg model.

New proposal for numerical simulations of theta-vacuum-like systems.

We propose a new approach to perform numerical simulations of theta-vacuum-like systems, test it in two analytically solvable models, and apply it to CP3. The main new ingredient in our approach is the method used to compute the probability distribution function of the topological charge at theta=0. We do not get unphysical phase transitions (flattening behavior of the free energy density) and reproduce the exact analytical results for the order parameter in the whole theta range within a few percent.

Interplay between double-exchange, superexchange, and Lifshitz localization in doped manganites

Considering the disorder caused in manganites by the substitution Mn→Fe or Ga, we accomplish a systematic study of doped manganites begun in previous papers. To this end, a disordered model is formulated and solved using the variational mean-field technique. The subtle interplay between double exchange, superexchange, and disorder causes similar effects on the dependence of T_(C) on the percentage of Mn substitution in the cases considered. Yet, in La_(2/3)Ca_(1/3)Mn_(1-y)Ga_(y)O_(3) our results suggest a quantum critical point (QCP) for y ≈ 0.1–0.2, associated to the localization of the electronic states of the conduction band. In the case of La_(x)Ca_(x)Mn_(1-y)Fe_(y)O_(3) (with x = 1/3,3/8) no such QCP is expected.

Variational mean-field approach to the double-exchange model

It has been recently shown that the double exchange Hamiltonian, with weak antiferromagnetic interactions, has a richer variety of first- and second-order transitions than previously anticipated, and that such transitions are consistent with the magnetic properties of manganites. Here we present a thorough discussion of the variational mean-field approach that leads to these results. We also show that the effect of the Berry phase turns out to be crucial to produce first-order paramagnetic-ferromagnetic transitions near half filling with transition temperatures compatible with the experimental situation. The computation relies on two crucial facts: the use of a mean-field ansatz that retains the complexity of a system of electrons with off-diagonal disorder, not fully taken into account by the mean-field techniques, and the small but significant antiferromagnetic superexchange interaction between the localized spins.

Finite density QCD: a new approach

We introduce a new approach to analyze the phase diagram of QCD at finite chemical potential and temperature, test it in the Gross-Neveu model at finite baryon density, and apply it to the study of the chemical potential-temperature phase diagram of QCD with four degenerate flavors of Kogut-Susskind type.

Phase diagram of QCD with four quark flavors at finite temperature and baryon density

Abstract We analyze the phase diagram of QCD with four staggered flavors in the ( μ , T ) plane using a method recently proposed by us. We explore the region T ≳ 0.7 T C and μ ≲ 1.4 T C , where T C is the transition temperature at zero baryon density, and find a first order transition line. Our results are quantitatively compatible with those obtained with the imaginary chemical potential approach and the double reweighting method, in the region where these approaches are reliable, T ≳ 0.9 T C and μ ≲ T C . But, in addition, our method allows us to extend the transition line to lower temperatures and higher chemical potentials.

Spontaneous symmetry breaking in fermion-gauge systems: a non standard approach

Abstract We propose a new method for the study of the chiral properties of the ground state in QFTs based on the computation of the probability distribution function of the chiral condensate. It can be applied directly in the chiral limit and therefore no mass extrapolations are needed. Furthermore this approach allows to write up equations relating the chiral condensate with quantities computable by standard numerical methods, the functional form of these relations depending on the broken symmetry group. As a check, we report some results for the compact Schwinger model.

Diquark condensation at strong coupling

The possibility of diquark condensation at sufficiently large baryon chemical potential and zero temperature is analyzed in QCD at strong coupling. In agreement with other strong coupling analysis, it is found that a first order phase transition separates a low density phase with chiral symmetry spontaneously broken from a high density phase where chiral symmetry is restored. In none of the phases diquark condensation takes place as an equilibrium state, but, for any value of the chemical potential, there is a metastable state characterized by a non-vanishing diquark condensate. The energy difference between this metastable state and the equilibrium state decreases with the chemical potential and is minimum in the high density phase. The results indicate that there is attraction in the quark-quark sector also at strong coupling, and that the attraction is more effective at high baryon density, but for infinite coupling it is not enough to produce diquark condensation. It is argued that the absence of diquark condensation is not a peculiarity of the strong coupling limit, but persists at sufficiently large finite couplings.

θ-vacuum systems via real action simulations

Abstract Inspired by the results of the Ising model within an imaginary external magnetic field, we introduce a transformation in quantum systems with a θ-vacuum term that amounts to a rescaling of z= cos θ 2 . Making use of this transformation we are able to determine the order parameter as a function of θ. The approach is successfully tested in models with both broken and unbroken CP symmetry at θ=π.

Chiral susceptibilities in noncompact QED: a new determination of the γ exponent and the critical couplings

Abstract We report the results of a measurement of susceptibilities in noncompact QED 4 in 8 4 , 10 4 and 12 4 lattices. Due to the potentialities of the MFA approach, we have done simulations in the chiral limit which are therefore free from arbitrary mass extrapolations. Our results in the Coulomb phase show unambiguously that the susceptibility critical exponent γ = 1 independently of the flavour symmetry group. The critical couplings extracted from these calculations are in perfect agreement with previous determinations based on the fermion effective action and plaquette energy, and outside the predictions of a logarithmically improved scalar mean field theory by eight standard deviations.