Monte Carlo Simulations of a Lattice Gas with Logarithmic Interactions

Phase transitions and their critical scaling properties, especially in systems with disorder, are important both for our theoretical understanding of our environment, but also for their practical use in applications and materials in our everyday life. This thesis presents results from finite size scaling analysis of critical phenomena in systems with disorder, using high-precision Monte Carlo simulations and state of the art numerical methods. Specifically, theoretical models suitable for simulations in the presence of uncorrelated or correlated disorder are studied. Uncorrelated strong disorder, as present in the two dimensional gauge glass model to study the vortex glass phase of high temperature superconductors in an applied magnetic field is shown to lack a finite temperature phase transition. Further, results from dynamic quantities, such as resistance and autocorrelation functions, indicate the existence of two distinct diverging correlation times, one associated with local relaxation and one associated with vortex phase slips. Correlated disorder is studied both in the superfluid transition of helium-4 and in the anisotropic critical scaling of a transverse Meissner-like transition in an experimental setup of a high temperature superconductor. For the superfluid helium transition, it is shown that the presence of fractally correlated disorder presumably alters the universality class of the pure model. Also, a comparison with experimental data suggests that the critical scaling theory describing the heat capacity of helium-4 may need to be modified in the presence of the disorder. In the case of superconductors, analyzing experimental data from resistance measurements in a system with columnar defects together with an anisotropy in the applied magnetic field, reveals a fully anisotropic scaling regime. Finally, a data analysis is presented from simulations of a charged particle gas system in three dimensions, where the normal Coulomb interaction between charges is changed into a logarithmic interaction. Previous work indicates the possibility of a transition similar to the Kosterlitz-Thouless transition in certain two dimensional systems. On the contrary, our simulations seem to favor a system whose critical scaling behavior is consistent with a transition occurring only at zero critical temperature. Overall, disorder in the model systems studied leads to important modifications of the critical scaling properties of pure systems, and thereby also to possible changes of the corresponding universality classes. This results in interesting predictions with experimentally relevant consequences.Phase transitions and their critical scaling properties, especially in systems with disorder, are important both for our theoretical understanding of our environment, but also for their practical use in applications and materials in our everyday life. This thesis presents results from finite size scaling analysis of critical phenomena in systems with disorder, using high-precision Monte Carlo simulations and state of the art numerical methods. Specifically, theoretical models suitable for simulations in the presence of uncorrelated or correlated disorder are studied. Uncorrelated strong disorder, as present in the two dimensional gauge glass model to study the vortex glass phase of high temperature superconductors in an applied magnetic field is shown to lack a finite temperature phase transition. Further, results from dynamic quantities, such as resistance and autocorrelation functions, indicate the existence of two distinct diverging correlation times, one associated with local relaxation and one associated with vortex phase slips. Correlated disorder is studied both in the superfluid transition of helium-4 and in the anisotropic critical scaling of a transverse Meissner-like transition in an experimental setup of a high temperature superconductor. For the superfluid helium transition, it is shown that the presence of fractally correlated disorder presumably alters the universality class of the pure model. Also, a comparison with experimental data suggests that the critical scaling theory describing the heat capacity of helium-4 may need to be modified in the presence of the disorder. In the case of superconductors, analyzing experimental data from resistance measurements in a system with columnar defects together with an anisotropy in the applied magnetic field, reveals a fully anisotropic scaling regime. Finally, a data analysis is presented from simulations of a charged particle gas system in three dimensions, where the normal Coulomb interaction between charges is changed into a logarithmic interaction. Previous work indicates the possibility of a transition similar to the Kosterlitz-Thouless transition in certain two dimensional systems. On the contrary, our simulations seem to favor a system whose critical scaling behavior is consistent with a transition occurring only at zero critical temperature. Overall, disorder in the model systems studied leads to important modifications of the critical scaling properties of pure systems, and thereby also to possible changes of the corresponding universality classes. This results in interesting predictions with experimentally relevant consequences.

