Hirshfeld Charge Research Articles

Why Are Some Pnictogen(III) Pincer Complexes Planar and Others Pyramidal?

Geometrically-constrained pnictogen pincer complexes have emerged in recent years as platforms for unique stoichiometric and catalytic chemical transformations. These complexes feature dynamic conformations ranging from fully planar at the pnictogen centre to distorted-pyramidal geometries, as well as variation between phases. Although the valued reactivity of pnictogen pincer complexes is ascribed to their geometries, there is no unified model to explain the observed conformational outcomes across different ligands and pnictogen centres. Here we propose such a model through computational analysis of more than 1300 structures across 64 complexes (16 ligands and 4 heavy pnictogens), explaining the experimental observations and making new predictions. By looking at signatures of bond stability (bond lengths, Wiberg bond indices) and delocalization (NPA charges, Hirshfeld charges), our framework posits a pnictogen-based σ-bonding effect that favours pyramidalization and exists in competition with a ligand-based π-bonding effect that favours planarity. Variations in structure as a function of pnictogen identity, ligand tethering, electronics, and aromaticity can be reconciled with reference to a balance between these two opposing forces. Careful consideration of the σ/π-bonding effects may aid in the rational design of future pnictogen pincer complexes with predictable geometries and reactivities.

Information-theoretic quantities as effective descriptors of electrophilicity and nucleophilicity in density functional theory.

Electrophilicity and nucleophilicity are two vastly important chemical concepts gauging the capability of atoms in molecules to accept and donate the maximal number of electrons. In our earlier studies, we proposed to simultaneously quantify them using the Kullback-Leibler divergence from the information-theoretic approach in density functional theory. However, several issues with this scheme remain to be clarified such as its general validity, predictability, and relationship with other information-theoretic quantities. In this work, we revisit the matter with bigger datasets and deeper theoretical insights. Five information-theoretic quantities including Kullback-Leibler divergence, Hirshfeld charge, Ghost-Berkowitz-Parr entropy, and second and third orders of relative Onicescu information energy are found to be reliable and robust descriptors of electrophilicity and nucleophilicity propensities. Employing these five descriptors, we design a list of new compounds and predict their electrophilicity and nucleophilicity scales. This work should markedly improve our confidence and capability in applying information-theoretic quantities to evaluate electrophilicity and nucleophilicity propensities and henceforth pave the route for more applications of these quantities from information-theoretic approach in density functional theory in the future. All structures were fully optimized at the M06-2X/6-311 + G(d) level of DFT functional using the Gaussian 16 package (version C01) with integration grids and tight self-consistent-field convergence. The solvent effect was taken into account by using the implicit solvent model (CPCM) in the CH2Cl2 solvent, and all 3D contour surfaces of Fukui function, local temperature, and ITA (information-theoretic approach) quantities were generated by GaussView. The Multiwfn 3.8 program was used to calculate the ITA indexes and atomic charges.

Co incorporated pentagraphene as an efficient single-atom catalyst for reduction of CO2: First-principles investigations

Hydrogenation of carbon dioxide (CO2) to formic acid is explored through density functional theory calculations over Co decorated penta-graphene (Co@PG) as single atom catalyst (SAC). Co is successfully stabilized over pentagraphene due to strong overlapping of Co-3d and C-2p states of pentagonal ring of pentagraphene near fermi level. The higher diffusion barrier for Co over pentagraphene evince that the clustering of Co atom is restricted. This suggests that Co@PG can be used as a catalyst for reducing CO2. Moreover, the Hirshfeld charge density, molecular electrostatic potential and spin density analysis indicates that the Co site is the highly active site for CO2 hydrogenation reaction. Both CO2 and H2 show strong interaction with Co@PG. In co-adsorbed state, the interactions are even larger due to greater mutual polarization of reactants and localization of positive charge on Co with adsorption energy of -2.57 eV. The hydrogenation can proceed at ambient temperature due to a small energy barrier of 0.29 eV, the competitive reaction i.e. CO hydrogenation with CO2 hydrogenations almost impossible over the catalyst surface due to larger energy barriers i.e. 1.98 eV. Thus, this work shows that Co@PG can be used as an efficient and promising catalyst for CO2 hydrogenation at room temperature.

The surface defect of Fe2Ge (001) induces a synergistic spin effect and regulates the surface magnetism

