Abstract

This study examines the influence of the number of nitrogen atoms in azo ligands on the geometric, thermodynamic, electronic, and spectroscopic properties of the trans isomers (γ and δ) of ruthenium (II) complexes of the formula RuCl₂L₂ (L = azben, azpy, or papm). The calculations were performed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT) at the B3LYP/Lanl2DZ level to determine the optimized structures, thermodynamic parameters, reactivity descriptors, electronic interactions, and electronic absorption properties. The results show that substituting pyridine with pyrimidine has little effect on the geometry of the complexes. However, the substitution of pyridine with phenyl, in the case of the γ-RuCl₂(azben)₂ isomer, profoundly alters the complex’s octahedral structure. Also, it significantly alters their stability and reactivity. At 298.15 K, the formation of the RuCl₂(papm)₂ and RuCl₂(azpy)₂ isomers is spontaneous and exothermic, whereas the formation of the RuCl₂(azben)₂ isomers is thermodynamically unfavorable under the same conditions. The RuCl₂(azpy)₂ complexes appear to be the most thermodynamically stable. Analysis of frontier orbitals and global reactivity descriptors indicates that replacing azopyridine with azobenzene increases the complexes’ electronic reactivity. In particular, the γ-RuCl₂(azben)₂ isomer exhibits the greatest affinity for interactions with DNA bases, making it a promising candidate for anticancer applications. Furthermore, TD-DFT calculations reveal that all complexes absorb in the visible spectrum via MLCT (Metal-to-Ligand Charge Transfer) transitions. These optical properties suggest that these compounds could serve as promising photosensitizers for dynamic cancer phototherapy.

Keywords

Ruthenium Complexes, Azo Ligands, DFT, TD-DFT, Electronic Reactivity, MLCT, Dynamic Phototherapy, Cancer

Share and Cite:

1. Introduction

Compounds containing the azo group (-N=N-) constitute an important family of chromophores with numerous applications in the food, cosmetics, dyes, textile, solar energy, and pharmaceutical industries [1] [2]. The interest in this chemical function stems from the physicochemical properties and biological activities that result from its combination with various aromatic or heteroaromatic rings, such as benzene, pyridine, or pyrimidine. Thus, azo derivatives represent a promising class of bioactive compounds with antifungal, anti-inflammatory, antitubercular, and anticancer properties [3].

Among these heterocycles, pyridine and pyrimidine play an important role in medicinal chemistry due to their strong ability to modulate the electronic and biological properties of the molecules that contain them. The bonding of the azo group (-N₁=N₂-) to benzene via the N₁ atom and, respectively, to benzene, pyridine, or pyrimidine via the N₂ atom results in the ligands azobenzene (azben), 2-phenylazopyridine (azpy) [4], and 2-phenylazopyrimidine (papm) [5]. The corresponding molecular structures are shown in Figure 1.

Figure 1. 2D structure of the azo ligands studied.

These three ligands differ primarily in the number of nitrogen atoms present in the aromatic ring bonded to the N₂ atom of the azo group. This structural difference is likely to significantly influence the electronic distribution, the coordination capacity, and, consequently, the physicochemical and biological properties of the compounds. Indeed, nitrogen atoms possess non-bonding electron pairs that can alter the ligand’s electron density, its affinity for metal centers, as well as its solubility, bioavailability, metabolism, and interaction with biological targets [6]. Many therapeutic molecules, in fact, exploit this structural feature. Azo ligands also constitute an important family of bidentate ligands capable of stabilizing ruthenium ions at oxidation states +II and +III [7]. Azo ligands are also an important class of bidentate ligands capable of stabilizing ruthenium ions in oxidation states +II and +III [7]. All RuCl₂L₂ complexes were treated as neutral Ru(II) species containing two anionic chloride ligands and two neutral azo ligands. The calculations were carried out using a singlet ground state (multiplicity = 1). The total charge and multiplicity of each molecular system were explicitly specified in the computational input files. If applicable: Alternative spin states were also examined and found to be higher in energy than the singlet state. The resulting RuCl₂L₂ complexes (L = azo ligand) exhibit particularly interesting electronic, photochemical, and biological properties, which are exploited in both dye-sensitized solar cells and the development of new anticancer agents [7]-[9]. In particular, several experimental and theoretical studies have shown that ruthenium complexes derived from 2-phenylazopyridine used as reference systems in this work exhibit promising activity against various types of cancer, notably colon, breast, and throat cancers [10] [11].

