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.

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.

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.

Insight into the generation of near infra-red (NIR) absorbing species in electrochromic surface-anchored metal-organic frameworks.

Due to their versatility and easy processing, Surface-Anchored Metal-Organic Frameworks (SurMOFs) have gained interest in recent times as promising electrochromic thin films. Herein a step forward in their use and characterization was achieved thanks to the integration of {Zn2(PDICl4)2} SurMOFs in a multi-layer electrochromic device (ECD), based on a membrane-like electrolyte. The optical and electrochemical properties of the ECD were fully characterized, revealing a two-step reduction process localized on the organic ligand and involving subsequent near infra-red (NIR) and cyan absorbing states, leading to optical modulation of the films. The species responsible for this absorption were isolated and identified in the reduced states. In parallel to experimental characterization, quantum chemistry was successfully used to investigate the structure-property relationship of the SurMOF, revealing additional information regarding the structure and the local environment of the electrochromic ligand.

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.

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.

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.

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.

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.

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.

Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties?

Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for ß and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.

Influence of the Nature of the Lipid Building Blocks on the Second-Order Nonlinear Optical Responses of an Embedded Di-8-ANEPPS Probe.

The extensive collection of lipids found in cell membranes is justified by the fact that each lipid contributes to their overall structure, dynamics, and properties and so to the biological processes taking place within them. It also showcases that, in order to deepen our understanding of membranes, we need to have a tool to differentiate lipid bilayers of varying composition. In this work, we investigate a suite of single-component saturated glycerophospholipids varying only in their headgroup structure by analyzing the second harmonic generation (SHG) nonlinear optical (NLO) response of a probe, di-8-ANEPPS, embedded into the membranes. The seven hydrophilic heads chosen (phosphatidylcholine (PC), phosphatidylethanolamine (PE), diaglycerol (GL), phosphatidylserine (PS), phosphatidylglycerol (PG), phosphatidylinositol (PI), and phosphatidyc acid (PA)) represent all the major headgroups that are part of mammalian plasma membranes and provide an assortment of neutral, zwiterrionic, and charged species. First, molecular dynamics simulations revealed that the lipidic arrangement is strongly sensitive to the nature of the hydrophilic head and less to the variety in the hydrophobic region. Membranes exhibiting drastically opposite structural properties can be pointed out: 1,2-dihexadecanoyl-rac-glycerol (DPGL) is the thickest and most ordered and aligned system, whereas 1,2-diacyl-sn-glycero-3-phospho-(1'-sn-glycerol) (DPPG) is thinnest and least ordered and aligned system. The structural analyses are then confronted with the molecular NLO responses, ß, computed at the time-dependent density functional theory (TDDFT) level. As the orientation of the chromophore is impacted by the various degrees of order within the lipid bilayers, the diagonal component of the ß tensor parallel to the bilayer normal, ßZZZ, is as well. In the end, this computational approach provides insights into the link between lipid building blocks and the NLO responses of the embedded dye.

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.

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.

Temperature-Dependent Structural Phase Transition in Rubrene Single Crystals: The Missing Piece from the Charge Mobility Puzzle?

Accurate structural models for rubrene, the benchmark organic semiconductor, derived from synchrotron X-ray data in the temperature range of 100-300 K, show that its cofacially stacked tetracene backbone units remain blocked with respect to each other upon cooling to 200 K and start to slip below that temperature. The release of the blocked slippage occurs at approximately the same temperature as the hole mobility crossover. The blocking between 200 and 300 K is caused by a negative correlation between the relatively small thermal expansion along the crystallographic b-axis and the relatively large widening of the angle between herringbone-stacked tetracene units. DFT calculations reveal that this blocked slippage is accompanied by a discontinuity in the variation with temperature of the electronic couplings associated with hole transport between cofacially stacked tetracene backbones.

Direct Triple Annulations: A Way to Design Large Triazastarphenes with Intertwined Hexagonal Packing.

A new straightforward synthetic strategy has been elaborated to achieve star-shaped triazatrinaphthylene and, for the first time, triazatrianthrylene derivatives. Their solution- and solid-state properties were thoroughly characterized by cyclic voltammetry, UV-vis absorption spectroscopy, X-ray diffraction, and density functional theory calculations. Original hexagonal molecular arrangements were found in the crystal phase, which opens a new pathway for designing materials with improved three-dimensional charge-transport properties.

Spontaneous exciton dissociation enables spin state interconversion in delayed fluorescence organic semiconductors.

Engineering a low singlet-triplet energy gap (ΔEST) is necessary for efficient reverse intersystem crossing (rISC) in delayed fluorescence (DF) organic semiconductors but results in a small radiative rate that limits performance in LEDs. Here, we study a model DF material, BF2, that exhibits a strong optical absorption (absorption coefficient = 3.8 × 105 cm-1) and a relatively large ΔEST of 0.2 eV. In isolated BF2 molecules, intramolecular rISC is slow (delayed lifetime = 260 µs), but in aggregated films, BF2 generates intermolecular charge transfer (inter-CT) states on picosecond timescales. In contrast to the microsecond intramolecular rISC that is promoted by spin-orbit interactions in most isolated DF molecules, photoluminescence-detected magnetic resonance shows that these inter-CT states undergo rISC mediated by hyperfine interactions on a ~24 ns timescale and have an average electron-hole separation of ≥1.5 nm. Transfer back to the emissive singlet exciton then enables efficient DF and LED operation. Thus, access to these inter-CT states, which is possible even at low BF2 doping concentrations of 4 wt%, resolves the conflicting requirements of fast radiative emission and low ΔEST in organic DF emitters.

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.

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.