The defect structure leads to the fracture and recombination of chemical bonds, lattice distortion, electron localization, etc., which has a significant improvement on the physical properties of materials, and give them new functions. The unique structure of Fe2Ge forms the surface structure of Ge end face and Fe end face in the direction of (001), but the surface is easy to form defect structure, and affecting its properties. It is found that the spin direction of Ge atom is opposite to the spin direction of Fe atom in perfect Fe2Ge (001) surface structure. When there is VGe defect, the magnetic moments of the perfect Fe2Ge (001) surface and Fe atoms are enhanced. When there are VFe, VFe-Ge and VGe-Fe defects, the magnetic moments of defects Fe2Ge (001) surface and Fe atoms are weakened. The reason is that the spin state electrons of Fe atoms induce Ge atoms to produce spin electrons, which forms a synergistic spin effect. The Hirshfeld Charge distribution and charge spin of Fe and Ge atoms are affected by the local potential field, Ge atoms are excited with spin direction of Fe atoms with VFe and VFe-Ge, and Ge atoms induce magnetic moments with spin direction of Fe atoms with VGe and VGe-Fe. The effects mechanism is that the five orbitals dzz, dxy, dzy, dzx and dxx-yy of Fe atom are greatly affected by the structure of the defect under the action of crystal field, while the s, px, py and pz orbitals are less affected. In VGe and VFe-Ge defects, the spin splitting of the d orbitals is larger to form a larger magnetic moment. In VFe and VGe-Fe defects, the d orbital spin splitting is smaller, and resulting in a smaller magnetic moment. The s, px, py and pz of Ge atoms are greatly affected, the spin splitting of pz is larger than that of px and py orbitals, and indicating that the spin of Ge atoms mainly comes from pz orbitals.

Computational exploring the potential of pure and Ag-decorated WTe2 for detecting volatile organic compounds (VOCs)

Recent research used a method called density functional theory (DFT) to examine how a thin layer of 2-dimensional Transition metal dichalcogenide (TMDC) WTe2 monolayer (ML) interacts with silver (Ag) and certain gases present in human breath. We explored how this interaction could assist in identifying particular volatile organic compounds (VOCs) such as acetone (CH3O6), ethanol (C2H5OH), methanol (CH3OH), toluene (C7H8), and formaldehyde (CH2O), Together with the gases causing interference, like (CO2) and (H2O), in exhaled breath. The WTe2 layer usually acts as a semiconductor, but when combined with Ag, it becomes metallic, enhancing its ability to detect these VOCs. They analyzed various properties such as structural, electronic, adsorption, Hirshfeld charge analysis, and sensing properties of the studied material. The findings highlight ethanol as particularly effective in interacting with the layer and remaining detectable for an extended period. In conclusion, they propose that this Ag-decorated WTe2 layer holds potential for creating sensitive sensors for breath analysis at room temperature.

Atomic-scale tribological behavior of pressure-induced MoO3(0 1 0)/Ag(1 1 1) heterogeneous sliding interface

The atomic-scale tribological properties of the MoO3(010)/Ag(111) heterogeneous sliding interface were studied by first-principles calculation method. The interface binding energy, sliding potential energy, friction force, charge density difference, partial density of states, Mulliken bond population and HiRshfeld charge were analyzed, and the essence of the friction was revealed. The results show that the sliding potential energy and friction exhibit anisotropy for the pressure-induced MoO3(010)/Ag(111) heterogeneous sliding interface. The friction force can be described by the change of the sliding potential energy, and the charge redistribution between interface atoms is the intrinsic reason for the change of tribological properties. This work provides atomic and electronic level insights for understanding the tribological characteristics between two-dimensional materials and soft metals.

Synthesis, single crystal X-ray, and DFT study of new hybrid-ligand complex [Cu(hfac)2(Me3TTF-CH=CH-Pyr)] and new mixed-valence radical ion salt (Me3TTF-CH=CH-Pyr)2(PF6)3

