Tuning the Photophysics of Two-Arm Bis[(dimethylamino)styryl]benzene Derivatives by Heterocyclic Substitution

The identification of novel molecular systems with high fluorescence and significant non-linear optical (NLO) properties is a hot topic in the continuous search for new emissive probes. Here, the photobehavior of three two-arm bis[(dimethylamino)styryl]benzene derivatives, where the central benzene was replaced by pyridine, furan, or thiophene, was studied by stationary and time-resolved spectroscopic techniques with ns and fs resolution. The three molecules under investigation all showed positive fluorosolvatochromism, due to intramolecular charge-transfer (ICT) dynamics from the electron-donor dimethylamino groups, and significant fluorescence quantum yields, because of the population of a planar and emissive ICT state stabilized by intramolecular hydrogen-bond-like interactions. The NLO properties (hyperpolarizability coefficient and TPA cross-section) were also measured. The obtained results allowed the role of the central heteroaromatic ring to be disclosed. In particular, the introduction of the thiophene ring guarantees high fluorescent quantum yields irrespective of the polarity of the medium, and the largest hyperpolarizability coefficient because of the increased conjugation. An important and structure-dependent involvement of the triplet state was also highlighted, with the intersystem crossing being competitive with fluorescence, especially in the thiophene derivative, where the triplet was found to significantly sensitize molecular oxygen even in polar environment, leading to possible applications in photodynamic therapy.

These push-pull systems show large charge displacement during the excitation and the presence of intramolecular charge transfer (ICT) in the excited state particularly favored in polar solvents.
Recently we have studied the photobehavior of distyrylbenzene analogs in which the central benzene was replaced by a heteroaromatic ring (furan, pyridine, or thiophene), and the lateral benzenes functionalized with strong electron-withdrawing nitro groups in the para position [25,26].

Solvent Effect on Spectral and Fluorescence Properties
The spectral characterization of the newly synthetized dimethylamino derivatives in many solvents of different polarity and polarizability is shown in Figure 1, while the molar absorption coefficients measured in Tol are reported in Table S1 and are all of the order of 30,000 M −1 cm −1 . The absorption spectrum shifted to the red at about 60 to 70 nm when replacing the π-deficient pyridine unit with the π-rich furan and thiophene rings, respectively, due to the increased degree of conjugation [25,28]. DMA-QP showed an absorption maximum at around 370 nm, with a symmetric bell-shaped spectral profile; while a more A detailed experimental study using conventional steady-state and time-resolved (with ns and fs resolution) spectroscopic techniques coupled with theoretical calculations at the TD-DFT level was performed. The combined experimental and theoretical study allowed the competitive processes deactivating the excited electronic states and their mechanisms to be pointed out, highlighting the effect given by the central heteroaromatic ring (furan, pyridine, and thiophene). The spectral characterization of the three compounds in solvents of different polarity to derive their hyperpolarizability coefficient through a solvatochromic method is also reported, together with their two-photon absorption (TPA) cross-sections measured by the two-photon excited fluorescence technique.
Particular attention was paid to the comparison of the photobehavior of the compounds investigated of the present study, with the symmetric dinitro-derivative analogs and the asymmetric ones that carry both the dimethylamino group on one edge and the nitro group on the opposite side, with the intent to investigate how the replacement of strong electron-withdrawing nitro substituents with strong electron-donor dimethylamino ones can tune the photobehavior and NLO properties in these symmetrical systems.

Solvent Effect on Spectral and Fluorescence Properties
The spectral characterization of the newly synthetized dimethylamino derivatives in many solvents of different polarity and polarizability is shown in Figure 1, while the molar absorption coefficients measured in Tol are reported in Table S1 and are all of the order of 30,000 M −1 cm −1 . The absorption spectrum shifted to the red at about 60 to 70 nm when replacing the π-deficient pyridine unit with the π-rich furan and thiophene rings, respectively, due to the increased degree of conjugation [25,28]. DMA-QP showed an absorption maximum at around 370 nm, with a symmetric bell-shaped spectral profile; while a more structured shape was found in the cases of DMA-QF and DMA-QT, for which a bathochromic shoulder was also observed. The absorption spectra of the investigated compounds were barely affected by the solvent polarity, whereas the polarizability was responsible for the small blue-shift of the spectrum when dealing with low refractive index media, with the most hypsochromic absorption being recorded in EtOAc.
Contrarily, a marked fluorosolvatochromism emerged from the fluorimetric investigation, suggesting a higher dipolar character of the excited state when going towards more polar media [46]. The different nature of the emissive excited state, acquiring an intramolecular charge transfer (ICT) character, is also supported by the clear change of the emission profile, which features a net vibronic structure in Tol, particularly for DMA-QF and DMA-QT, and turning into a bell-like fluorescence spectrum in more polar environments. The fluorosolvatochromic behavior was mainly evidenced in the case of the furan derivative, with the emission maximum shifting to the red at about 130 nm moving from Tol to EtOH, as compared to 50 and ca. 30 nm found for DMA-QP and DMA-QT, respectively (see Table 1 below). Simultaneously, an increase of the Stokes shift and full width at half maximum (FWHM) upon increasing the solvent polarity was underscored, accompanied for DMA-QP and DMA-QF by a not negligible fluorescence quantum yield (Φ F = 0.07 and 0.06 in Tol, respectively), which increased in polar media (Φ F = 0.20 and 0.37 in AcCN, respectively, see Figure 2 below and Table S1). The only exception to this Φ F trend was observed in the protic EtOH (Φ F = 0.06 for DMA-QP and 0.11 for DMA-QF), for which the most bathochromically shifted emission was also found. In fact, the fluorescence quantum yield was halved compared to its value in DMF and was more than three times smaller than in AcCN for both compounds. This result could be linked to the presence of intramolecular hydrogen-like bonds between the central heteroatom and the ethylene hydrogens, which stabilize a fluorescent planar compressed form (see below, Section 2.2). In a protic solvent, intermolecular H-bonds loosen the intramolecular interactions, favoring less compressed and more distorted geometries, in which non-radiative deactivation pathways predominate [47,48]. The presence of thermal equilibria between many slightly different conformations in fluid media was also in line with the much broader absorption spectrum recorded for DMA-QF in EtOH, where each conformation contributed to the absorption with its distinct absorption profile (7400 cm −1 for the FWHM vs. 5900 cm −1 in DMF). The highest fluorescence quantum yields (Φ F~0 .5) were measured in the case of DMA-QT, likely because of the more conjugated and planar structure conferred by the central less-aromatic thienyl ring. The emission capability of the thiophene-substrate, unlike the previous compounds, remained constant, being insensitive either to solvent polarity, polarizability, or even proticity, in line with the reduced tendency of the sulfur atom, as compared to nitrogen and oxygen, to form H-type bonds with solvent molecules. The weaker fluorescence of the first two systems in the series suggests that non-radiative pathways play a crucial role in their excited-state photodynamics (see below).  Contrarily, a marked fluorosolvatochromism emerged from the fluorimetric investigation, suggesting a higher dipolar character of the excited state when going towards more polar media [46]. The different nature of the emissive excited state, acquiring an intramolecular charge transfer (ICT) character, is also supported by the clear change of the emission profile, which features a net vibronic structure in Tol, particularly for DMA-QF and DMA-QT, and turning into a bell-like fluorescence spectrum in more polar environments. The fluorosolvatochromic behavior was mainly evidenced in the case of the furan derivative, with the emission maximum shifting to the red at about 130 nm moving

