Molecular Structure, Thermodynamic and Spectral Characteristics of Metal-Free and Nickel Complex of Tetrakis(1,2,5-thiadiazolo)porphyrazine

The Knudsen effusion method with mass spectrometric control of the vapor composition was used to study the possibility of a congruent transition to the gas phase and to estimate the enthalpy of sublimation of metal-free tetrakis(1,2,5-thiadiazolo)porphyrazine and its nickel complex (H2TTDPz and NiTTDPz, respectively). The geometrical and electronic structure of H2TTDPz and NiTTDPz in ground and low-lying excited electronic states were determined by DFT calculations. The electronic structure of NiTTDPz was studied by the complete active space (CASSCF) method, following accounting dynamic correlation by multiconfigurational quasi-degenerate second-order perturbation theory (MCQDPT2). A geometrical structure of D2h and D4h symmetry was obtained for H2TTDPz and NiTTDPz, respectively. According to data obtained by the MCQDPT2 method, the nickel complex possesses the ground state 1A1g, and the wave function of the ground state has the form of a single determinant. Electronic absorption and vibrational (IR and resonance Raman) spectra of H2TTDPz and NiTTDPz were studied experimentally and simulated theoretically.


Introduction
There is a growing interest in organic materials for application in optoelectronics due to their low-cost, high-throughput film manufacturing by solution-processing techniques, high-flex-stability, easy scaling up and integration in devices [1][2][3]. Compounds based on tetrapyrrole macrocycles have interesting spectral and non-linear optical properties [4] which are practically useful for optical communication, information storage, optical switching and processing of electro-optical signals [5][6][7][8][9][10].
The structural formulas of the molecules discussed here are presented in the Supplementary Materials. The unsubstituted porphyrazine (or tetra-azaporphyrin) and its metallic derivatives, including the Mg(II), Ni(II) and Cu(II) complexes (MPz), were synthesized in 1952 [11]. Their tetrabenzo fused derivatives, known as phthalocyanines are so far most widely studied. Substitution of benzene rings in phthalocyanines by aromatic heterocycles, e.g., by pyrazine [12,13] or 1,2,5-chalcogenadiazole [14], have a strong impact on the electronic properties of the central porphyrazine (Pz) core which is common to these systems.
The presence of a five-membered heterocycle containing nitrogen and sulfur atoms on the periphery of tetrakis(1,2,5-thiadiazolo)porphyrazine (TTDPz) considerably modulates the physico-chemical properties of the macrocycle and its intermolecular interaction as compared to phthalocyanines [14]. Unlike phthalocyanines, or their pyrazine-fused analogues, TTDPz have no H atoms on the periphery and their molecular packing during of H 2 TTDPz, we came to the conclusion that its sublimation is accompanied by partial decomposition. The ion H 2 TTDPz + has lower intensity than another heavy ion C 4 N 4 S + at m/z = 136 a.m.u. The C 4 N 4 S + ion corresponds to 1⁄4(TTDPz) and can be assigned to the C 4 N 4 S molecular species-1,2,5-thiadiazole-3,4-dicarbonitrile. The relative abundance of ions in mass spectrum of H 2 TTDPz recorded in the present work is similar to the relative abundance of ions in mass spectrum of ZnTTDPz recorded simultaneously with gas electron diffraction experiment [29]. However, the recording temperature of the ZnTTDPz mass spectrum was 861 K [29], which is 171 K higher than the recording temperature of the H 2 TTDPz mass spectrum. m/z = 136 a.m.u. The C4N4S + ion corresponds to 1⁄4(TTDPz) and can be assigned to the C4N4S molecular species-1,2,5-thiadiazole-3,4-dicarbonitrile. The relative abundance of ions in mass spectrum of H2TTDPz recorded in the present work is similar to the relative abundance of ions in mass spectrum of ZnTTDPz recorded simultaneously with gas electron diffraction experiment [29]. However, the recording temperature of the ZnTTDPz mass spectrum was 861 K [29], which is 171 K higher than the recording temperature of the H2TTDPz mass spectrum.  Mass-spectrometric studies performed in connection with the sublimation process show that NiTTDPz gives a stable stream of particles at temperatures T = 632-717 K, where a molecular ion (m/z = 602) dominates, followed by several ions of weaker (3-4%) intensity. No ions corresponding to oligomeric species were detected. According to the mass spectra analysis, we are inclined to conclude that saturated vapor consists only of the parent NiTTDPz.
The plot of ln(IT) = f(1000/T) for NiTTDPz + ion is shown in Figure 2. The set of the points is the result of measurements with a step-by-step increase in temperature, and then with a step-by-step decrease. Each point of the graph corresponds to the ion current measured after its stabilization at a given temperature. One can see that the hysteresis as the   52 25 Mass-spectrometric studies performed in connection with the sublimation process show that NiTTDPz gives a stable stream of particles at temperatures T = 632-717 K, where a molecular ion (m/z = 602) dominates, followed by several ions of weaker (3-4%) intensity. No ions corresponding to oligomeric species were detected. According to the mass spectra analysis, we are inclined to conclude that saturated vapor consists only of the parent NiTTDPz.
The plot of ln(IT) = f(1000/T) for NiTTDPz + ion is shown in Figure 2. The set of the points is the result of measurements with a step-by-step increase in temperature, and then with a step-by-step decrease. Each point of the graph corresponds to the ion current measured after its stabilization at a given temperature. One can see that the hysteresis as the temperature increases and decreases is practically absent. This allows us to conclude that points in plots are corresponding to the equilibrium states inside the effusion cell. The dependence ln(IT) = f(1000/T) could be closely approximated by a straight line, usually observed for vaporization in the considered temperature range without change of crystallographic modification of the solid phase and without significant change of the enthalpy of vaporization. The enthalpy of sublimation value ∆H s calculated by linear regression using the Clausius-Clayperon equation Ln(IT) = − ∆H RT + C was found to be 246(2) kJ·mol −1 .
temperature increases and decreases is practically absent. This allows us to conclude that points in plots are corresponding to the equilibrium states inside the effusion cell. The dependence ln(IT) = f(1000/T) could be closely approximated by a straight line, usually observed for vaporization in the considered temperature range without change of crystallographic modification of the solid phase and without significant change of the enthalpy of vaporization. The enthalpy of sublimation value ΔHs calculated by linear regression using the Clausius-Clayperon equation ln( ) IT H RT C = − + Δ was found to be 246 (2) kJ·mol −1 .