Jan 1, 20073 citations

From Order to Disorder in High Temperature Superconductors

Anders Vestergren

Phase transitions in a number of models related to hightemperature superconductors are investigated, using scalingmethods and Monte Carlo simulations. This thesis considers twomain topics. The first topic is phase transitions, phase diagrams, andvortex motion in high temperature superconductors at finitetemperature, subject to magnetic fields and disorder. We studya vortex glass model at finite temperature, with stronguncorrelated vortex pinning and a magnetic field. We find thatthe vortex glass exists at finite temperature and calculate thecritical exponents of the transition. We also investigate hightemperature superconductors with columnar disorder in zero andapplied magnetic fields. Some of these studies are alsorelevant for the superfluid to Mott insulator transition ofbosons in two dimensions. We find that the unscreened Boseglass transition belongs to a new universality class. Wecalculate the critical exponents of the superconductingtransition with columnar defects in zero applied magneticfield. The transverse Meissner transition is studied, and wefind an exotic universality class with a correlation volumethat is infinitely anisotropic in all directions. The second topic is confinement-deconfinement transitions incompact Abelian Higgs models. We develop a new order parameter,related to a large Wilson loop for fractionalized charges, anduse it to study the concept of topological order. Thesetransitions may be relevant for strongly correlated electronsin two dimensions.Phase transitions in a number of models related to hightemperature superconductors are investigated, using scalingmethods and Monte Carlo simulations. This thesis considers twomain topics. The first topic is phase transitions, phase diagrams, andvortex motion in high temperature superconductors at finitetemperature, subject to magnetic fields and disorder. We studya vortex glass model at finite temperature, with stronguncorrelated vortex pinning and a magnetic field. We find thatthe vortex glass exists at finite temperature and calculate thecritical exponents of the transition. We also investigate hightemperature superconductors with columnar disorder in zero andapplied magnetic fields. Some of these studies are alsorelevant for the superfluid to Mott insulator transition ofbosons in two dimensions. We find that the unscreened Boseglass transition belongs to a new universality class. Wecalculate the critical exponents of the superconductingtransition with columnar defects in zero applied magneticfield. The transverse Meissner transition is studied, and wefind an exotic universality class with a correlation volumethat is infinitely anisotropic in all directions. The second topic is confinement-deconfinement transitions incompact Abelian Higgs models. We develop a new order parameter,related to a large Wilson loop for fractionalized charges, anduse it to study the concept of topological order. Thesetransitions may be relevant for strongly correlated electronsin two dimensions.

Jan 1, 20040 citations

Phase transitions with critical loop disorder

Hannes Meier, Mats Wallin

