Understanding the Second-Order Nonlinear Optical Properties of One-Dimensional Ruthenium(II) Ammine Complexes

Abstract

First principles calculations of the linear and nonlinear optical properties are carried out for three series of one-dimensional ruthenium(II) ammine complexes with one pyridyl pyridinium ligand containing different numbers of C2H2 units or the (CHCH)Ph(CHCH) linker. The substitution effects upon introducing one pyridine or one N-methylimidazole moiety as donor instead of ammonia have also been studied in detail. These calculations employing density functional theory with different exchange-correlation functionals as well as the polarizable continuum model approach to describe the solvent effects show that these compounds are challenging for theoretical chemistry and that their nonlinear optical responses are complex and depend on many structural and electronic factors. Two major types of methods have been employed to calculate and analyze the first hyperpolarizabilities (β): (i) the summation-over-states scheme, applied to calculate the longitudinal β tensor component, truncated to one or two dominant excited states, of which the results are compared to Stark spectroscopy data, and (ii) the quadratic response schemes, providing the whole second harmonic generation response of the compound, which can be compared to the values derived from hyper-Rayleigh scattering measurements. Using response theory calculations, off-resonant first hyperpolarizability values display a monotonic increase as a function of the number of C2H2 units (from 0 to 3 units, corresponding to compounds 1–4, respectively) in the pyridyl pyridinium ligand (in the three series with NH3, pyridine, and N-methylimidazole substituent in opposite position to the large ligand). These calculations, which evidence the dominance of the longitudinal β tensor component and a consistent depolarization ratio close to 5.0, also predict a smaller β value when the (CHCH)3 linker of the pyridyl pyridinium ligand (compound 4) is replaced by the (CHCH)Ph(CHCH) linker (compound 5). To some extent, the calculations performed at a wavelength of 1064 nm confirm these trends, despite the presence of resonance enhancements. Calculations carried out within the summation-over-states scheme to analyze the origin of these β variations reproduce the monotonic increase of β in the 1–4 series, but compounds 4 and 5 are predicted to have similar β values. Upon relating these β variations to those of spectroscopic properties like the excitation energies, transition dipoles, and dipole moments of a pair of excited states (a low-energy metal to ligand charge transfer and a high-energy intraligand charge transfer excited state), the calculations point out the impact of including a second excited state in the treatment, since it leads to a decrease of β by about 50%. Though substituting the ammonia group by a pyridine has a limited impact on the linear and nonlinear optical properties of the complexes, the substitution by N-methylimidazole increases β by about 25%, which has been attributed to the contribution of N-methylimidazole to the HOMO and the subsequent batho...

Theoretical Modelling of Photoswitching of Hyperpolarisabilities in Ruthenium Complexes

“…On the other hand, the basis set LAN2TZ(f) for Ru is about the minimum necessary for a reliable description of both the GS and FDAES properties. This last factor explains the main differences between our data and those obtained by Inerbaev et al and by Zhang and Champagne, who used the smaller basis set LAN2DZ 23, 24…”

“…Selected simulated spectra are shown in Figure 2, together with the experimental spectra of complexes 3 and 4 as their PF 6 − salts. Comparisons with data published previously23, 24 show that although using a larger basis set for Ru does not affect the qualitative picture of the excitation properties, some significant quantitative differences are observed.…”

“…Previous gas-phase DFT calculations using the B3P86 functional with the LANL2DZ basis set proved qualitatively useful, but lack quantitative accuracy 21. A similar basis set combination LANL2DZ(Ru)/6-31G* or 6-31+G*(H,C,N,O) has been used by Inerbaev et al,23 and in a recent investigation by Zhang and Champagne 24. Here, we use a larger basis set (LANL2TZ(f)) for Ru in DFT calculations, as well as different larger basis sets for various high-level ab initio methods.…”

“…[g] Data taken from ref. 24. [h] For [ 2 ][ClO 4 ] 3 in H 2 O (a very weak NIR band at λ max = 855 nm is observed also),17e [ 3 ][PF 6 ] 3 ,19 and [ 4 ][PF 6 ] 3 ;21 λ max values in MeCN at room temperature; directly corresponding f os values are unavailable, so those quoted are in PrCN at 77 K (only slight variations due to changing the solvent and temperature are expected) 20…”

Spectroscopic and second-order nonlinear optical properties of Ruthenium(ii) complexes: a DFT/MRCI and ADC(2) study

“…Table 1 summarizes the results for some of the lowest excited states, i.e. 6. The gas phase TD-B3P86 results are only given for the FDAES, since the LDES is not accessible through linear-response (LR)-TD-DFT calculations.…”