3,4-Dimethylenecyclobutene: A Building Block for Design of Macrocycles with Excited State Aromatic Low-Lying High-Spin States.

3,4-Dimethylenecyclobutene (DMCB) is an unusual isomer of benzene. Motivated by recent synthetic progress to substituted derivatives of this scaffold, we carried out a theoretical and computational analysis with a particular focus on the extent of (anti)aromatic character in the lowest excited states of different multiplicities. We found that the parent DMCB is non-aromatic in its singlet ground state (S0), lowest triplet state (T1), and lowest singlet excited state (S1), while it is aromatic in its lowest quintet state (Q1) as this state is represented by a triplet multiplicity cyclobutadiene (CBD) ring and two uncoupled same-spin methylene radicals. Interestingly, the Q1 state, despite having four unpaired electrons, is placed merely 4.8 eV above S0, and there is a corresponding singlet tetraradical 0.16 eV above. The DMCB is potentially a highly useful structural motif for the design of larger molecular entities with interesting optoelectronic properties. Here, we designed macrocycles composed of fused DMCB units, and according to our computations, two of these have low-lying nonet states (i. e., octaradical states) at energies merely 2.40 and 0.37 eV above their S0 states as a result of local Hückel- and Baird-aromatic character of individual 6π- and 4π-electron monocycles.

On the third-order nonlinear optical responses of cis and trans stilbenes - a quantum chemistry investigation.

The second hyperpolarizabilities (γ) of the stilbene molecular switch in its trans and cis forms have been calculated using quantum chemistry methods to address their third-order nonlinear optical contrasts, to assess the reliability of lower-cost DFT methods, and to make comparisons with experiments. First, the reference CCSD(T) method shows that trans-stilbene presents a γ‖ value twice larger than its cis isomer (its γTHS value is 2.7 times larger). Among more cost-effective methods, reliable results are obtained at MP2 as well as with DFT, provided the CAM-B3LYP or ωB97X-D XCFs are employed. Supplementary DFT calculations have investigated the relationships between the accuracy of the exchange-correlation functionals, the fulfillment of Koopmans' theorem, and the delocalization error, and they demonstrated that satisfying Koopmans' theorem is not the condition for the best accuracy but that functionals with small delocalization errors are generally efficient. Using the selected CAM-B3LYP, large γ enhancements by about 70% (trans-stilbene) and 50% (cis-stilbene) have been evidenced when accounting for solvent effects using an implicit solvation model (IEFPCM), even for apolar solvents. Then, the frequency dispersion of the γ responses has been described using Bishop polynomial expansions, allowing comparisons with a broad set of experimental data. To a certain extent, no systematic agreement between the calculations and the measured values was found. On the one hand, the agreement is satisfactory for the γ(-ω;ω,-ω,ω) quantities, provided that the dominant vibrational contribution is taken into account. On the other hand, the agreement is poor for the γ(-2ω;ω,ω,0) and γ(-3ω;ω,ω,ω) quantities, while some inconsistencies between experimental values are also highlighted.

Second-order nonlinear optical properties of X-shaped pyrazine derivatives.

This work reports an investigation of the second-order NLO properties of two isomer series of X-shaped pyrazine derivatives, by means of HRS measurements and DFT calculations. The systems differ in the relative position of the donor and acceptor substituents with respect to the axis formed by the nitrogen atoms of the central pyrazine ring. Although the magnitude of the second harmonic signal is similar, HRS measurements revealed that the anisotropy of the NLO response strongly differs in the two chromophore series, the one of the 2,3-isomers being strikingly dipolar, while the one of the 2,6-isomers is mostly octupolar. The experimental observations are well supported by DFT calculations. In particular, the sum-over-states approach allows us to rationalize the different NLO anisotropies observed in the two isomer series through a detailed analysis of the symmetry of the low-lying excited states.

Dynamic effects on the nonlinear optical properties of donor acceptor stenhouse adducts: insights from combined MD + QM simulations.

