DFT Study of Molecular and Electronic Structure of Y, La and Lu Complexes with Porphyrazine and Tetrakis(1,2,5-thiadiazole)porphyrazine

Abstract

Electronic and geometric structures of Y, La and Lu complexes with porphyrazine (Pz) and tetrakis(1,2,5-thiadiazole)porphyrazine (TTDPz) were investigated by density functional theory (DFT) calculations and compared. The nature of the bonds between metal atoms and nitrogen atoms has been described using the analysis of the electron density distribution in the frame of Bader’s quantum theory of atoms in molecule (QTAIM). Simulation and interpretation of electronic spectra were performed with use of time-dependent density functional theory (TDDFT) calculations. Description of calculated IR spectra was carried out based on the analysis of the distribution of the potential energy of normal vibrations by natural vibrational coordinates.

1. Introduction

Macroheterocycles such as porphyrines, phthalocyanines and their analogues have found a number of applications [1,2,3,4]. The possibility of modification of peripheral substituents [5] or atoms in a central ring allows fine-tuning or modulation of their properties [5,6] such as stability, rigidity, extended conjugated π-electron cloud, conductivity, light absorption capability, long-lived fluorescence, possibility to host uncoupled electrons on the central metal ion or in the aromatic cloud, and so on [7].

Substitution of benzene rings in phthalocyanines by aromatic heterocycles, e.g., by pyrazine [8,9] or 1,2,5-chalcogenadiazole [10] have a strong impact of the electronic properties of the central porphyrazine (Pz) core which is common for these systems.

The presence of a five-membered heterocycle containing nitrogen and sulfur atoms on the periphery of tetrakis(1,2,5-thiadiazole)porphyrazine (TTDPz) considerably modulates the physico-chemical properties of the macrocycle and its ability for intermolecular interaction as compared to phthalocyanines [10]. Unlike phthalocyanines, or their pyrazine-fused analogues, TTDPz have no H atoms on the periphery and their molecular packing during crystal growth is determined by specific N…S interactions [11]. This is quite important for application of these phthalocyanine-type molecules as building blocks for novel functional materials in various fields [12].

Lanthanide complexes of phthalocyanines due to their rich structural diversity and peculiar spectral and electrochemical properties [13] are especially interesting as components of different devices [14]. Among the metal complexes of TTDPz, macrocycle complexes of p- and d-metals are known [MTTDPz]: M = Li [15], Mg [16], Ca [17], AlCl, GaCl, InOAc [18], Zn, Cu, Ni, Co, Fe, Mn [11,16,19]. Recently we have also obtained first representatives of the rare earth metal complexes [(acac)YTTDPz], [(acac)LuTTDPz] [20]. Among lanthanide complexes of unsubstituted porphyrazine [(acac)MPz] complexes with dysprosium, neodium and europium have been reported [21,22].

The optical and coordination properties of phthalocyanine-type complexes crucially depend on the electronic structure of the metal and macrocycle and can often be predicted by quantum-chemical calculations. Such investigations in the case of transition metal complexes s are often nontrivial due to the necessity to account for the multireference character of the wavefunction [23]. However the information about the ground-state properties can be readily obtained for closed-shell species. Density functional theory (DFT) provides a straightforward way to obtain quantitative and qualitative information about structural, electronic and spectral features. It is quite important to investigate the influence of a transition metal and a macrocyclic ligand on the properties of the porphyrine-type complexes, particularly, the peculiarities of chemical bonding and spectral properties. Comparison in the series of atoms of d¹ metals with a different radius (Y, La and Lu) allows the taking into account of only the size of atoms. The change of the periphery of the porphyrazine macrocycle (Pz) by fusion of electron-deficient 1,2,5-thiadiazole rings (TTDPz) allows the determination of the influence of the ligand on the coordination properties of central atoms.

Earlier in our laboratory, the Ca, Fe, Co and Zn complexes with porphyrazines and tetrathiadiazoporphyrazines ligands were investigated by quantum-chemical calculations [17,23]. Preliminary quantum-chemical calculations and interpretation of vibrational spectra were also carried out for complexes of tetra(1,2,5-thiadiazolo)porphyrazine with rare earth elements Y and Lu [20]. The main objective of the present study is to identify the influence of the molecular and electronic structures on the properties of macrocyclic complexes. It is quite important to compare the TTDPz complexes with the corresponding complexes of nonsubstituted porphyrazine in order to reveal the impact of annulation of the 1,2,5-thiadiazole rings. The nature of the chemical bonding between metal atoms and nitrogen atoms has been described using the analysis of the electron density distribution in the framework of Bader’s quantum theory of atoms in molecules (QTAIM) [24]. The lowest excited states were also calculated in order to explain the peculiarities and tendencies observed in the experimental electronic absorption spectra available for the Y and Lu tetra(1,2,5-thiadiazolo)porphyrazine complexes [20]. In addition, the vibrational absorption spectra were analyzed and interpreted.

2. Results and Discussion 2.1. Geometric Structure

Equilibrium structures of the Y, La and Lu macroheterocyclic complexes with porphyrazine (MClPz) and tetra(1,2,5-thiadiazolo)porphyrazine (MClTTDPz) bearing chloride as an axial ligand, were found to possess C_4v symmetry with doming distortion (Figure 1). Force-field calculations yielded no imaginary frequencies, indicating that the optimized configurations correspond to the minima on the potential energy hypersurfaces. The calculated molecular parameters are listed in Table 1.

