1
|
Maya J, Malloum A, Fifen JJ, Dhaouadi Z, Fouda HPE, Conradie J. Quantum cluster equilibrium theory applied to liquid ammonia. J Comput Chem 2024; 45:1279-1288. [PMID: 38353541 DOI: 10.1002/jcc.27327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 01/21/2024] [Accepted: 01/29/2024] [Indexed: 04/19/2024]
Abstract
Through this paper, the authors propose using the quantum cluster equilibrium (QCE) theory to reinvestigate ammonia clusters in the liquid phase. The ammonia clusters from size monomer to hexadecamer were considered to simulate the liquid ammonia in this approach. The clusterset used to model the liquid ammonia is an ensemble of different structures of ammonia clusters. After studious research of the representative configurations of ammonia clusters through the cluster research program ABCluster, the configurations have been optimized at the MN15/6-31++G(d,p) level of theory. These optimizations lead to geometries and frequencies as inputs for the Peacemaker code. The QCE study of this molecular system permits us to get the liquid phase populations in a temperature range of 190-260 K, covering the temperatures from the melting point to the boiling point. The results show that the population of liquid ammonia comprises mainly the ammonia hexadecamer followed by pentadecamer, tetradecamer, and tridecamer. We noted that the small-sized ammonia clusters do not contribute to the population of liquid ammonia. In addition, the thermodynamic properties, such as heat of vaporization, heat capacity, entropy, enthalpy, and free energies, obtained by the QCE theory have been compared to the experiment given some relatively good agreements in the gas phase and show considerable discrepancies in liquid phase except the density. Finally, based on the predicted population, we calculated the infrared spectrum of liquid ammonia at 215 K temperature. It comes out that the calculated infrared spectrum qualitatively agrees with the experiment.
Collapse
Affiliation(s)
- Josué Maya
- Department of Physics, Faculty of Science, University of Ngaoundere, Ngaoundere, Cameroon
- National Radiation Protection Agency, Yaounde, Cameroon
| | - Alhadji Malloum
- Department of Physics, Faculty of Science, University of Maroua, Maroua, Cameroon
- Department of Chemistry, University of the Free State, Bloemfontein, South Africa
| | - Jean Jules Fifen
- Department of Physics, Faculty of Science, University of Ngaoundere, Ngaoundere, Cameroon
| | - Zoubeida Dhaouadi
- Laboratoire de Spectroscopie Atomique Moléculaire et Application, Université de Tunis El Manar, Tunis, Tunisie
| | | | - Jeanet Conradie
- Department of Chemistry, University of the Free State, Bloemfontein, South Africa
| |
Collapse
|
2
|
Khire SS, Nakajima T, Gadre SR. REAlgo: Rapid and efficient algorithm for estimating MP2/CCSD energy gradients for large molecular clusters. J Chem Phys 2023; 159:184109. [PMID: 37955321 DOI: 10.1063/5.0174726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Accepted: 10/18/2023] [Indexed: 11/14/2023] Open
Abstract
This work reports the development of an algorithm for rapid and efficient evaluation of energy gradients for large molecular clusters employing correlated methods viz. second-order Møller-Plesset perturbation theory (MP2) theory and couple cluster singles and doubles (CCSD). The procedure segregates the estimation of Hartree-Fock (HF) and correlation components. The HF energy and gradients are obtained by performing a full calculation. The correlation energy is approximated as the corresponding two-body interaction energy. Correlation gradients for each monomer are approximated from the respective monomer-centric fragments comprising its immediate neighbours. The programmed algorithm is explored for the geometry optimization of large molecular clusters using the BERNY optimizer as implemented in the Gaussian suite of software. The accuracy and efficacy of the method are critically probed for a variety of large molecular clusters containing up to 3000 basis functions, in particular large water clusters. The CCSD level geometry optimization of molecular clusters containing ∼800 basis functions employing a modest hardware is also reported.
Collapse
Affiliation(s)
- Subodh S Khire
- RIKEN Center for Computational Science, Kobe 6500047, Japan
| | | | - Shridhar R Gadre
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune 411007, India
- Department of Chemistry, Savitribai Phule Pune University, Pune 411007, India
| |
Collapse
|
3
|
Da-yang TE, Lai CH. Potential energy surfaces of the Cu2+(NH3)n=1−10 clusters in solvent phase: A DFT study. Chem Phys 2023. [DOI: 10.1016/j.chemphys.2023.111902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
|
4
|
Patkar D, Bharati Ahirwar M, Deshmukh MM. A Tug of War between the Self- and Cross-associating Hydrogen Bonds in Neutral Ammonia-Water Clusters: Energetic Insights by Molecular Tailoring Approach. Chemphyschem 2022; 23:e202200476. [PMID: 36127809 DOI: 10.1002/cphc.202200476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Revised: 09/09/2022] [Indexed: 01/05/2023]
Abstract
In the present work, the energies of various types of individual HBs observed in neutral (NH3 )m (H2 O)n , (m+n=2 to 7) clusters were estimated using the molecular tailoring approach (MTA)-based method. The calculated individual HB energies suggest that the O-H…N HBs are the strongest (1.21 to 12.49 kcal mol-1 ). The next ones are the O-H…O (3.97 to 9.30 kcal mol-1 ) HBs. The strengths of N-H…N (1.09 to 5.29 kcal mol-1 ) and N-H…O (2.85 to 5.56 kcal mol-1 ) HBs are the weakest. The HB energies in dimers also follow this rank ordering. However, the HB energies in dimers are much smaller than those obtained by the MTA-based method due to the loss in cooperativity contribution in the dimers. Thus, the calculated cooperativity contributions, for different types of HBs, fall in the range 0.64 to 5.73 kcal mol-1 . We wish to emphasize based on the energetic rank ordering obtained by the MTA-based method that the O-H of water is a better HB donor than the N-H of ammonia. The reasons for the observed energetic rank ordering are two folds: (i) intrinsically stronger O-H…N HBs than the O-H…O ones as revealed by dimer energies and (ii) the higher cooperativity contribution in the former than the later ones. Indeed, the MTA-based method is useful in providing the missing energetic rank ordering of various type of HBs in neutral (NH3 )m (H2 O)n clusters, in the literature.
Collapse
Affiliation(s)
- Deepak Patkar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya, (A Central University), 470003, Sagar, India
| |
Collapse
|
5
|
Nakamura T, Fedorov DG. The catalytic activity and adsorption in faujasite and ZSM-5 zeolites: the role of differential stabilization and charge delocalization. Phys Chem Chem Phys 2022; 24:7739-7747. [PMID: 35293902 DOI: 10.1039/d1cp05851g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Adsorption and chemical reactions occurring on industrially important ZSM-5 and faujasite zeolite catalysts are investigated with the quantum-mechanical fragment molecular orbital method combined with periodic boundary conditions. Suitable fragmentation patterns are devised and tested providing important case studies of computing real materials with fragmentation methods. A good accuracy is demonstrated in comparison to full calculations, and a good agreement with the available experimental data is obtained. The full production cycle of p-xylene on faujasite zeolite is mapped. The catalytic role of the zeolite in the dehydration reaction, analyzed with the partition analysis, is attributed to the delocalization of the negative charge over the zeolite. On the other hand, an increase of the barrier in the Diels-Alder reaction by the zeolite is attributed to the preferential stabilization of the reactants over the transition state as demonstrated by the guest-zeolite interaction energy.
Collapse
Affiliation(s)
- Taiji Nakamura
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| |
Collapse
|
6
|
Malloum A, Conradie J. Non-covalent interactions in dimethylsulfoxide (DMSO) clusters and DFT benchmarking. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118522] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
7
|
Nakamura T, Yokaichiya T, Fedorov DG. Analysis of Guest Adsorption on Crystal Surfaces Based on the Fragment Molecular Orbital Method. J Phys Chem A 2022; 126:957-969. [PMID: 35080391 DOI: 10.1021/acs.jpca.1c10229] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For gaining insights into interactions in periodic systems, an analysis is developed based on the fragment molecular orbital method combined with periodic boundary conditions. The adsorption energy is decomposed into guest and surface polarization and deformation energy, guest-surface and guest-guest interactions, and the vibrational free energy. The analysis is applied to the adsorption of guest molecules to Ih (001) ice surface. The cooperativity effects result in a non-linear change in the adsorption energy with coverage due to many-body effects. The role of dispersion is found to be dominant for guests with long hydrophobic tails. A rule is proposed relating the length of the alkyl tail with the formation of the guest layer. The computed binding enthalpies are in good agreement with experimental values. For high coverage, adsorbed molecules can form an ordered layer known as self-assembled monolayer.
Collapse
Affiliation(s)
- Taiji Nakamura
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Tomoko Yokaichiya
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
8
|
|
9
|
Ahirwar MB, Patkar D, Yadav I, Deshmukh MM. Appraisal of individual hydrogen bond strengths and cooperativity in ammonia clusters via a molecular tailoring approach. Phys Chem Chem Phys 2021; 23:17224-17231. [PMID: 34369546 DOI: 10.1039/d1cp02839a] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we propose and test a method, based on the molecular tailoring approach (MTA), for the evaluation of individual hydrogen bond (HB) energies in ammonia (NH3)n clusters. This methodology was tested, in our earlier work, on water clusters. Liquid ammonia being a universal, non-aqueous ionizing solvent, such information of individual HB strength is indispensable in many studies. The estimated HB energies by an MTA-based method, in (NH3)n for n = 3-8, were calculated to be in the range of 0.65 to 5.54 kcal mol-1 with the cooperativity contribution falling between -0.54 and 1.88 kcal mol-1 both calculated at the MP2(full)/aug-cc-pVTZ level of theory. It is seen that the strong HBs in (NH3)n clusters were additionally strengthened by the large contribution of HB cooperativity. The accuracy of these estimated HB energies was validated by approximately estimating the molecular energy of a given cluster by adding the sum of HB energies to the sum of monomer energies. This approximately estimated molecular energy of a given cluster was found to be in excellent agreement with the actual calculated values. The negligibly small difference (less than 5.6 kcal mol-1) in these two values suggests that the estimated individual HB energies in ammonia clusters are quite reliable. Furthermore, these estimated HB energies by MTA are in excellent qualitative agreement with the other indirect measures of HB strength, such as HB bond distances and angles, N-H stretching frequency and the electron density values at the (3,-1) bond critical points.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr Harisingh Gour Vishwavidyalaya, (A Central University), Sagar, 470003, India.
| | | | | | | |
Collapse
|