Disorder-induced phase transitions in two-dimensional crystals

Effect of disorder on phases across two-dimensional thermal melting

We study melting in a two-dimensional system of classical particles with Gaussian-core interactions in disordered environments. The pure system validates the conventional two-step melting with a hexatic phase intervening between the solid and the liquid. This picture is modified in the presence of pinning impurities. A random distribution of pinning centers forces a hexaticlike low-temperature phase that transits into a liquid at a single melting temperature T_{m}^{RP}. In contrast, pinning centers located at randomly chosen sites of a perfect crystal anchor a solid at low temperatures which undergoes a direct transition to the liquid at T_{m}^{CP}. Thus, the two-step melting is lost in either case of disorder. We discuss the characteristics of melting depending on the nature of the impurities.

Structural transition in two-dimensional Hertzian spheres in the presence of random pinning

Using molecular dynamics simulation we have investigated the influence of random pinning on the phase diagram and melting scenarios of a two-dimensional system with the Hertz potential for α=5/2. It has been shown that random pinning can cardinally change the mechanism of first-order transition between the different crystalline phases (triangular and square) by virtue of generating hexatic and tetratic phases: a triangular crystal to hexatic transition is of the continuous Berezinskii-Kosterlitz-Thouless (BKT) type, a hexatic to tetratic transition is of first order, and finally, there is a continuous BKT-type transition from tetratic to the square crystal.

Dense hard disk ordering: influence of bidispersity and quenched disorder

Using Monte Carlo simulations, the impact on structural ordering in two-dimensional systems via the interplay of size bidispersity and quenched disorder in the form of an externally applied spatially random potential, is studied for a system of hard disks. By scanning across a wide range of dense packing fractions, size ratios and roughness of the applied potential, the phase diagram is constructed, which demonstrates that both quenched and size disorders shift the onset of translational order to higher packings, while maintaining the presence of the intermediate hexatic phase. At larger disorder strengths, the signatures of structural order are absent within the range of investigated packing fractions. Further, the dynamics with increasing potential strength is analysed for the mono-component system to obtain a spatio-temporal description of the melting process. Finally, the influence of the externally rough field on the Mermin-Wagner fluctuations, characteristic to two-dimensional systems, is investigated.

Melting of the Vortex Lattice through Intermediate Hexatic Fluid in an a-MoGe Thin Film

The hexatic fluid refers to a phase in between a solid and a liquid that has short-range positional order but quasi-long-range orientational order. In the celebrated theory of Berezinskii, Kosterlitz, and Thouless and subsequently refined by Halperin, Nelson, and Young, it was predicted that a two-dimensional hexagonal solid can melt in two steps: first, through a transformation from a solid to a hexatic fluid, which retains quasi-long-range orientational order; and then from a hexatic fluid to an isotropic liquid. In this Letter, using a combination of real space imaging and transport measurements, we show that the two-dimensional vortex lattice in an a-MoGe thin film follows this sequence of melting as the magnetic field is increased. Identifying the signatures of various transitions on the bulk transport properties of the superconductor, we construct a vortex phase diagram for a two-dimensional superconductor.

Disorder-Induced Topological State Transition in Photonic Metamaterials

The topological state transition has been widely studied based on the quantized topological band invariant such as the Chern number for the system without intense randomness that may break the band structures. We numerically demonstrate the disorder-induced state transition in the photonic topological systems for the first time. Instead of applying the ill-defined topological band invariant in a disordered system, we utilize an empirical parameter to unambiguously illustrate the state transition of the topological metamaterials. Before the state transition, we observe a robust surface state with well-confined electromagnetic waves propagating unidirectionally, immune to the disorder from permittivity fluctuation up to 60% of the original value. During the transition, a hybrid state composed of a quasiunidirectional surface mode and intensively localized hot spots is established, a result of the competition between the topological protection and Anderson localization.

Dewetting and spreading transitions for active matter on random pinning substrates

We show that sterically interacting self-propelled disks in the presence of random pinning substrates exhibit transitions among a variety of different states. In particular, from a phase separated cluster state, the disks can spread out and homogeneously cover the substrate in what can be viewed as an example of an active matter wetting transition. We map the location of this transition as a function of activity, disk density, and substrate strength, and we also identify other phases including a cluster state, coexistence between a cluster and a labyrinth wetted phase, and a pinned liquid. Convenient measures of these phases include the cluster size, which dips at the wetting-dewetting transition, and the fraction of sixfold coordinated particles, which drops when dewetting occurs.