This paper presents a comprehensive study that combines experimental techniques and Density Functional Theory (DFT) to investigate a newly synthesized mixed-ligand complex, [Cu(hfac)2(Me3TTF-CH=CH-Pyr)] (where hfac is hexafluoroacetylacetonate), and a novel mixed-valence radical ion salt (Me3TTF-CH=CH-Pyr)2(PF6)3. Both are constructed using ligands derived from the tetrathiafulvalene (TTF) donor molecule. These materials are meticulously characterized through single-crystal X-ray diffraction analysis. The complex crystallizes in the monoclinic system, specifically in the C2/c space group, with the following lattice parameters: a = 24.583(2) Å, b = 21.613(2) Å, c = 14.2066(13) Å, β = 99.789(5)°, and a volume of V = 7438.3(12) Å3. It exhibits five-coordinate behavior, forming a slightly distorted pyramid with a square base. The paper also provides a detailed account of the crystal structures and reports spectroscopic study results. Electrochemical behavior is explored through cyclic voltammetry in CH2Cl2 as the solvent. Additionally, the reactivity of these compounds is predicted using various theoretical tools, such as Local and Global reactivity descriptors, Hirshfeld charge analyses, and Molecular Electrostatic Potential. To elucidate charge transfer processes, Time-Dependent Density Functional Theory is utilized. Various types of interactions are analyzed through the Quantum Theory of Atoms in Molecules and Non-Covalent Interaction analyses. Finally, conductivity is estimated using frontier molecular orbitals (FMO) and gap energies as crucial parameters.

Ionic liquids intercalation in titanium carbide MXenes: A first-principles investigation.

Herein, we present a density functional theory with dispersion correction (DFT-D) calculations that focus on the intercalation of ionic liquids (ILs) electrolytes into the two-dimensional (2D) Ti3C2Tx MXenes. These ILs include the cation 1-ethyl-3-methylimidazolium (Emim+), accompanied by three distinct anions: bis(trifluoromethylsulfonyl)imide (TFSA-), (fluorosulfonyl)imide (FSA-) and fluorosulfonyl(trifluoromethanesulfonyl)imide (FTFSA-). By altering the surface termination elements, we explore the intricate geometries of IL intercalation in neutral, negative, and positive pore systems. Accurate estimation of charge transfer is achieved through five population analysis models, such as Hirshfeld, Hirshfeld-I, DDEC6 (density derived electrostatic and chemical), Bader, and VDD (voronoi deformation density) charges. In this work, we recommend the DDEC6 and Hirshfeld-I charge models, as they offer moderate values and exhibit reasonable trends. The investigation, aimed at visualizing non-covalent interactions, elucidates the role of cation-MXene and anion-MXene interactions in governing the intercalation phenomenon of ionic liquids within MXenes. The magnitude of this role depends on two factors: the specific arrangement of the cation, and the nature of the anionic species involved in the process.

Coordinate-Free and Low-Order Scaling Machine Learning Model for Atomic Partial Charge Prediction for Any Size of Molecules.

The atomic partial charge is of great importance in many fields, such as chemistry and drug-target recognition. However, conventional quantum-based computing of atomic charges is relatively slow, limiting further applications of atomic charge analysis. With the help of machine learning methods, various kinds of models appear to speed up atomic charge calculations. However, there are still some concerning problems. Some models based on geometric coordinates require high-accuracy geometry optimization as a preprocess, while other models have a limitation on the size of input molecules that narrow the applications of the model. Here, we propose a machine learning atomic charge model based on a message-passing featurizer. This preprocessing featurizer can quickly extract atomic environment information from a molecule according to the connectivity inside the molecule. The resulting descriptor can be used with a neural network to quickly predict the atomic partial charge. The model is able to automatically adapt to any size of molecule while remaining efficient and achieves a root-mean-square error in the Hirshfeld charge prediction of 0.018e, with an overall time complexity of O(n2). Thus, this model could enlarge the range of applications of atomic partial charge to more fields and cases.

Evaluation of hydrogen storage capacity of two-dimensional Sc2N MXene: A DFT study