Molecular Structure
The electronic configuration of Ni(II) is [Ar]3d 8 , and therefore it can form in ground state as either singlet or triplet complexes. Furthermore, the computational investigations in the case of Ni(II) complexes are often non-trivial due to the necessity to account for the multireference character of the wavefunction.
The electronic structure of NiTTDPz has been studied by the CASSCF method followed by an account of the dynamic electron correlation by multiconfigurational quasidegenerate second-order perturbation theory (MCQDPT2). The compositions of the wave functions are presented in Table 2 for the low-lying electronic states. According to data obtained by the MCQDPT2 method, the NiTTDPz complex possesses the ground states 1 A1g. The low-lying triplet state is higher by 88.4 kJ mol -1 in energy than the corresponding ground state (Table 1). It should be noted that, according to CASSCF calculations, NiTT-DPz possesses a triplet ground state. Such contradictory conclusions about the multiplicity of the ground state obtained using the CASSCF and MCQDPT2 methods are apparently due to the fact that the CASSCF calculations with a small active space do not practically take into account the dynamic correlation of electrons. An analysis of the data in Table 2 shows that the wave functions of the ground states and the most low-lying triplet states have the form of a single determinant. Therefore, for the D4h configuration in the electronic state 1 A1g, the geometry optimization, calculations of the force field, and vibrational and electronic spectra have been performed using the PBE0/pcseg-2 approach.

