Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTwo-dimensional melting in simple atomic systems: continuous vs. discontinous melting
358
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate the characteristics of two dimensional melting in simple atomic systems via isobaric-isothermal ($NPT$) and isochoric-isothermal ($NVT$) molecular dynamics simulations with special focus on the effect of the range of the potential on the melting. We find that the system with interatomic potential of longer range clearly exhibits a region (in the $PT$ plane) of (thermodynamically) stable hexatic phase. On the other hand, the one with shorter range potential exhibits a first-order melting transition both in $NPT$ and $NVT$ ensembles. Melting of the system with intermediate range potential shows a hexatic-like feature near the melting transition in $NVT$ ensemble, but it undergoes an unstable hexatic-like phase during melting process in $NPT$ ensemble, which implies existence of a weakly first order transition. The overall features represent a crossover from a continuous melting transition in the cases of longer-ranged potential to a discontinuous (first order) one in the systems with shorter and intermediate ranged potential. We also calculate the Binder cumulants as well as the susceptibility of the bond-orientational order parameter.
rate research
Read More
Motivated by the recently developed duality between elasticity of a crystal and a symmetric tensor gauge theory by Pretko and Radzihovsky, we explore its classical analog, that is a dual theory of the dislocation-mediated melting of a two-dimensional crystal, formulated in terms of a higher derivative vector sine-Gordon model. It provides a transparent description of the continuous two-stage melting in terms of the renormalization-group relevance of two cosine operators that control the sequential unbinding of dislocations and disclinations, respectively corresponding to the crystal-to-hexatic and hexatic-to-isotropic fluid transitions. This renormalization-group analysis compactly reproduces seminal results of the Coulomb gas description, such as the flows of the elastic couplings and of the dislocation and disclination fugacities, as well the temperature dependence of the associated correlation lengths.
In small confined systems predictions for the melting point strongly depend on the choice of quantity and on the way it is computed, even yielding divergent and ambiguous results. We present a very simple quantity which allows to control these problems -- the variance of the block averaged interparticle distance fluctuations.
The size and shape of a large variety of polymeric particles, including biological cells, star polymers, dendrimes, and microgels, depend on the applied stresses as the particles are extremely soft. In high-density suspensions these particles deform as stressed by their neighbors, which implies that the interparticle interaction becomes of many-body type. Investigating a two-dimensional model of cell tissue, where the single particle shear modulus is related to the cell adhesion strength, here we show that the particle deformability affects the melting scenario. On increasing the temperature, stiff particles undergo a first-order solid/liquid transition, while soft ones undergo a continuous solid/hexatic transition followed by a discontinuous hexatic/liquid transition. At zero temperature the melting transition driven by the decrease of the adhesion strength occurs through two continuous transitions as in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario. Thus, there is a range of adhesion strength values where the hexatic phase is stable at zero temperature, which suggests that the intermediate phase of the epithelial-to-mesenchymal transition could be hexatic type.
Two-dimensional systems may admit a hexatic phase and hexatic-liquid transitions of different natures. The determination of their phase diagrams proved challenging, and indeed those of hard-disks, hard regular polygons, and inverse power-law potentials, have been only recently clarified. In this context, the role of attractive forces is currently speculative, despite their prevalence at both the molecular and colloidal scale. Here we demonstrate, via numerical simulations, that attraction promotes a discontinuous melting scenario with no hexatic phase. At high-temperature, Lennard-Jones particles and attractive polygons follow the shape-dominated melting scenario observed in hard-disks and hard polygons, respectively. Conversely, all systems melt via a first-order transition with no hexatic phase at low temperature, where attractive forces dominate. The intermediate temperature melting scenario is shape-dependent. Our results suggest that, in colloidal experiments, the tunability of the strength of the attractive forces allows for the observation of different melting scenario in the same system.
The phase diagram of two-dimensional continuous particle systems is studied using Event-Chain Monte Carlo. For soft disks with repulsive power-law interactions $propto r^{-n}$ with $n gtrsim 6$, the recently established hard-disk melting scenario ($n to infty$) holds: a first-order liquid-hexatic and a continuous hexatic-solid transition are identified. Close to $n = 6$, the coexisting liquid exhibits very long orientational correlations, and positional correlations in the hexatic are extremely short. For $nlesssim 6$, the liquid-hexatic transition is continuous, with correlations consistent with the Kosterlitz-Thouless-Halperin-Nelson-Yong (KTHNY) scenario. To illustrate the generality of these results, we demonstrate that Yukawa particles likewise may follow either the KTHNY or the hard-disk melting scenario, depending on the Debye-Huckel screening length as well as on the temperature.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
