Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Exceptional three- to six-photon absorption at organometallic dendrimers

Ling Zhang a, Mahbod Morshedi a, Torsten Schwich a, Rika Kobayashi b and Mark G. Humphrey *a
aResearch School of Chemistry, Australian National University, Canberra, ACT 2601, Australia. E-mail: Mark.Humphrey@anu.edu.au
bNational Computational Infrastructure, Australian National University, Canberra, ACT 2601, Australia

Received 17th February 2024 , Accepted 13th May 2024

First published on 14th May 2024


Abstract

The light-intensity dependence of multi-photon absorption (MPA) affords outstanding spatial control. Furthermore, compared to the higher-energy photons needed for analogous linear absorption, the lower-energy photons involved in MPA often correspond to important wavelengths, such as those of the biological and telecommunications “windows”. It is therefore of crucial importance to develop molecules that exhibit outstanding MPA cross-sections. However, although progress has been made with two-photon absorption, there is currently a dearth of efficient instantaneous n-photon absorbers (n > 2), a key reason being the scarcity of structure–property studies required to understand higher-order MPA. We herein report systematically-varied metallodendrimers up to third-generation in size, together with their nonlinear absorptive responses over the spectral range 600–2520 nm. We show that the dendrimers exhibit exceptional instantaneous three- to six-photon absorption cross-sections, with maximal values increasing with dendrimer generation and installation of solubilizing group, and we report that changing the groups at the dendrimer periphery can shift the wavelengths of the nPA maxima. We also describe time-dependent DFT studies that have facilitated assignment of the key linear and nonlinear transitions and disclosed the crucial role of the metal in the outstanding MPA performance.


1 Introduction

Instantaneous multi-photon absorption (MPA) is a nonlinear optical (NLO) process involving the near simultaneous absorption of two or more photons. MPA has attracted ever-increasing interest from both the applied and fundamental perspectives.1 The light intensity dependence of MPA can afford exquisite three-dimensional control of interaction volume, because any subsequent “action” can be localized to the focal point of a laser, at which the light is sufficiently intense so as to manifest MPA effects. This 3D control has demonstrated or proposed uses in micromachining, data storage, photodynamic therapy, biological imaging, theranostics, and many other applications. Another practical advantage stems from the fact that the MPA excitation process involves n-photons of frequency (ν/n) for an overall transition energy . The use of such longer wavelength photons can result in decreased photo-damage, while simultaneously opening up technologically important spectroscopic regions (e.g. the windows where tissue and other biological material transparency are maximized: near-infrared I (NIR-I), 650–950 nm; NIR-II, 1000–1350 nm; NIR-III, 1550–1870 nm; and a newly proposed window, 2080–2340 nm).2,3 Essential to the development of MPA materials are wavelength-tunable laser systems to interrogate these spectral regions, with a particular need for ultra-short (femtosecond) low-repetition rate pulses to ensure instantaneous rather than stepwise effects are assayed.

The key to exploiting MPA is the development of materials with exceptional MPA merit. At the molecular level, MPA efficiency is usually quantified in terms of the n-photon absorption cross-sections σn. Extensive research has uncovered structure-2PA (two-photon absorption) activity relationships that have identified key molecular design criteria to maximize σ2.4–10 Less is known of higher-order effects, although crucial advances have been made, largely with 3PA/σ3;11,12 reports of 4PA/σ4 have been much less frequent,13,14 while descriptions of 5PA/σ5 and even higher-order MPA are sparse (Table S1).15,16 To some extent, this deficiency reflects the lack of routine access to the aforementioned low repetition rate, ultra-short pulse length long-wavelength-tunable lasers that are needed to drive the development of structure–activity relationships.

Dendrimers have attracted attention as putative NLO materials.17–19 Their monodisperse nature permits control of optical properties, and their hyperbranched structures facilitate maximization of both effective-chromophore density and molecular solubility. Organic dendrimers have commanded most attention, but organometallic dendrimers that exploit the flexibility in metal valence electron count and coordination geometry are also of interest. Square-planar platinum-containing dendrimers, for example, have been explored as optical limiters.20 Progression from 16 valence electron (VE) platinum to 18 VE metal centers is a logical approach to improve NLO responses because these properties are dependent on electron-richness as well as ease of polarization. With this in mind, we have previously carried out wavelength-dependence NLO studies of zero- and first-generation ruthenium alkynyl dendrimers,21–26 and reported their 2PA, 3PA and, in one case,26 4PA activity. These interesting results have prompted us to examine the effect of varying dendrimer content in a more rigorous fashion, to assess the impact of dendrimer generation and composition on MPA. We herein report the syntheses and characterization of systematically-varied zeroth- to third-generation ruthenium alkynyl dendrimers with varying dendrimer framework linker length, peripheral functionalization, and solubilizing group incorporation, together with wide-wavelength femtosecond Z-scan studies that have disclosed their exceptional instantaneous 2PA–6PA performance, identified record MPA coefficients, and aided in the development of structure–MPA property relationships (the MPA studies of a few of the dendrimers that lack peripheral substituents have been reported in a preliminary fashion).27 We also report computational studies of linear optical properties and 2PA that permit assignment of the major linear and nonlinear optical transitions, and thereby help to rationalize the experimental observations.

2 Results and discussion

2.1 Synthesis and characterization of ruthenium alkynyl dendrimers

The systematically varied dendrimers in Fig. 1 permit assessment of the influence of dendrimer generation and composition on MPA. These ruthenium alkynyl dendrimers are large π-electron-delocalized macromolecules for which the constituent wedges and dendrons have rigid conjugated sub-units, and for which certain design considerations needed to be taken into account prior to synthesis. The 1,4-phenyleneethynylene (PE)-based linkages must be a minimum length to avoid insuperable steric congestion at the periphery of the second- and third-generation dendrimers. The design of these PE linkages also needs to consider the “effective conjugation length”,28 which is likely to be a key factor in maximizing NLO response scaled by molecular size. While the branching inherent in the dendritic construction aids solubility in comparison to that of similar-size linear analogues,26 it may still be insufficient to ensure the solubility of higher-generation dendrimers, so linear alkyl substituents were incorporated in the larger examples of the present study to enhance solubility (linear alkyl groups were chosen in preference to alternatives such as alkoxy groups, to avoid perturbing the OPE π-systems that contribute to the NLO effects). With these considerations in mind, the resultant suite of dendrimers depicted in Fig. 1 permits assessment of the impact on MPA merit of varying generation (zeroth- to third-generation), peripheral substituent (chlorido, phenylethynyl, 4-nitrophenylethynyl), core linker (1PE or 2PE to ruthenium), and OPE linker (0PE to 3PE between ruthenium centers and dendrimer branching points), as well as the effect of introducing solubilizing substituents (ethyl groups).
image file: d4sc01127a-f1.tif
Fig. 1 Organoruthenium dendrimers in this study. [Ru]: trans-[Ru(κ2-dppe)2] (dppe = 1,2-bis(diphenylphosphino)ethane). Wavy lines: equivalent dendrons at dendrimer cores. Names follow the format: (i) nG (n-th generation dendrimer), (ii) xy (x and y phenyleneethynylene units between the branching points and the ligated Ru at each generation level, commencing at the core), (iii) a comma “,” separating the “xy” numbers of OPE units at each generation level, and in some cases (iv) –s (installation of solubilizing Et groups at the central phenylene of the OPE unit between the first- and second-generation branching points).

