Contrasting Photo-Switching Rates in Azobenzene Derivatives: How the Nature of the Substituent Plays a Role

A molecular design approach was used to create asymmetrical visible light-triggered azo-derivatives that can be good candidates for polymer functionalization. The specific electron–donor substituted molecules were characterized and studied by means of NMR analyses and UV-visible spectroscopy, comparing the results with Time Dependent Density Functional (TD-DFT) calculations. A slow rate of isomerization (ki = 1.5 × 10−4 s−1) was discovered for 4-((2-hydroxy-5methylphenyl) diazenyl)-3-methoxybenzoic acid (AZO1). By methylating this moiety, it was possible to unlock the isomerization mechanism for the second molecule, methyl 3-methoxy-4-((2-methoxy-5-methylphenyl) diazenyl)benzoate (AZO2), reaching promising isomerization rates with visible light irradiation in different solvents. It was discovered that this rate was heightened by one order of magnitude (ki = 3.1 × 10−3 s−1) for AZO2. A computational analysis using density functional (DFT/PBE0) and wavefunction (QD-NEVPT2) methodologies provided insight into the photodynamics of these systems. Both molecules require excitation to the second (S2) excited state situated in the visible region to initiate the isomerization. Two classic mechanisms were considered, namely rotation and inversion, with the former being energetically more favorable. These azo-derivatives show potential that paves the way for future applications as building blocks of functional polymers. Likewise, they could be really effective for the modification of existing commercial polymers, thus transferring their stimuli responsive properties to polymeric bulky structures, converting them into smart materials.


Introduction
A great attention during the past decades has been addressed to photochromic molecules in the field of materials science. These organic compounds can be included in many materials, changing the properties of the system when irradiated with light of the appropriate wavelength. The change can affect characteristics like: electronic or ionic conductivity, fluorescence, magnetism, shape, and conformation. The most studied photochromic molecules are azobenzenes [1] which have been widely investigated as promising photo-switchers [2][3][4][5][6] due to their simple preparation Rau has divided the azobenzenes into three spectroscopic classes based on their substitution patterns on the phenyl rings: azobenzene-type molecules, aminoazobenzene-type molecules, and pseudo-stilbenes, respectively [1]. The pseudo-stilbenes have presented the most efficient photoresponse. This behavior must be linked to a mixed photo-stationary state where the Z and E form are being constantly interconverted, caused by a shifted and overlapped absorption spectrum of the two isomers. Consequently, for these compounds, both the forward and reverse isomerization can be induced by a single wavelength of light in the visible region. On the other hand, the azobenzene-type molecules can isomerize very slowly from the Z configuration back into the E configuration when bulky substituents are introduced to delay the thermal relaxation [27,28]. Recently, the studies provided by Jerca et al. revealed that, using the proper solvent, it is possible to yield stable Z-isomers. This behavior can be enhanced by the nature and position of substituents on the aromatic rings [29]. This last feature of the azobenzene-type molecules is a very promising tool that allows fine-tuning of the optical response. Understanding the isomerization kinetics of different substituted azoderivatives in solution is of great importance in order to gain control over the response time to light irradiation and design photo responsive materials for precise applications. For instance, these azo systems can be used as moieties to bring stimuli-responsive polymers into a broad range of already Rau has divided the azobenzenes into three spectroscopic classes based on their substitution patterns on the phenyl rings: azobenzene-type molecules, aminoazobenzene-type molecules, and pseudo-stilbenes, respectively [1]. The pseudo-stilbenes have presented the most efficient photo-response. This behavior must be linked to a mixed photo-stationary state where the Z and E form are being constantly interconverted, caused by a shifted and overlapped absorption spectrum of the two isomers. Consequently, for these compounds, both the forward and reverse isomerization can be induced by a single wavelength of light in the visible region. On the other hand, the azobenzene-type molecules can isomerize very slowly from the Z configuration back into the E configuration when bulky substituents are introduced to delay the thermal relaxation [27,28]. Recently, the studies provided by Jerca et al. revealed that, using the proper solvent, it is possible to yield stable Z-isomers. This behavior can be enhanced by the nature and position of substituents on the aromatic rings [29]. This last feature of the azobenzene-type molecules is a very promising tool that allows fine-tuning of the optical response. Understanding the isomerization kinetics of different substituted azo-derivatives in solution is of great importance in order to gain control over the response time to light irradiation and design photo responsive materials for precise applications. For instance, these azo systems can be used as moieties to bring stimuli-responsive polymers into a broad range of already existing commercial polymeric families, thus transferring their stimuli responsive properties to polymeric bulky structures with the aim of enhancing materials performances. The idea of our work was to design, synthesize, and study novel asymmetrical azobenzene compounds with a precise substitutional pattern, bearing one reactive site. The specific properties derived from their substitutional pattern and the single reactive site could make them good candidates for the chemical modification of commercial polymers. The main goal was to create a proper photosensitive active group that, by means of its reactive site, could be grafted to a wide range of polymers, thus giving them stimuli-responsive properties. We focus our attention on the study of the E-Z-E isomerization kinetics of two related novel azo-derivatives (shown in Scheme 1). To the best of our knowledge, the synthesis of these two compounds has not been reported in the literature. These systems were well-designed to achieve an isomerization when irradiated by visible light and a thermal back isomerization when heated above room temperature. We investigated the conversion process in silico with a static approach (time dependent DFT) and the photodynamics of isomer conversion through a multi-configurational wavefunction methodology (QD-NEVPT2) in order to have a complete understanding of the properties of such systems. In the course of our investigation, some unusual and unexpected photochemistry of these asymmetrical ortho substituted derivatives were discovered. These findings can provide new insights regarding azobenzene isomerization as well as the possibility to reveal new applications for a well-known organic molecule. For instance, through a reactive site such as a carboxylic group, these photoactive moieties can be conveniently grafted to several classes of commercial polymers, such as polyvinyl alcohols (PVAs) and polyepichlorohydrin (PECH), thus conferring upon them light responsiveness [19]. In this way, commercial polymers could behave as novel smart materials.

Characterization
Subsequently, 1 H NMR and 13 C NMR spectra were recorded at 400 and 100.4 MHz, respectively, on a Varian Gemini 400 spectrometer (Varian, Inc. 3120 Hansen Way, Palo Alto, CA, USA) with proton noise decoupling for 13 C NMR. The chemical shifts were given in parts per million from TMS (Tetramethylsilane) in 1 H NMR spectra, while the central peak of the solvent was taken as the reference in the case of 13 C NMR spectra. The 13 C NMR spectra of the polymers were recorded at 30 • C, with a flip angle of 45 • , and the number of transients ranged from 20,000 to 40,000 with 10%-20% (w/v) sample solutions in deuterated dimethyl sulfoxide (DMSO-d6). A pulse delay time of 5 s for the 1 H NMR spectrum was used. Quantitative 13 C NMR spectra were performed using deuterated dimethyl sulfoxide (DMSO-d6) at 30 • C with a pulse delay time of 5 s. Photoisomerization was investigated using a UV-1800 Shimadzu Spectrophotometer (SHIMADZU EUROPA GmbH, Albert-Hahn-Strasse 6-10, 47269 Duisburg, Germany) with a scanning rate of 3000 to 2 nm/min. The kinetics of photoisomerization and thermal relaxation were followed by the same instrument. The samples for UV-visible spectroscopy were prepared via dilution of the sample into adequate solvent (DMSO, DMF, Dioxane) to give a final concentration of 7.4 × 10 −5 M in red light room (λ > 600 nm), to assure that they were predominantly in E form. After maintaining the azo dye solutions in the dark overnight, the spectra were taken in a  UV-grade quartz cell with optical beam path of 1 cm before and after irradiation with a desk lamp  equipped with a 11 W A E27 ID60 LED bulb (Jedi Lighting iDual bulb, TAO UK, 4 Deanfield Court,  Link 59 Business Park, Clitheroe, Lancs, UK). During the test the sample was placed at 2 cm from the light source. The optical signal output of white light was measured by an Ocean Optics USB2000 Miniature Fiber Optic Spectrometer (StellarNet, Inc., 14390 Carlson Circle, Tampa, FL, USA) in a dark room at room temperature, for more details the measured spectrum is shown in Figure S1 in the SI (3890 lux at 2 cm and 1030 lux at 8 cm, from the source, respectively). The solution was analyzed after 5 min irradiation. The process was repeated several times since the photostationary state (PPS) was observed. The Z-E thermal relaxation of solutions was studied after irradiating the samples reaching the photostationary state with the iDual lamp and then heating them in a TCC-100 Thermoelectrically Temperature Controlled cell holder (Aragonesses 13, Polig. Industrial 28100, Alcobendas, Madrid, Spain) for the analysis at different temperatures. The instruments collected the spectra each 2 min for 30 min. The temperatures used for the experiment were: 30, 40, 50, and 60 • C. Background was performed on each of the solvents used both before and after light irradiation and confirmed that the solvents absorption spectrum did not affect the region between 250 and 600 nm.