The second-order nonlinear optical (NLO) responses of a donor-acceptor stenhouse adduct (DASA) are investigated by using a computational approach combining molecular dynamics simulations and density functional theory (DFT) calculations. Specific force fields for the open and closed photoswitching forms are first parameterized and validated according to the Joyce protocol, in order to finely reproduce the geometrical features and potential energy surfaces of both isomers in chloroform solution. Then, DFT calculations are performed on structural snapshots extracted at regular time steps of the MD trajectories to address the influence of the thermalized conformational dynamics on the NLO responses related to hyper-Rayleigh scattering (HRS) experiments. We show that accounting for the structural dynamics largely enhances the HRS hyperpolarizability (ßHRS) compared to DFT calculations considering solely equilibrium geometries, and greatly improves the agreement with experimental measurements. Furthermore, we show that the NLO responses of the NLO-active open form are correlated with the bond order alternation along the triene bridge connecting the donor and acceptor moieties, which is rationalized using simple essential state models.

Nonlinear Optical Responses of Photoswitchable Donor-Acceptor Stenhouse Adducts.

This work combines hyper-Rayleigh scattering (HRS) experiments performed in the NIR range (1.30 and 1.60 µm) and quantum chemical calculations to provide a comprehensive description of the second harmonic generation (SHG) responses of donor-acceptor Stenhouse adducts (DASAs). Representative derivatives of the three generations of DASAs, which differ by the nature of their electron-donating and withdrawing moieties and also include clickable species, have been synthesized and their photoswitching behavior fully characterized. The HRS measurements allow us to establish relationships between the magnitude of the SHG response of open forms and the nature of the donor and acceptor groups. The largest SHG responses are obtained for derivatives incorporating either a barbituric acid or an indanedione acceptor unit, while N-methylaniline appears as the most efficient donor group. The calculations support well the experimental data and show that high hyperpolarizabilities are associated to low excitation energies and large extent of the photoinduced intramolecular charge transfer, which enhances the dipole moment variation between the ground and first dipole-allowed electronic excited state. In addition, a complete investigation of the photoswitching kinetics of DASAs in chloroform solution shows important differences, highlighting in particular the role of the donor group on the photoswitching efficiency.

Predicting the Second-Order Nonlinear Optical Responses of Organic Materials: The Role of Dynamics.

The last 30 years have witnessed an ever-growing application of computational chemistry for rationalizing the nonlinear optical (NLO) responses of organic chromophores. More specifically, quantum chemical calculations proved highly helpful in gaining fundamental insights into the factors governing the magnitude and character of molecular first hyperpolarizabilities (ß), be they either intrinsic to the chromophore molecular structure and arising from symmetry, chemical substitution, or π-electron delocalization, or induced by external contributions such as the laser probe or solvation and polarization effects. Most theoretical reports assumed a rigid picture of the investigated systems, the NLO responses being computed solely at the most stable geometry of the chromophores. Yet, recent developments combining classical molecular dynamics (MD) simulations and DFT calculations have evidenced the significant role of structural fluctuations, which may induce broad distributions of NLO responses, and even generate them in some instances.This Account presents recent case studies in which theoretical simulations have highlighted these effects. The discussion specifically focuses on the simulation of the second-order NLO properties that can be measured experimentally either from Hyper-Rayleigh Scattering (HRS) or Electric-Field Induced Second Harmonic Generation (EFISHG). More general but technical topics concerning several aspects of the calculations of hyperpolarizabilities are instead discussed in the Supporting Information.Selected examples include organic chromophores, photochromic systems, and ionic complexes in the liquid phase, for which the effects of explicit solvation, concentration, and chromophore aggregation are emphasized, as well as large flexible systems such as peptide chains and pyrimidine-based helical polymers, in which the relative variations of the responses were shown to be several times larger than their average values. The impact of geometrical fluctuations is also illustrated for supramolecular architectures with the examples of nanoparticles formed by organic dipolar dyes in water solution, whose soft nature allows for large shape variations translating into huge fluctuations in time of their NLO response, and of self-assembled monolayers (SAMs) based on indolino-oxazolidine or azobenzene switches, in which the geometrical distortions of the photochromic molecules, as well as their orientational and positional disorder within the SAMs, highly impact their NLO response and contrast upon switching. Finally, the effects of the rigidity and fluidity of the surrounding are evidenced for NLO dyes inserted in phospholipid bilayers.

