Abstract

Different density functional methods (DFT) have been used to optimize and study the chemistry of five potential anticancer complexes in terms of their electronic, conductive, and spectroscopic properties. Many of the computed properties in addition to the IR and QTAIM analysis of the NMR are dipole moment vector ( ), linear polarizability tensor ( ), first hyperpolarizability tensors ( ), polarizability exaltation index ( ), and chemical hardness ( ) of the complexes. Stable low energy geometries are obtained using basis set with effective core potential (ECP) approximation but, in the computation of atomic or molecular properties, the metal Ru atom is better treated with higher all electron basis set like DGDZVP. The spectroscopic features like the IR of the metal-ligand bonds and the isotropic NMR shielding tensor of the coordinated atoms are significantly influenced by the chemical environment of the participating atoms. The carboxylic and pyrazole units are found to significantly enhance the polarizabilities and hyperpolarizabilities of the complexes while the chloride only improves the polarity of the complexes. Fermi contacts (FC) have the highest effect followed by the PSO among all the four Ramsey terms which defined the total spin-spin coupling constant J (HZ) of these complexes.

1. Introduction

Ruthenium-based organometallic complexes are widely applied in medical research field as anticancer compounds [113]. Two of the ruthenium complexes, KP1019 [14, 15] and NAMI-A [1619], have passed through the phase two medical tests as potential anticancer drugs and might rival cisplatin which has been the most effective and widely used anticancer agent [2024]. Many other applications for ruthenium complexes in addition to study as anticancer agents are as photoluminescence, electrochemiluminescence, catalyst, and photochemical complexes [2530]. The complexes used in this work are designed to play a dual role as anticancer and nonlinear optical (NLO) materials. Docking study of some of these complexes has been found to be promising than many of the RAPTA complexes [31]. The most widely screened ruthenium complexes as anticancer are the half-sandwich complexes [2, 3, 3237]. RAPTA complexes as type of the half-sandwich ruthenium complexes under Dyson research teams have been experimentally proven as potential anticancer agents [2224, 3842]. The complexes (Figure 1) studied in this work have similar structures with another type of half-sandwich ruthenium anticancer complexes [4352]. The main differences in these models of complexes compared to the type of complexes from the Sandler research team are the pyrazole and carboxylic units of interest that are incorporated into the bidentate or tridentate ligands coordinating to the metal through the nitrogen atoms.

Our interest in this work is to compute the chemical properties of these complexes in relation to the unique features of the carboxylic or pyrazole units, total stability, conductivity, and reactivity either as anticancer or NLO materials. Besides the limitation of the unknown targets of ruthenium anticancer complexes [14, 20, 39, 5356] which is limiting their rational design and approval, other limitations in their application are the complexity of their reaction and instability [20]. In drug design there is need for a clear understanding of the physicochemical properties of the drug candidates [57] which will enhance their rational design. To understand the chemistry of these complexes, the IR spectroscopic differences, the isotropic NMR shielding ( ), magnetizabilities ( ), the hardness ( ), hyperpolarizability ( ), polarizability exaltation index ( ) dipole ( ), and anisotropy polarizations ( ) of these five complexes are computed. In addition to the spectroscopic features of these complexes, many factors computed in this work like the polarizabilities, hyperpolarizabilities [58], hardness [59, 60], and polarizability exaltation index [59], have been linked to the stability, reactivity, or selectivity [61] and conductivity of molecules. The main aim is to characterize and analyze the electric, conductive, and spectroscopic properties of the complexes in relation to the number and position of the carboxylic and pyrazole units. There is little known information about the hyperpolarizabilities of ruthenium metal complexes but many of the computed properties in this study have not been reported for these types of metal complexes to the best of our knowledge.

2. Computational Method