Analysis of the data in Table 1 shows that there is a significant difference (maximum is about 0.2 Å) in the internuclear distances M-Cl and M-N_p in the series of studied complexes, as well as in the parameters M-X1 and X1-X2 which characterize the distortion of the macrocyclic ligand. The change in the size of the coordination cavity (N_p…N_p)_opp of the ligand is also quite large in MClTTDPz over against MClPz. Additionally there is an increase of M-Cl, M-N_p, (N_p…N_p)_opp and (N_p…N_p)_adj distances in the LuClPz-YClPz-LaClPz series, corresponding to the increase in the ionic radius of metals [25]. A similar situation is observed for MClTTDPz. It should be noted that changes in the analogous parameters of MPz and MTTDPz complexes with a number of transition, alkali and alkali earth metals [15,17,23] are comparable to changes in the parameters of the complexes considered in this work. Obviously, such differences are due to the nature of the metals introduced into the macrocycle cavity. However, the nature of the metal has practically no effect on the structure of the periphery of the ligand, since the internuclear distances N_p-C_α, C_α-N_m, C_α-C_β, C_β-C_β, C_β-H (for MClPz), C_β-N_t, and N_t-S (for MClTTDPz) in the macrocyclic fragment change by no more than 0.002 Å, and the bond angles N_pC_αN_m, C_αN_mC_α, C_αN_pC_α and N_tSN_t (for MClTTDPz) change by no more than 0.5 degrees. Thus, we can conclude that the change in metal affects only the size of the macrocyclic cavity. It should also be noted that the introduction of thiadiazole as a substituent changes the bond lengths of pyrrole moieties by no more than 0.07 Å, and the values of bond angles change by no more than 4 degrees.

The value of lanthanide contraction Δr_Ln, which may be estimated as a difference between r(La-N_p) and r(Lu-N_p) distance values, is 0.202 Å for both MClPz and MClTTDPz. It should be noted that the obtained value of Δr_Ln is close to similar theoretical value in the series of compounds with tridentate ligand (0.208 Å) [26], experimental values in the series of ML₃ with monovalent ligands, such as the series of trichlorides (0.186(8) Å), tribromides (0.185(8) Å) and triiodides (0.193(8) Å) [27,28,29,30,31], as well as in the series of compounds with bidentate ligands, for example, in the series of tris-dipivaloylmethanates (0.186(6) Å) [32,33]. The similarity of Δr_Ln values in lanthanide compounds with three ligands and in complexes considered in the present study indicates that, with no regard of the size of the central atom in the lanthanide series, the geometric structure of the coordination cavity is adjusted to the size of the central atom with the formation of coordination bonds.

2.2. Electron Density Distribution and Chemical Bonding in MClPz and MClTTDPz

Selected net atomic charges and topological parameters of the electron density in bond critical points (BCP) are listed in Table 2. According to the Table 2 the incorporation of thiadiazole rings does not lead to a significant change in the charge of the metal and the ligands (the charge of the macroheterocyclic ligand is calculated as the sum of the net atomic charges on the ligand atoms). The presence of electron-withdrawing nitrogen atoms in the thiadiazole rings induces a shift of electron density from the central cavity and leads to a charge transfer in the row N_t←C_β←C_α. It should be noted that the charges on the atoms of the internal macrocycle (C_α, N_p, N_m) also do not change significantly, but the charge on the carbon atom C_β increases by 0.4 e (Table 2). Similar tendencies in the change of the charges were revealed in the complexes MPz and MTTDPz with transition metals Ca, Fe, Co and Zn [17,23]. The charge on the metal atoms is about +2.2 e, which is comparable to the charge on the Y, La and Lu atoms in the tridentate macroheterocyclic ligand [26]. In the studied complexes the charge of a ligand varies from −1.447 to −1.362 e, and in complexes with transition metals from −1.580 to −1.141 e. In addition, the charges on the atoms of the internal macrocycle (C_α, N_p, N_m) are close to each other in the case of all the complexes considered (the complexes studied in this work and the complexes with transition metals [17,23]). Based on the above-mentioned considerations, it can be concluded that the electron density distribution in the coordination cavity of metal complexes of porphyrazine and tetra(1,2,5-thiadiazole)porphyrazine metal does not change significantly regardless of metal nature and thiadiazole annelation.

The nature of a chemical bond can be determined by the value of the electron density, laplacian ∇²ρ. A positive value of the electron density laplacian ∇²ρ indicates ionic interaction. However, the positive values of ∇²ρ for M-Cl bonds, as well as the corresponding delocalization indices δ(M|Cl) representing the magnitudes of the electron exchange between the basins of the corresponding atoms, allow to argue that these bonds, along with an ionic component (Table 2), possess a significant covalent component (about 0.6). Moreover, it can be noted that the internuclear distances r(M-Cl) decrease by 0.03 Å and the delocalization indices increase by about 0.04 in the MClTTDPz complexes as compared to their MClPz analogues. Values of the parameters ∇²ρ, δ(M|N_p) and q(M|N_p) (Table 2) indicate mainly the ionic M-N_p bond with a noticeable covalent component (value of δ(M|N_p)) comparable to the corresponding value for the Pz and TTDPz complexes with Ca and Zn [17].

2.3. Electronic Absorption Spectra

