1
|
Prosmiti R, González-Lezana T. Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems. Chemphyschem 2024:e202400207. [PMID: 38837591 DOI: 10.1002/cphc.202400207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Indexed: 06/07/2024]
Abstract
The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.
Collapse
Affiliation(s)
- Rita Prosmiti
- Instituto de Fìsica Fundamental, IFF-CSIC, Serrano 123, 28006, Madrid, Spain
| | | |
Collapse
|
2
|
Alharzali N, Rodríguez-Segundo R, Prosmiti R. Modelling interactions of cationic dimers in He droplets: microsolvation trends in He nK 2+ clusters. Phys Chem Chem Phys 2021; 23:7849-7859. [PMID: 33220666 DOI: 10.1039/d0cp05406b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report the results of a detailed theoretical investigation of small K2+-doped He clusters. The structural characteristics and stabilities of such cations are determined from ab initio electronic structure calculations at the MRCI+Q level of theory. The underlying interactions show a multireference character and such effects are analyzed. The interaction potentials are constructed employing an interpolation technique within the inverse problem theory method, while the nuclear quantum effects are computed for the trimers, their spatial arrangements are discussed, and information was extracted on the orientational anisotropy of the forces. We found that energetically the most stable conformer corresponds to linear arrangements that are taking place under large amplitude vibrations, with high zero-point energy. We have further looked into the behavior of higher-order species with various He atoms surrounding the cationic dopant. By using a sum of potentials approach and an evolutionary programming method, we analyzed the structural stability of clusters with up to six He atoms in comparison with interactions energies obtained from MRCI+Q quantum chemistry computations. Structures containing Hen motifs that characterize pure rare gas clusters, appear for the larger K2+-doped He clusters, showing selective growth during the microsolvation process of the alkali-dimer cation surrounded by He atoms. Such results indicate the existence of local solvation microstructures in these aggregates, where the cationic impurity could get trapped for a short time, contributing to the slow ionic mobility observed experimentally in ultra-cold He-droplets.
Collapse
Affiliation(s)
- Nissrin Alharzali
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| | | | | |
Collapse
|
3
|
Borocci S, Grandinetti F, Sanna N, Antoniotti P, Nunzi F. Complexes of helium with neutral molecules: Progress toward a quantitative scale of bonding character. J Comput Chem 2020; 41:1000-1011. [PMID: 31960984 DOI: 10.1002/jcc.26146] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2019] [Revised: 12/07/2019] [Accepted: 01/01/2020] [Indexed: 11/12/2022]
Abstract
The complexes of helium with nearly 30 neutral molecules (M) were investigated by various techniques of bonding analysis and symmetry-adapted perturbation theory (SAPT). The main investigated function was the local electron energy density H(r), analyzed, in particular, so to estimate the degree of polarization (DoP) of He in the various He(M). As we showed recently (Borocci et al., J. Comput. Chem., 2019, 40, 2318-2328), the DoP is a quantitative index that is generally informative about the role of polarization (induction plus charge transfer [CT]) and dispersion in noncovalent noble gas complexes. As further evidence in this regard, we presently ascertained quantitative correlations between the DoP(He) of the He(M) and indices based on the electron density ρ(r), including the molecular electrostatic potential at the HeM bond critical point, as well as the percentage contributions of induction and dispersion to the SAPT binding energies. Based also on the explicit evaluation of the CT, accomplished through the study of the charge-displacement function, we derived a quantitative scale that ranks the He(M) according to their dispersive, inductive, and CT bonding character. Our taken approach could be conceivably extended to other types of noncovalent complexes.
Collapse
Affiliation(s)
- Stefano Borocci
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., Viterbo, Italy.,Istituto per i Sistemi Biologici del CNR, Monterotondo, RM Italy
| | - Felice Grandinetti
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., Viterbo, Italy.,Istituto per i Sistemi Biologici del CNR, Monterotondo, RM Italy
| | - Nico Sanna
- Dipartimento per la Innovazione nei sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, L.go dell'Università, s.n.c., Viterbo, Italy
| | | | - Francesca Nunzi
- Dipartimento di Chimica, Biologia e Biotecnologie (DCBB), Istituto CNR di Scienze e Tecnologie Chimiche "Giulio Natta" (CNR-SCITEC), Perugia, Italy
| |
Collapse
|
4
|
Papp D, Sarka J, Szidarovszky T, Császár AG, Mátyus E, Hochlaf M, Stoecklin T. Complex rovibrational dynamics of the Ar·NO + complex. Phys Chem Chem Phys 2017; 19:8152-8160. [PMID: 28225106 DOI: 10.1039/c6cp07731e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rotational-vibrational states of the Ar·NO+ cationic complex are computed, below, above, and well above the complex's first dissociation energy, using variational nuclear motion and close-coupling scattering computations. The HSLH potential energy surface used in this study (J. Chem. Phys., 2011, 135, 044312) is characterized by a first dissociation energy of D0 = 887.0 cm-1 and supports 200 bound vibrational states. The bound-state vibrational energies and the corresponding wave functions allow the interpretation of the scarcely available experimental results about the intermonomer vibrational motion of the complex. A very large number of long-lived quasibound combination states of the three vibrational modes, exhibiting a very similar energy-level structure as that of the bound states, are found embedded in the continuum. Additional short-lived resonance states are also identified and their properties are analyzed.
Collapse
Affiliation(s)
- Dóra Papp
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary.
| | - János Sarka
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary.
| | - Tamás Szidarovszky
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary.
| | - Attila G Császár
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary.
| | - Edit Mátyus
- Institute of Chemistry, Eötvös University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Majdi Hochlaf
- Université Paris-Est, Laboratoire Modélisation et Simulation Multi Echelle, MSME UMR 8208 CNRS, 5 bd Descartes, F-77454 Marne-la-Vallée, France
| | - Thierry Stoecklin
- Institut des Sciences Moléculaires, Université de Bordeaux, CNRS UMR 5255, 33405 Talence Cedex, France.
| |
Collapse
|
5
|
Li S, Zheng R, Chen SJ, Chen Y, Chen P. Investigations of the Rg-BrCl (Rg=He, Ne, Ar, Kr, Xe) binary van der Waals complexes: ab initio intermolecular potential energy surfaces, vibrational states and predicted pure rotational transition frequencies. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2017; 174:105-117. [PMID: 27888780 DOI: 10.1016/j.saa.2016.11.020] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2015] [Revised: 11/09/2016] [Accepted: 11/13/2016] [Indexed: 06/06/2023]
Abstract
The intermolecular potential energy surfaces (PESs) of the ground electronic state for the Rg-BrCl (Rg=He, Ne, Ar, Kr, Xe) van der Waals complexes have been constructed by using the coupled-cluster method in combination with the augmented quadruple-zeta correlation-consistent basis sets supplemented with an additional set of bond functions. The features of the anisotropic PESs for these complexes are remarkably similar, which are characterized by three minima and two saddle points between them. The global minimum corresponds to a collinear Rg-Br-Cl configuration. Two local minima, correlate with an anti-linear Rg-Cl-Br geometry and a nearly T-shaped structure, can also be located on each PES. The quantum bound state calculations enable us to investigate intermolecular vibrational states and rotational energy levels of the complexes. The transition frequencies are predicted and are fitted to obtain their corresponding spectroscopic constants. In general, the periodic trends are observed for this complex family. Comparisons with available experimental data for the collinear isomer of Ar-BrCl demonstrate reliability of our theoretical predictions, and our results for the other two isomers of Ar-BrCl as well as for other members of the complex family are also anticipated to be trustable. Except for the collinear isomer of Ar-BrCl, the data presented in this paper would be beneficial to improve our knowledge for these experimentally unknown species.
Collapse
Affiliation(s)
- Song Li
- School of Physics and Optoelectronic Engineering, Yangtze University, Jingzhou, People's Republic of China.
| | - Rui Zheng
- School of Mathematics and Information Science, North China University of Water Resources and Electric Power, Zhengzhou, People's Republic of China
| | - Shan-Jun Chen
- School of Physics and Optoelectronic Engineering, Yangtze University, Jingzhou, People's Republic of China
| | - Yan Chen
- School of Physics and Optoelectronic Engineering, Yangtze University, Jingzhou, People's Republic of China
| | - Peng Chen
- School of Physics and Optoelectronic Engineering, Yangtze University, Jingzhou, People's Republic of China
| |
Collapse
|
6
|
Theoretical investigation of potential energy surface and bound states for the van der Waals complex Ar–BrCl dimer. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.07.013] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
7
|
Valdés Á, Prosmiti R. Quantum mechanical characterization of the He4ICl weakly bound complex. J Phys Chem A 2013; 117:7217-23. [PMID: 23537513 DOI: 10.1021/jp4011263] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vibrational calculations are performed for the 12-dimensional He4ICl van der Waals complex using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface of the cluster is represented as a sum of the triatomic He-ICl ab initio parametrized terms plus the He-He interactions. The topology of the potential presents higher anisotropy compared to the one with a homonuclear dopant, and this is clearly reflected in the structure and energetics of the low-lying conformers of the system. In order to take advantage of the MCTDH method, natural potential fits are employed for the potential energy operator, and also, a mode combination scheme is introduced in order to speed up the computations. Zero-point energy, binding energies, and vibrationally averaged structures of different isomers of the He4ICl cluster are obtained. The present results predict that the (3,1,0) structure, involving three He atoms in the near T-shaped and one He atom in the linear configurations, to be the most stable one in accord with recent experimental findings. Comparisons with previous theoretical and experimental data are presented, and the stability of the high-order conformers is discussed in connection with the multiple minima (global and local) of the underlying potential surface.
Collapse
Affiliation(s)
- Álvaro Valdés
- Instituto de Física Fundamental (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain
| | | |
Collapse
|
8
|
McAllister LJ, Bruce DW, Karadakov PB. Quantum Chemical Investigation of Attractive Non-Covalent Interactions between Halomethanes and Rare Gases. J Phys Chem A 2012; 116:10621-8. [DOI: 10.1021/jp307869v] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Linda J. McAllister
- Department of Chemistry, University of York, Heslington, York, YO10 5DD, U.K
| | - Duncan W. Bruce
- Department of Chemistry, University of York, Heslington, York, YO10 5DD, U.K
| | - Peter B. Karadakov
- Department of Chemistry, University of York, Heslington, York, YO10 5DD, U.K
| |
Collapse
|
9
|
Kalemos A, Valdés Á, Prosmiti R. Theoretical investigation of the He-I2(E3Πg) ion-pair state: ab initio intermolecular potential and vibrational levels. J Chem Phys 2012; 137:034303. [PMID: 22830698 DOI: 10.1063/1.4733983] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.
Collapse
Affiliation(s)
- Apostolos Kalemos
- National and Kapodistrian University of Athens, Department of Chemistry, Laboratory of Physical Chemistry, Panepistimiopolis, Athens 15771, Greece.
| | | | | |
Collapse
|
10
|
Valdés Á, Prosmiti R, Villarreal P, Delgado-Barrio G. Theoretical investigation of the He4Br2 conformers. J Phys Chem A 2012; 116:7169-76. [PMID: 22667396 DOI: 10.1021/jp3026682] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Full dimensional quantum dynamics calculations of the three lowest isomers of the He(4)Br(2) van der Waals molecule in its ground electronic state are reported. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method and a realistic potential form that includes the sum of three body ab initio coupled-cluster single double triple [CCSD(T)] He-Br(2) interactions plus the He-He and Br-Br interactions. This potential exhibits several multiple minima, with the three lowest ones lying very close in energy, just within 2 cm(-1). Such small differences are also found in the calculated binding energies of the three most stable conformers, indicating the floppiness of the system and, thus, the need of accurate potential forms and quantum full dynamics methods to treat this kind of complexes. The 12 dimensional results reported in this work present benchmark data and, thus, can serve to evaluate approximate methods aiming to describe higher order rare gas-dihalogen (N > 4) complexes. A comparison with previous studies using different potential forms and approaches to the energetics for the He(4)Br(2) cluster is also presented.
Collapse
Affiliation(s)
- Álvaro Valdés
- Instituto de Física Fundamental (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| | | | | | | |
Collapse
|
11
|
Structure and dynamics of noble gas-halogen and noble gas ionic clusters: When theory meets experiment. Chem Phys 2012. [DOI: 10.1016/j.chemphys.2011.05.026] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
12
|
Theoretical study of spectroscopy, interaction, and dissociation of linear and T-shaped isomers of RgClF (Rg = He, Ne, and Ar) van der Waals complexes. Struct Chem 2011. [DOI: 10.1007/s11224-011-9914-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
13
|
Delgado-Tellez L, Valdés Á, Prosmiti R, Villarreal P, Delgado-Barrio G. Ab initio characterization of the Ne–I2 van der Waals complex: Intermolecular potentials and vibrational bound states. J Chem Phys 2011; 134:214304. [DOI: 10.1063/1.3596604] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Krishtal A, Vanommeslaeghe K, Olasz A, Veszprémi T, Van Alsenoy C, Geerlings P. Accurate interaction energies at density functional theory level by means of an efficient dispersion correction. J Chem Phys 2009; 130:174101. [DOI: 10.1063/1.3126248] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
García-Gutierrez L, Delgado-Tellez L, Valdés Á, Prosmiti R, Villarreal P, Delgado-Barrio G. Intermolecular Ab Initio Potential and Spectroscopy of the Ground State of HeI2 Complex Revisited. J Phys Chem A 2009; 113:5754-62. [DOI: 10.1021/jp901250u] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | | | - Álvaro Valdés
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Rita Prosmiti
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Pablo Villarreal
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | | |
Collapse
|
16
|
Becke AD, Johnson ER. A density-functional model of the dispersion interaction. J Chem Phys 2007; 123:154101. [PMID: 16252936 DOI: 10.1063/1.2065267] [Citation(s) in RCA: 903] [Impact Index Per Article: 53.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
We have recently introduced [J. Chem. Phys. 122, 154104 (2005)] a simple parameter-free model of the dispersion interaction based on the instantaneous in space, dipole moment of the exchange hole. The model generates remarkably accurate interatomic and intermolecular C6 dispersion coefficients, and geometries and binding energies of intermolecular complexes. The model involves, in its original form, occupied Hartree-Fock or Kohn-Sham orbitals. Here we present a density-functional reformulation depending only on total density, the gradient and Laplacian of the density, and the kinetic-energy density. This density-functional model performs as well as the explicitly orbital-dependent model, yet offers obvious computational advantages.
Collapse
Affiliation(s)
- Axel D Becke
- Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| | | |
Collapse
|
17
|
Diez-Pardos C, Valdés A, Prosmiti R, Villarreal P, Delgado-Barrio G. Additive intermolecular potentials from ab initio calculations: trends in Rg2–dihalogen van der Waals trimers. Theor Chem Acc 2007. [DOI: 10.1007/s00214-007-0368-5] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
18
|
Valdés A, Prosmiti R, Villarreal P, Delgado-Barrio G, Lemoine D, Lepetit B. Ab initiovibrational predissociation dynamics of He–I2(B) complex. J Chem Phys 2007; 126:244314. [PMID: 17614556 DOI: 10.1063/1.2748404] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Three-dimensional quantum mechanical calculations on the vibrational predissociation dynamics of HeI2 B state complex are performed using a potential energy surface accurately fitted to unrestricted open-shell coupled cluster ab initio data, further enabling extrapolation for large I2 bond lengths. A Lanczos iterative method with an optimized complex absorbing potential is used to determine energies and lifetimes of the vibrationally predissociating He,I2(B,v') complex for v'<or=26 of I2 vibrational excitations. The calculated predissociating state energies agree with recent experimental results within 0.5 cm(-1). This excellent agreement is remarkable since no adjustment was made with respect to the experiments. The present ab initio approach, however, shows its limitations in the fact that the computed lifetimes that are highly sensitive to subtle details of the potential energy surfaces such as anisotropy, are a factor of 1.5 larger than the available experimental data.
Collapse
Affiliation(s)
- Alvaro Valdés
- Instituto de Matemáticas y Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | | | | | | | |
Collapse
|
19
|
Sato T, Tsuneda T, Hirao K. Long-range corrected density functional study on weakly bound systems: Balanced descriptions of various types of molecular interactions. J Chem Phys 2007; 126:234114. [PMID: 17600411 DOI: 10.1063/1.2747243] [Citation(s) in RCA: 126] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The long-range correction scheme for the density functional theory, combined with a van der Waals functional, is examined for its applicability to a wide variety of weakly bonded complexes including dispersion, dipole-induced dipole, dipole-dipole, and hydrogen-bonded complexes. The present method is able to describe geometries and binding energies of all these complexes accurately. Explicit inclusion of the long-range exchange and dispersion interactions is found to be important for the balanced description of various kinds of weak interactions. The present method is a promising alternative for high-level ab initio methods in calculations of large and complex systems, because it gives equally correct descriptions for various types of molecular interactions with much less computational cost.
Collapse
Affiliation(s)
- Takeshi Sato
- Department of Applied Chemistry, School of Engineering, University of Tokyo, Tokyo 113-8656, Japan.
| | | | | |
Collapse
|
20
|
Johnson ER, Becke AD. A post-Hartree-Fock model of intermolecular interactions. J Chem Phys 2007; 123:24101. [PMID: 16050735 DOI: 10.1063/1.1949201] [Citation(s) in RCA: 613] [Impact Index Per Article: 36.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Intermolecular interactions are of great importance in chemistry but are difficult to model accurately with computational methods. In particular, Hartree-Fock and standard density-functional approximations do not include the physics necessary to properly describe dispersion. These methods are sometimes corrected to account for dispersion by adding a pairwise C6R6 term, with C6 dispersion coefficients dependent on the atoms involved. We present a post-Hartree-Fock model in which C6 coefficients are generated by the instantaneous dipole moment of the exchange hole. This model relies on occupied orbitals only, and involves only one, universal, empirical parameter to limit the dispersion energy at small interatomic separations. The model is extensively tested on isotropic C6 coefficients of 178 intermolecular pairs. It is also applied to the calculation of the geometries and binding energies of 20 intermolecular complexes involving dispersion, dipole-induced dipole, dipole-dipole, and hydrogen-bonding interactions, with remarkably good results.
Collapse
Affiliation(s)
- Erin R Johnson
- Department of Chemistry, Queen's University Kingston, Ontario, Canada, K7L 3N6
| | | |
Collapse
|
21
|
Valdés A, Prosmiti R, Villarreal P, Delgado-Barrio G. A theoretical study of He2ICl van der Waals cluster. J Chem Phys 2006; 125:014313. [PMID: 16863304 DOI: 10.1063/1.2208614] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The structure, energetics, and dynamics of He2ICl complex in its ground state are studied by means of ab initio electronic structure and quantum-mechanical calculations. Interaction energies for selected He2ICl configurations are calculated at the coupled-cluster [CCSD(T)] level of theory using a large-core pseudopotential for the I atom and the aug-cc-pVTZ and aug-cc-pV5Z basis sets for the Cl and He atoms, respectively. The surface is characterized around its lower five minima and the minimum energy pathways through them. The global minimum of the potential corresponds to a "police-nightstick (1)" configuration, the second one to a linear, the next one to tetrahedral configuration, and the following two to "bifork" and "police-nightstick (2)" structures, with well depths of -99.12, -97.42, -88.32, -85.84, and -78.54 cm(-1), respectively. An analytical form based on the sum of the three-body parametrized HeICl interactions plus the He-He interaction is found to represent very well the tetra-atomic CSSD(T) results. The present potential expression is employed to perform variational five-dimensional quantum-mechanical calculations to study the vibrational bound states of the van der Waals He2ICl complex. Results for total angular momentum J = 0 provide the binding energy D0 and the corresponding vibrationally averaged structure for different isomers of the cluster. Comparison of these results with recent experimental observations further justifies the potential used in this work.
Collapse
Affiliation(s)
- Alvaro Valdés
- Instituto de Matemáticas y Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | | | |
Collapse
|
22
|
Johnson ER, Becke AD. A post-Hartree-Fock model of intermolecular interactions: Inclusion of higher-order corrections. J Chem Phys 2006; 124:174104. [PMID: 16689564 DOI: 10.1063/1.2190220] [Citation(s) in RCA: 709] [Impact Index Per Article: 39.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have previously demonstrated that the dipole moment of the exchange hole can be used to derive intermolecular C(6) dispersion coefficients [J. Chem. Phys. 122, 154104 (2005)]. This was subsequently the basis for a novel post-Hartree-Fock model of intermolecular interactions [J. Chem. Phys. 123, 024101 (2005)]. In the present work, the model is extended to include higher-order dispersion coefficients C(8) and C(10). The extended model performs very well for prediction of intermonomer separations and binding energies of 45 van der Waals complexes. In particular, it performs twice as well as basis-set extrapolated MP2 theory for dispersion-bound complexes, with minimal computational cost.
Collapse
Affiliation(s)
- Erin R Johnson
- Department of Chemistry, Queen's University, Kingston, Ontario, Canada
| | | |
Collapse
|
23
|
|
24
|
Valdés A, Prosmiti R, Villarreal P, Delgado-Barrio G. Ab initio calculations, potential representation and vibrational dynamics of He2Br2 van der Waals complex. J Chem Phys 2005; 122:44305. [PMID: 15740246 DOI: 10.1063/1.1833352] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An intermolecular potential energy surface for He(2)Br(2) complex in the ground state is calculated at the levels of fourth-order (MP4) Moller-Plesset and coupled-cluster [CCSD(T)] approximations, using large-core pseudopotential for Br atoms and the aug-cc-pV5Z basis set for He. The surface is characterized by three minima and the minimum energy pathways through them. The global minimum corresponds to a linear He-Br(2)-He configuration, while the two other ones to "police-nightstick" and tetrahedral structures. The corresponding well depths are -90.39/-89.18, -81.23/-80.78 and -74.40/-74.02 cm(-1), respectively, at MP4/CCSD(T) levels of theory. It is found that results obtained by summing three-body parametrized HeBr(2) interactions and the He-He interaction are in very good accord with the corresponding MP4/CSSD(T) configuration energies of the He(2)Br(2). Variational calculations using a sum of three-body interactions are presented to study the bound states of the vdW He(2)Br(2) complex. The binding energy D(0) and the corresponding vibrationally averaged structure are determined for different isomers of the cluster and their comparison with the available experimental data is discussed.
Collapse
Affiliation(s)
- Alvaro Valdés
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | | | | | | |
Collapse
|
25
|
Valdés * Á, Prosmiti R, Villarreal P, Delgado-Barrio G. HeBr2complex: ground-state potential and vibrational dynamics fromab initiocalculations. Mol Phys 2004. [DOI: 10.1080/00268970412331290634] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
26
|
Prosmiti R, Valdés Á, Villarreal P, Delgado-Barrio G. Three-Dimensional ab Initio Potential and Ground State Dynamics of the HeI2 Complex. J Phys Chem A 2004. [DOI: 10.1021/jp048787t] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rita Prosmiti
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Álvaro Valdés
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Pablo Villarreal
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Gerardo Delgado-Barrio
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
27
|
Prosmiti R, López-López S, García-Vela A. Potential energy surface and rovibrational states of the ground Ar–HI complex. J Chem Phys 2004; 120:6471-7. [PMID: 15267536 DOI: 10.1063/1.1665467] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A potential energy surface for the ground electronic state of the Ar-HI van der Waals complex is calculated at the coupled-cluster with single and double excitations and a noniterative perturbation treatment of triple excitations [CCSD(T)] level of theory. Calculations are performed using for the iodine atom a correlation consistent triple-zeta valence basis set in conjunction with large-core Stuttgart-Dresden-Bonn relativistic pseudopotential, whereas specific augmented correlation consistent basis sets are employed for the H and Ar atoms supplemented with an additional set of bond functions. In agreement with previous studies, the equilibrium structure is found to be linear Ar-I-H, with a well depth of 205.38 cm(-1). Another two secondary minima are also predicted at a linear and bent Ar-H-I configurations with well depths of 153.57 and 151.57 cm(-1), respectively. The parametrized CCSD(T) potential is used to calculate rovibrational bound states of Ar-HI/Ar-DI complexes, and the vibrationally averaged structures of the different isomers are determined. Spectroscopic constants are also computed from the CCSD(T) surface and their comparison with available experimental data demonstrates the quality of the present surface in the corresponding configuration regions.
Collapse
Affiliation(s)
- Rita Prosmiti
- Instituto de Matemáticas y Física Fundamental, C.S.I.C., Serrano 123, 28006 Madrid, Spain.
| | | | | |
Collapse
|
28
|
Prosmiti R, Villarreal P, Delgado-Barrio G. Structure and bonding of ArClF: Intermolecular potential surface. Isr J Chem 2003. [DOI: 10.1560/kkrp-xqkm-hfe0-9tpa] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|