On the nature of homo- and hetero-dinuclear metal–metal quadruple bonds — Analysis of the bonding situation and benchmarking DFT against wave function methods

2010 ◽  
Vol 88 (11) ◽  
pp. 1079-1093 ◽  
Author(s):  
Nozomi Takagi ◽  
Andreas Krapp ◽  
Gernot Frenking

Homo- and hetero-dimetallic (d–d)8 analogues of the formally quadruply bonded [Re2Cl8]2– system with the general formula [MM′Cl8]x (M, M′ = Tc, Re, Ru, Os, Rh, Ir and x = –2, –1, 0, +1, +2) have been calculated with the density functional theory (DFT) functionals SVWN, BLYP, BP86, PBE, OLYP, OPBE, HCTH, B3LYP, O3LYP, X3LYP, BH&HLYP, TPSS, VSXC, TPPSh, and ab initio methods (CASPT2, CCSD(T)) using basis sets of triple-ζ quality. The performance of the functionals for the description of the metal–metal bond distance and the bond dissociation energy as well as the singlet–triplet gap was evaluated with respect to ab initio data at the CASPT2 level. Generally, the generalized gradient approximation (GGA) functionals, BLYP, BP86, and PBE, show good performance in the description of the metal–metal bond distance and for the dissociation energy. Hybrid functionals are not to be used for compounds of the type discussed here as they lead to increasingly too short and too weak bonds with the amount of exact exchange included. All functionals underestimate the singlet–triplet gap, with the GGA functionals BLYP, BP86, PBE being the closest to the CASPT2 values. The bonding situations of the [MM′Cl8]x compounds were analyzed at the DFT level (BP86) using the natural bond orbital (NBO) method and the energy decomposition analysis. The M–M bond in homodimetallic compounds, [MMCl8]x, becomes weaker from group 7 to group 8 to group 9 metals and the bond is weaker for 4d metal systems than for 5d transition metal compounds. The M–M bonds have approximately 50% covalent and 50% electrostatic character and the covalent contribution is dominated by the π orbitals, whereas the δ orbitals do not contribute significantly to the covalent bonding. Heterodimetallic systems, [MM′Cl8]x, have significantly stronger metal–metal bonds than the homodimetallic compounds. This comes from weaker Pauli repulsion and stronger electrostatic attraction. The most stable heterodimetallic bonds are observed for 5d–5d metal pairs.

2020 ◽  
Vol 2020 ◽  
pp. 1-9 ◽  
Author(s):  
P. O. Jomo ◽  
C. O. Otieno ◽  
P. W. O. Nyawere

We report the results of pressure-induced semiconductor-metal phase transition of the semiconducting chalcogenide compound KPSe6 under high pressure using the ab initio methods. The ground-state energy calculations were performed within density functional theory and the generalized gradient approximation using the pseudopotential method with plane-wave basis sets. The projector augmented-wave (PAW) pseudopotentials were used in our calculation. The optimized lattice parameters were found from total energy calculations as 13 Bohr, 1.6 Bohr, and 1.8 Bohr for cell dimensions one, two, and three, respectively, which are in good agreement with experimental calculations. At zero pressure, the material portrayed a semiconducting property with a direct bandgap of ≈1.7 eV. As we subjected the material to pressure, the band gap was observed to reduce until it disappeared. The phase transition from the semiconductor to metal was found to occur at ∼45 GPa, implying that the material underwent metallization as pressure was increased further.


Author(s):  
Yuming Zhao ◽  
Cody Marcus King-Poole

The noncovalent interactions between a redox-active molecule, phenyl-substituted dithiafulvene (Ph-DTF), and ten commonly encountered nitroaromatic compounds (NACs) were systematically investigated by means of density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations. Our modeling studies examined their 1:1 complexes in terms of equilibrium geometries, frontier molecular orbitals (FMOs), nature of noncovalent forces, intermolecular charge transfer (ICT), interaction energies and related energy decomposition analysis. The computational results indicate that Ph-DTF can form thermodynamically stable supramolecular complexes with trinitro-substituted benzenes (e.g., 2,4,6-trinisuchtrotoluene and picric acid), but its interactions with mono- and dinitrobenzenes do not exhibit such stability. The selective binding properties are further corroborated by AIMD simulations. Overall, this computational work establishes a comprehensive understanding of the nature of noncovalent interactions of Ph-DTF with various NACs, and the results can be used as theoretical guidance for the rational design of selective receptors and/or chemosensors for certain NACs that are of great concern in current industrial applications and environmental control.