The geometries of the complexes were optimized twice using PBE0 hybrid density functional [62] and combined basis set SBKJC VDZ [63] with effective core potential (ECP) (for ruthenium and chloride atoms where applicable) while other atoms are treated with basis set 6-31G in the first optimization that will be subsequently referred to as ECP(Ru,Cl)∣6-31G systems. In the second optimization, only the ruthenium atom is treated with SBKJC VDZ ECP basis set while other atoms are treated with improved basis set 6-31+G(d,p) and this will be referred to as ECP(Ru)∣6-31+G(d,p) systems. The PBE0 is obtained by casting the functional and correlation of Perdew, Burke, and Ernzerhof in a hybrid HF/DFT scheme with a fixed 1/4 ratio [64]. In application of SBKJC VDZ ECP basis set, 28 core electrons were removed from Ru (1s, 2s, 2p, 3s, 3p, and 3d) and 10 from Cl (1s, 2s, and 2p) atoms (where applicable) and were treated with pseudopotential while the valence electrons were treated with a double zeta quality functions. The choice of SBKJC VDZ ECP basis set is necessary for large systems of our type which also contain heavy metal like ruthenium. The choice combination of ECP basis set with PBE0 functional for the optimization is due to the past records of their effectiveness in computational study of metal clusters [65, 66] and because the ECP also incorporates relativistic corrections for the metal atoms [67]. This functional has been reported to improve the accuracy of excitation energies and charge transfer bands in metal complexes both in gas phase and in solution [68]. Also, the ECP basis set has been pointed out as a viable method for accurate calculations of transition metal polarizabilities [69]. All of the computed properties were done using Becke’s three-parameter exchange [70] and Lee-Yang-Parr’s correlation nonlocal functional which is usually denoted by B3LYP, combined with different combination of basis set. The properties of the systems optimized with PBE0/ECP(Ru)∣6-31+G(d,p) are computed with the combination of DGDZVP basis for ruthenium atom and 6-31+G(d,p) basis set for other atoms which will be further referred to as DGDZVP(Ru)∣6-31+G(d,p). The properties of the same set of the optimized systems are also computed with the same combinations of basis set as was used for the optimization (ECP(Ru)∣6-31+G(d,p)). The properties of all the systems optimized with lower basis sets ECP(Ru,Cl)∣6-31G were computed with the same combination of basis sets as used for their optimization (ECP(Ru,Cl)∣6-31G ) and also computed with minimal all electron basis set 3-21G [71] applied on all atoms. The external SBKJC VDZ ECP was obtained from EMSL Basis Set Library [72, 73] and incorporated into input files in a readable format for Gaussian 09 [74] suite of programs used for the computation. The NMR and one-bond NMR spin-spin coupling constants J(A,B) [75, 76] were computed in Gaussian package using the GIAO method. The wave function files (WFX) for the computation of QTAIM properties were obtained at the same level of theory through the fchk generated files using the AIMAll program package [77]. Intra- and interatomic properties were computed by the Proaim integration approach, as implemented in the AIMAll suite of programs. The accuracy of the integration process was guaranteed by keeping the atomic integral of the one-electron density Laplacian below the 10−4 a.u. within all atomic basins. In addition, the summation of all basin energies was compared with the total electronic energy of the molecule, calculated independently at the above-mentioned computational level, to check that the differences remain below the 1 kcal/mol.

The correlation of the computed properties was done with statistical package [78] while the RMSD of the geometries obtained at different optimization levels was built with UCSF Chimera package [79].

3. Results and Discussion

3.1. The Geometry and the IR Spectroscopic

The geometry is optimized twice using the combined basis set ECP(Ru, Cl)∣6-31G for the first optimization and ECP(Ru)∣6-31+G(d,p) for the second optimization. The zero imaginary obtained for the five complexes during the two methods of optimization is an indication that both methods detect the minimum geometry of the complexes except for complex 4 that was predicted during ECP(Ru,Cl)∣6-31G as transition with single imaginary frequency. There is no significant changes in the bond lengths (≫0.01 Å), angles (≫2.00), and the conformations (≫0.1) of the complexes when optimized with ECP(Ru,Cl)∣6-31G and with higher basis set ECP(Ru)∣6-31+G(d,p) (Table 1). The little significant change observed is in the conformational change of the carboxylic units of complex 1 which is responsible for it having the highest recorded RMSD of 0.068 (Table 1) when the geometries from the two optimizations are matched.