Two-Dimensional Disorder in Black Phosphorus and Monochalcogenide Monolayers

Ridged, orthorhombic two-dimensional atomic crystals with a bulk Pnma structure such as black phosphorus and monochalcogenide monolayers are an exciting and novel material platform for a host of applications. Key to their crystallinity, monolayers of these materials have a 4-fold degenerate structural ground state, and a single energy scale EC (representing the elastic energy required to switch the longer lattice vector along the x- or y-direction) determines how disordered these monolayers are at finite temperature. Disorder arises when nearest neighboring atoms become gently reassigned as the system is thermally excited beyond a critical temperature Tc that is proportional to EC/kB. EC is tunable by chemical composition and it leads to a classification of these materials into two categories: (i) Those for which EC ≥ kBTm, and (ii) those having kBTm > EC ≥ 0, where Tm is a given material's melting temperature. Black phosphorus and SiS monolayers belong to category (i): these materials do not display an intermediate order-disorder transition and melt directly. All other monochalcogenide monolayers with EC > 0 belonging to class (ii) will undergo a two-dimensional transition prior to melting. EC/kB is slightly larger than room temperature for GeS and GeSe, and smaller than 300 K for SnS and SnSe monolayers, so that these materials transition near room temperature. The onset of this generic atomistic phenomena is captured by a planar Potts model up to the order-disorder transition. The order-disorder phase transition in two dimensions described here is at the origin of the Cmcm phase being discussed within the context of bulk layered SnSe.

Random pinning changes the melting scenario of a two-dimensional core-softened potential system

In experiments two-dimensional systems are realized mainly on solid substrates, which introduce quenched disorder due to some inherent defects. The defects of substrates influence the melting scenario of the systems and have to be taken into account in the interpretation of experimental results. We present the results of molecular dynamics simulations of a two-dimensional system with a core-softened potential in which a small fraction of the particles is pinned, inducing quenched disorder. Ppotentials of this type are widely used for the qualitative description of systems with waterlike anomalies. In our previous publications it was shown that the system demonstrates an anomalous melting scenario: at low densities the system melts through two continuous transitions in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory with an intermediate hexatic phase, while at high densities the conventional first-order melting transition takes place. We find that the well-known disorder-induced widening of the hexatic phase occurs at low densities, while in the high-density part of the phase diagram random pinning transforms the first-order melting into two transitions: a continuous KTHNY-like solid-hexatic transition and a first-order hexatic-isotropic liquid transition.

Collective transport properties of driven Skyrmions with random disorder

We use particle-based simulations to examine the static and driven collective phases of Skyrmions interacting with random quenched disorder. We show that nondissipative effects due to the Magnus term reduce the depinning threshold and strongly affect the Skyrmion motion and the nature of the dynamic phases. The quenched disorder causes the Hall angle to become drive dependent in the moving Skyrmion phase, while different flow regimes produce distinct signatures in the transport curves. For weak disorder, the Skyrmions form a pinned crystal and depin elastically, while for strong disorder the system forms a pinned amorphous state that depins plastically. At high drives the Skyrmions can dynamically reorder into a moving crystal, with the onset of reordering determined by the strength of the Magnus term.

Numerical studies on critical properties of the Kosterlitz-Thouless phase for the gauge glass model in two dimensions

The critical exponents are estimated for the gauge glass model in two dimensions, in which only the Kosterlitz-Thouless (KT) phase appears in the low-temperature regime. The nonequilibrium relaxation method is applied to estimate the transition temperature and critical exponents: the static exponent η and the dynamical exponent z. Since the system exhibits criticality in the whole KT phase, we estimate the exponents on the boundary as well as inside the KT phase. The static exponent η depends on both the temperature and the strength of randomness, while the dynamical one z is almost constant throughout the KT phase, including the boundary.

Fluctuations of orientational order and clustering in a two-dimensional colloidal system under quenched disorder

Using both video microscopy of superparamagnetic colloidal particles confined in two dimensions and corresponding computer simulations of repulsive parallel dipoles, we study the formation of fluctuating orientational clusters and topological defects in the context of the KTHNY-like melting scenario under quenched disorder. We analyze cluster densities, average cluster sizes, and the population of noncluster particles, as well as the development of defects, as a function of the system temperature and disorder strength. In addition, the probability distribution of clustering and orientational order is presented. We find that the well-known disorder-induced widening of the hexatic phase can be traced back to the distinct development characteristics of clusters and defects along the melting transitions from the solid phase to the hexatic phase to the isotropic fluid.

