A Matter of Disorder: Monte Carlo Simulations of Phase Transitions in Strongly Disordered Systems

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, 20070 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 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

Superfluid Phase Transitions in Disordered Systems

Hannes Meier

This thesis presents results from large scale Monte Carlo simulations of systems subject to a superfluid phase transition in the presence of disorder. The simulations are performed by state-of-the-art, collective Monte Carlo algorithms treating phase degrees of freedom in effective models with amplitude fluctuations integrated out. In Paper I a model system for the possible solid to supersolid transition in 4He is presented.The Wolff cluster algorithm is used to study how the presence of linearly correlated random defects is able to alter the universality class of the 3-dimensional XY-model. In the pure case the superfluid density and heat capacity have singular onsets, which are not seen in the supersolid experiments where instead a smooth onset is obtained. Using finite size scaling of Monte Carlo data, we find a similar smooth onset in our simulations, governed by exponents ν=1 for the superfluid density and α=-1 for the heat capacity. These results are in qualitative agreement with experiments for the observed transition in solid 4He. In Paper II a systematic investigation of the scaling result z=d for the dynamic critical exponentat the Bose glass to superfluid quantum phase transition is performed. The result z=d has been believed to be exact for about 20 years, but although it has been questioned lately no accurate estimate of z has been available. An effective link current model of quantum bosons at T=0 with disorder in 2D is simulated using highly effective worm Monte Carlo simulations.The data analysis is based on a finite size scaling approach todetermine the quantum correlation time from simulationdata for boson world lines without any a priori assumption on the critical parameters. The resulting critical exponents are z=1.8 \pm 0.05, ν=1.15 \pm 0.03, and η=-0.3 \pm 0.1. This suggests that z=d is not satisfied.This thesis presents results from large scale Monte Carlo simulations of systems subject to a superfluid phase transition in the presence of disorder. The simulations are performed by state-of-the-art, collective Monte Carlo algorithms treating phase degrees of freedom in effective models with amplitude fluctuations integrated out. In Paper I a model system for the possible solid to supersolid transition in 4He is presented.The Wolff cluster algorithm is used to study how the presence of linearly correlated random defects is able to alter the universality class of the 3-dimensional XY-model. In the pure case the superfluid density and heat capacity have singular onsets, which are not seen in the supersolid experiments where instead a smooth onset is obtained. Using finite size scaling of Monte Carlo data, we find a similar smooth onset in our simulations, governed by exponents ν=1 for the superfluid density and α=-1 for the heat capacity. These results are in qualitative agreement with experiments for the observed transition in solid 4He. In Paper II a systematic investigation of the scaling result z=d for the dynamic critical exponentat the Bose glass to superfluid quantum phase transition is performed. The result z=d has been believed to be exact for about 20 years, but although it has been questioned lately no accurate estimate of z has been available. An effective link current model of quantum bosons at T=0 with disorder in 2D is simulated using highly effective worm Monte Carlo simulations.The data analysis is based on a finite size scaling approach todetermine the quantum correlation time from simulationdata for boson world lines without any a priori assumption on the critical parameters. The resulting critical exponents are z=1.8 \pm 0.05, ν=1.15 \pm 0.03, and η=-0.3 \pm 0.1. This suggests that z=d is not satisfied.

Jan 1, 20110 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}.