Evaluation of the nonlinear optical properties for an expanded porphyrin Hückel-Möbius aromaticity switch

Application of Inverse Design Approaches to the Discovery of Nonlinear Optical Switches

Molecular switches, in which a stimulus induces a large and reversible change in molecular properties, are of significant interest in the domain of photonics. Due to their commutable redox states with distinct nonlinear optical (NLO) properties, hexaphyrins have emerged as a novel platform for multistate switches in nanoelectronics. In this study, we employ an inverse design algorithm to find functionalized 26R→28R redox switches with maximal βHRS contrast. We focus on the role of core modifications, since a synergistic effect with meso-substitutions was recently found for the 30R-based switch. In contrast to these findings, the inverse design optima and subsequent database analysis of 26R-based switches confirm that core modifications are generally not favored when high NLO contrasts are targeted. Moreover, while push-pull combinations enhance the NLO contrast for both redox switches, they prefer a different arrangement in terms of electron-donating and electron-withdrawing functional groups. Finally, we aim at designing a three-state 26R→28R→ 30R switch with a similar NLO response for both ON states. Even though our best-performing three-state switch follows the design rules of the 30R-based component, our chemical compound space plots show that well-performing three-state switches can be found in regions shared by high-responsive 26R and 30R structures.

Designing hexaphyrins for high-potential NLO switches: the synergy of core-modifications and meso-substitutions

Due to the enormous size of the chemical compound space, usually only small regions are traversed with traditional direct molecular design approaches making the discovery for novel functionalized molecules for nonlinear optical applications challenging. By applying inverse molecular design algorithms, we aim to efficiently explore larger regions of the compound space in search of promising hexaphyrin-based molecular switches as measured by their first-hyperpolarizability (β_HRS) contrast. We focus on the 28R → 30R switch with a functionalization pattern allowing for centrosymmetric OFF states yielding zero β_HRS response. This switch is particularly challenging as full meso-substitution with a single type of functional group or core-modifications result in almost no contrast enhancement. We carried out four inverse design procedures during which two sets of core-modifications and three sets of meso-substitutions sites were systematically optimized. All 4 optimal switches are characterized by a mix of meso-substitutions and core-modifications, of which the best performing switch yields a 10-fold improvement over the parent macrocycle. Throughout the inverse design procedures, we collected and analyzed a database biased towards high NLO contrasts that contains 277 different patterns for hexaphyrin-based switches. We derived three design rules to obtain highly functional 28R → 30R NLO switches: (I) a combination of 2 strong EWG and 1 EDG group is the ideal recipe for increasing the NLO contrast, though their position also plays an important role. (II) The type of core-modification is less important when only the diagonal positions are core-modified. Switches with 4 core-modifications show a clear preference for oxygen. (III) Keeping centrosymmetry in the OFF state remains highly beneficial given the investigated functionalization pattern. Finally, we have demonstrated that combining meso-substitutions with core-modifications can synergistically improve the NLO contrast.

Effects of external field wavelength and solvation on the photophysical property and optical nonlinearity of 1,3-thiazolium-5-thiolates mesoionic compound

The photophysical property and optical nonlinearity of an electronic push-pull mesoionic compound, 2-(4-trifluoromethophenyl)-3-methyl-4-(4-methoxyphenyl)-1,3-thiazole-5-thiolate, were theoretically investigated with a reliable computing strategy. The essence of the optical properties were then explored through a variety of wavefunction analysis methods, including the natural transition orbital analysis, hole-electron analysis, (hyper)polarizability density analysis, decomposition of the (hyper)polarizability contribution, and (hyper)polarizability tensor analysis, at the level of electronic structure. The influences of the electric field and solvation on the absorption spectrum and (hyper)polarizability of the molecule are highlighted and clarified. The results are expected to provide guidance for people to understand the effects of external field wavelength and solvation on the optical properties of mesoscopic molecules.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Second-Order Nonlinear Optics Response of the Boron-Dipyrromethenes-Based Mislinked Expanded Porphyrins: Revealing the Role of the -BF2 Group