The synthetic procedures are described in the ESI and are depicted in Schemes S1–S7. The dendrimers were constructed by sequences of well-established procedures: vinylidene for chlorido ligand substitution, deprotonation of vinylidene ligands to form alkynyl ligands, protodesilylation of trialkylsilyl-protected alkynes to afford terminal alkynes, and Sonogashira palladium-catalyzed C–C coupling. Overall, the chlorido-terminated dendrimers were synthesized by divergent routes, while the 4-nitrophenyl-/phenyl-alkynyl-terminated dendrimers were accessed by a combination of divergent and convergent means. In all cases, the dendrimer syntheses exploited the steric control of the extent of metalation of a 1,3,5-trisubstituted arene that is possible with the sterically demanding ligated metal center trans-[Ru(κ2-dppe)2].29 This di- rather than tri-substitution, even in the presence of an excess of ruthenium reagent, affords rapid access to the dendrimer generation-defining 1,3,5-C6H3X2Y branching points, and thereby the key dendritic “wedge” intermediates, without the need for the protection/deprotection protocols required in analogous purely organic dendrimer synthesis.

The new complexes were initially characterized by a combination of 1H, 13C, and 31P NMR and IR spectroscopies, and satisfactory elemental analyses (Fig. S1–S45). The 31P NMR spectra, in particular, proved useful in confirming reaction completion and product composition because resonances arising from the trans-[Ru(C[triple bond, length as m-dash]CR)Cl(κ2-dppe)2] (49–50 ppm) and trans-[Ru(C[triple bond, length as m-dash]CR)(C[triple bond, length as m-dash]CR′)(κ2-dppe)2] (53–55 ppm) environments are well-separated. The dendrimers were further characterized by diffusion-ordered spectroscopy (DOSY: Fig. S46–S48), size-exclusion chromatography (SEC), transmission electron microscopy (TEM), and MS. DOSY confirmed the product purity; the Stejskal–Tanner plot-derived diffusion coefficients revealed that the high-generation dendrimers, as expected, diffuse more slowly than the low-generation dendrimers (Fig. S49) and, more importantly, that dendrimer diffusion is clearly distinguishable from that of possible impurities/by-products (such as branches, dendrons, wedges and other fragments). The Stokes–Einstein equation-derived hydrodynamic radii increase with increasing generation and installation of solubilizing ethyl groups, and upon proceeding from chlorido to phenylethynyl and then 4-nitrophenylethynyl peripheral ligands (Fig. S50). The hydrodynamic radii of the dendrimers are smaller than the corresponding radii of gyration calculated from molecular mechanics-optimized geometries (Fig. S51), with the disparity increasing with increasing generation, presumably due to the increasing proclivity of the dendrimer to rotate out of planarity. The SEC traces for all dendrimers exhibited a single peak, a narrow distribution for the number-averaged molecular weights (Mn), and dispersity values consistent with the presence of uniform ruthenium dendrimers for which Mn values coincide with their formula weights (Fig. S52). TEM micrographs of a representative example (2G22,03,01) revealed individual molecules with diameters of ca. 7.5 nm,27 consistent with that calculated by molecular modelling (Fig. S53). MS proved to be of limited utility due to the low ionization efficiency of the dendrimers, the highest molecular weight example for which it proved useful being the chlorido-functionalized 1G22,00–Cl. This complex afforded MS signals that were simulated as corresponding to [M − 6Cl + nMeCN]6+ (n = 4–6),27 behavior that is consistent with the experimentally-observed facile ligand substitution at trans-[Ru(C[triple bond, length as m-dash]CR)Cl(κ2-dppe)2].30,31

2.2 Experimental and computational studies of linear optical properties