Quantum Mechanical Calculations
Notoriously, these double-arm molecules could exist in fluid solution at room temperature as dynamic equilibria among conformers originating by rotation around the quasi-single bonds between the heteroaromatic ring and the ethenyl bridge. DFT calculations aided in carrying out a complete theoretical study on the conformational equilibrium of the three DMA derivatives, in order to distinguish between the most abundant conformers likely present in solution-either the compressed form (dipolar-like), the semielongated or the elongated species (quadrupolar-like). From the optimizations of the ground-and excited-state geometries, the compressed form was found to largely prevail (see Figure 3 and Table S3) [49], and was also experimentally attested by the lack of wavelength effect on emission and fluorescence excitation spectra (see Figure S1 as a representative example), being stabilized by the intramolecular H-like bonds between the hydrogens of the adjacent ethenyl bridges and the central electronegative heteroatom (Figure 3). As a matter of fact, the intramolecular N⋯H/O⋯H distances of 2.54 and 2.56 Å, respectively, and the NĤC-OĤC angles (97.8° and 96.9°, respectively) demonstrate the Hbond type nature of this interaction and, thus, the prevalence of the bent structures (Table  S4). For this reason, the pyridine and furan derivatives will be considered from now on to be dipolar-like push-pull systems with the dipole moment oriented along the y-axis. Adversely, the larger S-H distance (2.83 Å) and wider CĤS angle (102.9°), jointly with the more extended optimized structure, led to the interpretation of DMA-QT as a quadrupolar molecule (Table S4). In addition, the predicted excited-state dipole moments (Table  S5), calculated through the Mulliken atomic charges and increasing upon excitation, showed the largest difference (ΔμCT) for the pyridine derivative and the lowest values for the ground and excited state in the case of the thiophene-substituted compound. In contrast, the QXX component of the quadrupole moment was found to slightly decrease upon

Quantum Mechanical Calculations
Notoriously, these double-arm molecules could exist in fluid solution at room temperature as dynamic equilibria among conformers originating by rotation around the quasi-single bonds between the heteroaromatic ring and the ethenyl bridge. DFT calculations aided in carrying out a complete theoretical study on the conformational equilibrium of the three DMA derivatives, in order to distinguish between the most abundant conformers likely present in solution-either the compressed form (dipolar-like), the semielongated or the elongated species (quadrupolar-like). From the optimizations of the ground-and excited-state geometries, the compressed form was found to largely prevail (see Figure 3 and Table S3) [49], and was also experimentally attested by the lack of wavelength effect on emission and fluorescence excitation spectra (see Figure S1 as a representative example), being stabilized by the intramolecular H-like bonds between the hydrogens of the adjacent ethenyl bridges and the central electronegative heteroatom (Figure 3). As a matter of fact, the intramolecular N· · · H/O· · · H distances of 2.54 and 2.56 Å, respectively, and the NĤC-OĤC angles (97.8 • and 96.9 • , respectively) demonstrate the H-bond type nature of this interaction and, thus, the prevalence of the bent structures (Table S4). For this reason, the pyridine and furan derivatives will be considered from now on to be dipolar-like push-pull systems with the dipole moment oriented along the y-axis. Adversely, the larger S-H distance (2.83 Å) and wider CĤS angle (102.9 • ), jointly with the more extended optimized structure, led to the interpretation of DMA-QT as a quadrupolar molecule (Table S4). In addition, the predicted excited-state dipole moments (Table S5), calculated through the Mulliken atomic charges and increasing upon excitation, showed the largest difference (∆µ CT ) for the pyridine derivative and the lowest values for the ground and excited state in the case of the thiophene-substituted compound. In contrast, the Q XX component of the quadrupole moment was found to slightly decrease upon excitation in the case of DMA-QP and DMA-QF, while it halved in the case of DMA-QT (Table S5). TD-DFT calculations (see Tables S6-S17) managed to successfully reproduce the experimental absorption spectra, whose main transition (S 0 → S 1 ) is significantly described by the π-π* HOMO-LUMO configuration (>70% for DMA-QP and >80% for DMA-QF and DMA-QT). The corresponding molecular orbitals (see Tables S7, S11, S14 and S16) showed a higher degree of ICT character for the relaxed excited state, which corroborates the occurrence, during the excited-state lifetime, of some charge displacement involving the central heterocycle and the ethenyl groups. This charge movement was more distinct in the case of the pyridine derivative, being a π-poor heterocycle in comparison to the furan and thiophene units. In addition, in the electron-density difference maps ( Figures S2-S5), a consistent weakening of the double bonds is observed in the compressed DMA-QP when going from the S 0 to S 1 state, which likely favors trans-cis photoisomerization in the singlet manifold. As for DMA-QF, electron density depletion on the furan unit, with a consequent decrease in its aromaticity, is instead observed.

