1
|
Prampolini G, Greff da Silveira L, Vilhena JG, Livotto PR. Predicting Spontaneous Orientational Self-Assembly: In Silico Design of Materials with Quantum Mechanically Derived Force Fields. J Phys Chem Lett 2022; 13:243-250. [PMID: 34968058 DOI: 10.1021/acs.jpclett.1c03517] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
De novo design of self-assembled materials hinges upon our ability to relate macroscopic properties to individual building blocks, thus characterizing in such supramolecular architectures a wide range of observables at varied time/length scales. This work demonstrates that quantum mechanical derived force fields (QMD-FFs) do satisfy this requisite and, most importantly, do so in a predictive manner. To this end, a specific FF, built solely based on the knowledge of the target molecular structure, is employed to reproduce the spontaneous transition to an ordered liquid crystal phase. The simulations deliver a multiscale portrait of such self-assembly processes, where conformational changes within the individual building blocks are intertwined with a 200 ns ensemble reorganization. The extensive characterization provided not only is in quantitative agreement with the experiment but also connects the time/length scales at which it was performed. Realizing QMD-FF predictive power and unmatched accuracy stands as an important leap forward for the bottom-up design of advanced supramolecular materials.
Collapse
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91 501-970 Porto Alegre, Brazil
| | - J G Vilhena
- Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91 501-970 Porto Alegre, Brazil
| |
Collapse
|
2
|
Vilhena JG, Greff da Silveira L, Livotto PR, Cacelli I, Prampolini G. Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation. J Chem Theory Comput 2021; 17:4449-4464. [PMID: 34185536 DOI: 10.1021/acs.jctc.1c00213] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.
Collapse
Affiliation(s)
- J G Vilhena
- Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-CNR, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
3
|
Sasaki R, Takahashi Y, Hayashi Y, Kawauchi S. Atomistic Mechanism of Anisotropic Heat Conduction in the Liquid Crystal 4-Heptyl-4'-cyanobiphenyl: All-Atom Molecular Dynamics. J Phys Chem B 2020; 124:881-889. [PMID: 31880459 DOI: 10.1021/acs.jpcb.9b08158] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
All-atom molecular dynamics simulations were performed on 4-heptyl-4'-cyanobiphenyl (7CB) to study the mechanism of heat conduction in this nematic liquid crystal atomistically. To describe 7CB properly, the AMBER-type force field was optimized for the dihedral parameter of biphenyl and the Lennard-Jones parameters. The molecular dynamics calculation using the optimized force field well reproduced the experimental values of the isotropic-nematic phase transition temperature, density, and anisotropy of the thermal conductivity. Furthermore, the contributions of convection, intramolecular interaction, and intermolecular interaction to the thermal conductivity were determined by performing thermal conductivity decomposition analysis. According to the analysis, the contributions of convection, bond stretching, and bond bending interactions were higher in the direction parallel to the nematic director than that perpendicular to the director, which is the origin of the anisotropy in the nematic phase. This result indicates that the anisotropy is caused by well-aligned covalent bonds and high mobility parallel to the director. This quantitative description of the mechanism of heat conduction of 7CB is foreseen to provide new insights toward designing highly thermally conductive liquid-crystalline materials.
Collapse
Affiliation(s)
- Ryoma Sasaki
- Department of Chemical Science and Engineering, School of Materials and Chemical Technology , Tokyo Institute of Technology , 2-12-1-E4-6 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan.,Tokyo Tech Academy for Convergence of Materials and Informatics (TAC-MI) , Tokyo Institute of Technology , 2-12-1 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan
| | - Yuki Takahashi
- Department of Chemical Science and Engineering, School of Materials and Chemical Technology , Tokyo Institute of Technology , 2-12-1-E4-6 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan
| | - Yoshihiro Hayashi
- Department of Chemical Science and Engineering, School of Materials and Chemical Technology , Tokyo Institute of Technology , 2-12-1-E4-6 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan.,Research Institute of Polymer Science and Technology (RIPST) , Tokyo Institute of Technology , 2-12-1 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan
| | - Susumu Kawauchi
- Department of Chemical Science and Engineering, School of Materials and Chemical Technology , Tokyo Institute of Technology , 2-12-1-E4-6 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan.,Tokyo Tech Academy for Convergence of Materials and Informatics (TAC-MI) , Tokyo Institute of Technology , 2-12-1 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan.,Research Institute of Polymer Science and Technology (RIPST) , Tokyo Institute of Technology , 2-12-1 Ookayama , Meguro-ku, Tokyo 152-8552 , Japan
| |
Collapse
|
4
|
Affiliation(s)
- Michael P. Allen
- Department of Physics, University of Warwick, Coventry, UK
- H. H. Wills Physics Laboratory, Royal Fort, Bristol, UK
| |
Collapse
|
5
|
Jana PK, Chen W, Alava MJ, Laurson L. Nanoscale liquid crystal lubrication controlled by surface structure and film composition. Phys Chem Chem Phys 2018; 20:18737-18743. [PMID: 29961781 DOI: 10.1039/c8cp03353f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Liquid crystals have emerged as potential candidates for next-generation lubricants due to their tendency to exhibit long-range ordering. Here, we construct a full atomistic model of 4-cyano-4-hexylbiphenyl (6CB) nematic liquid crystal lubricants mixed with hexane and confined by mica surfaces. We explore the effect of the surface structure of mica, as well as lubricant composition and thickness, on the nanoscale friction in the system. Our results demonstrate the key role of the structure of the mica surfaces, specifically the positions of potassium (K+) ions, in determining the nature of sliding friction with monolayer lubricants, including the presence or absence of stick-slip dynamics. With the commensurate setup of confining surfaces, when the grooves created between the periodic K+ ions are parallel to the sliding direction we observe a lower friction force as compared to the perpendicular situation. Random positions of ions exhibit even smaller friction forces with respect to the previous two cases. For thicker lubrication layers the surface structure becomes less important and we observe a good agreement with the experimental data on bulk viscosity of 6CB and the additive hexane. In case of thicker lubrication layers, friction may still be controlled by tuning the relative concentrations of 6CB and hexane in the mixture.
Collapse
Affiliation(s)
- Pritam Kumar Jana
- COMP Centre of Excellence, Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Espoo, Finland.
| | | | | | | |
Collapse
|
6
|
Frechette L, Stratt RM. The inherent dynamics of isotropic- and nematic-phase liquid crystals. J Chem Phys 2017; 144:234505. [PMID: 27334177 DOI: 10.1063/1.4953618] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The geodesic (shortest) pathways through the potential energy landscape of a liquid can be thought of as defining what its dynamics would be if thermal noise were removed, revealing what we have called the "inherent dynamics" of the liquid. We show how these inherent paths can be located for a model liquid crystal former, showing, in the process, how the molecular mechanisms of translation and reorientation compare in the isotropic and nematic phases of these systems. These mechanisms turn out to favor the preservation of local orientational order even under macroscopically isotropic conditions (a finding consistent with the experimental observation of pseudonematic domains in these cases), but disfavor the maintenance of macroscopic orientational order, even in the nematic phase. While the most efficient nematic pathways that maintain nematic order are indeed shorter than those that do not, it is apparently difficult for the system to locate these paths, suggesting that molecular motion in liquid-crystal formers is dynamically frustrated, and reinforcing the sense that there are strong analogies between liquid crystals and supercooled liquids.
Collapse
Affiliation(s)
- Layne Frechette
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Richard M Stratt
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
7
|
Nozawa T, Takahashi KZ, Narumi T, Yasuoka K. Comparison of the accuracy of periodic reaction field methods in molecular dynamics simulations of a model liquid crystal system. J Comput Chem 2015; 36:2406-11. [PMID: 26525311 DOI: 10.1002/jcc.24222] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2015] [Revised: 07/22/2015] [Accepted: 09/07/2015] [Indexed: 11/09/2022]
Abstract
A periodic reaction field (PRF) method is a technique to estimate long-range interactions. The method has the potential to effectively reduce the computational cost while maintaining adequate accuracy. We performed molecular dynamics (MD) simulations of a model liquid-crystal system to assess the accuracy of some variations of the PRF method in low-charge-density systems. All the methods had adequate accuracy compared with the results of the particle mesh Ewald (PME) method, except for a few simulation conditions. Furthermore, in all of the simulation conditions, one of the PRF methods had the same accuracy as the PME method.
Collapse
Affiliation(s)
- Takuma Nozawa
- Department of Mechanical Engineering, Keio University, Yokohama, 223-8522, Japan
| | - Kazuaki Z Takahashi
- Department of Mechanical Engineering, Keio University, Yokohama, 223-8522, Japan.,Research Institute for Sustainable Chemistry, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, 305-8568, Japan
| | - Tetsu Narumi
- Department of Computer Science, University of Electro-Communications, Tokyo, 182-8585, Japan
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, Yokohama, 223-8522, Japan
| |
Collapse
|
8
|
Nozawa T, Takahashi KZ, Kameoka S, Narumi T, Yasuoka K. Application of isotropic periodic sum method for 4-pentyl-4′-cyanobiphenyl liquid crystal. MOLECULAR SIMULATION 2015. [DOI: 10.1080/08927022.2014.998210] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
9
|
Zhang J, Guo H. Transferability of Coarse-Grained Force Field for nCB Liquid Crystal Systems. J Phys Chem B 2014; 118:4647-60. [DOI: 10.1021/jp411615f] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Jianguo Zhang
- Beijing National Laboratory
for Molecular Sciences, State Key Laboratory of Polymer Physics and
Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Hongxia Guo
- Beijing National Laboratory
for Molecular Sciences, State Key Laboratory of Polymer Physics and
Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
10
|
Barone V, Cacelli I, De Mitri N, Licari D, Monti S, Prampolini G. Joyce and Ulysses: integrated and user-friendly tools for the parameterization of intramolecular force fields from quantum mechanical data. Phys Chem Chem Phys 2013; 15:3736-51. [PMID: 23389748 DOI: 10.1039/c3cp44179b] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Joyce program is augmented with several new features, including the user friendly Ulysses GUI, the possibility of complete excited state parameterization and a more flexible treatment of the force field electrostatic terms. A first validation is achieved by successfully comparing results obtained with Joyce2.0 to literature ones, obtained for the same set of benchmark molecules. The parameterization protocol is also applied to two other larger molecules, namely nicotine and a coumarin based dye. In the former case, the parameterized force field is employed in molecular dynamics simulations of solvated nicotine, and the solute conformational distribution at room temperature is discussed. Force fields parameterized with Joyce2.0, for both the dye's ground and first excited electronic states, are validated through the calculation of absorption and emission vertical energies with molecular mechanics optimized structures. Finally, the newly implemented procedure to handle polarizable force fields is discussed and applied to the pyrimidine molecule as a test case.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | | | | | | | | | | |
Collapse
|
11
|
Prampolini G, Bellina F, Biczysko M, Cappelli C, Carta L, Lessi M, Pucci A, Ruggeri G, Barone V. Computational Design, Synthesis, and Mechanochromic Properties of New Thiophene-Based π-Conjugated Chromophores. Chemistry 2012; 19:1996-2004. [DOI: 10.1002/chem.201203672] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2012] [Indexed: 12/25/2022]
|
12
|
Cacelli I, Cimoli A, Livotto PR, Prampolini G. An automated approach for the parameterization of accurate intermolecular force-fields: pyridine as a case study. J Comput Chem 2012; 33:1055-67. [PMID: 22410966 DOI: 10.1002/jcc.22937] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2011] [Revised: 12/23/2011] [Accepted: 12/27/2011] [Indexed: 11/09/2022]
Abstract
An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts.
Collapse
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università degli Studi di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | | | | | | |
Collapse
|
13
|
Ilk Capar M, Cebe E, Zakharov A. Dynamic phenomena and viscous properties in a liquid crystal: A theoretical treatment and molecular dynamic simulations. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.08.037] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
14
|
Zhang J, Su J, Guo H. An Atomistic Simulation for 4-Cyano-4′-pentylbiphenyl and Its Homologue with a Reoptimized Force Field. J Phys Chem B 2011; 115:2214-27. [DOI: 10.1021/jp111408n] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Jianguo Zhang
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Jiaye Su
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Hongxia Guo
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
15
|
Pizzirusso A, Savini M, Muccioli L, Zannoni C. An atomistic simulation of the liquid-crystalline phases of sexithiophene. ACTA ACUST UNITED AC 2011. [DOI: 10.1039/c0jm01284j] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Xue XG, Zhao L, Lu ZY, Li MH, Li ZS. Molecular dynamics simulation study on the isomerization and molecular orientation of liquid crystals formed by azobenzene and (1-cyclohexenyl)phenyldiazene. Phys Chem Chem Phys 2011; 13:11951-7. [DOI: 10.1039/c0cp02094j] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Sarman S, Laaksonen A. Flow alignment phenomena in liquid crystals studied by molecular dynamics simulation. J Chem Phys 2009; 131:144904. [DOI: 10.1063/1.3238549] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
18
|
Cacelli I, Cimoli A, De Gaetani L, Prampolini G, Tani A. Chemical Detail Force Fields for Mesogenic Molecules. J Chem Theory Comput 2009; 5:1865-76. [DOI: 10.1021/ct900002p] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Antonella Cimoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Luca De Gaetani
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Giacomo Prampolini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Alessandro Tani
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| |
Collapse
|