Transition metals induce control of enhanced NLO
properties of functionalized organometallic complexes
under laser modulations

S Taboukhat, N Kichou, Jean-Luc Fillaut, K Waszkowska, A Zawadzka, A
Migalska-Zalas, Olivier Alévêque, Abdelkrim El-Ghayoury, Bouchta Sahraoui

► To cite this version:

S Taboukhat, N Kichou, Jean-Luc Fillaut, K Waszkowska, A Zawadzka, et al.. Transition metals induce control of enhanced NLO properties of functionalized organometallic complexes under laser modulations. Scientific Reports, 2020, 10 (1), pp.15292. 10.1038/s41598-020-71769-2. hal-02959802

HAL Id: hal-02959802

https://hal.science/hal-02959802v1

Submitted on 7 Oct 2020

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

OPEN

# Transition metals induce control of enhanced NLO properties of functionalized organometallic complexes under laser modulations

S. Taboukhat¹, N. Kichou²,³, J.-L. Fillaut², O. Alévêque¹, K. Waszkowska¹, A. Zawadzka⁴, A. El-Ghayoury¹, A. Migalska-Zalas⁵ & B. Sahraoui¹

The molecular engineering of organometallic complexes has recently attracted renewed interest on account of their potential technological applications for optoelectronics in general and optical data storage. The transition metal which induces control of enhanced nonlinear optical properties of functionalized organometallic complexes versus not only the intensity but also the polarization of the incident laser beam is original and important for all optical switching. This makes organometallic complexes valuable and suitable candidates for nonlinear optical applications. In the present work, we report the synthesis and full characterization of four organometallic complexes consisting of N, N-dibutylamine and azobenzene fragments but differ by auxiliary alkynyl ligands or metal cations. Thus, a ferrocenyl derivative 1 and three ruthenium complexes 2–4 have been prepared. The nonlinear optical properties of the four new azo-based ruthenium and iron organometallic complexes in the solid state, using polymethylmethacrylate as matrix, have been thoroughly studied. This concept is extended to computing the HOMO and LUMO energy levels of the considered complexes, dipole moment, first and second order hyperpolarizabilities using the 6–31+G(d,p) + LANL2DZ mixed basis set. The second and third nonlinear optical properties of the resulting polymer composites were obtained by measuring SHG and THG response by means of the Maker fringe technique using a laser generating at 1,064 nm with a 30 ps pulse duration. The values of the second and third order NLO susceptibilities of the four organometallic complexes were found to be higher than the common references used. Theoretical calculation shows that the large first and second order hyperpolarizabilities are caused by strong intramolecular charge transfer between the transition metal parts and the ligands though a conjugated transmitter. These results indicate that the present organometallic complexes are valuable candidates for optoelectronic and photonic applications.

The organometallic complexes are among the most studied materials in nonlinear optics (NLO)¹,². Their opto-electronic properties are highly dependent on the metal-d orbital occupation and energy levels, which can be easily modulated through the choice of the metal, its oxidation state and the surrounding ligands³. Among these complexes, the family of alkynyl-metal complexes that can be defined as carbon-rich transition metal complexes which have a nearly linear M–C≡C–R structure demonstrated particularly suitable characteristics for NLO studies due to an efficient electronic coupling between the metal and distant groups via the π-conjugate pathway. This electronic coupling, known as the metal-to-ligand charge transfer (MLCT), is highly dependent on the degree of overlap between the filled metal d orbitals and the unoccupied π-orbitals of the alkynyl moiety⁴⁻⁶. The physical and optoelectronic properties of the metal alkynyl complexes can be optimized through a careful design, related to the choice of both the organic and the metallic end groups, for example by using iron or ruthenium, and of

¹MOLTECH-Anjou, UMR 6200, CNRS, Univ Angers, 2 Bd Lavoisier, 49045 Angers, France. ²Univ Rennes, CNRS, ISCR-UMR 6226, 35000 Rennes, France. ³Department of Chemistry, Faculty of Sciences, University of Mouloud Mammeri, 15000 Tizi-ouzou, Algeria. ⁴Department of Automatic and Measurement Systems, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland. ⁵Faculty of Science and Technology, Jan Dlugosz University in Czestochowa, Al. Armii Krajowej 13/15, 42-200 Czestochowa, Poland. ⁶email: bouchta.sahraoui@univ-angers.fr
the organic bridges of the π-conjugated entity<sup>7-12</sup>. This diversity could be used for different applications such as optoelectronics, quantum electronics, and optical communications and in the field of molecular engineering of functional materials<sup>13-15</sup>. Moreover it could be also useful for tuning the properties of organometallic based magnetics for spintronic and quantum computing devices<sup>16</sup> and also electronic and optical properties for photovoltaic and optoelectronic applications<sup>17</sup>.

Much research has been carried out on the study of azobenzene derivatives which, thanks to their reversible and stable *trans*-cis photoisomerisation over many cycles, are among the most studied photochromic compounds in recent years. Azobenzene systems possess also obvious advantages in the field of NLO, due to their high intra-molecular charge transfer and their high hyperpolarizabilities<sup>18-20</sup> which can be coupled with the photoinduced *trans*-cis isomerization. Hence, they could function as NLO molecular switches by applying light of different wavelengths to obtain varying amounts of cis and trans isomers. On the other hand, donor-acceptor substituted azo dyes, which are molecules with easily polarizable electrons, show large second-order nonlinearities<sup>21,22</sup>. A major advantage is also the possibility of inducing a noncentrosymmetric structure within the material, for example by orientation by the Corona poling technique, because obtaining macroscopic quadratic NLO properties requires breaking the centrosymmetry of the medium by orientation of the chromophores. These properties of azobenzene derivatives give the possibilities of many applications in optical switching and data storage<sup>23,24</sup>.

We report herein the synthesis and full characterization of four new azo-based ruthenium and iron organometallic complexes. The nonlinear optical (NLO) responses of these complexes, related to metal acetylide complexes, have been evaluated by means of second and third harmonic generations (SHG, THG) experimental techniques and also quantum chemical calculations. We undertook this study to better understand the physical contribution of the metal center and its coordination sphere to the NLO susceptibilities of the complexes. The parameters responsible for second and third-order nonlinear-optical effects such as second and third harmonic generation are the first and second order hyperpolarizabilities. A conventional key to achieving high values of these parameters in metal acetylide complexes is to lengthen the π-conjugate system and strengthen the donor-acceptor character<sup>21,25-30</sup>. In this work the effect of the metal transition complexation on the intramolecular charge transfer as well as the role of auxiliary ligands have been investigated. The relations between theoretically calculated NLO properties and the experimental obtained second and third order susceptibilities are presented. Both results indicate that these molecular materials are valuable candidates for application in optoelectronics and photonics<sup>31</sup>.

## Experimental studies

All new compounds are air stable in the solid state and reasonably air stable in solution, but the vinylidene intermediates are not air stable and all reactions were performed under an argon atmosphere. All the solvents were of HPLC grade, further purified in a solvent system containing drying columns or dried over 4 Å molecular sieves. All commercially available reagents were used without further purification. The starting materials, ferrocenylacetylene<sup>32</sup>, cis-[Ru(dppm)<sub>2</sub>Cl<sub>2</sub>]<sup>33</sup>, cis-[RuCl<sub>2</sub>(PPh<sub>3</sub>)<sub>2</sub>(4,4'-dimethyl-2,2'-bipyridine)]<sup>34</sup> and RuCl(pppe) η<sup>5</sup>-Cp<sup>35</sup>, (E)-N,N-dibutyl-4-((4-bromophenyl)diazenyl)aniline<sup>5,36</sup> and (E)-N,N-dibutyl-4-((4-ethynylphenyl) diazenyl)aniline<sup>21,37</sup>, 6 were prepared by literature methods.

Flash column chromatography was performed on silica gel (high-purity grade, 230–400 mesh, 40–63 µm, Sigma-Aldrich) according to a standard technique. Nuclear magnetic resonance spectra (¹H, ¹³C and ³¹P) were recorded on a Bruker spectrometer (400 MHz). Chemical shifts for ¹H and ¹³C spectra are recorded in parts per million and are calibrated to solvent residual peaks (for example: CHCl₃: ¹H 7.26 ppm, ¹³C 77.16 ppm) according to H. E. Gottlieb, V. Kotlyar and A. Nudelman<sup>38</sup>. Multiplicities are indicated by s (singlet), bs (broad singlet), d (doublet), t (triplet), q (quadruplet), quint (quintuplet) and m (multiplet). Coupling constants, J, are reported in Hertz. Exact mass was obtained through Matrix Assisted Laser Desorption Ionization Time Of Flight Mass Spectrometry (MALDI-TOF MS).

