Computational simulation of vibrationally resolved spectra for spin-forbidden transitions

A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials

Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.

The design and synthesis of green emissive iridium(III) complexes guided by calculations of the vibrationally-resolved emission spectra

A key challenge in developing emissive materials for organic light-emitting diodes is to optimize their colour saturation, which means targeting narrowband emitters. In this combined theoretical and experimental study, we investigate the use of heavy atoms in the form of trimethylsilyl groups as a tool to reduce the intensity of the vibrations in the 2-phenylpyridinato ligands of emissive iridium(III) complexes that contribute to the vibrationally coupled modes that broaden the emission profile. An underutilised computational technique, Frank-Condon vibrationally coupled electronic spectral modelling, was used to identify the key vibrational modes that contribute to the broadening of the emission spectra in known benchmark green-emitting iridium(III) complexes. Based on these results, a family of eight new green-emitting iridium complexes containing trimethylsilyl groups substituted at different positions of the cyclometalating ligands has been prepared to explore the impact that these substituents have on reducing the intensity of the vibrations and the resulting reduction in the contribution of vibrationally coupled emission modes to the shape of the emission spectra. We have demonstrated that locating a trimethylsilyl group at the N4 or N5 position of the 2-phenylpyridine ligand damps the vibrational modes of the iridium complex and provides a modest narrowing of the emission spectrum of 8-9 nm (or 350 cm^-1). The strong correlation between experimental and calculated emission spectra highlights the utility of this computational method to understand how the vibrational modes contribute to the profile of the emission spectra in phosphorescent iridium(III) emitters.

Low-Temperature Luminescence in Organic Helicenes: Singlet versus Triplet State Circularly Polarized Emission

The experimental measurement of the photophysical and chiroptical properties of helicene-based π-conjugated emitters with electron-accepting (-CN, -py, -NO₂) or donating (TMS, NMe₂, NH₂) moieties is reported at low temperature (77 K). The samples exhibit strong circularly polarized phosphorescence in frozen solution of 2-MeTHF, with a luminescence dissymmetry factor reaching 1.6 × 10^-2 and a lifetime of over 0.46 s for the most active molecule, the nitro compound. The theoretical investigation shows that although the singlet (S₁) and triplet (T₁) excited-state emissions mainly arise from the helicene core, the rotatory strengths of the spin-allowed versus spin-forbidden emission have opposite signs. Further analysis of the spin-orbit coupling matrix elements shows that there is no strong mixing between S₁ and T₁, justifying the different signs of the rotatory strengths. In the case of the nitro compound, the enhanced phosphorescence emission is due to an efficient intersystem crossing.

Optical Activity of Spin-Forbidden Electronic Transitions in Metal Complexes from Time-Dependent Density Functional Theory with Spin-Orbit Coupling

The calculation of magnetic transition dipole moments and rotatory strengths was implemented at the zeroth-order regular approximation (ZORA) two-component relativistic time-dependent density functional theory (TDDFT) level. The circular dichroism of the spin-forbidden ligand-field transitions of tris(ethylenediamine)cobalt(III) computed in this way agrees very well with available measurements. Phosphorescence dissymmetry factors g lum and the corresponding lifetimes are evaluated for three N-heterocyclic-carbene-based iridium complexes, two of which contain helicene moieties, and for two platinahelicenes. The agreement with experimental data is satisfactory. The calculations reproduce the signs and order of magnitude of g lum , and the large variations of phosphorescence lifetimes among the systems. The electron spin contribution to the magnetic transition dipole moment is shown to be important in all of the computations.

Theoretical Investigation of the Circularly Polarized Luminescence of a Chiral Boron Dipyrromethene (BODIPY) Dye

Over the last decade, molecules capable of emitting circularly polarized light have attracted growing attention for potential technological and biological applications. The efficiency of such devices depend on multiple parameters, in particular the magnitude and wavelength of the peak of emitted light, and also on the dissymmetry factor for chiral applications. In light of these considerations, molecular systems with tunable optical properties, preferably in the visible spectral region, are particularly appealing. This is the case of boron dipyrromethene (BODIPY) dyes, which exhibit large molecular absorption coefficients, have high fluorescence yields, are very stable, both thermally and photochemically, and can be easily functionalized. The latter property has been extensively exploited in the literature to produce chromophores with a wide range of optical properties. Nevertheless, only a few chiral BODIPYs have been synthetized and investigated so far. Using a recently reported axially chiral BODIPY derivative where an axially chiral BINOL unit has been attached to the chromophore unit, we present a comprehensive computational protocol to predict and interpret the one-photon absorption and emission spectra, together with their chiroptical counterparts. From the physico-chemical properties of this molecule, it will be possible to understand the origin of the circularly polarized luminescence better, thus helping to fine-tune the properties of interest. The sensitivity of such processes require accurate results, which can be achieved through a proper account of the vibrational structure in optical spectra. Methodologies to compute vibrationally-resolved electronic spectra can now be applied on relatively large chromophores, such as BODIPYs, but require more extensive computational protocols. For this reason, particular attention is paid in the description of the different steps of the protocol, and the potential pitfalls. Finally, we show how, by means of appropriate tools and approaches, data from intermediate steps of the simulation of the final spectra can be used to obtain further insights into the properties of the molecular system under investigation and the origin of the visible bands.

The ONIOM/PMM Model for Effective Yet Accurate Simulation of Optical and Chiroptical Spectra in Solution: Camphorquinone in Methanol as a Case Study

This paper deals with the development and first validation of a composite approach for the simulation of chiroptical spectra in solution aimed to strongly reduce the number of full QM computations without any significant accuracy loss. The approach starts from the quantum mechanical computation of reference spectra including vibrational averaging effects and taking average solvent effects into account by means of the polarizable continuum model. Next, the snapshots of classical molecular dynamics computations are clusterized and one reference configuration from each cluster is used to compute a reference spectrum. Local fluctuation effects within each cluster are then taken into account by means of the perturbed matrix model. The performance of the proposed approach is tested on the challenging case of the optical and chiroptical spectra of camphorquinone in methanol solution. Although further validations are surely needed, the results of this first study are quite promising also taking into account that agreement with experimental data is reached by just a couple of full quantum mechanical geometry optimizations and frequency computations.

Vibrationally resolved electronic circular dichroism and circularly polarized luminescence spectra of a boron-fused double helicene: A theoretical study

In this work, we present the theoretical study of the vibrationally resolved absorption (ABS), electronic circular dichroism (ECD), emission (EMI), and circularly polarized luminescence (CPL) spectra of a boron-fused double helicene, with a detailed and complete discussion of the alternative possible approximate methods. Given the fact that few examples of CPL calculations exist, the potential energy surfaces (PESs) have been constructed and compared with Adiabatic (AH) and Vertical Hessian (VH) models. All the vibronic calculations have accounted for Duschinsky mixings, Franck-Condon (FC) effect and Herzberg-Teller (HT) contribution. Moreover, different HT expansions have been checked and compared, by computing the derivatives of the electric and magnetic dipole transition moments around the equilibrium geometries of the initial and final states. Our results show that both AH and VH models have well reproduced the experimental vibronic structures and VH model shows a better performance in the simulation of spectral lineshapes. They also show that HT effects dominate the shapes of EMI and CPL, tuning the relative heights of the different vibronic peaks, improving the agreement with the experiment for EMI. Moreover, HT effects are the main reason for the mirror-symmetry breaking between ECD and CPL spectra. Furthermore, interesting interference effects between FC and HT contributions have also been addressed.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?