The complexation of ruthenium derived from the RuCl₃∙3H₂O precursor by azopyridine ligands leads to the formation of five geometric isomers designated α-Cl, β-Cl, γ-Cl, δ-Cl, and ε-Cl [12] [13]. The α-Cl, β-Cl, and ε-Cl isomers belong to the cis configuration, while the γ-Cl and δ-Cl isomers exhibit a trans configuration, defined by the relative arrangement of the two chloride ligands around the metal center. Previous studies have shown that trans isomers generally possess greater stability and more promising biological properties. For this reason, only the γ-Cl and δ-Cl isomers will be considered in this study.

The main objective of this work is to evaluate, using density functional theory (DFT), the effect of replacing the pyridine ring in azopyridine with a benzene or pyrimidine ring on the geometric, electronic, and spectroscopic properties of the γ-RuCl₂(azpy)₂ and δ-RuCl₂(azpy)₂ complexes. More specifically, this study aims to highlight the influence of the number of nitrogen atoms present in the aromatic ring of the azo ligand on the stability, reactivity, and optical properties of the ruthenium complexes.

To achieve this objective, density functional theory (DFT) is used as the computational method. Recognized for its excellent balance between accuracy and computational cost, this approach is now one of the most powerful tools for modeling molecular and organometallic systems [14]. Based on solving the Schrödinger equation using the electron density, DFT allows for the determination, with good accuracy of the optimized structures, electronic properties, reactivity descriptors, and spectroscopic properties of the compounds under study.

2. Materials and Methods

2.1. Geometric Optimization and Frequency Analysis

The molecular systems studied in this work are obtained by substituting the pyridine ring of the γ-Cl and δ-Cl isomers of the reference complex RuCl₂(azpy)₂ with a benzene or pyrimidine ring. This strategy led to a comparative study of six ruthenium complexes, with the aim of evaluating the influence of the number of nitrogen atoms present in the ligand’s aromatic ring on their structural, electronic, and spectroscopic properties.

Quantum chemistry calculations are performed within the framework of density functional theory (DFT), using the Becke three-parameter hybrid functional combined with the Lee-Yang-Parr correlation functional (B3LYP) [15]. This choice is motivated by the fact that B3LYP offers an excellent balance between accuracy and computational cost. This functional has been extensively validated for the study of transition metal complexes, particularly those of ruthenium, and satisfactorily reproduces the geometric parameters, electronic properties, and relative energies of the various isomers.

Geometric optimizations and vibrational frequency calculations are performed using the B3LYP/LanL2DZ theory [16]. The choice of the LanL2DZ basis set is justified by the presence of the ruthenium atom, a transition metal for which this basis set incorporates an Effective Core Potential (ECP). This approach allows for the relativistic effects of inner-shell electrons to be taken into account while significantly reducing computational cost, without compromising the quality of the results.

All calculations were performed using the Gaussian software [17]. Following geometric optimization, a frequency calculation was systematically performed to verify the nature of the obtained structures. The absence of imaginary frequencies confirms that each optimized structure corresponds to a local minimum on the potential energy surface, i.e., a stable ground state.

The optimized structures then serve as the basis for determining the main geometric parameters (bond lengths and angles around the metal center), thermodynamic quantities (enthalpy, entropy, and Gibbs free energy), as well as electronic properties such as the energies of the boundary orbitals (HOMO and LUMO). These parameters enable the analysis of the complexes’ stability, their chemical reactivity, and their potential ability to interact with biological systems or exhibit photochemical properties.

