1
|
Milchev A, Binder K. Surface enrichment and interdiffusion in blends of semiflexible polymers of different stiffness. SOFT MATTER 2022; 18:3781-3792. [PMID: 35514321 DOI: 10.1039/d2sm00036a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A model for a mixture of two kinds of semiflexible polymers (A and B) with the same chain length (NA = NB = 32), but different persistence lengths, confined between parallel planar repulsive walls in a common good solvent is studied by molecular dynamics simulations. In the isotropic phase at low polymer concentrations, both polymers are repelled by the walls, and the system is anisotropic near the walls over a range controlled by the polymer linear dimensions. Close to the concentrations where in the bulk nematic order sets in, precursors of thick nematic layers at the walls are observed, strongly enriched by a stiffer component, which hence is depleted in the center of the slit pore. At larger concentrations, where in the bulk a uniformly mixed nematic phase occurs, the enrichment of B-chains at the walls is rather minor, extending over the scale of the transverse correlation length of concentration fluctuations, which is of the order of a few monomeric diameters only for the present model. In this ordered phase, both self-diffusion and interdiffusion of chains (in the direction perpendicular to the director) are found to be significantly slowed down in comparison to dilute solutions. Since equilibration times scale with the square of the slit thickness, incomplete equilibration is predicted when polymeric coatings on substrate containing polymers differing in stiffness are produced.
Collapse
Affiliation(s)
- Andrey Milchev
- Institute for Physical Chemistry, Bulgarian Academia of Sciences, 1113, Sofia, Bulgaria.
| | - Kurt Binder
- Institute of Physics, Johannes Gutenberg University Mainz, Staudingerweg 7, 55128 Mainz, Germany
| |
Collapse
|
2
|
Takemoto K, Ishii Y, Washizu H, Kim K, Matubayasi N. Simulating the nematic-isotropic phase transition of liquid crystal model via generalized replica-exchange method. J Chem Phys 2022; 156:014901. [PMID: 34998348 DOI: 10.1063/5.0073105] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in the gREM, allowing for the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.
Collapse
Affiliation(s)
- Kengo Takemoto
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Yoshiki Ishii
- Graduate School of Information Science, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Hitoshi Washizu
- Graduate School of Information Science, University of Hyogo, Kobe, Hyogo 650-0047, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
3
|
Kowaguchi A, Brumby PE, Yasuoka K. Phase Transitions and Hysteresis for a Simple Model Liquid Crystal by Replica-Exchange Monte Carlo Simulations. Molecules 2021; 26:1421. [PMID: 33808039 PMCID: PMC7961897 DOI: 10.3390/molecules26051421] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Revised: 03/01/2021] [Accepted: 03/01/2021] [Indexed: 11/17/2022] Open
Abstract
In this work, the advantages of applying the temperature and pressure replica-exchange method to investigate the phase transitions and the hysteresis for liquid-crystal fluids were demonstrated. In applying this method to the commonly used Hess-Su liquid-crystal model, heat capacity peaks and points of phase co-existence were observed. The absence of a smectic phase at higher densities and a narrow range of the nematic phase were reported. The identity of the crystalline phase of this system was found to a hexagonal close-packed solid. Since the nematic-solid phase transition is strongly first order, care must be taken when using this model not to inadvertently simulate meta-stable nematic states at higher densities. In further analysis, the Weighted Histogram Analysis Method was applied to verify the precise locations of the phase transition points.
Collapse
Affiliation(s)
| | - Paul E. Brumby
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan;
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan;
| |
Collapse
|
4
|
Skutnik RA, Eichler JC, Mazza MG, Schoen M. The temperature dependence of the helical pitch in a cholesteric liquid crystal. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1881638] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Robert A. Skutnik
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakulät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Berlin, Germany
| | - Jan-Christoph Eichler
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakulät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Berlin, Germany
| | - Marco G. Mazza
- Interdisciplinary Centre for Mathematical Modelling and Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire, UK
- Max Planck Institute for Dynamics and Self-Organization, Göttingen, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakulät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Berlin, Germany
- Department of Chemical Engineering, Imperial College London, London, UK
| |
Collapse
|
5
|
Geier IS, Wandrei SM, Skutnik RA, Schoen M. Molecular theory of a ferromagnetic nematic liquid crystal. Phys Rev E 2019; 100:022702. [PMID: 31574629 DOI: 10.1103/physreve.100.022702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2019] [Indexed: 06/10/2023]
Abstract
We employ a version of classical density functional theory to study the phase behavior of a simple model liquid crystal in an external field. The uniaxially symmetric molecules have a spherically symmetric core with superimposed orientation-dependent attractions. The interaction between the cores consists of a hard-sphere repulsion plus an isotropic square-well attraction. The anisotropic part of the interaction potential allows for the formation of a uniaxially symmetric nematic phase. The orientation of the molecules couples to an external polar field. The external field is capable of rotating the nematic director n[over ̂] in the x-z plane. The field is also capable of changing the topology of the phase diagram in that it suppresses the phase coexistence between an isotropic liquid and a nematic phase observed in the absence of the field. We study the transition from an unpolar to a polar nematic phase in terms of the orientation-distribution function (odf), nematic and polar order parameters, and components of n[over ̂]. If represented suitably the odf allows us to study orientational changes during the switching process between nonpolar and polar nematic phases. We also give a simple argument that explains why nematic order is lost whereas polar order persists up to the gas-liquid critical point along the coexistence curve. We also discuss the relevance of our theory for future experimental studies.
Collapse
Affiliation(s)
- Immanuel S Geier
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Stefanie M Wandrei
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Robert A Skutnik
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
- Department of Chemical and Biomolecular Engineering, Engineering Building I, Box 7905, North Carolina State University, 911 Partners Way, Raleigh, North Carolina 27695, USA
| |
Collapse
|
6
|
Affiliation(s)
- Michael P. Allen
- Department of Physics, University of Warwick, Coventry, UK
- H. H. Wills Physics Laboratory, Royal Fort, Bristol, UK
| |
Collapse
|
7
|
Skutnik RA, Lehmann L, Püschel-Schlotthauer S, Jackson G, Schoen M. The formation of biaxial nematic phases in binary mixtures of thermotropic liquid-crystals composed of uniaxial molecules. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1581292] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Robert A. Skutnik
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Berlin, Germany
| | - Louis Lehmann
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Berlin, Germany
| | - Sergej Püschel-Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Berlin, Germany
| | - George Jackson
- Department of Chemical Engineering, South Kensington Campus, Imperial College London, London, UK
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Berlin, Germany
- Department of Chemical Engineering, South Kensington Campus, Imperial College London, London, UK
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
8
|
Zhao Y, Stratt RM. Measuring order in disordered systems and disorder in ordered systems: Random matrix theory for isotropic and nematic liquid crystals and its perspective on pseudo-nematic domains. J Chem Phys 2018; 148:204501. [PMID: 29865812 DOI: 10.1063/1.5024678] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Surprisingly long-ranged intermolecular correlations begin to appear in isotropic (orientationally disordered) phases of liquid crystal forming molecules when the temperature or density starts to close in on the boundary with the nematic (ordered) phase. Indeed, the presence of slowly relaxing, strongly orientationally correlated, sets of molecules under putatively disordered conditions ("pseudo-nematic domains") has been apparent for some time from light-scattering and optical-Kerr experiments. Still, a fully microscopic characterization of these domains has been lacking. We illustrate in this paper how pseudo-nematic domains can be studied in even relatively small computer simulations by looking for order-parameter tensor fluctuations much larger than one would expect from random matrix theory. To develop this idea, we show that random matrix theory offers an exact description of how the probability distribution for liquid-crystal order parameter tensors converges to its macroscopic-system limit. We then illustrate how domain properties can be inferred from finite-size-induced deviations from these random matrix predictions. A straightforward generalization of time-independent random matrix theory also allows us to prove that the analogous random matrix predictions for the time dependence of the order-parameter tensor are similarly exact in the macroscopic limit, and that relaxation behavior of the domains can be seen in the breakdown of the finite-size scaling required by that random-matrix theory.
Collapse
Affiliation(s)
- Yan Zhao
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Richard M Stratt
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
9
|
Schoen M, Haslam AJ, Jackson G. Perturbation Theory versus Thermodynamic Integration. Beyond a Mean-Field Treatment of Pair Correlations in a Nematic Model Liquid Crystal. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11345-11365. [PMID: 28772076 DOI: 10.1021/acs.langmuir.7b01849] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The phase behavior and structure of a simple square-well bulk fluid with anisotropic interactions is described in detail. The orientation dependence of the intermolecular interactions allows for the formation of a nematic liquid-crystalline phase in addition to the more conventional isotropic gas and liquid phases. A version of classical density functional theory (DFT) is employed to determine the properties of the model, and comparisons are made with the corresponding data from Monte Carlo (MC) computer simulations in both the grand canonical and canonical ensembles, providing a benchmark to assess the adequacy of the DFT results. A novel element of the DFT approach is the assumption that the structure of the fluid is dominated by intermolecular interactions in the isotropic fluid. A so-called augmented modified mean-field (AMMF) approximation is employed accounting for the influence of anisotropic interactions. The AMMF approximation becomes exact in the limit of vanishing density. We discuss advantages and disadvantages of the AMMF approximation with respect to an accurate description of isotropic and nematic branches of the phase diagram, the degree of orientational order, and orientation-dependent pair correlations. The performance of the AMMF approximations is found to be good in comparison with the MC data; the AMMF approximation has clear advantages with respect to an accurate and more detailed description of the fluid structure. Possible strategies to improve the DFT are discussed.
Collapse
Affiliation(s)
- Martin Schoen
- Department of Chemical Engineering, Imperial College London , South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Andrew J Haslam
- Department of Chemical Engineering, Imperial College London , South Kensington Campus, London SW7 2AZ, United Kingdom
| | - George Jackson
- Department of Chemical Engineering, Imperial College London , South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
10
|
Püschel-Schlotthauer S, Meiwes Turrión V, Hall CK, Mazza MG, Schoen M. The Impact of Colloidal Surface-Anchoring on the Smectic A Phase. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:2222-2234. [PMID: 28165243 DOI: 10.1021/acs.langmuir.6b03941] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Liquid-crystalline phases are known for their unique properties, i.e., the combination of fluidity and long-range orientational and/or positional order. The presence of a colloidal particle gives rise to perturbations of this order locally. These perturbations are the origin of intercolloidal forces driving the colloidal self-assembly in a directed manner. Hence, the understanding of these perturbations is the first step in understanding and controlling the self-assembly process. Here, we perform Monte Carlo simulations to investigate the perturbations of orientational and positional order in a smectic A phase caused by a spherical colloid. We model the host phase via an interaction potential that reproduces characteristic features of phase behavior, structure, dynamics, and elasticity [S. Püschel-Schlotthauer et al. J. Chem. Phys. 2016, 145, 164903]. For strong homeotropic anchoring conditions, we find a Saturn ring defect and an onion structure in the smectic A phase; the latter has never been reported for colloids so far. For strong planar anchoring conditions, we find Boojum defects that become elongated at low temperature, similar to what is observed in experiments. However, for weak planar anchoring conditions, a double surface ring defect is exhibited in the smectic A phase.
Collapse
Affiliation(s)
- Sergej Püschel-Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin , Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Victor Meiwes Turrión
- Max Planck Institute for Dynamics and Self-Organization (MPIDS) , Am Faßberg 17, 37077 Göttingen, Germany
| | - Carol K Hall
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Engineering Building I, Box 7905, 911 Partners Way, Raleigh, North Carolina 27695, United States
| | - Marco G Mazza
- Max Planck Institute for Dynamics and Self-Organization (MPIDS) , Am Faßberg 17, 37077 Göttingen, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin , Straße des 17. Juni 135, 10623 Berlin, Germany
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Engineering Building I, Box 7905, 911 Partners Way, Raleigh, North Carolina 27695, United States
| |
Collapse
|
11
|
Püschel-Schlotthauer S, Meiwes Turrión V, Stieger T, Grotjahn R, Hall CK, Mazza MG, Schoen M. A novel model for smectic liquid crystals: Elastic anisotropy and response to a steady-state flow. J Chem Phys 2016; 145:164903. [PMID: 27802618 DOI: 10.1063/1.4965711] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
By means of a combination of equilibrium Monte Carlo and molecular dynamics simulations and nonequilibrium molecular dynamics we investigate the ordered, uniaxial phases (i.e., nematic and smectic A) of a model liquid crystal. We characterize equilibrium behavior through their diffusive behavior and elastic properties. As one approaches the equilibrium isotropic-nematic phase transition, diffusion becomes anisotropic in that self-diffusion D⊥ in the direction orthogonal to a molecule's long axis is more hindered than self-diffusion D∥ in the direction parallel to that axis. Close to nematic-smectic A phase transition the opposite is true, D∥ < D⊥. The Frank elastic constants K1, K2, and K3 for the respective splay, twist, and bend deformations of the director field n̂ are no longer equal and exhibit a temperature dependence observed experimentally for cyanobiphenyls. Under nonequilibrium conditions, a pressure gradient applied to the smectic A phase generates Poiseuille-like or plug flow depending on whether the convective velocity is parallel or orthogonal to the plane of smectic layers. We find that in Poiseuille-like flow the viscosity of the smectic A phase is higher than in plug flow. This can be rationalized via the velocity-field component in the direction of the flow. In a sufficiently strong flow these smectic layers are not destroyed but significantly bent.
Collapse
Affiliation(s)
- Sergej Püschel-Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Victor Meiwes Turrión
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), Am Faßberg 17, 37077 Göttingen, Germany
| | - Tillmann Stieger
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Robin Grotjahn
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Carol K Hall
- Department of Chemical and Biomolecular Engineering, Engineering Building I, Box 7905, North Carolina State University, 911 Partners Way, Raleigh, North Carolina 27695, USA
| | - Marco G Mazza
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), Am Faßberg 17, 37077 Göttingen, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
12
|
Püschel-Schlotthauer S, Stieger T, Melle M, Mazza MG, Schoen M. Coarse-grained treatment of the self-assembly of colloids suspended in a nematic host phase. SOFT MATTER 2016; 12:469-480. [PMID: 26477506 DOI: 10.1039/c5sm01860a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The complex interplay of molecular scale effects, nonlinearities in the orientational field and long-range elastic forces makes liquid-crystal physics very challenging. A consistent way to extract information from the microscopic, molecular scale up to the meso- and macroscopic scale is still missing. Here, we develop a hybrid procedure that bridges this gap by combining extensive Monte Carlo (MC) simulations, a local Landau-de Gennes theory, classical density functional theory, and finite-size scaling theory. As a test case to demonstrate the power and validity of our novel approach we study the effective interaction among colloids with Boojum defect topology immersed in a nematic liquid crystal. In particular, at sufficiently small separations colloids attract each other if the angle between their center-of-mass distance vector and the far-field nematic director is about 30°. Using the effective potential in coarse-grained two-dimensional MC simulations we show that self-assembled structures formed by the colloids are in excellent agreement with experimental data.
Collapse
Affiliation(s)
- Sergej Püschel-Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | | | | | | |
Collapse
|
13
|
Stieger T, Püschel-Schlotthauer S, Schoen M, Mazza MG. Flow-induced deformation of closed disclination lines near a spherical colloid immersed in a nematic host phase. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1096973] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Tillmann Stieger
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| | - Sergej Püschel-Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
- Department of Chemical and Biomolecular Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA
| | - Marco G. Mazza
- Max-Planck-Institut für Dynamik und Selbstorganisation, Am Faßberg 17, 37077 Göttingen, Germany
| |
Collapse
|
14
|
Schlotthauer S, Skutnik RA, Stieger T, Schoen M. Defect topologies in chiral liquid crystals confined to mesoscopic channels. J Chem Phys 2015; 142:194704. [DOI: 10.1063/1.4920979] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sergej Schlotthauer
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, Berlin 10623, Germany
| | - Robert A. Skutnik
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, Berlin 10623, Germany
| | - Tillmann Stieger
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, Berlin 10623, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, Berlin 10623, Germany
- Department of Chemical and Biomolecular Engineering, North Carolina State University, 911 Partners Way, Raleigh, North Carolina 27695, USA
| |
Collapse
|
15
|
Zhang CY, Jian XL, Lu W. Structure and percolation of one-patch spherocylinders. SOFT MATTER 2015; 11:1362-1368. [PMID: 25575168 DOI: 10.1039/c4sm02402h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
When the volume fraction exceeds the threshold, the colloidal particles would form a spanning cluster to realize percolation, which is affected by the shape of the particles, interaction between particles, etc. In this paper, we use the Monte Carlo method to study the structure and percolation of a system of one-patch spherocylinders which have been fabricated recently [Chaudhary et al., J. Am. Chem. Soc., 2012, 134, 12901]. With strong adsorption, one-patch spherocylinders self-assemble into multipods which further make contact with each other to form a percolation network at a high volume fraction, while the percolation network is inhibited by the local structures in a system of one-patch spheres. The main multipods are dipods when the patch angle equals π/3, while they are tetrapods and pentapods when the patch angle equals 2π/3. With enhancing the adsorption, the bigger the patch angle, the more the percolation threshold drops. The orientational order parameter, the distribution of the relative orientation between the nearest neighbors and the probabilities of a spherocylinder owning n adsorbing neighbors have been calculated to analyze the formation and transition of the structures.
Collapse
Affiliation(s)
- Cheng-Yu Zhang
- College of Engineering, Nanjing Agricultural University, Nanjing 210031, China.
| | | | | |
Collapse
|
16
|
Stieger T, Schoen M, Mazza MG. Effects of flow on topological defects in a nematic liquid crystal near a colloid. J Chem Phys 2014; 140:054905. [PMID: 24511977 DOI: 10.1063/1.4862953] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
We perform molecular dynamics simulations of a nematic liquid crystal flowing around a colloidal particle. We study the flow-induced modifications of the topological defects in the liquid crystal due to the presence of the colloid. We show that flow distorts Boojum defects into an asymmetrically larger downstream lobe, and that Saturn ring defects are convected downstream along the flow direction, which is in agreement with experimental observations. Additionally, for a Janus colloid with both parallel and perpendicular patches, exhibiting a Boojum defect and a Saturn ring defect, we find that the Boojum defect facing the upstream direction is destroyed and the Saturn ring is convected downstream.
Collapse
Affiliation(s)
- Tillmann Stieger
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Marco G Mazza
- Max-Planck-Institut für Dynamik und Selbstorganisation (MPIDS), Am Faßberg 17, 37077 Göttingen, Germany
| |
Collapse
|
17
|
Melle M, Schlotthauer S, Hall CK, Diaz-Herrera E, Schoen M. Disclination lines at homogeneous and heterogeneous colloids immersed in a chiral liquid crystal. SOFT MATTER 2014; 10:5489-5502. [PMID: 24954626 DOI: 10.1039/c4sm00959b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
In the present work we perform Monte Carlo simulations in the isothermal-isobaric ensemble to study defect topologies formed in a cholesteric liquid crystal due to the presence of a spherical colloidal particle. Topological defects arise because of the competition between anchoring at the colloidal surface and the local director. We consider homogeneous colloids with either local homeotropic or planar anchoring to validate our model by comparison with earlier lattice Boltzmann studies. Furthermore, we perform simulations of a colloid in a twisted nematic cell and discuss the difference between induced and intrinsic chirality on the formation of topological defects. We present a simple geometrical argument capable of describing the complex three-dimensional topology of disclination lines evolving near the surface of the colloid. The presence of a Janus colloid in a cholesteric host fluid reveals a rich variety of defect structures. Using the Frank free energy we analyze these defects quantitatively indicating a preferred orientation of the Janus colloid relative to the cholesteric helix.
Collapse
Affiliation(s)
- Michael Melle
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany.
| | | | | | | | | |
Collapse
|
18
|
Giura S, Schoen M. Density-functional theory and Monte Carlo simulations of the phase behavior of a simple model liquid crystal. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:022507. [PMID: 25215749 DOI: 10.1103/physreve.90.022507] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2014] [Indexed: 06/03/2023]
Abstract
We consider the phase behavior of a simple model of a liquid crystal by means of modified mean-field density-functional theory (MMF DFT) and Monte Carlo simulations in the grand canonical ensemble (GCEMC). The pairwise additive interactions between liquid-crystal molecules are modeled via a Lennard-Jones potential in which the attractive contribution depends on the orientation of the molecules. We derive the form of this orientation dependence through an expansion in terms of rotational invariants. Our MMF DFT predicts two topologically different phase diagrams. At weak to intermediate coupling of the orientation dependent attraction, there is a discontinuous isotropic-nematic liquid-liquid phase transition in addition to the gas-isotropic liquid one. In the limit of strong coupling, the gas-isotropic liquid critical point is suppressed in favor of a fluid- (gas- or isotropic-) nematic phase transition which is always discontinuous. By considering three representative isotherms in parallel GCEMC simulations, we confirm the general topology of the phase diagram predicted by MMF DFT at intermediate coupling strength. From the combined MMF DFT-GCEMC approach, we conclude that the isotropic-nematic phase transition is very weakly first order, thus confirming earlier computer simulation results for the same model [see M. Greschek and M. Schoen, Phys. Rev. E 83, 011704 (2011)].
Collapse
Affiliation(s)
- Stefano Giura
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany and Department of Chemical and Biomolecular Engineering, Engineering Building I, Box 7905, North Carolina State University, 911 Partners Way, Raleigh, North Carolina 27695, USA
| |
Collapse
|
19
|
Schoen M, Giura S, Klapp SHL. Phase behavior of an amphiphilic fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:012310. [PMID: 24580230 DOI: 10.1103/physreve.89.012310] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2013] [Indexed: 06/03/2023]
Abstract
We invoke mean-field density functional theory (DFT) to investigate the phase behavior of an amphiphilic fluid composed of a hard-sphere core plus a superimposed anisometric Lennard-Jones perturbation. The orientation dependence of the interactions consists of a contribution analogous to the interaction potential between a pair of "spins" in the classical, three-dimensional Heisenberg fluid and another one reminiscent of the interaction between (electric or magnetic) point dipoles. At fixed orientation both contributions are short-range in nature decaying as r-6 (r being the separation between the centers of mass of a pair of amphiphiles). Based upon two mean-field-like approximations for the pair correlation function that differ in the degree of sophistication we derive expressions for the phase boundaries between various isotropic and polar phases that we solve numerically by the Newton-Raphson method. For sufficiently strong coupling between the Heisenberg "spins" both mean-field approximations generate three topologically different and generic types of phase diagrams that are observed in agreement with earlier work [see, for example, Tavares et al., Phys. Rev. E 52, 1915 (1995)]. Whereas the dipolar contribution alone is incapable of stabilizing polar phases on account of its short-range nature it is nevertheless important for details of the phase diagram such as location of the gas-isotropic liquid critical point, triple, and tricritical points. By tuning the dipolar coupling constant suitably one may, in fact, switch between topologically different phase diagrams. Employing also Monte Carlo simulations in the isothermal-isobaric ensemble the general topology of the DFT phase diagrams is confirmed.
Collapse
Affiliation(s)
- Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany and Department of Chemical and Biomolecular Engineering, Engineering Building I, Box 7905, North Carolina State University, 911 Partners Way, Raleigh, North Carolina 27695, USA
| | - Stefano Giura
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 115, 10623 Berlin, Germany
| | - Sabine H L Klapp
- Institut für Theoretische Physik, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Hardenbergstr. 36, 10623 Berlin, Germany
| |
Collapse
|
20
|
Melle M, Theile M, Hall CK, Schoen M. Nanoconfinement-induced structures in chiral liquid crystals. Int J Mol Sci 2013; 14:17584-607. [PMID: 23989605 PMCID: PMC3794743 DOI: 10.3390/ijms140917584] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2013] [Revised: 08/08/2013] [Accepted: 08/16/2013] [Indexed: 01/15/2023] Open
Abstract
We employ Monte Carlo simulations in a specialized isothermal-isobaric and in the grand canonical ensemble to study structure formation in chiral liquid crystals as a function of molecular chirality. Our model potential consists of a simple Lennard-Jones potential, where the attractive contribution has been modified to represent the orientation dependence of the interaction between a pair of chiral liquid-crystal molecules. The liquid crystal is confined between a pair of planar and atomically smooth substrates onto which molecules are anchored in a hybrid fashion. Hybrid anchoring allows for the formation of helical structures in the direction perpendicular to the substrate plane without exposing the helix to spurious strains. At low chirality, we observe a cholesteric phase, which is transformed into a blue phase at higher chirality. More specifically, by studying the unit cell and the spatial arrangement of disclination lines, this blue phase can be established as blue phase II. If the distance between the confining substrates and molecular chirality are chosen properly, we see a third structure, which may be thought of as a hybrid, exhibiting mixed features of a cholesteric and a blue phase.
Collapse
Affiliation(s)
- Michael Melle
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Straße des 17. Juni 135, Berlin 10623, Germany; E-Mails: (M.T.); (M.S.)
| | - Madlona Theile
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Straße des 17. Juni 135, Berlin 10623, Germany; E-Mails: (M.T.); (M.S.)
| | - Carol K. Hall
- Department of Chemical and Biomolecular Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA; E-Mail: (C.K.H.)
| | - Martin Schoen
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Technische Universität Berlin, Straße des 17. Juni 135, Berlin 10623, Germany; E-Mails: (M.T.); (M.S.)
- Department of Chemical and Biomolecular Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA; E-Mail: (C.K.H.)
| |
Collapse
|
21
|
Giura S, Márkus BG, Klapp SHL, Schoen M. Isotropic-polar phase transitions in an amphiphilic fluid: density functional theory versus computer simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:012313. [PMID: 23410337 DOI: 10.1103/physreve.87.012313] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2012] [Indexed: 06/01/2023]
Abstract
We investigate the critical line separating isotropic from polar phases in an amphiphilic bulk fluid by means of density functional theory (DFT) and Monte Carlo (MC) simulations in the isothermal-isobaric ensemble. The intermolecular interactions are described by a Lennard-Jones potential in which the attractive contribution is modified by an orientation-dependent function. The latter consists of two terms: The first one has the orientation dependence of a classical three-dimensional Heisenberg interaction, whereas, the second one has the orientation dependence of a classical dipole-dipole interaction. However, both contributions are short range. Employing DFT together with a modified mean-field (MMF) approximation for the orientation-dependent pair correlation function, we derive an analytical expression for the critical line separating isotropic from polar liquidlike phases. In parallel MC simulations, we locate the line of critical points through an analysis of Binder's second-order cumulant of the polar-order parameter. Comparison with DFT shows that the dipolelike contribution is irrelevant for the isotropic-polar phase transition. As far as the Heisenberg contribution is concerned, the MC data are in semiquantitative agreement with the DFT predictions for sufficiently strong coupling between molecular orientations. For weaker coupling, the variation in the ratio of critical density and temperature ρ(c)/T(c) with the Heisenberg coupling constant ε(H) is underestimated by the MMF treatment. The MC results suggest that this is because ρ(c) increases with decreasing ε(H) such that the assumption on which the MMF approach rests becomes less applicable in the weaker-coupling limit.
Collapse
Affiliation(s)
- Stefano Giura
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Strasse des 17. Juni 115, Berlin 10623, Germany
| | | | | | | |
Collapse
|
22
|
Cartwright JHE, Piro O, Sánchez PA, Sintes T. Ice polyamorphism in the minimal Mercedes-Benz model of water. J Chem Phys 2012; 137:244503. [PMID: 23277941 DOI: 10.1063/1.4772202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate ice polyamorphism in the context of the two-dimensional Mercedes-Benz model of water. We find a first-order phase transition between a crystalline phase and a high-density amorphous phase. Furthermore, we find a reversible transformation between two amorphous structures of high and low density; however, we find this to be a continuous and not an abrupt transition, as the low-density amorphous phase does not show structural stability. We discuss the origin of this behavior and its implications with regard to the minimal generic modeling of polyamorphism.
Collapse
Affiliation(s)
- Julyan H E Cartwright
- Instituto Andaluz de Ciencias de la Tierra, CSIC-Universidad de Granada, E-18071 Granada, Spain.
| | | | | | | |
Collapse
|
23
|
Greschek M, Gubbins KE, Schoen M. Imprinting substrate structures onto a nematic liquid crystal. J Chem Phys 2012; 137:144703. [DOI: 10.1063/1.4757391] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
|
24
|
Melle M, Schlotthauer S, Mazza MG, Klapp SHL, Schoen M. Defect topologies in a nematic liquid crystal near a patchy colloid. J Chem Phys 2012; 136:194703. [PMID: 22612105 DOI: 10.1063/1.4717619] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Using isothermal-isobaric Monte Carlo simulations we investigate defect topologies due to a spherical colloidal particle immersed in a nematic liquid crystal. Defects arise because of the competition between the preferential orientation at the colloid's surface and the far-field director ̂n(0). Considering a chemically homogeneous colloid as a special case we observe the well-known surface and saturn ring defect topologies for weak and strong perpendicular anchoring, respectively; for homogeneous, strong parallel anchoring we find a boojum defect topology that has been seen experimentally [see P. Poulin and D. A. Weitz, Phys. Rev. E 57, 626 (1998)] but not in computer simulations. We also consider a heterogeneous, patchy colloid where the liquid-crystal molecules anchor either preferentially planar or perpendicular at the surface of the colloid. For a patchy colloid we observe a boojum ring defect topology in agreement with recent experimental studies [see M. Conradi, M. Ravnik, M. Bele, M. Zorko, S. Žumer, and I. Muševič, Soft Matter 5, 3905 (2009)]. We also observe two other novel defect topologies that have not been reported thus far neither experimentally nor theoretically.
Collapse
Affiliation(s)
- Michael Melle
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Strasse des 17. Juni 135, 10623 Berlin, Germany
| | | | | | | | | |
Collapse
|
25
|
Diffusivity maximum in a reentrant nematic phase. Int J Mol Sci 2012; 13:7854-7871. [PMID: 22837730 PMCID: PMC3397562 DOI: 10.3390/ijms13067854] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2012] [Revised: 06/11/2012] [Accepted: 06/13/2012] [Indexed: 12/04/2022] Open
Abstract
We report molecular dynamics simulations of confined liquid crystals using the Gay–Berne–Kihara model. Upon isobaric cooling, the standard sequence of isotropic–nematic–smectic A phase transitions is found. Upon further cooling a reentrant nematic phase occurs. We investigate the temperature dependence of the self-diffusion coefficient of the fluid in the nematic, smectic and reentrant nematic phases. We find a maximum in diffusivity upon isobaric cooling. Diffusion increases dramatically in the reentrant phase due to the high orientational molecular order. As the temperature is lowered, the diffusion coefficient follows an Arrhenius behavior. The activation energy of the reentrant phase is found in reasonable agreement with the reported experimental data. We discuss how repulsive interactions may be the underlying mechanism that could explain the occurrence of reentrant nematic behavior for polar and non-polar molecules.
Collapse
|
26
|
Daoulas KC, Rühle V, Kremer K. Simulations of nematic homopolymer melts using particle-based models with interactions expressed through collective variables. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:284121. [PMID: 22738833 DOI: 10.1088/0953-8984/24/28/284121] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
We develop a hybrid Monte Carlo approach for modelling nematic liquid crystals of homopolymer melts. The polymer architecture is described with a discrete worm-like chain model. A quadratic density functional accounts for the limited compressibility of the liquid, while an additional quadratic functional of the local orientation tensor of the segments captures the nematic ordering. The approach can efficiently address large systems parametrized according to volumetric and conformational properties, representative of real polymeric materials. The results of the simulations regarding the influence of the molecular weight on the isotropic-nematic transition are compared to predictions from a Landau-de Gennes free energy expansion. The formation of the nematic phase is addressed within Rouse-like dynamics, realized using the current model.
Collapse
Affiliation(s)
- Kostas Ch Daoulas
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| | | | | |
Collapse
|
27
|
Greschek M, Schoen M. Orientational prewetting of planar solid substrates by a model liquid crystal. J Chem Phys 2012; 135:204702. [PMID: 22128947 DOI: 10.1063/1.3660377] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present grand canonical ensemble Monte Carlo simulations of prewetting transitions in a model liquid crystal at structureless solid substrates. Molecules of the liquid crystal interact via anisometric Lennard-Jones potentials and can be anchored planar or homeotropically at the substrates. Fluid-substrate attraction is modeled by a Yukawa potential of variable range. By monitoring the grand-potential density and the nematic order parameter as functions of the chemical potential μ, several discontinuous prewetting, wetting, and isotropic-nematic phase transitions are observed. These transitions depend on both the range of the fluid-substrate attraction and the specific anchoring at the substrate. Our results show that at substrates characterized by degenerate anchoring prewetting occurs at lower μ compared with cases in which the anchoring is monostable. This indicates that prewetting transitions are driven by orientational entropy because degenerate anchoring allows for more orientationally distinct configurations of molecules compared with monostable anchoring. In addition, by analyzing local density and various local order parameters, a detailed picture of the structure of various phases emerges from our simulations.
Collapse
Affiliation(s)
- Manuel Greschek
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany.
| | | |
Collapse
|
28
|
Melle M, Giura S, Schlotthauer S, Schoen M. Finite-size scaling analysis of isotropic-polar phase transitions in an amphiphilic fluid. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:035103. [PMID: 22173064 DOI: 10.1088/0953-8984/24/3/035103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
We present Monte Carlo simulations of the isotropic-polar (IP) phase transition in an amphiphilic fluid carried out in the isothermal-isobaric ensemble. Our model consists of Lennard-Jones spheres where the attractive part of the potential is modified by an orientation-dependent function. This function gives rise to an angle dependence of the intermolecular attractions corresponding to that characteristic of point dipoles. Our data show a substantial system-size dependence of the dipolar order parameter. We analyze the system-size dependence in terms of the order-parameter distribution and a cumulant involving its first and second moments. The order parameter, its distribution, and susceptibility observe the scaling behavior characteristic of the 3D Ising universality class. Because of this scaling behavior and because all cumulants have a common intersection irrespective of system size we conclude that the IP phase transition is continuous. Considering pressures 1.3 ≤ P ≤ 3.0 we demonstrate that a line of continuous phase transitions exists which is analogous to the Curie line in systems exhibiting a ferroelectric transition. Our results are qualitatively consistent with Landau's theory of continuous phase transitions.
Collapse
Affiliation(s)
- Michael Melle
- Stranski-Laboratorium für Physikalische und Theoretische Chemie, Fakultät für Mathematik und Naturwissenschaften, Technische Universität Berlin, Berlin, Germany
| | | | | | | |
Collapse
|