Molecular Structure
The electronic configuration of Ni(II) is [Ar]3d 8 , and therefore it can form in ground state as either singlet or triplet complexes. Furthermore, the computational investigations in the case of Ni(II) complexes are often non-trivial due to the necessity to account for the multireference character of the wavefunction.
The electronic structure of NiTTDPz has been studied by the CASSCF method followed by an account of the dynamic electron correlation by multiconfigurational quasi-degenerate second-order perturbation theory (MCQDPT2). The compositions of the wave functions are presented in Table 2 for the low-lying electronic states. According to data obtained by the MCQDPT2 method, the NiTTDPz complex possesses the ground states 1 A 1g . The lowlying triplet state is higher by 88.4 kJ mol −1 in energy than the corresponding ground state (Table 1). It should be noted that, according to CASSCF calculations, NiTTDPz possesses a triplet ground state. Such contradictory conclusions about the multiplicity of the ground state obtained using the CASSCF and MCQDPT2 methods are apparently due to the fact that the CASSCF calculations with a small active space do not practically take into account the dynamic correlation of electrons. An analysis of the data in Table 2 shows that the wave functions of the ground states and the most low-lying triplet states have the form of a single determinant. Therefore, for the D 4h configuration in the electronic state 1 A 1g, the geometry optimization, calculations of the force field, and vibrational and electronic spectra have been performed using the PBE0/pcseg-2 approach.
In ref [23] it was shown that crystal field theory (CFT) can be used to describe the sequence of electronic states of MPz and MTTDPz (M=Fe, Co) complexes. However, in the case of the singlet state of nickel hemi-porphyrazine [34], it is impossible to describe the sequence of electronic states using the crystal field theory. In the framework of this theory, the most energetically favorable states are those with the least repulsion between the electrons occupying the d-shell of the metal and orbitals of the macrocycle. From this point of view, the occupation of the b 2g , e g , and a 1g MOs are the most favorable, but not b 1g .

State
Contributions ∆E, kJ/mol 1 A 1g 96/(a 2u ) 2 (a 1u ) 2 (e g ) 4 (b 1g ) 0 / 0.0 1 E g 99/(a 2u ) 2 (a 1u ) 2 (e g ) 3 (b 1g ) 1 / 263.0 1 E u 97/(a 2u ) 2 (a 1u ) 2 (e g ) 2 (b 1g ) 2 / 578.1 1 A 1g 97/(a 2u ) 2 (a 1u ) 2 (e g ) 2 (b 1g ) 2 / 581.5 1 E u 96/(a 2u ) 2 (a 1u ) 2 (e g ) 2 (b 1g ) 2 / 737.2 3 B 1g 100/(b 2g ) 2 (e g ) 4 (a 1g ) 1 (b 1g ) 1 / 82.3 3 E g 93/(b 2g ) 2 (e g ) 3 (a 1g ) 2 (b 1g ) 1 / 92.6 3 A 2g 30/(b 2g ) 2 (e g ) 2 (a 1g ) 2 (b 1g ) 2 /+70/(b 2g ) 1 (e g ) 4 (a 1g ) 2 (b 1g ) 1 / 174.8 3 E g 32/(b 2g ) 1 (e g ) 3 (a 1g ) 2 (b 1g ) 2 /+47/(b 2g ) 1 (e g ) 3 (a 1g ) 2 (b 1g ) 2 /+16/(b 2g ) 2 (e g ) 3 (a 1g ) 1 (b 1g ) 2 / 331.0 3 B 2g 100/(b 2g ) 1 (e g ) 4 (a 1g ) 1 (b 1g ) 2 / 338.5 The crystal field theory (CFT, [35,36]) can be used to describe the sequence of NiTTDPz electronic states (Table 1) despite the fact that the shapes of two CASSCF active molecular orbitals (a 2u and a 1u , Figure 3) in the singlet NiTTDPz state comprise atoms of the macrocycle rather than the metal atom. This conclusion is confirmed by the fact that the b 1g orbital is unoccupied in the ground state and this orbital has significantly greater energy than the other three orbitals ( Figure 4).   Shapes of CASSCF active molecular orbitals ( Figure 3) of triplet NiTTDPz state and their composition analysis show that the corresponding components of the d-orbitals of the metal atom make a dominant contribution. The orbitals of the macrocycle atoms are almost not involved in the formation of these molecular orbitals. The orbital of b 1g symmetry is an exception, since according to Figure 3 the contribution of the macrocycle orbitals can be visually observed. It should be noted that no noticeable interaction of metal d-orbitals and macrocycle orbitals was found. Thus, the crystal field theory (CFT) can be used to describe the sequence of electronic states. A diagram of the energies of active in the CASSCF calculations molecular orbitals ( Figure 4) confirms this conclusion.  According to DFT calculations the molecules under consideration have a planar structure with symmetry D2h and D4h for H2TTDPz and NiTTDPz, respectively ( Figure 5). Note that the singlet and triplet states of the NiTTDPz molecule have significant differences in geometric structure (Table 3). In the case of triplet NiTTDPz, a significant (about 0.1 Å) increase in the (Np…Np)opp and, accordingly, in Np-M and (Np…Np)adj bond length (about 0.05 Å) distance is observed. When analyzing the data in Table 3, it can be noted that the smallest size of coordination cavity is observed for the NiTTDPz molecule in the singlet state, and the largest is for the H2TTDPz molecule. This confirms the conclusion of the influence of the metal on the size of the coordination cavity [23,27,32]. According to DFT calculations the molecules under consideration have a planar structure with symmetry D 2h and D 4h for H 2 TTDPz and NiTTDPz, respectively ( Figure 5). Note that the singlet and triplet states of the NiTTDPz molecule have significant differences in geometric structure (Table 3). In the case of triplet NiTTDPz, a significant (about 0.1 Å) increase in the (N p . . . N p ) opp and, accordingly, in N p -M and (N p . . . N p ) adj bond length (about 0.05 Å) distance is observed. When analyzing the data in Table 3, it can be noted that the smallest size of coordination cavity is observed for the NiTTDPz molecule in the singlet state, and the largest is for the H 2 TTDPz molecule. This confirms the conclusion of the influence of the metal on the size of the coordination cavity [23,27,32].