Synthesis of 4-((2-hydroxy-5-methylphenyl) diazenyl)-3-methoxybenzoic acid (AZO1) and
Methyl-3-methoxy-4-((2-methoxy-5-methylphenyl)diazinyl)benzoate (AZO2) For the synthesis of 4-((2-hydroxy-5-methylphenyl) diazenyl)-3-methoxybenzoic acid (AZO1), 12 mmol of 4-Amino-3-methoxybenzoic acid, 80 mL of water, and 7 mL of 35% hydrochloric acid were introduced in a 100 mL beaker and mildly heated at 40 • C to get the amine chlorohydrate (Scheme 1). The resulting suspension was then cooled to ≈ 0 • C in an ice bath. Meanwhile, 1.35 g (19.6 mmol) sodium nitrite were dissolved in 2.8 mL of water and added dropwise to the first suspension to obtain the corresponding diazonium salt (solution A). Solution B was prepared by dissolving 1.616 g (15 mmol) p-Cresol in 14 mL of water with 2.10 g (52.5 mmol) of sodium hydroxide. Solution B was then carefully added dropwise to solution A at 2-3 • C. A dark red precipitate was obtained. The crude product was then re-dissolved in 200 mL of a 10% aqueous solution of sodium hydroxide, the solution was diluted with 800 mL of water and then 20 mL 37% hydrochloric acid were added carefully. The product precipitated as a pure red solid which was then filtered and dried in oven at 45 • C; yield 58%, m.p.: 262.7 • C (deg.) Polymers 2019, 11, x FOR PEER REVIEW  4 of 21   overnight, the spectra were taken in a UV-grade quartz cell with optical beam path of 1 cm before  and after irradiation with a desk lamp equipped with a 11 W A E27 ID60 LED bulb (Jedi Lighting  iDual bulb,TAO UK, 4 Deanfield Court, Link 59 Business Park, Clitheroe, Lancs, UK). During the test the sample was placed at 2 cm from the light source. The optical signal output of white light was measured by an Ocean Optics USB2000 Miniature Fiber Optic Spectrometer (StellarNet, Inc., 14390 Carlson Circle, Tampa, FL, USA) in a dark room at room temperature, for more details the measured spectrum is shown in Figure S1 in the SI (3890 lux at 2 cm and 1030 lux at 8 cm, from the source, respectively). The solution was analyzed after 5 min irradiation. The process was repeated several times since the photostationary state (PPS) was observed. The Z-E thermal relaxation of solutions was studied after irradiating the samples reaching the photostationary state with the iDual lamp and then heating them in a TCC-100 Thermoelectrically Temperature Controlled cell holder (Aragonesses 13, Polig. Industrial 28100, Alcobendas, Madrid, Spain) for the analysis at different temperatures. The instruments collected the spectra each 2 min for 30 min. The temperatures used for the experiment were: 30, 40, 50, and 60 °C. Background was performed on each of the solvents used both before and after light irradiation and confirmed that the solvents absorption spectrum did not affect the region between 250 and 600 nm.

Synthesis of 4-((2-hydroxy-5-methylphenyl) diazenyl)-3-methoxybenzoic acid (AZO1) and Methyl-3methoxy-4-((2-methoxy-5-methylphenyl)diazinyl)benzoate (AZO2)
For the synthesis of 4-((2-hydroxy-5-methylphenyl) diazenyl)-3-methoxybenzoic acid (AZO1), 12 mmol of 4-Amino-3-methoxybenzoic acid, 80 mL of water, and 7 mL of 35% hydrochloric acid were introduced in a 100 mL beaker and mildly heated at 40 °C to get the amine chlorohydrate (Scheme 1). The resulting suspension was then cooled to ≈ 0 °C in an ice bath. Meanwhile, 1.35 g (19.6 mmol) sodium nitrite were dissolved in 2.8 mL of water and added dropwise to the first suspension to obtain the corresponding diazonium salt (solution A). Solution B was prepared by dissolving 1.616 g (15 mmol) p-Cresol in 14 mL of water with 2.10 g (52.5 mmol) of sodium hydroxide. Solution B was then carefully added dropwise to solution A at 2-3 °C. A dark red precipitate was obtained. The crude product was then re-dissolved in 200 mL of a 10% aqueous solution of sodium hydroxide, the solution was diluted with 800 mL of water and then 20 mL 37% hydrochloric acid were added carefully. The product precipitated as a pure red solid which was then filtered and dried in oven at 45°C; yield 58%, m.p.: 262.7 °C (deg.)  For the methylation of the hydroxyl groups of the previous molecule and the production of methyl-3-methoxy-4-((2-methoxy-5-methylphenyl) diazinyl) benzoate (AZO2), to 7.4 mmol of AZO1 in 50 mL of dry THF, 0.83 g of NaH were added (60% dispersion in mineral oil) followed by 1.73 mL (27 mmol) of methyl iodide. The mixture was stirred at 80 • C for 24 h. The solvent was evaporated with a rotavapor and the product was recrystallized from 2-propanol, obtaining an orange fine powder. Yield: 49%, m.p.: 138.4 • C.
The spectral characterization of both compounds is summarized below.

Computational Methods
Density functional calculations were performed using the ADF [30] 2017 103 modelling suite, employing the Perdew, Burke, and Ernzerhof hybrid functional (PBE0) [31][32][33]. The basis sets consisted of all electron triple zeta Slater type orbitals with one polarization function (TZP). The COnductor like Screening MOdel (COSMO) dielectric continuum solvation scheme was employed in the calculations with default appropriate parameters for dimethylformamide, dimethyl sulfoxide and dichloromethane [34]. Geometry optimizations did not involve any symmetry constraints. Time dependent (TD-PBE0) runs for the calculation of UV-visible transitions required 60 singlet-singlet excitations [35].
Additional multiconfigurational self-consistent field calculations [36] were performed using the ORCA 4.0.1 software package [37,38]. The PES scans were carried out by relaxing the geometries in vacuo at the state average (six singlets) CASSCF (4,3) level to a minimum at every scan point. At each of the scan points the wavefunction subsequently underwent a perturbation through the Quasi-Degenerate N-Electron Valence Perturbation Theory (QD-NEVPT2) approach [39,40] in the Nakano formulation [41]. The Ahlrichs def2-TZVP contracted Gaussian type basis set [42] was employed for every element making use of the resolution of the identity (RI) Coulomb integral approximation scheme with the def2-TZVP/C correlation density fitting basis set [43]. Using the state average CASSCF wavefunction, a non-SCF CI run was performed for the ground (S 0 ) and second excited (S 2 ) states which then underwent a density matrix analysis using the natural bond orbital scheme (NBO 6 program [44]). The natural localized molecular orbitals were generated, and their bond orders computed using the procedure devised by Schleyer et al. [45].