The UV-vis-NIR spectra of the dendrimers are shown in Fig. S54 and S55, with important data for the dendrimers and their precursors being listed in Tables 1 and S2, respectively. Proceeding from chlorido- to phenylethynyl-terminated dendrimer results in little change in the lowest-energy bands, but an increase in absorptivity at wavelengths corresponding to higher-energy absorptions, while installation of nitro acceptor groups at the peripheral phenylethynyl ligands results in the appearance of a lower-energy band. Calculations on related monometallic model complexes suggest that the lowest energy band in the spectra of the chlorido- and phenylalkynyl-containing dendrimers corresponds to charge transfer from Ru to the OPE-containing ligand,32 while that in the spectra of the nitro-containing examples corresponds to ligated ruthenium to nitrophenyl charge transfer.28,33 Dendrimer generation increase while otherwise maintaining dendrimer composition results in a general increase in absorptivity scaling roughly with molecular size, but no change in profile.
Table 1 Linear optical absorption and two- to four-photon absorption cross-section maximaa
Complex λ max [ε]c σ 2 (λmax)b σ 2 /M (λmax)b σ 3 (λmax)b σ 3 /M (λmax)b σ 4 (λmax)b σ 4 /M (λmax)b Ref.
a Solvent CH2Cl2. b nm. c 104 L mol−1 cm−1. d GM = 10−50 cm4 s−1 photon−1. e GM mol g−1. f 10−80 cm6 s−2 photon−2. g 10−80 cm6 s−2 photon−2 mol−1 g−1. h 10−110 cm8 s−3 photon−3. i 10−110 cm8 s−3 photon−3 mol−1 g−1.
3G22,03,02,01–NO2 427 [82] 163[thin space (1/6-em)]600 (750) 3.11 (750) 38[thin space (1/6-em)]400 (1200) 0.73 (1200) 3800 (1600) 0.072 (1600) This work
123[thin space (1/6-em)]000 (850) 2.34 (850) 3000 (1900) 0.057 (1900)
3G22,03,02,01 418 [106] 113[thin space (1/6-em)]200 (725) 2.20 (725) 22[thin space (1/6-em)]200 (1200) 0.43 (1200) 5850 (1650) 0.11 (1650) 27
45[thin space (1/6-em)]900 (900) 0.89 (900)
3G22,03,02,01–s 424 [81] 111[thin space (1/6-em)]500 (700) 2.15 (700) 23[thin space (1/6-em)]200 (1250) 0.45 (1250) 7300 (1650) 0.14 (1650) 27
50[thin space (1/6-em)]600 (900) 0.98 (900)
3G22,03,02,00–Cl 412 [110] 113[thin space (1/6-em)]700 (750) 2.26 (750) 23[thin space (1/6-em)]100 (1200) 0.46 (1200) 4250 (1800) 0.085 (1800) This work
13[thin space (1/6-em)]600 (1400) 0.27 (1400)
2G22,03,01–NO2 430 [62] 42[thin space (1/6-em)]200 (750) 1.68 (750) 15[thin space (1/6-em)]600 (1200) 0.62 (1200) 2100 (1850) 0.083 (1850) This work
7600 (1450) 0.30 (1450)
2G22,03,01 415 [49] 43[thin space (1/6-em)]600 (700) 1.78 (700) 13[thin space (1/6-em)]750 (1250) 0.56 (1250) 2950 (1650) 0.12 (1650) 27
20[thin space (1/6-em)]000 (900) 0.81 (900)
2G22,03,01–s 421 [52] 42[thin space (1/6-em)]450 (750) 1.71 (750) 15[thin space (1/6-em)]000 (1250) 0.60 (1250) 3700 (1650) 0.15 (1650) 27
16[thin space (1/6-em)]000 (875) 0.64 (875)
2G22,03,00–Cl 421 [45] 38[thin space (1/6-em)]500 (700) 1.60 (700) 13[thin space (1/6-em)]600 (1250) 0.56 (1250) 1000 (1750) 0.041 (1750) This work
1G22,01–NO2 430 [33] 22[thin space (1/6-em)]000 (680) 2.05 (680) 11[thin space (1/6-em)]400 (1200) 1.06 (1200) 900 (1750) 0.084 (1750) This work
650 (1900) 0.060 (1900)
1G22,01 427 [19] 13[thin space (1/6-em)]700 (725) 1.31 (725) 4800 (1200) 0.46 (1200) 1200 (1650) 0.11 (1650) 27
8800 (900) 0.84 (900)
1G22,00–Cl 421 [20] 12[thin space (1/6-em)]700 (680) 1.26 (680) 5800 (1200) 0.58 (1200) 300 (1850) 0.030 (1850) This work
2900 (1350) 0.29 (1350)
0G21–NO2 418 [11.6] 6000 (875) 1.55 (875) 5000 (1200) 1.29 (1200) 180 (1750) 0.046 (1750) 45
0G21 422 [12.6] 3200 (900) 0.85 (900) 2300 (1200) 0.61 (1200) 45
0G20–Cl 424 [13.5] 9500 (650) 2.68 (650) 1800 (1200) 0.51 (1200) 45
6000 (875) 1.69 (875)
2G12,02,01–NO2 424 [63] 31[thin space (1/6-em)]200 (850) 1.29 (850) 14[thin space (1/6-em)]400 (1250) 0.60 (1250) 1300 (1850) 0.054 (1850) This work
600 (2000) 0.025 (2000)
2G12,02,01 340 [96] 34[thin space (1/6-em)]700 (750) 1.47 (750) 7950 (1250) 0.34 (1250) 2700 (1650) 0.11 (1650) 27
409 [46] 18[thin space (1/6-em)]000 (900) 0.76 (900)
2G12,02,00–Cl 409 [46] 35[thin space (1/6-em)]600 (700) 1.56 (700) 13[thin space (1/6-em)]700 (1200) 0.60 (1200) 2200 (1750) 0.096 (1750) This work
1G12,01–NO2 424 [32] 23[thin space (1/6-em)]200 (850) 2.22 (850) 6200 (1200) 0.59 (1200) 400 (1950) 0.038 (1950) This work
1G12,01 343 [38] 14[thin space (1/6-em)]000 (725) 1.38 (725) 2500 (1250) 0.25 (1250) 900 (1650) 0.089 (1650) 27
412 [18] 6300 (875) 0.62 (875)
1G12,00–Cl 412 [18] 14[thin space (1/6-em)]400 (750) 1.47 (750) 3200 (1200) 0.32 (1200) 650 (1600) 0.066 (1600) This work
0G11–NO2 403 [11.2] 2500 (670) 0.70 (670) 740 (1290) 0.21 (1290) 45
1100 (800) 0.31 (800)
0G11 412 [11.6] 1490 (650) 0.43 (650) 100 (1240) 0.029 (1240) 45
370 (810) 0.11 (810)
0G10–Cl 413 [9.9] 1500 (680) 0.46 (680) 190 (1290) 0.059 (1290) 45
1050 (845) 0.32 (845)


To shed more light on the UV-vis-NIR spectra, and in particular understand charge transfer within the dendritic interior, further computational studies were undertaken in the present work. Once again, due to the size of the dendrimers, it was necessary to approach assignment of the underlying transitions responsible for the linear absorption bands by undertaking calculations on monometallic models of the key components of the dendrimers, specifically, the trans-[Ru(C[triple bond, length as m-dash]CR)(C[triple bond, length as m-dash]CR′)(κ2-dppe)2] linear units (R = Ph, R′ = (1,4-C6H4C[triple bond, length as m-dash]C)0–3Ph; R = 1,4-C6H4C[triple bond, length as m-dash]CPh, R′ = (1,4-C6H4C[triple bond, length as m-dash]C)1-2Ph; R = R′ = (1,4-C6H4C[triple bond, length as m-dash]C)2Ph) that terminate at the 1,3,5-C6H3X2Y branching points in the dendritic structures. The phenyl groups of the experimental κ2-dppe ligands were replaced by H in the model complexes (Fig. 2: dHpe ligands), an approximation that previous studies on monometallic complexes had confirmed was an acceptable compromise to mitigate computational expense.34 The geometry-optimized coordinates are given in Tables S3–S9 (PBE0 (ref. 35)/6-31G(d),36 no symmetry constraints, D3BJ dispersion correction,37 PCM CH2Cl2 (ref. 38 and 39)). The low-energy absorption bands were calculated using both the hybrid functional PBE0 and the long-range corrected functional CAM-B3LYP (6G(d,p)/SDD non-metal/Ru basis sets)40–43 (Table S10). The latter was found to significantly overestimate the transition energies and, as the former was found to better approximate the experimental data, PBE0 was employed for the subsequent linear absorption computations.


