Benoît Champagne

Design and Characterization of Molecular Nonlinear Optical Switches

Nanoscale structures, including molecules, supramolecules, polymers, functionalized surfaces, and crystalline/amorphous solids, can commute between two or more forms, displaying contrasts in their nonlinear optical (NLO) properties. Because of this property, they have high potential for applications in data storage, signal processing, and sensing. As potential candidates for integration into responsive materials, scientists have been intensely studying organic and organometallic molecules with switchable first hyperpolarizability over the past two decades. As a result of this, researchers have been able to synthesize and characterize several families of molecular NLO switches that differ by the stimulus used to trigger the commutation. These stimuli can include light irradiation, pH variation, redox reaction, and ion recognition, among others. The design of multistate (including several switchable units) and multifunctional (triggered with different stimuli) systems has also motivated a large amount of work, aiming at the improvement of the storage capacity of optical memories or the diversification of the addressability of the devices. In complement to the synthesis of the compounds and the characterization of their NLO responses by means of hyper-Rayleigh scattering, quantum chemical calculations play a key role in the design of molecular switches with high first hyperpolarizability contrasts. Through the latter, we can gain a fundamental understanding of the various factors governing the efficiency of the switches. These are not easily accessible experimentally, and include donor/acceptor contributions, frequency dispersion, and solvent effects. In this Account, we illustrate the similarities of the experimental and theoretical tools to design and characterize highly efficient NLO switches but also the difficulties in comparing them. After providing a critical overview of the different theoretical approaches used for evaluating the first hyperpolarizabilities, we report two case studies in which theoretical simulations have provided guidelines to design NLO switches with improved efficiencies. The first example presents the joint theoretical/experimental characterization of a new family of multi-addressable NLO switches based on benzazolo-oxazolidine derivatives. The second focuses on the photoinduced commutation in merocyanine-spiropyran systems, where the significant NLO contrast could be exploited for metal cation identification in a new generation of multiusage sensing devices. Finally, we illustrate the impact of environment on the NLO switching properties, with examples based on the keto-enol equilibrium in anil derivatives. Through these representative examples, we demonstrate that the rational design of molecular NLO switches, which combines experimental and theoretical approaches, has reached maturity. Future challenges consist in extending the investigated objects to supramolecular architectures involving several NLO-responsive units, in order to exploit their cooperative effects for enhancing the NLO responses and contrasts.

Electron correlation effects on the first hyperpolarizability of push–pull π-conjugated systems

The first hyperpolarizability (β) of representative push-pull π-conjugated compounds has been calculated at several levels of approximation to assess the effects of electron correlation. First, the 6-31+G(d) basis set has been shown to give the best balance between accuracy and computational resources for a polyene linker whereas for polyyne linker, the 6-31G(d) basis set is already an optimal choice. As a result of cancellations between higher order contributions, the MP2 method turns out to be the method of choice to predict β of push-pull π-conjugated systems since it closely reproduces the values obtained with the reference CCSD(T) scheme. Moreover, the SDQ-MP4 and CCSD approaches provide rarely improved estimates over MP2 while the MP4 method does not represent an improvement over MP4-SDQ or the SCS-MP2 method, over MP2. Among density functional theory exchange-correlation functionals, LC-BLYP is reliable when characterizing the changes of first hyperpolarizability upon enlarging the π-conjugated linker or upon changing the polyyne linker into a polyene segment. Nevertheless, its reliability is very similar to what can be achieved with the Hartree-Fock method and the MP2 scheme is by far more accurate. On the other hand, the BLYP, B3LYP, and BHandHLYP functionals perform quantitatively better in a number of cases but the trends are poorly described. This is also the case of the B2-PLYP and mPW2-PLYP functionals, which are often the most accurate, though they underestimate the increase of β when going from polyyne to polyene linkers and overestimate the enhancement of β with chain length.

Basis set and electron correlation effects on the polarizability and second hyperpolarizability of model open-shell π-conjugated systems