In this work, the efficiency of hydrogen storage on Sc2N MXenes monolayer is systematically investigated using first-principles calculations. After determining the structural stability of designed monolayer, adsorption sites for the storage of hydrogen atoms (H) and molecules (H2) are identified. The calculated binding energies have indicated that the hydrogen atoms have favorable adsorption at the center or H site of the Sc atoms and hydrogen molecules are found to be the most stable when adsorbed on the top/T site of Sc atoms. The binding energy assessment identified the chemisorption for hydrogen atoms and Kubas-type interactions for H2. To further ascertain the suitability of the Sc2N monolayer for hydrogen storage various analyses including binding energy, partial density of states (PDOS), electron localization function (ELF), Mullikan charge analysis and Hirshfeld charge analysis are conducted. In order to determine the thermal stability of the Sc2N monolayer, desorption temperature with respect to pressure and molecular dynamics study has also been conducted and found the monolayer thermally stable. Desorption isotherms condition with increasing pressure confirms the gravimetric storage uptake. The obtained results have revealed that the hydrogen storage capacity considered monolayer has reached up to 10.8 wt%, which is reasonably high. Hence, the obtained results of hydrogen atom and molecules adsorption have clearly established that Sc2N monolayer is as an exemplary, stable, reliable and efficient medium for hydrogen storage applications.

Quantum chemical calculations and molecular dynamics simulations to investigate the mechanism of interaction of six dye levelers with copper surface

Six dye levelers are researched by quantum chemical calculations and molecular dynamics simulations. The electrostatic potential (ESP) and Hirshfeld charge of the levelers were calculated to compare the physical factors affecting their adsorption capacity (electrostatic effects) and to analyse the possible physical adsorption sites for each leveler. The average localised ionisation energy (ALIE) and the condensed Fukui function (CFF) of the levelers were calculated to compare the chemical factors (reactivity) affecting their adsorption capacity and to analyse the possible chemisorption sites for each leveler. On this basis, the adsorption capacity of the six levelers on the surface of Cu(111) at different temperatures was investigated through molecular dynamics simulations, and the adsorption characteristics of the six levelers were summarised.

Nickel-Catalyzed Three-Component Unsymmetrical Bis-Allylation of Alkynes with Alkenes: A Density Functional Theory Study.

Density functional theory (DFT) characterizations were employed to resolve the structural and energetic aspects and product selectivities along the mechanistic reaction paths of the nickel-catalyzed three-component unsymmetrical bis-allylation of alkynes with alkenes. Our putative mechanism initiated with the in situ generation of the active catalytic species [Ni(0)L2] (L = NHC) from its precursors [Ni(COD)2, NHC·HCl] to activate the alkyne and alkene substrates to form the final skipped trienes. This proceeds via the following five sequential steps: oxidative addition (OA), β-F elimination, ring-opening complexation, C-B cleavage and reductive elimination (RE). Both the OA and RE steps (with respective free energy barriers of 24.2 and 24.8 kcal·mol-1) contribute to the observed reaction rates, with the former being the selectivity-controlling step of the entire chemical transformation. Electrophilic/nucleophilic properties of selected substrates were accurately predicted through dual descriptors (based on Hirshfeld charges), with the chemo- and regio-selectivities being reasonably predicted and explained. Further distortion/interaction and interaction region indicator (IRI) analyses for key stationary points along reaction profiles indicate that the participation of the third component olefin (allylboronate) and tBuOK additive played a crucial role in facilitating the reaction and regenerating the active catalyst, ensuring smooth formation of the skipped triene product under a favorably low dosage of the Ni(COD)2 catalyst (5 mol%).

Triboelectric Charging Properties of the Functional Groups of Common Pharmaceutical Materials Using Density Functional Theory Calculations.

Triboelectrification is a ubiquitous and poorly understood phenomenon in powder processing, particularly for pharmaceutical powders. Charged particles can adhere to vessel walls, causing sheeting; they can also cause agglomeration, threatening the stability of powder formulations, and in extreme cases electrostatic discharges, which present a serious fire and explosion hazard. Triboelectrification is highly sensitive to environmental and material conditions, which makes it very difficult to compare experimental results from different publications. In this work, density functional theory (DFT) is used to investigate the charge transfer characteristics of several functional groups of paracetamol in order to better understand the mechanisms of charging at the nanoscale and the influence of the environmental and material properties on charge transfer. This is achieved by studying the structure and electronic properties at the molecule-substrate interface. Using this molecule-substrate approach, the charging contributions of individual functional groups are explored by examining the Hirschfeld charges, the charge density difference between the molecule and substrate, the density of states, and the location of the frontier orbitals (HOMO and LUMO) of a paracetamol molecule. Charge density difference calculations indicate a significant transfer of charge from the molecule to the surface. Observable regions of electron density enrichment and depletion are evident around the electron-donating and -withdrawing groups, respectively. The density of states for the paracetamol molecule evolves as it approaches the surface, and the band gap disappears upon contact with the substrate. Hirshfeld charge analysis reveals asymmetry in the charge redistribution around the molecule, highlighting the varying charging tendencies of different atoms.