NLO Properties
The solvatochromic behavior showed by the three compounds under investigation was exploited in order to estimate second-order non-linear optical properties appealing for optoelectronic and photonics applications, quantified through the hyperpolarizability coefficients. The solvatochromic method has been successfully applied for push-pull compounds, proving to be a valid alternative and simplified way to explore those NLO properties without relying on more sophisticated techniques such as EFISH (electron-field induced second-harmonic) and HRS (hyper-Rayleigh scattering) [12,50,51]. From the slope given by the linear plot of the Stokes shift vs. the solvent properties ( Figure 4), namely the dielectric constant and refractive index, derived by Equation (1), the experimental dipolemoment difference along the charge-transfer direction, ΔμCT, between the first excited singlet state and the ground state was determined. The Stokes shift recorded in EtOH was excluded from the plot for both DMA-QP and DMA-QF because of the intramolecular Htype bonds with the solvent affecting spectral positions beside solvatochromism. The ΔμCT values match the fluorosolvatochromic behavior of the DMA derivatives; the largest difference was found for the furan-derivative, followed by the pyridine-and thiophene-substituted compounds. This information was then manipulated using Oudar's formula (Equation (2)) for the pyridine and furan derivatives and in the modified equation (Equation (3)) for the quadrupolar DMA-QT compound, to give the dynamic hyperpolarizability coefficients (βCT) in Tol (see Table 2). The sizeable βCT values (248, 924, and 1030 esu −1 cm 5 for DMA-QP, DMA-QF, and DMA-QT, respectively) were found to be comparable to other NLO materials and their trend reproduces the increased conjugation when replacing the central pyridine with the furan and, particularly, with the thienyl unit, as found for other push-pull molecules [51].

NLO Properties
The solvatochromic behavior showed by the three compounds under investigation was exploited in order to estimate second-order non-linear optical properties appealing for optoelectronic and photonics applications, quantified through the hyperpolarizability coefficients. The solvatochromic method has been successfully applied for push-pull compounds, proving to be a valid alternative and simplified way to explore those NLO properties without relying on more sophisticated techniques such as EFISH (electron-field induced second-harmonic) and HRS (hyper-Rayleigh scattering) [12,50,51]. From the slope given by the linear plot of the Stokes shift vs. the solvent properties ( Figure 4), namely the dielectric constant and refractive index, derived by Equation (1), the experimental dipole-moment difference along the charge-transfer direction, ∆µ CT , between the first excited singlet state and the ground state was determined. The Stokes shift recorded in EtOH was excluded from the plot for both DMA-QP and DMA-QF because of the intramolecular H-type bonds with the solvent affecting spectral positions beside solvatochromism. The ∆µ CT values match the fluorosolvatochromic behavior of the DMA derivatives; the largest difference was found for the furan-derivative, followed by the pyridine-and thiophene-substituted compounds. This information was then manipulated using Oudar's formula (Equation (2)) for the pyridine and furan derivatives and in the modified equation (Equation (3)) for the quadrupolar DMA-QT compound, to give the dynamic hyperpolarizability coefficients (β CT ) in Tol (see Table 2). The sizeable β CT values (248, 924, and 1030 esu −1 cm 5 for DMA-QP, DMA-QF, and DMA-QT, respectively) were found to be comparable to other NLO materials and their trend reproduces the increased conjugation when replacing the central pyridine with the furan and, particularly, with the thienyl unit, as found for other push-pull molecules [51].
values match the fluorosolvatochromic behavior of the DMA derivatives; the largest difference was found for the furan-derivative, followed by the pyridine-and thiophene-substituted compounds. This information was then manipulated using Oudar's formula (Equation (2)) for the pyridine and furan derivatives and in the modified equation (Equation (3)) for the quadrupolar DMA-QT compound, to give the dynamic hyperpolarizability coefficients (CT) in Tol (see Table 2). The sizeable βCT values (248, 924, and 1030 esu −1 cm 5 for DMA-QP, DMA-QF, and DMA-QT, respectively) were found to be comparable to other NLO materials and their trend reproduces the increased conjugation when replacing the central pyridine with the furan and, particularly, with the thienyl unit, as found for other push-pull molecules [51].    Table 2. Predicted and experimental data required for the determination of the dynamic (β CT ) and static (β 0 ) hyperpolarizability coefficients of the investigated compounds in Tol: ground state dipole moment (µ GS ); excited-state dipole moment (µ ES ); Onsager's cavity ray (a); dipole moment difference (∆µ CT = µ ES − µ GS ); absorption maximum frequency ( v eg ); theoretical oscillator strength (f ). In the case of DMA-QT, the distance between the barycenter of the two poles (d = 8.00 Å) is also used as the double dipole model is considered. The considerable emission capability of the three molecules under study enabled the application of the two-photon excited fluorescence (TPEF) technique to probe the third order NLO properties and, thus, measure the two-photon cross-sections (σ TPA ). The twophoton excitation (TPE) spectra ( Figure 5) were obtained by exciting the sample (DMA-QF and DMA-QT in DMF) with two simultaneous photons of halved frequency, tuning the excitation radiation in the 830-1200 nm range, and acquiring the fluorescence excitation spectrum to be compared with the one-photon excitation (OPE) spectral profile. Given its hypsochromic absorption spectrum and the accessible spectral window in the employed experimental setup, the same experiment could not be carried out in the case of DMA-QP. The obtained σ TPA values were on the order of tens of GM. The TPE spectrum, in red, was found to differ from the black OPE profile in both investigated compounds (see Figure 5). However, this result is not surprising as distinct selection rules govern OPA and TPA electronic transitions, modifying transition probabilities, which are also sensitive to molecular symmetry. For this reason, the blue edge of the TPE spectra may envisage the onset of one-photon partially forbidden but two-photon allowed S 0 → S 2 transitions predicted at OPA = 317 and 323 nm (see Tables S10 and S15). The mismatch between the TPE and OPE maxima may also imply that different conformations might be preferentially excited with the TPEF technique, resulting in a slightly red-shifted fluorescence-excitation profile. As a matter of fact, the TPA process can be influenced by conformation changes as it correlates to the direction and magnitude of the transition dipole moment, which may differ when dealing with distinct conformations [52,53].
predicted at OPA = 317 and 323 nm (see Tables S10 and S15). The mismatch between the TPE and OPE maxima may also imply that different conformations might be preferentially excited with the TPEF technique, resulting in a slightly red-shifted fluorescenceexcitation profile. As a matter of fact, the TPA process can be influenced by conformation changes as it correlates to the direction and magnitude of the transition dipole moment, which may differ when dealing with distinct conformations [52,53].