Electric-field induced second harmonic generation responses of push-pull polyenic dyes: experimental and theoretical characterizations.

The second-order nonlinear optical properties of four series of amphiphilic cationic chromophores involving different push-pull extremities and increasingly large polyenic bridges have been investigated both experimentally, by means of electric field induced second harmonic (EFISH) generation, and theoretically, using a computational approach combining classical molecular dynamics (MD) and quantum chemical (QM) calculations. This theoretical methodology allows to describe the effects of structural fluctuations on the EFISH properties of the complexes formed by the dye and its iodine counterion, and provides a rationale to EFISH measurements. The good agreement between experimental and theoretical results proves that this MD + QM scheme constitutes a useful tool for a rational, computer-aided, design of SHG dyes.

A quantum chemical investigation of the second hyperpolarizability of p-nitroaniline.

Recent measurements of the third harmonic scattering responses of molecules have given a new impetus for computing molecular second hyperpolarizabilities (γ) and for deducing structure-property relationships. This paper has employed a variety of wavefunction and density functional theory methods to evaluate the second hyperpolarizability of the p-nitroaniline prototypical push-pull π-conjugated molecule, addressing also numerical aspects, such as the selection of an integration grid and the impact of the order of differentiation vs the achievable accuracy by using the Romberg quadrature. The reliability of the different methods has been assessed by comparison to reference Coupled-Cluster Singles and Doubles with perturbative treatment of the Triples results. On the one hand, among wavefunction methods, the MP2 scheme offers the best accuracy/cost ratio for computing the static γ. On the other hand, using density functional theory, γ remains a challenging property to compute because all conventional, global hybrid or range-separated hybrid, exchange-correlation functionals underestimate static γ values by at least 15%. Even tuning the range-separating parameter to minimize the delocalization errors does not enable to improve the γ values. Nevertheless, the original double-hybrid B2-PLYP functional, which benefits from 27% of PT2 correlation and 53% Hartree-Fock exchange, provides accurate estimates of static γ values. Unfortunately, the best performing exchange-correlation functionals for γ are not necessarily reliable for the first hyperpolarizability, ß, and vice versa. In fact, the ß of p-nitroaniline (pNA) could be predicted, with a good accuracy, with several hybrid exchange-correlation functionals (including by tuning the range-separating parameter), but these systematically underestimate γ. As for γ, the MP2 wavefunction method remains the best compromise to evaluate the first hyperpolarizability of pNA at low computational cost.

First hyperpolarizability of the di-8-ANEPPS and DR1 nonlinear optical chromophores in solution. An experimental and multi-scale theoretical chemistry study.

The solvent effects on the linear and second-order nonlinear optical properties of an aminonaphtylethenylpyridinium (ANEP) dye are investigated by combining experimental and theoretical chemistry methods. On the one hand, deep near infrared (NIR) hyper-Rayleigh scattering (HRS) measurements (1840-1950 nm) are performed on solutions of di-8-ANEPPS in deuterated chloroform, dimethylformamide, and dimethylsulfoxide to determine their first hyperpolarizablity (ßHRS). For the first time, these HRS experiments are carried out in the picosecond regime in the deep NIR with very moderate (≤3 mW) average input power, providing a good signal-to-noise ratio and avoiding solvent thermal effects. Moreover, the frequency dispersion of ßHRS is investigated for Disperse Red 1 (DR1), a dye commonly used as HRS external reference. On the other hand, these are compared with computational chemistry results obtained by using a sequential molecular dynamics (MD) then quantum mechanics (QM) approach. The MD method allows accounting for the dynamical nature of the molecular structures. Then, the QM part is based on TDDFT/M06-2X/6-311+G* calculations using solvation models ranging from continuum to discrete ones. Measurements report a decrease of the ßHRS of di-8-ANEPPS in more polar solvents and these effects are reproduced by the different solvation models. For di-8-ANEPPS and DR1, comparisons show that the use of a hybrid solvation model, combining the description of the solvent molecules around the probe by point charges with a continuum model, already achieves quasi quantitative agreement with experiment. These results are further improved by using a polarizable embedding that includes the atomic polarizabilities in the solvent description.