Here, the mislinked expanded porphyrins singly (labeled A) and doubly (labeled B) neo-confused [22]smaragdyrin, the boron-dipyrromethenes-based mislinked expanded porphyrins singly (labeled C) and doubly (labeled D) neo-confused [22]smaragdyrin, where both C and D include a -BF₂ group, are chosen to serve as the study objects, and theoretical calculations are carried out to study the role of the -BF₂ group in the second-order nonlinear optics (NLO) behaviors. Results highlighted that the -BF₂ group plays an important role for the second-order behaviors in mislinked expanded porphyrins; namely, embedding the -BF₂ group well enhanced the hyper-Rayleigh scattering (HRS) value {β_HRS(0;0,0)}, C{βHRS(0;0,0)}A{βHRS(0;0,0)} = 2.0 and D{βHRS(0;0,0)}B{βHRS(0;0,0)} = 2.9, main owning to the fact that installing -BF₂ increases the electron delocalization degree and decreases the excited energy of the crucial excited state.

Fine-Tuning of Nonlinear Optical Contrasts of Hexaphyrin-Based Molecular Switches Using Inverse Design

In the search for new nonlinear optical (NLO) switching devices, expanded porphyrins have emerged as ideal candidates thanks to their tunable chemical and photophysical properties. Introducing meso-substituents to these macrocycles is a successful strategy to enhance the NLO contrasts. Despite its potential, the influence of meso-substitution on their structural and geometrical properties has been scarcely investigated. In this work, we pursue to grasp the underlying pivotal concepts for the fine-tuning of the NLO contrasts of hexaphyrin-based molecular switches, with a particular focus on the first hyperpolarizability related to the hyper-Rayleigh scattering (β_HRS). Building further on these concepts, we also aim to develop a rational design protocol. Starting from the (un)substituted hexaphyrins with various π-conjugation topologies and redox states, structure-property relationships are established linking aromaticity, photophysical properties and β_HRS responses. Ultimately, inverse molecular design using the best-first search algorithm is applied on the most favorable switches with the aim to further explore the combinatorial chemical compound space of meso-substituted hexaphyrins in search of high-contrast NLO switches. Two definitions of the figure-of-merit of the switch performance were used as target objectives in the optimization problem. Several meso-substitution patterns and their underlying characteristics are identified, uncovering molecular symmetry and the electronic nature of the substituents as the key players for fine-tuning the β_HRS values and NLO contrasts of hexaphyrin-based switches.

Theoretical Framework of 1,3-Thiazolium-5-Thiolates Mesoionic Compounds: Exploring the Nature of Photophysical Property and Molecular Nonlinearity

Inspired by a previous experimental study on the first-order hyperpolarizabilities of 1,3-thiazolium-5-thiolates mesoionic compounds using the hyper-Rayleigh scattering technique, we theoretically investigated the UV-vis absorption spectra and every-order polarizabilities of these mesoionic molecules. Based on the fact that the photophysical and nonlinear properties observed in the experiment can be perfectly replicated, our theoretical calculations explored the essential characteristics of the optical properties of the mesoionic compounds with different electron-donating groups at the level of electronic structures through various wave function analysis methods. The influence of the electron-donating ability of the donor on the optical properties of the molecules and the contribution of the mesoionic ring moiety to their optical nonlinearity are clarified, which have not been reported by any research so far. This work will help people understand the nature of optical properties of mesoionic-based molecules and provide guidance for the rational design of molecules with excellent photoelectric performance in the future.

Chalcone-based molecules: Experimental and theoretical studies on the two-photon absorption and molecular first hyperpolarizability

Five chalcone-based molecules denominated by C-3 ((E)-1-(4-methoxyphenyl)-3-phenylprop-2-en-1-one), C-4 ((E)-1,3-bis(4-methoxyphenyl)prop-2-en-1-one), C-5 ((E)-1-(benzo[d][1,3]dioxol-5-yl)-3-(4-methoxyphenyl)prop-2-en-1-one), C-6 ((E)-3-(naphthalen-1-yl)-1-phenylprop-2-en-1-one) and C-7 ((E)-1-(4-methoxyphenyl)-3-(naphthalen-1-yl)prop-2-en-1-one) were synthesized by Claisen-Schmidt reaction in solution of NaOH in water/ethanol 2:1. The aldehydes used were benzaldehyde, anisaldehyde, and β-naphthaldehyde, while the used ketones were acetophenone, p-methoxyacetophenone, and 3,4-methylenedioxyacetophenone. Z-scan and hyper-Rayleigh scattering techniques were used to study the nonlinear optical properties of these compounds in dichloromethane medium. By using Z-scan technique with femtosecond pulses, two-photon absorption cross-sections (σ^TPA) were determined, while the first molecular electronic hyperpolarizabilities (β^HRS) were evaluated by the hyper-Rayleigh scattering technique, with picosecond pulses. From the recorded two-photon absorption spectra, it was identified that compound C-7 presented the highest σ^TPA, regarding the HOMO-LUMO transition, with a value of 40 GM, while C-6 achieved the lowest value for the same transition with 13 GM. Concerning the values of the first molecular hyperpolarizability, compound C-4 presented the highest value, 38 × 10^-30 cm⁴ statvolt^-1, while C-3 presented the lowest β^HRS value of about 16 × 10^-30 cm⁴ statvolt^-1. Time-dependent density functional theory calculations were used to simulate the one- and two-photon absorption spectra, as well to predict the theoretical value of β^HRS in dichloromethane and vacuum medium.