On the prospects of using B4C3 as a potential electrode material for lithium-ion batteries

Lithium-ion batteries (LIBs) served as promising energy-storage solutions for various applications, from portable electronics to electric-vehicles. There is a growing demand for high-performance electrode materials with improved storage capacity, high cyclic stability, low energy barriers, reduced material toxicity, and cost-effectiveness. This study explores the electrochemical feasibility of B4C3 monolayer as an anode material for LIBs using density functional theory (DFT) computations and molecular dynamic (MD) simulations. The optimum anodic nature of B4C3 is reported based on its structural, electronic, diffusion, storage behaviour, and lithiation studies. The adsorption sites for Li atoms were systematically explored. The lithium saturation requires the adsorption of 29 lithium atoms on B4C3 while preserving stability confirmed by Ab-initio molecular dynamics (AIMD) simulations and ensuring favorable transport characteristics during the lithiation process based on Hirshfeld charge analysis and steric interaction calculations. The storage capacity of the host B4C3 is calculated as 2770.2 mAhg−1 which is notably higher than that of graphite anode. The computed value of open circuit voltage is 0.62 V per Li. The value of diffusion coefficient and ionic conductivity σ is 0.7 × 10−8 m2/s and 0.78 S/m respectively which predicts excellent overall anodic operation of the material. The transition barrier of the Li ions in the material is determined through Cl-NEB, which shows the lowest value of 0.06 eV along the x-axis.

2D Ni2P/N-doped graphene heterostructure as a Novel electrocatalyst for hydrogen evolution reaction: A computational study

The main prerequisite for designing electrocatalysts with favorable performance is to examine the links between electronic structural features and catalytic activity. In this work, Ni2P as a model electrocatalyst and one of the most potent catalysts for hydrogen evolution reaction (HER) was utilized to develop various Ni2P and carbon-based (graphene and N-doped graphene) heterostructures. The characteristics of such structures (Ni2P, graphene, N-doped graphene, Ni2P/graphene, and Ni2P/N-doped graphene), including binding energies, the projected density of states (PDOS), band structure, charge density difference, charge transfer, Hirshfeld charge analysis, and minimum-energy path (MEP) towards HER were calculated and analyzed by density functional theory (DFT) approach. The coupling energy values of hybrid systems were correlated with the magnitude of charge transfer. The main factors driving a promising water-splitting reaction were explained by the data of PDOS, band structures, and charge analysis, including the inherent electronegativity of the N species alongside shifting the Fermi level toward the conductive band. It was also shown that a significant drop occurs in the HER energy barrier on Ni2P/graphene compared to the pristine Ni2P due to N doping on the graphene layer in the Ni2P/N-doped graphene heterostructure.

Density functional theory study on the enhanced adsorption mechanism of VOCs by Cu4 cluster decorated UiO-66

Volatile organic compounds (VOCs) are a class of air pollutants that greatly harm human health and the ecological environment, convenient and effective treatment of VOCs has become a challenge. This work proposed a transition metal cluster decorated UiO-66 (a typical type of metal-organic framework) as the effective adsorbent for the purification of VOCs with typical functional groups based on density functional theory (DFT) calculations. The results show that the VOCs adsorption energy on the modified UiO-66 increased significantly, especially for Cu4 cluster-modified UiO-66, and the adsorption energy of trichloroethene is greater than that of the other considered VOCs, which may be induced by halogens or carbon-carbon double bond in trichloroethene. We compared chloromethane, chloroethylene, and ethylene and found that VOCs with carbon-carbon double bond have stronger adsorption energy, increasing with the number of chlorine atoms. The density of states (DOS), Hirshfeld charge analysis, and orbital distribution further confirm the charge distribution in UiO-66 after the decoration of the Cu4 cluster. Cu4 cluster can act as a bridge to promote the formation of chemical bonds between UiO-66 and VOCs molecules. Therefore, Cu4 cluster-modified UiO-66 can be a promising VOCs adsorption material, especially for vinyl chloride molecules.