2.2. Thermodynamic Quantities of Reactions

The complexes under study are assumed to be obtained experimentally by the reaction between ruthenium trichloride trihydrate, RuCl₃∙3H₂O, and two azo ligands L (L = azpy, azben, or papm), according to the following overall reaction:

${\text{RuCl}}_3\cdot 3{\text{H}}_2\text{O}+2\text{L}\to {\text{RuCl}}_2{\text{L}}_2+\frac{1}{2}{\text{Cl}}_2+3{\text{H}}_2\text{O}$

The thermodynamic parameters associated with the formation reactions were calculated from the electronic energies obtained at the optimized geometries, including zero-point energy and thermal corrections derived from frequency calculations at 298.15 K. The Gibbs free energies were subsequently obtained by combining the corresponding enthalpic and entropic contributions.

All thermodynamic quantities were calculated for the gas phase at 298.15 K and under the standard-state convention adopted in the calculations. To evaluate the thermodynamic feasibility of this reaction, the standard thermodynamic properties of the reactants and products are determined based on calculations performed using the Gaussian software. The values for the standard enthalpy of formation (Δ_fH˚), the standard Gibbs free energy of formation (Δ_fG˚), and the standard entropy of formation (Δ_fS˚) are extracted from the calculation output files.

These data are then used to calculate the changes in reaction enthalpy (Δ_rH˚), Gibbs free energy of reaction (Δ_rG˚), and reaction entropy (Δ_rS˚), in accordance with Hess’s law. These parameters allow us to assess whether the reaction is endothermic or exothermic, its spontaneity at the given temperature, and the evolution of disorder during the formation of the ruthenium complexes under study.

2.3. Global Reactivity Descriptors

Boundary orbitals play a crucial role in interpreting the chemical reactivity of molecular systems [18]. In particular, the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) provide essential information on a molecule’s ability to donate or accept electrons [19]. These parameters thus provide a relevant basis for evaluating the stability and reactivity of the complexes under study.

In this work, several global reactivity descriptors, derived from the hard and soft acid-base (HSAB) theory developed by Pearson [20], are calculated and analyzed. These include the energy gap between the HOMO and LUMO orbitals (energy gap, ΔE), the chemical potential (μ), the chemical hardness (η), and the electrophilicity index (ω) [21]. These descriptors are used to characterize the electronic stability of the complexes, their ability to transfer charge, and their behavior with respect to nucleophilic or electrophilic species.

The values of these various descriptors are determined from the energies of the boundary orbitals according to the following relationships:

$\Delta \text{E}={\text{E}}_{\text{LUMO}}-{\text{E}}_{\text{HOMO}}$ ; $\mu =- \frac{\text{PI}+\text{AE}}{2}$ ; $\eta =\left(\text{PI}-\text{AE}\right)$ ; $\omega = \frac{{\mu }^2}{2\eta }$

where $\text{AE}=-{\text{E}}_{\text{LUMO}}$ is the electron affinity, which characterizes a molecule’s susceptibility to nucleophilic attack, and $\text{PI}=-{\text{E}}_{\text{HOMO}}$ is the ionization potential, which characterizes a molecule’s susceptibility to electrophilic attack.

2.4. Analysis of the Natural NBO Population

The Natural Bond Orbital (NBO) method is a widely used post-processing analysis tool for interpreting the results of quantum chemistry calculations [22]. It describes the electronic structure of molecules using localized orbitals, which are similar to the classical representation of chemical bonds and non-bonding doublets.

This method highlights donor-acceptor interactions between occupied and unoccupied orbitals, which are responsible for phenomena such as electronic delocalization, hyperconjugation, and charge transfer. NBO analysis also provides the natural populations, atomic charges, and stabilization energies associated with these interactions, thereby offering a detailed description of the nature of chemical bonds.