The electronic absorption spectra of MClPz and MClTTDPz molecules simulated on the basis of time-dependent density functional theory (TDDFT) calculations are shown in Figure 2. The calculated electronic spectra of complexes MClPz and MClTTDPz (M=Y, La, Lu) are quite similar. The main difference in the electronic absorption spectra in the MClPz series is the bathochromic shift of the B_x (λ₂) band. In the case of MClTTDPz molecules, the difference in the electronic spectra appears in the change of positions and relative intensities of the B_x (λ₂) and B_y (λ₃) bands, as well as in the relatively small bathochromic shift of the Q-band (λ₁). In general, the spectra can be described by the model of Gouterman [34,35,36], who suggested electronic absorption spectra of porphyrazine molecules, to possess four intense bands B_x, B_y, Q_x, and Q_y, and three bands (B_x, B_y, and Q) in the case of their metal complexes. Interpretation of the electronic spectra was carried out on the basis of the results of TDDFT calculations. The calculated oscillator strengths (f) for the lowest excited states along with their composition (in terms of one-electron transitions) are given in Table 3 (full list of the most active transitions can be found in Supplementary Materials). Analysis of the data in Table 3 demonstrates that there are three excited electronic states of MClPz molecules characterized by a high oscillator strength. These transitions correspond to three intense bands—B_x, B_y, and Q. The low-lying calculated excited state 1¹E for MClPz complexes is formed predominantly by the electronic transitions from HOMO and HOMO-1 to LUMO, which can be assigned to the Q-band. It should be noted that the compositions of the wave functions and the oscillator strengths corresponding to the electronic state 1¹E are very close to each other for all MClPz (M = Y, La, Lu) complexes. The excited state 4¹E of MClPz is formed by the similar electronic transitions and corresponds to the B_x band, which is shifted to larger wavelengths in the case of LaClPz. Transitions in the higher excited states are predicted to have a highest oscillator strength and can be associated with the B_y band. These excited states (6¹E for LaClPz and 5¹E for YClPz and LuClPz complexes) are composed mainly of the electronic transitions from a₁ and a₂ orbitals to the LUMO e*. The changes of electronic spectra are more significant in complexes of the MClTTDPz (Figure 2, Table 3), where the B_x and B_y bands are different in wavelengths and relative intensities. Thus, it can be concluded that the incorporation of a substituent can enhance the effect of the metal on the spectral properties of the complexes.

2.4. Molecular Orbitals

The energy diagram of molecular orbitals for MClPz and MClTTDPz is shown in Figure 3. The HOMO energies are close in the series LaClPz-YClPz-LuClPz and LaClTTDPz-YClTTDPz-LuClTTDPz, which allows us to conclude that the influence of a metal is insignificant. At the same time, the fusion of thiadiazole rings to the porphyrazine macrocycle leads to a decrease in the HOMO energy, which indicates an increase in the acceptor properties of the compound. In addition, the thiadiazole ring presence significantly reduces the HOMO–LUMO gap from approximately 2.6 to 2.3 eV. A similar situation was observed when considering the Ca(II) and Zn(II) complexes with porphyrazine and tetrakis(1,2,5-thiadiazole)porphyrazine [17] and the HOMO–LUMO gap of magnesium(II) tetra(1,2,5-chalcogenadiazolo)porphyrazines, [TXDPzMg] (X=O, S, Se, Te) is close and varies from ~2.0 to ~2.2 eV [37].

The shapes of molecular orbitals (MOs) that participate in electronic transitions with large oscillator strengths are shown in Figure 4 and Figure 5. Figure 4 and Figure 5 show that LUMO and HOMO are linear combinations of π-bonding and π-antibonding molecular orbitals of the macrocycle, while HOMO-1 and HOMO-6 contain both atomic orbitals (AOs) of the macrocycle and AOs of metal and axial ligand Cl. It should be noted that the shapes of HOMO and LUMO for all considered complexes are similar. However, HOMO exhibits a bonding character for the C_α-C_β bond in the case of MClPz and an antibonding character for MClTTDPz complexes. Differences in the molecular orbitals of the complexes MClPz are observed in the compositions of the HOMO-1 and HOMO-6 MOs (Figure 4). Thus, in the LaClPz molecule, the electron density is shifted towards the Cl atom, while in LuClPz and YClPz the Cl atom has practically no electron density (HOMO-1). In MClTTDPz molecules, the main contribution to transitions to excited electronic states with the highest oscillator strength is made by the LUMO, HOMO, HOMO-2, HOMO-3, and HOMO-5 (HOMO-4 for LaClTTDPz) orbitals. The shapes of these orbitals are shown in Figure 5. A comparative analysis of the selected MOs reveals differences in the shape of 6e orbitals. In the case of LaClTTDPz, the electron density belongs to the macrocycle, while in the La and Y complexes, the orbitals of the chlorine atom contribute. Of further note is the shape of the 4a₁ orbitals, which are a linear combination of bonding and antibonding MOs of the macrocycle and metal. Thus, it can be observed that, in the LaClTTDPz molecule, it includes a bonding π-orbital at the C_β-C_β bond, while these MOs of the LuClTTDPz and YClTTDPz molecules can be considered as completely antibonding.

2.5. Vibrational Spectra

Recently, a description of the absorption bands of the vibrational spectrum of YClTTDPz and LuClTTDPz complexes was published in [20]. Detailed interpretation of IR-spectra for MClTTDPz and MClPz is performed in the present work. The IR spectra were simulated on the basis of the normal mode frequencies and band intensities, which have been calculated by the DFT (B3LYP/pcseg-2) method in a harmonic approximation. Description of the main vibrations is presented in Table 4 (full list of the most active vibrations can be found in Supplementary Materials) and simulated spectra are shown in Figure 6.