Excited-State Dynamics
Ultrafast transient absorption (TA) spectroscopy was applied to shed light on the excited-state deactivation dynamics of the three investigated molecules in solvents of different polarities and to draw an overall picture of the competitive decay channels involved. The experimental results were subsequently rationalized through global analysis in order to identify all the transient species, together with their lifetimes and evolution-associated spectra (EAS). The relative results for the three molecules under study are compiled in Tables 3 and S19. Generally, the time-resolved transient spectra of the three DMA derivatives (Figures 6, S6 and S7) featured two kinds of footprints in the explored spectral window (450-780 nm); first, a broad positive excited-state absorption (ESA) band, referring to the S 1 →S n transitions; second, a negative signal perfectly matching the steady-state emission spectrum, thus attributed to the stimulated emission (SE). The fs-TA spectrum of DMA-QT in Tol is reported as a case study in Figure 6. Here, the global analysis pointed out five transients (Panel C); the first four spectral profiles share many common features, a broad peak at ca. 700 nm, formed in <1 ps, and a structured stimulated emission signal at ca. 575 nm. A closer look reveals that the emission minimum is slightly red-shifted when going from the black transient to the gray profile. This evolution is the well-known solvation (inertial, black, and diffusional, light green) dynamics [54,55], followed by further stabilization due to the structural relaxation (SR) of the molecule in the new excited configuration. The fourth transient, resembling the previous ones, is the more stabilized, emissive S 1 state, overlapping the steady-state fluorescence (green-shaded area in Figure 6, left). The violet transient, centered at 610 nm and appearing at later delays, was instead ascribed to the long-lived T 1 state. This assignment is justified by the following evidence: the transient overlaps the spectral shape resulting from nanosecond transient absorption measurement (purple-shaded area in Figure 6, left, see below); from the representative kinetics ( Figure 6, inset of Panel B), the rise time of this transient (kinetic at 615 nm) is concomitant with the S 1 decay (kinetic at 710 nm); and no trace of decay was found, meaning that its lifetime far exceeds the experimental time window of 3.2 ns. The rise-decay kinetic can be therefore interpreted as the occurrence of an ISC process. Fluorescence up-conversion (FUC) experiments corroborate the results given by the complementary transient absorption. As a matter of fact, all the aforementioned transients with almost equal lifetimes were identified by the FUC measurements, except for the absence of the dark triplet state (Figures S8-S10). The global analysis of these experimental data then showed that DMA-QT displayed the longest S 1 lifetimes (τ S1 = 840 ps in Tol), around one order of magnitude longer than those obtained for the pyridine and furan-analogs (τ S1 = 89 and 93 ps, respectively), in line with the highest fluorescence quantum yield in the DMA-derivative series.  When it comes to a more polar environment, the photobehavior deduced from the TA experiment is described by the same transients listed above, which are the two solvation components and the structural relaxation (hundreds of picoseconds), with the only exception of a new species, reported in red (see Figure 6, right, and Figures S6 and S7). This species not only presented a blue-shifted ESA band relative to the Franck-Condon singlet excited state, namely the locally excited (LE) state, but also a red-shifted SE signal, resembling the steady-state fluorescence spectrum in polar media (orange-shaded area in the relative figures). Its emissive nature was also confirmed by the fluorescence up-conversion measurements (see Figure 7 for DMA-QF in DMF as a reference example and Figures S8 and S10 for the remaining compounds). The marked reddening (ca. 100 nm) of the time-resolved emission EAS spectra (Figure 7, Panel C) over time in polar DMF, showing the stabilization of the emissive state, draws attention to the enhanced ICT character of relaxed S1. The change in the singlet excited-state nature, S1,LE → S1,ICT, for all the investigated compounds of the series occurred during the first ultrafast transient species (<2 ps), When it comes to a more polar environment, the photobehavior deduced from the TA experiment is described by the same transients listed above, which are the two solvation components and the structural relaxation (hundreds of picoseconds), with the only exception of a new species, reported in red (see Figure 6, right, and Figures S6 and S7). This species not only presented a blue-shifted ESA band relative to the Franck-Condon singlet excited state, namely the locally excited (LE) state, but also a red-shifted SE signal, resembling the steady-state fluorescence spectrum in polar media (orange-shaded area in the relative figures). Its emissive nature was also confirmed by the fluorescence upconversion measurements (see Figure 7 for DMA-QF in DMF as a reference example and Figures S8 and S10 for the remaining compounds). The marked reddening (ca. 100 nm) of the time-resolved emission EAS spectra (Figure 7, Panel C) over time in polar DMF, showing the stabilization of the emissive state, draws attention to the enhanced ICT character of relaxed S 1 . The change in the singlet excited-state nature, S 1,LE → S 1,ICT , for all the investigated compounds of the series occurred during the first ultrafast transient species (<2 ps), that is during the inertial solvation. Consequently, the newly formed and more polar singlet state was progressively stabilized until the relaxed structure was reached, resulting in the low-lying fluorescence emission at longer wavelengths. Analogous to what was observed in the steady-state fluorimetric measurements, the largest red-shift of the S 1,ICT was found for DMA-QF, while it was only hinted at in the case of the other two compounds. As observed for other symmetrical push-pull compounds [26,41,56], the singlet lifetimes in the cases of DMA-QP and DMA-QF increased one order of magnitude when going from a non-polar solvent such as Tol towards a polar medium (τ S1 = 970 and 830 ps in DMF, respectively), implying the involvement of a planar ICT state and matching the parallel increase of the Φ F . This result goes along with the decrease of triplet production in a polar solvent (see below), where the stabilization of the S 1,ICT is detrimental to the population of the triplet state. The lifetimes of the emissive excited states obtained from the FUC measurements were employed to determine the fluorescence kinetic constants (k F ) in both non-polar and polar media (Table S2). The k F values were found to be~10 8 s −1 in both Tol, as also inferred by the large oscillator strengths calculated for the S 1 →S 0 transition (Tables S6-S18) and DMF, highlighting a fully-allowed process even when the excited state turns into an ICT state. However, a small drop in k F was observed for both DMA-QP and DMA-QF due to the gain in internal conversion originating from the stabilized ICT state. in a polar solvent (see below), where the stabilization of the S1,ICT is detrimental to the population of the triplet state. The lifetimes of the emissive excited states obtained from the FUC measurements were employed to determine the fluorescence kinetic constants (kF) in both non-polar and polar media (Table S2). The kF values were found to be ~10 8 s −1 in both Tol, as also inferred by the large oscillator strengths calculated for the S1→S0 transition (Tables S6-S18) and DMF, highlighting a fully-allowed process even when the excited state turns into an ICT state. However, a small drop in kF was observed for both DMA-QP and DMA-QF due to the gain in internal conversion originating from the stabilized ICT state. Nanosecond laser-flash photolysis measurements in aerated and de-aerated Tol (Figures 8 and S12) and DMF ( Figure S11) solutions provide details on the triplet excited-state properties of the DMA derivatives under investigation, particularly the transient absorption spectra, triplet lifetimes, and quantum yields. For all the investigated compounds, the triplet formation was relatively fast (kT ~ 10 8 -10 9 s −1 , Table 4). TD-DFT calculations indeed revealed for DMA-QP an S0 → T3 transition described by a HOMO-2 → LUMO + 1 configuration almost isoenergetic with the S0 → S1 transition but showing lower charge- Nanosecond laser-flash photolysis measurements in aerated and de-aerated Tol (Figures 8 and S12) and DMF ( Figure S11) solutions provide details on the triplet excited-state properties of the DMA derivatives under investigation, particularly the transient absorption spectra, triplet lifetimes, and quantum yields. For all the investigated compounds, the triplet formation was relatively fast (k T~1 0 8 -10 9 s −1 , Table 4). TD-DFT calculations indeed revealed for DMA-QP an S 0 → T 3 transition described by a HOMO-2 → LUMO + 1 configuration almost isoenergetic with the S 0 → S 1 transition but showing lower charge-transfer character (Tables S6 and S7, and Figure S2). The presence of this triplet excited state in close proximity to S 1 may favor ISC and thus be the reason for the significant triplet kinetic constant (k T = 1.2 × 10 9 s −1 ). As for the furan and thiophene derivatives, for which the abovementioned transition was also predicted but at much higher energies (Tables S10 and S15), the noticeable values obtained for the k T process may be explained by the presence of the heavier central atoms. The ns-TA spectrum recorded in Tol is constituted by two distinct signals ( Figure 8): a positive ESA band, peaked at 460 nm for DMA-QP, 650 nm for DMA-QF, and 610 nm for DMA-QT in N 2 -purged solution, which is barely affected by solvent polarity (with the only exception of a 10 nm-blue shift in the case of the pyridine-derivative in DMF); and a hypsochromic, negative signal (<350 nm for DMA-QP, centered at 430 nm for DMA-QF and 440 nm for DMA-QT) overlapping the steady-state absorption. The latter was hence assigned to ground-state bleaching (GSB), while the positive ESA band was interpreted as a T 1 → T n signal given that its formation took place during the laser pulse, it followed a first-order decay kinetic, and molecular oxygen quenched its lifetime. This assignment was also validated by sensitization experiments from high-energy triplet donors (DTK) or energy transfer to low-energy triplet acceptors (D2TO). An additional positive signal was detected at later delays in the case of DMA-QF at ca. λ T = 420 nm in Tol and ca. λ T = 470 nm in DMF (Figures 8 and S11). Considering the rise-decay kinetics, with the rise time matching the very short decay time (τ T,N2 = 0.41 µs in Tol, see Table 4) of the ESA at 650 nm, and given its blue-shifted position relative to the GSB, the delayed signal might be tentatively assigned to the cis-trans isomer of the furan derivative, thus formed by rotation of one double bond from a highly reactive triplet state. DMA-QF faced the most marked solvent effect; the triplet production was 0.41 in Tol, but then dropped below the detection limit in DMF. The ultrashort τ T combined with a very low ∆OD limited the quantitative characterization of the triplet state in this solvent. Similarly, the modest triplet production yield measured for DMA-QP in a non-polar solvent (Φ T = 0.11 in Tol) and the short triplet lifetime, together with low fluorescence efficiency, suggest a dominant role of reactivity even in the excited-state deactivation dynamics of the pyridine substrate. In fact, when the non-radiative kinetic (k nr = 1−Φ F −Φ T τ F ) constants are considered (Table S2), high values of the order of~10 9 -10 10 s −1 are found for DMA-QF and DMA-QP in both non-polar and polar media, pointing to a major contribution of non-radiative decay channels (i.e., isomerization and internal conversion) for their excited singlet states. Table 4. Triplet properties of investigated compounds in solvents of different polarity: absorption maximum (λ T ) and triplet lifetime in the presence (τ T,AIR ) and absence (τ T,N2 ) of molecular oxygen with the relative quenching constant (k OX ); molar absorption coefficient (ε T ) and triplet quantum yield (Φ T ), kinetic constant (given by k T = Φ T /τ S1 ); singlet oxygen quantum yield (Φ ∆ ). (ФT = 0.11 in Tol) and the short triplet lifetime, together with low fluorescence efficiency, suggest a dominant role of reactivity even in the excited-state deactivation dynamics of the pyridine substrate. In fact, when the non-radiative kinetic ( = ) constants are considered (Table S2), high values of the order of ~10 9 -10 10 s −1 are found for DMA-QF and DMA-QP in both non-polar and polar media, pointing to a major contribution of nonradiative decay channels (i.e., isomerization and internal conversion) for their excited singlet states.  On the contrary, the quadrupolar-like thiophene analog showed longer triplet lifetimes (tens of microseconds) with high triplet quantum yields both in Tol (Φ T = 0.80) and in DMF (Φ T = 0.43). In the latter solvent, the sum of fluorescence and intersystem-crossing quantum yields add up to almost the entirety, pointing to the absence of other deactivation channels for DMA-QT in a polar medium. When it comes to Tol, this value exceeded 100% (Φ F + Φ T = 1.3), thus entailing some peculiar phenomenon that was worth analyzing. The effect of concentration on the emission spectra were thus evaluated up to the conditions reached in the flash photolysis experiments, which require about one order of magnitude more highly concentrated solutions. In fact, self-absorption phenomena were revealed ( Figure S13) and the fluorescence quantum yield was found to decrease from 0.51 to 0.27 with concentration. This implies that in a concentrated solution the compound is caught in some sort of radiative trapping, where the absorbed energy is transferred radiatively to a nearby molecule contributing to the quantitative population of the triplet state. In addition, the DMA-QT excited triplet state proved capable of efficiently sensitizing singlet molecular oxygen (Φ ∆~7 0%, Figure S14), making this compound alluring for PDT application.