image file: d4sc01127a-f2.tif
Fig. 2 Computational models of the key dendrimer building blocks 1-M-n (n = 1–4), 2-M-n (n = 2, 3), and 3-M-3.

The transitions with calculated oscillator strengths greater than 0.3 are listed in Table S11, isodensity plots of the key orbitals involved in the low-energy transitions are depicted in Fig. S56–S62, and natural transition orbitals for the low-energy transitions are shown in Fig. S63. A number of observations can be made. Contributions to these frontier orbitals are dominated by the OPE ligands and the metal, the exceptions being the smallest complexes for which there is charge-transfer involving the Ru(κ2-dHpe)2 moiety (1-M-2: LUMO+1; 1-M-3: LUMO+2). Increasing the number of phenylethynyl units in unsymmetrical complexes while otherwise maintaining complex composition (proceeding from 1-M-2 to 1-M-3 and finally 1-M-4) results in a diminishing LUMO ← HOMO contribution to, and a diminishing red-shift in, the lowest-energy band, which is predominantly C2RuC2Ph to OPE charge-transfer in character; a similar trend is seen proceeding from 1-M-2 to 2-M-3. Increasing the number of phenylethynyl units in symmetrical complexes (proceeding from 1-M-1 to 2-M-2 and then 3-M-3) results in diminishing LUMO ← HOMO contributions to, and diminishing red-shifts for, the lowest-energy bands, which correspond to symmetric charge transfer from the core C2RuC2 unit to the OPE ligands. Comparison of the positional isomers 1-M-3 and 2-M-2, and 1-M-4 and 2-M-3, reveals that the isomers with the phenylethynyl ligand (1-M-n (n = 3, 4)) exhibit the lower energy bands. The LUMOs of the 1-M-n complexes have significant contributions from all phenyl groups of the OPE ligand with the exception of 1-M-4, for which the proximal and distal rings of the OPE make comparatively small contributions (Fig. S59), and consistent with prior reports suggesting limits in π-delocalization in such complexes.28

2.3 Experimental studies of nonlinear optical properties

The complexes in the present study are insufficiently emissive to employ multiphoton-excited fluorescence to measure their nonlinear absorption, so nonlinear absorptive and nonlinear refractive data for the dendrimers and their key precursors were obtained using a combination of open-aperture and closed-aperture Z-scan experiments over the spectral range 600–2520 nm, and employing low-repetition rate ca. 130 femtosecond light pulses. The open-aperture experiments revealed very large absorptive nonlinearities that have been replotted as the corresponding multi-photon absorption cross-sections in Fig. 3, 4 and S64–S84. The multi-photon absorption plots also contain the linear optical absorption spectra plotted at twice to six-times the wavelength, to highlight the coincidence or near-coincidence of multiples of certain bands to MPA bands. In particular, the wavelengths of the maximal values of the intensity dependence-assigned MPA correspond closely to the appropriate multiples of the wavelengths of the low-energy MLCT bands (assigned by analogy with the low-energy bands of the model complexes which are predominantly LUMO ← HOMO and MLCT in nature), rather than multiples of the higher-energy, more intense, MLCT/ILCT/LLCT admixtures. This correlation with the low-energy MLCT bands emphasizes the importance of ligated metal modules and of certain types of charge transfer in the MPA performance. The closed-aperture experiments show negative or zero γreal values for all compounds, with the negative maximal values occurring in regions of significant MPA; this is consistent with the anticipated dependence of γreal on all nonlinear absorption processes through a nonlinear Kramers–Krönig relationship.44
image file: d4sc01127a-f3.tif
Fig. 3 Wavelength dependence of the nonlinear absorption of 3G22,03,02,01–NO2. Plots of σ2 (blue scatter), σ3 (red scatter), and σ4 (green scatter), overlaid on the UV/Vis spectrum (black line), and including plots of the UV/Vis spectrum as a function of twice (blue line), three times (red line), and four times (green line) the wavelength.

image file: d4sc01127a-f4.tif
Fig. 4 Wavelength dependence of the nonlinear absorption of 3G22,03,02,01. Plots of σ2 (blue scatter), σ3 (red scatter), σ4 (green scatter), σ5 (purple scatter), and σ6 (orange scatter) overlaid on the UV/Vis spectrum (black line), and including plots of the UV/Vis spectrum as a function of twice (blue line), three times (red line), four times (green line), five times (purple line), and six times (orange line) the wavelength.

The nonlinear absorption maximal values of the dendrimers are collected in Table 1, organized by core OPE linker (Gx, x = 1, 2), dendrimer generation (nG, n = 0–3), and peripheral substituent (–NO2 = 4-nitrophenylethynyl, “blank” = phenylethynyl, –Cl = chlorido). The phenylethynyl-terminated dendrimers exhibit nonlinear absorption maxima at 700–750 nm, 875–900 nm, 1200–1250 nm, and 1650 nm, confirmed to be 2PA, 2PA, 3PA, and 4PA in nature, respectively, from the intensity dependencies of the corresponding open-aperture fs Z-scan traces (3G22,03,02,01–s: Fig. S85–S89), and supported by the location of the nPA maxima at wavelengths close to the corresponding multiples of that of the linear absorption maxima. The chlorido-terminated dendrimers exhibit nonlinear absorption maxima at similar wavelength ranges of ca. 650–750 nm (2PA), 1200–1300 nm (3PA), and 1600–1850 nm (4PA), while the 4-nitrophenyl-terminated dendrimers show nonlinear absorption maxima at the red-shifted wavelengths of ca. 650–875 nm (2PA), 1200–1450 nm (3PA), and 1600–2000 nm (4PA); the red-shift following installation of peripheral nitro substituent mimics the corresponding multiple of the 1PA profile with its longer wavelength MLCT to the nitro group, and is of a sufficient extent that 3G22,03,02,01–NO2 exhibits appreciable 4PA at wavelengths beyond 2000 nm (Fig. 3). In general, the σn maximal values increase with generation increase and installation of nitro group; the increased activity accompanying generation increase results in measurable 4PA being observed for all nG (n > 0) examples. Focusing on the σ4 data, for which few precedent structure–activity guidelines exist, the 4PA cross-sections increase with dendrimer generation, installation of alkyl solubilizing groups, and progression from peripheral chlorido and nitrophenylalkynyl to phenylalkynyl ligand.