The numbers of the vibrational modes observed for each of the complexes are 102, 114, 120, 108, and 129 and are directly proportional to the number of atoms 36, 40, 42, 38, and 45 in each of the respective complexes 1, 2, 3, 4, and 5. The IR spectra of the five complexes are shown in Figure 2 and the possible assignments of the vibration to the bonds stretching and bending are shown in Table 2. Differences in the complexes and their respective IR spectra were used to suggest the corresponding bond vibrations. The Ru-nitrogen (Ru–N) bonds are assigned to indicate the differences in the vibration of ruthenium bonding to pyrazole or pyridine or phenanthroline units. The peculiarity of weak vibrations around 700 for the complexes 1, 2, and 3 with chloride ligands which are not obvious in other complexes without chloride is an indication of Ru–Cl vibration around this region. The Ru-carbon (Ru–C) bonds are assigned to vibration at 797 which is common to all the complexes. The assigned Ru–C vibrations are within the possible Metal-C vibrations (548–829) reported for carbonyl metal bonds [83] and of that which was specifically reported for the Ru–C vibration (838) [84]. The common vibrations around 1419 and that peculiar to complex 4 and at 1477 for other complexes are assigned to the possible C–C and C–N vibrations and can be equally suggested to be a result of Ru-ligand (Ru-L) charge transfer [82]. All the vibrations at 3600–3625, 1758–1784, 1050–1150, and 735–745 are all peculiar to the complexes with carboxylic unit and are assigned to different modes of vibrations as indicated in Table 2.

4. Spin-Spin Coupling (J (HZ)) of Ru-Ligand Bonds

The four Ramsey terms, Fermi contact (FC), spin dipole (SD), diamagnetic spin orbit (DSO), and paramagnetic spin orbit (PSO), were used to determine the total nuclear spin-spin coupling (J) [85]. These Ramsey terms are computed with different combinations of basis sets. The properties computed with DGDZVP(Ru)∣6-21+G(d,p) and 3-21G are shown in Supplementary Table S1, available online at http://dx.doi.org/10.1155/2014/656830, while those computed with ECP(Ru)∣6-21+G(d,p) and ECP(Ru,Cl)∣6-21G are shown in Table S2. The FC and SD represent the spin polarization densities while DSO and PSO represent orbital current densities [85]. Also, the FC and SD describe the interaction of the electronic spin with the nuclear magnetic field which consists of two parts, namely, extended dipole field outside the nucleus and a strongly localized field inside the nucleus [86]. The SD term corresponds to the dipole field, whereas the FC term describes a contact interaction between electron spin and nuclear magnetic field at the surface of the nucleus (i.e., hyperfine interaction between nucleus and electron) [85]. Long-range coupling in p systems should contain substantial SD contributions if the coupling nuclei are incorporated into the p system [86].

The Ramsey terms computed with only minimal basis set 3-21G are in a close relationship with the one computed with higher basis set DGDZVP(Ru)∣6-21+G(d,p) but it underestimates the values of the FC and SD and consequentially the values of the total spin-spin coupling constant (J (HZ)) (Table S1) while the magnitude of PSO and DMO remains virtually the same. However, the same trends of values obtained with 3-21G and DGDZVP(Ru)∣6-21+G(d,p) ( -square = 0.9864; value = 2.20E-016; -statistic = 3192) are a clear indication that the minimal basis set can equally reproduce a good feature of the Ramsey terms of these complexes. All the Ramsey properties computed with basis sets either when only the Ru atom is treated with ECP (ECP(Ru)∣6-21+G(d,p)) or when Ru and Cl atoms are treated with ECP (ECP(Ru,Cl)∣6-21G ) significantly underestimate the Ramsey terms except the DSO which appears to be unaffected with core electron approximation with ECP (Table S2). The Ramsey term that was most underestimated when ECP was used is the FC term and there was no significant improvement even when the ECP was limited to Ru metal with application of higher basis set on all other atoms in the complexes (ECP(Ru)∣6-21+G(d,p)) (Table S2). This further supports the reported characteristics of the Ramsey term which pointed out that the SD coupling mechanism requires occupied and unoccupied non-s orbitals while FC operator probes the -electrons at the sites of coupling nuclei [86, 87] which are predominantly the ECP approximated core orbitals. However, the trend and the values of the total spin-spin coupling constant (J (HZ)) of the ECP(Ru)∣6-21+G(d,p) systems improved ( -square = 0.7781; value = 7.28E-016; -statistic = 155.3) compared to the ECP(Ru,Cl)∣6-2G systems ( -square = −0.007818; value = 0.4215; -statistic = 0.6587).