In the Y-La-Lu series a high-frequency band shift occurs both in MClTTDPz and in MClPz. The main differences between MClTTDPz and MClPz spectra occur in the positions and characters of the most intensive band, as well as the larger number of intense peaks in MClPz especially noticeable in high-frequency region. The most intensive peak is located in the 982 cm⁻¹ region and corresponds predominantly to C_α-C_β stretching for MClPz complexes. In the case of MClTTDPz the strongest band is observed at ~1290 cm⁻¹ (ω₈₇-ω₈₈) and mainly represented by stretching and bending in the macrocyclic fragment. Spectra of MClPz contain more intensive peaks as compared to MClTTDPz. The largest contribution to the ω₈₅-ω₈₆ bands in MClPz is C_β-C_β stretching. The N_m-C_α stretching vibrations contribute strongly to relatively weak bands ω₈₂-ω₈₃. In the case of the LaClPz ω₇₆ band representing mainly in-plane bending/stretching, it possesses a higher intensity than the corresponding one in LuClPz and YClPz. The out-of-plane deformations provide a major contribution to the medium band ω₄₈. The M-Cl stretching vibrations make a significant contribution to the weak bands in the low frequencies region for MClPz and for MClTTDPz too. MClTTDPz have similar IR-spectra in the 500–1600 cm⁻¹ region where the influence of the metal is almost absent. In this region vibrations of the macrocyclic core are observed. The ω₅₀-ω₅₁ and ω₈₀-ω₈₁ bands in MClTTDPz spectra are composed by N_p-C_α and N_m-C_α stretching vibrations combined with in-plane deformations of the rings. Despite the noted significant covalent character of the bonding of the metal with the macrocycle, the stretching vibrations involving metal are detected only in the low frequencies region or have a low intensity. A performed description of calculated IR spectra shows that the most observed bands cannot be assigned to any single vibrational mode and are the result of the some vibrational modes superpositions. 3. Materials and Methods Computational Details

The DFT-based investigation of MClPz and MClTTDPz included geometry optimizations and computations of the harmonic vibrations followed by TDDFT calculations of the electronic absorption spectrum. The number of the calculated excited states was 30. The calculations were performed using the B3LYP functional with the use of pcseg-2 basis set [38] taken from the EMSL BSE library [39,40] for describing the electron shells of C, N, S, Cl and H atoms. For describing the 28 core electron shells and 11 valence electrons of the yttrium atom pseudopotential combined with corresponding basis set cc-pVTZ-PP [41] were used. In case of the lanthanum atom 46 core electron shells were described by pseudopotential [42] and 11 valence electrons were described by Def2-TZVPP [42] basis sets. For the description of 11 valence electrons of the Lu atom triple-ζ valence basis sets [43] were used and 60 core shells were described by an effective core pseudopotential [44] developed by the Stuttgart group. The Firefly QC [45] package, which is partially based on the GAMESS(US) [46] source code was used in all the calculations.

Optimized Cartesian coordinates of MPz and MTTDPz are available from Supplementary Materials.

Description of the vibrational modes is carried out based on the analysis of the distribution of the potential energy of normal vibrations by natural vibrational coordinates. This analysis was performed using the VibModule program [47].

The quantum theory of atoms in molecules (QTAIM) analysis [48] was performed using the AIMAll [49] software package. Topological parameters of ρ(r) in bond critical points and charges on atoms are collected in Supplementary Materials.

The molecular models and orbitals demonstrated in the paper were visualized by means of the Chemcraft program [50].

4. Conclusions The influence of the nature of the metal (Y, La, Lu) and the ligand (either porphyrazine or thiadiazole-annelated porphyrazine) on the geometry, electronic structure and electronic and vibrational absorption spectra of the macroheterocyclic complex was studied by DFT calculations at the B3LYP/pcseg-2 level. It has been shown that the ligand cavity can adjust itself to the size of the central atom in a way that results in doming distortion of the ligand without changing the structure of the periphery of the macrocycle. It has been found that the electron density distribution in the coordination cavity of metal complexes of porphyrazine and tetra(1,2,5-thiadiazole)porphyrazine slightly depends on the nature of the ligand and metal. The value of the lanthanide contraction Δr_Ln in the studied complexes has been determined, which turned out to be close to Δr_Ln in the series of lanthanide halogenides with three monovalent ligands and in the series of lanthanide compounds with three bidentate ligands. The incorporation of a substituent can enhance the effect of the metal on the spectral properties of the complexes. The electronic spectra of investigated complexes can be described by the model of Gouterman. It has been shown that in the calculated IR spectra the most intense bands cannot be assigned to any single vibrational mode and are the result of the superposition of vibrational modes.

Figure 1. Molecular models of porphyrazine (MClPz) (a) and tetra(1,2,5-thiadiazolo)porphyrazine (MClTTDPz) (b) complexes with atom labeling (M = Y, La, Lu).

Figure 2. Calculated time-dependent density functional theory (TDDFT) electronic absorption spectra for MClPz (a) and MClTTDPz (b) complexes.

Figure 3. Molecular orbital (MO) level diagram for MClPz and MClTTDPz complexes. The values of highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gaps are given in eV.

Figure 4. Influence of the metal (Y, La, Lu) on the molecular orbitals of MClPz complexes.

Figure 5. Influence of the metal (Y, La, Lu) on the molecular orbitals of MClTTDPz complexes.

Figure 6. Simulated vibrational spectra of the MClPz (a) and MClTTDPz (b) complexes.

	LaClPz	LaClTTDPz	LuClPz	LuClTTDPz	YClPz	YClTTDPz
M-Cl	2.711	2.680	2.493	2.468	2.523	2.497
M-Np	2.462	2.486	2.260	2.284	2.291	2.313
Np-Cα	1.369	1.377	1.370	1.379	1.370	1.379
Cα-Nm	1.330	1.321	1.329	1.319	1.329	1.319
Cα-Cβ	1.456	1.460	1.454	1.458	1.454	1.458
Cβ-Cβ	1.355	1.423	1.355	1.422	1.355	1.422
Cβ-H	1.077		1.077		1.077
Cβ-Nt		1.316		1.317		1.317
Nt-S		1.644		1.643		1.643
(Np…Np)opp	4.023	4.139	3.979	4.096	3.987	4.101
(Np…Np)adj	2.844	2.927	2.813	2.897	2.819	2.900
∠ (NpMCl)	125.2	123.6	118.3	116.3	119.5	117.6
∠ (MNpCα)	124.7	122.4	125.1	123.1	125.2	122.9
∠ (NpCαNm)	127.6	128.1	127.3	127.9	127.4	128.0
∠ (CαNmCα)	125.0	127.1	124.5	126.6	124.6	126.7
∠ (CαNpCα)	107.5	111.6	107.6	111.5	107.6	111.6
∠ (NtSNt)		100.3		100.4		100.4
M-X1 1	1.419	1.377	1.072	1.011	1.129	1.071
X1-X2 2		1.028		0.775		0.767