Discussion
The study of the trans,trans isomers of the three bis[(dimethylamino)styryl]benzene derivatives under investigation proved very informative about the role of the central heteroatoms in tuning the spectral, photophysical, and NLO properties of this molecular series. In fact, the π-rich furan and thiophene rings endow the system with an increased degree of conjugation, shifting the absorption spectra towards longer wavelengths as opposed to the absorption of DMA-QP, which instead features a π-deficient pyridine ring. Conjugation also plays a leading role in determining the entity of the hyperpolarizability coefficients, which were found to be large and comparable when dealing with DMA-QF and DMA-QT, due to the great electron-transport ability of furan and thiophene [32], and about four times lower in the case of the pyridine derivative. On the other hand, this behavior revealed a modest contribution offered by charge transfer to second-order NLO properties; fluorosolvatochromism and excited-state dynamics indeed revealed a higher ICT degree (higher ∆µ CT ) for the relaxed excited state of the DMA-QF derivative, but this did not increase its β value as compared to that of DMA-QT (β = 920 and 1030 10 −30 esu −1 cm 5 , while ∆µ CT = 27.6 and 21.1 D for DMA-QF and DMA-QT, respectively). Interestingly, the close dependence of the β value on conjugation has also been corroborated by previous results on push-pull cationic chromophores [34,51], as well as on symmetric and asymmetric styrylbenzene derivatives [25,57], analogs of those reported in this work, whose hyperpolarizabilities were always found to be on the order of 1000 × 10 −30 esu −1 cm 5 either when symmetrically substituted with two nitro groups or when bearing nitro and dimethylamino substituents at the two opposite sides of the molecule, without altering the conjugation of the central structure. These values were instead halved in the case of an asymmetric substitution with a nitro and a methoxy group (460 and 420 esu −1 cm 5 for the furan and thiophene derivatives, respectively), which indeed boosts the ICT degree (∆µ CT = 34.6 and 31 D for the furan and thiophene derivatives, respectively) but reduces conjugation, thus affecting the β value [57].
Conversely, the ICT character is fundamental for enhancing third order NLO properties, such as TPA abilities, with DMA-QF exhibiting higher TPA cross-sections relative to DMA-QT. The σ TPA were, however, modest and the TPE spectra did not overlap the OPE profiles. This finding could be explained by taking into account the conformational equilibria between compressed and elongated forms. On the one hand, the compressed and bent configuration is favored because of intramolecular hydrogen-like bonds; on the other hand, the elongated form is expected to show a higher TPA cross-section because of the larger transition dipole moment associated with a straight configuration, as already reported for other systems [53], thus contributing to the redshift of the TPE spectrum. In fact, the existence of other configurations also becomes apparent when DMA-QF is studied in a protic solvent like EtOH, where intermolecular hydrogen bonds can be introduced between the heteroatom and solvent molecules, resulting in the notable broadening of the absorption spectrum.
The presence of ICT dynamics for the three molecules was evidenced by the important fluorosolvatochromism, especially marked for DMA-QF, and the significant spectral evolution probed by fs-resolved ultrafast spectroscopy. However, the lack of excited-state lifetime quenching with solvent polarity, which was instead found to increase for DMA-QP and DMA-QF when passing from a scarcely to a highly polar medium because of the competition with intersystem crossing, points to a planar ICT state which remains significantly fluorescent. The emissive capability in a polar environment differed radically from what has been previously observed for asymmetric styrylbenzene analogs, whose highly dipolar push-pull A-π-D structure leads to the formation of an ICT state that becomes non-emissive upon increasing polarity [29]. As for the DMA-QP derivative, its non-negligible emission efficiencies even in a highly polar water-AcCN mixture have already been exploited to unravel an interesting acido(fluoro)chromic behavior, upon a first protonation on the central pyridine followed by a second protonation involving the lateral dimethylamino substituent [58]. As a rule, the fluorescence quantum yields of the three bis[(dimethylamino)styryl]benzene derivatives in polar media were always higher than those of the bis(nitrostyryl)benzene analogs (Φ F = 0.20, 0.37, and 0.53 for DMA-QP, DMA-QF, and DMA-QT, respectively, vs. Φ F = 0.004, 0.10, and 0.21 for the corresponding dinitro compounds) [26], making the former better emissive probes with a view to possible fluorescence imaging applications in a biological environment. In this sense, the triplet properties also awaken deep interest in the DMA-QT molecule. While DMA-QP and DMA-QF experienced a drop for the triplet quantum yield in polar DMF, where the stabilization of the ICT state put the triplet out of reach, the intersystem crossing in DMA-QT remained very efficient with increasing polarity, with Φ T = 0.43 and very significant production of singlet oxygen by sensitization from the triplet. These properties make the dimethylaminothiophene derivative even more attractive than its nitro-substituted counterpart, where both fluorescence and triplet production were found to reduce with polarity [26].
DMA-QT, therefore, features a very interesting photobehavior in a polar environment, where its excited state deactivates entirely by fluorescence and intersystem crossing (followed by molecular oxygen sensitization) in an almost 50/50 ratio, making the molecule promising for applications as both an imaging probe and a PDT agent.