Electronic Absorption Spectra
Analyzing the electronic absorption spectra of molecules simulated by the TDDFT method, one can notice significant differences ( Figure 6). The Q-band is located at about 564 nm in the spectrum of NiTTDPz, while in the spectrum of the H 2 TTDPz, splitting into two bands, Q y (584 nm) and Q x (554 nm), is observed. In addition to significant changes in the region of the Q-band, one can also observe a change in relative intensities in the region of the B-band (≈300-330 nm). It is clearly seen that the B 2 -band in the spectrum of NiTTDPz has a lower intensity in comparison with the B 1 -band. In the case of H 2 TTDPz, these bands are close in intensity. The obtained spectra of H 2 TTDPz and NiTTDPz can be described using the four-orbital model of Gouterman [37][38][39]; both spectra in general are quite typical for this class of compounds. 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 4 300 400 500 600 700 800 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 Tables 4  and 5 for H 2 TTDPz and NiTTDPz, respectively. Analysis of the data in Tables 4 and 5 demonstrates that, in the case of the NiTTDPz complex, it can be seen that the Q-band corresponds to the transition between the frontier orbitals and the formation of the 1 1 E u state. The Q y and Q x bands in the H 2 TTDPz spectrum correspond to transitions from the ground state to the excited electronic states 1 1 B 2u and 1 1 B 3u , respectively. The 1 1 B 3u state is formed due to the HOMO-1→LUMO and HOMO→LUMO+1 transitions and the 1 1 B 2u state is formed due to the transition between the frontier orbitals. The excited states with a strong contribution to the electronic transition from the Gouterman type a 2u (NiTTDPz) and b 1u (H 2 TTDPz) orbitals, localized mainly on the inner nitrogen atoms of porphyrazine, have the highest intensity in the calculated spectra and are denoted as B band. Table 4. Calculated composition of the lowest excited states and corresponding oscillator strengths for H 2 TTDPz molecule.