2017 ◽  
Vol 114 (48) ◽  
pp. 12649-12656 ◽  
Author(s):  
Daniel S. Levine ◽  
Martin Head-Gordon

An energy decomposition analysis (EDA) for single chemical bonds is presented within the framework of Kohn–Sham density functional theory based on spin projection equations that are exact within wave function theory. Chemical bond energies can then be understood in terms of stabilization caused by spin-coupling augmented by dispersion, polarization, and charge transfer in competition with destabilizing Pauli repulsions. The EDA reveals distinguishing features of chemical bonds ranging across nonpolar, polar, ionic, and charge-shift bonds. The effect of electron correlation is assessed by comparison with Hartree–Fock results. Substituent effects are illustrated by comparing the C–C bond in ethane against that in bis(diamantane), and dispersion stabilization in the latter is quantified. Finally, three metal–metal bonds in experimentally characterized compounds are examined: a MgI–MgI dimer, the ZnI–ZnI bond in dizincocene, and the Mn–Mn bond in dimanganese decacarbonyl.


2014 ◽  
Vol 69 (7) ◽  
pp. 385-395 ◽  
Author(s):  
Nicole Holzmann ◽  
Gernot Frenking

Quantum chemical calculations using density functional theory at the BP86 level in conjunction with triple-zeta polarized basis sets have been carried out for the title compounds. The nature of the bonding between the diatomic fragment and the NHC ligands is investigated with an energy decomposition analysis. The chemical bonds in the [(NHCMe)2(E2)] complexes can be discussed in terms of donor-acceptor interactions which consist of two NHCMe→E2←NHCMe donor components and two weaker components of the NHCMe←E2→NHCMe π backdonation. The out-of-phase (+)=(-) contribution of the s donation is always stronger than the in-phase (+)=(+) contribution. The electronic reference state of N2 in the dinitrogen complex [(NHCMe)2(N2)] is the highly excited 11Γg state which explains the anti-periplanar arrangement of the ligands. The gauche arrangement of the ligands in the heavier homologues [(NHCMe)2(E2)] (E = P-Bi) may be discussed using either the excited 11Γg state or the X1Σg+ ground state of E2 as reference states for the donor-acceptor bonds. The EDA-NOCV calculations suggest that the latter bonding model is better suited for the complexes where E = As-Bi while the phosphorus complex is a borderline case.


2020 ◽  
Author(s):  
Luis Vasquez ◽  
Agnieszka Dybala-Defratyka

<p></p><p>Very often in order to understand physical and chemical processes taking place among several phases fractionation of naturally abundant isotopes is monitored. Its measurement can be accompanied by theoretical determination to provide a more insightful interpretation of observed phenomena. Predictions are challenging due to the complexity of the effects involved in fractionation such as solvent effects and non-covalent interactions governing the behavior of the system which results in the necessity of using large models of those systems. This is sometimes a bottleneck and limits the theoretical description to only a few methods.<br> In this work vapour pressure isotope effects on evaporation from various organic solvents (ethanol, bromobenzene, dibromomethane, and trichloromethane) in the pure phase are estimated by combining force field or self-consistent charge density-functional tight-binding (SCC-DFTB) atomistic simulations with path integral principle. Furthermore, the recently developed Suzuki-Chin path integral is tested. In general, isotope effects are predicted qualitatively for most of the cases, however, the distinction between position-specific isotope effects observed for ethanol was only reproduced by SCC-DFTB, which indicates the importance of using non-harmonic bond approximations.<br> Energy decomposition analysis performed using the symmetry-adapted perturbation theory (SAPT) revealed sometimes quite substantial differences in interaction energy depending on whether the studied system was treated classically or quantum mechanically. Those observed differences might be the source of different magnitudes of isotope effects predicted using these two different levels of theory which is of special importance for the systems governed by non-covalent interactions.</p><br><p></p>


2019 ◽  
Author(s):  
Brandon B. Bizzarro ◽  
Colin K. Egan ◽  
Francesco Paesani