The maximal nonlinear absorption values of the precursors 16–18, 20, 24, and 35 were also determined and are collected in Table S2, together with those of the key intermediates 22, 26, and 36. The data reveal significant 2PA and 3PA activity for all precursors and intermediates, and with the largest examples 35 and 36 also displaying measurable 4PA, outcomes consistent with the behavior of the dendrimers.

Various procedures to compare 2PA data for different molecules have been examined, including scaling the experimental data by molecular volume, cost of production,46 number of constituent chromophore units,47 and number of “effective” electrons,48,49 but the most commonly accepted approach to compare data is to scale by molecular weight, a procedure we have implemented for σn (n = 2–4) in Table 1. For the chlorido-terminated dendrimers, the σ3/M data reveal a super-linear increase on generation increase for examples with a 1,4-C2C6H4C2 core-[Ru] linkage, but the σ3/M data are relatively invariant for examples with a longer 1,4-C2C6H4C2–1,4-C6H4C2 core-[Ru] linkage (for example 0G10–Cl 0.059 (1290), 1G12,00–Cl 0.32 (1200), 2G12,02,00–Cl 0.60 (1200), c.f.0G20–Cl 0.51 (1200), 1G22,00–Cl 0.58 (1200), 2G22,03,00–Cl 0.56 (1250)). In general, the σ3/M data increase on introduction of nitro substituent while the σ4/M data are greater in the absence of the nitro group on the peripheral phenylalkynyl ligand.

The molecular weight-scaled performances of certain dendrimer synthesis intermediates also deserve comment (Table S2). The hexaruthenium wedge precursor 26, with a strongly dipolar composition, exhibits the largest σ2/M value (4.44 GM mol−1 g−1 at 900 nm). The Sonogashira coupling precursors 35 and 36 are, effectively, organometallic dendrimers that are peripherally hexaiodo- and dodecaiodo-functionalized, respectively, and thereby benefit from a strong heavy-atom effect, with striking σ2/M parameters, the largest σ3/M parameter from the present study (36), significant σ4 activity, and the second largest σ4/M parameter (35).

As mentioned above, the maximal values increase with dendrimer generation and, indeed, the third-generation dendrimers and some of the second-generation examples are sufficiently active as to also show measurable higher-order MPA at longer wavelengths, the maximal data being collected in Table 2 and the wavelength dependencies being shown in Fig. 4, S65, S67 and S69, with the assignments again being supported by intensity dependencies of the corresponding open-aperture fs Z-scan traces (3G22,03,02,01–s: Fig. S85–S89). Note that the small red-shift in σn maxima seen on replacing phenylalkynyl by chlorido peripheral ligand and the larger red-shift that results from introduction of nitro substituent results in possible 6PA activity (chlorido) and 5PA and 6PA activity (nitrophenylalkynyl) being outside our measurement window. Reports of molecular 5PA are rare; excluding our dendrimer study, only two σ5 data are extant [(E)-3-(4-(2-(1-hexyl-4-methyl-1H-imidazole-5-yl)vinyl)pyridinium-1-yl)propyl sulfate] (0.00192 × 10−140 cm10 s−4 photon−4 at 2100 nm, fs nonlinear transmission)16 and a spiro-fused ladder-type oligo(p-phenylene) (9320 × 10−140 cm10 s−4 photon−4 at 1540 nm, 120 fs Z-scan and five-photon excited fluorescence).15 The present data reveal an increase in σ5 with generation increase, replacement of phenylalkynyl by chlorido ligand, and installation of solubilizing alkyl groups, and with the σ5/M data revealing that the dendrimer generation increases are super-linear in nature. Dendrimers 3G22,03,02,01–s and 3G22,03,02,01 display the first quantifiable molecular 6PA data for organometallics and only the second and third overall, the precedent being the aforementioned spiro-fused ladder-type oligo(p-phenylene) which is active at a much shorter incident wavelength (86.7 × 10−170 cm12 s−5 photon−5 at 1820 nm, 120 fs Z-scan and six-photon excited fluorescence).15 As with σ4 and σ5, the σ6 and σ6/M data increase with introduction of the ethyl solubilizing groups.

Table 2 Linear optical absorption and five- and six-photon absorption cross-section maximaa
Complex λ max [ε]c σ 5 (λmax)b σ 5 /M (λmax)b σ 6 (λmax)b σ 6 /M (λmax)b
a Solvent CH2Cl2. b nm. c 104 L mol−1 cm−1. d 10−140 cm10 s−4 photon−4. e 10−140 cm10 s−4 photon−4 mol−1 g−1. f 10−170 cm12 s−5 photon−5. g 10−170 cm12 s−5 photon−5 mol−1 g−1.
3G22,03,02,01 (ref. 27) 418 [106] 350 (2100) 0.0068 (2100) 25 (2420) 0.00049 (2420)
3G22,03,02,01–s27 424 [81] 500 (2050) 0.0096 (2050) 53 (2470) 0.0010 (2470)
3G22,03,02,00–Cl 412 [110] 500 (2300) 0.0100 (2300)
2G22,03,01 (ref. 27) 415 [49] 100 (2160) 0.0041 (2160)
2G22,03,01–s27 421 [52] 170 (2100) 0.0068 (2100)
2G22,03,00–Cl 421 [45] 110 (2160) 0.0046 (2160)


2.4 Computational studies of two-photon absorption