Spectroscopic and Computational Study of ZnCl2-Methanol Low-Melting-Temperature Mixtures.

Alcoholic electrolyte mixtures have wide applications in industries. In this study, a series of mixtures composed of ZnCl2 and methanol (MeOH) with ZnCl2 mol % from 6.7 to 25 were prepared, and their spectral, structural, and thermodynamic properties were studied using infrared (IR) spectroscopy, differential scanning calorimetry (DSC), and density functional theory (DFT) calculations. The DFT-assisted analysis of excess spectra, supported by 2D-correlation spectroscopy, led to the identification of the major constituents of ZnCl2-MeOH mixtures, namely, MeOH monomer, MeOH dimer, and ZnCl2-3MeOH complex, produced after dissociation of MeOH trimer which represents the bulk methanol. The Hirshfeld charge analysis showed that in the competition between the O-H···Cl hydrogen bond (H-bond) and Zn ← O coordination bond to transfer charge in ZnCl2-MeOH complexes, the latter always dominates, making MeOH positively charged. The phase diagram of the binary system showed the presence of V-shaped glass transition temperatures (Tg), characteristic of low-melting mixture solvents (LoMMSs). The present study provides insights into the microscopic properties of the system and sheds light on the understanding of the general principles to prepare deep-eutectic solvents (DESs) or LoMMSs using inorganic salts and alcoholic compounds.

Quantum chemical study of the effect of oxygen-containing functional groups on N2O adsorption on char

This work aims to investigate the effect of oxygen-containing functional groups (OFGs) on nitrous oxide (N2O) adsorption on char based on density-functional theory (DFT) at M06-2X/6-311G(d) level using pristine and defective zigzag char models. The char model configurations were geometrically optimized and single-point energy was calculated. The adsorption energy, bond length, Mayer bond order (MBO), atomic dipole moment corrected Hirshfeld (ADCH) charges, electrostatic potential (ESP), partial density of states (PDOS) and interaction region indicator (IRI) were analyzed. Results show that the presence of OFGs has little influence on the structure of char model but increases its activity. In addition, it was found that the effect of OFGs on N2O adsorption on char surface is related to char model and adsorption configurations. Further char configurations analysis after adsorbing N2O molecules suggests that the OFGs addition can significantly affect adsorption energies by changing chemical bond length and MBO. Moreover, the ESP analysis indicates that the OFGs addition increases the overall level of ESP on char surface, leading to the formation of orbital hybridization and electrons transfer and further promoting N2O adsorption on char surface. ADCH analysis indicates that OFGs affect the active site charge and facilitate N2O adsorption. IRI analysis reveals a variety of weak interactions in the structure after the incorporation of OFGs. Therefore, this study provides a theoretical foundation for a better understanding of the heterogenous adsorption of N2O on char.

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment.

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges on 1-particle basis set and (for WFT methods) -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.

Comparative Analysis of Nickel-Phosphine Complexes with Cumulated Double Bond Ligands: Structural Insights and Electronic Interactions via ETS-NOCV and QTAIM Approaches.

This study presents a comprehensive analysis of nickel-phosphine complexes, specifically Ni(PH3)2(OCCH2), Ni(PH3)2(H2CCO), Ni(PH3)2(H2CCCH2), Ni(PH3)2(NNCH2), and Ni(PH3)2(η1-H2CNN). Utilizing ETS-NOCV analysis, we explored orbital energy decomposition and the Hirshfeld charges of the ligands, providing insights into the electronic structures and donor-acceptor interactions within these complexes. The interactions in the ketene and allene complexes exhibit similar deformation densities and NOCV orbital shapes to those calculated for Ni(PH3)2(NNCH2), indicating consistent interaction characteristics across these complexes. The total interaction energy for all η2 complexes is observed to be over 60 kcal/mol, slightly exceeding that of the analogous carbon dioxide complex reported earlier. Furthermore, the study highlights the stronger back-donation as compared to donor interactions across all η2 complexes. This is further corroborated by Hirshfeld analysis, revealing the charge distribution dynamics within the ligand fragments. The research offers new perspectives on the electron distribution and interaction energies in nickel-phosphine complexes, contributing to a deeper understanding of their catalytic and reactive behaviors.