Two-dimensional melting under quenched disorder

We study the influence of quenched disorder on the two-dimensional melting behavior of superparamagnetic colloidal particles, using both video microscopy and computer simulations of repulsive parallel dipoles. Quenched disorder is introduced by pinning a fraction of the particles to an underlying substrate. We confirm the occurrence of the Kosterlitz-Thouless-Halperin-Nelson-Young scenario and observe an intermediate hexatic phase. While the fluid-hexatic transition remains largely unaffected by disorder, the hexatic-solid transition shifts to lower temperatures with increasing disorder. This results in a significantly broadened stability range of the hexatic phase. In addition, we observe spatiotemporal critical(like) fluctuations, which are consistent with the continuous character of the phase transitions. Characteristics of first-order transitions are not observed.

Effects of confinement and external fields on structure and transport in colloidal dispersions in reduced dimensionality

In this work, we focus on low-dimensional colloidal model systems, via simulation studies and also some complementary experiments, in order to elucidate the interplay between phase behavior, geometric structures and transport properties. In particular, we try to investigate the (nonlinear!) response of these very soft colloidal systems to various perturbations: uniform and uniaxial pressure, laser fields, shear due to moving boundaries and randomly quenched disorder. We study ordering phenomena on surfaces or in monolayers by Monte Carlo computer simulations of binary hard-disk mixtures, the influence of a substrate being modeled by an external potential. Weak external fields allow a controlled tuning of the miscibility of the mixture. We discuss the laser induced de-mixing for the three different possible couplings to the external potential. The structural behavior of hard spheres interacting with repulsive screened Coulomb or dipolar interaction in 2D and 3D narrow constrictions is investigated using Brownian dynamics simulations. Due to misfits between multiples of the lattice parameter and the channel widths, a variety of ordered and disordered lattice structures have been observed. The resulting local lattice structures and defect probabilities are studied for various cross sections. The influence of a self-organized order within the system is reflected in the velocity of the particles and their diffusive behavior. Additionally, in an experimental system of dipolar colloidal particles confined by gravity on a solid substrate we investigate the effect of pinning on the dynamics of a two-dimensional colloidal liquid. This work contains sections reviewing previous work by the authors as well as new, unpublished results. Among the latter are detailed studies of the phase boundaries of the de-mixing regime in binary systems in external light fields, configurations for shear induced effects at structured walls, studies on the effect of confinement on the structures and defect densities in three-dimensional systems, the effect of confinement and barriers on two-dimensional flow and diffusion, and the effect of pinning sites on the diffusion.

Glassy phases and driven response of the phase-field-crystal model with random pinning

We study the structural correlations and the nonlinear response to a driving force of a two-dimensional phase-field-crystal model with random pinning. The model provides an effective continuous description of lattice systems in the presence of disordered external pinning centers, allowing for both elastic and plastic deformations. We find that the phase-field crystal with disorder assumes an amorphous glassy ground state, with only short-ranged positional and orientational correlations, even in the limit of weak disorder. Under increasing driving force, the pinned amorphous-glass phase evolves into a moving plastic-flow phase and then, finally, a moving smectic phase. The transverse response of the moving smectic phase shows a vanishing transverse critical force for increasing system sizes.

Effects of discreteness on gauge glass models in two and three dimensions

The q-state clock gauge glass model is studied to see the effect of discreteness on the Kosterlitz-Thouless (KT) transition and the ferromagnetic (FM) critical phenomenon in random systems. The nonequilibrium relaxation analysis is applied. In two dimensions, the successive transitions of paramagnetic (PM), KT, and FM phases are investigated along the Nishimori line for q=6, 8, 10, 12, 14, 16, and 1024 (recognized as infinity) cases. For the upper critical temperature, it is found that the transition temperature is almost the same as in the continuous case for all q values. The lower transition temperature is found to be proportional to 1/q2. In three dimensions, the critical behavior of the PM-FM transition is studied along the Nishimori line for q=6, 8, 16, and 1024 cases. It is found that the spin discreteness is irrelevant, and the transition belongs to the same universality class as in the (continuous) XY case.

Disordering transitions and peak effect in polydisperse particle systems