Molecular Design of the Azo-derivatives:
Azo compounds are usually synthesized with procedures like: azo coupling reaction (coupling of diazonium salts with activated aromatic compounds), Mills reaction (reaction between aromatic nitroso derivatives and anilines) and Wallach reaction (transformation of azoxybenzenes into 4-hydroxy substituted azo-derivatives in acid media) [46]. They are mostly produced with symmetric structures that have two potential reactive sites that offer the possibility to combine them with other monomers. In this way, polymers capable to react when irradiated by an electromagnetic irradiation can be produced. They are typically sensitive to UV light irradiation and they can have unstable Z form that can convert back to the more stable E form within picoseconds.
Polymers 2020, 12, 1019 6 of 20 A preliminary study on the possible substitution pattern for the novel azobenzene derivatives synthesized in this work was performed. The objective was to achieve a well-designed asymmetric structure to: • Produce a molecule with only one reactive site for further modifications and be used as a side-chain group.

•
Move the excitation wavelength of the E-Z isomerization from the UV region of the absorbance spectrum to the visible light one.

•
Enhance the stability of the Z form once obtained after irradiation with white light.
To successfully synthesize the first compound AZO1, an azo coupling reaction procedure was chosen with p-Cresol and 4-amino-3-methoxybenzoic acid as precursors. In this specific electrophilic aromatic substitution reaction, a phenol with a methyl group in para position was selected as the activated arene. The substitution occurs stereotypically in para to the activating group when a bulky substituent is selected as the attacking species. However, ortho substitution takes place when the para position is already occupied by another substituent. The 4-amino-3-methoxybenzoic acid was chosen to give the molecule the desired asymmetry and a reactive site for further modifications. In addition, with the presence of an electron donor group on the aromatic ring, like the methoxy one, we aimed to move the isomerization wavelength to the visible region of the absorption spectrum to trigger the isomerization with white light irradiation [47]. Lastly, it was reported that the methyl group could enhance the stability of the Z-isomer prolonging its half-life in the darkness [48].
The conditions of the reaction were carefully controlled both in pH and temperature due to the instability of the diazonium salts formed during the process. A mildly acidic or neutral solution is usually needed when an aromatic amine is involved in this reaction: nevertheless, if the solution is too acidic, the concentration of free amine becomes too small and the reaction does not occur. On the other hand, the phenol is usually not sufficiently active for the reaction, but, the more reactive phenoxide ion is a better initiator to generate the target azobenzene, when slightly alkaline conditions are present. Nonetheless, neither amines nor phenols can react in moderately alkaline solution because the diazonium ions are converted into the corresponding diazo hydroxide [49].
From the account given by Bandara et al. [50] regarding a specific set of azobenzene derivatives, the hydroxyl groups can hinder the isomerization by forming hydrogen bonds with the azo group, or other substituent on the opposite ring. We therefore envisaged to apply this same idea for the present work. The AZO1 compound was thus methylated to give AZO2 (Scheme 1). This modification generates a structure with two electron donor groups in ortho position on both the phenyl rings, changing the conjugation of the electrons on the -N=N-double bound. This modification was needed given the small isomerization conversion registered for AZO1 E-isomer when irradiated by white light, as shown in the following sections.
The novel azobenzene compounds were obtained with reasonably good yields and purities. Their structures were characterized by 1 H-NMR, 13 C-NMR, and UV-Vis spectroscopy.

H-NMR and 13 C-NMR Characterizations:
In Figure 2 1 H-NMR spectra of AZO1 and AZO2 molecules dissolved in deuterated DMSO are shown. The proton spectra exhibited some similarities. The peaks corresponding to the aromatic hydrogens are highlighted in the insets showing the δ region between 7.0 ppm and 8.0 ppm. The peaks related to the hydrogen atoms of the methoxy and methyl groups are present in the region from 2.0 to 4.0 ppm. A peak corresponding to the -OH of the carboxylic group is visible at 12.18 ppm (Figure 2A). After the conversion to methyl ester, this peak disappeared as it is evident in the Figure 2B. On the other hand, a peak corresponding to the methyl ester group is present at 3.91 ppm. All the peaks of the aromatic rings were shifted after the methylation, but their number remained unchanged as expected. The small peak at 3.33 ppm is characteristic of the residual water in d-DMSO. From these evidences, one can conclude that the methylation of the AZO1 molecule was successfully accomplished. AZO2 The 13 C-NMR chemical shifts of the molecules AZO1 and AZO2 are given in Table 1. From the data shown, it is evident that two new signals appear after the methylation of AZO1. They can be related to the two additional carbon atoms present in the AZO2 structure after AZO1 methylation. The 13 C-NMR chemical shifts of the molecules AZO1 and AZO2 are given in Table 1. From the data shown, it is evident that two new signals appear after the methylation of AZO1. They can be related to the two additional carbon atoms present in the AZO2 structure after AZO1 methylation.

UV/Vis Spectra Interpretation:
The UV-visible absorption spectra of the azo-derivatives exhibit two characteristic absorption bands. The first characteristic absorption band is present around 330-320 nm, while the second

UV/Vis Spectra Interpretation:
The UV-visible absorption spectra of the azo-derivatives exhibit two characteristic absorption bands. The first characteristic absorption band is present around 330-320 nm, while the second

UV/Vis Spectra Interpretation:
The UV-visible absorption spectra of the azo-derivatives exhibit two characteristic absorption bands. The first characteristic absorption band is present around 330-320 nm, while the second absorption band is present around 412-325 nm.
Comparing the absorption spectra, it is evident that the methoxy groups attached on both the phenyl rings of AZO2 are blue shifting the λ 1 of 26 nm in DMF (Table 2). This can be explained by the lower conjugation effect that they induce on the -N=N-bond electrons. Table 2. Wavelength maxima and experimentally measured rate constants for isomerization (k i ), reverse isomerization (k r ) and the content of the Z-isomer at the PSS (α Z ). AZO1 and AZO2 optical properties comparison. * = the obtained value is not reliable because it is falling within the experimental uncertainty of the instrument (see text).