In this work, NBO analysis is used to characterize the electronic interactions within the studied ruthenium complexes, to evaluate the strength of metal-ligand bonds, and to highlight the effect of azo ligand substitution on the electronic distribution and stability of the complexes. It thus complements the geometric, energetic, and orbital analyses obtained by density functional theory (DFT), providing a more refined chemical interpretation of the results.

2.5. Method for Obtaining UV-Visible Spectra

Time-Dependent Density Functional Theory (TD-DFT) is an extension of Density Functional Theory (DFT) to systems subject to time-dependent perturbations. Based on the Runge-Gross theorem, it establishes that the time-dependent electron density uniquely determines the external potential applied to the system. This approach thus allows for the description of excited electronic states and transitions between different energy levels.

In practice, TD-DFT relies on solving the time-dependent Kohn-Sham equations, in which electrons evolve under the influence of an external perturbation, such as an electromagnetic field. This method is now widely used in quantum chemistry to predict excitation energies, absorption wavelengths, oscillator forces, and the optical properties of molecules and coordination complexes [23].

In this work, the electronic properties of excited states are studied using the TD-B3LYP/LanL2DZ theory, implemented in the Gaussian software package. The calculations allow us to determine the electronic transition energies, absorption wavelengths, oscillator forces, and the nature of the transitions involved. These parameters are then used to simulate and analyze the UV-visible absorption spectra of the ruthenium complexes under study, in order to evaluate the effect of azo ligand substitution on their optical and photochemical properties.

3. Results and Discussion

3.1. Comparison of the Geometry of RuCl₂L₂ Complexes

The structures of the compounds studied are derived from those of the γ-RuCl₂(azpy)₂ and δ-RuCl₂(azpy)₂ isomers by replacing pyridine with pyrimidine or benzene. These compounds differ in the number of nitrogen atoms. Pyridine and pyrimidine contain one and two, respectively. As for benzene, its structure contains no nitrogen atoms. Figure 2 shows the ground-state structures of the γ-RuCl₂L₂ and δ-RuCl₂L₂ isomers (L = azben, azpy, or papm). These structures were obtained following geometric optimization and frequency calculations. Furthermore, these minimum-energy geometries are characterized by the absence of imaginary frequencies. Table 1 summarizes the geometric parameters of the compounds in this minimum-energy state.

Figure 2. Molecular geometries of the complexes calculated at the B3LYP/Lanl2DZ level.

Table 1. Bond lengths (in Å) and bond angles (in ˚) involving ruthenium, determined using the B3LYP/Lanl2DZ theory [6] [13].

Analysis of the geometric data shows that the calculated bond lengths and bond angles are very close to those obtained experimentally. The geometry of the compounds is therefore accurately reproduced by the theoretical calculations. The theoretical approach used is thus suitable for predicting the structure of this type of ruthenium complex. The bond lengths between ruthenium and the two ligand binding sites (N2, Np, or C) range from 2.00 Å to 2.11 Å. In general, the substitution of pyrimidine for pyridine has little effect on the geometry of the ruthenium azopyridine complexes. However, the substitution of pyridine with phenyl, in the case of the γ-RuCl₂(azben)₂ isomer, profoundly alters the complex’s octahedral structure (Np-Ru-Np angle). This difference can be explained, on the one hand, by the Yann-Teller effect and, on the other hand, by repulsive interactions between the free electron pairs on the nitrogen atoms of pyridine and pyrimidine.

3.2. Enthalpy and Spontaneity of the Reaction Leading to the Formation of the Studied Complexes

The spontaneity and enthalpy of a chemical reaction involving the isomers of the studied complexes are quantified, respectively, by determining the change in reaction enthalpy (Δ_rH˚) and the change in reaction free energy (Δ_rG˚). A positive value of Δ_rH˚ and a positive value of Δ_rG˚ indicate, respectively, that the formation reaction under study is endothermic and non-spontaneous at the temperature of study. Conversely, negative changes in these quantities allow these formation reactions to be characterized as exothermic and spontaneous. In addition to these two quantities, the change in entropy (Δ_rS˚) accompanying the formation reaction is also considered. A positive value of Δ_rS˚ indicates an increase in disorder, whereas a negative value indicates a decrease in disorder. Table 2 presents the values for the change in enthalpy and the change in free enthalpy of the trans isomers of the ruthenium complexes studied.

