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
|
Yang S, Zhang Z, Zhang DH. A full-dimensional ab initio potential energy and dipole moment surfaces for (NH 3) 2. J Chem Phys 2021; 155:164306. [PMID: 34717358 DOI: 10.1063/5.0072063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) for the ammonia dimer (NH3)2 are reported. The database of the PES consists of 27 736 ab initio energy points and all of these points were calculated at the UCCSD(T)-F12a/AVTZ level. The PES was fitted by using the fundamental invariant neural network (FI-NN) method that satisfies the permutational symmetry of identical atoms, and the root mean square fitting error for the PES is very small as low as 0.562 meV. The geometries for the (NH3)2 DMS are the same as those used for the PES and are calculated at the XYG3/AVTZ level. This PES can describe a variety of internal floppy motions, including all kinds of vibrational modes no matter intermolecular or intramolecular. The CCSD(T)-PES can dissociate correctly to two NH3 monomers, with De = 1135.55 cm-1 (13.58 kJ/mol) which agrees accurately with the 13.5 ± 0.3 kJ/mol predicted by previous work.
Collapse
Affiliation(s)
- Shuo Yang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Zhaojun Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| |
Collapse
|
3
|
Malloum A, Conradie J. Hydrogen bond networks of ammonia clusters: What we know and what we don’t know. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116199] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
4
|
Ninomiya M, Doi H, Matsumoto Y, Mochizuki Y, Komeiji Y. Ab Initio Fragment Molecular Orbital-Based Molecular Dynamics (FMO-MD) Simulations of (NH 3) 32 Cluster: Effects of Electron Correlation. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2020. [DOI: 10.1246/bcsj.20190320] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Moeko Ninomiya
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Hideo Doi
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology, AIST, Tsukuba Central 2, Tsukuba, Ibaraki 305-8568, Japan
| | - Yoshiteru Matsumoto
- Department of Chemistry, Faculty of Science, Shizuoka University, 836 Ohya, Suruga, Shizuoka 422-8529, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Biomedical Research Institute, National Institute of Advanced Industrial Science and Technology, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
| |
Collapse
|
5
|
Malloum A, Fifen JJ, Conradie J. Large-Sized Ammonia Clusters and Solvation Energies of the Proton in Ammonia. J Comput Chem 2020; 41:21-30. [PMID: 31568565 DOI: 10.1002/jcc.26071] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2019] [Revised: 08/19/2019] [Accepted: 08/28/2019] [Indexed: 01/10/2023]
Abstract
The absolute solvation energies (free energies and enthalpies) of the proton in ammonia are used to compute the pKa of species embedded in ammonia. They are also used to compute the solvation energies of other ions in ammonia. Despite their importance, it is not possible to determine experimentally the solvation energies of the proton in a given solvent. We propose in this work a direct approach to compute the solvation energies of the proton in ammonia from large-sized neutral and protonated ammonia clusters. To undertake this investigation, we performed a geometry optimization of neutral and protonated ammonia 30-mer, 40-mer, and 50 mer to locate stable structures. These structures have been fully optimized at both APFD/6-31++g(d,p) and M06-2X/6-31++g(d,p) levels of theory. An infrared spectroscopic study of these structures has been provided to assess the reliability of our investigation. Using these structures, we have computed the absolute solvation free energy and the absolute solvation enthalpy of the proton in ammonia. It comes out that the absolute solvation free energy of the proton in ammonia is calculated to be -1192 kJ mol-1 , whereas the absolute solvation enthalpy is evaluated to be -1214 kJ mol-1 . © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Alhadji Malloum
- Department of Chemistry, University of the Free State, Bloemfontein, 9300, South Africa
| | - Jean J Fifen
- Department of Physics, Faculty of Science, The University of Ngaoundere, 454, Ngaoundere, Cameroon
| | - Jeanet Conradie
- Department of Chemistry, University of the Free State, Bloemfontein, 9300, South Africa
| |
Collapse
|
6
|
Malloum A, Fifen JJ, Conradie J. Structures and spectroscopy of the ammonia eicosamer, (NH3)n=20. J Chem Phys 2018; 149:024304. [DOI: 10.1063/1.5031790] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Alhadji Malloum
- Faculty of Science, Department of Physics, University of Ngaoundere, P.O. Box 454, Ngaoundere, Cameroon
| | - Jean Jules Fifen
- Faculty of Science, Department of Physics, University of Ngaoundere, P.O. Box 454, Ngaoundere, Cameroon
| | - Jeanet Conradie
- Department of Chemistry, University of the Free State, P.O. Box 339, Bloemfontein 9300, South Africa
| |
Collapse
|
7
|
Mella M, Curotto E. Quest for Inexpensive Hydrogen Isotopic Fractionation: Do We Need 2D Quantum Confining in Porous Materials or Are Rough Surfaces Enough? The Case of Ammonia Nanoclusters. J Phys Chem A 2016; 120:8148-8159. [PMID: 27704841 DOI: 10.1021/acs.jpca.6b08005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We study the adsorption energetics and quantum properties of the molecular hydrogen isotopes H2, D2, and T2 onto the surface of rigid ammonia nanoclusters with quantum simulations and accurate model potential energy surfaces (PES). A highly efficient diffusion Monte Carlo (DMC) algorithm for rigid rotors allowed us to accurately define zero-point adsorption energies for the three isotopes, as well as the degree of translational and rotational delocalization that each affords on the surface. From the data emerges that the quantum adsorption energy (Eads) of T2 can be up to twice the one of H2 at 0 K, suggesting the possibility of exploiting some form of solid ammonia to selectivity separate hydrogen isotopes at low temperatures (≃20 K). This is discussed by focusing on the structural motif that may be more effective for the task. The analysis of the contributions to Eads, however, surprisingly indicates that the average kinetic energy (Ekin) and rotation energy (Erotkin) of T2 can also be, respectively, 2 times and 20 times higher than those of H2; this finding markedly deviates from what is predicted for hydrogen molecules inside carbon nanotubes (CNT) or metallic-organic frameworks (MOF), where Ekin and Erotkin is higher for H2 due to the unavoidable effects of confinement and hindrance to its rotational motion. The rationale for these differences is provided by the geometrical distributions for the rigid rotors, which reveal an increasingly stronger coupling between rotational and translational degrees of freedom upon increasing the isotopic mass. This effect has never been observed before on adsorbing surfaces (e.g., graphite) and is induced by a strongly anisotropic and anharmonic bowl-like potential experienced by the rotors.
Collapse
Affiliation(s)
- Massimo Mella
- Dipartimento di Scienza ed Alta Tecnologia, Università degli Studi dell'Insubria , via Valleggio 9, 22100 Como, Italy
| | - E Curotto
- Department of Chemistry and Physics, Arcadia University , Glenside, Pennsylvania 19038-3295, United States
| |
Collapse
|
8
|
|
9
|
Malloum A, Fifen JJ, Dhaouadi Z, Engo SGN, Jaidane NE. Structures and relative stabilities of ammonia clusters at different temperatures: DFT vs. ab initio. Phys Chem Chem Phys 2015; 17:29226-42. [DOI: 10.1039/c5cp03374h] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The global minimum energy structures of (NH3)n=2–10are pointed out for the first time at a given temperature.
Collapse
Affiliation(s)
- Alhadji Malloum
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere 454
- Cameroon
| | - Jean Jules Fifen
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere 454
- Cameroon
| | - Zoubeida Dhaouadi
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications
- Faculté des Sciences de Tunis
- Université de Tunis El Manar
- Tunis
- Tunisie
| | - Serge Guy Nana Engo
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere 454
- Cameroon
| | - Nejm-Eddine Jaidane
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications
- Faculté des Sciences de Tunis
- Université de Tunis El Manar
- Tunis
- Tunisie
| |
Collapse
|
10
|
|
11
|
|
12
|
Jenkins R, Curotto E, Mella M. Replica exchange with Smart Monte Carlo and Hybrid Monte Carlo in manifolds. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.10.053] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
13
|
Mella M, Curotto E. Quantum simulations of the hydrogen molecule on ammonia clusters. J Chem Phys 2013; 139:124319. [PMID: 24089779 DOI: 10.1063/1.4821648] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Mixed ammonia-hydrogen molecule clusters [H2-(NH3)n] have been studied with the aim of exploring the quantitative importance of the H2 quantum motion in defining their structure and energetics. Minimum energy structures have been obtained employing genetic algorithm-based optimization methods in conjunction with accurate pair potentials for NH3-NH3 and H2-NH3. These include both a full 5D potential and a spherically averaged reduced surface mimicking the presence of a para-H2. All the putative global minima for n ≥ 7 are characterized by H2 being adsorbed onto a rhomboidal ammonia tetramer motif formed by two double donor and two double acceptor ammonia molecules. In a few cases, the choice of specific rhombus seems to be directed by the vicinity of an ammonia ad-molecule. Diffusion Monte Carlo simulations on a subset of the species obtained highlighted important quantum effects in defining the H2 surface distribution, often resulting in populating rhomboidal sites different from the global minimum one, and showing a compelling correlation between local geometrical features and the relative stability of surface H2. Clathrate-like species have also been studied and suggested to be metastable over a broad range of conditions if formed.
Collapse
Affiliation(s)
- Massimo Mella
- Dipartimento di Scienze ed Alta Tecnologia, Università degli Studi dell'Insubria, via Valleggio 11, 22100 Como, Italy
| | | |
Collapse
|
14
|
Unn-Toc W, Halberstadt N, Meier C, Mella M. Exploring the importance of quantum effects in nucleation: the archetypical Ne(n) case. J Chem Phys 2012; 137:014304. [PMID: 22779645 DOI: 10.1063/1.4730033] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effect of quantum mechanics (QM) on the details of the nucleation process is explored employing Ne clusters as test cases due to their semi-quantal nature. In particular, we investigate the impact of quantum mechanics on both condensation and dissociation rates in the framework of the microcanonical ensemble. Using both classical trajectories and two semi-quantal approaches (zero point averaged dynamics, ZPAD, and Gaussian-based time dependent Hartree, G-TDH) to model cluster and collision dynamics, we simulate the dissociation and monomer capture for Ne(8) as a function of the cluster internal energy, impact parameter and collision speed. The results for the capture probability P(s)(b) as a function of the impact parameter suggest that classical trajectories always underestimate capture probabilities with respect to ZPAD, albeit at most by 15%-20% in the cases we studied. They also do so in some important situations when using G-TDH. More interestingly, dissociation rates k(diss) are grossly overestimated by classical mechanics, at least by one order of magnitude. We interpret both behaviours as mainly due to the reduced amount of kinetic energy available to a quantum cluster for a chosen total internal energy. We also find that the decrease in monomer dissociation energy due to zero point energy effects plays a key role in defining dissociation rates. In fact, semi-quantal and classical results for k(diss) seem to follow a common "corresponding states" behaviour when the proper definition of internal and dissociation energies are used in a transition state model estimation of the evaporation rate constants.
Collapse
Affiliation(s)
- Wesley Unn-Toc
- Laboratoire Collisions Agrégats Réactivité-IRSAMC, Université Paul Sabatier, 118 Route de Narbonne, 31062 Toulouse, France
| | | | | | | |
Collapse
|
15
|
Roberts K, Sebsebie R, Curotto E. A rare event sampling method for diffusion Monte Carlo using smart darting. J Chem Phys 2012; 136:074104. [PMID: 22360233 DOI: 10.1063/1.3685453] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We identify a set of multidimensional potential energy surfaces sufficiently complex to cause both the classical parallel tempering and the guided or unguided diffusion Monte Carlo methods to converge too inefficiently for practical applications. The mathematical model is constructed as a linear combination of decoupled Double Wells [(DDW)(n)]. We show that the set (DDW)(n) provides a serious test for new methods aimed at addressing rare event sampling in stochastic simulations. Unlike the typical numerical tests used in these cases, the thermodynamics and the quantum dynamics for (DDW)(n) can be solved deterministically. We use the potential energy set (DDW)(n) to explore and identify methods that can enhance the diffusion Monte Carlo algorithm. We demonstrate that the smart darting method succeeds at reducing quasiergodicity for n ≫ 100 using just 1 × 10(6) moves in classical simulations (DDW)(n). Finally, we prove that smart darting, when incorporated into the regular or the guided diffusion Monte Carlo algorithm, drastically improves its convergence. The new method promises to significantly extend the range of systems computationally tractable by the diffusion Monte Carlo algorithm.
Collapse
Affiliation(s)
- K Roberts
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | | | | |
Collapse
|
16
|
Abbaspour M. Transport, thermodynamic, and structural properties of fluid ammonia using a new intermolecular potential: The inversion method and molecular dynamics simulation. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.09.017] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
17
|
Patrone M, Mella M. Sequential growth simulation of (NH3)n clusters (n=2–8) in ultracold superfluid environment. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.07.086] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
18
|
Curotto E, Mella M. Quantum Monte Carlo simulations of selected ammonia clusters (n = 2–5): Isotope effects on the ground state of typical hydrogen bonded systems. J Chem Phys 2010; 133:214301. [DOI: 10.1063/1.3506027] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
|