To shed further light on the outstanding nonlinear absorption behavior of the dendrimers, the 2PA cross-sections of the aforementioned dendritic fragment model complexes 1-M-n (n = 1–3) were calculated through a residue of the quadratic response function50 implemented in Dalton 2020.1 (ref. 51 and 52) (CAM-B3LYP, dHpe co-ligands, PCM53 CH2Cl2, stuttgart_rsc_1997_ecp54 (Ru) and 6+G(d)55 (non-metal atoms) basis sets; further details are given in the ESI). Calculation of 2PA spectra is computationally demanding, so the two-photon transition strengths were computed using the ten lowest excited states only. Despite this necessary simplification, the calculations are generally in good agreement with the experimental results; the 2PA profiles are similar, with one (1-M-1/1-M-1(dppe)), two (1-M-2/1-M-2(dppe)), and three (1-M-3/1-M-3(dppe)) maxima in the 450 nm range, and the significant red shift in the lowest-energy 2PA maximum that is seen experimentally upon lengthening the OPE alkynyl ligand being reproduced computationally (Table 3 and Fig. 5). We note that the shortest wavelength 2PA maxima coincide with the onset of one-photon absorption and deconvoluting experimental data in this wavelength range is known to be difficult,56 so comments are necessarily cautious.
Table 3 Experimentala and calculatedb 2PA maximal values
Complex σ 2,exp (GM) (λmax (nm)) Model σ 2,calc[thin space (1/6-em)](GM) (λmax (nm))
a Ph2PCH2CH2PPh2 (dppe) ligands. b H2PCH2CH2PH2 (dHpe) ligands.
1-M-1(dppe) 630 ± 80 (600) 1-M-1 290 (574)
1-M-2(dppe) 860 ± 200 (650) 1-M-2 850 (540)
1000 ± 170 (750) 470 (704)
1-M-3(dppe) 1400 ± 170 (600) 1-M-3 370 (572)
1560 ± 180 (700) 1100 (623)
900 ± 70 (850) 800 (763)



image file: d4sc01127a-f5.tif
Fig. 5 Calculated 2PA cross-sections of 1-M-n (n = 1–3) (red traces), experimental 2PA cross-sections of 1-M-n(dppe) (n = 1–3) (red circles), linear absorption spectra of 1-M-n(dppe) (n = 1–3) (black traces), and plots of the UV/Vis spectrum as a function of twice the wavelength (light gray traces). Insets: the metal-to-ligand charge transfer (MLCT) corresponding to the calculated low-energy 2PA maxima.

The 2PA peak of 1-M-1 at 574 nm is attributed to the largest transition moment (4.32 eV (287 nm): Table S12), corresponding to the linear MLCT transition at 285 nm (1-M-1 calc.) and 328 nm (1-M-1(dppe) exp.). Similarly, the lower-lying 2PA maximum of 1-M-2 and the lowest-lying 2PA maximum of 1-M-3 correspond to metal-to-alkynyl-ligand charge transfer, confirming the key role the metal plays in the MPA merit of these complexes.

The calculated 2PA cross-sections for 1-M-n (n = 1–3) are collected in Table 3, together with those of the experimental complexes 1-M-n(dppe) (n = 1–3). The increase in maximal values seen on proceeding from 1-M-1 to 1-M-3 reproduces the experimental trend seen with 1-M-n(dppe) (n = 1–3). The clear increase in σ2 (both maximal value and integrated intensity) upon OPE alkynyl ligand lengthening highlights the desirability of incorporating the longer OPE units in MPA materials, and confirms that the OPE-rich dendrimer composition in the present study is ideal for MPA efficiency.

3 Conclusions

The present study has identified optimal compositions for MPA-efficient molecules. The trans-[Ru(κ2-dppe)2] modules incorporate polarizable low oxidation state pseudo-octahedral electron-rich metal centers, for which stepwise chemistry at the two axial sites permits the designed construction of π-electron-delocalizable architectures in which the metal engages in strong low-energy MLCT interactions. The trans-[Ru(κ2-dppe)2] units impart greatly superior solubility and optical nonlinearity to those of analogous purely organic aryleneethynylene-based structures. The resultant metallodendrimers have been definitively identified by a combination of small molecule and macromolecule characterization techniques.

Multi-photon absorption studies have identified exceptional performance across a spectral range extending deep into the NIR region. The dendrimers exhibit outstanding 2PA, 3PA, and 4PA cross-sections. Installation of peripheral nitro substituents red-shifts the MPA maxima, such that exceptional 4PA is seen at wavelengths beyond 2000 nm. Structural modifications that enhance MPA merit have been identified: increasing the length of the OPE linkage to the core, installation of alkyl solubilizing groups, and peripheral phenylalkynyl ligand incorporation are all beneficial compositional changes. Examination of the synthesis intermediates has identified peripheral iodination as a further beneficial compositional change. Computational studies of linear and two-photon absorption suggest the use of the 3PE linker as optimum for π-delocalization and maximizing MPA merit.

In addition to record values of (2–4)PA, the present study has also afforded molecules exhibiting higher-order MPA. Rare examples of molecular 5PA have been found with the second- and third-generation metallodendrimers. Generation increase, installation of peripheral chlorido ligands, and incorporation of alkyl solubilizing groups have been identified as molecular modifications that improve 5PA. The second and third examples of molecular 6PA have been found with the third-generation metallodendrimers; these exhibit activity at much longer wavelengths than the only extant literature example. The σ6 and σ6/M parameters of the dendrimers increase on incorporation of alkyl substituents at the OPE group between the first- and second-generation branching points.

Data availability

The data that support the findings of this study are available in the ESI of this article.

Author contributions

Conceptualization: M. G. H. Data curation: L. Z., M. M. Formal analysis: L. Z., M. M., T. S. Funding acquisition: M. G. H. Investigation: L. Z., M. M., T. S. Methodology: L. Z., M. M., R. K., M. G. H. Project administration: M. G. H. Resources: M. G. H. Software: L. Z., M. M., R. K. Supervision: M. G. H. Validation: L. Z., M. M., R. K. Visualization: L. Z., M. M. Writing – original draft: L. Z., M. M., M. G. H. Writing – review and editing: L. Z., M. M., R. K., M. G. H.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

We thank the Australian Research Council for support (DP220100111 to M. G. H.). L. Z. thanks the China Scholarship Council and the Australian National University for an ANU-CSC Postgraduate Scholarship. This work was supported by computational resources provided by the Australian Government through the National Computational Infrastructure (NCI) under the ANU Merit Allocation Scheme (project eo34: M. G. H.).