Table 2. Values of thermodynamic reaction parameters.

The values of the enthalpy changes and free enthalpy for the isomers of the RuCl₂(papm)₂ and RuCl₂(azpy)₂ complexes are all negative, indicating spontaneous and exothermic formation reactions. In the case of the RuCl₂(azben)₂ isomers, the calculated changes are all positive, reflecting non-spontaneous and endothermic reactions at 298 K. Furthermore, examination of the various ∆_rG˚ values in Table 2 allows us to establish the following ascending order of reaction spontaneity:

δ-RuCl₂(azben)₂ ˂ γ-RuCl₂(azben)₂ ˂ γ-RuCl₂(papm)₂ ˂ δ-RuCl₂(papm)₂ ˂ γ-RuCl₂(azpy)₂ ˂ δ-RuCl₂(azpy)₂

This ranking highlights the greater stability of the isomers of ruthenium azopyridine complexes RuCl₂(azpy)₂. With regard to the entropy change ∆_rS˚, the recorded values range from 9.88 cal∙mol⁻¹∙K⁻¹ to 15.54 cal∙mol⁻¹∙K⁻¹. These positive values of ∆_rS˚ indicate an increase in disorder during the formation of the compounds studied.

3.3. Comparison of the Overall Reactivity of the Studied Complex Isomers

Frontier molecular orbitals (FMOs) play a key role in the chemical reactivity of compounds. The molecular orbitals of interest are the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). The distribution and energies of these various frontier orbitals for the studied compounds are illustrated and compared in Figure 3. Additionally, Table 3 provides information on the percentage contribution of the various atomic entities (ruthenium, chlorine, and ligand) to the formation of the frontier orbitals.

Figure 3. Energy levels and distribution of the frontier orbitals of the compounds studied.

Table 3. Composition of the frontier orbitals of the various complexes.

Analysis of the distribution of the isodensity surface and the various contributions reveals that ruthenium’s contribution to the formation of the compounds’ HOMO ranges from 30% to 62%. This significant contribution of ruthenium’s d orbitals to the formation of the HOMO shows that the compounds’ HOMO is predominantly localized on ruthenium. It is therefore characterized by the ruthenium d orbitals. As for the LUMO of the compounds containing the azpy and papm ligands, it consists of between 1% and 2% of the ruthenium orbitals. The LUMO of these compounds is therefore localized on the ligand. However, the LUMO of the γ and δ isomers of RuCl₂(azben)₂ shows contributions from ruthenium’s d orbitals of 53% and 27%, respectively. The LUMO of γ-RuCl₂(azben)₂ has a metallic character. Thus, for this isomer, ligand-metal transitions (LMCTs) could be observed between its HOMO and LUMO. However, for this compound, MLCT transitions are observed during the electron promotion of the HOMO-1 and HOMO-2 orbitals to the LUMO+1 and LUMO+2 orbitals. The composition of the HOMO and LUMO boundary orbitals of the other compounds favors metal-ligand charge transfer (MLCT) transitions of the t_2g → π* type.

The energies of these orbitals are listed in the table. They are used to determine the values of the various global reactivity parameters.

Table 4. Global reactivity parameters (in eV).

