1
|
Windom ZW, Datta M, Huda MM, Sabuj MA, Rai N. Understanding speciation and solvation of glyphosate from first principles simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
2
|
Huda MM, Saha C, Jahan N, Wilson WN, Rai N. Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores. J Phys Chem B 2022; 126:1352-1364. [PMID: 35119855 DOI: 10.1021/acs.jpcb.1c10572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (∼20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. The dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule's local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.
Collapse
Affiliation(s)
- Md Masrul Huda
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Chinmoy Saha
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Nusrat Jahan
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Woodrow N Wilson
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, United States
| |
Collapse
|
3
|
Huda MM, Rai N. Probing Early-Stage Aggregation of Low Molecular Weight Gelator in an Organic Solvent. J Phys Chem B 2020; 124:2277-2288. [PMID: 32105082 DOI: 10.1021/acs.jpcb.9b11865] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Md Masrul Huda
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State 39762, Mississippi, United States
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State 39762, Mississippi, United States
| |
Collapse
|
4
|
Goel H, Windom ZW, Jackson AA, Rai N. CO2 sorption in triethyl(butyl)phosphonium 2-cyanopyrrolide ionic liquid via first principles simulations. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111323] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
5
|
Goel H, Ling S, Ellis BN, Taconi A, Slater B, Rai N. Predicting vapor liquid equilibria using density functional theory: A case study of argon. J Chem Phys 2018; 148:224501. [PMID: 29907054 DOI: 10.1063/1.5025726] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Sanliang Ling
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Breanna Nicole Ellis
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Anna Taconi
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| |
Collapse
|
6
|
Goel H, Rai N. Transferable Potentials for Chloroethenes: Insights into Nonideal Solution Behavior of Environmental Contaminants. ACS OMEGA 2018; 3:3646-3654. [PMID: 31458615 PMCID: PMC6641527 DOI: 10.1021/acsomega.8b00044] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Accepted: 03/19/2018] [Indexed: 06/10/2023]
Abstract
Predicting the nonideal phase behavior of binary and multicomponent systems remains a significant challenge for particle-based simulations. Here, we develop a transferable force field for chloroethenes, common environmental contaminants, that can accurately model the vapor liquid phase equilibria including azeotrope formation. The new all-atom force field reproduces saturated liquid densities, saturated vapor pressures, boiling points, and critical properties within 1, 10, 1, and 1% of the experiment data, respectively. Furthermore, the vapor liquid equilibria of trichloroethylene and 1-propanol binary mixture, which forms a minimum boiling point azeotrope, is predicted with a reasonable accuracy. The microstructure of neat and binary systems is explored using pair correlation functions and spatial distribution functions. As the new force field is consistent with transferable potentials for phase equilibria (TraPPE) force field, it expands the applicability of TraPPE force field to chloroethenes.
Collapse
Affiliation(s)
| | - Neeraj Rai
- E-mail: . Phone: +1
(662) 3250790. Fax: +1 (662) 3252482 (N.R.)
| |
Collapse
|
7
|
Goel H, Windom ZW, Jackson AA, Rai N. Performance of density functionals for modeling vapor liquid equilibria of CO 2 and SO 2. J Comput Chem 2018; 39:397-406. [PMID: 29164642 DOI: 10.1002/jcc.25123] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Revised: 10/12/2017] [Accepted: 11/07/2017] [Indexed: 01/16/2023]
Abstract
Vapor liquid equilibria (VLE) and condensed phase properties of carbon dioxide and sulfur dioxide are calculated using first principles Monte Carlo (FPMC) simulations to assess the performance of several density functionals, notably PBE-D3, BLYP-D3, PBE0-D3, M062X-D3, and rVV10. GGA functionals were used to compute complete vapor liquid coexistence curves (VLCCs) to estimate critical properties, while the hybrid and nonlocal van der Waals functionals were used only for computing density at a single state point due to the high computational cost. Our results show that the BLYP-D3 functional performs well in predicting VLE properties for both molecules when compared with other functionals. In the liquid phase, pair correlation functions reveal that there is not a significant difference in the location of the peak for the first solvation shell while the peak heights are different for different functionals. Overall, the BLYP-D3 functional is a good choice for modeling VLE of acidic gases with significant environmental implications such as CO2 and SO2 . © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Zachary W Windom
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Amber A Jackson
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS, 39762
| |
Collapse
|
8
|
Goel H, Windom ZW, Butler CL, Rai N. Phase Equilibria and Condensed Phase Properties of Fluorinated Alkanes via First Principles Simulations. ChemistrySelect 2017. [DOI: 10.1002/slct.201701972] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems Mississippi State University, MS
| | - Zachary W. Windom
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems Mississippi State University, MS
| | - Charles L. Butler
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems Mississippi State University, MS
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering and Center for Advanced Vehicular Systems Mississippi State University, MS
| |
Collapse
|