Sample
Solvent To understand the optical properties of the two compounds, an in-silico analysis was carried out with the time dependent density functional (TD-DFT) approach. A DFT geometry optimisation (see computational methods) was performed for the ground state for both structural isomers of AZO1 and AZO2. The simulated electronic vertical transitions of both molecules are displayed on Figure 3 and these will occur in the frontier orbital region mapped out in Figure 4. The calculated bands (TD-PBE0) compare favourably well to the experimental ones, presenting a redshift in the region of 20-50 nm. The relative intensities are incorrectly described, likely due to solvent interaction effects.
The band maxima corresponding to the electronic excitations are listed in Tables 3 and 4 for AZO1 and AZO2, respectively. The most intense transition in the visible region can be assigned to the second excited state for both E-AZO1 and E-AZO2, each showing oscillator strengths of 0.57 and 0.91 respectively. The conventional assignment of distinctly separate n→ π* and π→ π* transitions is not rigorously valid in the case of E-AZO1 since there is an admixture of both transitions in the first two bands.   The band maxima corresponding to the electronic excitations are listed in Tables 3 and 4 for AZO1 and AZO2, respectively. The most intense transition in the visible region can be assigned to the second excited state for both E-AZO1 and E-AZO2, each showing oscillator strengths of 0.57 and 0.91 respectively. The conventional assignment of distinctly separate n→ π* and π→ π* transitions is not rigorously valid in the case of E-AZO1 since there is an admixture of both transitions in the first two bands.   Table 4. Calculated excitation energies, band maxima assignments and oscillator strengths (f osc ) of E-AZO2 (experimental values from UV/Vis analysis are shown in brackets).

Excited State # Main Character
The first and second electronic transitions of E-AZO1 correspond to two maxima calculated to be at 489 and 441 nm and have a majority of π 2 → π* character. The other band with significant oscillator strength in this compound is due to a π 1 → π* transition occurring at 338 nm.
The spectral features of E-AZO2 are altogether significantly different in that n→ π* transitions do not admix with the π→ π* type excitations. Since the former transitions are symmetry forbidden, they have a very low oscillator strength while the latter show up as very strong intensity bands. To broaden the study on the optical properties of the two moieties, the samples were dissolved in different solvents, according to the solubility of the different molecules. Figure 5 shows the dipole moments (in Debye) of each solvent compared to the wavelengths of the main peaks (λ 1 , λ 2 , λ 3 ) of the absorption spectra of the azo-derivatives.
For AZO1 molecule, a broader range of solvents could be tested. The first peculiar aspect is that AZO1 main peak (λ 1 ) is splitting its signal passing from a more polar solvent like dimethyl sulfoxide to a less polar one like dichloromethane. This can be ascribed to the different mechanisms of excitation followed by the electrons involved in the azo double bond when irradiated by light. In more polar solvents the signal linked to this excitation are overlapped (λ 1 ), then their wavelength split (λ 1 and λ 3 ) being better resolved when the polarity of the solvent decreases. The second crucial aspect is highlighted by the hypsochromic shift of λ 1 and the bathochromic shift of λ 2 of AZO1 molecule in polar solvents. This can be related to the stabilization of the electronic cloud of AZO1, that is favoured by the presence of a less polar solvent like dioxane. This, in fact, can decrease the conjugation effect of the electrons of the azo double bond, shifting lambda to lower wavelengths. AZO2 λ 1 is slightly influenced by the solvent polarity, a bathochromic shift is registered passing from a less polar solvent to a more polar one [51].
The UV-Vis spectrum was calculated for E-AZO1 in various implicit solvent dielectric media to determine how the effect that polarity shifts the visible band wavelengths. As can be seen from the values in Supplementary Table S4, the values of λ 1 tend to increase when passing from vacuum to the solvated environment.
Neglecting dispersion interactions between solvent and solute represents the limiting case of solvation in non-polar media. However, the value of λ 1 changed only a few for the used selection of solvents, with varying degree of polarity differences. The observed trend for DMSO is not reproduced by the calculation and in fact the value is slightly higher with respect to CH 2 Cl 2 and DMF, but the change is so small that barely any significance can be inferred from it. The failure to account for the experimental trend is a sign that explicit interactions between the solute and the solvent, which are not considered at this level of theory, are essential to explain the subtle band shifts. For AZO1 molecule, a broader range of solvents could be tested. The first peculiar aspect is that AZO1 main peak (λ1) is splitting its signal passing from a more polar solvent like dimethyl sulfoxide to a less polar one like dichloromethane. This can be ascribed to the different mechanisms of excitation followed by the electrons involved in the azo double bond when irradiated by light. In more polar solvents the signal linked to this excitation are overlapped (λ1), then their wavelength split (λ1 and λ3) being better resolved when the polarity of the solvent decreases. The second crucial aspect is highlighted by the hypsochromic shift of λ1 and the bathochromic shift of λ2 of AZO1 molecule in polar solvents. This can be related to the stabilization of the electronic cloud of AZO1, that is favoured by the presence of a less polar solvent like dioxane. This, in fact, can decrease the conjugation effect of the electrons of the azo double bond, shifting lambda to lower wavelengths. AZO2 λ1 is slightly influenced by the solvent polarity, a bathochromic shift is registered passing from a less polar solvent to a more polar one [51].
The UV-Vis spectrum was calculated for E-AZO1 in various implicit solvent dielectric media to determine how the effect that polarity shifts the visible band wavelengths. As can be seen from the values in Supplementary Table S4, the values of λ1 tend to increase when passing from vacuum to the solvated environment.
Neglecting dispersion interactions between solvent and solute represents the limiting case of solvation in non-polar media. However, the value of λ1 changed only a few for the used selection of solvents, with varying degree of polarity differences. The observed trend for DMSO is not reproduced by the calculation and in fact the value is slightly higher with respect to CH2Cl2 and DMF, but the change is so small that barely any significance can be inferred from it. The failure to account for the experimental trend is a sign that explicit interactions between the solute and the solvent, which are not considered at this level of theory, are essential to explain the subtle band shifts.

E-Z Photoisomerization
The photoisomerization kinetics of azo-derivatives was established by monitoring the

E-Z Photoisomerization
The photoisomerization kinetics of azo-derivatives was established by monitoring the absorbance of S 2 absorption band. Figure 6 shows the UV-Vis spectra of the AZO1 and AZO2 molecules after irradiation at different times. Under white light irradiation the intensity of the S 2 absorption band gradually decreases with the irradiation time, while the S 4 absorption band increases, until the photo-stationary state is reached. Both the solvent polarity and substituent nature influence the composition of the PSS. In fact, comparing the isomerization of both molecules, there are big differences that are easy to point out. First, only a slight change is visible in AZO1 absorption spectra ( Figure 6A). The E-Z conversion rate is about 3.3 % after 4 h 30 min of irradiation compared to 45.7% of AZO2 after 50 min of irradiation. The kinetic constant of the isomerization for the two molecules differ by an order of magnitude (Table 2).
influence the composition of the PSS. In fact, comparing the isomerization of both molecules, there are big differences that are easy to point out. First, only a slight change is visible in AZO1 absorption spectra ( Figure 6A). The E-Z conversion rate is about 3.3 % after 4 h 30 min of irradiation compared to 45.7% of AZO2 after 50 min of irradiation. The kinetic constant of the isomerization for the two molecules differ by an order of magnitude (Table 2).  Table 5 lists the photoisomerization rate constants calculated from UV-Vis experiments on AZO2 in different solvents at 303K. It can be observed that AZO2 ki increases going from a less polar solvent like dioxane to a more polar one, such as DMSO. This evidence can be related to the photoisomerization mechanism. When rotation or inversion-assisted rotation pathway is preferred, an acceleration of the photoisomerization rate in polar solvents is noticed, due to the stabilization of the dipolar transition state. A change in the isomerization mechanism can be responsible for the lower percentage of Z fraction registered when AZO2 is dissolved in dioxane [51].