We show numerically that in a binary system of Yukawa particles, a dispersity-driven disordering transition occurs. In the presence of quenched disorder this disordering transition coincides with a marked increase in the depinning threshold, known as a peak effect. We find that the addition of poorly pinned particles can increase the overall pinning in the sample by increasing the amount of topological disorder present. If the quenched disorder is strong enough to create a significant amount of topological disorder in the monodisperse system, addition of a poorly pinned species generates further disorder but does not produce a peak in the depinning force. Our results indicate that for binary mixtures, optimal pinning occurs for topological defect fraction densities from 0.2 to 0.25. For defect densities below this range, the system retains orientational order. We determine the effect of the pinning density, strength, and radius on the depinning peak and find that the peak effect is more pronounced in weakly pinning systems.

Phase transition and critical dynamics in site-diluted Josephson-junction arrays

Measurements of the IV characteristics of site-diluted Josephson-junction arrays have revealed intriguing effects of percolative disorder on the phase transition and the vortex dynamics in a two-dimensional XY system. Different from other types of phase transitions, the Kosterlitz-Thouless transition was eliminated with the introduction of percolative disorder far below the percolation threshold. Even after the Kosterlitz-Thouless order had been removed, the system remained superconducting at low temperatures by establishing a different type of order. Near the percolation threshold, evidence was found that, as a consequence of the underlying fractal structure, the critical dynamics of the phase degrees of freedom persisted down to zero temperature.

Two-dimensional diffusion in the presence of topological disorder

How topological defects affect the dynamics of particles hopping between lattice sites of a distorted, two-dimensional crystal is addressed. Perturbation theory and numerical simulations show that weak, short-ranged topological disorder leads to a finite reduction of the diffusion coefficient. Renormalization group theory and numerical simulations suggest that longer-ranged disorder, such as that from randomly placed dislocations or random disclinations with no net disclinicity, leads to subdiffusion at long times.

Crystal structure in nematic emulsion

We describe the experimental observation of a crystal structure formed by glycerol droplets suspended in a nematic liquid crystal. The structure exhibits a high density hexagonal ordering. We have experimentally observed a noticeable interaction between droplets with tangential boundary conditions. Within the scope of known models we discuss the nature of appropriate mechanisms of the interaction.

Freezing of dynamical exponents in low dimensional random media

A particle in a random potential with logarithmic correlations in dimensions d = 1,2 is shown to undergo a dynamical transition at T(dyn)>0. In d = 1 exact results show T(dyn) = T(c), the static glass transition temperature, and that the dynamical exponent changes from z(T) = 2+2(T(c)/T)(2) at high T to z(T) = 4T(c)/T in the glass phase. The same formulas are argued to hold in d = 2. Dynamical freezing is also predicted in the 2D random gauge XY model and related systems. In d = 1 a mapping between dynamics and statics is unveiled and freezing involves barriers as well as valleys. Anomalous scaling occurs in the creep dynamics, relevant to dislocation motion experiments.

Glass transition of a particle in a random potential, front selection in nonlinear renormalization group, and entropic phenomena in Liouville and sinh-Gordon models

We study via renormalization group (RG), numerics, exact bounds, and qualitative arguments the equilibrium Gibbs measure of a particle in a d-dimensional Gaussian random potential with translationally invariant logarithmic spatial correlations. We show that for any d>/=1 it exhibits a transition at T=T(c)>0. The low-temperature glass phase has a nontrivial structure, being dominated by a few distant states (with replica symmetry breaking phenomenology). In finite dimension this transition exists only in this "marginal glass" case (energy fluctuation exponent straight theta=0) and disappears if correlations grow faster (single ground-state dominance straight theta>0) or slower (high-temperature phase). The associated extremal statistics problem for correlated energy landscapes exhibits universal features which we describe using a nonlinear Kolmogorov (KPP) RG equation. These include the tails of the distribution of the minimal energy (or free energy) and the finite-size corrections, which are universal. The glass transition is closely related to Derrida's random energy models. In d=2, the connection between this problem and Liouville and sinh-Gordon models is discussed. The glass transition of the particle exhibits interesting similarities with the weak- to strong-coupling transition in Liouville (c=1 barrier) and with a transition that we conjecture for the sinh-Gordon model, with correspondence in some exact results and RG analysis. Glassy freezing of the particle is associated with the generation under RG of new local operators and of nonsmooth configurations in Liouville. Applications to Dirac fermions in random magnetic fields at criticality reveal a peculiar "quasilocalized" regime (corresponding to the glass phase for the particle), where eigenfunctions are concentrated over a finite number of distant regions, and allow us to recover the multifractal spectrum in the delocalized regime.