Asymmetric Allenylation of Alkynes mediated by Chiral Organoselenated Reagents under Oxidative Conditions.

The reactivity of novel chiral lactamide-substituted diselenide-based reagents under oxidative conditions was exploited to develop a metal-free method for the preparation of enantioenriched allenylamides from simple alkynes in good yields, and with enantiomeric excesses up to 99 %. The key of the success in this method is attributed to the hydrogen-bonded lactamide appendages that ensure configurational stability of chiral vinyl selenoxide intermediates for an optimal enantiotopic ß-syn-elimination step.

Multi-State Second-Order Nonlinear Optical Switches Incorporating One to Three Benzazolo-Oxazolidine Units: A Quantum Chemistry Investigation.

This contribution employs quantum chemistry methods to describe the variations of the second nonlinear optical responses of molecular switches based on benzazolo-oxazolidine (BOX) units, connected by π-linkers, along their successive opening/closing. Under the fully closed forms, all of them display negligible first hyperpolarizability (ß) values. When one BOX is opened, which is sketched as C→O, a push-pull π-conjugated segment is formed, having the potential to enhance ß and to set the depolarization ratio (DR) to its one-dimensional-like value (DR = 5). This is observed when only one BOX is open, either for the monoBOX species (C→O) or for the diBOX (CC→CO) and triBOX (CCC→CCO) compounds, i.e., when the remaining BOXs stay closed. The next BOX openings have much different effects. For the diBOXs, the second opening (CO→OO) is associated with a decrease of ß, and this decrease is tuned by controlling the conformation of the π-linker, i.e., the centrosymmetry of the whole compound because ß vanishes in centrosymmetric compounds. For the triBOXs, the second opening gives rise to a Λ-shape compound, with a negligible change of ß, but a decrease of the DR whereas, along the third opening, ß remains similar and the DR decreases to the typical value of octupolar systems (DR = 1.5).

Stealth Luminescent Organic Nanoparticles Made from Quadrupolar Dyes for Two-Photon Bioimaging: Effect of End-Groups and Core.

Molecular-based Fluorescent Organic Nanoparticles (FONs) are versatile light-emitting nano-tools whose properties can be rationally addressed by bottom-up molecular engineering. A challenging property to gain control over is the interaction of the FONs' surface with biological systems. Indeed, most types of nanoparticles tend to interact with biological membranes. To address this limitation, we recently reported on two-photon (2P) absorbing, red to near infrared (NIR) emitting quadrupolar extended dyes built from a benzothiadiazole core and diphenylamino endgroups that yield spontaneously stealth FONs. In this paper, we expand our understanding of the structure-property relationship between the dye structure and the FONs 2P absorption response, fluorescence and stealthiness by characterizing a dye-related series of FONs. We observe that increasing the strength of the donor end-groups or of the core acceptor in the quadrupolar (D-π-A-π-D) dye structure allows for the tuning of optical properties, notably red-shifting both the emission (from red to NIR) and 2P absorption spectra while inducing a decrease in their fluorescence quantum yield. Thanks to their strong 1P and 2P absorption, all FONs whose median size varies between 11 and 28 nm exhibit giant 1P (106 M-1.cm-1) and 2P (104 GM) brightness values. Interestingly, all FONs were found to be non-toxic, exhibit stealth behaviour, and show vanishing non-specific interactions with cell membranes. We postulate that the strong hydrophobic character and the rigidity of the FONs building blocks are crucial to controlling the stealth nano-bio interface.