Z-E Thermal Relaxation
The kinetic Z-E thermal relaxation of AZO1 and AZO2 were investigated at several temperatures in order to determine its thermodynamic parameters. However, in the case of AZO1, thermal relaxation was impossible to study due to the hindered isomerization. A rough estimation of the kinetic constant relative to Z-E thermal relaxation of AZO1 was attempted from the data obtained in DMF: unfortunately, the obtained value is not reliable because it is falling within the experimental uncertainty of the instrument (for sake of explanation the tentative fitting is shown in the SI with a value R 2 = 0.040). The rate of Z-E thermal relaxation of AZO2 molecule was monitored in the dark at the absorption maximum of its E-isomer. The thermal Z to E isomerization rate, summarized in Table  Table 5 lists the photoisomerization rate constants calculated from UV-Vis experiments on AZO2 in different solvents at 303K. It can be observed that AZO2 k i increases going from a less polar solvent like dioxane to a more polar one, such as DMSO. This evidence can be related to the photoisomerization mechanism. When rotation or inversion-assisted rotation pathway is preferred, an acceleration of the photoisomerization rate in polar solvents is noticed, due to the stabilization of the dipolar transition state. A change in the isomerization mechanism can be responsible for the lower percentage of Z fraction registered when AZO2 is dissolved in dioxane [51]. Table 5. AZO2 E-Z photoisomerization rate constants at 303K. In this order: k i , direct isomerization kinetic constant; k r , thermal relaxation kinetic constant; α Z , content of the Z-isomer at the PSS (photo-stationary state); τ 1/2 , half-life of Z-isomer in the dark.