State
Composition 333 0.02 311 0.55 362 294 0.19 (89) 257 0.26 (88) 244 0.07 (7) 240 0.34 231 0.05 The shapes of molecular orbitals (MOs) that participate in electronic transitions with large oscillator strengths are shown in Figure 7 and the energy diagram of molecular orbitals for H 2 TTDPz and NiTTDPz is shown in the Figure 8. The symmetry of the frontier MOs in the NiTTDPz complex is also typical for metal complexes of porphyrazines and tetrakis(1,2,5-thiadiazolo)porphyzarines-a 1u for HOMO and a pair of degenerate orbitals e g * for LUMO [26,27,32,40,41]. The shapes of HOMO orbitals for H 2 TTDPz and NiTTDPz are similar, despite different types of symmetry. The b 2g * orbital (LUMO) of the H 2 TTDPz molecule is a linear combination of bonding orbitals along the C α -C β and N m -C α bonds, as well as antibonding orbitals predominantly belonging to N t atoms and pairs of N p and S atoms. In the case of NiTTDPz, the shape of the 1e g * orbital can be characterized in the same way, but the contribution of the d-orbital of the nickel atom is noticeable. In addition to the frontier orbitals, HOMO-3 (H 2 TTDPz), HOMO-4 (NiTTDPz), and HOMO-7 of both compounds are involved in a large number of electronic transitions. Analyzing the composition of the HOMO-3 and HOMO-4 orbitals of the corresponding molecules, one can note the fact that, despite different symmetry (b 3u and a 2u , respectively), these MOs are practically identical in composition and only slightly differ in the regions of the N p atoms.
atoms. In the case of NiTTDPz, the shape of the 1eg * orbital can be characterized in the same way, but the contribution of the d-orbital of the nickel atom is noticeable. In addition to the frontier orbitals, HOMO-3 (H2TTDPz), HOMO-4 (NiTTDPz), and HOMO-7 of both compounds are involved in a large number of electronic transitions. Analyzing the composition of the HOMO-3 and HOMO-4 orbitals of the corresponding molecules, one can note the fact that, despite different symmetry (b3u and a2u, respectively), these MOs are practically identical in composition and only slightly differ in the regions of the Np atoms.