The Ramsey terms which have the greatest magnitude are the FC term followed by the PSO while the DSO is the smallest. The least magnitude of DSO further supports the report that DSO is the least interesting and also the least investigated mechanism [86]. A large DSO contribution can only be expected if the charge distribution around one of the coupling nuclei strongly deviates from isotropy as can be found in certain transition metal atoms [85] which is not the case in these complexes.

The chloride atoms in complexes 1, 2, and 3 have the lowest magnitude of the FC and DSO but the highest SD and PSO which consequentially lead to the highest magnitude of the J (HZ) compare to other coordinated atoms of the ligands. The coordinated nitrogen atoms have the highest magnitude of FC but lowest magnitude of PSO and DSO compared to other coordinated atoms in the complexes but their J (HZ) is lesser than that of the chloride atoms. The features of the PSO that show the chloride atoms as having the highest magnitude followed by the arene atoms and nitrogen atom being the least are an indication that the high electronegativity effect of chloride and a high -strength of the arene strongly determine the magnitude of PSO. This is in agreement with the reports that show the possible relation of PSO coupling mechanism to bond order, -strength, electronegativity, and the magnetizability of a bond ( ) [86]. The relatively same values of the spin-spin coupling of the chloride atoms in complexes 1, 2, and 3 are an indication that there is no serious effect on the coupling strength of the chloride to the ruthenium metal atom due to changes in their chemical environment (Table S1). The coupling spin-spin constant of the coordinated nitrogen atoms of the phenanthroline unit of complex 3 has the lowest spin-spin coupling constant among all the bidentate Ru–N couplings. The coordinated nitrogen atom of the central pyridine unit of tridentate ligand of complex 4 has the lowest spin-spin coupling among all its coordinated nitrogen of the ligands while no significant variation is observed for the three coordinated nitrogen atoms of the terpyridine unit of complex 5. In all the atoms of the ligands that are coordinated to the Ru metal atom, the carbon atoms of the arene unit have the lowest magnitude of the total spin-spin coupling and a significant variation in the coupling strength of the carbon atoms of the same arene unit. There appear to be significant changes in the magnitude of the SD and consequential significant change in the coupling J (HZ) of the coordinated arene carbon atoms as the chemical environment changes from one complex to another.

In complexes 1, 2, and 3 where there is chloride atom, a very high correlation of the SD with the PSO (0.93, 0.89, and 0.98, resp.) is observed but no such high correlation is observed for complexes 4 (0.32) and 5 (0.08) without the chloride atoms. The effect of SD on the total spin-spin coupling (J (HZ)) was also high on the three complexes with chloride (Table 3) while the effect of the FC on the spin-spin coupling constant significantly improved in the complexes 4 and 5 without the chloride atoms. Generally, the Ramsey terms which have the highest effects on the total spin-spin coupling constants are FC and PSO while SD and DSO are average and depend on the type of the complexes (Table 3). Even though the magnitude of DSO is the least in all the complexes (Table S1), yet it has a significantly improved effect in complexes 4 and 5 (Table 3) that are without chloride ligand.

5. NMR Properties of the Complexes

The computed NMR properties using the combined basis set DGDZVP(Ru)∣6-31+G(d,p) (Table S3) and a single minimal basis set 3-21G are in a closed range (Table S4) but when the ECP(Ru)∣6-31+G(d,p) basis was used, the of Ru atom was underestimated which equally affects the values obtained for other atoms (Table S5). There have been reported applications of ECP to study only the chemical shifts of the ligands and it was not expected to perform well for the shielding of the heavy atom since NMR parameters depend critically on the electron density in the nuclear regions which may not be accurately modelled by the RECPs approach [88]. It is also difficult to measure accurately, and the experimental value often includes large errors. The magnitudes of the , ( ), and many other atomic properties computed with the minimal basis set 3-21G and DGDZVP(Ru)∣6-31+G(d,p) are relatively close. This observation further supports the view that the is generally insensitive to electron correlation corrections and vibrational corrections [89]. There are significant variations in the computed properties with ECP as the intra-atomic magnetizabilities ( ) of the ruthenium atom that was characterized with ECP basis set were underestimated while its or was overestimated to compensate for the underestimated intra which consequentially lead to overestimation of the of ruthenium (compare Table S3 with Table S5). This observation does not completely agree with the reported view that relativistic effect in is much smaller, since the magnetizability does not require an accurate electronic structure in the vicinity of nuclei [89]. Magnetizability is expressed as a second derivative of the total electronic energy with respect to the external magnetic field [89]. The effect of ECP is much more prominent in the computed as there was no compensation for the underestimated intra-atomic isotropic shielding ( ) which consequentially lead to the complete underestimation of the of the ruthenium atom.

