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
|
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
|
4
|
Cacelli I, Prampolini G. Parametrization and Validation of Intramolecular Force Fields Derived from DFT Calculations. J Chem Theory Comput 2015; 3:1803-17. [PMID: 26627623 DOI: 10.1021/ct700113h] [Citation(s) in RCA: 87] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The energy and its first and second geometrical derivatives obtained by DFT calculations for a number of conformations of a single molecule are used to parametrize intramolecular force fields, suitable for computer simulations. A systematic procedure is proposed to adequately treat either fully atomistic or more simplified force fields, as within the united atom approach or other coarse grained models. The proposed method is tested and validated by performing molecular dynamics simulations on several different molecules, comparing the results with literature force fields and relevant experimental data. Particular emphasis is given to the united atom approach for flexible molecules characterized by "soft" torsional potentials which are known to retain a high degree of chemical specificity.
Collapse
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Giacomo Prampolini
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| |
Collapse
|
5
|
Sarman S, Wang YL, Laaksonen A. Non-Newtonian rheological properties of shearing nematic liquid crystal model systems based on the Gay–Berne potential. Phys Chem Chem Phys 2015; 17:16615-23. [DOI: 10.1039/c5cp02468d] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The logarithm of the viscosity of a nematic liquid crystal is a linear function of the square root of the shear rate in the non-Newtonian regime.
Collapse
Affiliation(s)
- Sten Sarman
- Department of Materials and Environmental Chemistry
- Arrhenius Laboratory
- Stockholm University
- 106 91 Stockholm
- Sweden
| | - Yong-Lei Wang
- Department of Materials and Environmental Chemistry
- Arrhenius Laboratory
- Stockholm University
- 106 91 Stockholm
- Sweden
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry
- Arrhenius Laboratory
- Stockholm University
- 106 91 Stockholm
- Sweden
| |
Collapse
|
6
|
Prampolini G, Yu P, Pizzanelli S, Cacelli I, Yang F, Zhao J, Wang J. Structure and Dynamics of Ferrocyanide and Ferricyanide Anions in Water and Heavy Water: An Insight by MD Simulations and 2D IR Spectroscopy. J Phys Chem B 2014; 118:14899-912. [DOI: 10.1021/jp511391b] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|
7
|
Pizzirusso A, Di Pietro ME, De Luca G, Celebre G, Longeri M, Muccioli L, Zannoni C. Order and Conformation of Biphenyl in Cyanobiphenyl Liquid Crystals: A Combined Atomistic Molecular Dynamics and1H NMR Study. Chemphyschem 2014; 15:1356-67. [DOI: 10.1002/cphc.201400082] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2014] [Indexed: 11/07/2022]
|
8
|
Sarman S, Laaksonen A. Twist viscosities and flow alignment of biaxial nematic liquid crystal phases of a soft ellipsoid-string fluid studied by molecular dynamics simulation. Phys Chem Chem Phys 2012; 14:11999-2013. [DOI: 10.1039/c2cp41456b] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Affiliation(s)
- Sten Sarman
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, 106 91 Stockholm, Sweden.
| | | |
Collapse
|
9
|
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]
|
10
|
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
|
11
|
Cacelli I, Lami CF, Prampolini G. Force-field modeling through quantum mechanical calculations: Molecular dynamics simulations of a nematogenic molecule in its condensed phases. J Comput Chem 2009; 30:366-78. [DOI: 10.1002/jcc.21062] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
12
|
Tiberio G, Muccioli L, Berardi R, Zannoni C. Towardsin SilicoLiquid Crystals. Realistic Transition Temperatures and Physical Properties forn-Cyanobiphenyls via Molecular Dynamics Simulations. Chemphyschem 2009; 10:125-36. [DOI: 10.1002/cphc.200800231] [Citation(s) in RCA: 147] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
13
|
Cifelli M, De Gaetani L, Prampolini G, Tani A. Atomistic Computer Simulation and Experimental Study on the Dynamics of the n-Cyanobiphenyls Mesogenic Series. J Phys Chem B 2008; 112:9777-86. [DOI: 10.1021/jp802935q] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Mario Cifelli
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Luca De Gaetani
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Giacomo Prampolini
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | - Alessandro Tani
- Dipartimento di Chimica e Chimica Industriale, Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| |
Collapse
|
14
|
De Gaetani L, Prampolini G, Tani A. Subdiffusive dynamics of a liquid crystal in the isotropic phase. J Chem Phys 2008; 128:194501. [DOI: 10.1063/1.2916681] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
|
15
|
Cacelli I, Gaetani LD, Prampolini G, Tani A. Liquid Crystal Properties of the n-Alkyl-cyanobiphenyl Series from Atomistic Simulations with Ab Initio Derived Force Fields. J Phys Chem B 2007; 111:2130-7. [PMID: 17288473 DOI: 10.1021/jp065806l] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Lengthy molecular dynamics (MD) simulations were performed at constant atmospheric pressure and different temperatures for the series of the 4-n-alkyl-4'-cyanobiphenyls (nCB) with n = 6, 7, and 8. The accurate atomistic force field (Bizzarri, M.; Cacelli, I.; Prampolini, G; Tani, A. J. Phys. Chem. A 2004, 108, 10336), successfully employed to reproduce thermodynamic and transport properties of the 5CB molecule, has here been extended to higher homologues. Nematic and isotropic phases were found for all members of the series, and also, a smectic phase was (tentatively) identified for 8CB at 1 atm and 300 K. Transition temperatures reproduce the experimental values within +/-10 K. Also, structural properties as second and fourth rank orientational order parameters are in good agreement with the corresponding experimental quantities. This means that the well-known odd-even effect, observed for many properties along the nCB series, is well reproduced, despite the narrow range of oscillations, e.g., in clearing temperatures. A detailed analysis of the correlation between molecular properties and odd-even effects is presented.
Collapse
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale Universita di Pisa, via Risorgimento 35, I-56126 Pisa, Italy.
| | | | | | | |
Collapse
|
16
|
Peláez J, Wilson M. Molecular orientational and dipolar correlation in the liquid crystal mixture E7: a molecular dynamics simulation study at a fully atomistic level. Phys Chem Chem Phys 2007; 9:2968-75. [PMID: 17551620 DOI: 10.1039/b614422e] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular dynamics simulations are reported for the four component nematic liquid crystal mixture E7, which is used commercially. We are able to show the growth of a nematic phase directly from an isotropic liquid over a 100 ns period for an all-atom model, and study orientational and dipole order within the nematic phase. The simulations show that the cyanoterphenyl component of the mixture, 5CT, is more ordered than the three cyanobiphenyl components. The simulations show also that both parallel and anti-parallel dipole correlation take place in E7 but that the strong anti-parallel dipole correlation is localised to particular arrangements of molecules. It is possible to identify two key preferred configurations for molecular pairs in the fluid, which explain the form of the dipole correlation function, g(1)(r).
Collapse
Affiliation(s)
- Jorge Peláez
- Department of Chemistry, Durham University, South Road, Durham, UKDH1 3LE.
| | | |
Collapse
|
17
|
Amovilli C, Cacelli I, Cinacchi G, De Gaetani L, Prampolini G, Tani A. Structure and dynamics of mesogens using intermolecular potentials derived from ab initio calculations. Theor Chem Acc 2006. [DOI: 10.1007/s00214-006-0209-y] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
18
|
Prampolini G. Parametrization and Validation of Coarse Grained Force-Fields Derived from ab Initio Calculations. J Chem Theory Comput 2006; 2:556-67. [DOI: 10.1021/ct050328o] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Giacomo Prampolini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| |
Collapse
|