The basis set and electron correlation effects on the static polarizability (alpha) and second hyperpolarizability (gamma) are investigated ab initio for two model open-shell pi-conjugated systems, the C(5)H(7) radical and the C(6)H(8) radical cation in their doublet state. Basis set investigations evidence that the linear and nonlinear responses of the radical cation necessitate the use of a less extended basis set than its neutral analog. Indeed, double-zeta-type basis sets supplemented by a set of d polarization functions but no diffuse functions already provide accurate (hyper)polarizabilities for C(6)H(8) whereas diffuse functions are compulsory for C(5)H(7), in particular, p diffuse functions. In addition to the 6-31G(*)+pd basis set, basis sets resulting from removing not necessary diffuse functions from the augmented correlation consistent polarized valence double zeta basis set have been shown to provide (hyper)polarizability values of similar quality as more extended basis sets such as augmented correlation consistent polarized valence triple zeta and doubly augmented correlation consistent polarized valence double zeta. Using the selected atomic basis sets, the (hyper)polarizabilities of these two model compounds are calculated at different levels of approximation in order to assess the impact of including electron correlation. As a function of the method of calculation antiparallel and parallel variations have been demonstrated for alpha and gamma of the two model compounds, respectively. For the polarizability, the unrestricted Hartree-Fock and unrestricted second-order Møller-Plesset methods bracket the reference value obtained at the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples level whereas the projected unrestricted second-order Møller-Plesset results are in much closer agreement with the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples values than the projected unrestricted Hartree-Fock results. Moreover, the differences between the restricted open-shell Hartree-Fock and restricted open-shell second-order Møller-Plesset methods are small. In what concerns the second hyperpolarizability, the unrestricted Hartree-Fock and unrestricted second-order Møller-Plesset values remain of similar quality while using spin-projected schemes fails for the charged system but performs nicely for the neutral one. The restricted open-shell schemes, and especially the restricted open-shell second-order Møller-Plesset method, provide for both compounds gamma values close to the results obtained at the unrestricted coupled cluster level including singles and doubles with a perturbative inclusion of the triples. Thus, to obtain well-converged alpha and gamma values at low-order electron correlation levels, the removal of spin contamination is a necessary but not a sufficient condition. Density-functional theory calculations of alpha and gamma have also been carried out using several exchange-correlation functionals. Those employing hybrid exchange-correlation functionals have been shown to reproduce fairly well the reference coupled cluster polarizability and second hyperpolarizability values. In addition, inclusion of Hartree-Fock exchange is of major importance for determining accurate polarizability whereas for the second hyperpolarizability the gradient corrections are large.

Two‐Way Molecular Switches with Large Nonlinear Optical Contrast

Molecular switches: Highly efficient acido- and photoswitchable frequency doublers (see scheme) based on the indolinooxazolidine core are studied by means of hyper-Rayleigh experiments and quantum-chemical calculations.To optimize the nonlinear optical (NLO) contrast, a series of indolinooxazolidine derivatives with electron-withdrawing substituents in the para position on the indolinic residue have been synthesized. Their linear and nonlinear optical properties have been characterized by UV-visible absorption and hyper-Rayleigh scattering measurements, as well as by ab initio calculations. The two-way photo- or pH-triggered switching mechanism has been demonstrated by comparing the absorption spectra of the zwitterionic and protonated open forms (POF). Hyper-Rayleigh measurements have revealed that the second-order NLO contrast between the closed indolinooxazolidine and the open pi-conjugated colored forms remain very large upon substitution. Theory and measurements show that for the POFs the amplitude of the first hyperpolarizability follows the Hammett parameters of the withdrawing groups. However, because the measurements are performed in resonance, to recover this behavior, elaborate procedures including homogeneous and inhomogeneous broadenings, as well as single-mode vibronic structures are necessary to extrapolate to the static limit.

Long-range corrected density functional theory study on static second hyperpolarizabilities of singlet diradical systems

