Density functional theory (DFT) and atoms in molecules theory (AIM) were used to study the characteristic of the noncovalent interactions in complexes formed between Lewis bases (NH(3), H(2)O, and H(2)S) and Lewis acids (ClF, BrF, IF, BrCl, ICl, and IBr). In order to compare halogen and hydrogen bonds interactions, this study included hydrogen complexes formed by some Lewis bases and HF, HCl, and HBr Lewis acids. Ab initio, wave functions were generated at B3LYP/6-311++G(d,p) level with optimized structures at the same level. Criteria based on a topological analysis of the electron density were used in order to characterize the nature of halogen interactions in Lewis complexes. The main purpose of the present work is to provide an answer to the following questions: (a) why can electronegative atoms such as halogens act as bridges between two other electronegative atoms? Can a study based on the electron charge density answer this question? Considering this, we had performed a profound study of halogen complexes in the framework of the AIM theory. A good correlation between the density at the intermolecular bond critical point and the energy interaction was found. We had also explored the concentration and depletion of the charge density, displayed by the Laplacian topology, in the interaction zone and in the X-Y halogen donor bond. From the atomic properties, it was generally observed that the two halogen atoms gain electron population in response to its own intrinsic nature. Because of this fact, both atoms are energetically stabilized.

In this work we investigate the nature of the Cl···N interactions in complexes formed between substituted ammonium [NHn(X3-n) (with n = 0, 1, 2, 3 and X = -CH3, -F] as Lewis bases and F-Cl molecule as Lewis acid. They have been chosen as a study case due to the wide range of variation of their binding energies, BEs. Møller-Plesset [MP2/6-311++G(2d,2p)] calculations show that the BEs for this set of complexes lie in the range from 1.27 kcal/mol (in F-Cl···NF3) to 27.62 kcal/mol [in F-Cl···N(CH3)3]. The intermolecular distribution of the electronic charge density and their L(r) = -¼∇(2)ρ(r) function have been investigated within the framework of the atoms in molecules (AIM) theory. The intermolecular interaction energy decomposition has also been analyzed using the reduced variational space (RVS) method. The topological analysis of the L(r) function reveals that the local topological properties measured at the (3,+1) critical point [in L(r) topology] are good descriptors of the strength of the halogen bonding interactions. The results obtained from energy decomposition analysis indicate that electrostatic interactions play a key role in these halogen bonding interactions. These results allow us to establish that, when the halogen atom is bonded to a group with high electron-withdrawing capacity, the electrostatic interaction between the electron cloud of the Lewis base and the halogen atom unprotected nucleus of the Lewis acid produces the formation and determines the geometry of the halogen bonded complexes. In addition, a good linear relationship has been established between: the natural logarithm of the BEs and the electrostatic interaction energy between electron charge distribution of N atom and nucleus of Cl atom, denoted as V e-n(N,Cl) within the AIM theory.

We report here an exhaustive and complete conformational study on the conformational potential energy hypersurface (PEHS) of dopamine (DA) interacting with the dopamine D2 receptor (D2-DR). A reduced 3D model for the binding pocket of the human D2-DR was constructed on the basis of the theoretical model structure of bacteriorhodopsin. In our reduced model system, only 13 amino acids were included to perform the quantum mechanics calculations. To obtain the different complexes of DA/D2-DR, we combined semiempirical (PM6), DFT (B3LYP/6-31G(d)), and QTAIM calculations. The molecular flexibility of DA interacting with the D2-DR was evaluated from potential energy surfaces and potential energy curves. A comparative study between the molecular flexibility of DA in the gas phase and at D2-DR was carried out. In addition, several molecular dynamics simulations were carried out to evaluate the molecular flexibility of the different complexes obtained. Our results allow us to postulate the complexes of type A as the "biologically relevant conformations" of DA. In addition, the theoretical calculations reported here suggested that a mechanistic stepwise process takes place for DA in which the protonated nitrogen group (in any conformation) acts as the anchoring portion, and this process is followed by a rapid rearrangement of the conformation allowing the interaction of the catecholic OH groups.

No abstract

The nature of F-BrX-R interactions (with X = F, Cl, Br, I and R = -H, -F) has been investigated through theoretical calculation of molecular potential electrostatic (MEP), molecular polarizability, atoms in molecules (AIM) analysis and energetic decomposition analysis (EDA). A detailed analysis of the MEPs reveals that considering only the static electrostatic interactions is not sufficient to explain the nature of these interactions. The molecular polarizabilities of X-R molecules suggest that the deformation capacity of the electronic cloud of the lone pairs of the X atom plays an important role in the stability of these complexes. The topological analysis of the L(r) = -¼∇(2)ρ(r) function and the detailed analysis of the atomic quadrupole moments reveal that the BrX interactions are electrostatic in nature. The electron acceptor Br atom causes a polarization of the electronic cloud (electronic induction) on the valence shell of the X atom. Finally, the electrostatic forces and charge transfer play an important role not only in the stabilization of the complex, but also in the determination of the molecular geometry of equilibrium. The dispersive and polarization forces do not influence the equilibrium molecular geometry.

In this work, the intermolecular distribution of the electronic charge density in the aromatic hydrogen/halogen bonds is studied within the framework of the atoms in molecules (AIM) theory and the molecular electrostatic potentials (MEP) analysis. The study is carried out in nine complexes formed between benzene and simple lineal molecules, where hydrogen, fluorine and chlorine atoms act as bridge atoms. All the results are obtained at MP2 level theory using cc-pVTZ basis set. Attention is focused on topological features observed at the intermolecular region such as bond, ring and cage critical points of the electron density, as well as the bond path, the gradient of the density maps, molecular graphs and interatomic surfaces. The strength of the interaction increases in the following order: F[Symbol: see text]pi < Cl[Symbol: see text]pi < H[Symbol: see text]pi. Our results show that the fluorine atom has the capability to interact with the pi-cloud to form an aromatic halogen bond, as long as the donor group is highly electron withdrawing. The Laplacian topology allows us to state that the halogen atoms can act as nucleophiles as well as electrophiles, showing clearly their dual character.

In this work we studied the host–guest interactions between confined molecules and zeolites and their relationship with the energies involved in the reaction of methylation of benzene by methanol in H-ZSM-5 and H-Beta zeolites employing density functional theory (DFT) methods and the quantum theory of atoms in molecules. Results show that the strength of the interactions related to adsorption and coadsorption processes is higher in the catalyst with the larger cavity; however, the confinement effects are higher in the smaller zeolite, explaining, from an electronic viewpoint, the reason why the stabilization energy is higher in H-ZSM-5 than in H-Beta. The confinement effects of the catalyst on the confined species for methanol adsorption, benzene coadsorption, and the formed intermediates dominate this stabilization. For the transition state (TS), the stability of the TS is achieved due to the stabilizing effect of the surrounding zeolite framework on the formed carbocationic species (CH3 +) which is higher in H-ZSM-5 than in H-Beta. In both TSs the methyl cation is multicoordinated forming the following H2O···CH3 +···CB concerted bonds. It is demonstrated that, through the electron density analysis, the criteria can be defined to discriminate between interactions related to the confinement effects and the reaction itself (adsorption, coadsorption, and bond-breaking and bond-forming processes) and, thus, to discriminate the relative contributions of the degree of confinement to the reaction energies for two zeolite catalysts with different topologies.

In the present work, the distribution of the electronic charge density in the ethene protonation reaction by a zeolite acid site is studied within the framework of the density functional theory and the atoms in molecules (AIM) theory. The key electronic effects such as topological distribution of the charge density involved in the reaction are presented and discussed. The results are obtained at B3LYP/6-31G(**) level theory. Attention is focused on topological parameters such as electron density, its Laplacian, kinetic energy density, potential energy density, and electronic energy density at the bond critical points (BCP) in all bonds involved in the interaction zone, in the reactants, pi-complex, transition state, and alkoxy product. In addition, the topological atomic properties are determined on the selected atoms in the course of the reaction (average electron population, N(Omega), atomic net charge, q(Omega), atomic energy, E(Omega), atomic volume, v(Omega), and first moment of the atomic charge distribution, M(Omega)) and their changes are analyzed exhaustively. The topological study clearly shows that the ethene interaction with the acid site of the zeolite cluster, T5-OH, in the ethene adsorbed, is dominated by a strong O-H...pi interaction with some degree of covalence. AIM analysis based on DFT calculation for the transition state (TS) shows that the hydrogen atom from the acid site in the zeolitic fragment is connected to the carbon atom by a covalent bond with some contribution of electrostatic interaction and to the oxygen atom by closed shell interaction with some contribution of covalent character. The C-O bond formed in the alkoxy product can be defined as a weaker shared interaction. Our results show that in the transition state, the dominant interactions are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the concentration and accumulation of the charge density along the bond path between the nuclei linked increases.