Impact of van der Waals interactions on the structural and nonlinear optical properties of azobenzene switches.

The geometrical structures, relative Z-E energies, and second-order nonlinear responses of a collection of azobenzene molecules symmetrically substituted in the meta-position with functional groups of different bulkiness are investigated using various ab initio and density functional approximations. We show that RI-MP2 and RI-CC2 approximations provide very similar geometries and relative energies and evidence that London dispersion interactions existing between bulky meta-substituents stabilize the Z conformer. The ωB97X-D exchange-correlation functional provides an accurate description of these effects and gives a good account of the nonlinear optical response of the molecules. We show that density functional approximations should include no less than 50% of Hartree-Fock exchange to provide accurate hyperpolarizabilities. A property-structure analysis of the azobenzene derivatives reveals that the main contribution to the first hyperpolarizability comes from the azo bond, but phenyl meso-substituents can enhance it.

Self-assembling, structure and nonlinear optical properties of fluorescent organic nanoparticles in water.

Owing to their intense emission, low toxicity and solubility in aqueous medium, fluorescent organic nanoparticles (FONs) have emerged as promising alternatives to inorganic ones for the realization of exogenous probes for bioimaging applications. However, the intimate structure of FONs in solution, as well as the role played by intermolecular interactions on their optical properties, remains challenging to study. Following a recent Second-Harmonic Scattering (SHS) investigation led by two of us [Daniel et al., ACS Photonics, 2015, 2, 1209], we report herein a computational study of the structural organization and second-order nonlinear optical (NLO) properties of FONs based on dipolar chromophores incorporating a hydrophobic triphenylamine electron-donating unit and a slightly hydrophilic aldehyde electron-withdrawing unit at their extremities. Molecular dynamics simulations of the FON formation in water are associated with quantum chemical calculations, to provide insight into the molecular aggregation process, the molecular orientation of the dipolar dyes within the nanoparticles, and the dynamical behavior of their NLO properties. Moreover, the impact of intermolecular interactions on the NLO responses of the FONs is investigated by employing the tight-binding version of the recently developed simplified time-dependent density functional theory (sTD-DFT) approach, allowing the all-atom quantum mechanics treatment of nanoparticles.

Bispericyclic Diels-Alder Dimerization of ortho-Quinols in Natural Product (Bio)Synthesis: Bioinspired Chemical 6-Step Synthesis of (+)-Maytenone.

Many natural products of plant or microbial origins are derived from enzymatic dearomative oxygenation of 2-alkylphenolic precursors into 6-alkyl-6-hydroxycyclohexa-2,4-dienones. These so-called ortho-quinols cyclodimerize via a remarkably selective bispericyclic Diels-Alder reaction. Whether or not the intervention of catalytic or dirigent proteins is involved during this final step of the biosynthesis of these natural products, this cyclodimerization of ortho-quinols can be chemically reproduced in the laboratory with the same strict level of site-specific regioselectivity and stereoselectivity. This unique yet unified process, which finds its rationale in the inherent chemical reactivity of those ortho-quinols, is illustrated herein by an efficient and bioinspired first chemical synthesis of one of the most structurally complex and synthetically challenging examples of such natural cyclodimers, the bisditerpenoid (+)-maytenone.

Photophysics, Electrochemistry and Efficient Electrochemiluminescence of Trigonal Truxene-Core Dyes.