Condensed matter systems undergoing phase transitions rarely allow exact solutions. The presence of disorder renders the situation even worse but collective Monte Carlo methods and parallel algorithms allow numerical descriptions. This thesis considers classical phase transitions in disordered spin systems in general and in effective models of superfluids with disorder and novel interactions in particular. Quantum phase transitions are considered via a quantum to classical mapping. Central questions are if the presence of defects changes universal properties and what qualitative implications follow for experiments. Common to the cases considered is that the disorder maps out correlated structures. All results are obtained using large-scale Monte Carlo simulations of effective models capturing the relevant degrees of freedom at the transition.Considering a model system for superflow aided by a defect network, we find that the onset properties are significantly altered compared to the $\lambda$ -transition in $^{4}$ He. This has qualitative implications on expected experimental signatures in a defect supersolid scenario.For the Bose glass to superfluid quantum phase transition in 2D we determine the quantum correlation time by an anisotropic finite size scaling approach. Without a priori assumptions on critical parameters, we find the critical exponent $z=1.8 \pm 0.05$ contradicting the long standing result $z=d$ .Using a 3D effective model for multi-band type-1.5 superconductors we find that these systems possibly feature a strong first order vortex-driven phase transition. Despite its short-range nature details of the interaction are shown to play an important role.Phase transitions in disordered spin models exposed to correlated defect structures obtained via rapid quenches of critical loop and spin models are investigated. On long length scales the correlations are shown to decay algebraically. The decay exponents are expressed through known critical exponents of the disorder generating models. For cases where the disorder correlations imply the existence of a new long-range-disorder fixed point we determine the critical exponents of the disordered systems via finite size scaling methods of Monte Carlo data and find good agreement with theoretical expectations.Condensed matter systems undergoing phase transitions rarely allow exact solutions. The presence of disorder renders the situation even worse but collective Monte Carlo methods and parallel algorithms allow numerical descriptions. This thesis considers classical phase transitions in disordered spin systems in general and in effective models of superfluids with disorder and novel interactions in particular. Quantum phase transitions are considered via a quantum to classical mapping. Central questions are if the presence of defects changes universal properties and what qualitative implications follow for experiments. Common to the cases considered is that the disorder maps out correlated structures. All results are obtained using large-scale Monte Carlo simulations of effective models capturing the relevant degrees of freedom at the transition.Considering a model system for superflow aided by a defect network, we find that the onset properties are significantly altered compared to the $\lambda$-transition in $^{4}$He. This has qualitative implications on expected experimental signatures in a defect supersolid scenario.For the Bose glass to superfluid quantum phase transition in 2D we determine the quantum correlation time by an anisotropic finite size scaling approach. Without a priori assumptions on critical parameters, we find the critical exponent $z=1.8 \pm 0.05$ contradicting the long standing result $z=d$.Using a 3D effective model for multi-band type-1.5 superconductors we find that these systems possibly feature a strong first order vortex-driven phase transition. Despite its short-range nature details of the interaction are shown to play an important role.Phase transitions in disordered spin models exposed to correlated defect structures obtained via rapid quenches of critical loop and spin models are investigated. On long length scales the correlations are shown to decay algebraically. The decay exponents are expressed through known critical exponents of the disorder generating models. For cases where the disorder correlations imply the existence of a new long-range-disorder fixed point we determine the critical exponents of the disordered systems via finite size scaling methods of Monte Carlo data and find good agreement with theoretical expectations.

Jan 1, 20151 citation

Phase transitions in effective lattice models for strongly correlated systems

Steinar Kragset