Notes and references

  1. G. S. He, L.-S. Tan, Q. Zheng and P. N. Prasad, Chem. Rev., 2008, 108, 1245 CrossRef CAS PubMed.
  2. Z. Feng, T. Tang, T. Wu, X. Yu, Y. Zhang, M. Wang, J. Zheng, Y. Ying, S. Chen, J. Zhou, X. Fan, D. Zhang, S. Li, M. Zhang and J. Qian, Light: Sci. Appl., 2021, 10, 197 CrossRef CAS PubMed.
  3. C. J. Liu, A. Roy, A. A. Simons, D. M. Farinella and P. Kara, Sci. Rep., 2020, 10, 16351 CrossRef CAS PubMed.
  4. M. Albota, D. Beljonne, J. L. Bredas, J. Ehrlich, J. Y. Fu, A. A. Heikal, S. E. Hess, T. Kogej, M. D. Levin, S. R. Marder, D. McCord-Maughon, J. W. Perry, H. Rockel, M. Rumi, G. Subramaniam, W. W. Webb, X. L. Wu and C. Xu, Science, 1998, 281, 1653 CrossRef CAS PubMed.
  5. M. P. Joshi, J. Swiatkiewicz, F. Xu, P. N. Prasad, B. A. Reinhardt and R. Kannan, Opt. Lett., 1998, 23, 1742 CrossRef CAS PubMed.
  6. B. A. Reinhardt, L. L. Brott, S. J. Clarson, A. G. Dillard, J. C. Bhatt, R. Kannan, L. Yuan, G. S. He and P. N. Prasad, Chem. Mater., 1998, 10, 1863 CrossRef CAS.
  7. B. R. Cho, K. H. Son, S. H. Lee, Y. S. Song, Y. K. Lee, S. J. Jeon, J. H. Choi, H. Lee and M. H. Cho, J. Am. Chem. Soc., 2001, 123, 10039 CrossRef CAS PubMed.
  8. L. Xu, J. Zhang, L. Yin, X. Long, W. Zhang and Q. Zhang, J. Mater. Chem. C, 2020, 8, 6342 RSC.
  9. L. Xu, W. Lin, B. Huang, J. Zhang, X. Long, W. Zhang and Q. Zhang, J. Mater. Chem. C, 2021, 9, 1520 RSC.
  10. S. J. Weishäupl, D. C. Mayer, Y. Cui, P. Kumar, H. Oberhofer, R. A. Fischer, J. Hauer and A. Pöthig, J. Mater. Chem. C, 2022, 10, 6912 RSC.
  11. M. Chołuj, R. Behera, E. F. Petrusevich, W. Bartkowiak, Md. M. Alam and R. Zaleśny, J. Phys. Chem. A, 2022, 126, 752 CrossRef PubMed.
  12. F. Zhou and W. Ji, Laser Photonics Rev., 2017, 11, 1700021 CrossRef.
  13. F. E. Hernández, K. D. Belfield, I. Cohanoschi, M. Balu and K. J. Schafer, Appl. Opt., 2004, 43, 5394 CrossRef PubMed.
  14. J. Szeremeta, M. Nyk, D. Wawrzynczyk and M. Samoc, Nanoscale, 2013, 5, 2388 RSC.
  15. Y. Jiang, K. F. Li, K. Gao, H. Lin, H. L. Tam, Y. Y. Liu, Y. Shu, K. L. Wong, W. Y. Lai, K. W. Cheah and W. Huang, Angew. Chem., Int. Ed., 2021, 60, 10007 CrossRef CAS PubMed.
  16. Q. Zheng, H. Zhu, S.-C. Chen, C. Tang, E. Ma and X. Chen, Nat. Photonics, 2013, 7, 234 CrossRef CAS.
  17. W. Wu, C. Wang, Q. Li, C. Ye, J. Qin and Z. Li, Sci. Rep., 2014, 4, 6101 CrossRef CAS PubMed.
  18. R. Tang, S. Zhou, Z. Cheng, G. Yu, Q. Peng, H. Zeng, G. Guo, Q. Li and Z. Li, Chem. Sci., 2017, 8, 340 RSC.
  19. W. Wu, C. Ye, J. Qin and Z. Li, ACS Appl. Mater. Interfaces, 2013, 5, 7033 CrossRef CAS PubMed.
  20. R. Vestberg, R. Westlund, A. Eriksson, C. Lopes, M. Carlsson, B. Eliasson, E. Glimsdal, M. Lindgren and E. Malmström, Macromolecules, 2006, 39, 2238 CrossRef CAS.
  21. C. E. Powell, J. P. Morrall, S. A. Ward, M. P. Cifuentes, E. G. A. Notaras, M. Samoc and M. G. Humphrey, J. Am. Chem. Soc., 2004, 126, 12234 CrossRef CAS PubMed.
  22. M. P. Cifuentes, C. E. Powell, J. P. Morrall, A. M. McDonagh, N. T. Lucas, M. G. Humphrey, M. Samoc, S. Houbrechts, I. Asselberghs, K. Clays, A. Persoons and T. Isoshima, J. Am. Chem. Soc., 2006, 128, 10819 CrossRef CAS PubMed.
  23. M. Samoc, J. P. Morrall, G. T. Dalton, M. P. Cifuentes and M. G. Humphrey, Angew. Chem., Int. Ed., 2007, 46, 731 CrossRef CAS PubMed.
  24. R. L. Roberts, T. Schwich, T. C. Corkery, M. P. Cifuentes, K. A. Green, J. D. Farmer, P. J. Low, T. B. Marder, M. Samoc and M. G. Humphrey, Adv. Mater., 2009, 21, 2318 CrossRef CAS.
  25. K. A. Green, M. P. Cifuentes, M. Samoc and M. G. Humphrey, Coord. Chem. Rev., 2011, 255, 2025 CrossRef CAS.
  26. P. V. Simpson, L. A. Watson, A. Barlow, G. Wang, M. P. Cifuentes and M. G. Humphrey, Angew. Chem., Int. Ed., 2016, 55, 2387 CrossRef CAS PubMed.
  27. L. Zhang, M. Morshedi and M. G. Humphrey, Angew. Chem., Int. Ed., 2022, 61, e202116181 CrossRef CAS PubMed.
  28. B. Babgi, L. Rigamonti, M. P. Cifuentes, T. C. Corkery, M. D. Randles, T. Schwich, S. Petrie, R. Stranger, A. Teshome, I. Asselberghs, K. Clays, M. Samoc and M. G. Humphrey, J. Am. Chem. Soc., 2009, 131, 10293 CrossRef CAS PubMed.
  29. C. E. Powell, S. K. Hurst, J. P. Morrall, R. L. Roberts, M. P. Cifuentes, M. Samoc and M. G. Humphrey, Organometallics, 2007, 26, 4456 CrossRef CAS.
  30. D. Touchard, C. Morice, V. Cadierno, P. Haquette, L. Toupet and P. H. Dixneuf, J. Chem. Soc., Chem. Commun., 1994, 859 RSC.
  31. P. J. West, M. P. Cifuentes, T. Schwich, M. D. Randles, J. P. Morrall, E. Kulasekera, S. Petrie, R. Stranger and M. G. Humphrey, Inorg. Chem., 2012, 51, 10495 CrossRef CAS PubMed.
  32. C. E. Powell, M. P. Cifuentes, J. P. L. Morrall, R. Stranger, M. G. Humphrey, M. Samoc, B. Luther-Davies and G. A. Heath, J. Am. Chem. Soc., 2003, 125, 602 CrossRef CAS PubMed.
  33. B. A. Babgi, M. S. Kodikara, M. Morshedi, H. Wang, C. Quintana, T. Schwich, G. J. Moxey, N. Van Steerteghem, K. Clays, R. Stranger, M. P. Cifuentes and M. G. Humphrey, ChemPlusChem, 2018, 83, 630 CrossRef CAS PubMed.
  34. L. Zhang, M. Morshedi, M. Kodikara and M. G. Humphrey, Angew. Chem., Int. Ed., 2022, 62, e202208168 Search PubMed.
  35. C. Adamo and V. Barone, J. Chem. Phys., 1999, 110, 6158 CrossRef CAS.
  36. W. J. Hehre, R. Ditchfield and J. A. Pople, J. Chem. Phys., 1969, 56, 2257 CrossRef.
  37. S. Grimme, S. Ehrlich and L. Goerigk, J. Comput. Chem., 2011, 32, 1456 CrossRef CAS PubMed.
  38. S. Miertus, E. Scrocco and J. Tomasi, Chem. Phys., 1981, 55, 117 CrossRef CAS.
  39. J. Tomasi, B. Mennucci and R. Cammi, Chem. Rev., 2005, 105, 2999 CrossRef CAS PubMed.
  40. T. Yanai, D. Tew and N. Handy, Chem. Phys. Lett., 2004, 393, 51 CrossRef CAS.
  41. M. J. G. Peach, T. Helgaker, P. Sazek, T. W. Keal, O. B. Lutnæs, D. J. Tozer and N. C. Handy, Phys. Chem. Chem. Phys., 2006, 8, 558 RSC.
  42. R. Krishnan, J. S. Binkley, R. Seeger and J. A. Pople, J. Chem. Phys., 1980, 72, 650 CrossRef CAS.
  43. D. Andrae, U. Haeussermann, M. Dolg, H. Stoll and H. Preuss, Theor. Chem. Acc., 1990, 77, 123 Search PubMed.
  44. M. Sheik-Bahae, D. J. Hagan and E. W. van Stryland, Phys. Rev. Lett., 1990, 65, 96 CrossRef CAS PubMed.
  45. T. Schwich, A. Barlow, M. P. Cifuentes, J. Szeremeta, M. Samoc and M. G. Humphrey, Chem.–Eur. J., 2017, 23, 8395 CrossRef CAS PubMed.
  46. T. Schwich, M. P. Cifuentes, P. A. Gugger, M. Samoc and M. G. Humphrey, Adv. Mater., 2011, 23, 1433 CrossRef CAS.
  47. M. Drobizhev, A. Karotki, M. Kruk and A. Rebane, Chem. Phys. Lett., 2002, 355, 175 CrossRef CAS.
  48. J. P. Moreno and M. G. Kuzyk, J. Chem. Phys., 2005, 123, 194101 CrossRef PubMed.
  49. M. G. Kuzyk, J. Pérez-Moreno and S. Shafei, Phys. Rep., 2013, 529, 297 CrossRef.
  50. J. Olsen and P. Jørgensen, J. Chem. Phys., 1985, 82, 3235 CrossRef CAS.
  51. K. Aidas, C. Angeli, K. L. Bak, V. Bakken, R. Bast, L. Boman, O. Christiansen, R. Cimiraglia, S. Coriani, P. Dahle, E. K. Dalskov, U. Ekström, T. Enevoldsen, J. J. Eriksen, P. Ettenhuber, B. Fernández, L. Ferrighi, H. Fliegl, L. Frediani, K. Hald, A. Halkier, C. Hättig, H. Heiberg, T. Helgaker, A. C. Hennum, H. Hettema, E. Hjertenæs, S. Høst, I.-M. Høyvik, M. F. Iozzi, B. Jansik, H. J. Aa. Jensen, D. Jonsson, P. Jørgensen, J. Kauczor, S. Kirpekar, T. Kjærgaard, W. Klopper, S. Knecht, R. Kobayashi, H. Koch, J. Kongsted, A. Krapp, A. Kristensen, A. Ligabue, O. B. Lutnæs, J. I. Melo, K. V. Mikkelsen, R. H. Myhre, C. Neiss, C. B. Nielsen, P. Norman, J. Olsen, J. M. H. Olsen, A. Osted, M. J. Packer, F. Pawlowski, T. B. Pedersen, P. F. Provasi, S. Reine, Z. Rinkevicius, T. A. Ruden, K. Ruud, V. Rybkin, P. Salek, C. C. M. Samson, A. Sánchez de Merás, T. Saue, S. P. A. Sauer, B. Schimmelpfennig, K. Sneskov, A. H. Steindal, K. O. Sylvester-Hvid, P. R. Taylor, A. M. Teale, E. I. Tellgren, D. P. Tew, A. J. Thorvaldsen, L. Thøgersen, O. Vahtras, M. A. Watson, D. J. D. Wilson, M. Ziolkowski and H. Ågren, WIREs Comput. Mol. Sci., 2014, 4, 269 CrossRef CAS PubMed , Dalton, a Molecular Electronic Structure Program, Release Dalton2020.1, 2022: see http://daltonprogram.org/.
  52. H. Ågren, O. Vahtras, H. Koch, P. Jørgensen and T. Helgaker, J. Chem. Phys., 1993, 98, 6417 CrossRef.
  53. L. Frediani, H. Ågren, L. Ferrighi and K. Ruud, J. Chem. Phys., 2005, 123, 144117 CrossRef PubMed.
  54. D. Andrae, U. Haeussermann, M. Dolg, H. Stoll and H. Preuss, Theor. Chem. Acc., 1990, 77, 123 Search PubMed.
  55. R. Krishnan, J. S. Binkley, R. Seeger and J. A. Pople, J. Chem. Phys., 1980, 72, 650 CrossRef CAS.
  56. K. Kamada, K. Ohta, Y. Iwase and K. Kondo, Chem. Phys. Lett., 2003, 372, 386 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available: Methods, instrumentation, Fig. S1–S81 and Tables S1–S12. See DOI: https://doi.org/10.1039/d4sc01127a

This journal is © The Royal Society of Chemistry 2024