1 X1 is dummy atom located in center between Np atoms. 2 X2 is dummy atom located in center between S atoms.

	YClPz	LaClPz	LuClPz	YClTTDPz	LaClTTDPz	LuClTTDPz
∇2ρ (M-Np), a.u.	0.205	0.160	0.235	0.193	0.152	0.221
q(M|Np)	0.348	0.356	0.340	0.354	0.362	0.347
δ(M|Np)	0.349	0.366	0.357	0.337	0.353	0.346
q(M|Cl)	0.779	0.793	0.763	0.757	0.770	0.740
δ(M|Cl)	0.536	0.564	0.552	0.574	0.606	0.588
q(M)	+2.171	+2.21	+2.125	+2.173	+2.218	+2.127
q(Np)	−1.191	−1.179	−1.187	−1.164	−1.153	−1.16
q(Nm)	−1.121	−1.123	−1.121	−1.124	−1.125	−1.123
q(Cα)	0.928	0.920	0.929	0.963	0.958	0.963
q(Cβ)	0.017	0.014	0.017	0.526	0.525	0.526
q(ligand) 1	−1.391	−1.423	−1.362	−1.416	−1.447	−1.386

1 The sum of the net atomic charges on the ligand atoms.

State	Composition (%)	Λ (nm)	f
YClPz
11 E	2a1→e*(16)a2→e* (82)	510	0.15	λ1
41 E	1a1→e*(53)2a1→e*(38)a2→e* (5)	334	0.13	λ2
51 E	1a1→e*(41)2a1→e*(42)a2→e* (10)	312	0.60	λ3
LaClPz
11 E	3a1→e*(15)a2→e* (82)	510	0.15	λ1
41 E	2a1→e*(47)3a1→e*(44)a2→e* (4)	346	0.08	λ2
61 E	1a1→e* (12)2a1→e*(38)3a1→e*(36)a2→e* (11)	311	0.62	λ3
LuClPz
11E	2a1→e*(16)a2→e* (82)	510	0.16	λ1
41E	1a1→e*(56)2a1→e*(35)a2→e* (5)	332	0.12	λ2
51E	1a1→e*(39)2a1→e*(45)a2→e* (11)	311	0.62	λ3
YClTTDPz
11 E	5a1→1e*(5)3a2→1e* (91)	588	0.27	λ1
51 E	4a1→1e*(43)5a1→1e*(45)2b1→1e* (10)	342	0.23	λ2
61 E	3a1→1e*(5)4a1→1e*(41)5a1→1e*(36)3a2→1e*(6)3a2→2e* (7)	319	0.67	λ3
LaClTTDPz
11 E	5a1→1e*(6)3a2→1e* (91)	582	0.27	λ1
51 E	4a1→1e*(80)5a1→1e*(10)3b1→1e* (10)	344	0.08
61 E	3a1→1e*(21)4a1→1e*(7)5a1→1e*(59)3a2→1e*(5)3a2→2e* (6)	329	0.48	λ2
71 E	2a1→1e*(13) 3a1→1e*(61)5a1→1e* (14)	308	0.36	λ3
LuClTTDPz
11 E	5a1→1e*(5)3a2→1e* (91)	590	0.27	λ1
51 E	4a1→1e*(39)5a1→1e*(49)3b1→1e* (10)	342	0.26	λ2
61 E	3a1→1e*(5)4a1→1e*(45)5a1→1e*(32)3a2→1e*(6)3a2→2e* (7)	318	0.67	λ3