In three research articles we have studied the critical properties of effective lattice models for strongly correlated electron systems by Monte Carlo simulations. A similar model is used in a fourth article for investigating thermal fluctuations of vortices in a rotating Bose–Einstein condensate. In the first part of this thesis we review the necessary background and introduce the models one by one. The last part is a collection of the papers.\n\nPaper I [1]: We consider the scaling of the mean square dipole moment in a plasma with logarithmic interactions in a two- and three-dimensional system. In both cases, we establish the existence of a low-temperature regime where the mean square dipole moment does not scale with system size and a hightemperature regime does scale with system size. Thus, there is a nonanalytic change in the polarizability of the system as a function of temperature, and hence a metal-insulator transition in both cases. The relevance of this transition in three dimensions to quantum phase transitions in 2 + 1-dimensional systems is briefly discussed.\n\nPaper II [2]: The existence of a discontinuity in the inverse dielectric constant of the two-dimensional Coulomb gas is demonstrated on purely numerical grounds. This is done by expanding the free energy in an applied twist and performing a finite-size scaling analysis of the coefficients of higher-order terms. The phase transition, driven by unbinding of dipoles, corresponds to the Kosterlitz-Thouless transition in the 2D XY model. The method developed is also used for investigating the possibility of a Kosterlitz-Thouless phase transition in a threedimensional system of point charges interacting with a logarithmic pair-potential, a system related to effective theories of low-dimensional strongly correlated systems. We also contrast the finite-size scaling of the fluctuations of the dipole moments of the two-dimensional Coulomb gas and the three-dimensional logarithmic system to those of the three-dimensional Coulomb gas.\n\nPaper III [3]: We perform large-scale Monte Carlo simulations on an effective gauge theory for an easy plane quantum anti-ferromagnet, including a Berry phase term that projects out the S = 1/2 sector. Without a Berry phase term, the model exhibits a phase transition in the 3DXY universality class associated with proliferation of gauge-charge neutral U(1) vortices. The instantons that eliminate the phase transition in the gauge-charged sector are cancelled by the Berry phases. The result is a first order phase transition. This gauge theory therefore does not exhibit deconfined criticality.\n\nPaper IV [4]: We perform Monte Carlo studies of vortices in three dimensions in a cylindrical confinement, with uniform and nonuniform density. The former is relevant to rotating 4He, the latter is relevant to a rotating trapped Bose condensate. In the former case we find dominant angular thermal vortex fluctuations close to the cylinder wall. For the latter case, a novel effect is that at low temperatures the vortex solid close to the center of the trap crosses directly over to a tension-less vortex tangle near the edge of the trap. At higher temperatures an intermediate tensionful vortex liquid located between the vortex solid and the vortex tangle, may exist.In three research articles we have studied the critical properties of effective lattice models for strongly correlated electron systems by Monte Carlo simulations. A similar model is used in a fourth article for investigating thermal fluctuations of vortices in a rotating Bose–Einstein condensate. In the first part of this thesis we review the necessary background and introduce the models one by one. The last part is a collection of the papers.\n\nPaper I [1]: We consider the scaling of the mean square dipole moment in a plasma with logarithmic interactions in a two- and three-dimensional system. In both cases, we establish the existence of a low-temperature regime where the mean square dipole moment does not scale with system size and a hightemperature regime does scale with system size. Thus, there is a nonanalytic change in the polarizability of the system as a function of temperature, and hence a metal-insulator transition in both cases. The relevance of this transition in three dimensions to quantum phase transitions in 2 + 1-dimensional systems is briefly discussed.\n\nPaper II [2]: The existence of a discontinuity in the inverse dielectric constant of the two-dimensional Coulomb gas is demonstrated on purely numerical grounds. This is done by expanding the free energy in an applied twist and performing a finite-size scaling analysis of the coefficients of higher-order terms. The phase transition, driven by unbinding of dipoles, corresponds to the Kosterlitz-Thouless transition in the 2D XY model. The method developed is also used for investigating the possibility of a Kosterlitz-Thouless phase transition in a threedimensional system of point charges interacting with a logarithmic pair-potential, a system related to effective theories of low-dimensional strongly correlated systems. We also contrast the finite-size scaling of the fluctuations of the dipole moments of the two-dimensional Coulomb gas and the three-dimensional logarithmic system to those of the three-dimensional Coulomb gas.\n\nPaper III [3]: We perform large-scale Monte Carlo simulations on an effective gauge theory for an easy plane quantum anti-ferromagnet, including a Berry phase term that projects out the S = 1/2 sector. Without a Berry phase term, the model exhibits a phase transition in the 3DXY universality class associated with proliferation of gauge-charge neutral U(1) vortices. The instantons that eliminate the phase transition in the gauge-charged sector are cancelled by the Berry phases. The result is a first order phase transition. This gauge theory therefore does not exhibit deconfined criticality.\n\nPaper IV [4]: We perform Monte Carlo studies of vortices in three dimensions in a cylindrical confinement, with uniform and nonuniform density. The former is relevant to rotating 4He, the latter is relevant to a rotating trapped Bose condensate. In the former case we find dominant angular thermal vortex fluctuations close to the cylinder wall. For the latter case, a novel effect is that at low temperatures the vortex solid close to the center of the trap crosses directly over to a tension-less vortex tangle near the edge of the trap. At higher temperatures an intermediate tensionful vortex liquid located between the vortex solid and the vortex tangle, may exist.

Jan 1, 20060 citations

Phase transitions in systems with critical cluster defects