Z-E Thermal Relaxation
The kinetic Z-E thermal relaxation of AZO1 and AZO2 were investigated at several temperatures in order to determine its thermodynamic parameters. However, in the case of AZO1, thermal relaxation was impossible to study due to the hindered isomerization. A rough estimation of the kinetic constant relative to Z-E thermal relaxation of AZO1 was attempted from the data obtained in DMF: unfortunately, the obtained value is not reliable because it is falling within the experimental uncertainty of the instrument (for sake of explanation the tentative fitting is shown in the SI with a value R 2 = 0.040). The rate of Z-E thermal relaxation of AZO2 molecule was monitored in the dark at the absorption maximum of its E-isomer. The thermal Z to E isomerization rate, summarized in Table 6, was determined for AZO2 at different temperatures and solvents. The process was found to be first order in all the studied solvents. 3.14 × 10 −6 2.14 × 10 −6 313. 15 1.18 × 10 −5 1.02 × 10 −5 6.23 × 10 −6 323. 15 3.22 × 10 −5 3.01 × 10 −5 1.73 × 10 −5 333. 15 5.59 × 10 −5 6.93 × 10 −5 3.81 × 10 −5 Calculation of the first-order rate constants at various temperatures allowed us to estimate thermodynamic activation parameters such as the activation energy (E a ), activation enthalpy (∆H ‡ ), and activation entropy (∆S ‡ ) for thermal Z to E reaction using Arrhenius and Eyring equations (equations shown in Appendix A). We used the linear correlation between the enthalpy and the entropy of activation to determine the mechanism of AZO2 thermal relaxation. In all the solvents, the Arrhenius plots were linear, proving that there is no significant change in the reaction pathway over the examined temperature range. AZO2 compound has a positive enthalpy of activation and negative entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ∆S ‡ absolute value and the lowest ∆H ‡ , and therefore exhibits the fastest thermal relaxation [52,53]. This is consistent with the calculation that will be shown below. In order to understand the dynamic photochemical and thermal processes involved in the isomerisation of AZO1 and AZO2 a choice of appropriate computational strategy was paramount. The employment of density functional theory is unsuitable in this instance because there is no guarantee, that at every stage in the transformation within the specific excited state hypersurface, the electron wavefunction of the ground state remains singly determinantal. The best methodology is to adopt multi-configurational self-consistent field theory [36,54] with a posteriori perturbational corrections for a good description of dynamical electron correlation such as the NEVPT2 approach [40,55] (see Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer change involve the either the twisting of the  6, was determined for AZO2 at different temperatures and solvents. The process was found to be first order in all the studied solvents. Calculation of the first-order rate constants at various temperatures allowed us to estimate thermodynamic activation parameters such as the activation energy (Ea), activation enthalpy (ΔH ‡ ), and activation entropy (ΔS ‡ ) for thermal Z to E reaction using Arrhenius and Eyring equations (equations shown in Appendix A). We used the linear correlation between the enthalpy and the entropy of activation to determine the mechanism of AZO2 thermal relaxation. In all the solvents, the Arrhenius plots were linear, proving that there is no significant change in the reaction pathway over the examined temperature range. AZO2 compound has a positive enthalpy of activation and negative entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute value and the lowest ΔH ‡ , and therefore exhibits the fastest thermal relaxation [52,53]. This is consistent with the calculation that will be shown below. In order to understand the dynamic photochemical and thermal processes involved in the isomerisation of AZO1 and AZO2 a choice of appropriate computational strategy was paramount. The employment of density functional theory is unsuitable in this instance because there is no guarantee, that at every stage in the transformation within the specific excited state hypersurface, the electron wavefunction of the ground state remains singly determinantal. The best methodology is to adopt multiconfigurational self-consistent field theory [36,54] with a posteriori perturbational corrections for a good description of dynamical electron correlation such as the NEVPT2 approach [40,55] (see Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer change involve the either the twisting of the (N=N-C) angle (inversion) or a rotation of the planes that make up the (C-N=N-C) dihedral (rotation). Only the latter was found to exhibit more favourable energies for the photochemical conversion. This contradicts the hydrogen bonding inhibition scenario as presented by Bandara et al. [50], since rotation of the molecule around the -CNNC-dihedral planes does not necessarily require the termination of the intramolecular hydrogen bond with the adjacent hydroxyl group. An alternative explanation is given below to account for the slower rates of photo-switching.
The alternative twisting scenario may be found in the supporting information section. More elaborate reaction coordinates combining angular and dihedral changes have also been studied for unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S2) presents the highest oscillator strength for both AZO1 and AZO2: this will be the most efficient conversion pathway.
A relaxed potential energy surface (PES) scan of the -CNNC-dihedral at the QD-NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of theory is presented in Figure 7 and in Figure 8 for AZO1 and AZO2, respectively. Therein the energies of each of the lowest six singlets were monitored throughout this reaction coordinate. The main composition of the states at the beginning of the (N=N-C) angle (inversion) or a rotation of the planes that make up the Polymers 2019, 11, x FOR PEER REVIEW 6, was determined for AZO2 at different temperatu order in all the studied solvents. Calculation of the first-order rate constants thermodynamic activation parameters such as the and activation entropy (ΔS ‡ ) for thermal Z to E (equations shown in Appendix A). We used the entropy of activation to determine the mechanism Arrhenius plots were linear, proving that there is n the examined temperature range. AZO2 compound entropy for all the solvents used (see Table 7).  [52,53]. This is consisten order to understand the dynamic photochemical an of AZO1 and AZO2 a choice of appropriate compu of density functional theory is unsuitable in this in stage in the transformation within the specific exc of the ground state remains singly determinan configurational self-consistent field theory [36,54] good description of dynamical electron correlat Computational Methods).
Two conventional mechanisms usually adopte the twisting of the (N=N-C) angle (inversion) o N=N-C) dihedral (rotation). Only the latter was f photochemical conversion. This contradicts the hy by Bandara et al. [50], since rotation of the molecu necessarily require the termination of the intramol group. An alternative explanation is given below t The alternative twisting scenario may be fou elaborate reaction coordinates combining angular unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excite for both AZO1 and AZO2: this will be the most eff A relaxed potential energy surface (PES) NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of th AZO1 and AZO2, respectively. Therein the energie throughout this reaction coordinate. The main c (C-N=N-C) dihedral (rotation). Only the latter was found to exhibit more favourable energies for the photochemical conversion. This contradicts the hydrogen bonding inhibition scenario as presented by Bandara et al. [50], since rotation of the molecule around the -CNNC-dihedral planes does not necessarily require the termination of the intramolecular hydrogen bond with the adjacent hydroxyl group. An alternative explanation is given below to account for the slower rates of photo-switching.
The alternative twisting scenario may be found in the supporting information section. More elaborate reaction coordinates combining angular and dihedral changes have also been studied for unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S 2 ) presents the highest oscillator strength for both AZO1 and AZO2: this will be the most efficient conversion pathway.
A relaxed potential energy surface (PES) scan of the -CNNC-dihedral at the QD-NEVPT2(4,3)// CASSCF (4,3)/def2-TZVP level of theory is presented in Figure 7 and in Figure 8 for AZO1 and AZO2, respectively. Therein the energies of each of the lowest six singlets were monitored throughout this reaction coordinate. The main composition of the states at the beginning of the process (E-isomer) is qualitatively described with different colours. The curves are qualitatively similar to what was previously reported for azobenzene [57]. process (E-isomer) is qualitatively described with different colours. The curves are qualitatively similar to what was previously reported for azobenzene [57].  The photochemical process originates with the promotion of S0S2 (π 1 π* 1 ) which can allow the system enough flexibility to distort the dihedral once the π bond strength has a formal bond order of zero. However, due to orbital polarisation there remains a marginally bonding interaction and some thermal energy is still required to break it. This barrier is computed to be 36.0 kJ mol −1 in the case of Calculation of the first-order rate constants at various temperatures allowed us to estimate thermodynamic activation parameters such as the activation energy (Ea), activation enthalpy (ΔH ‡ ), and activation entropy (ΔS ‡ ) for thermal Z to E reaction using Arrhenius and Eyring equations (equations shown in Appendix A). We used the linear correlation between the enthalpy and the entropy of activation to determine the mechanism of AZO2 thermal relaxation. In all the solvents, the Arrhenius plots were linear, proving that there is no significant change in the reaction pathway over the examined temperature range. AZO2 compound has a positive enthalpy of activation and negative entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute value and the lowest ΔH ‡ , and therefore exhibits the fastest thermal relaxation [52,53]. This is consistent with the calculation that will be shown below. In order to understand the dynamic photochemical and thermal processes involved in the isomerisation of AZO1 and AZO2 a choice of appropriate computational strategy was paramount. The employment of density functional theory is unsuitable in this instance because there is no guarantee, that at every stage in the transformation within the specific excited state hypersurface, the electron wavefunction of the ground state remains singly determinantal. The best methodology is to adopt multiconfigurational self-consistent field theory [36,54] with a posteriori perturbational corrections for a good description of dynamical electron correlation such as the NEVPT2 approach [40,55] (see Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer change involve the either the twisting of the (N=N-C) angle (inversion) or a rotation of the planes that make up the (C-N=N-C) dihedral (rotation). Only the latter was found to exhibit more favourable energies for the photochemical conversion. This contradicts the hydrogen bonding inhibition scenario as presented by Bandara et al. [50], since rotation of the molecule around the -CNNC-dihedral planes does not necessarily require the termination of the intramolecular hydrogen bond with the adjacent hydroxyl group. An alternative explanation is given below to account for the slower rates of photo-switching.
The alternative twisting scenario may be found in the supporting information section. More elaborate reaction coordinates combining angular and dihedral changes have also been studied for unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S2) presents the highest oscillator strength for both AZO1 and AZO2: this will be the most efficient conversion pathway.
A relaxed potential energy surface (PES) scan of the -CNNC-dihedral at the QD-NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of theory is presented in Figure 7 and in Figure 8 for AZO1 and AZO2, respectively. Therein the energies of each of the lowest six singlets were monitored throughout this reaction coordinate. The main composition of the states at the beginning of the (CNNC) = 139 • and some 8.6 kcal/mol above the S 2 minimum. process (E-isomer) is qualitatively described with different colours. The curves are qualitatively similar to what was previously reported for azobenzene [57].  The photochemical process originates with the promotion of S0S2 (π 1 π* 1 ) which can allow the system enough flexibility to distort the dihedral once the π bond strength has a formal bond order of zero. However, due to orbital polarisation there remains a marginally bonding interaction and some thermal energy is still required to break it. This barrier is computed to be 36.0 kJ mol −1 in the case of  Calculation of the first-order rate constants at various thermodynamic activation parameters such as the activation and activation entropy (ΔS ‡ ) for thermal Z to E reaction u (equations shown in Appendix A). We used the linear corr entropy of activation to determine the mechanism of AZO2 the Arrhenius plots were linear, proving that there is no significan the examined temperature range. AZO2 compound has a posit entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute value and the fastest thermal relaxation [52,53]. This is consistent with the ca order to understand the dynamic photochemical and thermal p of AZO1 and AZO2 a choice of appropriate computational stra of density functional theory is unsuitable in this instance beca stage in the transformation within the specific excited state hy of the ground state remains singly determinantal. The b configurational self-consistent field theory [36,54] with a pos good description of dynamical electron correlation such a Computational Methods).
Two conventional mechanisms usually adopted to describ the twisting of the (N=N-C) angle (inversion) or a rotation N=N-C) dihedral (rotation). Only the latter was found to ex photochemical conversion. This contradicts the hydrogen bon by Bandara et al. [50], since rotation of the molecule around necessarily require the termination of the intramolecular hydr group. An alternative explanation is given below to account fo The alternative twisting scenario may be found in the elaborate reaction coordinates combining angular and dihedr unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S2) for both AZO1 and AZO2: this will be the most efficient conve A relaxed potential energy surface (PES) scan of (CNNC) = 150 • and 6.6 kcal/mol above the S 2 minimum.
The photochemical process originates with the promotion of S 0 →S 2 (π 1 π* 1 ) which can allow the system enough flexibility to distort the dihedral once the π bond strength has a formal bond order of zero. However, due to orbital polarisation there remains a marginally bonding interaction and some thermal energy is still required to break it. This barrier is computed to be 36.0 kJ mol −1 in the case of AZO1 and 27.4 kcal mol −1 for AZO2 and takes place at Polymers 2019, 11, x FOR PEER REVIEW 6, was determined for AZO2 at different temperatures and solvents. The proce order in all the studied solvents. Calculation of the first-order rate constants at various temperatures a thermodynamic activation parameters such as the activation energy (Ea), acti and activation entropy (ΔS ‡ ) for thermal Z to E reaction using Arrhenius (equations shown in Appendix A). We used the linear correlation between entropy of activation to determine the mechanism of AZO2 thermal relaxation Arrhenius plots were linear, proving that there is no significant change in the the examined temperature range. AZO2 compound has a positive enthalpy of a entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute value and the lowest ΔH ‡ , and fastest thermal relaxation [52,53]. This is consistent with the calculation that w order to understand the dynamic photochemical and thermal processes involv of AZO1 and AZO2 a choice of appropriate computational strategy was param of density functional theory is unsuitable in this instance because there is no g stage in the transformation within the specific excited state hypersurface, the of the ground state remains singly determinantal. The best methodolog configurational self-consistent field theory [36,54] with a posteriori perturba good description of dynamical electron correlation such as the NEVPT2 Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer ch the twisting of the (N=N-C) angle (inversion) or a rotation of the planes t N=N-C) dihedral (rotation). Only the latter was found to exhibit more favou photochemical conversion. This contradicts the hydrogen bonding inhibition by Bandara et al. [50], since rotation of the molecule around the -CNNCdi necessarily require the termination of the intramolecular hydrogen bond with group. An alternative explanation is given below to account for the slower rat  Calculation of the first-order rate constants at v thermodynamic activation parameters such as the activ and activation entropy (ΔS ‡ ) for thermal Z to E reac (equations shown in Appendix A). We used the linea entropy of activation to determine the mechanism of AZ Arrhenius plots were linear, proving that there is no sig the examined temperature range. AZO2 compound has entropy for all the solvents used (see Table 7).  [52,53]. This is consistent wit order to understand the dynamic photochemical and th of AZO1 and AZO2 a choice of appropriate computation of density functional theory is unsuitable in this instanc stage in the transformation within the specific excited s of the ground state remains singly determinantal. configurational self-consistent field theory [36,54] with good description of dynamical electron correlation Computational Methods).
Two conventional mechanisms usually adopted to the twisting of the (N=N-C) angle (inversion) or a r N=N-C) dihedral (rotation). Only the latter was found photochemical conversion. This contradicts the hydrog by Bandara et al. [50], since rotation of the molecule ar necessarily require the termination of the intramolecula group. An alternative explanation is given below to acc respectively. One metric to measure the π symmetry bond strength from the electronic density is the natural localized molecular orbital (NLMO) bond order taken from the CASCI(4,3) density matrix. A comparative analysis of the azo bond orders in ground (S 0 ) and excited (S 2 ) states is displayed on Table 8. It may be seen that the calculated π bond order is identical from the outset (S 0 ). The S 0 →S 2 excitation suppresses the π bond strength differently in AZO1 and AZO2. In AZO1, the π bond remnant is stronger than in AZO2 thereby justifying the additional energetic requirement for the dihedral rotation and rendering AZO1 kinetically more inert in the S 2 state. Disregarding thermal contributions to the transition state energies, i.e., ∆E ‡ (S 2 ) ≈ ∆G ‡ (S 2 ) and assuming a high quantum yield for the S 2 excitation, the Eyring equation yields a k i (AZO1)/k i (AZO2) ratio of 3.1 × 10 −2 using the above calculated values which is reasonably close to experimental one (4.8 × 10 −2 ). Overall, a 9 kJ mol −1 energy difference leads to a decrease in the reaction rate constant close to two orders of magnitude.
Towards more acute angles one hallmark feature [57] of azobenzene derivatives comes to the fore, which is the conical intersection (CI) allowing S 2 to configurationally admix with S 0 leading to the cis isomer, in the region 80 • < 333. 15 5.59 × 10 6.93 × 10 3.81 × 10 Calculation of the first-order rate constants at various temperatures allowed us to estimat thermodynamic activation parameters such as the activation energy (Ea), activation enthalpy (ΔH ‡ and activation entropy (ΔS ‡ ) for thermal Z to E reaction using Arrhenius and Eyring equation (equations shown in Appendix A). We used the linear correlation between the enthalpy and th entropy of activation to determine the mechanism of AZO2 thermal relaxation. In all the solvents, th Arrhenius plots were linear, proving that there is no significant change in the reaction pathway ove the examined temperature range. AZO2 compound has a positive enthalpy of activation and negativ entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute value and the lowest ΔH ‡ , and therefore exhibits th fastest thermal relaxation [52,53]. This is consistent with the calculation that will be shown below. I order to understand the dynamic photochemical and thermal processes involved in the isomerisatio of AZO1 and AZO2 a choice of appropriate computational strategy was paramount. The employmen of density functional theory is unsuitable in this instance because there is no guarantee, that at ever stage in the transformation within the specific excited state hypersurface, the electron wavefunctio of the ground state remains singly determinantal. The best methodology is to adopt multi configurational self-consistent field theory [36,54] with a posteriori perturbational corrections for good description of dynamical electron correlation such as the NEVPT2 approach [40,55] (se Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer change involve the eithe the twisting of the (N=N-C) angle (inversion) or a rotation of the planes that make up the (C N=N-C) dihedral (rotation). Only the latter was found to exhibit more favourable energies for th photochemical conversion. This contradicts the hydrogen bonding inhibition scenario as presented by Bandara et al. [50], since rotation of the molecule around the -CNNC-dihedral planes does no necessarily require the termination of the intramolecular hydrogen bond with the adjacent hydroxy group. An alternative explanation is given below to account for the slower rates of photo-switching The alternative twisting scenario may be found in the supporting information section. Mor elaborate reaction coordinates combining angular and dihedral changes have also been studied fo unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S2) presents the highest oscillator strengt for both AZO1 and AZO2: this will be the most efficient conversion pathway.
A relaxed potential energy surface (PES) scan of the -CNNC-dihedral at the QD NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of theory is presented in Figure 7 and in Figure 8 fo AZO1 and AZO2, respectively. Therein the energies of each of the lowest six singlets were monitored throughout this reaction coordinate. The main composition of the states at the beginning of th (-CNNC-) < 100 • as the lone pairs (n) become degenerate with the π orbitals. The energy gap in this nonadiabatic coupling (H 02 ) is 98.3 kJ mol −1 in AZO1 and 83.1 kJ mol −1 in AZO2. From the Landau-Zener relation [58], it is known that the probability of transition between the two adiabatic surfaces S 0 and S 2 can be given by e −H 2 02 /α where α is a collection of constants incorporating nuclear velocity, and the difference in reaction coordinate slope in each adiabatic state. H 2 02 is the squared energy gap between the two states at the avoided crossing. As the values of H 02 are considerable for both molecules (98.3 kJ mol −1 in AZO1 and 83.1 kJ mol −1 in AZO2), the probability of radiationless transition is practically zero. Thus, the S 0 ←S 2 transition will require the emission of one photon at 1217 nm (AZO1) and 1437 nm (AZO2) wavelengths (infra-red region).
The potential energy landscapes also afford an overview of the thermal back-isomerization which has the values +74.9 and +74.4 kJ mol −1 . These values are in good agreement with the relaxation activation energies determined experimentally even neglecting the necessary thermal corrections (zero point, vibrational, rotational and translational).
In light of the current data, the presence of the -OH group in AZO1 cannot be considered to play a direct role in the rapid back isomerization of AZO1 if one considers this neutral tautomer alone as alluded to in the literature [50]. In view of the acidic properties of AZO1, all the possible tautomeric zwitterionic isomers were optimized at the PBE0/TZP level and their optical properties assessed (see Supporting Information section).
Our results show that a low energy tautomer may play a role in the isomerization process if the following equilibrium occurs (Scheme 2). This zwitterionic isomer should be more favored in a polar solvent medium [∆E(DMF) = 3.6 kcal mol −1 ]: Our results show that the presence of the hydroxy group does indeed influence the conversion process in AZO1, but not by creating strong hydrogen bonds with nitrogen but by instead allowing for the existence of diazo-phenolate that has distinctly unique photodynamic properties (Supplementary Figure S6).
alluded to in the literature [50]. In view of the acidic properties of AZO1, all the possible tautomeric zwitterionic isomers were optimized at the PBE0/TZP level and their optical properties assessed (see Supporting Information section).
Our results show that a low energy tautomer may play a role in the isomerization process if the following equilibrium occurs (Scheme 2). This zwitterionic isomer should be more favored in a polar solvent medium [ΔE(DMF) = 3.6 kcal mol −1 ]: Scheme 2. Tautomeric equilibrium between two structural isomers (A, C) of AZO1. Scheme 2. Tautomeric equilibrium between two structural isomers (A,C) of AZO1.

Conclusions
In this work, two novel asymmetrical azobenzene derivatives were successfully synthesized, characterized, and studied. A molecular design was done in order to plan and create moieties with one reactive suitable site for polymer functionalization. The molecules were set to have electron-donor groups that permit isomerization when irradiated by visible light irradiation and a sufficiently stable Z-isomer with a thermal relaxation at temperatures above 303K. A wide-ranging study on the E-Z-E isomerization mechanism of the azo-derivatives was performed. The definite substitutional pattern was found to change the specific photochemistry of the molecules. A hindered isomerization was discovered in the case of AZO1 molecule, subsequently analyzed by the TD-DFT calculations.
On the other hand, by methylating AZO1 to AZO2, the photochemistry of AZO2 was found to follow the bibliography trends. AZO2 isomerization occurs under visible light irradiation, giving a stable Z-isomer at room temperature, with thermal relaxation rate constants increasing at higher temperature, as expected. The interaction between solvent and chromophore plays a role in determining the rate of this process and therefore should be considered in future applications.
By means of TD-DFT calculations, the molecules were studied and compared. The calculated bands (TD-PBE0) relate well to the experimental ones, presenting a minor systematic redshift. The conventional assignment of distinctly separated n→ π* and π→ π* transitions was found to be not rigorously valid in the case of E-AZO1 since an admixture of both transitions in the first two bands was found. The spectral features of E-AZO2 are different in that n→ π* transitions do not admix with the π→ π* type excitations.
The dynamic photochemical and thermal processes involved in the isomerisation of AZO1 and AZO2 was studied by means of multi-configurational self-consistent field theory with a posteriori perturbational corrections. The rotation mechanism of the planes that make up the Polymers 2019, 11, x FOR PEER REVIEW 6, was determined for AZO2 at different temperatur order in all the studied solvents. Calculation of the first-order rate constants thermodynamic activation parameters such as the and activation entropy (ΔS ‡ ) for thermal Z to E (equations shown in Appendix A). We used the l entropy of activation to determine the mechanism o Arrhenius plots were linear, proving that there is n the examined temperature range. AZO2 compound entropy for all the solvents used (see Table 7). AZO2 in DMF has the highest ΔS ‡ absolute val fastest thermal relaxation [52,53]. This is consistent order to understand the dynamic photochemical an of AZO1 and AZO2 a choice of appropriate computa of density functional theory is unsuitable in this ins stage in the transformation within the specific excit of the ground state remains singly determinan configurational self-consistent field theory [36,54] good description of dynamical electron correlati Computational Methods).
Two conventional mechanisms usually adopte the twisting of the (N=N-C) angle (inversion) or N=N-C) dihedral (rotation). Only the latter was fo photochemical conversion. This contradicts the hy by Bandara et al. [50], since rotation of the molecu necessarily require the termination of the intramole group. An alternative explanation is given below to The alternative twisting scenario may be fou elaborate reaction coordinates combining angular unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited for both AZO1 and AZO2: this will be the most effi A relaxed potential energy surface (PES) NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of th AZO1 and AZO2, respectively. Therein the energies throughout this reaction coordinate. The main co (C-N=N-C) dihedral angle was found to be the most favourable for the E-Z isomer change. The photo-chemical process originates with the electronic promotion of S 0 →S 2 and an energetic barrier in the S 2 surface was found along the reaction coordinate valued at +36.0 and +27.4 kJ mol −1 for AZO1 and AZO2 respectively. This difference in thermal barriers for the S 2 state explains why AZO1 is kinetically more inert to isomerisation with respect to AZO2. Another route for AZO1 conversion is through dihedral rotation of the proton tautomer but the ground state of the Z isomer in this instance is thermodynamically unstable and rapidly switches back to the E isomer. This is consistent with the spectral evidence.
While the topology of the hydrogen bonds in the ortho position may be relevant in the photo-switching of certain azo derivatives, this feature is entirely reliant on the nature of the most favourable reaction coordinate. In light of previous results [53], this shows that subtle interplays between substituent groups and the aryl rings pose difficulties in the modular design of tuneable azo compounds since their interconversion shows up as system specific.
The photo-switching in AZO1 is inhibited due to a stronger π bond order hindering rotation around the N-N bond axis. Further studies will be required to attempt to draw systematic functional group trends within the same mechanistic scenario.
To conclude, in this work a molecular design approach was used to plan and produce two novel moieties, suitable candidates for applications that aim to take advantage of visible light irradiation. Combined with TD-DFT calculation, a comprehensive study of their properties was achieved giving also more information about the important role of the substituents that must be selected during the planning phase of these molecules. These findings can provide new insights regarding substituted azobenzene isomerization as well as the possibility to reveal new applications for a well-known organic molecule.