We synthesized and characterized a series of dyes built from a spirofluorene or truxene core. The quadrupolar spirofluorene system is the initial building unit for the design and preparation of more complex star-shaped dyes consisting of a truxene core bearing three di- or triphenylamine moieties with or without a thiophene connector. Their photophysical, electrochemical, and electrochemiluminescence (ECL) properties were first investigated in solution. Structure/activity relationships were derived and rationalized by comparing the quadrupolar system and trigonal truxene-core derivatives using computational studies. The photophysical and redox characteristics are drastically tuned by the introduction of a thiophene bridge and electron-donor substituents at their terminal branches. These comparative studies show the essential role of the stability of both radical cations and anions to obtain efficient ECL dyes. The stabilization of the radicals is directly related to the charge delocalization due to the π-conjugation by the thiophene bridge. The brightest ECL is achieved by annihilation and coreactant (benzoyl peroxide) pathways with the blue-emitting truxene dye, which is 2- and 4.5-times greater than that of the quadrupolar compound and reference [Ru(bpy)3 ]2+ emitter, respectively. Such an extensive study on these extended π-conjugated molecules presenting different core structures may guide the design and synthesis of new ECL dyes with a strong efficiency.

Exciton Interactions, Excimer Formation, and [2π+2π] Photodimerization in Nonconjugated Curcuminoid-BF2 Dimers.

We describe the synthesis of a series of covalently linked dimers of quadrupolar curcuminoid-BF2 dyes and the detailed investigation of their solvent-dependent spectroscopic and photophysical properties. In solvents of low polarity, intramolecular folding induces the formation of aggregated chromophores, the UV/Vis absorption spectra of which display the optical signature characteristic of weakly-coupled H-aggregates. The extent of folding and, in turn, of ground-state aggregation is strongly dependent on the nature of the flexible linker. Steady-state and time-resolved fluorescence emission spectroscopies show that the Frenkel exciton relaxes into a fluorescent symmetrical excimer state with a long lifetime. Furthermore, our in-depth studies show that a weakly emitting excimer lies on the pathway toward a photocyclomer. Two-dimensional 1 H NMR spectroscopy and density functional theory (DFT) allowed the structure of the photoproduct to be established. To our knowledge, this represents the first example of a [2π+2π] photodimerization of the curcuminoid chromophore.

Flavylium Fluorophores as Near-Infrared Emitters.

This theoretical work rationalizes the absorption and fluorescence emission properties of conjugated dyes composed of dimethylamino flavylium heterocycles linked by a polymethine chain, which were recently reported to act as efficient shortwave infrared emitters. Density functional theory is used to characterize the electronic structure of the low-lying excited states as a function of the polymethine chain length. Decomposition of the computed excitations in terms of diabatic states is also performed to deconvolute the excited states wavefunction into charge-transfer intramolecular excitations. Based on these results, chemical substitution patterns consisting in enhancing the electron-withdrawing strength of the polymethine bridge and the electron-donating ability of the lateral flavylium fragments, are proposed to further redshift the photoluminescence of the fluorophores.

Performance of DFT functionals for calculating the second-order nonlinear optical properties of dipolar merocyanines.

The second-order nonlinear optical responses of a series of recently designed dipolar merocyanines are investigated using the 2006 Minnesota family of hybrid exchange-correlation functionals (XCFs), as well as the LC-BLYP, ωB97XD and CAM-B3LYP long-range (LR) corrected XCFs. The performance of these different levels of approximation is discussed in regard to reference second-order Møller-Plesset calculations and experimental data obtained from Hyper-Rayleigh Scattering (HRS) measurements. Particular focus is given to the influence of the amount of exact Hartree-Fock exchange included in the XCF on the magnitude of the static HRS responses, as well as to the impact of tuning the range-separation parameter in LR-XCFs, according to a system-specific nonempirical procedure. Frequency dispersion effects are also investigated, as well as their crucial role in the comparison between theoretical and experimental data.

Nonlinear optical properties of spirocyclohexadine photochromes: insights from DFT calculations.

The second-order nonlinear optical properties of 1,3-indandione-derived spirocyclohexadine compounds, a new class of photochromes showing sensitivity to both UV and visible lights, are investigated by means of density functional theory. The calculations predict that these compounds should display a large first hyperpolarizability variation upon commutation between their cyclic spiropyran and their merocyanine form. We also show that significant enhancement of the NLO response of the merocyanine form, as well as of the NLO contrast, can be achieved by means of appropriate chemical substitution, which makes these derivatives highly appealing for potential use in photo-responsive NLO devices.