Synthesis
The synthesis procedure of the DMA-QP compound has been already described in a previous paper [58], while those regarding the furan and thiophene derivatives are reported below.
DMA-QF (4,4 -((1E,1 E)-furan-2,5-diylbis(ethene-2,1-diyl))bis (N,N-dimethylaniline). To a stirred solution of the 2,5-dimethylfuran in acetic anhydride (3 eq), pdimethylaminobenzaldehyde (2 eq, 0.45 mmol), KOAc (1 eq), and acetic anhydride (3 eq) were added. In addition, a catalytic amount of I 2 was added and the reaction mixture was left on reflux for 24 h. The reaction mixture was poured into the ice and saturated with aqueous NaOH. The precipitate was formed, filtered, and washed with water. After removal of the solvent, the residue was worked up with water and toluene and dried over MgSO 4 . The crude reaction product was chromatographed and the pure E,E-isomer of DMA-furan was obtained in the last fractions using petroleum ether-diethylether (5%) mixture as eluent. To a stirred solution of 2,5-dimethylthiophene in tetrachloromethane, N-bromosuccinimide (0.6 eq) and azobisisobutyronitrile (0.1 eq) were added. The reaction mixture was heated on reflux and irradiated with a halogen lamp (75 W) overnight. After cooling down to RT, the reaction mixture was filtrated to remove succinimide and evaporated. After the removal of the solvent, an extraction with dichloromethane and water was carried out. The extract was dried and concentrated. The crude product was dissolved in benzene and triphenylphosphine (PPh 3 ) in benzene was added. After stirring overnight at RT the precipitate was filtered off and used in the next step after drying. To a stirred solution of obtained phosphonium salt (0.45 mmol) and p-dimethylaminobenzaldehyde

Synthesis Characterization
The 1 H and 13 C NMR spectra were recorded on a spectrometer at 600 MHz (Figures S15-S18). All NMR spectra were measured in CDCl 3 using tetramethylsilane as a reference. High-resolution mass spectra (HRMS) were obtained on a matrix-assisted laser desorption/ionization time-of-flight MALDI-TOF/TOF mass spectrometer (4800 Plus MALDI-TOF/TOF analyzer, Applied Biosystems Inc., Foster City, CA, USA) equipped with Nd:YAG laser operating at 355 nm with a firing rate of 200 Hz in the positive ion reflector mode, described in detail in Ref. [58]. A total of 1600 shots per spectrum were taken with a mass range of 100-1000 Da, a focus mass of 500 Da, and a delay time of 100 ns. Nicotinamide and azithromycin were used for external mass calibration in positive ion mode.
Singlet oxygen quantum yields (Φ ∆ ) were measured in air-equilibrated Tol and DMF solutions using the DMA derivatives as sensitizers. The 1 O 2 phosphorescence spectra were detected using a spectrofluorimeter FS5 (Edinburgh Instrument) equipped with an InGaAs detector. Phenalenone (Φ ∆ = 0.99 and 1.00 in Tol and DMF, respectively) was employed as a reference compound for comparison purposes [61].
Even though approximated, the solvatochromic method managed to give a valid estimation of the second order hyperpolarizability coefficients of β CT values measured by the more sophisticated EFISH (electric-field-induced harmonic generation) or HRS (hyper-Rayleigh scattering) techniques, with the advantages of simplicity and the application of conventional steady-state stationary spectroscopy only. As for the case of the EFISH method, the solvatochromic method provides the β CT dominant contribution, which corresponds to the β XXX component of the β tensor describing the CT transition and being referred to as the frequency of the exciting laser of EFISH allows for direct comparison and good agreement between the results of the two approaches.
This method is based on the results of the fluorosolvatochromic behavior. The dependence of the Stokes shift (∆ ν) from the solvent properties (dielectric constants, ε, and refractive index, n) expressed as f (ε, n) = ε−1 ε+2 − n 2 −1 n 2 +2 retains the information about the dipolemoment difference between the ICT excited state and the ground state (∆µ CT = µ E − µ G ) according to Equation (1): where ∆ ν = ν abs − ν em is the Stokes shift (in cm −1 ), a is the cavity radius within Onsager's model, taken as 60% of the calculated diameter along the CT direction resulting from the optimized geometry [62] in cm, h is the Planck constant and c is the speed of light in a vacuum.
In the case of DMA-QT, as also reported for other quadrupolar push-pull systems [12,25], another equation was applied that relates the stokes shift to the quadrupole moment difference (∆Q = Q E − Q G ) as follows: where Q = 2µd and thus ∆Q CT = 4d 2 (∆µ CT ) 2 . In this case, the quadrupole is considered as two opposite dipoles, separated by a distance d between their barycenter, sharing the central heteroatomic ring. The dynamic hyperpolarizability coefficient was then derived through Oudar's formula [13] where r eg is the transition dipole moment related to the oscillator strength (f ) obtained by the integrated absorption band as f = 4.32 × 10 −9 ε(v)dv [63] or predicted by quantum mechanical calculations, ν eg is the maximum frequency (in cm −1 ) of the bathochromic CT transition, and ν L is the frequency of the incident radiation, chosen as 1064 nm of Nd:YAG laser, to which the β values would be referred. The static, frequency-independent hyperpolarizability coefficient β 0 has been calculated considering Equation (4) [64]: The experimental setup for the femtosecond transient absorption and fluorescence up-conversion measurements have been widely described elsewhere [3,24,27]. Briefly, the 400-nm excitation pulses of about 60 fs were generated by an amplified Ti-Sapphire laser system (Spectra Physics, Mountain View, CA, USA). The transient absorption spectrometer (Helios, Ultrafast Systems, Sarasota, FL, USA) was characterized by 150 fs time and 1.5 nm spectral resolutions. Probe pulses were produced in the 450-750 nm range by passing a small fraction of the 800 nm excitation radiation through an optical delay line (time window of 3200 ps) and focusing it onto a Sapphire crystal (2 mm thick) to generate a white-light continuum. In the fluorescence up-conversion setup (Halcyone, Ultrafast System), the 400-nm pulse excites the sample while the fundamental laser beam acts as the "gate" light. After passing through the delay line, the "gate" beam sums to the sample emission promoting the up-conversion process. The time resolution is about 200 fs while the spectra resolution is 1.5 nm. All the measurements were carried out at ca. 2.0 × 10 −3 J/cm 2 of energy fluence and under the magic angle condition, stirring the solution in a 2 mm cuvette (0.5 < A < 1.0 at λ exc = 400 nm) during the experiments to avoid photoproduct interferences. Photodegradation was checked by recording the absorption spectra before and after the time-resolved measurement, where no significant change was observed. The experimental data matrixes were first analyzed using the Surface Xplorer PRO (Ultrafast Systems) software, where it was possible to perform SVD of the 3D matrix to derive the principal components (spectra and kinetics). Successively, a global analysis using GloTarAn software was performed in order to obtain the lifetimes and the evolution-associated spectra (EAS) of the detected transient.
The same Nd:YAG pump coupled with an optical parametric oscillator (OPO-Surelite P/N 996-0210, Continuum, Newton, MA, USA) was used to perform the two-photon excited fluorescence experiments. The OPO can be tuned to produce radiation in the 410-2200 nm range (signal between 410 and 750 nm; idler between 820 and 2200 nm). The fluorescence intensity is detected by a PTM tube (Hamamatsu R2257-Y001), powered by a high-voltage power supply (PS-310, SRS) and given in mV by an oscilloscope LeCroy-Wave Runner-LT322 (500 MHz, 200 MS s −1 , DSO). Solutions of known concentration (about 1 × 10 −5 M) in DMF were prepared together with fluorescein in buffered water at pH 11 (σ = 26 GM at 930 nm) [69] as reference compounds to determine the TPA cross-sections through a comparative approach.

Quantum Mechanical Calculations
Quantum mechanical calculations were performed by using the Gaussian 16 package [70]. DFT with the B3LYP functional was chosen as the method to optimize the ground state geometry of these small organic push-pull systems and to derive their properties [71]. In contrast the lowest singlet excited states were investigated using TD-DFT excited-state calculations with WB97XD functional. Every calculation was submitted using 6-31 g + G(d) as the basis set including the solvent effect (Tol) according to the conductor-like polarizable continuum model (CPCM) [72].