Aromaticity as a Guiding Concept for Spectroscopic Features and Nonlinear Optical Properties of Porphyrinoids

With their versatile molecular topology and aromaticity, porphyrinoid systems combine remarkable chemistry with interesting photophysical properties and nonlinear optical properties. Hence, the field of application of porphyrinoids is very broad ranging from near-infrared dyes to opto-electronic materials. From previous experimental studies, aromaticity emerges as an important concept in determining the photophysical properties and two-photon absorption cross sections of porphyrinoids. Despite a considerable number of studies on porphyrinoids, few investigate the relationship between aromaticity, UV/vis absorption spectra and nonlinear properties. To assess such structure-property relationships, we performed a computational study focusing on a series of Hückel porphyrinoids to: (i) assess their (anti)aromatic character; (ii) determine the fingerprints of aromaticity on the UV/vis spectra; (iii) evaluate the role of aromaticity on the NLO properties. Using an extensive set of aromaticity descriptors based on energetic, magnetic, structural, reactivity and electronic criteria, the aromaticity of [4n+2] π-electron porphyrinoids was evidenced as was the antiaromaticity for [4n] π-electron systems. In agreement with previous studies, the absorption spectra of aromatic systems display more intense B and Q bands in comparison to their antiaromatic homologues. The nature of these absorption bands was analyzed in detail in terms of polarization, intensity, splitting and composition. Finally, quantities such as the average polarizability and its anisotropy were found to be larger in aromatic systems, whereas first and second hyperpolarizability are influenced by the interplay between aromaticity, planarity and molecular symmetry. To conclude, aromaticity dictates the photophysical properties in porphyrinoids, whereas it is not the only factor determining the magnitude of NLO properties.

New electron delocalization tools to describe the aromaticity in porphyrinoids

There are several possible pathways in the macrocycle of large porphyrinoids and, among aromaticity indices, only AV_minis capable of recognizing the most aromatic one.

The novel link between planar möbius aromatic and third order nonlinear optical properties of metal-bridged polycyclic complexes

Metal–bridged polcyclic aromatic complexes, exhibiting unusual optical effects such as near-infrared photoluminescence with particularly large Stokes shifts, long lifetimes and aggregation enhancement, have been established as unique “carbonloong chemistry”. Herein, the electronic structures, aromaticities, absorption spectra and third order nonlinear optical (NLO) responses of metal–bridged polcyclic aromatic complexes (M = Fe, Re, Os and Ir) are investigated using the density functional theory computations. It is found that the bridge–head metal can stabilize and influence rings, thus creating π–, σ– and metalla–aromaticity in an extended, π–conjugated framework. Interestingly, metal radius greatly influence the bond, aromaticity, liner and third order NLO properties, which reveals useful information to develop new applications of metal regulatory mechanism in NLO materials field. Significantly, the novel relationship between the aromaticity and third order NLO response has firstly been proposed, that the metal-bridged polycyclic complex with larger aromaticity will exhibit larger third order nonlinear optical response. It is our expectation that the novel link between aromaticity and NLO response could provide valuable information for scientists to develop the potential NLO materials on the basis of metal–bridged polycyclic complexes.

Computational Studies on Optoelectronic and Nonlinear Properties of Octaphyrin Derivatives

The electronic and nonlinear optical (NLO) properties of octaphyrin derivatives were studied by employing the DFT/TDFT at CAM-B3LYP/6-311++G (2d, 2p) level of the theory. Thiophene, phenyl, methyl and cyano moieties were substituted on the molecular framework of octaphyrin core, in order to observe the change in optoelectronic and nonlinear response of these systems. The frontier molecular orbital studies and values of electron affinity reveals that the studied compounds are stable against the oxygen and moisture present in air. The calculated ionization energies, adiabatic electron affinity and reorganization energy values indicate that octaphyrin derivatives can be employed as effective n-type material for Organic Light Emitting Diodes (OLEDs). This character shows an enhancement with the introduction of an electron withdrawing group in the octaphyrin framework. The polarizability and hyperpolarizability values of octaphyrin derivatives demonstrate that they are good candidates for NLO devices. The nonlinear response of these systems shows enhancement on the introduction of electron donating groups on octaphyrin moiety. However, these claims needs further experimental verification.