Frequency, cm−1	Irel, %	Symmetry	Assignment 1	Exp, cm−1
YClTTDPz
1101.6(ω80-ω81)	59	E	r(Np-Cα) (24), r(Nm-Cα) (28), r(Cα-Cβ) (15), r(Cβ-Nt) (9), φ(CβCβNt) (6), φ(CβNtS) (7)	1090 acac [20]
1280.0(ω87-ω88)	100	E	r(Np-Cα) (27), r(Cβ-Cβ) (10), r(Cβ-Nt) (10), φ(CαNpCα) (6), φ(NpCαNm) (10), φ(NmCαCβ) (9), φ(CαCβNt) (6)	1260 acac [20]
LaClTTDPz
1100.2(ω80-ω81)	61	E	r(Np-Cα) (28), r(Nm-Cα) (27), r(Cα-Cβ) (14), r(Cβ-Nt) (8), φ(CβCβNt) (5), φ(CβNtS) (6)
1274.8(ω87-ω88)	100	E	r(Np-Cα) (29), r(Cβ-Cβ) (8), r(Cβ-Nt) (9), φ(CαNpCα) (6), φ(NpCαNm) (10), φ(CαNmCα) (5), φ(NmCαCβ) (9), φ(CαCβNt) (6)
LuClTTDPz
1102.1(ω80-ω81)	59	E	r(Np-Cα) (24), r(Nm-Cα) (28), r(Cα-Cβ) (15), r(Cβ-Nt) (9), φ(CβCβNt) (6), φ(CβNtS) (7)	1090 acac [20]
1281.5(ω87-ω88)	100	E	r(Np-Cα) (26), r(Cβ-Cβ) (10), r(Cβ-Nt) (10), φ(CαNpCα) (6), φ(NpCαNm) (10), φ(NmCαCβ) (9), φ(CαCβNt) (8)	1262 acac [20]
YClPz
839.81 (ω48)	46	A1	OPB(Cβ-Np-Nm-Cα) (35), OPB(H-Cα-Cβ-Cβ) (44), θ(Cα-Nm-Cα-Np) (5), θ(Nm-Cα-Cβ-Cβ) (11)
982.30 (ω53-ω54)	100	E	r(Np-Y) (5), r(Np-Cα) (14), r(Cα-Cβ) (53), r(Cα-Cβ) (5)
LaClPz
838.13 (ω48)	44	A1	OPB(Cβ-Np-Nm-Cα) (35), OPB(H-Cα-Cβ-Cβ) (44), θ(Cα-Nm-Cα-Np) (5), θ(Nm-Cα-Cβ-Cβ) (11)
975.65(ω53-ω54)	100	E	r(Np-La) (5), r(Np-Cα) (15), r(Nm-Cα) (5), r(Cα-Cβ) (56)
LuClPz
839.62 (ω48)	46	A1	OPB(Cβ-Np-Nm-Cα) (35), OPB(H-Cα-Cβ-Cβ) (45), θ(Cα-Nm-Cα-Np) (5), θ(Nm-Cα-Cβ-Cβ) (11)
983.45 (ω53-ω54)	100	E	r(Np-Lu) (5), r(Np-Cα) (7), r(Np-Cα) (7), r(Cα-Cβ) (26), r(Cα-Cβ) (27)

1 Coordinates are listed provided that their contributions (shown in parentheses) are greater than ~5%. Assignment of vibrational modes based on potential energy distribution. The following designations of the coordinates are used: r—stretching of the bond; φ—bending, a change in the angle; OPB—out-of-plane bending; θ—a change in the dihedral angle.

Supplementary Materials

Cartesian coordinates of YClPz optimized B3LYP/pcseg-2 level of theory, Cartesian coordinates of LaClPz optimized B3LYP/pcseg-2 level of theory, Cartesian coordinates of LuClPz optimized B3LYP/pcseg-2 level of theory, Cartesian coordinates of YClTTDPz optimized B3LYP/pcseg-2 level of theory, Cartesian coordinates of LaClTTDPz optimized B3LYP/pcseg-2 level of theory, Cartesian coordinates of LuClTTDPz optimized B3LYP/pcseg-2 level of theory, Table S1: Calculated composition of the lowest excited states and corresponding oscillator strengths for MClPz and MClTTDPz complexes, Table S2: Assignment of the IR vibrations of the MClPz and MClTTDPz complexes, Table S3: Bond lengths and topological parameters of ρ(r) in bond critical points of the MClPz, Table S4: Bond lengths and topological parameters of ρ(r) in bond critical points of the MClTTDPz, Table S5: Charge on atoms in MClPz and MClTTDPz.

Author Contributions

Conceptualization, P.A.S.; Methodology, Y.A.Z.; Investigation, I.A.K. and A.V.E.; Resources, Y.A.Z.; Data Curation, I.V.R.; Writing-original draft preparation, I.V.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Russian Science Foundation (grant № 19-73-00256).

Data Availability Statement

The data presented in this study are available in supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Sorokin, A.B. Recent progress on exploring µ-oxo bridged binuclear porphyrinoid complexes in catalysis and material science. Coord. Chem. Rev. 2019, 389, 141-160.

2. Yu, Z.; Hagfeldt, A.; Sun, L. The application of transition metal complexes in hole-transporting layers for perovskite solar cells: Recent progress and future perspectives. Coord. Chem. Rev. 2020, 406, 213143.

3. Almeida-Marrero, V.; Van De Winckel, E.; Anaya-Plaza, E.; Torres, T.; De La Escosura, A. Porphyrinoid biohybrid materials as an emerging toolbox for biomedical light management. Chem. Soc. Rev. 2018, 47, 7369-7400.

4. Wöhrle, D.; Schnurpfeil, G.; Makarov, S.G.; Kazarin, A.; Suvorova, O.N. Practical Applications of Phthalocyanines-From Dyes and Pigments to Materials for Optical, Electronic and Photo-electronic Devices. Macroheterocycles 2012, 5, 191-202.

5. Lukyanets, E.A.; Nemykin, V.N. The key role of peripheral substituents in the chemistry of phthalocyanines and their analogs. J. Porphyr. Phthalocyanines 2010, 14, 1-40.

6. Longevial, J.; Clément, S.; Wytko, J.A.; Ruppert, R.; Weiss, J.; Richeter, S. Peripherally Metalated Porphyrins with Applications in Catalysis, Molecular Electronics and Biomedicine. Chem. Eur. J. 2018, 24, 15442-15460.

7. Bettini, S.; Valli, L.; Giancane, G. Applications of Photoinduced Phenomena in Supramolecularly Arranged Phthalocyanine Derivatives: A Perspective. Molecules 2020, 25, 3742.

8. Donzello, M.P.; Ercolani, C.; Novakova, V.; Zimcik, P.; Stuzhin, P.A. Tetrapyrazinoporphyrazines and their metal derivatives. Part I: Synthesis and basic structural information. Coord. Chem. Rev. 2016, 309, 107-179.

9. Novakova, V.; Donzello, M.P.; Ercolani, C.; Zimcik, P.; Stuzhin, P.A. Tetrapyrazinoporphyrazines and their metal derivatives. Part II: Electronic structure, electrochemical, spectral, photophysical and other application related properties. Coord. Chem. Rev. 2018, 361, 1-73.