**Syntheses.** Synthesis of 1: 250 mg (1.2 mmol) of ferrocenylacetylene<sup>32</sup> and 466 mg (1.2 mmol) of 5<sup>36</sup> were dissolved in 20 mL of THF and 20 mL of Et<sub>3</sub>N in a round bottom flask. This flask was purged three times with argon. 11 mg (60 µmol) of copper iodide and 41 mg (60 µmol) of bis(triphenylphosphine)palladium(II) dichloride were added and the reaction mixture was stirred overnight at 50 °C. After the solvent was removed, the residue was extracted 3 times with 25 mL of CH<sub>2</sub>Cl<sub>2</sub> and washed with 50 mL of water, and then dried over Na<sub>2</sub>SO<sub>4</sub>. The resulting mixture was concentrated under reduced pressure. The crude product was purified by chromatography on silica gel eluted with ether/petroleum ether (1/4, v/v) to give 1. Red-orange crystalline powder, yield 366 mg, 59%. ¹H NMR (400 MHz, CDCl<sub>3</sub>) δ 7.89 (d, 7.8 Hz, 2H), 7.84 (d, 7.8 Hz, 2H), 7.62 (d, 8.0 Hz, 2H), 6.72, (d, 8.0 Hz, 2H), 4.56 (s, 2H), 4.29 (br. s., 7H), 3.39 (t, 7.8 Hz, 4H), 1.66 (m, 4H), 1.43 (m, 4H), 1.01 (t, 7.7 Hz, 6H). ¹³C NMR (101 MHz, CDCl<sub>3</sub>) δ 152.22, 150.69, 143.21, 132.04, 125.37, 124.59, 122.19, 111.13, 90.25, 86.12, 71.49, 70.03, 68.96, 65.25, 50.99, 29.53, 20.34, 14.01. MALDI-TOF MS calculated for C<sub>32</sub>H<sub>35</sub>FeN<sub>3</sub> m/z: 517.22, found: 517.56.

Preparation of 2: cis-[Ru(dppm)<sub>2</sub>Cl<sub>2</sub>]<sup>33</sup>, (300 mg, 0.22 mmol) was added to 6<sup>21,37</sup> (74 mg, 0.22 mmol) and KPF<sub>6</sub> (0.6 mmol) in CH<sub>2</sub>Cl<sub>2</sub> (25 mL) and stirred for 16 h. The dark red solution was filtered and the solvent removed in vacuo. The solid residue was washed with deoxygenated diethylether (3 × 10 mL) to remove any excess of 6 and then redissolved in CH<sub>2</sub>Cl<sub>2</sub>, and K<sub>2</sub>CO<sub>3</sub> (0.6 mmol) was added to the vinylidene solution (³¹P<sup>1</sup>H)NMR: δ (ppm) -16.7 (s, PPh<sub>2</sub>)) and the stirring continued for a further 2 h. The resulting red-orange solution was filtered before 20 mL of heptane was added. The solid thus obtained was isolated by filtration, washed with pentane, and dried under vacuum to give the product as a red-orange crystalline powder. Complex 2, yield 163 mg, 60%. ³¹P NMR (162 MHz, CDCl<sub>3</sub>) δ -6.7 ppm. ¹H NMR (400 MHz, CDCl<sub>3</sub>) δ 7.81 (d, 7.8 Hz, 2H), 7.53-7.44 (m, 16H) 7.32 (t, 8.0 Hz, 4H), 7.30 (t, 8.0 Hz, 4H), 7.28 (t, 8.0 Hz, 8H), 7.21(t, 8.0 Hz, 8H), 7.11 (t, 7.2 Hz, 2H), 6.72 (d, 7.8 Hz, 2H), 6.14 (d, 7.8 Hz, 2H), 4.96 (m., 4H), 3.38 (d, 7.8 Hz, 2H), 1.66 (m., 4H) 1.42 (m., 4H), 1.02 (t., 7.0 Hz, 6H).
¹³C NMR (101 MHz, CDCl₃) δ 149.8, 148.6, 143.5, 134.2 (quint, J<sub>CP</sub> = 11 Hz), 133.8 (quint, J<sub>CP</sub> = 11 Hz), 133.4, 132.3, 129.28, 127.81, 127.56, 124.54, 121.22, 114.59, 111.13, 50.97, 50.32 (t, J<sub>CP</sub> = 10 Hz), 29.55, 20.36, 14.02. Anal. Calcd for C₇₂H₇₀ClN₃P₄Ru: C, 69.87; H, 5.70; N, 3.39. Found: C, 69.05; H, 5.93, N, 3.14.