Characterization of the binding of six actinyls AnO22+/+(An = U/Np/Pu) with three expanded porphyrins by density functional theory

Oxasapphyrin has stronger binding affinity than amethyrin and grandephyrin in complexation with six hydrated actinyl cations, AnO₂(H₂O)₅^2+/+(An = U/Np/Pu).

Understanding conductivity in molecular switches: a real space approach in octaphyrins

In recent years, expanded porphyrins have emerged as a promising class of π-conjugated switches whose conductance is studied from the electron density.

Understanding the molecular switching properties of octaphyrins

Triggering Hückel–Möbius topological and aromaticity switches in octaphyrins by protonation and redox reactions.

Exploring the chemical bonding, infrared and UV-vis absorption spectra of OH radicals adsorption on the smallest fullerene

In the present work, the density-functional theory calculations were performed on C20 hydroxylated fullerene. B3LYP functionals with 6-31G(d,p) basis set were utilized to gain insight into the bonding characters and intramolecular interactions of hydroxyl groups adsorbed on the cage. Interestingly, we observed that the C20 cage has the bonding patterns with spherical orbitals configuration [1S(2)1P(6)1D(10)1F(2)], and the adsorbed hydroxyl groups significantly affect the chemical bonding of the cage surface. Analysis of vertical electron affinities and vertical ionization potentials indicates that the polyhydroxylated derivative with eight hydroxyl groups is more stable than others. The intramolecular interaction of these derivatives considered here reveals that the more the hydroxyl groups in derivatives, the stronger the interaction in stabilizing structures. On the basis of theoretical studies, the hydroxyl groups largely enhance the infrared intensities, especially for the polyhydroxylated derivatives.

Design and synthesis of the first triply twisted Möbius annulene

As long as 50 years ago theoretical calculations predicted that Möbius annulenes with only one π surface and one edge would exhibit peculiar electronic properties and violate the Hückel rules. Numerous synthetic attempts notwithstanding, the first singly twisted Möbius annulene was not prepared until 2003. Here we present a general, rational strategy to synthesize triply or even more highly twisted cyclic π systems. We apply this strategy to the preparation of a triply twisted [24]dehydroannulene, the structure of which was confirmed by X-ray analysis. Our strategy is based on the topological transformation of 'twist' into 'writhe'. The advantage is twofold: the product exhibits a lower degree of strain and precursors can be designed that inherently include the writhe, which, after cyclization, ends up in the Möbius product. With our strategy, triply twisted systems are easier to prepare than their singly twisted counterparts.

De Proft F. Exploring the structure–aromaticity relationship in Hückel and Möbius N-fused pentaphyrins using DFT

Fine-tuning of the molecular topology and macrocyclic aromaticity of N-fused pentaphyrins by changing the number of π electrons and meso-substituents.

Performance of density functional theory in computing nonresonant vibrational (hyper)polarizabilities

A set of exchange-correlation functionals, including BLYP, PBE0, B3LYP, BHandHLYP, CAM-B3LYP, LC-BLYP, and HSE, has been used to determine static and dynamic nonresonant (nuclear relaxation) vibrational (hyper)polarizabilities for a series of all-trans polymethineimine (PMI) oligomers containing up to eight monomer units. These functionals are assessed against reference values obtained using the Møller-Plesset second-order perturbation theory (MP2) and CCSD methods. For the smallest oligomer, CCSD(T) calculations confirm the choice of MP2 and CCSD as appropriate for assessing the density functionals. By and large, CAM-B3LYP is the most successful, because it is best for the nuclear relaxation contribution to the static linear polarizability, intensity-dependent refractive index second hyperpolarizability, static second hyperpolarizability, and is close to the best for the electro-optical Pockels effect first hyperpolarizability. However, none of the functionals perform satisfactorily for all the vibrational (hyper)polarizabilities studied. In fact, in the case of electric field-induced second harmonic generation all of them, as well as the Hartree-Fock approximation, yield the wrong sign. We have also found that the Pople 6-31+G(d) basis set is unreliable for computing nuclear relaxation (hyper)polarizabilities of PMI oligomers due to the spurious prediction of a nonplanar equilibrium geometry.