NiTTDPz H2TTDPz
LUMO, 1 * LUMO, 1 * The HOMO-LUMO gap is 2.40 eV for H2TTDPz and 2.56 eV for NiTTDPz, respectively. It should be noted that this value is typical for unsubstituted porphyrazines, but falls outside the range of previously studied thiadiazol annulated porphyrazines (2.24-

Vibrational Spectra
Theoretical results were also used for interpretation of the experimental vibrational spectra of H2TTDPz and NiTTDPz. The IR and Raman spectra were simulated on the basis of the normal mode frequencies and band intensities, which have been calculated by the DFT (PBE0/pcseg-2) method in a harmonic approximation. Description of the main IR active vibrations is presented in Table 6. As can be seen from Figure 9, the simulated IR spectra have fairly good correspondence with the experimental spectra and theoretical data can be used for assignment of the most intense bands.

Vibrational Spectra
Theoretical results were also used for interpretation of the experimental vibrational spectra of H 2 TTDPz and NiTTDPz. The IR and Raman spectra were simulated on the basis of the normal mode frequencies and band intensities, which have been calculated by the DFT (PBE0/pcseg-2) method in a harmonic approximation. Description of the main IR active vibrations is presented in Table 6. As can be seen from Figure 9, the simulated IR spectra have fairly good correspondence with the experimental spectra and theoretical data can be used for assignment of the most intense bands.
The two most intensive bands, located for NiTTDPz in the 1100-1350 cm −1 region, correspond to skeletal vibrations of the macrocycle with a predominant contribution of the N p -C α stretching vibrations. For H 2 TTDPz four intense bands are present in this region and have comparable contribution from the stretching vibrations of the N p -C α and C α -C β bonds in the pyrrole and pyrrolenine rings and in-plane deformation modes. With increasing frequency, a decrease of the out-of-plane vibrations contribution is observed. The ratio of stretching vibrations at the periphery also increases. One strong band is present for NiTTDPz and three less intense bands for H 2 TTDPz in the 1500-1650 cm −1 region. These bands make a considerable contribution to the stretching vibrations of the bridging N m -C α bonds, C α -C β bonds in pyrrole rings and C β -N t bonds in the fused heterocycle. The out-of-plane deformations of 1,2,5-thiadiazole rings appear as medium-strong bands at 500-550 cm −1 . Calculations predict a medium peak of N p -H stretching at 3554 cm −1 . In the experimental spectrum, this vibration is observed at lower frequencies 3291 cm −1 , evidencing existence of strong intramolecular hydrogen bonding the center of macrocycle.   The resonance Raman spectra of NiTTDPz and H 2 TTDPz obtained on rotating KBr disks using various excitation wavelengths at 80 K are displayed in Figures 10 and 11 along with the theoretical spectra representing Raman active vibrational modes. The description of the Raman vibrations is presented in Table 7. The calculated spectral patterns are in reasonable agreement with the experimental spectra, although the position of some bands varies by up to 100 cm −1 . Nevertheless, the theoretical data on the contribution of the vibrations of different fragments of the macrocycle to the normal modes are quite useful in the assignment of the experimental spectra.  (Table 5), which should be present near λex = 514.5 nm.
The metal free macrocycle H2TTDPz, due to the presence of two pyrrole and two pyrrolenine type fragments, has lower D2h symmetry and its resonance Raman spectra (see Figure 10) are richer than for NiTTDPz (D4h). It can be seen that intense vibrations of the macrocyclic skeleton appearing as single bands for NiTTDPz are split into two components in the experimental spectra of H2TTDPz (1572 → 1566 and 1555; 1534 → 1532 and 1517; 1182 → 1177 and 1158; 709 → 711 and 697 cm −1 ). The intensity ratio of the components is dependent on the excitation wavelength. When excitation is shifted from the Qto Soret band region, the bands containing a considerable contribution from the vibrational modes of the fused 1,2,5-thiadiazole rings are enhanced. This is not surprising since the Q-band transition is mainly localized on the atoms of the central macrocyclic core, while molecules with strong participation from the fused heterocycles participate in the electronic transitions in the Soret band region. Another remarkable feature in the experimental spectra of H2TTDPz is high intensity of the bands at 1430-1440, 1334 and 563 cm −1 , especially at the excitation at 454.5 nm.        In the Raman spectra of NiTTDPz recording using excitation in the virtual state (λ ex = 1064 nm), two intense bands are observed at 1572 and 1182 cm −1 . They are also dominant in the resonance Raman spectra recorded upon excitation at 454.5 and 647.2 nm, i.e., in the region of the absorption bands of the two lowest electronic ππ* transitions. These depolarized bands can be assigned to the non-totally symmetric B 1g and B 2g modes which are calculated at 1681 and 1242 cm −1 and involve skeletal vibrations of the central porphyrazine core with strong contribution from the N m -C α and N p -C α bonds formed by meso-and pyrrolic nitrogen atoms, respectively. In the spectrum obtained using excitation between the Soret and Q bands (λ ex = 514.5 nm), along with these depolarized bands several additional bands are enhanced. The polarized bands at 708, 865 and 1266 cm −1 should originate from totally symmetric modes. Indeed, calculations predict A 1g vibrations involving atoms constituting fused 1,2,5-thiadiazole rings at 735, 844 and 1413 cm −1 . The appearance of the anomalously polarized band at 1534 cm −1 might indicate the enhancement of the A 2g type vibration. Evidently, A 1g and A 2g vibrations are enhanced due to vibronic coupling with weak electronic transitions 2 1 E u and 3 1 E u (Table 5), which should be present near λ ex = 514.5 nm.
The metal free macrocycle H 2 TTDPz, due to the presence of two pyrrole and two pyrrolenine type fragments, has lower D 2h symmetry and its resonance Raman spectra (see Figure 10) are richer than for NiTTDPz (D 4h ). It can be seen that intense vibrations of the macrocyclic skeleton appearing as single bands for NiTTDPz are split into two components in the experimental spectra of H 2 TTDPz (1572 → 1566 and 1555; 1534 → 1532 and 1517; 1182 → 1177 and 1158; 709 → 711 and 697 cm −1 ). The intensity ratio of the components is dependent on the excitation wavelength. When excitation is shifted from the Q-to Soret band region, the bands containing a considerable contribution from the vibrational modes of the fused 1,2,5-thiadiazole rings are enhanced. This is not surprising since the Q-band transition is mainly localized on the atoms of the central macrocyclic core, while molecules with strong participation from the fused heterocycles participate in the electronic transitions in the Soret band region. Another remarkable feature in the experimental spectra of H 2 TTDPz is high intensity of the bands at 1430-1440, 1334 and 563 cm −1 , especially at the excitation at 454.5 nm.

Experimental Details
H 2 TTDPz was synthesized from the lithium(I) complex Li 2 TTDPz as described earlier [40]. NiTTDPz was prepared from H 2 TTDPz and nickel(II) acetate in DMSO at 100 • C following previously published methodology [42].
UV/Vis absorption spectra were recorded on a Cary 60 spectrometer. Resonance Raman spectra were recorded using Dilor XY multi-channel spectrometer and excitation by Ar + and Kr + lasers (Spectra Physics) for samples in KBr pellets at 80 K. FT Raman spectra were measured using IFX 66 CS/FRA 107 Bruker interferometer and excitation by NdYAG Atlas laser (Type 300, 1064 nm) at 10 K.
The sublimation of H 2 TTDPz and NiTTDPz has been investigated by Knudsen method using the MI-1201 commercial magnetic sector mass spectrometer adapted to thermodynamic studies and described in detail in [43,44]. The solid samples were sublimated from a stainless steel effusion cell. The ratio of the cross-sectional square of the cell to the square of the effusion orifice was about 1000, which made it possible to practically eliminate the violation of thermodynamic equilibrium due to the efflux of vapors from the effusion orifice in the experiments performed. The cell temperature was measured by a tungsten-rhenium thermocouple W-Re 5/20. The energy of ionizing electrons was 70 eV, accelerating voltage 5 kV, and cathode emission current I emis of 0.5 mA.

Computational Details
The electronic structure of NiTTDPz has been studied by the CASSCF method followed by accounting for dynamic electron correlation by multiconfigurational quasi-degenerate second-order perturbation theory (MCQDPT2). Eight electrons in five molecular orbitals consisting mainly of the 3d orbitals of Ni atom were selected for the active space. The doubly occupied orbitals corresponding to the 1s orbitals of C, N, S, and Ni and the 2s and 2p orbitals of S and Ni were frozen in the MCQDPT2 calculations. The triple-zeta basis sets pcseg-2 [45] from the Basis Set Exchange database [46,47] were used in all calculations.
DFT/PBE0-based investigations of H 2 TTDPz and NiTTDPz included geometry optimizations followed by computations of harmonic vibrations and TDDFT calculations of the electronic absorption spectra. The number of the calculated excited states was 30. The applicability of PBE0 functional for first row transition metal systems was shown by Jensen et al. [48]. All calculations were performed using the Firefly QC package [49], which is partially based on the GAMESS (US) [50] source code. The molecular models and orbitals demonstrated in the paper were visualized by means of the Chemcraft program [51]. Optimized Cartesian coordinates of H 2 TTDPz and NiTTDPz are available in the 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 [52].

Conclusions
Based on mass spectrometric investigations, it was established that NiTTDPz forms a stable stream of particles and enthalpy of sublimation was estimated using the Clausius-Clayperon equation. However, analysis of the mass spectrum showed that sublimation of H 2 TTDPz is accompanied by partial decomposition.
It has been found that H 2 TTDPz and NiTTDPz have a planar macrocycle structure of D 2h and D 4h symmetry, respectively. The complexes of nickel NiTTDPz possess the ground state 1 A 1g and the wave function has the form of a single determinant. It has been shown that for the studied complex NiTTDPz the crystal field theory (CFT) can be used to describe the sequences of the electronic states.
The effect of the introduction of the metal was studied on the basis of the results of TDDFT calculations. The electronic spectra of investigated molecules can be described by the four-orbital Gouterman model.
It has been shown that the theoretical data on the contribution of the vibrations of different fragments of the macrocycle to the normal modes are useful in the assignment of the experimental IR and Raman spectra.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.