The analysis of the global reactivity descriptors highlights a significant influence of the nature of the ligand and the isomerism on the electronic properties of the studied complexes. Table 4 shows the value of these parameters. In general, a small HOMO-LUMO energy gap (ΔE) and low chemical hardness (η) are associated with a greater ease of electron-density redistribution and, consequently, enhanced reactivity. Conversely, high values of ΔE and η indicate greater resistance to electronic perturbations and are generally characteristic of more stable systems. Among the studied complexes, γ-RuCl₂(azben)₂ appears to be the most reactive, exhibiting the smallest energy gap (ΔE = 1.19 eV) and the lowest chemical hardness (η = 0.59 eV). This high reactivity is further supported by its high electrophilicity index (ω = 19.67 eV), which is the largest value in the series. In contrast, δ-RuCl₂(azben)₂ exhibits the largest energy gap (ΔE = 2.78 eV) and the highest chemical hardness (η = 1.39 eV), indicating greater electronic stability and lower reactivity. The other complexes display intermediate behavior. Thus, the substitution of pyridine by benzene in the azben ligand leads to electronic properties that strongly depend on the isomeric form: the γ isomer is characterized by enhanced reactivity, whereas the δ isomer exhibits significantly greater electronic stability. These results suggest that γ-RuCl₂(azben)₂ has a particular propensity to participate in interactions involving charge transfer, whereas δ-RuCl₂(azben)₂ is distinguished by its greater electronic stability, particularly associated with its high chemical hardness.

By mainly considering ΔE, η, and ω, the following ranking can be proposed:

Increasing reactivity:

δ-RuCl₂(azben)₂ < δ-RuCl₂(papm)₂ < γ-RuCl₂(azpy)₂ < γ-RuCl₂(papm)₂ < δ-RuCl₂(azpy)₂ < γ-RuCl₂(azben)₂

Conversely, in terms of electronic stability:

γ-RuCl₂(azben)₂ < δ-RuCl₂(azpy)₂ < γ-RuCl₂(papm)₂ < γ-RuCl₂(azpy)₂ < δ-RuCl₂(papm)₂ < δ-RuCl₂(azben)₂

3.4. Comparison of Electrostatic Potential Distribution

Electrostatic potential maps show the three-dimensional distribution of charges within molecules and allow us to visualize the distribution of charges in a molecule. Understanding this charge distribution helps determine how the molecule interacts with another chemical system. Figure 4 shows the charge distribution within the various complexes.

A negative electrostatic potential corresponds to an electrophilic site (proton attraction) due to a concentration of electron density (red color on the ESP surface). Conversely, a positive electrostatic potential corresponds to a nucleophilic site (proton repulsion) due to low electron density (shown in blue on the ESP surface). The electrostatic potential increases in the following order of colors: red - orange - yellow - green - blue.

Examination of the surfaces shows that the regions containing the two chlorine atoms and the ruthenium, on the one hand, and the nitrogen atoms not bonded to ruthenium, on the other, are the electrophilic sites. This observation is explained by the presence of free electron pairs on these atoms. The carbon atoms of the various aromatic groups constitute potential nucleophilic attack sites, as indicated by the blue color of these regions.

Figure 4. Molecular electrostatic potential surfaces of the complexes studied.

3.5. NBO Analysis of the Structure of the Studied Complexes

The electronic density of the studied complexes is analyzed using the Natural Bond Orbital (NBO) method. The NBO second-order perturbation analysis was carefully reassessed. Only donor-acceptor interactions involving occupied donor orbitals and low-lying unoccupied acceptor orbitals were considered in the discussion. The corresponding stabilization energies were used to identify the main charge-transfer and hyperconjugative interactions contributing to complex stabilization. This provides information on the distribution of electrons among the various atomic orbitals that make up the chemical system. Table 5 illustrates the electronic configuration of the central metal within the studied ruthenium complex isomers. The electron population of ruthenium ranges from 44.00843 to 44.11269 electrons. This deviation from the electron population of the ruthenium atom on its own (44 electrons) is explained by the delocalization of electrons from the nitrogen atoms of the ligand to ruthenium during complex formation. This delocalization is more pronounced in the RuCl₂(azben)₂ complexes. These electrons are distributed equally among the core orbitals, valence orbitals, and Rydberg orbitals. In fact, in these compounds, 36 electrons do not participate in the chemical reaction because they occupy the core orbitals, and eight electrons make up the orbitals that may be involved in an interaction or reaction with another chemical system. As for the virtual orbitals, their electron population ranges from 0.03328 to 0.04741. This low population of Rydberg electrons highlights the accuracy and rigor of the theoretical framework used in the distribution of electrons among the various molecular orbitals and, in particular, atomic orbitals. These different orbitals are described using the same atomic orbitals of ruthenium.