Within the spin-unrestricted density functional theory (DFT) the long-range correction (LC) scheme combined with the Becke-Lee-Yang-Parr exchange-correlation functional, referred to as LC-UBLYP method, has been applied to the calculation of the second hyperpolarizability (gamma) of open-shell singlet diradical systems of increasing complexity and has demonstrated good performance: (i) for the simplest H(2) dissociation model, the gamma values calculated by the LC-UBLYP method significantly overshoot the full configuration interaction result but reproduce qualitatively the evolution of gamma as a function of the diradical character, (ii) for small singlet diradical 1,3-dipole systems, the diradical character dependence of gamma determined by the UCCSD and UCCSD(T) reference methods is reproduced semiquantitatively by the LC-UBLYP method except in the small diradical character region, where the spin-unrestricted solutions coincide with spin-restricted solutions, (iii) the LC-UBLYP method also closely reproduces the UCCSD(T) results on the diradical character dependence of gamma of the p-quinodimethane model system, particularly in the intermediate and large diradical character regions, whereas it shows an abrupt change for a diradical character (y) close to 0.2 originating from the triplet instability, (iv) the reliability of LC-UBLYP to reproduce reference coupled cluster results on open-shell singlet systems with intermediate and large diradical characters has also been substantiated in the case of gamma of 1,4-bis-(imidazol-2-ylidene)-cyclohexa-2,5-diene (BI2Y), then (v), for real systems built from a pair of phenalenyl radicals separated by a conjugated linker, the LC-UBLYP results have been found to closely match the UBHandHLYP values-which, for small systems are in good agreement with those obtained using correlated molecular orbital methods-whereas the UB3LYP results can be much different. These results are not only important from the viewpoint of an efficient determination of the nonlinear optical properties of open-shell singlet systems, but also from the viewpoint of defining new challenges for elaborating improved exchange-correlation functionals.

Time dependent density functional theory investigation of the resonance Raman properties of the julolidinemalononitrile push-pull chromophore in various solvents

The absorption and resonance Raman spectra have been investigated for the first excited state of the julolidinemalononitrile push-pull chromophore in cyclohexane, dichloromethane, and acetonitrile by means of time dependent density functional theory calculations. The effect of nonlocal exchange has been considered by using three different hybrid exchange-correlation functionals containing 20%, 35%, and 50% of exact Hartree-Fock exchange. The interactions with the solvent have been described by the polarizable continuum model. The short-time approximation expression has been used to evaluate the resonance Raman intensities, while the vibronic theory of resonance Raman spectroscopy has been employed to determine both the intensities and the excitation profiles. It is shown that a consistent description of the vibronic structure of the excited state and resonance Raman spectra can be obtained provided that an adequate amount, close to 35%, of exact exchange is included in the exchange-correlation functional. The effect of increasing the polarity of the solvent is well represented by the polarizable continuum model, both for the absorption spectra and resonance Raman intensities. In particular, these simulations can reproduce the observed variations of the 1560 cm(-1) band intensity and attribute them to elongations of a CC double bond upon electronic excitation. Moreover, the short-time approximation has been found sufficient to reproduce most of the results of the more evolved vibronic theory of resonance Raman spectroscopy, which includes summations over vibrational excited states, for both the spectral signatures and their solvent dependencies.

Assessment of time-dependent density functional schemes for computing the oscillator strengths of benzene, phenol, aniline, and fluorobenzene.

In present study the relevance of using the time-dependent density functional theory (DFT) within the adiabatic approximation for computing oscillator strengths (f) is assessed using different LDA, GGA, and hybrid exchange-correlation (XC) functionals. In particular, we focus on the lowest-energy valence excitations, dominating the UV/visible absorption spectra and originating from benzenelike HOMO(pi)-->LUMO(pi(*)) transitions, of several aromatic molecules: benzene, phenol, aniline, and fluorobenzene. The TDDFT values are compared to both experimental results obtained from gas phase measurements and to results determined using several ab initio schemes: random phase approximation (RPA), configuration interaction single (CIS), and a series of linear response coupled-cluster calculations, CCS, CC2, and CCSD. In particular, the effect of the amount of Hartree-Fock (HF) exchange in the functional is highlighted, whereas a basis set investigation demonstrates the need of including diffuse functions. So, the hybrid XC functionals--and particularly BHandHLYP--provide f values in good agreement with the highly correlated CCSD scheme while these can be strongly underestimated using pure DFT functionals. These results also display systematic behaviors: (i) larger f and squares of the transition dipole moments (mid R:mumid R:(2)) are associated with larger excitation energies (DeltaE); (ii) these relationships present generally a linear character with R>0.9 in least-squares fit procedures; (iii) larger amounts of HF exchange in the XC functional lead to larger f, R:mumid R:(2), as well as DeltaE values; (iv) these increases in f, mid R:mumid R:(2), and DeltaE are related to increased HOMO-LUMO character; and (v) these relationships are, however, not universal since the linear regression parameters (the slopes and the intercepts at the origin) depend on the system under investigation as well as on the nature of the excited state.