Solvent Dipole Moment (D)
A (s −1 ) Ea (kJ mol −1 ) ΔH ‡ (kJ mol −1 ) ΔS ‡ (J mol −1 K −1 ) AZO2 in DMF has the highest ΔS ‡ absolute value and the lowest ΔH ‡ , and therefore exhibits the fastest thermal relaxation [52,53]. This is consistent with the calculation that will be shown below. In order to understand the dynamic photochemical and thermal processes involved in the isomerisation of AZO1 and AZO2 a choice of appropriate computational strategy was paramount. The employment of density functional theory is unsuitable in this instance because there is no guarantee, that at every stage in the transformation within the specific excited state hypersurface, the electron wavefunction of the ground state remains singly determinantal. The best methodology is to adopt multiconfigurational self-consistent field theory [36,54] with a posteriori perturbational corrections for a good description of dynamical electron correlation such as the NEVPT2 approach [40,55] (see Computational Methods).
Two conventional mechanisms usually adopted to describe the isomer change involve the either the twisting of the (N=N-C) angle (inversion) or a rotation of the planes that make up the (C-N=N-C) dihedral (rotation). Only the latter was found to exhibit more favourable energies for the photochemical conversion. This contradicts the hydrogen bonding inhibition scenario as presented by Bandara et al. [50], since rotation of the molecule around the -CNNC-dihedral planes does not necessarily require the termination of the intramolecular hydrogen bond with the adjacent hydroxyl group. An alternative explanation is given below to account for the slower rates of photo-switching.
The alternative twisting scenario may be found in the supporting information section. More elaborate reaction coordinates combining angular and dihedral changes have also been studied for unsubstituted azobenzene [56].
As seen from Tables 3 and 4 the second excited state (S2) presents the highest oscillator strength for both AZO1 and AZO2: this will be the most efficient conversion pathway.
A relaxed potential energy surface (PES) scan of the -CNNC-dihedral at the QD-NEVPT2(4,3)//CASSCF (4,3)/def2-TZVP level of theory is presented in Figure 7 and in Figure 8 for AZO1 and AZO2, respectively. Therein the energies of each of the lowest six singlets were monitored throughout this reaction coordinate. The main composition of the states at the beginning of the (NNC) angle twisting to the linear position, Figure S3-Tight grid scan of AZO1 (left) and AZO2 (right) evidencing the two (avoided) crossings of S 2 and S 3 , Table S1-Relative energies of the zwitterionic isomers of AZO1 as calculated through PBE0/TZP (DMF); Table S2-Calculated UV/Vis absorption spectrum of E-AZO1-C at the PBE0/TZP (DMF) level of theory; Figure S4-Calculated UV-Vis spectra of the two tautomers (convoluted with a Gaussian function at FWHM=60); Figure S5 Table S4-Calculated band absorption wavelengths for the AZO1 system; Figure S8