The hydrogen atoms that are involved in HB are characterized with higher , lower , and lower than any other hydrogen atoms in the complex (Figure 3). The QTAIM properties obtained from the highest combination of basis set DGDZVP(6-31+G(d,p)) (Table S3) will be considered for further discussion. The magnitude of of Ru atom is higher than other atoms in the complexes but it is not the case for its total bonding isotropic NMR shielding contribution. The highest among each of the of Ru to all the coordinated atoms is the arene carbon atoms (Ru–C bonding) except in complex 3 where the highest contribution is towards the Ru–Cl bonding. The amount of contribution of Ru atom to any coordinated atom is in proportion to the amount contributed by the coordinating atoms and is also higher except for the Ru–Cl bonding where it is lower. This type of proportion does not occur in the case of the of Ru atom to the coordinated atoms. It is Ru atom Sigma_Iso contribution that is far either higher as in the case of Ru–C bonding with arene atoms or lower as in the case of Ru–N or Ru–Cl atoms bonding. The nature of each Sigma_Iso shows that the arene unit is electron deficient -ligand while the Ru, Cl, and N unit are electron rich. It has been shown that in ionic and polar systems current density flow is localized in atomic basins and the contribution of the to the total of the atom is minimal which is reverse in covalently bonded molecules where the contribution of to the is considerable and in some cases, like aromatic molecules, even surpasses the [90]. Therefore, the lowest values obtained for the Cl atom (Table S3) is an indication that it could be ionic bond. It has been demonstrated that is able to verify the exact nature of aromaticity/antiaromaticity among different molecules and also distinguish the correct aromaticity order among the sets of aromatic/antiaromatic molecules [90].

The selected correlation of the computed QTAIM NMR atomic properties with the and NMR are shown in Table 4. Besides those shown in Table 4, the values of anisotropy is highly correlation with with the correlation values of 0.69, 0.62, 0.68, 0.96, and 0.89 for the respective complexes while it is inverse correlation with with correlation values of −0.80, −0.79, −0.83, −0.86, and −0.80 for the respective complexes. It is interesting to point out that there is a very strong correlation of atomic charges with the as was previously observed in the literature [90] but with the it is low (Table 4). Also, the correlation of the atomic charges with the is very low. The of the atom is found to be highly correlated with the electronic kinetic energy of an atom , the approximation to virial-based total energy of an atom ( _Scaled( )), and the number of electrons localized in an atom (LI ) but relatively lower correlated with the number of electron delocalization from an atom (DI )) and the total dipole moment contribution of an atom (μ ).

The high correlation obtained from the computed atomic properties makes us propose a model for easy computation of NMR properties of molecules.

The proposed model equation is where , , , , and .

The values, the level of the correlation, and the significance of the QTAIM properties ( ) used to model the are shown in Table 5. From the values of the , there is a clear indication that there is a very high correlation of the with the computed electronic kinetic energy of atoms (K(A)), LI , approximation to virial-based total energy of atoms ( _Scaled ), and DI ). However, despite the high correlation, wide residual ranges are observed for the complexes 1, 2, 3, and 5 with carboxylic unit. The carboxylic group is responsible for a wide margin as the oxygen and the hydroxyl units are the most poorly fitted among all the atoms in each of the complexes.

6. Polarizability and Hyperpolarizabilities of the Complexes

The energy ( ) of a molecule that is subjected to a static electric field ( ) can be expressed as where is the energy of the molecule in the absence of an electronic field, is the component of the dipole moment vector, is the linear polarizability tensor, and are first and second hyperpolarizability tensors, where , , and label are the , , and components, respectively [68].

The terms in the energy of molecule due to the effect of the external electric fields are computed except the second hyperpolarizability tensor ( ).