An analytical derivative procedure for the calculation of vibrational Raman optical activity spectra

We present an analytical time-dependent Hartree-Fock algorithm for the calculation of the derivatives of the electric dipole-magnetic dipole polarizability with respect to atomic Cartesian coordinates. Combined with analogous procedures to determine the derivatives of the electric dipole-electric dipole and electric dipole-electric quadrupole polarizabilities, it enables a fully analytical evaluation of the three frequency-dependent vibrational Raman optical activity (VROA) invariants within the harmonic approximation. The procedure employs traditional non-London atomic orbitals, and the gauge-origin dependence of the VROA intensities has, therefore, been assessed for the commonly used aug-cc-pVDZ and rDPS:3-21G basis sets.

Resonance Raman scattering of rhodamine 6G as calculated by time-dependent density functional theory: vibronic and solvent effects.

The geometries, UV-vis absorption spectra, and resonance Raman (RR) intensities have been determined for the S1 and S3 excited states of rhodamine 6G (R6G) in vacuum and ethanol by means of DFT/TDDFT methodologies with the aim of better understanding the structures and properties of the excited states. The RR spectra have been simulated from the vibronic theory of RR scattering as well as within the short-time approximation, while the solvent effects have been modeled using the polarizable continuum model. The S1 and S3 states of R6G present UV-vis absorption bands with similar vibronic structure, i.e., a shoulder at smaller wavelengths, although this shoulder is relatively more intense and more sensitive to the solvent in the case of S3. These differences are corroborated by the larger geometry relaxations upon excitation for S3 and the fact that the charge transfer of S3 is reduced in ethanol. Moreover, the differences between S1 and S3 are magnified when considering the RR spectra. On one hand, the RR spectrum of R6G in resonance with the S0 --> S1 transition presents many transitions of which the relative intensities strongly vary when the excitation wavelength gets closer to the maximum of absorption. The RR spectrum of R6G in resonance with S1 is however little influenced by the solvent. On the other hand, the RR spectrum of R6G in resonance with the S0 --> S3 transition displays only a few bands, strongly depends on the solvent, and is little affected when changing the excitation wavelength within the limits of the absorption band. As a consequence, the short-time approximation is suitable to reproduce the RR spectrum of R6G in resonance with S3 for a broad range of excitation wavelengths, whereas the vibronic theory approach is needed for describing the RR spectrum of R6G in resonance with S1 close to resonance.

Vibrational polarizability of polyacetylene chains

Within the double harmonic oscillator approximation, ab initio vibrational contributions to the static electric dipole polarizability have been computed for the all‐trans polyacetylene chains homologous series. Atomic basis set and electron correlation effects have been studied on ethylene and trans‐butadiene as well as the use of semiempirical Hamiltonians. By using the 6‐31G* atomic basis set within the Hartree–Fock and the second order Mo/ller–Plesset procedures as well as by using semiempirical Austin model 1 Hamiltonians, the evolution with chain length of the vibrational polarizability per structural unit has been investigated and compared to the electronic contribution. Although smaller than the electronic contributions (αvibration≊10% αelectronic), the longitudinal component to the polarizability presents a similar exaltation as chain length grows, but a slower saturation to an asymptotic value per unit cell. Inclusion of electron correlation via the second order Mo/ller–Plesset technique turns ou...