Research Article  Open Access
The Spectroscopic and Conductive Properties of Ru(II) Complexes with Potential Anticancer Properties
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 metalligand 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 spinspin coupling constant J (HZ) of these complexes.
1. Introduction
Rutheniumbased organometallic complexes are widely applied in medical research field as anticancer compounds [1–13]. Two of the ruthenium complexes, KP1019 [14, 15] and NAMIA [16–19], 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 [20–24]. Many other applications for ruthenium complexes in addition to study as anticancer agents are as photoluminescence, electrochemiluminescence, catalyst, and photochemical complexes [25–30]. 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 halfsandwich complexes [2, 3, 32–37]. RAPTA complexes as type of the halfsandwich ruthenium complexes under Dyson research teams have been experimentally proven as potential anticancer agents [22–24, 38–42]. The complexes (Figure 1) studied in this work have similar structures with another type of halfsandwich ruthenium anticancer complexes [43–52]. 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, 53–56] 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 631G in the first optimization that will be subsequently referred to as ECP(Ru,Cl)∣631G 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 631+G(d,p) and this will be referred to as ECP(Ru)∣631+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 threeparameter exchange [70] and LeeYangParr’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)∣631+G(d,p) are computed with the combination of DGDZVP basis for ruthenium atom and 631+G(d,p) basis set for other atoms which will be further referred to as DGDZVP(Ru)∣631+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)∣631+G(d,p)). The properties of all the systems optimized with lower basis sets ECP(Ru,Cl)∣631G were computed with the same combination of basis sets as used for their optimization (ECP(Ru,Cl)∣631G) and also computed with minimal all electron basis set 321G [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 onebond NMR spinspin 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 oneelectron 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 abovementioned 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)∣631G for the first optimization and ECP(Ru)∣631+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)∣631G as transition with single imaginary frequency. There is no significant changes in the bond lengths (≫0.01 Å), angles (≫2.0^{0}), and the conformations (≫0.1) of the complexes when optimized with ECP(Ru,Cl)∣631G and with higher basis set ECP(Ru)∣631+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 values for the complexes 4 and 5 with superscript “#” indicated their association with the middle nitrogen of their tridentate. 
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 Runitrogen (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 Rucarbon (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 MetalC 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 Ruligand (RuL) 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. SpinSpin Coupling (J (HZ)) of RuLigand 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 spinspin coupling (J) [85]. These Ramsey terms are computed with different combinations of basis sets. The properties computed with DGDZVP(Ru)∣621+G(d,p) and 321G are shown in Supplementary Table S1, available online at http://dx.doi.org/10.1155/2014/656830, while those computed with ECP(Ru)∣621+G(d,p) and ECP(Ru,Cl)∣621G 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]. Longrange 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 321G are in a close relationship with the one computed with higher basis set DGDZVP(Ru)∣621+G(d,p) but it underestimates the values of the FC and SD and consequentially the values of the total spinspin 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 321G and DGDZVP(Ru)∣621+G(d,p) (square = 0.9864; value = 2.20E016; 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)∣621+G(d,p)) or when Ru and Cl atoms are treated with ECP (ECP(Ru,Cl)∣621G) 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)∣621+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 nons 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 spinspin coupling constant (J (HZ)) of the ECP(Ru)∣621+G(d,p) systems improved (square = 0.7781; value = 7.28E016; statistic = 155.3) compared to the ECP(Ru,Cl)∣62G 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 spinspin 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 spinspin constant of the coordinated nitrogen atoms of the phenanthroline unit of complex 3 has the lowest spinspin 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 spinspin 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 spinspin 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 spinspin coupling (J (HZ)) was also high on the three complexes with chloride (Table 3) while the effect of the FC on the spinspin 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 spinspin 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)∣631+G(d,p) (Table S3) and a single minimal basis set 321G are in a closed range (Table S4) but when the ECP(Ru)∣631+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 321G and DGDZVP(Ru)∣631+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 intraatomic 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 intraatomic 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(631+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].
(a)
(b)
(c)
(d)
(e)
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 virialbased 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 virialbased 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 quasiPythagorean 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 HOMOLUMO 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)∣631G and ECP (Ru)∣631+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)∣631+G(d,p) give better results that follow the same trend with the basis set DGDZVP(Ru)∣631+G(d,p). The trend of the results obtained from 321G is very similar to the best result obtained from DGDZVP(Ru)∣631+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 redoxswitching noninnocent ligands (NILs) [68] and many of the h5monocyclopentadienylnitrilecobalt 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 backdonation 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 rutheniumbased 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)∣631G and ECP(Ru)∣631+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)∣631G 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 321G are very similar to the one computed with the combination of higher basis set DGDZVP(Ru)∣631+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 metalligand 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 spinspin 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 RuLigand 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 virialbased 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 metalligand spin spin coupling of the complexes from two different systems treated with DGDZVP(Ru)631+G(d,p) and 321G respectively
Table S2: The values of the Ramsey terms that define all the metalligand spin spin coupling of the complexes from two different systems treated with ECP(Ru)631+G(d,p) and ECP(Ru,Cl)631G* respectively
Supplementary Table S3: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis set DGDZVP(Ru)631+G(d,p)
Supplementary Table S4: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis sets 321G
Supplementary Table S5: The atomic NMR magnetazability and shielding for the selected atoms in the complexes using combined basis sets ECP(Ru)631+G(d,p)
References
 J. Ruiz, C. Vicente, C. de Hero, and D. Bautista, “A novel ruthenium(II) arene based intercalator with potent anticancer activity,” Dalton Transactions, vol. 26, pp. 5071–5073, 2009. View at: Publisher Site  Google Scholar
 A. Martinez, J. Suarez, T. Sand, R. S. Magliozzi, and R. A. SanchezDelgado, “Interactions of areneRu(II)chloroquine complexes of known antimalarial and antitumor activity with human serum albumin (HSA) and transferring,” Journal of Inorganic Biochemistry, vol. 105, no. 1, pp. 39–45, 2011. View at: Google Scholar
 C. Gossans, I. Tavernelli, and U. Rothlisberger, “DNA structural distortions induced by rutheniumarene anticancer compounds,” Journal of the American Chemical Society, vol. 130, pp. 10921–10928, 2008. View at: Publisher Site  Google Scholar
 N. Gligorijević, S. Aranđelović, L. Filipović et al., “Picolinate ruthenium(II)–arene complex with in vitro antiproliferative and antimetastatic properties: Comparison to a series of ruthenium(II)–arene complexes with similar structure,” Journal of Inorganic Biochemistry, vol. 108, pp. 53–61, 2012. View at: Publisher Site  Google Scholar
 E. Alessio, G. Mestroni, A. Bergamo, and G. Sava, “Ruthenium antimetastatic agents,” Current Topics in Medicinal Chemistry, vol. 4, no. 15, pp. 1525–1535, 2004. View at: Google Scholar
 A. Bergamo, C. Gaiddon, J. H. Schellens, J. H. Beijnen, and G. Sava, “Approaching tumour therapy beyond platinum drugs: status of the art and perspectives of ruthenium drug candidates,” Journal of Inorganic Biochemistry, vol. 106, pp. 90–99, 2012. View at: Publisher Site  Google Scholar
 I. Bratsos, D. Urankar, E. Zangrando et al., “1(2Picolyl)substituted 1,2,3triazole as novel chelating ligand for the preparation of ruthenium complexes with potential anticancer activity,” Dalton Transactions, vol. 40, pp. 5188–5199, 2011. View at: Publisher Site  Google Scholar
 M. Brindell, E. Kuliś, S. K. Elmroth, K. Urbańska, and G. Stochel, “Lightinduced anticancer activity of [RuCl_{2}(DMSO)_{4}] complexes,” Journal of Medicinal Chemistry, vol. 48, no. 23, pp. 7298–7304, 2005. View at: Publisher Site  Google Scholar
 V. Rajendiran, M. Murali, E. Suresh, M. Palaniandavar, V. S. Periasamy, and M. A. Akbarsha, “Noncovalent DNA binding and cytotoxicity of certain mixedligand ruthenium(II) complexes of 2,2′dipyridylamine and diimines,” Dalton Transactions, no. 16, pp. 2157–2170, 2008. View at: Publisher Site  Google Scholar
 A. Romerosa, T. CamposMalpartida, C. Lidrissi et al., “Synthesis, characterization, and DNA binding of new watersoluble cyclopentadienyl ruthenium(II) complexes incorporating phosphines,” Inorganic Chemistry, vol. 45, no. 3, pp. 1289–1298, 2006. View at: Publisher Site  Google Scholar
 W. F. Schmid, R. O. John, G. Mühlgassner et al., “Metalbased paullones as putative CDK inhibitors for antitumor chemotherapy,” Journal of Medicinal Chemistry, vol. 50, no. 25, pp. 6343–6355, 2007. View at: Google Scholar
 S. K. Singh, S. Joshi, A. R. Singh, J. K. Saxena, and D. S. Pandey, “DNA binding and topoisomerase II inhibitory activity of watersoluble ruthenium(II) and rhodium(III) complexes,” Inorganic Chemistry, vol. 46, pp. 10869–10876, 2007. View at: Publisher Site  Google Scholar
 C. Tan, J. Liu, L. Chen, S. Shi, L. Ji, and J. Inorg, “Synthesis, structural characteristics, DNA binding properties and cytotoxicity studies of a series of Ru(III) complexes,” Journal of Inorganic Biochemistry, vol. 102, no. 8, pp. 1644–1653, 2008. View at: Google Scholar
 P. Heffeter, K. Bock, B. Atil et al., “Intracellular protein binding patterns of the anticancer ruthenium drugs KP_{1019}and KP_{1339},” Journal of Biological Inorganic Chemistry, vol. 15, pp. 737–748, 2010. View at: Publisher Site  Google Scholar
 G. Zhao and H. Lin, “Metal complexes with aromatic Ncontaining ligands as potential agents in cancer treatment,” Current Medicinal Chemistry—AntiCancer Agents, vol. 5, pp. 137–147, 2005. View at: Publisher Site  Google Scholar
 M. Bacac, A. C. Hotze, K. van der et al., “The hydrolysis of the anticancer ruthenium complex NAMIA affects its DNA binding and antimetastatic activity: an NMR evaluation.,” Journal of Inorganic Biochemistry, vol. 98, pp. 402–412, 2004. View at: Publisher Site  Google Scholar
 A. H. Velders, A. Bergamo, E. Alessio et al., “Synthesis and chemical−pharmacological characterization of the antimetastatic NAMIAtype Ru(III) complexes (Hdmtp)[transRuCl_{4}(dmsoS)(dmtp)], (Na)[transRuCl_{4}(dmsoS)(dmtp)], and [merRuCl_{3}(H_{2}O)(dmsoS)(dmtp)] (dmtp = 5,7Dimethyl[1,2,4]triazolo[1,5a]pyrimidine),” Journal of Medicinal Chemistry, vol. 47, no. 5, pp. 1110–1121, 2004. View at: Publisher Site  Google Scholar
 A. V. Vargiu, A. Robertazzi, A. Magistrato, P. Ruggerone, and P. Carloni, “The hydrolysis mechanism of the anticancer ruthenium drugs NAMIA and ICR investigated by DFTPCM calculations,” The Journal of Physical Chemistry B, vol. 112, pp. 4401–4409, 2008. View at: Publisher Site  Google Scholar
 P. Mura, M. Camalli, L. Messori, F. Piccioli, P. Zanello, and M. Corsini, “Synthesis, structural characterization, solution chemistry, and preliminary biological studies of the ruthenium(III) complexes [TzH][transRuCl_{4}(Tz)_{2}] and [TzH][transRuCl_{4}(DMSO)(Tz)]·(DMSO), the thiazole analogues of antitumor ICR and NAMIA,” Inorganic Chemistry, vol. 43, pp. 3863–3670, 2004. View at: Publisher Site  Google Scholar
 C. S. Allardyce, A. Dorcier, C. Scolaro, and P. J. Dyson, “Development of organometallic (organotransition metal) pharmaceuticals,” Applied Organometallic Chemistry, vol. 19, no. 1, pp. 1–10, 2005. View at: Publisher Site  Google Scholar
 P. J. Dyson and G. Sava, “Metalbased antitumour drugs in the post genomic era,” Dalton Transactions, vol. 23, no. 16, pp. 1929–1933, 2006. View at: Google Scholar
 S. Chatterjee, S. Kundu, A. Bhattacharyya, C. G. Hartinger, and P. J. Dyson, “The ruthenium(II)arene compound RAPTAC induces apoptosis in EAC cells through mitochondrial and p53JNK pathways,” Journal of Biological Inorganic Chemistry, vol. 13, pp. 1149–1155, 2008. View at: Publisher Site  Google Scholar
 M. Groessl, Y. O. Tsybin, C. G. Hartinger, B. K. Keppler, and P. J. Dyson, “Ruthenium versus platinum: interactions of anticancer metallodrugs with duplex oligonucleotides characterised by electrospray ionisation mass spectrometry,” Journal of Biological Inorganic Chemistry, vol. 15, no. 5, pp. 677–688, 2009. View at: Google Scholar
 A. E. Egger, C. G. Hartinger, A. K. Renfrew, and P. J. Dyson, “Metabolization of [Ru(η^{6}C_{6}H_{5}CF_{3})(pta)Cl_{2}]: a cytotoxic RAPTAtype complex with a strongly electron withdrawing arene ligand,” Journal of Biological Inorganic Chemistry, vol. 15, pp. 919–927, 2010. View at: Publisher Site  Google Scholar
 P. Byabartta, S. K. Jasimuddin, G. Mostafa, T.H. Lu, and C. Sinha, “The synthesis, spectral studies and electrochemistry of 1,10(phenanthroline)bis{1alkyl2(arylazo)imidazole}ruthenium(II) perchlorate. Single crystal Xray structure of [Ru(phen)(HaaiMe)_{2}](ClO_{4})_{2} [phen=1,10phenanthroline, HaaiMe=1methyl2(phenylazo)imidazole],” Polyhedron, vol. 22, pp. 849–859, 2003. View at: Publisher Site  Google Scholar
 H. A. J. AlLawati, F. E. O. Suliman, S. M. Z. Al Kindy, A. M. AlLawati, G. B. Varma, and I. E. M. Nour, “Enhancement of on chip chemiluminescence signal intensity of tris(1,10phenanthroline)ruthenium(II) peroxydisulphate system for analysis of chlorpheniramine maleate in pharmaceutical formulations,” Talanta, vol. 82, pp. 1999–2002, 2010. View at: Publisher Site  Google Scholar
 X. M. Chen, G. H. Wu, J. M. Chen et al., “A novel electrochemiluminescence sensor based on bis(2,2 ′bipyridine)5amino1,10phenanthroline ruthenium(II) covalently combined with graphite oxide,” Biosensors and Bioelectronics, vol. 26, no. 2, pp. 872–876, 2010. View at: Google Scholar
 H. Wei, J. Yin, and E. Wang, “Bis(2,2'bipyridine)(5,6epoxy5,6dihydro[1,10] phenanthroline)ruthenium: synthesis and electrochemical and electrochemiluminescence characterization,” Analytical Chemistry, vol. 80, pp. 5635–5639, 2008. View at: Publisher Site  Google Scholar
 E. Rajalakshmanan and V. Alexander, “Synthesis, luminescence, and electrochemical studies of tris (homoleptic) ruthenium (II) and osmium (II) complexes of 6′tolyl2, 2′: 4′, 2′′terpyridine,” Inorganic Chemistry, vol. 46, pp. 6252–6260, 2007. View at: Google Scholar
 H. J. Bolink, L. Cappelli, E. Coronado, and P. Gaviña, “Observation of electroluminescence at room temperature from a ruthenium(II) bisterpyridine complex and its use for preparing lightemitting electrochemical cells,” Inorganic Chemistry, vol. 44, no. 17, pp. 5966–5968, 2005. View at: Publisher Site  Google Scholar
 A. A. Adeniyi and P. A. Ajibade, “Inhibitory activities and possible anticancer targets of Ru(II)based complexes using computational docking method,” Journal of Molecular Graphics and Modelling, vol. 38, pp. 60–69, 2012. View at: Publisher Site  Google Scholar
 S. GrgurićSipka, I. Ivanović, G. Rakić et al., “Ruthenium(II)–arene complexes with functionalized pyridines: synthesis, characterization and cytotoxic activity,” European Journal of Medicinal Chemistry, vol. 45, pp. 1051–1058, 2010. View at: Publisher Site  Google Scholar
 B. Biersack, M. Zoldakova, K. Effenberger, and R. Schobert, “Synthesis of [1,2,4]triazolo[1,5a]pyridines of potential PGE2 inhibitory properties,” European Journal of Medicinal Chemistry, vol. 44, no. 5, pp. 1972–1975, 2010. View at: Google Scholar
 F. Linares, M. A. Galindo, S. Galli, M. A. Romero, J. A. Navarro, and E. Barea, “Tetranuclear coordination assemblies based on halfsandwich ruthenium(II) complexes: noncovalent binding to DNA and cytotoxicity,” Inorganic Chemistry, vol. 48, no. 15, pp. 7413–7420, 2009. View at: Publisher Site  Google Scholar
 S. Das, S. Sinha, R. Britto, K. Somasundaram, and A. G. Samuelson, “Cytotoxicity of half sandwich ruthenium(II) complexes with strong hydrogen bond acceptor ligands and their mechanism of action,” Journal of Inorganic Biochemistry, vol. 104, pp. 93–104, 2010. View at: Publisher Site  Google Scholar
 T. Gianferrara, I. Bratsos, E. Iengo et al., “Synthetic strategies towards rutheniumporphyrin conjugates for anticancer activity,” Dalton Transactions, vol. 38, pp. 10742–10756, 2009. View at: Google Scholar
 E. Meggers, G. E. AtillaGokcumen, K. Gründler, C. Frias, and A. Prokop, “Inert ruthenium halfsandwich complexes with anticancer activity,” Dalton Transactions, no. 48, pp. 10882–10888, 2009. View at: Publisher Site  Google Scholar
 A. Casini, C. Gabbiani, F. Sorrentino et al., “Emerging protein targets for anticancer metallodrugs: inhibition of thioredoxin reductase and cathepsin B by antitumor ruthenium (II)arene compounds,” Journal of Medicinal Chemistry, vol. 51, no. 21, pp. 6773–6781, 2008. View at: Publisher Site  Google Scholar
 W. H. Ang, E. Daldini, C. Scolaro, R. Scopelliti, L. JuilleratJeannerat, and P. J. Dyson, “Development of organometallic ruthenium−arene anticancer drugs that resist hydrolysis,” Inorganic Chemistry, vol. 45, pp. 9006–9013, 2006. View at: Publisher Site  Google Scholar
 C. A. Vock, W. H. Ang, C. Scolaro et al., “Development of ruthenium antitumor drugs that overcome multidrug resistance mechanisms,” Journal of Medicinal Chemistry, vol. 50, no. 9, pp. 2166–2175, 2007. View at: Publisher Site  Google Scholar
 W. H. Ang, E. Daldini, L. JuilleratJeanneret, and P. J. Dyson, “Strategy to tether organometallic rutheniumarene anticancer compounds to recombinant human serum albumin,” Inorganic Chemistry, vol. 46, no. 22, pp. 9048–9050, 2007. View at: Publisher Site  Google Scholar
 W. H. Ang, A. de Luca, C. ChapuisBernasconi, L. JuilleratJeanneret, M. Lo Bello, and P. J. Dyson, “Organometallic ruthenium inhibitors of glutathioneStransferase P11 as anticancer drugs,” ChemMedChem, vol. 2, no. 12, pp. 1799–1806, 2007. View at: Publisher Site  Google Scholar
 F. Barragán, P. LópezSenín, and L. Salassa, “Photocontrolled DNA binding of a receptortargeted organometallic ruthenium(II) complex,” Journal of the American Chemical Society, vol. 133, pp. 14098–14108, 2011. View at: Publisher Site  Google Scholar
 T. Bugarcic, A. Habtemariam, J. Stepankova et al., “The contrasting chemistry and cancer cell cytotoxicity of bipyridine and bipyridinediol ruthenium(II) arene complexes,” Inorganic Chemistry, vol. 47, no. 24, pp. 11470–11486, 2008. View at: Google Scholar
 M. CastellanoCastillo, H. Kostrhunova, V. Marini et al., “Binding of mismatch repair protein MutS to mispaired DNA adducts of intercalating ruthenium(II) arene complexes,” Journal of Biological Inorganic Chemistry, vol. 13, no. 6, pp. 993–999, 2008. View at: Google Scholar
 S. J. Dougan, A. Habtemariam, S. E. McHale, S. Parsons, and P. J. Sadler, “Catalytic organometallic anticancer complexes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, pp. 11628–11633, 2008. View at: Publisher Site  Google Scholar
 S. J. Dougan, M. Melchart, A. Habtemariam, S. Parsons, and P. J. Sadler, “A Neutral spirocyclic hexacoordinated germanium(IV) complex: hypervalent germanium compounds with sulfurcontaining eightmembered rings,” Inorganic Chemistry, vol. 45, pp. 10882–10894, 2006. View at: Publisher Site  Google Scholar
 R. Fernández, M. Melchart, A. Habtemariam, S. Parsons, and P. J. Sadler, “Use of chelating ligands to tune the reactive site of halfsandwich ruthenium(II)–arene anticancer complexes,” Chemistry, vol. 10, no. 20, pp. 5173–5179, 2004. View at: Publisher Site  Google Scholar
 Y. Fu, A. Habtemariam, A. M. Basri, D. Braddick, G. J. Clarkson, and P. J. Sadler, “Structure–activity relationships for organometallic osmium arene phenylazopyridine complexes with potent anticancer activity,” Dalton Transactions, no. 40, pp. 10553–10562, 2011. View at: Publisher Site  Google Scholar
 A. Habtemariam, M. Melchart, R. Fernandez et al., “Structureactivity relationships for cytotoxic ruthenium(II) arene complexes containing N,N, N,O, and O,Ochelating ligands,” Journal of Medicinal Chemistry, vol. 49, no. 23, pp. 6858–6868, 2006. View at: Google Scholar
 H. Kostrhunova, J. Florian, O. Novakova et al., “DNA interactions of monofunctional organometallic osmium(II) antitumor complexes in cellfree media,” Journal of Medicinal Chemistry, vol. 51, pp. 3635–3643, 2008. View at: Google Scholar
 H. K. Liu, F. Wang, J. A. Parkinson, J. Bella, and P. J. Sadler, “Ruthenation of duplex and singlestranded d(CGGCCG) by organometallic anticancer complexes,” Chemistry, vol. 12, pp. 6151–6165, 2006. View at: Publisher Site  Google Scholar
 G. Sava, A. Bergamoa, and P. J. Dyson, “Metalbased antitumour drugs in the postgenomic era: what comes next?” Dalton Transactions, vol. 40, no. 36, pp. 9069–9075, 2011. View at: Publisher Site  Google Scholar
 G. Gasser, I. Ott, and N. MetzlerNolte, “Organometallic anticancer compounds,” Journal of Medicinal Chemistry, vol. 54, no. 1, pp. 3–25, 2011. View at: Publisher Site  Google Scholar
 S. Page, “Education in chemistry,” 2012, http://www.rsc.org/eic/. View at: Google Scholar
 C. S. Allardyce and P. J. Dyson, “Ruthenium in medicine: current clinical uses and future prospects,” Platinum Metals Review, vol. 45, no. 2, pp. 62–69, 2001. View at: Google Scholar
 J. Pratuangdejkul, P. Jaudon, and O. C. Ducrocq, “Cation$\pi $ interactions in serotonin: conformational, electronic distribution, and energy decomposition analysis,” Journal of Chemical Theory and Computation, vol. 2, pp. 746–760, 2006. View at: Google Scholar
 P. S. Liyanage, R. M. de Silva, and K. M. N. de Silva, “Nonlinear optical (NLO) properties of novel organometallic complexes: high accuracy density functional theory (DFT) calculations,” Journal of Molecular Structure, vol. 639, pp. 195–201, 2003. View at: Publisher Site  Google Scholar
 A. Alparone, “Dipole (hyper)polarizabilities of fluorinated benzenes: an ab initio investigation,” Journal of Fluorine Chemistry, vol. 144, pp. 94–101, 2012. View at: Publisher Site  Google Scholar
 M. Karabacak and M. Cinar, “FTIR, FTRaman, UV spectra and DFT calculations on monomeric and dimeric structure of 2amino5bromobenzoic acid,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 86, pp. 590–599, 2012. View at: Publisher Site  Google Scholar
 E. E. Ebenso, T. Arslan, F. Kandemirli, N. Caner, and I. Love, “Quantum chemical studies of some rhodanine azosulpha drugs as corrosion inhibitors for mild steel in acidic medium,” International Journal of Quantum Chemistry, vol. 110, no. 5, pp. 1003–1018, 2010. View at: Google Scholar
 C. Adamo and V. Barone, “Toward reliable density functional methods without adjustable parameters: the PBE0 model,” The Journal of Chemical Physics, vol. 110, no. 13, pp. 6158–6170, 1999. View at: Publisher Site  Google Scholar
 W. J. Stevens, M. Krauss, H. Basch, and P. G. Jasien, “Relativistic compact effective potentials and efficient, sharedexponent basis sets for the third, fourth, and fifthrow atoms,” Canadian Journal of Chemistry, vol. 70, pp. 612–630, 1992. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, pp. 3865–3868, 1996. View at: Publisher Site  Google Scholar
 B. Marchal, P. Carbonniére, D. Begue, C. Pouchan, and P. Carbonniére, “Structural and vibrational determination of small galliumarsenide clusters from CCSD(T) and DFT calculations,” Chemical Physics Letters, vol. 453, no. 1–3, pp. 49–54, 2008. View at: Google Scholar
 R. Marchal, P. Carbonnière, and C. Pouchan, “Structural and vibrational properties prediction of Sn_{n}Te_{n} clusters (n = 2–8) using the GSAM approach,” Computational and Theoretical Chemistry, vol. 990, pp. 100–105, 2012. View at: Publisher Site  Google Scholar
 A. Bagno and M. Bonchio, “DFT calculations of ^{99}Ru chemical shifts with allelectron and effective core potential basis sets,” European Journal of Inorganic Chemistry, vol. 2002, no. 6, pp. 1475–1483, 2002. View at: Publisher Site  Google Scholar
 Y. Liu, C.G. Liu, S.L. Sun, G.C. Yang, and Y.Q. Qiu, “Redoxswitching secondorder nonlinear optical responses of N^{∧}N^{∧}N ruthenium complexes,” Computational and Theoretical Chemistry, vol. 979, pp. 112–118, 2012. View at: Publisher Site  Google Scholar
 N. P. Labello, A. M. Ferreira, and H. A. Kurtz, “Utilizing relativistic effective core potentials for accurate calculations of molecular polarizabilities on transition metal compounds,” The Journal of Physical Chemistry A, vol. 110, pp. 13507–13513, 2006. View at: Publisher Site  Google Scholar
 A. D. Becke, “Density‐functional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, pp. 5648–5652, 1993. View at: Publisher Site  Google Scholar
 K. D. Dobbs and W. J. Hehre, “Molecular orbital theory of the properties of inorganic and organometallic compounds 5. Extended basis sets for firstrow transition metals,” Journal of Computational Chemistry, vol. 6, pp. 861–879, 1987. View at: Google Scholar
 D. J. Feller, “The role of databases in support of computational chemistry calculations,” Journal of Computational Chemistry, vol. 17, pp. 1571–1586, 1996. View at: Publisher Site  Google Scholar
 K. L. Schuchardt, B. T. Didier, T. Elsethagen et al., “Basis set exchange: a community database for computational sciences,” Journal of Chemical Information and Modeling, vol. 47, no. 3, pp. 1045–1052, 2007. View at: Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Gaussian, Inc., Wallingford, Conn, USA, 2009.
 T. Helgaker, M. Watson, and N. C. Handy, “Analytical calculation of nuclear magnetic resonance indirect spin–spin coupling constants at the generalized gradient approximation and hybrid levels of densityfunctional theory,” Journal of Chemical Physics, vol. 113, pp. 9402–9409, 2000. View at: Publisher Site  Google Scholar
 V. Sychrovsky, J. Grafenstein, and D. Cremer, “Nuclear magnetic resonance spinspin coupling constants from coupled perturbed density functional theory,” The Journal of Chemical Physics, vol. 113, no. 9, pp. 3530–3547, 2000. View at: Publisher Site  Google Scholar
 T. A. Keith, “AIMAll (Version 12.06.03),” TK Gristmill Software, Overland Park, Kan, USA, 2012, http://aim.tkgristmill.com/. View at: Google Scholar
 R Development Core Team, “A language and environment for statistical computing,” R Foundation for Statistical Computing, Vienna, Austria, 2009, http://www.Rproject.org. View at: Google Scholar
 E. F. Pettersen, T. D. Goddard, C. C. Huang, G. S. Couch, E. C. Meng, and T. E. Ferrin, “UCSF Chimera: a visualization system for exploratory research and analysis,” Journal of Computational Chemistry, vol. 25, pp. 1605–1612, 2004. View at: Google Scholar
 F. Caruso, M. Rossi, A. Benson et al., “Ruthenium–arene complexes of curcumin: xray and density functional theory structure, synthesis, and spectroscopic characterization, in vitro antitumor activity, and DNA docking studies of (pCymene)ru(curcuminato)chloro,” Journal of Medicinal Chemistry, vol. 55, pp. 1072–1081, 2012. View at: Publisher Site  Google Scholar
 C. S. K. Rajapakse, A. Martinez, B. Naoulou et al., “Synthesis, characterization, and in vitro antimalarial and antitumor activity of new ruthenium(II) complexes of chloroquine,” Inorganic Chemistry, vol. 48, pp. 1122–1131, 2009. View at: Publisher Site  Google Scholar
 K. M. Omberg, J. R. Schoonover, J. A. Treadway, R. M. Leasure, R. B. Dyer, and T. J. Meyer, “Midinfrared spectrum of [Ru(bpy)3]^{2+},” Journal of the American Chemical Society, vol. 119, no. 30, pp. 7013–7018, 1997. View at: Google Scholar
 G. Bor, U. K. Dietler, P. L. Stanghellini et al., “Structural and infrared spectroscopic characterization of Co_{6}C(CO)_{12}S_{2}: a highnuclearity carbido carbonyl cluster spontaneously formed from dicobalt octacarbonyl and carbon disulphide,” Journal of Organometallic Chemistry, vol. 213, no. 1, pp. 277–292, 1981. View at: Publisher Site  Google Scholar
 P. L. Stanghellini, L. Cognolato, G. Bor, and S. F. A. Kettle, “Identification of the IR frequencies of the interstitial carbon atom in some M_{6}C (M = Fe, Ru) carbonyl clusters by^{13} Cenrichment,” Journal of Crystallographic and Spectroscopic Research, vol. 13, no. 2, pp. 127–133, 1983. View at: Publisher Site  Google Scholar
 D. Cremer and J. Grafenstein, “Calculation and analysis of NMR spin–spin coupling constants,” Physical Chemistry Chemical Physics, vol. 9, pp. 2791–2816, 2007. View at: Publisher Site  Google Scholar
 R. Gaur and L. Mishra, “Synthesis and characterization of Ru(II)–DMSO–Cl–chalcone complexes: DNA binding, nuclease, and topoisomerase II inhibitory activity,” Inorganic Chemistry, vol. 51, pp. 3059–3070, 2012. View at: Publisher Site  Google Scholar
 J. Pratuangdejkul, P. Jaudon, O. C. Ducrocq et al., “Cationπ interactions in serotonin: conformational, electronic distribution, and energy decomposition analysis,” Journal of Chemical Theory and Computation, vol. 2, pp. 746–760, 2006. View at: Publisher Site  Google Scholar
 T. Helgaker and M. Jaszunski, “Ab initio methods for the calculation of NMR shielding and indirect spinspin coupling constants,” Chemical Reviews, vol. 99, pp. 293–352, 1999. View at: Publisher Site  Google Scholar
 T. Yoshizawa and M. Hada, “Relativistic quantumchemical calculations of magnetizabilities of noble gas atoms using the Douglas Kroll Hess method,” Chemical Physics Letters, vol. 458, pp. 223–226, 2008. View at: Publisher Site  Google Scholar
 C. ForoutanNejad, “Interatomic magnetizability: a QTAIMbased approach toward deciphering magnetic aromaticity,” The Journal of Physical Chemistry A, vol. 115, pp. 12555–12560, 2011. View at: Publisher Site  Google Scholar
 T. Mineva, E. Sicilia, and N. Russo, “Densityfunctional approach to hardness evaluation and its use in the study of the maximum hardness principle,” Journal of the American Chemical Society, vol. 120, no. 35, pp. 9053–9058, 1998. View at: Google Scholar
 G. de Luca, E. Sicilia, N. Russo, and T. Mineva, “On the hardness evaluation in solvent for neutral and charged systems,” Journal of the American Chemical Society, vol. 124, no. 7, pp. 1494–1499, 2002. View at: Publisher Site  Google Scholar
 T. Mineva, V. Parvanov, I. Petrov, N. Neshev, and N. Russo, “Fukui indices from perturbed Kohn–Sham orbitals and regional softness from Mayer atomic valences,” The Journal of Physical Chemistry A, vol. 105, no. 10, pp. 1959–1967, 2001. View at: Publisher Site  Google Scholar
 R. G. Pearson and J. Chem, “Recent advances in the concept of hard and soft acids and bases,” Journal of Chemical Education, vol. 64, pp. 561–567, 1987. View at: Publisher Site  Google Scholar
 A. A. Adeniyi and P. A. Ajibade, “Theoretical study of the electronic and spectroscopic properties of some Ru(II) anticancer complexes,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 105, pp. 456–465, 2013. View at: Publisher Site  Google Scholar
 P. J. Mendesa, J. P. P. Ramalho, A. J. E. Candeias, M. P. Robalod, and M. H. Garcia, “Density functional theory calculations on η^{5}monocyclopentadienylnitrilecobalt complexes concerning their secondorder nonlinear optical properties,” Journal of Molecular Structure: THEOCHEM, vol. 729, pp. 109–113, 2005. View at: Publisher Site  Google Scholar
 I. C. de Silva, R. M. de Silva, and K. M. N. de Silva, “Investigations of nonlinear optical (NLO) properties of Fe, Ru and Os organometallic complexes using high accuracy density functional theory (DFT) calculations,” Journal of Molecular Structure: THEOCHEM, vol. 728, no. 1–3, pp. 141–145, 2005. View at: Publisher Site  Google Scholar
 E. E. Ebenso, D. A. Isabirye, and N. O. Eddy, “Adsorption and quantum chemical studies on the inhibition potentials of some thiosemicarbazides for the corrosion of mild steel in acidic medium,” International Journal of Molecular Sciences, vol. 11, pp. 2473–2498, 2010. View at: Google Scholar
Copyright
Copyright © 2014 Adebayo A. Adeniyi and Peter A. Ajibade. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.