1
|
Gvozden K, Novak Ratajczak S, Orellana AG, Kentzinger E, Rücker U, Dhont JKG, De Michele C, Stiakakis E. Self-Assembly of All-DNA Rods with Controlled Patchiness. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2104510. [PMID: 34837474 DOI: 10.1002/smll.202104510] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Revised: 09/15/2021] [Indexed: 05/23/2023]
Abstract
Double-stranded DNA (dsDNA) fragments exhibit noncovalent attractive interactions between their tips. It is still unclear how DNA liquid crystal self-assembly is affected by such blunt-end attractions. It is demonstrated that stiff dsDNA fragments with moderate aspect ratio can specifically self-assemble in concentrated aqueous solutions into different types of smectic mesophases on the basis of selectively screening of blunt-end DNA stacking interactions. To this end, this type of attractions are engineered at the molecular level by constructing DNA duplexes where the attractions between one or both ends are screened by short hairpin caps. All-DNA bilayer and monolayer smectic-A type of phases, as well as a columnar phase, can be stabilized by controlling attractions strength. The results imply that the so far elusive smectic-A in DNA rod-like liquid crystals is a thermodynamically stable phase. The existence of the bilayer smectic phase is confirmed by Monte-Carlo simulations of hard cylinders decorated with one attractive terminal site. This work demonstrates that DNA blunt-ends behave as well-defined monovalent attractive patches whose strength and position can be potentially precisely tuned, highlighting unique opportunities concerning the stabilization of nonconventional DNA-based lyotropic liquid crystal phases assembled by all-DNA patchy particles with arbitrary geometry and composition.
Collapse
Affiliation(s)
- Katarina Gvozden
- Biomacromolecular Systems and Processes, Institute of Biological Information Processing (IBI-4), Forschungszentrum Jülich, D-52425, Jülich, Germany
| | - Sanja Novak Ratajczak
- Biomacromolecular Systems and Processes, Institute of Biological Information Processing (IBI-4), Forschungszentrum Jülich, D-52425, Jülich, Germany
| | - Alberto G Orellana
- Dipartimento di Fisica, Sapienza Universita di Roma, Piazzale A. Moro 5, Roma, 00185, Italy
| | - Emmanuel Kentzinger
- Jülich Centre for Neutron Science JCNS and Peter Grünberg Institut PGI, JARA-FIT, Forschungszentrum Jülich, D-52425, Jülich, Germany
| | - Ulrich Rücker
- Jülich Centre for Neutron Science JCNS and Peter Grünberg Institut PGI, JARA-FIT, Forschungszentrum Jülich, D-52425, Jülich, Germany
| | - Jan K G Dhont
- Biomacromolecular Systems and Processes, Institute of Biological Information Processing (IBI-4), Forschungszentrum Jülich, D-52425, Jülich, Germany
| | - Cristiano De Michele
- Dipartimento di Fisica, Sapienza Universita di Roma, Piazzale A. Moro 5, Roma, 00185, Italy
| | - Emmanuel Stiakakis
- Biomacromolecular Systems and Processes, Institute of Biological Information Processing (IBI-4), Forschungszentrum Jülich, D-52425, Jülich, Germany
| |
Collapse
|
2
|
Nikoubashman A. Ordering, phase behavior, and correlations of semiflexible polymers in confinement. J Chem Phys 2021; 154:090901. [DOI: 10.1063/5.0038052] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Arash Nikoubashman
- Institute of Physics, Johannes Gutenberg University Mainz, Staudingerweg 7, 55128 Mainz, Germany
| |
Collapse
|
3
|
Peters VFD, Vis M, Wensink HH, Tuinier R. Algebraic equations of state for the liquid crystalline phase behavior of hard rods. Phys Rev E 2020; 101:062707. [PMID: 32688562 DOI: 10.1103/physreve.101.062707] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2020] [Accepted: 06/01/2020] [Indexed: 05/24/2023]
Abstract
Based on simplifications of previous numerical calculations [H. Graf and H. Löwen, Phys. Rev. E 59, 1932 (1999)1063-651X10.1103/PhysRevE.59.1932], we propose algebraic free energy expressions for the smectic-A liquid crystal phase and the crystal phases of hard spherocylinders. Quantitative agreement with simulations is found for the resulting equations of state. The free energy expressions can be used to straightforwardly compute the full phase behavior for all aspect ratios and to provide a suitable benchmark for exploring how attractive interrod interactions mediate the phase stability through perturbation approaches such as free-volume or van der Waals theory.
Collapse
Affiliation(s)
- V F D Peters
- Laboratory of Physical Chemistry, Department of Chemical Engineering and Chemistry & Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
| | - M Vis
- Laboratory of Physical Chemistry, Department of Chemical Engineering and Chemistry & Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
- Laboratoire de Chimie, École Normale Supérieure de Lyon, 69364 Lyon Cedex 07, France
| | - H H Wensink
- Laboratoire de Physique des Solides - UMR 8502, CNRS & Université Paris-Saclay, Orsay, France
| | - R Tuinier
- Laboratory of Physical Chemistry, Department of Chemical Engineering and Chemistry & Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
- Van t Hoff Laboratory for Physical and Colloid Chemistry, Department of Chemistry & Debye Institute for Nanomaterials Science, Utrecht University, 3584 CH Utrecht, The Netherlands
| |
Collapse
|
4
|
Wensink HH. Polymeric Nematics of Associating Rods: Phase Behavior, Chiral Propagation, and Elasticity. Macromolecules 2019. [DOI: 10.1021/acs.macromol.9b01421] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Henricus H. Wensink
- Laboratoire de Physique des Solides—UMR 8502, CNRS, Université Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
5
|
Abstract
Many nanoparticle-based chiral liquid crystals are composed of polydisperse rod-shaped particles with considerable spread in size or shape, affecting the mesoscale chiral properties in, as yet, unknown ways. Using an algebraic interpretation of Onsager-Straley theory for twisted nematics, we investigate the role of length polydispersity on the pitch of nanorod-based cholesterics with a continuous length polydispersity, and find that polydispersity enhances the twist elastic modulus, K 2 , of the cholesteric material without affecting the effective helical amplitude, K t . In addition, for the infinitely large average aspect ratios considered here, the dependence of the pitch on the overall rod concentration is completely unaffected by polydispersity. For a given concentration, the increase in twist elastic modulus (and reduction of the helical twist) may be up to 50% for strong size polydispersity, irrespective of the shape of the unimodal length distribution. We also demonstrate that the twist reduction is reinforced in bimodal distributions, obtained by doping a polydisperse cholesteric with very long rods. Finally, we identify a subtle, non-monotonic change of the pitch across the isotropic-cholesteric biphasic region.
Collapse
|
6
|
Avendaño C, Jackson G, Wensink HH. Nanorings in planar confinement: the role of repulsive surfaces on the formation of lacuna smectics. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1484950] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Carlos Avendaño
- School of Chemical Engineering and Analytical Science, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - George Jackson
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, UK
| | - Henricus H. Wensink
- Laboratoire de Physique des Solides UMR 8502, CNRS, Universite Paris-Sud, Universite Paris-Saclay, Orsay, France
| |
Collapse
|
7
|
Hvozd M, Patsahan T, Holovko M. Isotropic-Nematic Transition and Demixing Behavior in Binary Mixtures of Hard Spheres and Hard Spherocylinders Confined in a Disordered Porous Medium: Scaled Particle Theory. J Phys Chem B 2018; 122:5534-5546. [PMID: 29385800 DOI: 10.1021/acs.jpcb.7b11834] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We develop the scaled particle theory to describe the thermodynamic properties and orientation ordering of a binary mixture of hard spheres (HS) and hard spherocylinders (HSC) confined in a disordered porous medium. Using this theory, the analytical expressions of the free energy, the pressure, and the chemical potentials of HS and HSC have been derived. The improvement of obtained results is considered by introducing the Carnahan-Starling-like and Parsons-Lee-like corrections. Phase diagrams for the isotropic-nematic transition are calculated from the bifurcation analysis of the integral equation for the orientation singlet distribution function and from the conditions of thermodynamic equilibrium. Both the approaches correctly predict the isotropic-nematic transition at low concentrations of hard spheres. However, the thermodynamic approach provides more accurate results and is able to describe the demixing phenomena in the isotropic and nematic phases. The effects of porous medium on the isotropic-nematic phase transition and demixing behavior in a binary HS/HSC mixture are discussed.
Collapse
Affiliation(s)
- M Hvozd
- Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine , 1 Svientsitskii Street , 79011 Lviv , Ukraine
| | - T Patsahan
- Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine , 1 Svientsitskii Street , 79011 Lviv , Ukraine
| | - M Holovko
- Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine , 1 Svientsitskii Street , 79011 Lviv , Ukraine
| |
Collapse
|
8
|
Langenbach K. Co-Oriented Fluid Functional Equation for Electrostatic interactions (COFFEE). Chem Eng Sci 2017. [DOI: 10.1016/j.ces.2017.08.025] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
9
|
Brumby PE, Wensink HH, Haslam AJ, Jackson G. Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11754-11770. [PMID: 28885848 DOI: 10.1021/acs.langmuir.7b02254] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.
Collapse
Affiliation(s)
- Paul E Brumby
- Department of Mechanical Engineering, Keio University , 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Henricus H Wensink
- Laboratoire de Physique des Solides, Université Paris Sud & CNRS , 91405 Orsay Cedex, France
| | - Andrew J Haslam
- Department of Chemical Engineering and Qatar Carbonates and Carbon Storage Research Centre, Imperial College London , South Kensington Campus, London SW7 2AZ, United Kingdom
| | - George Jackson
- Department of Chemical Engineering and Qatar Carbonates and Carbon Storage Research Centre, Imperial College London , South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
10
|
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
|
11
|
Rickayzen G, Heyes DM. Isotropic-nematic phase transition of uniaxial variable softness prolate and oblate ellipsoids. J Chem Phys 2017; 146:164505. [PMID: 28456193 DOI: 10.1063/1.4981887] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Onsager's theory of the isotropic-nematic phase separation of rod shaped particles is generalized to include particle softness and attractions in the anisotropic interparticle force field. The procedure separates a scaled radial component from the angular integral part, the latter being treated in essentially the same way as in the original Onsager formulation. Building on previous treatments of more idealised hard-core particle models, this is a step toward representing more realistic rod-like systems and also allowing temperature (and in principle specific chemical factors) to be included at a coarse grained level in the theory. The focus of the study is on the coexisting concentrations and associated coexistence properties. Prolate and oblate ellipsoids are considered in both the small and very large aspect ratio limits. Approximations to the terms in the angular integrals derived assuming the very large (prolate) and very small (oblate) aspect ratios limits are compared with the formally exact treatment. The approximation for the second virial coefficient matches the exact solution for aspect ratios above about 20 for the prolate ellipsoids and less than ca. 0.05 for the oblate ellipsoids from the numerical evaluation of the angular integrals. The temperature dependence of the coexistence density could be used to help determine the interaction potential of two molecules. The method works at temperatures above a certain threshold temperature where the second virial coefficient is positive.
Collapse
Affiliation(s)
- G Rickayzen
- School of Physical Sciences, University of Kent, Canterbury, Kent CT2 7NH, United Kingdom
| | - D M Heyes
- Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 OEX, United Kingdom
| |
Collapse
|
12
|
Ferreiro-Córdova C, Wensink HH. Spinodal instabilities in polydisperse lyotropic nematics. J Chem Phys 2016; 145:244904. [DOI: 10.1063/1.4972523] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Affiliation(s)
- C. Ferreiro-Córdova
- Laboratoire de Physique des Solides - UMR 8502, CNRS, Université Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| | - H. H. Wensink
- Laboratoire de Physique des Solides - UMR 8502, CNRS, Université Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
13
|
Meyer H, van der Schoot P, Schilling T. Percolation in suspensions of polydisperse hard rods: Quasi universality and finite-size effects. J Chem Phys 2015; 143:044901. [PMID: 26233158 DOI: 10.1063/1.4926946] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a study of connectivity percolation in suspensions of hard spherocylinders by means of Monte Carlo simulation and connectedness percolation theory. We focus attention on polydispersity in the length, the diameter, and the connectedness criterion, and we invoke bimodal, Gaussian, and Weibull distributions for these. The main finding from our simulations is that the percolation threshold shows quasi universal behaviour, i.e., to a good approximation, it depends only on certain cumulants of the full size and connectivity distribution. Our connectedness percolation theory hinges on a Lee-Parsons type of closure recently put forward that improves upon the often-used second virial approximation [T. Schilling, M. Miller, and P. van der Schoot, e-print arXiv:1505.07660 (2015)]. The theory predicts exact universality. Theory and simulation agree quantitatively for aspect ratios in excess of 20, if we include the connectivity range in our definition of the aspect ratio of the particles. We further discuss the mechanism of cluster growth that, remarkably, differs between systems that are polydisperse in length and in width, and exhibits non-universal aspects.
Collapse
Affiliation(s)
- Hugues Meyer
- Research Unit for Physics and Materials Science, Université du Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Paul van der Schoot
- Department of Applied Physics, Eindhoven University of Technology, P.O. Box 513, 3500 MB Eindhoven, The Netherlands
| | - Tanja Schilling
- Research Unit for Physics and Materials Science, Université du Luxembourg, L-1511 Luxembourg, Luxembourg
| |
Collapse
|
14
|
Chatterjee AP. Connectedness percolation in isotropic systems of monodisperse spherocylinders. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:375302. [PMID: 26325334 DOI: 10.1088/0953-8984/27/37/375302] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The dependence upon aspect ratio and connectedness range of the percolation threshold for monodisperse spherocylinders as calculated from an analogy between lattice and continuum percolation is compared with findings from: (i) recent Monte Carlo (MC) simulations, and: (ii) a recent elaboration of an integral equation-based approach that uses the connectedness Ornstein-Zernike equation. The results from each of these approaches are found to be in reasonable semi-quantitative agreement, especially for the case of particles that have a finite hard core diameter.
Collapse
Affiliation(s)
- Avik P Chatterjee
- Department of Chemistry, SUNY College of Environmental Science and Forestry, One Forestry Drive, Syracuse, NY 13210, USA
| |
Collapse
|
15
|
Chatterjee AP, Grimaldi C. Random geometric graph description of connectedness percolation in rod systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:032121. [PMID: 26465440 DOI: 10.1103/physreve.92.032121] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/22/2015] [Indexed: 06/05/2023]
Abstract
The problem of continuum percolation in dispersions of rods is reformulated in terms of weighted random geometric graphs. Nodes (or sites or vertices) in the graph represent spatial locations occupied by the centers of the rods. The probability that an edge (or link) connects any randomly selected pair of nodes depends upon the rod volume fraction as well as the distribution over their sizes and shapes, and also upon quantities that characterize their state of dispersion (such as the orientational distribution function). We employ the observation that contributions from closed loops of connected rods are negligible in the limit of large aspect ratios to obtain percolation thresholds that are fully equivalent to those calculated within the second-virial approximation of the connectedness Ornstein-Zernike equation. Our formulation can account for effects due to interactions between the rods, and many-body features can be partially addressed by suitable choices for the edge probabilities.
Collapse
Affiliation(s)
- Avik P Chatterjee
- Department of Chemistry, SUNY College of Environmental Science and Forestry, One Forestry Drive, Syracuse, New York 13210, USA
| | - Claudio Grimaldi
- Laboratory of Physics of Complex Matter, Ecole Polytechnique Fédérale de Lausanne, Station 3, CH-1015 Lausanne, Switzerland
| |
Collapse
|
16
|
Wu L, Malijevský A, Jackson G, Müller EA, Avendaño C. Orientational ordering and phase behaviour of binary mixtures of hard spheres and hard spherocylinders. J Chem Phys 2015; 143:044906. [DOI: 10.1063/1.4923291] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Liang Wu
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Alexandr Malijevský
- Department of Physical Chemistry, ICT Prague, 166 28 Praha 6, Czech Republic, and Institute of Chemical Process Fundamentals of ASCR, 16502 Praha 6, Czech Republic
| | - George Jackson
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Carlos Avendaño
- School of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, United Kingdom
| |
Collapse
|
17
|
van Westen T, Oyarzún B, Vlugt TJH, Gross J. An analytical equation of state for describing isotropic-nematic phase equilibria of Lennard-Jones chain fluids with variable degree of molecular flexibility. J Chem Phys 2015; 142:244903. [DOI: 10.1063/1.4922921] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
|
18
|
Oyarzún B, van Westen T, Vlugt TJH. Isotropic-nematic phase equilibria of hard-sphere chain fluids—Pure components and binary mixtures. J Chem Phys 2015; 142:064903. [DOI: 10.1063/1.4907639] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
19
|
Cuetos A, Martínez-Haya B. Liquid crystal phase diagram of soft repulsive rods and its mapping on the hard repulsive reference fluid. Mol Phys 2015. [DOI: 10.1080/00268976.2014.996191] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
20
|
Wensink HH, Trizac E. Generalized Onsager theory for strongly anisometric patchy colloids. J Chem Phys 2014; 140:024901. [PMID: 24437905 DOI: 10.1063/1.4851217] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The implications of soft "patchy" interactions on the orientational disorder-order transition of strongly elongated colloidal rods and flat disks is studied within a simple Onsager-van der Waals density functional theory. The theory provides a generic framework for studying the liquid crystal phase behaviour of highly anisometric cylindrical colloids which carry a distinct geometrical pattern of repulsive or attractive soft interactions localized on the particle surface. In this paper, we apply our theory to the case of charged rods and disks for which the local electrostatic interactions can be described by a screened-Coulomb potential. We consider infinitely thin rod like cylinders with a uniform line charge and infinitely thin discotic cylinders with several distinctly different surface charge patterns. Irrespective of the backbone shape the isotropic-nematic phase diagrams of charged colloids feature a generic destabilization of nematic order at low ionic strength, a dramatic narrowing of the biphasic density region, and a reentrant phenomenon upon reducing the electrostatic screening. The low screening regime is characterized by a complete suppression of nematic order in favor of positionally ordered liquid crystal phases.
Collapse
Affiliation(s)
- H H Wensink
- Laboratoire de Physique des Solides - UMR 8502, Université Paris-Sud and CNRS, 91405 Orsay Cedex, France
| | - E Trizac
- Laboratoire de Physique Théorique et Modèles Statistiques - UMR 8626, Université Paris-Sud and CNRS, 91405 Orsay Cedex, France
| |
Collapse
|
21
|
Chatterjee AP. Percolation thresholds for polydisperse circular disks: A lattice-based exploration. J Chem Phys 2014; 141:034903. [DOI: 10.1063/1.4890280] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Wu L, Müller EA, Jackson G. Understanding and Describing the Liquid-Crystalline States of Polypeptide Solutions: A Coarse-Grained Model of PBLG in DMF. Macromolecules 2014. [DOI: 10.1021/ma401230x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Liang Wu
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Erich A. Müller
- 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
|
23
|
van Westen T, Vlugt TJH, Gross J. The isotropic-nematic and nematic-nematic phase transition of binary mixtures of tangent hard-sphere chain fluids: An analytical equation of state. J Chem Phys 2014; 140:034504. [DOI: 10.1063/1.4860980] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
24
|
Wu L, Jackson G, Müller EA. Liquid crystal phase behaviour of attractive disc-like particles. Int J Mol Sci 2013; 14:16414-42. [PMID: 23965962 PMCID: PMC3759919 DOI: 10.3390/ijms140816414] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2013] [Revised: 07/24/2013] [Accepted: 07/25/2013] [Indexed: 11/17/2022] Open
Abstract
We employ a generalized van der Waals-Onsager perturbation theory to construct a free energy functional capable of describing the thermodynamic properties and orientational order of the isotropic and nematic phases of attractive disc particles. The model mesogen is a hard (purely repulsive) cylindrical disc particle decorated with an anisotropic square-well attractive potential placed at the centre of mass. Even for isotropic attractive interactions, the resulting overall inter-particle potential is anisotropic, due to the orientation-dependent excluded volume of the underlying hard core. An algebraic equation of state for attractive disc particles is developed by adopting the Onsager trial function to characterize the orientational order in the nematic phase. The theory is then used to represent the fluid-phase behaviour (vapour-liquid, isotropic-nematic, and nematic-nematic) of the oblate attractive particles for varying values of the molecular aspect ratio and parameters of the attractive potential. When compared to the phase diagram of their athermal analogues, it is seen that the addition of an attractive interaction facilitates the formation of orientationally-ordered phases. Most interestingly, for certain aspect ratios, a coexistence between two anisotropic nematic phases is exhibited by the attractive disc-like fluids.
Collapse
Affiliation(s)
- Liang Wu
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, UK; E-Mails: (L.W.); (G.J.)
| | - George Jackson
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, UK; E-Mails: (L.W.); (G.J.)
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, UK; E-Mails: (L.W.); (G.J.)
| |
Collapse
|
25
|
van Westen T, Oyarzún B, Vlugt TJH, Gross J. The isotropic-nematic phase transition of tangent hard-sphere chain fluids--pure components. J Chem Phys 2013; 139:034505. [PMID: 23883045 DOI: 10.1063/1.4815916] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An extension of Onsager's second virial theory is developed to describe the isotropic-nematic phase transition of tangent hard-sphere chain fluids. Flexibility is introduced by the rod-coil model. The effect of chain-flexibility on the second virial coefficient is described using an accurate, analytical approximation for the orientation-dependent pair-excluded volume. The use of this approximation allows for an analytical treatment of intramolecular flexibility by using a single pure-component parameter. Two approaches to approximate the effect of the higher virial coefficients are considered, i.e., the Vega-Lago rescaling and Scaled Particle Theory (SPT). The Onsager trial function is employed to describe the orientational distribution function. Theoretical predictions for the equation of state and orientational order parameter are tested against the results from Monte Carlo (MC) simulations. For linear chains of length 9 and longer, theoretical results are in excellent agreement with MC data. For smaller chain lengths, small errors introduced by the approximation of the higher virial coefficients become apparent, leading to a small under- and overestimation of the pressure and density difference at the phase transition, respectively. For rod-coil fluids of reasonable rigidity, a quantitative comparison between theory and MC simulations is obtained. For more flexible chains, however, both the Vega-Lago rescaling and SPT lead to a small underestimation of the location of the phase transition.
Collapse
Affiliation(s)
- Thijs van Westen
- Process and Energy Laboratory, Delft University of Technology, Leeghwaterstraat 44, 2628 CA Delft, The Netherlands
| | | | | | | |
Collapse
|
26
|
Franco-Melgar M, Haslam AJ, Jackson G. A generalisation of the Onsager trial function approach: Describing nematic liquid crystals with an algebraic equation of state. Mol Phys 2012. [DOI: 10.1080/00268976.2012.750060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
27
|
van Westen T, Vlugt TJH, Gross J. An analytical approximation for the orientation-dependent excluded volume of tangent hard sphere chains of arbitrary chain length and flexibility. J Chem Phys 2012; 137:044906. [DOI: 10.1063/1.4737663] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Wu L, Wensink H, Jackson G, Müller E. A generic equation of state for liquid crystalline phases of hard-oblate particles. Mol Phys 2012. [DOI: 10.1080/00268976.2011.649794] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
29
|
Ilg P. Enhanced Landau-de Gennes potential for nematic liquid crystals from a systematic coarse-graining procedure. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:061709. [PMID: 23005116 DOI: 10.1103/physreve.85.061709] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2011] [Revised: 04/27/2012] [Indexed: 06/01/2023]
Abstract
The macroscopic theory of nematics is conveniently described in terms of the phenomenological Landau-de Gennes free energy. Here we show how such an effective free energy can be obtained explicitly from a microscopic model via the help of a systematic coarse-graining procedure. We test our approach for the two- and three-dimensional Lebwohl-Lasher model of nematics. The effective free energy that we obtain is consistent with the phenomenological Landau-de Gennes form for weak orientational ordering and the Maier-Saupe theory of the isotropic-nematic transition. For strong orientational ordering, however, the effective free energy increases rapidly and diverges logarithmically near the fully oriented state. The explicit form for the regularized Landau-de Gennes potential proposed here restricts the order parameter to physical admissible values and reproduces our numerical data accurately.
Collapse
Affiliation(s)
- Patrick Ilg
- ETH Zürich, Department of Materials, Polymer Physics, HCI H541, CH-8093 Zürich, Switzerland
| |
Collapse
|
30
|
You XM, Vlasov AY, Anton L, Masters AJ. Isotropic and nematic phases of the hard spheroid fluid: a virial approach. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:061706. [PMID: 23005113 DOI: 10.1103/physreve.85.061706] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/01/2012] [Indexed: 06/01/2023]
Abstract
We use a high level virial expansion to investigate the properties of the isotropic and nematic phases of the hard spheroid fluid. We use the Monte Carlo techniques described previously to calculate the virial coefficients up to seventh order and we represent the dependence of these coefficients on particle orientations via a spherical harmonic expansion. The expansion coefficients are determined using Lebedev quadrature which carries out the angular integration required exactly. For fairly spherical spheroids (1/3<a/b<3, where a and b are the semiaxes of the spheroid) the virial expansion appears to be convergent in both the isotropic and nematic phases up to high density. For more aspherical particles, the predictions varied in a more erratic fashion upon increasing the number of virials. The use of a finite angular basis set to represent the orientationally dependent virial coefficients could be a contributory cause to this.
Collapse
Affiliation(s)
- X-M You
- School of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | | | | | | |
Collapse
|
31
|
Wensink HH, Jackson G. Cholesteric order in systems of helical Yukawa rods. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2011; 23:194107. [PMID: 21525561 DOI: 10.1088/0953-8984/23/19/194107] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We consider the interaction potential between two chiral rod-like colloids which consist of a thin cylindrical backbone decorated with a helical charge distribution on the cylinder surface. For sufficiently slender helical rods a simple scaling expression is derived which relates the chiral 'twisting' potential to the microscopic properties of the particles, such as the internal helical pitch, charge density and electrostatic screening parameter. To predict the behaviour of the macroscopic cholesteric pitch of the fluid bulk phase we invoke a simple second-virial theory generalized to treat anisotropic states with weakly twisted director fields. It is shown that, while particles with weakly coiled helices always form a cholesteric phase whose helical sense is commensurate with that of the internal helix, more strongly coiled rods lead to the formation of a cholesteric state of opposite sense. The correlation between the helical symmetry at the microscopic and macroscopic scale is found to be very sensitive to the pitch of the Yukawa helix. Mixing helical particles of sufficiently disparate length and internal pitch may give rise to a demixing of the uniform cholesteric phase into two fractions with a different macroscopic pitch. Our findings could be relevant to the interpretation of experimental observations in systems of cellulose and chitin microfibres, DNA and fd virus rods.
Collapse
Affiliation(s)
- H H Wensink
- Department of Chemical Engineering, Imperial College London, London, UK.
| | | |
Collapse
|
32
|
Varga S, Jackson G. A study of steric chirality: the chiral nematic phase of a system of chiral two-site HGO molecules. Mol Phys 2011. [DOI: 10.1080/00268976.2011.556577] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
33
|
Zheng YX, Yu YX, Li YF. An Equation of State for the Isotropic–Nematic Phase Transition of Semiflexible Polymers. Ind Eng Chem Res 2011. [DOI: 10.1021/ie102379h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yuan-Xiang Zheng
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Yang-Xin Yu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Ying-Feng Li
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
34
|
Brumby PE, Haslam AJ, de Miguel E, Jackson G. Subtleties in the calculation of the pressure and pressure tensor of anisotropic particles from volume-perturbation methods and the apparent asymmetry of the compressive and expansive contributions. Mol Phys 2011. [DOI: 10.1080/00268976.2010.530301] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
35
|
Kumar R, Audus D, Fredrickson GH. Phase Separation in Symmetric Mixtures of Oppositely Charged Rodlike Polyelectrolytes. J Phys Chem B 2010; 114:9956-76. [DOI: 10.1021/jp101413a] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rajeev Kumar
- Materials Research Laboratory, University of California at Santa Barbara, Santa Barbara, California, 93106-5080
| | - Debra Audus
- Materials Research Laboratory, University of California at Santa Barbara, Santa Barbara, California, 93106-5080
| | - Glenn H. Fredrickson
- Materials Research Laboratory, University of California at Santa Barbara, Santa Barbara, California, 93106-5080
| |
Collapse
|
36
|
Llovell F, Galindo A, Blas FJ, Jackson G. Classical density functional theory for the prediction of the surface tension and interfacial properties of fluids mixtures of chain molecules based on the statistical associating fluid theory for potentials of variable range. J Chem Phys 2010; 133:024704. [DOI: 10.1063/1.3449143] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
|
38
|
|
39
|
Cheung DL, Schmidt M. Quenched-annealed density functional theory for interfacial behavior of hard rods at a hard rod matrix. J Chem Phys 2009; 131:214705. [DOI: 10.1063/1.3267728] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
40
|
Avendaño C, Müller EA. Liquid crystalline and antinematic behavior of shape-persistent macrocycles from molecular-dynamics simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:061702. [PMID: 20365180 DOI: 10.1103/physreve.80.061702] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2009] [Indexed: 05/29/2023]
Abstract
In this work we present a molecular-dynamics study of a coarse-grained (CG) model for a system of planar shape-persistent macrocycles (SPMs). SPMs are synthetic organic rigid macromolecules typically comprised of meta- and para-aromatics groups connected by acetylene and/or diacetylene units. In the CG model, each SPM is represented as a rigid hexagonal arrangement of 24 soft-repulsive spheres, resembling a large ring or hoop. The supramolecular arrangement of these macrocycles at high pressures is studied using N-P-T molecular-dynamics simulation both by expansion of an initial hexagonal lattice structure and also by compression of an isotropic phase. In both cases, systems under consideration exhibit an isotropic-smectic-A phase transition, which is detected by monitoring relevant order parameters and analyzing snapshots of equilibrium configurations. The smectic-A phase is unique; although the molecules form layers, the system presents antinematic order where the orientation of the molecular axes is perpendicular to the direction of the layers themselves. Due to their planar geometry, the SPM molecules would be expected to form columnar or nematic phases. On the contrary, these phases seem suppressed by a novel smectic-A phase, formed by the mutual interpenetration of the cycles. These results are a unique example of how molecular nonconvexity can, by itself, induce mesomorphism in anisotropic systems.
Collapse
Affiliation(s)
- Carlos Avendaño
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | | |
Collapse
|
41
|
Franco-Melgar M, Haslam AJ, Jackson G. Advances in generalised van der Waals approaches for the isotropic–nematic fluid phase equilibria of thermotropic liquid crystals–an algebraic equation of state for attractive anisotropic particles with the Onsager trial function. Mol Phys 2009. [DOI: 10.1080/00268970903352335] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
42
|
Holmqvist P, Ratajczyk M, Meier G, Wensink HH, Lettinga MP. Supersaturated dispersions of rodlike viruses with added attraction. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:031402. [PMID: 19905111 DOI: 10.1103/physreve.80.031402] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/14/2009] [Indexed: 05/20/2023]
Abstract
The kinetics of isotropic-nematic (I-N) and nematic-isotropic (N-I) phase transitions in dispersions of rodlike fd viruses are studied. Concentration quenches were applied using pressure jumps in combination with polarization microscopy, birefringence, and turbidity measurements. The full biphasic region could be accessed, resulting in the construction of an experimental analog of the bifurcation diagram. The N-I spinodal points for dispersions of rods with varying concentrations of depletion agent (dextran) were obtained from orientation quenches using cessation of shear flow in combination with small-angle light scattering. We found that the location of the N-I spinodal point is independent of the attraction, which was confirmed by theory. Surprisingly, the experiments showed that also the absolute induction time, the critical nucleus, and the growth rate are insensitive of the attraction if the concentration is scaled to the distance to the phase boundaries.
Collapse
Affiliation(s)
- P Holmqvist
- Institut für Festkörperforschung, Forschungszentrum Jülich, D-52425 Jülich, Germany
| | | | | | | | | |
Collapse
|
43
|
Wensink HH, Jackson G. Generalized van der Waals theory for the twist elastic modulus and helical pitch of cholesterics. J Chem Phys 2009; 130:234911. [DOI: 10.1063/1.3153348] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Malijevský A, Jackson G, Varga S. Many-fluid Onsager density functional theories for orientational ordering in mixtures of anisotropic hard-body fluids. J Chem Phys 2008; 129:144504. [DOI: 10.1063/1.2982501] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|