Table 5. Electronic configuration of the ruthenium nucleus.

NBO analysis makes it possible to assess the intramolecular interactions that contribute to the stabilization of the molecular structure. This contribution to stability is measured using the second-order perturbation energy E². Furthermore, a higher value of the perturbation energy indicates a greater contribution to the stabilization of the compound. Table 6 presents the values of the second-order perturbation energies associated with stabilizing interactions involving ruthenium.

Table 6. Interactions involving ruthenium and associated second-order perturbation energy values in kcal∙mol⁻¹.

An examination of the E² energy values for the RuCl₂(papm)₂ isomers indicates that the interactions between the lone pair LP(1)N and the antibonding LP*(4) Ru orbital are associated with the highest second-order perturbation energy values (99.74 kcal∙mol⁻¹, 99.67 kcal∙mol⁻¹, and 94.06 kcal∙mol⁻¹). The recorded E² values show that, regardless of the ligand, the interactions between the ligand’s bonds or free doublets and the ruthenium orbitals have higher energy. They reflect electron donation from the ligand (papm, azpy, or azben) to ruthenium. Indeed, within the complexes, the delocalization of the π electrons from the azo group has the effect of stabilizing the complex. This greater energy transfer is indicative of a ligand-metal transition (LMCT) at low wavelengths. However, stabilizing interactions between the ruthenium d orbitals and the π*(N=N) orbitals of the azo bond or the free doublets of the chlorine atoms are observed. These interactions result in electron back-donation. The papm, azpy, and azben ligands are therefore electron acceptors. This delocalization of ruthenium’s electron density—which is less significant—stabilizes the resulting complex. This type of interaction is responsible for electron transfers from ruthenium to the ligand (MLTC), which can be exploited in photodynamic therapy.

3.6. Effect of Substitution on Spectroscopic Properties

Azo compounds are known for their photochromic properties. The spectroscopic properties of these compounds were studied in vacuum and in acetonitrile using time-dependent DFT at the TD-B3LYP/Lanl2DZ level. Figure 5 compares the absorption spectra obtained in the solvent.

Figure 5. UV-visible spectra in acetonitrile of the various complexes studied.

Table 7. Characteristic parameters (excitation energy E_exc, wavelength λ, and oscillator strength f) of the UV-visible spectra of the complexes studied [6].

The wavelength values obtained experimentally and those obtained by calculation are very close (as shown in Table 7). The B3LYP functional and the Lanl2DZ pseudopotential used are therefore appropriate for simulating and predicting the UV-visible absorption spectra of ruthenium complexes. The interaction of these compounds with light gives rise to two types of electronic transitions. The first type of transition, observed at wavelengths below 400 nm, consists of intra-ligand transitions (LLCTs). The transitions observed between 400 nm and 800 nm are metal-ligand transitions (MLCTs). This latter type of transition is responsible for the photochemical properties of ruthenium complexes.