The dipole was computed by The anisotropy polarization was calculated using

Besides the mean polarizability ( ) and polarizability anisotropies ( , , and ), we also evaluated the polarizability exaltation index ( ), which is determined as , where (mol) is the mean polarizability of a molecule and is the summation of the atomic polarizability of the atoms which constitute the molecule. According to literature [59], calculated values have been employed to estimate relative aromaticity of furan homologues and stabilities of atomic clusters. A large negative value denotes a very stable structure.

The accurate single values of first static hyperpolarizability ( ) of the complexes were computed using quasi-Pythagorean problem [58] on the 10 components of the 3 × 3 × 3 matrix as , , , , , , , , , and , respectively, from G03 output as where , and .

Since the values of in G09 are reported in atomic units (a.u.), the calculated total values were converted into electrostatic units (esu) (1 a.u. =8.6393 × 10−33 esu).

Hardness ( ) is a molecular property commonly employed to characterize the relative stability and reactivity of series of compounds. An operative expression of is given in terms of ionization energy (IE) and electron affinity (EA), which can be approximated by the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies ( ), respectively, The calculated values are included in Table 6 with all other conductive properties of the complexes. In agreement with the relative energies, the most stable isomer is predicted to be the hardest one [59]. However, the molecules having a small energy gap known as soft and having a large energy gap are known as hard molecules [60]. The hardness has been found to be a good indicator of more stable molecular isomers [91, 92]. It is also pointed out to be a measure of the degree of electron localization, molecular polarizability, and resonance energy [93]. Pearson observed that a large HOMO-LUMO gap increases stability which suggests that molecules seem to follow a rule of nature to arrange them so as to be as hard as possible [94].

In both cases of ECP basis sets ECP(Ru,Cl)∣6-31G and ECP (Ru)∣6-31+G(d,p), the computed conductive properties are overestimated beside the dipole (μ), polarizability exaltation index (Γ), band gap, and hardness ( ) that were underestimated. However, ECP(Ru)∣6-31+G(d,p) give better results that follow the same trend with the basis set DGDZVP(Ru)∣6-31+G(d,p). The trend of the results obtained from 3-21G is very similar to the best result obtained from DGDZVP(Ru)∣6-31+G(d,p) basis sets but there is little overestimation of the values of Γ and underestimation of other properties. The computed hyperpolarizability values of these complexes are found to be lower than that observed for RAPTA complexes [95] but the NLO properties of complex 2 are higher than some of the reported high hyperpolarizability materials of ruthenium complexes with redox-switching noninnocent ligands (NILs) [68] and many of the h5-monocyclopentadienylnitrilecobalt complexes [96]. This implies that these complexes will have intense intramolecular CT between metal and ligand (MLCT or LMCT) which is responsible for their relatively large value of as a result of back-donation in the complexes [96, 97]. Complex 2 has the highest value of the first hyperpolarizabilities ( ) while complex 5 has the least. The trend of the magnitude of computed conductivity properties of the complexes for dipole, , , , , , , and band gap (or ) are 3 < 1 = 2 < 5 < 4; 5 < 3 < 2 < 4 < 1; 2 < 5 < 3 < 4 < 1; 5 < 2 < 3 < 4 < 1; 4 ≪ 5 < 1 < 3 < 2; 2 < 3 < 5 < 4 < 1; 2 ≪ 3 < 1 < 4 < 5; and 1 < 4 < 5 < 2 < 3, respectively. The computed value of using different basis sets does not follow the same trend which is an indication that the values of are very sensitive to the type of the basis sets. The polarizabilities properties of complexes 2 and 5 are relatively high except for the which is an indication that the pyridine as bidentate and tridentate in the respective complexes contributes to the polarizabilities. The complex 4 has the lowest dipole moment because there is no carboxylic unit which would have improved its polarity which is responsible also for its poor polarizabilities, first static hyperpolarizabilities ( ), and higher band gap. Since the values of Γ can be used to predict the stability of the complexes [59], therefore, complexes 2 and 3 may probably be the most stable and coupling with their lowest band gap (very soft) will make them very reactive [98] which may eventually help in their anticancer activities. Also, complexes 2 and 3 have the highest values of which is an indication that they are the best NLO materials which may eventually help in their anticancer activities.