10. Donzello, M.P.; Ercolani, C.; Stuzhin, P.A. Novel families of phthalocyanine-like macrocycles-Porphyrazines with annulated strongly electron-withdrawing 1,2,5-thia/selenodiazole rings. Coord. Chem. Rev. 2006, 250, 1530-1561.

11. Suzuki, Y.; Fujimori, M.; Yoshikawa, H.; Awaga, K. Packing motifs and magneto-structural correlations in crystal structures of metallo-tetrakis(1,2,5-thiadiazole)porphyrazine series, MTTDPz (M=H 2, Fe, Co, Ni, Cu, Zn). Chem. Eur. J. 2004, 10, 5158-5164.

12. Koifman, O.I.; Ageeva, T.A.; Beletskaya, I.P.; Averin, A.D.; Yakushev, A.A.; Tomilova, L.G.; Dubinina, T.V.; Tsivadze, A.Y.; Gorbunova, Y.G.; Stuzhin, P.A.; et al. Macroheterocyclic Compounds-A Key Building Block in New Functional Materials and Molecular Devices. Macroheterocycles 2020.

13. Weiss, R.; Fischer, J. Lanthanide Phthalocyanine Complexes. In The Porphyrin Handbook: Phthalocyanines: Spectroscopic and Electrochemical Characterization; Elsevier Inc.: Cambridge, MA, USA, 2003; Volume 16, pp. 171-246. ISBN 9780123932266.

14. Bouvet, M.; Gaudillat, P.; Suisse, J.-M. Lanthanide macrocyclic complexes: From molecules to materials and from materials to devices. J. Porphyr. Phthalocyanines 2013, 17, 628-635.

15. Stuzhin, P.A.; Ivanova, S.S.; Hamdoush, M.; Kirakosyan, G.A.; Kiselev, A.; Popov, A.; Sliznev, V.; Ercolani, C. Tetrakis(1,2,5-thiadiazolo)porphyrazines. 9. Synthesis and spectral and theoretical studies of the lithium(i) complex and its unusual behaviour in aprotic solvents in the presence of acids. Dalton Trans. 2019, 48, 14049-14061.

16. Stuzhin, P.A.; Bauer, E.M.; Ercolani, C. Tetrakis(thiadiazole)porphyrazines. 1. Syntheses and Properties of Tetrakis(thiadiazole)porphyrazine and Its Magnesium and Copper Derivatives. Inorg. Chem. 1998, 37, 1533-1539.

17. Otlyotov, A.A.; Ryzhov, I.V.; Kuzmin, I.A.; Zhabanov, Y.A.; Mikhailov, M.S.; Stuzhin, P.A. DFT Study of Molecular and Electronic Structure of Ca(II) and Zn(II) Complexes with Porphyrazine and tetrakis(1,2,5-thiadiazole)porphyrazine. Int. J. Mol. Sci. 2020, 21, 2923.

18. Donzello, M.P.; Agostinetto, R.; Ivanova, S.S.; Fujimori, M.; Suzuki, Y.; Yoshikawa, H.; Shen, J.; Awaga, K.; Ercolani, C.; Kadish, K.M.; et al. Tetrakis(thiadiazole)porphyrazines. 4. Direct template synthesis, structure, general physicochemical behavior, and redox properties of Al III, GaIII, and InIII complexes. Inorg. Chem. 2005, 44, 8539-8551.

19. Bauer, E.M.; Cardarilli, D.; Ercolani, C.; Stuzhin, P.A.; Russo, U. Tetrakis(thiadiazole)porphyrazines. 2. Metal complexes with Mn(II), Fe(II), Co(II), Ni(II), and Zn(II). Inorg. Chem. 1999, 38, 6114-6120.

20. Tarakanova, E.N.; Hamdoush, M.; Eroshin, A.V.; Ryzhov, I.V.; Zhabanov, Y.A.; Stuzhin, P.A. Tetra(1,2,5-thiadiazolo)porphyrazines. 10. Synthesis, spectral characterization and DFT study of complexes with yttrium(III) and lutetium(III). Polyhedron 2021, 193, 114877.

21. Khelevina, O.G.; Voinov, A.A. Synthesis and kinetic stability of dysprosium(III), europium(III), and neodymium(III) complexes with tetraazaporphine in acetic acid solutions. Russ. J. Gen. Chem. 2000, 70, 778-783.

22. Khelevina, O.G.; Voinov, A.A. Kinetic stability of the dysprosium(III) complex with tetraazaporphine in water-acetic acid and methanol-acetic acid mixtures. Russ. J. Coord. Chem. Khimiya 1999, 25, 445-449.

23. Zhabanov, Y.A.; Sliznev, V.V.; Ryzhov, I.V.; Stuzhin, P.A. Peculiarities of electronic structure and chemical bonding in iron and cobalt metal complexes of porphyrazine and tetra(1,2,5-thiadiazole)porphyrazine. J. Porphyr. Phthalocyanines 2020, 24, 1146-1154.

24. Bader, R.F.W.; Slee, T.S.; Cremer, D.; Kraka, E. Description of conjugation and hyperconjugation in terms of electron distributions. J. Am. Chem. Soc. 1983, 105, 5061-5068.

25. Shannon, R.D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Crystallogr. Sect. A 1976, 32, 751-767.

26. Zhabanov, Y.A.; Giricheva, N.I.; Zakharov, A.V.; Islyaikin, M.K. Distinctive features of the structure of hemihexaphyrazine complexes with Y, La, and Lu according to quantum chemical data. J. Mol. Struct. 2017, 1132.