The spectra of the δ-RuCl₂(papm)₂ and δ-RuCl₂(azpy)₂ isomers exhibit a very intense band centered at 400 nm and a second, much weaker band centered at 700 nm. The transitions associated with these bands are of the intra-ligand LLCT and metal-ligand MLCT types, respectively. The substitution of pyrimidine for pyridine results in a slight bathochromic effect (shift toward longer wavelengths) accompanied by a more or less significant hypochromic effect (decrease in absorption intensity). In the case of the isomers γ-RuCl₂(papm)₂ and γ-RuCl₂(azpy)₂, the absorption spectra have the same profile and exhibit a single intense band centered at 589.674 nm and 463.634 nm, respectively. These bands reveal MLCT transitions with an absorption intensity greater than that of the δ-RuCl₂L₂ isomers. Additionally, shoulders are observed on either side of the intense band. The left shoulder, located at a wavelength below 400 nm, is characteristic of an intra-ligand transition. Furthermore, the substitution of pyrimidine for pyridine also produces a slight bathochromic and hypsochromic effect. However, the substitution of pyridine with benzene affects the shape of the complex’s absorption spectrum, which now exhibits only a single intense absorption band centered at 382.209 nm for δ-RuCl₂(azben)₂ and at 436.905 nm for γ-RuCl₂(azben)₂. These bands are characteristic of intra-ligand transitions. Furthermore, metal-ligand transitions in these compounds, observed at 450.003 nm for γ-RuCl₂(azben)₂ and at 519.230 nm for δ-RuCl₂(azben)₂, show that cyclomethylation has a significant hypsochromic effect (shifts toward shorter wavelengths) on this transition. Overall, the results reveal a structure-property relationship in which the nature of the aromatic ring and the geometric arrangement of the ligands around ruthenium significantly modulate the absorption properties. The MLCT transitions located in the visible range especially pronounced for certain γ isomers give these complexes potential interest for photophysical and photochemical applications.

4. Conclusions

The objective of this work was to evaluate the influence of the nature of the azo ligand’s aromatic ring on the geometric, thermodynamic, electronic, and spectroscopic properties of ruthenium complexes of the general formula RuCl₂L₂ (L = papm, azpy, or azben), using density functional theory (DFT) and its time-dependent extension (TD-DFT) at the B3LYP/LanL2DZ theory level.

The results show that replacing the ligand’s aromatic ring has little effect on the coordination geometry around the metal center. The characteristic bond lengths and angles of ruthenium remain largely unchanged, reflecting the robustness of the octahedral environment of the complexes studied.

The thermodynamic study reveals that the reactions leading to the formation of isomers of the RuCl₂(papm)₂ and RuCl₂(azpy)₂ complexes are spontaneous and exothermic at 298.15 K, unlike those of the RuCl₂(azben)₂ complexes, which become favorable only at higher temperatures. The azopyridine-derived complexes thus appear to be the most thermodynamically stable.

Analysis of frontier orbitals and global reactivity descriptors reveals a significant influence of the ligand’s nature on the electronic properties. Replacing azopyridine with azobenzene increases the complexes’ electronic reactivity and enhances their ability to transfer charge. Among the compounds studied, the γ-RuCl₂(azben)₂ isomer exhibits the electronic characteristics most conducive to interactions with DNA bases, suggesting interesting biological potential.

The spectroscopic study conducted using TD-DFT shows that all of the complexes absorb in the visible region. The observed electronic transitions are primarily of the Metal-to-Ligand Charge Transfer (MLCT) type, characteristic of ruthenium complexes. These optical properties confirm their potential utility as photosensitizers in dynamic phototherapy applications for cancer treatment.

Overall, this study highlights that modifying the nature of the ligand’s aromatic ring is a relevant strategy for modulating the electronic properties, reactivity, and optical properties of ruthenium complexes. The results thus provide useful insights for the rational design of new organometallic complexes intended for biomedical and photochemical applications.

Looking ahead, it would be interesting to:

• study the interactions of these complexes with biomolecules of therapeutic interest (DNA, target proteins, or enzymes) using molecular docking and molecular dynamics approaches to better understand their mechanisms of action;

• to examine the effect of increasing the number of nitrogen atoms or extending the conjugation of the ligands on the electronic, spectroscopic, and biological properties of the complexes;

• to explicitly evaluate the photophysical properties (fluorescence yield, excited-state lifetimes, singlet oxygen production) in order to assess their potential in dynamic phototherapy;

• finally, to compare theoretical predictions with experimental results from synthesis, characterization, and biological evaluation in order to validate the trends revealed by quantum calculations.

Author Contributions

N’Guessan Kouakou Nobel: Initiating Author; Kangah Niamké Jean Baptiste: Initiating Author; Koné Mamadou Guy-Richard: Initiating Author.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References