7. Conclusions

The electronic, conductive, and spectroscopic properties of five ruthenium-based complexes are computed using density functional theory (DFT). The properties computed are related to the unique features of their carboxylic or pyrazole units, the total stability, conductivity, and reactivity of these complexes as either anticancer or NLO materials. Both applications of the ECP basis sets ECP(Ru,Cl)∣6-31G and ECP(Ru)∣6-31+G(d,p) with the hybrid DFT functional PBE0 locate the minimum geometries of the complexes except the complex 4 that was predicted as transition when optimized with ECP(Ru,Cl)∣6-31G basis sets. There are no significant changes in the geometries when comparing the geometries from the two optimization methods. Though application of the ECP basis sets give optimized stable geometries of the complexes but computing the properties like the , polarizabilities, and hyperpolarizabilities with ECP gives either overestimated or underestimated values. The properties computed with the minimal basis set 3-21G are very similar to the one computed with the combination of higher basis set DGDZVP(Ru)∣6-31+G(d,p) which is an indication that minimum basis set can be an advantage in the computation of the large molecules of these types. The difference in the IR features of these complexes shows that there is either bathochromic or hypsochromic shift in the IR vibration modes of the metal-ligand bonds like Ru–N as a result of the change in the chemical environment of participating atoms.

The Ramsey term which has the greatest magnitude and effect on the computed spin-spin coupling constant J (HZ) is the FC term followed by the PSO and SD while the DSO is least. The J (HZ) properties of the Ru–Cl are the highest among all the present Ru-Ligand bonds and do not significantly change due to change in the chemical environment of the complexes. Next to it is the Ru–N, but the magnitude of its J (HZ) changes with the change in the chemical environment as the Ru–N of complex 3 with phenanthroline environment is characterized with the lowest magnitude among all the computed J (HZ) of Ru–N in the complexes.

The hydrogen atoms that are involve in hydrogen bonding (HB) are characterized with higher and and lower than any other hydrogen atom in the complexes. The isotropic bonding contribution of Ru atom to the coordinated atoms ( ) is far higher as in the case of Ru–C bonding with arene atoms or lower as in the case of Ru–N or Ru–Cl atoms bonding. The nature of each shows that the arene unit is electron deficient π-ligand while the Ru, Cl, and N units are electron rich. Also the lowest values obtained for the Cl atom are an indication that it could be ionic bond. The values of anisotropy show high correlation with but inverse correlation with . The is found to be highly correlated with the electronic kinetic energy of an atom , the approximation to virial-based total energy of an atom ( _Scaled ), and the number of electrons localized in an atom (Loc ).

Other computed properties show that complexes 2 and 3 with two units of carboxylic are better NLO materials characterized with lower band gap and higher hyperpolarizabilities. Also, complexes 2 and 3 are predicted as the most stable and characterized with the highest value of Γ and the most reactive with lowest band gap (very soft). Also, the pyrazole units appear to contribute to the NLO features of complexes 1 and 4 which are found to be better NLO materials than complex 5 despite the absence of carboxylic unit in complex 4. The presence of the chloride ligand in complexes 1, 2, and 3 does not have any effect on the conductive properties computed except that it enhances the polarity of the complexes resulting in their higher dipole above complexes 4 and 5 without the chloride ligand.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledged the financial support of Govan Mbeki Research and Development Centre, University of Fort Hare, South Africa. The CHPC in South Africa is gratefully acknowledged for providing the computing facilities and some of the software programs that are used for the computation.

Supplementary Materials

Supplementary Table S1: The values of the Ramsey terms that define all the metal-ligand spin spin coupling of the complexes from two different systems treated with DGDZVP(Ru)|6-31+G(d,p) and 3-21G respectively

Table S2: The values of the Ramsey terms that define all the metal-ligand spin spin coupling of the complexes from two different systems treated with ECP(Ru)|6-31+G(d,p) and ECP(Ru,Cl)|6-31G* respectively

Supplementary Table S3: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis set DGDZVP(Ru)|6-31+G(d,p)

Supplementary Table S4: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis sets 3-21G

Supplementary Table S5: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis sets ECP(Ru)|6-31+G(d,p)

  1. Supplementary Materials