27. Giricheva, N.I.; Zakharov, A.V.; Shlykov, S.A.; Girichev, G.V. Molecular structure of GdCl3. Nuclear dynamics of the trichlorides of Gd, Tm, and Lu. J. Chem. Soc. Dalton Trans. 2000, 19, 3401-3403.

28. Zakharov, A.V.; Vogt, N.; Shlykov, S.A.; Giricheva, N.I.; Galanin, I.E.; Girichev, G.V.; Vogt, J. The molecular structures of LaCl3 and LaBr3 reinvestigated by gas-phase electron diffraction. J. Mol. Struct. 2004, 707, 147-152.

29. Shlykov, S.A.; Giricheva, N.I.; Lapykina, E.A.; Girichev, G.V.; Oberhammer, H. The molecular structure of TbI3, DyI3, HoI3 and ErI3 as determined by synchronous gas-phase electron diffraction and mass spectrometric experiment assisted by quantum chemical calculations. J. Mol. Struct. 2010, 978, 170-177.

30. Giricheva, N.I.; Shlykov, S.A.; Lapykina, E.A.; Oberhammer, H.; Girichev, G.V. The molecular structure of PrI3 and GdI3 as determined by synchronous gas-phase electron diffraction and mass spectrometric experiment assisted by quantum chemical calculations. Struct. Chem. 2011, 22, 385-392.

31. Ischenko, A.A.; Girichev, G.V.; Tarasov, Y.I. Difraktciya Elektronov: Struktura i Dinamika Svobodnykh Molekul i Kondensirovannogo Sostoyaniya Veshestva; Fizmatlit: Moscow, Russia, 2012.

32. Giricheva, N.I.; Belova, N.V.; Shlykov, S.A.; Girichev, G.V.; Vogt, N.; Tverdova, N.V.; Vogt, J. Molecular structure of tris(dipivaloylmethanato)lanthanum(III) studied by gas electron diffraction. J. Mol. Struct. 2002, 605, 171-176.

33. Belova, N.V.; Girichev, G.V.; Hinchley, S.L.; Kuzmina, N.P.; Rankin, D.W.H.; Zaitzeva, I.G. Molecular structure of tris(dipivaloylmethanato)lutetium(iii) studied by gas electron diffraction and ab initio and DFT calculations. Dalton Trans. 2004, 1715-1718.

34. Gouterman, M. Spectra of porphyrins. J. Mol. Spectrosc. 1961, 6, 138-163.

35. Gouterman, M.; Wagnière, G.H.; Snyder, L.C. Spectra of porphyrins. Part II. Four orbital model. J. Mol. Spectrosc. 1963, 11, 108-127.

36. Weiss, C.; Kobayashi, H.; Gouterman, M. Spectra of porphyrins. Part III. Self-consistent molecular orbital calculations of porphyrin and related ring systems. J. Mol. Spectrosc. 1965, 16, 415-450.

37. Zhabanov, Y.A.; Tverdova, N.V.; Giricheva, N.I.; Girichev, G.V.; Stuzhin, P.A. DFT Study of molecular and electronic structure of magnesium (II) tetra(1,2,5-chalcogenadiazolo) porphyrazines, [TXDPzMg] (X = O, S, Se, Te). J. Porphyr. Phthalocyanines 2017, 21, 439-452.

38. Jensen, F. Unifying general and segmented contracted basis sets. Segmented polarization consistent basis sets. J. Chem. Theory Comput. 2014, 10, 1074-1085.

39. Pritchard, B.P.; Altarawy, D.; Didier, B.; Gibson, T.D.; Windus, T.L. New Basis Set Exchange: An Open, Up-to-Date Resource for the Molecular Sciences Community. J. Chem. Inf. Model. 2019.

40. Schuchardt, K.L.; Didier, B.T.; Elsethagen, T.; Sun, L.; Gurumoorthi, V.; Chase, J.; Li, J.; Windus, T.L. Basis set exchange: A community database for computational sciences. J. Chem. Inf. Model. 2007, 47, 1045-1052.

41. Peterson, K.A.; Figgen, D.; Dolg, M.; Stoll, H. Energy-consistent relativistic pseudopotentials and correlation consistent basis sets for the 4d elements Y-Pd. J. Chem. Phys. 2007, 126, 124101.

42. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. Des. Assess. Accuracy 2005, 7, 3297-3305.

43. Weigand, A.; Cao, X.; Yang, J.; Dolg, M. Quasirelativistic f-in-core pseudopotentials and core-polarization potentials for trivalent actinides and lanthanides: Molecular test for trifluorides. Theor. Chem. Acc. 2009, 126, 117-127.

44. Dolg, M.; Stoll, H.; Savin, A.; Preuss, H. Energy-adjusted pseudopotentials for the rare earth elements. Theor. Chim. Acta 1989, 75, 173-194.

45. Granovsky, A.A. Firefly Version 8. Available online: http://classic.chem.msu.su/gran/firefly/index.html (accessed on 28 December 2020).

46. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; et al. General atomic and molecular electronic structure system. J. Comput. Chem. 1993, 14, 1347-1363.

47. Vishnevskiy, Y.V.; Zhabanov, Y.A. New implementation of the first-order perturbation theory for calculation of interatomic vibrational amplitudes and corrections in gas electron diffraction. J. Phys. Conf. Ser. 2015, 633, 12076.

48. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, UK; Clarendon Press: Oxford, UK, 1990; ISBN 9780198558651.

49. Some References Related to AIMAll. Available online: http://aim.tkgristmill.com/references.html (accessed on 13 March 2020).

50. Zhurko, G.A.; Zhurko, D.A. ChemCraft Version 1.6 (Build 312); Version 1.6 (Build 312) Ed. Available online: http://www.chemcraftprog.com/index.html (accessed on 28 December 2020).