<div> <div> <div> <p>Interaction energies of halide-water dimers, X<sup>-</sup>(H<sub>2</sub>O), and trimers, X<sup>-</sup>(H<sub>2</sub>O)<sub>2</sub>, with X = F, Cl, Br, and I, are investigated using various many-body models and exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. Analysis of the results obtained with the many-body models demonstrates the need to capture important short-range interactions in the regime of large inter-molecular orbital overlap, such as charge transfer and charge penetration. Failure to reproduce these effects can lead to large deviations relative to reference data calculated at the coupled cluster level of theory. Decompositions of interaction energies carried out with the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) method demonstrate that permanent and inductive electrostatic energies are accurately reproduced by all classes of XC functionals (from generalized gradient corrected (GGA) to hybrid and range-separated functionals), while significant variance is found for charge transfer energies predicted by different XC functionals. Since GGA and hybrid XC functionals predict the most and least attractive charge transfer energies, respectively, the large variance is likely due to the delocalization error. In this scenario, the hybrid XC functionals are then expected to provide the most accurate charge transfer energies. The sum of Pauli repulsion and dispersion energies are the most varied among the XC functionals, but it is found that a correspondence between the interaction energy and the ALMO EDA total frozen energy may be used to determine accurate estimates for these contributions. </p> </div> </div> </div>


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Samaneh Pasban ◽  
Heidar Raissi

AbstractHexakis (m-phenylene ethynylene) (m-PE) macrocycles, with aromatic backbones and multiple hydrogen-bonding side chains, had a very high propensity to self-assemble via H-bond and π–π stacking interactions to form nanotubular structures with defined inner pores. Such stacking of rigid macrocycles is leading to novel applications that enable the researchers to explored mass transport in the sub-nanometer scale. Herein, we performed density functional theory (DFT) calculations to examine the drug delivery performance of the hexakis dimer as a novel carrier for doxorubicin (DOX) agent in the chloroform and water solvents. Based on the DFT results, it is found that the adsorption of DOX on the carrier surface is typically physisorption with the adsorption strength values of − 115.14 and − 83.37 kJ/mol in outside and inside complexes, respectively, and so that the essence of the drug remains intact. The negative values of the binding energies for all complexes indicate the stability of the drug molecule inside and outside the carrier's cavities. The energy decomposition analysis (EDA) has also been performed and shown that the dispersion interaction has an essential role in stabilizing the drug-hexakis dimer complexes. To further explore the electronic properties of dox, the partial density of states (PDOS and TDOS) are calculated. The atom in molecules (AIM) and Becke surface (BS) methods are also analyzed to provide an inside view of the nature and strength of the H-bonding interactions in complexes. The obtained results indicate that in all studied complexes, H-bond formation is the driving force in the stabilization of these structures, and also chloroform solvent is more favorable than the water solution. Overall, our findings offer insightful information on the efficient utilization of hexakis dimer as drug delivery systems to deliver anti-cancer drugs.


2015 ◽  
Vol 11 ◽  
pp. 2727-2736 ◽  
Author(s):  
Diego M Andrada ◽  
Nicole Holzmann ◽  
Thomas Hamadi ◽  
Gernot Frenking

Fifteen cyclic and acylic carbenes have been calculated with density functional theory at the BP86/def2-TZVPP level. The strength of the internal X→p(π) π-donation of heteroatoms and carbon which are bonded to the C(II) atom is estimated with the help of NBO calculations and with an energy decomposition analysis. The investigated molecules include N-heterocyclic carbenes (NHCs), the cyclic alkyl(amino)carbene (cAAC), mesoionic carbenes and ylide-stabilized carbenes. The bonding analysis suggests that the carbene centre in cAAC and in diamidocarbene have the weakest X→p(π) π-donation while mesoionic carbenes possess the strongest π-donation.


2003 ◽  
Vol 793 ◽  
Author(s):  
Daniel I Bilc ◽  
S.D. Mahanti ◽  
M.G. Kanatzidis

ABSTRACTComplex quaternary chalcogenides (AgSb)xPbn-2xTen (0<x<n/2) are thought to be narrow band-gap semiconductors which are very good candidates for room and high temperature thermoelectric applications. These systems form in the rock-salt structure similar to the well known two component system PbTe (x=0). In these systems Ag and Sb occupy Pb sites randomly although there is some evidence of short-range order. To gain insights into the electronic structure of these compounds, we have performed electronic structure calculations in AgSbTe2 (x=n/2). These calculations were carried out within ab initio density functional theory (DFT) using full potential linearized augmented plane wave (LAPW) method. The generalized gradient approximation (GGA) was used to treat the exchange and correlation potential. Spinorbit interaction (SOI) was incorporated using a second variational procedure. Since it is difficult to treat disorder in ab initio calculations, we have used several ordered structures for AgSbTe2. All these structures show semimetallic behavior with a pseudogap near the Fermi energy. Te and Sb p orbitals, which are close in energy, hybridize rather strongly indicating a covalent interaction between Te and Sb atoms.


Sign in / Sign up

Export Citation Format

Share Document