**Preparation of 3:** A deoxygenated solution of the complex cis-[RuCl₂(PPh₃)₂(4,4'-dimethyl-2,2'-bipyridine)]³⁴ (200 mg, 0.24 mmol), 6 (135 mg, 0.4 mmol) and KPF₆ (43.1 mg, 0.23 mmol) in 20 mL of CH₂Cl₂/MeOH (1/1) was stirred for 48 h at room temperature. The mixture was filtered off in celite column, concentrated to ca. 1 mL and precipitated with heptane. The solid obtained was washed with pentane, dried under vacuum and dissolved in CH₂Cl₂ (15 mL) (³¹P{¹H} NMR (162 MHz, CDCl₃): 20.9 ppm (singlet, PPh₃), ~143.92 (septet, PF₆, ¹Jₚf = 710 Hz). K₂CO₃ (2 mmol) was added and the reaction mixture was stirred for another 6 h. The crude product was purified by chromatography on silica gel eluted with CH₂Cl₂/petroleum ether (1/2, v/v) to give a red-orange solid identified as 3. Yield 130 mg, 45%. ³¹P NMR (162 MHz, CDCl₃) δ, 29.5 ppm. ¹H NMR (400 MHz, CDCl₃) δ 8.90 (d, 5.5 Hz, 1H), 8.22 (d, 5.5 Hz, 1H), 7.81 (d, 7.8 Hz, 2H), 7.72-7.46 (unresolved m., 20H), 7.16 - 6.93 (unresolved m.,16H), 6.72 (d, 7.8 Hz, 2H), 6.63 (d, 5.5 Hz, 1H), 5.95 (d, 5.5 Hz, 1H), 3.38 (t, 7.2 Hz, 4H), 2.33 (s, 3H), 2.21 (s, 3H), 1.65 (s, 4H), 1.42 (s, 4H), 1.00 (s, 7.8 Hz, 6H). ¹³C NMR (101 MHz, CDCl₃) δ 156.59, 155.53, 154.88, 151.93, 149.72, 144.73, 143.51, 142.60, 134.18 (t, 4.2 Hz), 133.95, 133.03, 132.14, 132.05, 131.97, 131.94, 128.57, 128.45, 128.04, 127.05 (t, 4.0 Hz), 124.55, 120.63, 120.37, 111.14, 50.95, 29.56, 21.03, 20.65, 20.35, 14.02. Anal. Calcd for C₇₀H₆₈ClN₅P₂Ru: C, 71.38; H, 5.82; N, 5.95. Found: C, 71.96; H, 5.58; N, 5.74.

**Preparation of 4:** A solution of RuCl(pppe)η²-Cp³⁵, (380 mg, 0.5 mmol), 6 (200 mg, 0.6 mmol), and KPF₆ (2.5 mmol) in 15 mL of CH₂Cl₂ was stirred for 48 h at room temperature under argon. Deoxygenated diethyl ether (2×50 mL) was added and the mixture filtered. The filtrate was dissolved in CH₂Cl₂ (15 mL) ³¹P{¹H} NMR (singlet at 42 ppm (PPh₃), ~143.87 (septet, PF₆, ¹Jₚf = 709 Hz). Then K₂CO₃ (2 mmol) was added and the reaction mixture was stirred for another 6 h. The crude product was purified by chromatography on silica gel eluted with diethylether/petroleum ether (1/2, v/v) to give 4. Red solid, yield 250 mg, 48%. ³¹P NMR (162 MHz, CDCl₃) δ, 50.2 ppm (s, PPh₃). ¹H NMR (400 MHz, CDCl₃) δ 7.83 (d, 7.8 Hz, 2H), 7.70 (d, 7.8 Hz, 2H), 7.51 (m, 10H), 7.23 (m, 8H), 7.16 (m, 12H), 6.71 (d, 7.8 Hz, 2H), 4.37 (s, 5H) 3.38 (t, 7.0 Hz, 4H), 1.66 (m, 4H), 1.42 (m, 4H), 1.01 (t, 7.1 Hz, 6H). ¹³C NMR (101 MHz, CDCl₃) δ 149.88, 149.19, 143.47, 138.77 (m, ¹Jₚ₄ + ³Jₚ₄ = 42 Hz), 133.82, 130.93, 128.81, 128.47, 127.27, 124.66, 122.57, 122.02, 116.33, 111.13, 85.32, 50.96, 29.51, 20.35, 14.01. Anal. Calcd for C₆₃H₆₁N₃P₄Ru: C, 71.78; H, 5.83; N, 3.99. Found: C, 71.96; H, 5.58; N, 3.74.

**Electrochemical experiments.** Electrochemical experiments were carried out with a Biologic SP-150 potentiostat driven by the EC-Lab software including ohmic drop compensation. Cyclic Voltammetry (CV) was performed in a three-electrode cell controlled at a temperature of 293 K in a glove box containing dry, oxygen-free (< 1 ppm) argon. Working electrodes were glassy carbon planar disk electrodes (Ø = 3 mm). Counter electrodes were platinum wires. Reference electrodes were Ag/AgNO₃ (0.01 M CH₃CN). Experiments were recorded in dry HPLC-grade acetonitrile with tetrabutylammonium hexafluorophosphate (Bu₄NPF₆, electrochemical grade, Fluka) as supporting electrolyte. All the potential reported were calibrated versus ferricinium/ferrocene couple (Fc⁺/Fc) (IUPAC Recommendation)³⁹. Based on repetitive measurements, absolute errors on potentials were found to be around ± 5 mV.

**Preparation of host–guest films of complexes and optical absorption measurements.** For the host–guest films preparation, the various chromophores and PMMA (Sigma-Aldrich, Mw = 15,000 g/mol) have been dissolved in dichloromethane at concentration of 50 g/L. The concentration of the compounds was 100 µmol towards 1 g of PMMA. The solutions were deposited using the spin-coater (SCS G3) at 1,000 rpm on BK 7 glass plates substrates with 1 mm thickness (which were cleaned in distilled water using ultrasonic bath, acetone, and ethanol) and then were dried. Obtained guest–host polymer films were kept at room temperature during two days in order to eliminate any remaining of solvent. The thickness of deposited films was calculated by the profilometer (Dektak 6M, Veeco).

We used Lambda 950 UV/Vis/NIR spectrophotometer (PerkinElmer) with the range 300–1,200 nm in order to measure the absorption spectra of the samples.

## NLO studies

**SHG and THG measurements.** After the preparation of the thin films, we used the corona poling technique in order to orientate the chromophores at the thin layer surface. First, they were heated on a hot plate at selected poling temperature of 100 °C. Then an external electric field was provided by applying a voltage of +5 kV to a tungsten needle fixed at 1 cm above the polymer surface, while the electrode under the glass substrate was grounded. With the remaining electric field, the heater was switched off, and the sample was cooled down to room temperature, the corona field was turned off.

In order to evaluate the nonlinear optical response of our metal complexes 1, 2, 3 and 4, we use the Third and Second Harmonic Generation (SHG & THG), which are based on the Maker Fringe technique⁴⁰. In this Harmonic Generation technique, an incident laser beam at the frequency ω interacts with a nonlinear medium in order to generate another beam with double and triple frequencies.

We used a Nd:YAG laser working at 1,064 nm with 30 ps pulse duration and we employed 10 Hz repetition rate (Fig. 1). The input energy of laser pulses was controlled by laser power/energy meter (LabMax TOP, COHERENT) to be 95 µJ for SHG and 140 µJ for THG measurements respectively. The thin films were mounted on a rotating stage and the beam was focalized on them. A fast photodiode measured a portion of the input beam in order to synchronize the rotation and the fundamental beam. The interference filter at 532 nm (or 355 nm) was used to select the desired wavelength of light²⁰,⁴¹⁻⁴⁴ and was used to cut the pump beam before the photomultiplier, that allowed the third harmonic or the second harmonic generated signal to be detected. The polarization
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Experimental setup for the THG and SHG measurements: (BS) the beam splitters, (P) the Glan polarizers, (PhD) the photodiode, (λ/2) the half wave plate, (L) the lens, (S) the sample, (RS) the rotation stage, (F) the neutral density filters, (KG3) the KG3 filter, (IF) the interference filter, (PMT) the photomultiplier tube.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Synthesis of the ferrocenyl derivative 1.</figcaption></figure>

of second and third harmonic was controlled by polarizer placed before the photomultiplier. We used established reference materials, quartz glass plate and silica glass plate respectively, for the second harmonic generation and third harmonic generation measurements$^{45}$.

**Theoretical quantum chemical calculations of NLO properties.** For nonlinear optical properties calculation of the first and second order hyperpolarizabilities (β, γ), the time-dependent Hartree-Fock (TDHF) level of theory with functional LanL2DZ for transition metals Ru Fe and 6-31+G(d,p) functional for the other concerned atoms (C, H, N, P, Cl) have been used. Indeed, in order to obtain good theoretical results we combined two methods (LANL2DZ functional$^{46}$ for transition metals and all-electron basis sets for all other non-transition-metal atoms), which is prevalent in computations for transition-metal-containing complexes. In order to get a precise estimation of the hyperpolarizabilities one should take into account d polarization functions on the carbon and nitrogen atoms and the addition of p functions on hydrogen atoms and diffuse functions. Calculations were performed for structures optimized by density functional theory DFT/B3LYP at the 6-31+G(d,p) + LANL2DZ mixed basis set. The hybrid-GGA functional B3LYP$^{47}$ gives the best performance in predicting the good structure and HOMO, LUMO properties. All calculations were done using GAUSSIAN 09 program package.

## Results and discussion

**Synthesis and characterization.** The catalytic coupling of ferrocenylacetylene$^{32}$ with (E)-N,N-dibutyl-4-((4-bromophenyl)diazenyl)aniline$^{36}$ in $NEt_3$ at 50 °C by $PdCl_2(PPh_3)_2/CuI$ catalysts gave the red-orange crystalline ferrocenyl derivative **1** (59%; Fig. 2). The acetylenic moiety was characterized by $^{13}C\{^1H\}$ NMR (δ 90.2 ppm (~$C_6H_4-C$; C-) and 86.1 (C; $C$-Cp)).

The synthesis of the ruthenium acetylide derivatives **2-4** was attempted by the activation of the azo-containing alkyne **6** (*E*)-N,N-dibutyl-4-((4-ethynylphenyl)diazenyl)aniline$^{21,37}$, in a two-steps procedure, via the formation of vinylidene intermediates (Fig. 3) followed by deprotonation, as was shown to occur for $RuCl_2(dppe)_2$ analogues$^{21,37}$. This approach is in contrast with classical metal-carbon bond formation via a metal-halide unit and an organometallic$^{48}$ and was preferred here because it supports easily functional groups that the alkyne bears. This route was first used to prepare the $\eta^1$-ethynyl trans-[RuCl(-C;C-R)(dppm)$_2$] **2**. cis-RuCl$_2$(dppm)$_2$$^{33}$ was reacted with **6** in the presence of $KPF_6$ in $CH_2Cl_2$ to give a pale green intermediate, likely to be the vinylidene $[Cl(dppm)Ru = C = CH-R]^+$ salt ($^{31}P\{^1H\}$NMR: δ (ppm) -16.7 (s, PPh$_2$)). On deprotonation by potassium carbonate, this intermediate led to the red ruthenium complex **2** in good yield (60%; Fig. 3) ($^{31}P\{^1H\}$NMR: δ (ppm) -6.7 (s, PPh$_2$)). A singlet resonance for the Ru-C;C carbon nuclei (δ 114.6 ppm) confirms the presence of the alkynyl ligand.

cis-[RuCl$_2$(PPh$_3$)$_2$(4,4'-dimethyl-2,2'-bipyridine)]$^{34}$ was used to prepare the $\eta^1$-ethynyl cis-[RuCl(-C;C-R)(PPh$_3$)$_2$(4,4'-dimethyl-2,2'-bipyridine)], **3**. In the first step, a vinylidene complex (see Fig. 3) is obtained by chloride displacement in a $CH_2Cl_2/MeOH$ 1:1, v:v solution with $PF_6^-$ salt and activation of the azo-containing
<figure><img src="image_5.png" /><figcaption>**Figure 3.** General synthetic scheme of ruthenium acetylenes derivatives.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Synthesis of the ruthenium organometallic complexes 2–4.</figcaption></figure>

alk-1-yne 6 H-C≡C-R, at 40 °C. The formation of this vinylidene complex is confirmed by ³¹P{¹H} NMR spectra that shows a singlet at 22.1 ppm as expected due to the equivalence of the opposite position of phosphorus atoms in a *trans* isomer. This vinylidene species decomposes slowly in air, due to reaction between the vinylidene group and atmospheric oxygen as observed for similar compounds⁴⁹. It is quickly deprotonated upon addition of potassium carbonate in dichloromethane to give the expected η¹-ethynyl derivative 3 (Fig. 4)⁴⁹.

Finally, the reaction between RuCl(PPh₃)₂(η⁵-C₅H₅) and the alkyne 6, in the presence of KPF₆, has given the corresponding cationic η¹-vinylidene complex [Ru(C=CHR)(PPh₃)₂(η⁵-C₅H₅)]⁺, the formation was attested by ³¹P{¹H} NMR (singlet at 42 ppm). This cation is readily deprotonated to give the corresponding η¹-ethynyl derivative Ru(-C≡C-R) (PPh₃)₂(η⁵-C₅H₅), 4 (³¹P{¹H} NMR (singlet at 50.2 ppm).

These new alkynyl complexes were characterized by MS or satisfactory microanalyses, UV-vis spectroscopy and ¹H-, ³¹P- and ¹³C-NMR spectroscopy (see part 1.1 and Supporting Information).

**Cyclic voltammetry.** In order to gain insight into the electronic environment of the metal atom in the new alkynyl complexes, the electrochemical properties of these compounds were measured by means of cyclic voltammetry in 0.1 M tetra-butylammonium hexafluorophosphate ([n-Bu₄N]PF₆) methylene chloride solutions. The results are summarized in Table 1. The cyclic voltammetric scan of compound 1 exhibits one reversible oxidation wave that corresponds to the ferrocene oxidation with a half wave potential of E₁/₂ = 0.10 V (vs. FeCp₂/FeCp₂⁺). This value can be compared to that of free ferrocene (i.e. internal reference used for calibration), which clearly reveals an electron transfer from donor ferrocene to the ethynyl unit, which behaves as a modest electron with-drawing group: the oxidation of the iron center becomes more anodic compared to free ferrocene due to the removal of electron density. A second anodic wave is observed at E₁/₂ = 0.50 V. This wave is reversible
<table><caption>**Table 1.** Cyclic voltammetry and linear optical data. Measured in CH₂Cl₂, (vs. FeCp/FeCp⁺). aPt disc working, Pt wire auxiliary, Ag reference electrode; ferrocene/ferrocenium redox couple as internal reference, CH₂Cl₂ containing 0.1 M [n-Bu₄N][PF₆], 20 °C, scan rate 100 mV s⁻¹ at 5.10⁻⁴ M.</caption><thead><tr><th>Complex</th><th>E1/2 [V] [ΔE, (I<sub>pc</sub>/I<sub>pa</sub>)]<sup>a</sup></th><th>λ<sub>max</sub> [nm] (ε [10<sup>3</sup> M<sup>-1</sup> cm<sup>-1</sup>])</th></tr></thead><tbody><tr><td>1</td><td>0.10 [0.70, 1]; 0.50 [0.70, 1]</td><td>461 (54)</td></tr><tr><td>2</td><td>-0.06 [0.80, 1]; 0.27 [0.70, 1]</td><td>500 (47)</td></tr><tr><td>3</td><td>-0.30 [0.70, 1]; 0.30 [0.70, 1]</td><td>514 (23)</td></tr><tr><td>4</td><td>-0.05 [0.60, 1]; 0.25 [0.70, 1]; 0.50 [0.80, 1]</td><td>481 (15)</td></tr></tbody></table>

<table><caption>**Table 2.** The values of χ<sup>(2)</sup><sub>eff</sub>, χ<sup>(3)</sup><sub>elec</sub>, and γ<sub>elec</sub> obtained for 1, 2, 3 and 4 guest-host films.</caption><thead><tr><th>Compounds</th><th>χ<sup>(2)</sup><sub>eff</sub>, s-p (pm V<sup>-1</sup>)</th><th>χ<sup>(2)</sup><sub>eff</sub>, p-p (pm V<sup>-1</sup>)</th><th>χ<sup>(3)</sup><sub>elec</sub>, s-s (10<sup>-20</sup> m<sup>2</sup>V<sup>-2</sup>)</th><th>γ<sub>elec</sub> (10<sup>-47</sup> m<sup>5</sup>V<sup>-2</sup>)</th></tr></thead><tbody><tr><td>1</td><td>0.15</td><td>1.46</td><td>0.61</td><td>5.55</td></tr><tr><td>2</td><td>0.09</td><td>0.69</td><td>0.35</td><td>3.22</td></tr><tr><td>3</td><td>0.05</td><td>0.35</td><td>0.25</td><td>2.26</td></tr><tr><td>4</td><td>0.99</td><td>5.41</td><td>0.27</td><td>2.44</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** UV-visible absorption spectroscopy of compounds 1–4.</figcaption></figure>

($I_{pa}/I_{pc} \cong 1$) ($I_{pc}$ and $I_{pa}$ are the peak cathodic and anodic currents), and is tentatively attributed to the oxidation of the azo-dialkylamino fragment.

Cyclic voltammetric scans of complexes 2 and 3 using window potentials of –1.0 to 1.0 V each show two electrochemically reversible waves with half-wave potentials of –0.06 V and +0.27 V for 2, and –0.30 V and +0.30 V for 3, (vs. FeCp₂/FeCp₂⁺). The lowest potential oxidation waves are attributed to the one electron oxidation of the ruthenium centers. These observations are consistent with previous results⁵⁰,⁵¹. Conversely, the highest oxidation potentials are attributed to the azo-dialkylamino fragment.

The cyclic voltammograms of the acetylide compound 4 shows a reversible oxidation wave ($I_{pa}/I_{pc} \cong 1$), corresponding to the Ru<sup>II/III</sup> couple at a potential [$E_{1/2} = -0.05$ V vs. FeCp₂/FeCp₂⁺] which is concordant with the data mentioned for RuX(PPh₃)₂($η^5-C_5H_5$) series⁵². A second reversible oxidation wave is observed at $E_{1/2} = 0.25$ V, that would correspond to the waves observed in 2 and 3, that we attributed to the oxidation of the azo-dialkylamino fragment. Surprisingly, a third wave, quasi reversible, was observed at 0.5 V vs. FeCp₂/FeCp₂⁺. This wave was tentatively attributed to a further oxidation of the azo-dialkylamino fragment to two-electron oxidation products⁵¹,⁵³.

**Photochemical studies in solution and in solid state.** The UV-vis absorption spectra of the azo-ethyl-nyl complexes 1–4 were recorded in dichloromethane, and the absorption data of these compounds are listed in Table 2. The ferrocenyl compound 1 (Fig. 5) shows a wide and strong absorption band between 390–530 nm, presumably due to superimposed π–π* and d–d transition arising from the ferrocene moiety⁵⁴⁻⁵⁷, together with the n–π* and π–π* absorption bands of the azobenzene chromophore.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Normalized UV-vis spectra of the compounds 1, 2, 3 and 4 incorporated in PMMA films at the concentration 100 µmol/g.</figcaption></figure>

The UV-Vis spectra for the acetylide derivatives **1–4** are summarized in Table **1** and presented in Fig. **5**. They are very similar with those of parent compounds mentioned in the literature.

The UV-vis spectrum of complex **2** contains a broad low energy non-symmetrical absorption band with a maximum centred around 514 nm. This low energy band is accompanied by a smooth shoulder at 400 to 430 nm with can be consistent with the $\pi-\pi^*$ absorption bands of the ethynyl-azo unit as observed in *trans*-[RuCl(-C:C-R)(dppe)₂]¹⁴, that overlaps with MLCT bands of the metal-acetylide unit.

Like the parent *cis*-[RuCl₂(PPh₃)₂(4,4'-dimethyl-2,2'-bipyridine)] complex³⁴,⁵¹, **3** has a major absorption band in the 450–550 nm range (λmax = 514 nm), that can been attributed to ruthenium-to-bipyridyl MLCT transitions⁵¹. This broad band is accompanied by a distinct shoulder with an apparent maximum at 400 nm which can be consistent with the absorption bands of the azobenzene chromophore.

The UV-vis spectra for the alkynyl-metal complex **4** shows a broad absorption band from 380 to 550 nm, that would correspond to the superimposition of both a MLCT band around 400 nm, that is concordant with the data mentioned for other Ru(-C:C-R) (PPh₃)₂(η⁵-C₅H₅) (R = H, CHO) complexes⁵², and the $\pi-\pi^*$ absorption bands of the azobenzene chromophore, as observed in the previous compounds.

The UV-Visible absorption spectra of the four investigated metal complexes dissolved in dichloromethane and then embedded in PMMA films at the concentration 100 µmol/g are given in Fig. 6. The spectra of the obtained films exhibit wide absorption band with maximum at 464 nm, 500 nm, 514 nm, and 481 nm respectively for the films of complexes **1, 2, 3** and **4** which correspond to ligand centered (LC, $\pi-\pi^*$ and n-$\pi^*$) transitions⁵⁸. They are absolutely similar to the absorption spectra obtained for the four investigated metal complexes dissolved in dichloromethane. More interestingly for the measurements of SHG and THG, we observed that the samples show high optical transparency, at the wavelengths higher than 700 nm.

**SHG and THG responses.** The measurements of SHG and THG of compounds **1, 2, 3** and **4** in PMMA films at the concentration 100 µmol/g, were obtained by using the Maker fringe technique for s- and p-polarized fundamental beam. We used the corona poling just before starting the measurements in order to increase a uniaxial orientation of the molecules and consequently enhance the 2nd order NLO properties of the compounds in the polymer films. The obtained SHG response is due to the break of the centrosymmetry of the studied organometallic complexes. Figure 7 represents the dependence of the second harmonic intensity generated as a function of the incident angle. As shown in this figure, the four films have a maximum signal between 60° and 65° and zero intensity at normal incidence of the fundamental beam. The intensity depends on the polarization p or s. In our case, the polarization of the second harmonic signal was found to be always p-polarized regardless the incident polarization.

For the calculation of the quadratic NLO susceptibilities, we used the comparative model of Lee and coworkers which compares the intensities of SHG for films and for quartz plate which is expressed by the following equation and which takes into account the linear optical absorption⁵⁹.

$$ \chi^{(2)} = \chi_{\text{Quartz}}^{(2)} \left(\frac{2}{\pi}\right) \left(\frac{I_{\text{Quartz}}^{\text{coh}}}{d}\right) \left(\frac{\frac{\alpha d}{2}}{1 - e^{-\frac{\alpha d}{2}}}\right) \left(\frac{I_{\text{Quartz}}^{2\omega}}{I_{\text{Quartz}}^{2\omega}}\right)^{1/2} \quad (1) $$

where $\chi_{\text{Quartz}}^{(2)} = 1.0\text{pm/v}$ is the quadratic susceptibility of quartz, $I_{\text{Quartz}}^{\text{coh}}$ is the coherent length of quartz, $\alpha$ is the absorption coefficient at doubled laser wavelength, $d$ is the film thickness, $I_{\text{Quartz}}^{2\omega}$ and $I_{\text{Quartz}}^{2\omega}$ are the SHG intensities for the sample and quartz respectively under the same conditions of measurement. The obtained quadratic NLO susceptibilities for **1–4** films are presented in Table **2**. For the input-output values of $\chi^{(2)}$ were found to be
<figure><img src="image_9.png" /><figcaption>**Figure 7.** SHG intensity as a function of incident angle in 1, 2, 3 and 4 guest-host films at the s-, and p-polarized fundamental beam.</figcaption></figure>

higher at polarization p-p than polarization s-p, which is caused by the symmetry peculiarities of the guest-host polymeric films after poling. On the basis of these results (Table 2) the film of complex 4 was found to exhibit higher second-order optical nonlinearity than the three other films. This result might be due to different polarity and hence different charge transfer from the metal to the π-conjugated acetylide azobenzene system (MLCT) in complex 4 as compared with the other two ruthenium metal complexes. This is intrinsic to the presence of bipyridine (3) or cyclopentadienyl (4) pi-acceptor ligands, whereas the dppm (2) and PPh₃ ligands are sigma-donor
<figure><img src="image_10.png" /><figcaption>Figure 8. THG intensity as a function of incident angle in 1, 2, 3 and 4 guest-host films.</figcaption></figure>

ligands. The structure of 1 is different since there is no metal-carbon bond between iron and acetylide. The charge transfer will be weak there, as the dibutylamino moiety acts as a strong electron donor unit. In fact in complex 3 the CV and the UV-visible measurements indicate a charge transfer from ruthenium to the bipyridyl unit⁵⁰.

For the calculation of cubic NLO susceptibilities, we used the model of Kubodera and Kobayashi. This model compares directly the maximum amplitudes of third harmonic light intensities of the studied material with those of flat silica which is taken a reference (1 mm in thickness). This model takes into account the optical absorption of the studied compounds⁶⁰:

$$ \chi^{(3)} = \chi_{\text{Silica}}^{(3)} \left(\frac{2}{\pi}\right) \left(\frac{I_{\text{Silica}}^{\text{coh}}}{d}\right) \left(\frac{\frac{\alpha d}{2}}{1 - e^{-\frac{\alpha d}{2}}}\right) \left(\frac{I_{\text{Silica}}^{3\omega}}{I_{\text{Silica}}^{3\omega}}\right)^{1/2} \quad (2) $$

where $\chi_{\text{Silica}}^{(3)}$ is the cubic susceptibility of silica, $I_{\text{Silica}}^{\text{coh}}$ is the coherent length of silica, $\alpha$ is the absorption coefficient at triple laser wavelength, $d$ is the film thickness, $I_{\text{Silica}}^{3\omega}$ and $I_{\text{Silica}}^{3\omega}$ are the THG intensities of the samples and silica at the same conditions, respectively. The Fig. 8 represents the dependences of the third harmonic intensity generated as function of the incident angle.

The calculated values of $\chi^{(3)}$ for 1–4 films are given in Table 2. We have the same results before and after corona poling which is expected as the orientation has no significant influence on third-order NLO properties. Due to the extended $\pi$-conjugated chain complexes and coupling between the metal⁶¹ and the remote groups through the $\pi$-conjugated path of the investigated metal complexes, the third-order nonlinear optical susceptibilities values of $\chi^{(3)}$ (Table 2) are stronger than the one of the silica⁶². The complex 1 is characterized with a higher third order NLO susceptibility than 2, 3 and 4 due to more effective electronic delocalization which is caused by charge transfer from the (N,N-dibutylamine) and ferrocenyl donor groups to the azo accepting unit. from third order susceptibility (Table 2) and taking into account the concentration of the complexes and the local field factor, we calculated the electronic contribution of second hyperpolarizability $\gamma_{elec}$ and we found that third order susceptibility, $\gamma_{elec}$ for 1 is also much higher than the other complexes⁶²,⁶³.

Although the experimental methods used in the literature to determine the $\chi^{(2)}$ and $\chi^{(3)}$ susceptibilities are slightly different, a comparison between the performances of our four organometallic complexes 1–4 and the ones reported for the state-of-the-art SHG/THG compounds could be made. However we have to underline here that for the third harmonic generation as we investigate thin film so what we are measuring is electronic contributions which is not the case for example for degenerate four waves mixing but comparison could be done. From Table 3, we can clearly see that the $\chi^{(3)}$ susceptibility values are slightly lower in the case of our organometallic complexes while the $\chi^{(2)}$ susceptibility values are generally much higher for our complexes 1–4. Thus, the sizeable SHG response for the complexes 1 and 4 indicates that these compounds are valuable candidates for potential applications in NLO and more specifically, for the evolution towards applications in molecular electronics and photonics and in related fields including second-order optical effects in organometallic compounds induced by an acoustic field²⁶,⁶⁴. Considering third harmonic generation, electronic contributions were measured as we investigated thin films. This is obviously a difference with results previously reported by using the degenerate four waves mixing approach⁶⁵⁻⁷¹.

**Hyperpolarizabilities and HOMO, LUMO analysis.** To engineer optimal nonlinear optical properties, the origins of the nonlinear optical phenomena must be understood. This requires knowledge of the bonding properties of the atoms in the molecules. To understand the connection between molecular structure and NLO property, we have to expand the scope of our investigation to include the computing the first and second order
<table><caption>**Table 3.** Comparison of the SHG/THG NLO properties of 1–4 compounds with the values of some selected compound in the literature (AISHG: acoustically induced second harmonic generation/DFWM: degenerate four wave mixing/AIOSHG: acoustically induced optical second harmonic generation/EFISH: electric field induced second harmonic).</caption><thead><tr><th>Compounds</th><th>λ(nm)</th><th>Technique</th><th>χ<sup>(2)</sup> (pm V<sup>-1</sup>)</th><th>χ<sup>(3)</sup> (10<sup>-20</sup> m<sup>2</sup> V<sup>-2</sup>)</th></tr></thead><tbody><tr><td>1</td><td>1,064</td><td>SHG/THG</td><td>1.46</td><td>0.61</td></tr><tr><td>2</td><td>1,064</td><td>SHG/THG</td><td>0.69</td><td>0.35</td></tr><tr><td>3</td><td>1,064</td><td>SHG/THG</td><td>0.35</td><td>0.25</td></tr><tr><td>4</td><td>1,064</td><td>SHG/THG</td><td>5.41</td><td>0.27</td></tr><tr><td>Trans Ru(2,5-C≡CC-th-CHO)Cl(dppe)<sub>2</sub><sup>26</sup></td><td>1,760/532</td><td>AISHG</td><td>0.4</td><td>–</td></tr><tr><td>Trans Ru(2,2',5,5'-C≡C-th-th-CHO)Cl(dppe)<sub>2</sub><sup>12,21,26</sup></td><td>1,760/532</td><td>AISHG/DFWM</td><td>0.8</td><td>1.2</td></tr><tr><td>Trans Ru(2,2',5,5',2'',5''-C≡C-th-th-th-CHO)Cl(dppe)<sub>2</sub><sup>26</sup></td><td>1,760/532</td><td>AISHG/DFWM</td><td>0.27</td><td>2.10</td></tr><tr><td>Trans Ru(-C≡C-th-(E)CH=CH-th-CHO)Cl(dppe)<sub>2</sub><sup>26</sup></td><td>1,760/532</td><td>AISHG/DFWM</td><td>0.64</td><td>2.50</td></tr><tr><td>C<sub>30</sub>H<sub>28</sub>N<sub>4</sub>O<sub>4</sub>Sn<sup>72</sup></td><td>1,300</td><td>THG</td><td>–</td><td>1.99</td></tr><tr><td>Pb<sub>3</sub>Ge<sub>5</sub>O<sub>12</sub><sup>26</sup></td><td>1,760</td><td>AIOSHG</td><td>0.61</td><td>–</td></tr><tr><td>Ru-X(CuCR)(dppe)<sub>2</sub><sup>21</sup></td><td>1,064</td><td>EFISH</td><td>1.4</td><td>–</td></tr><tr><td>LiNbO<sub>3</sub><sup>26</sup></td><td>1,760</td><td>AIOSHG</td><td>0.12</td><td>–</td></tr></tbody></table>

<table><caption>**Table 4.** Some selected components of the frequency-dependent β (−2ω;ω,ω) values at ω = 0.042827 atomic unit (1,064 nm) for compounds 1–4.</caption><thead><tr><th>Compounds</th><th>β<sub>x</sub> × 10<sup>-30</sup> esu</th><th>β<sub>y</sub> × 10<sup>-30</sup> esu</th><th>β<sub>z</sub> × 10<sup>-30</sup> esu</th><th>β<sub>tot</sub> × 10<sup>-30</sup> esu</th></tr></thead><tbody><tr><td>1</td><td>0.803</td><td>-2.759</td><td>-212.401</td><td>212.420</td></tr><tr><td>2</td><td>-3.786</td><td>-1.414</td><td>-66.023</td><td>66.146</td></tr><tr><td>3</td><td>51.038</td><td>1.3107</td><td>-18.692</td><td>54.369</td></tr><tr><td>4</td><td>70.634</td><td>-1.444</td><td>-58.023</td><td>91.422</td></tr><tr><td>TransRu(4-CCC6H4CHO)Cl(dppe)<sub>2</sub></td><td></td><td></td><td></td><td>40<sup>74</sup></td></tr><tr><td>TransRu(4-CCC6H<sub>4</sub>NO<sub>2</sub>)Cl(dppe)<sub>2</sub></td><td></td><td></td><td></td><td>55<sup>74</sup></td></tr><tr><td>Organometallic complexes</td><td></td><td></td><td></td><td>10-1,300<sup>74</sup><br>5.7-140<sup>30</sup></td></tr></tbody></table>

<table><caption>**Table 5.** Some selected components of the frequency-dependent γ (−2ω;ω,ω,0) values at ω = 0.042827 atomic unit (1,064 nm) for compounds 1–4.</caption><thead><tr><th>Compounds</th><th>γ<sub>xxxx</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>yyyy</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>zzzz</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>xxyy</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>xxzz</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>yyzz</sub> × 10<sup>-36</sup> esu</th><th>γ<sub>tot</sub> × 10<sup>-36</sup> esu</th></tr></thead><tbody><tr><td>1</td><td>4.59</td><td>6.32</td><td>1,600.73</td><td>1.18</td><td>-11.54</td><td>8.57</td><td>321.59</td></tr><tr><td>2</td><td>52.59</td><td>49.03</td><td>1507.85</td><td>9.73</td><td>14.74</td><td>15.11</td><td>337.73</td></tr><tr><td>3</td><td>951.55</td><td>61.87</td><td>106.64</td><td>5.17</td><td>279.13</td><td>3.69</td><td>339.21</td></tr><tr><td>4</td><td>1,338.12</td><td>31.63</td><td>24.50</td><td>17.23</td><td>46.12</td><td>8.48</td><td>307.59</td></tr><tr><td>TransRu(4-CCC6H4CHO)Cl(dppe)<sub>2</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td>300<sup>75</sup></td></tr><tr><td>TransRu(4-CCC6H<sub>4</sub>NO<sub>2</sub>)Cl(dppe)<sub>2</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td>320<sup>75</sup></td></tr><tr><td>Organometallic complexes</td><td></td><td></td><td></td><td></td><td></td><td></td><td>9-360<sup>65</sup><br>5.5-100<sup>30</sup></td></tr></tbody></table>

hyperpolarizabilities (β, γ). The obtained theoretical results are collected in Tables 4, 5. The β<sub>tot</sub> value has been calculated using the following expression<sup>73</sup>.

$$ \beta_{tot} = \sqrt{\beta_{x}^{2} + \beta_{y}^{2} + \beta_{z}^{2}} \qquad (3) $$

Generally, an increase in the β value occurs together with the hypochromic shift effect (see Table 4) what is seen for compounds 1 and 4. Theoretical results suggest that better second order hyperpolarizabilities possess the ferrocenyl compound 1 and the ruthenium metal complex 4. These results are in good agreement with the experimental results, where the best λ<sup>(2)</sup> were obtained for molecules 1 and 4. The HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital) are very important aspects to consider the NLO properties. In general if the HOMO–LUMO gap decreases (molecule 2–3), thus causing a bathochromic
<figure><img src="image_11.png" /><figcaption>Figure 9. The electric dipole moment vector of 1–4.</figcaption></figure>

shift the NLO properties should increase. In our case, there is no such relationship because absorption bands for molecules 2, 3 are close to resonance and the second harmonic signal can be partially absorbed.

The calculated hyperpolarizabilities of these molecules were compared with first order hyperpolarizabilities of other organometallic compounds (see Table 4 30,74). The study reveals that the investigated complexes have the same β range as other organometallics complexes known in literature hence in general may have potential applications in the development of non-linear optical materials.

The average second-order hyperpolarizability γtot values have been calculated using the following formula:

$$ \gamma_{tot} = \left(\frac{1}{5}\right) [\gamma_{xxxx} + \gamma_{yyyy} + \gamma_{zzzz} + 2(\gamma_{xxyy} + \gamma_{xxzz} + \gamma_{yyzz})] \quad (4) $$

The calculated second order hyperpolarizability tensor components are presented in Table 5. The hyperpolarizability of γtot obtained for molecules 1–4 are quite large comparable to the values obtained for other organometallic materials as presented in Table 5 30,61. The second order frequency-dependent hyperpolarizability for compound 1 and 2 is dominated by the longitudinal component of γzzzz, but for the molecule 3 and 4 the maximal tensor component is achieved for γxxxx. The same conclusions can be drawn in relation to the components obtained for the calculation of the first hyperpolarizability. Domination of particular component indicates a substantial delocalization of charges in these directions. This electron distribution can be skewed by substituents, the enlarge of this redistribution is measured by dipole moment, and facility of redistribution in response to an externally applied electric field by hyperpolarizability.

From Fig. 9 we can see the direction of the electric dipole moment. It can be deduced also from Table 6 that for the molecules 1 and 2, the dipole moment vectors are directed mainly along the longest part of the molecule μx while for molecules 3, 4, the dipole moment vectors in the direction of phosphine ligands μs. If we take into account the maximum values of γzzzz for 1 and 2 and γxxxx for 3, 4 a good agreement is obtained between theoretical and experimental results.

HOMO, LUMO and HOMO–LUMO-energy gaps have been also calculated in order to analyze the relationship of NLO properties with the molecular structure of the studied compounds. Interpretations of orbital energies give useful forecasting. As usually, the HOMO and LUMO energy levels depend on electron donating strength of donor and electron – withdrawing strength of acceptor, respectively.

The energy gap (ΔE$_{L-H}$ = E$_{LUMO}$ − E$_{HOMO}$), explains the charge transfer interaction taking place within the ligands. The HOMO represents the ability to donate an electron and the LUMO represents the ability to obtain an electron. The HOMO and LUMO energies and HOMO–LUMO energy gaps are given in Table 6. The width of the energy gap obtained for investigated materials falls within limits 2–3 eV. Which correlates well with the maximum in the absorption spectrum, see Table 1. The reduction in the HOMO–LUMO energy gap explains eventually charge transfer interaction taking place within the molecules. Generally, if conjugation bonds in a molecule increases, the HOMO–LUMO gap decreases and leads to a bathochromic shift. The HOMO–LUMO
<table><caption>**Table 6.** The calculated electric dipole moments (Debye) and dipole moment components for complexes 1–4. (DFT/B3LYP levels in 6-31G (d,p) + lanl2dz basis set).</caption><thead><tr><th>Compounds</th><th>μ<sub>x</sub> (Debye)</th><th>μ<sub>y</sub> (Debye)</th><th>μ<sub>z</sub> (Debye)</th><th>μ<sub>tot</sub> (Debye)</th></tr></thead><tbody><tr><td>1</td><td>0.15</td><td>0.65</td><td>-4.48</td><td>4.53</td></tr><tr><td>2</td><td>0.03</td><td>-0.33</td><td>5.21</td><td>5.21</td></tr><tr><td>3</td><td>9.70</td><td>0.44</td><td>5.10</td><td>10.96</td></tr><tr><td>4</td><td>4.90</td><td>-3.50</td><td>-1.53</td><td>6.42</td></tr></tbody></table>

<table><caption>**Table 7.** HOMO, LUMO energies (E<sub>HOMO</sub>, E<sub>LUMO</sub>) and HOMO–LUMO gaps (E<sub>LUMO</sub>–E<sub>HOMO</sub>) calculated at the DFT/B3LYP levels in 6-31G (d) and lanl2dz basis set.</caption><thead><tr><th>Compound</th><th>E<sub>HOMO</sub> (eV)</th><th>E<sub>LUMO</sub> (eV)</th><th>E<sub>LUMO</sub>-E<sub>HOMO</sub> (eV)</th></tr></thead><tbody><tr><td>Ligand</td><td>-7.24111</td><td>-1.70347</td><td>5.53764</td></tr><tr><td>1</td><td>-5.06633</td><td>-2.16036</td><td>2.90597</td></tr><tr><td>2</td><td>-4.48454</td><td>-1.74456</td><td>2.73998</td></tr><tr><td>3</td><td>-4.18738</td><td>-2.13342</td><td>2.05396</td></tr><tr><td>4</td><td>-5.519867</td><td>-2.51445</td><td>3.00541</td></tr></tbody></table>

energy gap decreases also upon metal complexation (see Table 7) E<sub>LUMO</sub>–E<sub>HOMO</sub> for free ligand is higher than for complexes 1–4 what should be related to increasing of NLO properties.

Based on the obtained experimental and theoretical results (see Supporting Information pages S10-S11) we can deduce that in our case the large second and third order optical nonlinearity is caused by strong intramolecular charge transfer between donors and acceptors through a conjugated transmitter.

The theoretical results show that the geometry of the molecule can determine nonlinear properties. The choice of the auxiliary ligands allows the direction of dipole moment of the complex to be modified and it can change the linear and nonlinear optical properties the whole system. Comparison of theoretical results with the values published in the literature indicates that all these complexes exhibit large microscopic and macroscopic second and third-order NLO properties.

## Conclusion

In summary, we have reported herein the synthesis and full characterization of four specific organometallic complexes (1–4). The UV–Visible absorption spectra as well as the cyclic voltammetry measurements indicated a different metal to ligand charge transfer (MLCT). Thus, the second and third-order nonlinear optical properties of these four organometallic complexes have been evaluated by means of second and third harmonic generation (SHG and THG) techniques. The values of the second order NLO susceptibility of our four complexes were found to be the highest for complex 4 while complex 1 showed the strongest third harmonic generation. The difference of obtained values of second and third nonlinear susceptibilities is due to modification of the energy states and the introduction of a charge transfer state upon complexation which are of different nature within the four metal complexes. Calculations correctly reproduce HOMO, LUMO energies and hyperpolarizabilities of the organometallic metal complexes 1–4. The calculated HOMO and LUMO energies show that charge transfer occurs within the molecules. From theoretical calculation, it is clear that the presence of two Ph<sub>3</sub>P acting as the most electron-acceptor groups in 3 and 4 complexes skew the direction of the dipole moment causing electron delocalization in other direction than for the molecule 1 and 2. The possible modulation of the nonlinear optical susceptibility by the nature of the metal atom used and its coordination sphere is an important phenomenon which opens new perspectives for both photonic and opto-electronic applications.

Received: 7 May 2020; Accepted: 10 August 2020

Published online: 17 September 2020

## References

1. Cole, J. M. & Ashcroft, C. M. Generic classification scheme for second-order dipolar nonlinear optical organometallic complexes that exhibit second harmonic generation. *J. Phys. Chem. A* https://doi.org/10.1021/acs.jpca.8b11687 (2018).

2. Popczyk, A. *et al*. Selected organometallic compounds for third order nonlinear optical application. *Nanomaterials* (Basel, Switzerland) **9** (2019).

3. Ashcroft, C. M. & Cole, J. M. 5—Molecular engineering of organic and organometallic second-order nonlinear optical materials. in *Handbook of Organic Materials for Electronic and Photonic Devices* (2nd Edn.) (ed. Ostroverkhova, O.) 139–176 (Woodhead Publishing, 2019). https://doi.org/10.1016/B978-0-08-102284-9.00005-X.

4. Kodikara, M. S., Stranger, R. & Humphrey, M. G. Long-range corrected DFT calculations of first hyperpolarizabilities and excitation energies of metal alkynyl complexes. *ChemPhysChem* **19**, 1537–1546 (2018).

5. Humphrey, M. G. Ruthenium alkynyl complexes in non-linear optics*. *Aust. J. Chem.* **71**, 731–742 (2018).
6. Zawadzka, A. et al. Diagnostic and control of linear and nonlinear optical effects in selected self-assembled metallophthalocyanine chlorides nanostructures. *Dyes Pigm.* **157**, 151–162 (2018).

7. Fillaut, J.-L. Design of ruthenium-alkynyl complexes as nonlinear optical chromophores. *Display Imaging* **2**, 115 (2016).

8. Amar, A. et al. Photoisomerisation in aminoazobenzene-substituted ruthenium(II) tris(bipyridine) complexes: Influence of the conjugation pathway. *Chem. Eur. J.* **21**, 8262–8270 (2015).

9. Dragonetti, C. et al. Novel fullerene platinum alkynyl complexes with high second-order nonlinear optical properties as a springboard for NLO-active polymer films. *Organometallics* **35**, 1015–1021 (2016).

10. Dragonetti, C. et al. The role of substituents on functionalized 1,10-phenanthroline in controlling the emission properties of cationic iridium(III) complexes of interest for electroluminescent devices. *Inorg. Chem.* **46**, 8533–8547 (2007).

11. Boigel, J. et al. Second-order NLO switches from molecules to polymer films based on photochromic cyclometalated platinum(II) complexes. *J. Am. Chem. Soc.* **136**, 5367–5375 (2014).

12. Colombo, A. et al. Unexpectedly high second-order nonlinear optical properties of simple Ru and Pt alkynyl complexes as an analytical springboard for NLO-active polymer films. *Chem. Commun.* **50**, 7986–7989 (2014).

13. Luc, J., Fillaut, J.-L., Niziol, J. & Sahraoui, B. Large third-order nonlinear optical properties of alkynyl ruthenium chromophore thin films using third harmonic generation. *J. Optoelectron. Adv. Mater.* **9**, 2826–2832 (2007).

14. Xu, L. et al. 4-Diphenylamino-phenyl substituted pyrazine: Nonlinear optical switching by protonation. *J. Mater. Chem. C* **3**, 9191–9196 (2015).

15. Si, J., Mitsuyu, T., Ye, P., Shen, Y. & Hirao, K. Optical poling and its application in optical storage of a polyimide film with high glass transition temperature. *Appl. Phys. Lett.* **72**, 762–764 (1998).

16. Sakurai, M., Koley, P. & Aono, M. Tunable magnetism of organometallic nanoclusters by graphene oxide on-surface chemistry. *Sci. Rep.* **9**, 1–7 (2019).

17. Roknuzzaman, M. et al. Electronic and optical properties of lead-free hybrid double perovskites for photovoltaic and optoelectronic applications. *Scientific Reports* **9**, 718 (2019).

18. Tonnelé, C., Champagne, B., Muccioli, L. & Castet, F. Nonlinear optical contrast in azobenzene-based self-assembled monolayers. *Chem. Mater.* **31**, 6759–6769 (2019).

19. Xie, Q. et al. Novel photo-controllable third-order nonlinear optical (NLO) switches based on azobenzene derivatives. *Dyes Pigm.* **170**, 107599 (2019).

20. Bouchouit, M. et al. Synthesis, spectral, theoretical calculations and optical properties performance of substituted-azobenzene dyes. *Opt. Quantum Electron.* **48** (2016).

21. Colombo, A. et al. Nonlinear optical properties of intriguing Ru σ-acetylide complexes and the use of a photocrosslinked polymer as a springboard to obtain SHG active thin films. *Dalton Trans.* **45**, 11052–11060 (2016).

22. Virkki, M., Tuominen, O., Kauranen, M. & Priimagi, A. Photoinduced nonlinear optical response in azobenzene-functionalized molecular glass. *Opt. Express* **24**, 4964–4971 (2016).

23. Kreicberga, J., Laipniece, L., Bērziņa, G. & Kampars, V. Synthesis of pyridinium betaine azo chromophores. *Chem. Heterocycl. Compd.* **46**, 438–444 (2010).

24. Liaros, N. et al. NLO response of photoswitchable azobenzene-based materials. *ChemPhysChem* **14**, 2961–2972 (2013).

25. Bella, S. D. Second-order nonlinear optical properties of transition metal complexes. *Chem. Soc. Rev.* **30**, 355–366 (2001).

26. Fillaut, J.-L. et al. Design and synthesis of ruthenium oligothienylacetylide complexes. New materials for acoustically induced nonlinear optics. *Organometallics* **24**, 687–695 (2005).

27. Gauthier, N. et al. Electron-rich iron/ruthenium arylalkynyl complexes for third-order nonlinear optics: Redox-switching between three states. *Chem. Eur. J.* **17**, 5561–5577 (2011).

28. Green, K. A., Cifuentes, M. P., Samoc, M. & Humphrey, M. G. Syntheses and NLO properties of metal alkynyl dendrimers. *Coord. Chem. Rev.* **255**, 2025–2038 (2011).

29. Grelaud, G., Cifuentes, M. P., Paul, F. & Humphrey, M. G. Group 8 metal alkynyl complexes for nonlinear optics. *J. Organomet. Chem.* **751**, 181–200 (2014).

30. Hurst, S. K. et al. Organometallic complexes for nonlinear optics. 22.1 Quadratic and cubic hyperpolarizabilities of trans-bis(bidentate phosphine)ruthenium σ-arylvinlidene and σ-arylalkynyl complexes. *Organometallics* **20**, 4664–4675 (2001).

31. Kulyk, B. et al. Optimization and diagnostic of nonlinear optical features of π-conjugated benzodifuran-based derivatives. *RSC Adv.* **6**, 14439–14447 (2016).

32. Chen, W. et al. Alkyne-protected ruthenium nanoparticles. *J. Phys. Chem. C* **114**, 18146–18152 (2010).

33. Chaudret, B., Delavaux, B. & Poilblanc, R. Bisdiphenylphosphinomethane in dinuclear complexes. *Coord. Chem. Rev.* **86**, 191–243 (1988).

34. Adams, C. J. Synthesis and reactions of dimethyltin dithiooxalate: A convenient dithiooxalate transfer reagent. *J. Chem. Soc. Dalton Trans.* **1545–1550** (2002) https://doi.org/10.1039/B109520J.

35. Gamasa, M. P. et al. Phosphine substitution in indenyl- and cyclopentadienylruthenium complexes. Effect of the η5 ligand in a dissociative pathway. *Organometallics* **15**, 302–308 (1996).

36. Lemercier, G., Alexandre, M., Mulatier, J.-C. & Andraud, C. Synthesis of a pentaerythritol derivative bearing azo functions. *J. Chem. Res.* **2003**, 542–543 (2003).

37. Gherab, K. N. et al. Design and photoinduced surface relief grating formation of photoresponsive azobenzene based molecular materials with ruthenium acetyliides. *J. Mater. Chem.* **20**, 2858 (2010).

38. Gottlieb, H. E., Kotlyar, V. & Nudelman, A. NMR chemical shifts of common laboratory solvents as trace impurities. *J. Organ. Chem.* **62**, 7512–7515 (1997).

39. Gritzner, G. & Kūta, J. Recommendations on reporting electrode potentials in nonaqueous solvents: IUPC commission on electrochemistry. *Electrochim. Acta* **29**, 869–873 (1984).

40. Maker, P. D., Terhune, R. W., Nisenoff, M. & Savage, C. M. Effects of dispersion and focusing on the production of optical harmonics. *Phys. Rev. Lett.* **8**, 21 (1962).

41. Arroudj, S. et al. Synthesis, spectral, optical properties and theoretical calculations on schiff bases ligands containing o-tolidine. *Opt. Mater.* **56**, 116–120 (2016).

42. Kulyk, B. et al. Second and third order nonlinear optical properties of nanostructured ZnO thin films deposited on α-BBO and LiNbO3. *Opt. Commun.* **281**, 6107–6111 (2008).

43. Taboukhat, S. et al. Pyrene-based iminopyridine ligand and its metal complexes for nonlinear optical performance. in *2019 21st International Conference on Transparent Optical Networks (ICTON)* 1–4 (IEEE, 2019). https://doi.org/10.1109/ICTON.2019.8840532.

44. Taboukhat, S. et al. Nonlinear optical properties of some selected highly conjugated molecules based on TTF for optoelectronics applications. in *2017 19th International Conference on Transparent Optical Networks (ICTON)* 1–4 (IEEE, 2017). https://doi.org/10.1109/ICTON.2017.8024902.

45. Kulyk, B. et al. Tuning the nonlinear optical properties of BODIPYs by functionalization with dimethylaminostyryl substituents. *Dyes Pigm.* **137**, 507–511 (2017).

46. Yang, Y., Weaver, M. N. & Kenneth M. Merz, J. Assessment of the “6–31 + G** + LANL2DZ” mixed basis set coupled with density functional theory methods and the effective core potential: prediction of heats of formation and ionization potentials for first-row-transition-metal complexes. https://pubs.acs.org/doi/pdf/10.1021/jp807643p (2009) https://doi.org/10.1021/jp807643p.
47. Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. *J. Chem. Phys.* **98**, 5648–5652 (1993).

48. Long, N. J. & Williams, C. K. Metal alkynyl σ complexes: Synthesis and materials. *Angew. Chem. Int. Ed.* **42**, 2586–2617 (2003).

49. Adams, C. J. & Pope, S. J. A. Diimine–acetylide compounds of ruthenium: The structural and spectroscopic effects of oxidation. *Inorg. Chem.* **43**, 3492–3499 (2004).

50. Batista, A. A. et al. Electrochemical and spectroscopic studies on RuCl2(PPh3)2(N)2 and RuCl2(PPh3)2(N–N) complexes (N = pyridine derivatives and N–N = phenanthroline or bipyridine derivatives). X-ray structure of RuCl2(PPh3)2(phen). *Polyhedron* **20**, 2123–2128 (2001).

51. McDonagh, A. M. et al. Organometallic complexes for nonlinear optics Part 20. Syntheses and molecular quadratic hyperpolarizabilities of alkynyl complexes derived from (E)-4,4'-HC CC6H4NNC6H4NO2. *J. Organomet. Chem.* **9** (2000).

52. Hurst, S. K. et al. Organometallic complexes for nonlinear optics: Part 23. Quadratic and cubic hyperpolarizabilities of acetylide and vinylidene complexes derived from protected and free formylphenylacetylenes. *J. Organomet. Chem.* **633**, 114–124 (2001).

53. Stiborová, M., Asfaw, B., Frei, E. & Schmeiser, H. H. Oxidation of azo dyes by peroxidase: Additional evidence of a one-electron mechanism of oxidation of dimethylaminoazobenzene and Sudan I (solvent yellow 14). *Collect. Czech. Chem. Commun.* **61**, 962–972 (1996).

54. Roberts, D. J. et al. The synthesis and characterisation of novel ferrocenyl polyphenylenes. *Dalton Trans.* **41**, 8850–8860 (2012).

55. Kulyk, B., Guichaoua, D., Ayadi, A., El-Ghayoury, A. & Sahraoui, B. Metal-induced efficient enhancement of nonlinear optical response in conjugated azo-based iminopyridine complexes. *Org. Electron.* **36**, 1–6 (2016).

56. Ayadi, A., Korchi, K. E., Guichaoua, D., Taboukhat, S. & El-Ghayoury, A. Azo-Based Ligands and Metal Complexes for NLO Applications. in *2019 21st International Conference on Transparent Optical Networks (ICTON)* 1–4 (IEEE, 2019). https://doi.org/10.1109/ICTON.2019.8840333.

57. Kulyk, B. et al. Nonlinear refraction and absorption activity of dimethylaminostyryl substituted BODIPY dyes. *RSC Adv.* **6**, 84854–84859 (2016).

58. Arroudj, S. et al. Effect of the complexation on the NLO electronic contribution in film based conjugated quinoline ligand. *Phys. B* **516**, 1–6 (2017).

59. Lee, G. J., Cha, S. W., Jeon, S. J., Jin, J. I. & Yoon, J. S. Second-order nonlinear optical properties of unpoled bent molecules in powder and in vacuum-deposited film. *J. Korean Phys. Soc.* **39**, 912–915 (2001).

60. Kubodera, K. & Kobayashi, H. Determination of third-order nonlinear optical susceptibilities for organic materials by third-harmonic generation. *Mol. Cryst. Liq. Cryst. Incorp. Nonlinear Opt.* **182**, 103–113 (1990).

61. Niziol, J., Fillaut, J.-L., Sniechowski, M., Khammar, F. & Sahraoui, B. Alkynyl-ruthenium based azo-organometallics: Thermal evolution of thin films microstructure. *Opt. Mater.* **34**, 1670–1676 (2012).

62. Gubler, U. & Bosshard, C. Optical third-harmonic generation of fused silica in gas atmosphere: Absolute value of the third-order nonlinear optical susceptibility $\chi^{\Lambda}(3)$. *Phys. Rev. B* **61**, 10702–10710 (2000).

63. Stegeman, G. I. & Stegeman, R. A. *Nonlinear Optics: Phenomena, Materials and Devices*. (Wiley, 2012).

64. Migalska-Zalas, A. et al. Second-order optical effects in organometallic nanocomposites induced by an acoustic field. *Phys. Rev. B* **71**, 035119 (2005).

65. Kajzar, F., Messier, J. & Rosilio, C. Nonlinear optical properties of thin films of polysilane. *J. Appl. Phys.* **60**, 3040–3044 (1986).

66. van Beek, J. B., Kajzar, F. & Albrecht, A. C. Resonant third-harmonic generation in all-trans β-carotene: The vibronic origins of the third-order nonlinear susceptibility in the visible region. *J. Chem. Phys.* **95**, 6400–6412 (1991).

67. Burland, D. M., Walsh, C. A., Kajzar, F. & Sentein, C. Comparison of hyperpolarizabilities obtained with different experimental methods and theoretical techniques. *J. Opt. Soc. Am. B* **8**, 2269–2281 (1991).

68. Sahraoui, B. & Rivoire, G. Degenerate four-wave mixing in absorbing isotropic media. *Opt. Commun.* **138**, 109–112 (1997).

69. Sahraoui, B., Phu, X. N., Sallé, M. & Gorgues, A. Electronic and nuclear contributions to the third-order nonlinear optical susceptibilities of new *p-N*, *N*-dimethylaniline tetrathaifulvalene derivatives. *Opt. Lett.* **23**, 1811–1813 (1998).

70. Iliopoulos, K. et al. Physical origin of the third order nonlinear optical response of orthogonal pyrrolo-tetrathiafulvalene derivatives. *Appl. Phys. Lett.* **97**, 101104 (2010).

71. Migalska-Zalas, A. et al. $\chi^{(2)}$ grating in Ru derivative chromophores incorporated within the PMMA polymer matrices. *J. Phys. Chem. B* **108**, 14942–14947 (2004).

72. García-López, M. C. et al. Microwave-assisted synthesis, third-order nonlinear optical properties, voltammetry cyclic and theoretical calculations of organotin compounds bearing push-pull Schiff bases. *J. Organomet. Chem.* **806**, 68–76 (2016).

73. Thanthiriwatte, K. S. & Nalin de Silva, K. M. Non-linear optical properties of novel fluorenyl derivatives-ab initio quantum chemical calculations. *J. Mol. Struct. Theochem.* **617**, 169–175 (2002).

74. Migalska-Zalas, A. Theoretical study of the effect of $\pi$-conjugated transmitter of D–$\pi$–A ruthenium systems on the quadratic NLO properties. *Opt. Quant. Electron.* **48**, 183 (2016).

75. Migalska-Zalas, A. Theoretical simulation of the third order nonlinear optical properties of some selected organometallics complexes. *Digest. J. Nanomater. Biostruct.* **3**, 1–8 (2008).

## Acknowledgements

The presented results are part of a project IMAGE that has received funding from European Union's Horizon 2020 research and innovation program under Marie Skłodowska-Curie grant agreement No. 778156. Said Taboukhat would like to acknowledge the French Embassy in Morocco for partly financial support of this study. Calculations have been carried out in Wroclaw Centre for Networking and Supercomputing (http://www.wcss.pl), grant no. 282.

## Author contributions

All authors reviewed the manuscript and participate equally to this research work.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-71769-2.

Correspondence and requests for materials should be addressed to B.S.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020