1
|
Bhendale M, Srivastava A, Singh JK. Insights into the Phase Diagram of Pluronic L64 Using Coarse-Grained Molecular Dynamics Simulations. J Phys Chem B 2022; 126:4731-4744. [PMID: 35708274 DOI: 10.1021/acs.jpcb.2c02429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate the concentration-dependent phase diagram of pluronic L64 in aqueous media at 300 and 320 K using coarse-grained (CG) molecular dynamics (MD) simulations. The CG model is derived by adapting the Martini model for nonbonded interactions along with the Boltzmann inversion (BI) of bonded interactions from all-atom (AA) simulations. Our derived CG model successfully captures the experimentally observed micellar-, hexagonal-, lamellar-, and polymer-rich solution phase. The end-to-end distance reveals the conformational change from an open-chain structure in the micellar phase to a folded-chain structure in the lamellar phase, increasing the orientational order. An increase in temperature leads to expulsion of water molecules from the L64 moiety, suggesting an increase in L64 hydrophobicity. Thermodynamic analysis using the two-phase thermodynamics (2PT) method suggests the entropy of the system decreases with increasing L64 concentration and the decrease in free energy (F) with temperature is mainly driven by the entropic factor (-TS). Further, the increase in aggregation number at lower concentrations and self-assembly at very high concentrations is energetically driven, whereas the change from the micellar phase to the lamellar phase with increasing L64 concentration is entropically driven. Our model provides molecular insights into L64 phases which can be further explored to design functionality-based suprastructures for drug delivery and tissue engineering applications.
Collapse
Affiliation(s)
- Mangesh Bhendale
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Arpita Srivastava
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India.,Prescience Insilico Private Limited, Fifth Floor, Novel MSR Building, Marathahalli, Bengaluru, Karnataka 560037, India
| |
Collapse
|
2
|
Fuller J, Moon J, Zhang J, Chidambaram D, An Q. Coordination and thermophysical properties of select trivalent lanthanides in LiCl-KCl. Phys Chem Chem Phys 2022; 24:13102-13109. [PMID: 35588243 DOI: 10.1039/d1cp03056f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The coordination chemistry of various fission and decay products, such as actinides and lanthanides, are crucial to the commercial deployment of molten salt reactors as they can affect the thermophysical properties. Here, we examined the structure, coordination environment, and physical properties such as the density and the vibrational density of states for three lanthanide species, namely Ce, Eu, and Sm in the LiCl-KCl eutectic system using a combination of quantum mechanics simulations and spectroscopic experiments. Quantum mechanics molecular dynamics (QM-MD) modelling was employed to determine the physical properties of each system resulting in accurate local coordination of each species. Then, the vibrational density of states (DOS) was determined using a two-phase thermodynamic modelling which was then compared to the experimentally obtained Raman spectra of the species in molten LiCl-KCl having the eutectic composition. We find that Ce3+, Eu3+ and Sm3+ all adopt octahedral local coordination environments in the eutectic salt composition in good agreement with experimental results. Ce3+ is found to fluctuate between an octahedral six-coordinated and a seven-coordinated structure due to the increased local proximity of Cl in the eutectic salt, resulting in a lower fluidicity/diffusivity than the other trivalent lanthanides studied. The thermophysical properties of the eutectic composition with trivalent lanthanides were not significantly different from the pure eutectic salt composition, but several changes were noted.
Collapse
Affiliation(s)
- Jon Fuller
- Department of Chemical and Materials Engineering, University of Nevada - Reno, Nevada 89577, USA.
| | - Jeremy Moon
- Department of Chemical and Materials Engineering, University of Nevada - Reno, Nevada 89577, USA.
| | - Jing Zhang
- Department of Chemical and Materials Engineering, University of Nevada - Reno, Nevada 89577, USA.
| | - Dev Chidambaram
- Department of Chemical and Materials Engineering, University of Nevada - Reno, Nevada 89577, USA.
| | - Qi An
- Department of Chemical and Materials Engineering, University of Nevada - Reno, Nevada 89577, USA.
| |
Collapse
|
3
|
Li S, Pokuri BSS, Ryno SM, Nkansah A, De'Vine C, Ganapathysubramanian B, Risko C. Determination of the Free Energies of Mixing of Organic Solutions through a Combined Molecular Dynamics and Bayesian Statistics Approach. J Chem Inf Model 2020; 60:1424-1431. [PMID: 31935097 DOI: 10.1021/acs.jcim.9b01113] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
As new generations of thin-film semiconductors are moving toward solution-based processing, the development of printing formulations will require information pertaining to the free energies of mixing of complex mixtures. From the standpoint of in silico material design, this move necessitates the development of methods that can accurately and quickly evaluate these formulations in order to maximize processing speed and reproducibility. Here, we make use of molecular dynamics (MD) simulations, in combination with the two-phase thermodynamic (2PT) model, to explore the free energy of mixing surfaces for a series of halogenated solvents and high-boiling point solvent additives used in the development of thin-film organic semiconductors. Although the combined methods generally show good agreement with available experimental data, the computational cost to traverse the free-energy landscape is considerable. Hence, we demonstrate how a Bayesian optimization scheme, coupled with the MD and 2PT approaches, can drastically reduce the number of simulations required, in turn shrinking both the computational cost and time.
Collapse
Affiliation(s)
- Shi Li
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | | | - Sean M Ryno
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | - Asare Nkansah
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | - Camron De'Vine
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| | | | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research (CAER), University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
4
|
Tiwari A, Honingh C, Ensing B. Accurate calculation of zero point energy from molecular dynamics simulations of liquids and their mixtures. J Chem Phys 2019; 151:244124. [PMID: 31893925 DOI: 10.1063/1.5131145] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
The two-phase thermodynamic (2PT) method is used to compute the zero point energy (ZPE) of several liquids and their mixtures. The 2PT method uses the density of states (DoS), which is computed from the velocity autocorrelation (VAC) function obtained from a short classical molecular dynamics trajectory. By partitioning the VAC and the DoS of a fluid into solid and gaslike components, quantum mechanical corrections to thermodynamical properties can be computed. The ZPE is obtained by combining the partition function of the quantum harmonic oscillator with the vibrational part of the solidlike DoS. The resulting ZPE is found to be in excellent agreement with both experimental and ab initio results. Solvent effects such as hydrogen bonding and polarization can be included by the utilization of ab initio density functional theory based molecular dynamics simulations. It is found that these effects significantly influence the DoS of water molecules. The obtained results demonstrate that the 2PT model is a powerful method for efficient ZPE calculations, in particular, to account for solvent effects and polarization.
Collapse
Affiliation(s)
- A Tiwari
- Van't Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - C Honingh
- Van't Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - B Ensing
- Van't Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
5
|
|
6
|
Dhabal D, Nguyen AH, Singh M, Khatua P, Molinero V, Bandyopadhyay S, Chakravarty C. Excess entropy and crystallization in Stillinger-Weber and Lennard-Jones fluids. J Chem Phys 2015; 143:164512. [DOI: 10.1063/1.4933420] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi 110016, India
| | - Andrew Huy Nguyen
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Murari Singh
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Prabir Khatua
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Valeria Molinero
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Sanjoy Bandyopadhyay
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi 110016, India
| |
Collapse
|
7
|
Chen M, Pendrill R, Widmalm G, Brady JW, Wohlert J. Molecular Dynamics Simulations of the Ionic Liquid 1-n-Butyl-3-Methylimidazolium Chloride and Its Binary Mixtures with Ethanol. J Chem Theory Comput 2014; 10:4465-79. [DOI: 10.1021/ct500271z] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Mo Chen
- Department
of Food Science, Cornell University, Ithaca, New York 14853, United States
| | - Robert Pendrill
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| | - Göran Widmalm
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| | - John W. Brady
- Department
of Food Science, Cornell University, Ithaca, New York 14853, United States
| | - Jakob Wohlert
- Wallenberg
Wood Science Center, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
| |
Collapse
|
8
|
Yeh KY, Huang SN, Chen LJ, Lin ST. Diffusive and quantum effects of water properties in different states of matter. J Chem Phys 2014; 141:044502. [DOI: 10.1063/1.4890572] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Kuan-Yu Yeh
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
- Refining & Manufacturing Research Institute, CPC Corporation, Chia-Yi 60051, Taiwan
| | - Shao-Nung Huang
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
| | - Li-Jen Chen
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
| | - Shiang-Tai Lin
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|