Hannes Meier, Mats Wallin

Jan 1, 20151 citation

Unscreened universality class for superconductors with columnar disorder

Anders Vestergren, Jack Lidmar, Mats Wallin

The phase transition in a model for vortex lines in high temperature superconductors with columnar defects, i.e., linearly correlated quenched random disorder, is studied with finite size scaling and Monte Carlo simulations. Previous studies of critical properties have mainly focused on the limit of strongly screened vortex line interactions. Here the opposite limit of weak screening is considered. The simulation results provide evidence for a distinct universality class, with values of the critical exponents that differ from the case of strong screening of the vortex interaction. In particular, scaling is anisotropic and characterized by a nontrivial value of the anisotropy exponent $\ensuremath{\zeta}={\ensuremath{\nu}}_{\ensuremath{\Vert}}/{\ensuremath{\nu}}_{\ensuremath{\perp}}.$ The exponents we find, $\ensuremath{\zeta}=1.25\ifmmode\pm\else\textpm\fi{}0.1,{\ensuremath{\nu}}_{\ensuremath{\perp}}=1.0\ifmmode\pm\else\textpm\fi{}0.1,z=1.95\ifmmode\pm\else\textpm\fi{}0.1,$ are similar to certain experimental results for ${\mathrm{YBa}}_{2}{\mathrm{Cu}}_{3}{\mathrm{O}}_{7}.$ The phase transition in a model for vortex lines in high temperature superconductors with columnar defects, i.e., linearly correlated quenched random disorder, is studied with finite size scaling and Monte Carlo simulations. Previous studies of critical properties have mainly focused on the limit of strongly screened vortex line interactions. Here the opposite limit of weak screening is considered. The simulation results provide evidence for a distinct universality class, with values of the critical exponents that differ from the case of strong screening of the vortex interaction. In particular, scaling is anisotropic and characterized by a nontrivial value of the anisotropy exponent $\ensuremath{\zeta}={\ensuremath{\nu}}_{\ensuremath{\Vert}}/{\ensuremath{\nu}}_{\ensuremath{\perp}}.$ The exponents we find, $\ensuremath{\zeta}=1.25\ifmmode\pm\else\textpm\fi{}0.1,{\ensuremath{\nu}}_{\ensuremath{\perp}}=1.0\ifmmode\pm\else\textpm\fi{}0.1,z=1.95\ifmmode\pm\else\textpm\fi{}0.1,$ are similar to certain experimental results for ${\mathrm{YBa}}_{2}{\mathrm{Cu}}_{3}{\mathrm{O}}_{7}.$

arXiv:cond-mat/0208225Physical review. B, Condensed matterMar 4, 20038 citations

\lambda

-transition in

^{4}

He. This has qualitative implications on expected experimental signatures in a defect supersolid scenario.For the Bose glass to superfluid quantum phase transition in 2D we determine the quantum correlation time by an anisotropic finite size scaling approach. Without a priori assumptions on critical parameters, we find the critical exponent

z=1.8 \pm 0.05

contradicting the long standing result

z=d

.Using a 3D effective model for multi-band type-1.5 superconductors we find that these systems possibly feature a strong first order vortex-driven phase transition. Despite its short-range nature details of the interaction are shown to play an important role.Phase transitions in disordered spin models exposed to correlated defect structures obtained via rapid quenches of critical loop and spin models are investigated. On long length scales the correlations are shown to decay algebraically. The decay exponents are expressed through known critical exponents of the disorder generating models. For cases where the disorder correlations imply the existence of a new long-range-disorder fixed point we determine the critical exponents of the disordered systems via finite size scaling methods of Monte Carlo data and find good agreement with theoretical expectations.

\lambda

-transition in

^{4}

z=1.8 \pm 0.05

contradicting the long standing result

z=d

{\mathrm{YBa}}_{2}{\mathrm{Cu}}_{3}{\mathrm{O}}_{7}.