Search

The evaluation of four high-level composite methods based on the modification of Gaussian-3 (G3) theory for radicals and 18 exchange-correlation density functionals, including modern long-range and dispersion-corrected functionals, in the modelization of singlet diradicals has been performed in this work. Structural parameters and properties such as singlet-triplet gaps, electron affinities, ionization potentials, dipole moments, enthalpies of formation, and bond dissociation energies have been calculated in a set of six well-characterized singlet diradicals, and benchmarked against experimental data and wavefunction-based CASSCF/CASPT2 calculations. The complexity of the open-shell singlet ground state is revealed in the difficulties to properly represent the diradical character reported by some DFT functionals, specially those that do not comprise a certain amount of Hartree-Fock exchange in their formulation. We find that STGs, EAs, dipole moments, and thermochemical properties are, in general, satisfactorily calculated, while for IPs larger deviations with respect to the experiments are found in all cases. The best overall performance is accounted for by hybrid functionals, including some of the long-range corrected functionals, but also pure functionals, comprising the kinetic energy density in their formulation, are found to be competent. Composite methods perform satisfactorily, especially G3(MP2)-RAD and G3X(MP2)-RAD, which calculate singlet-triplet gaps and electron affinities more accurately. On the other hand, G3-RAD and G3X-RAD provide better ionization potentials. This study emphasizes that the use of recently developed functionals, within the broken symmetry approximation, is an appropriate tool for the simulation of organic singlet diradicals, with similar accuracy compared to more expensive composite methods. Nevertheless, suitable selection of the methodology is still crucial for the accomplishment of accurate results.

The triplet-triplet absorption spectra of three newly synthesized N-substituted 4,5,6,7-tetrachlorophthalimides (TCP) were measured experimentally and calculated with density functional theory. The heavy atom effect increases the intersystem crossing rate, and the transient triplet absorbance could be measured. Fluorescence emission was not observed. The transient absorption spectra show two peaks in the region of 385-410 nm and 770-830 nm having a red shift with increasing size of N substituent. The singlet and triplet states and their vertical transitions were investigated theoretically by NEVPT2/CASSCF(12,9) and CAM-B3LYP methods for a benchmark molecule N-Phenyl-maleimide, and by CAM-B3LYP for N-Phenyl-TCP. It is shown that the singlet excited S1 state and the pi-pi* triplet state have similar geometries, thus intersystem crossing is most likely for N-Phenyl-TCP. No crossing between the singlet ground state and the singlet excited states could be found on the linear synchronous transit paths of S0-S1-T1-S0 states using CAM-B3LYP method. The computed singlet ground state absorption and triple-triplet absorption spectra agree well with the experimental ones, in both spectra the first state is a pi-pi* state while the second state is the n-pi* state with an energetic difference of 0.1 eV.

The results of STO-6G valence bond calculations with Lewis-type valence bond structures are reported for the S=0 spin ground states of the linear NNO (C∞v) and NON (D∞h) isomers of N2O. It is calculated that the singlet diradical character of the NON isomer is substantially larger than for the NNO isomer. The results of B3LYP/6-31G(d) density functional calculations are also reported for the C∞v and D∞h ground states, the single-determinant approximations to the lowest-energy S=1 spin excited states of these two isomers, and for the triangular C2v isomer. The calculated bond lengths are in accord with qualitative valence bond considerations.

Cyclic organic peroxides are a broad and highly sought-after class of peroxide compounds that present high reactivity and even explosive character. The unusually high reactivity of these peroxides can generally be attributed to the rupture of O–O bonds. Cyclic diperoxides are a very interesting series of substituted compounds in which tetroxane is the most prominent member. Gas-phase thermolysis of the simplest substituted member of the series [3-methyl-1,2,4,5-tetroxane or methylformaldehyde diperoxide (MFDP)] has been observed to yield one acetaldehyde, one formaldehyde, and one oxygen molecule as reaction products. DFT at the 6-311 + G** level of theory using the BHANDHLYP correlation–exchange functional was applied via the Gaussian09 program to calculate the critical points of the potential energy surface (PES) of this reaction. Equatorial and axial isomers were studied. The singlet state PES of MFDP was calculated, and an open diradical structure was found to be the first intermediate in a stepwise reaction. Two PESs were subsequently obtained: singlet state (S) and triplet state (T) PESs. After that, two alternative stepwise reactions were found to be possible: 1) one in which either an acetaldehyde, or 2) formaldehyde molecule is initially formed. For second one, exothermic reactions were observed for both the S and T PESs. The reaction products include a oxygen molecule in either S or T state, with the T reaction being the most exothermic. When calculations were performed at the CASSCF(10,10)/6-311 + G** level, spin–orbit coupling permitted S to T crossing at the open diradical intermediate stage, a non-adiabatic reaction was observed, and lower activation energies and higher exothermicity were generally seen for the T PES than for the S PES. These results were compared with the corresponding results for tetroxane. The spin–orbit coupling of MFDP and tetroxane yielded identical values, so it appears that the methyl substituent does not have any effect on this coupling.

The geometries, energies, and electronic properties of the two possible configurations of bis-[dibenzo[a.i]fluorenylidene] were investigated theoretically by density functional theory DFT B3LYP at the UB3LYP/6-311 + G(2d,p) // UB3LYP/6-31 + G(d,p) level of theory. According to the performed calculations, it was found that the singlet is 3.4 kcal mol(-1) lower in energy compared to triplet state at room temperature. This gap is compared with those of other alkenes like ethylene, (61.9 kcal mol(-1)) tetra-tert-butyethylene, (6.4 kcal mol(-1)) and bis-fluorenylidene (19.5 kcal mol(-1)). These results confirm the experimental findings of the paramagnetic properties determined by Franzen and Joschek. The low singlet-triplet gap in the case of bis-[dibenzo[a.i]fluorenylidene] is the result of a steric destabilization of the singlet due to strain and stabilization of the triplet electronic state by delocalization of each free electron within each aromatic moiety. This correlates with the special electronic structure of the triplet state of this compound, where facial interaction of two hydrogen atoms lying close to the lobes of each p-orbital occupied with a single electron at the distorted double bond in the triplet electronic state.

Three fluoranthenes and one substituted fluoranthene, 2,2-dimethyl-2H-dibenzo[cd,k]fluoranthene, were investigated using the unrestricted symmetry-broken and complete active space methods. It was shown that four Kekulean hydrocarbons are diradicals, implying that their ground state is a triplet. In the energetically less favorable singlet state these hydrocarbons exhibit pronounced diradical character. This occurance is explained with the tendency of the investigated molecules to delocalize their π-electrons. This leads to aromatic stabilization which is stronger than destabilization due to unpaired electrons. Our results for 2,2-dimethyl-2H-dibenzo[cd,k]fluoranthene are in excellent accord with experimental findings of McMaster et al. concerning this compound.

The polycyclic p-quinodimethanes are proposed to be the novel candidates of the high-performance nonlinear optical (NLO) materials because of their large third order polarizabilities (γ). We investigate the switchable NLO responses of a series of polycyclic p-quinodimethanes with redox properties by employing the density functional theory (DFT). The polycyclic p-quinodimethanes are forecasted to exhibit obvious pure diradical characters because of their large y0 index (the y0 index is a value between 0 [closed-shell state] and 1 [pure biradical state]). The γ values of these polycyclic p-quinodimethanes and their corresponding one-electron and two-electron reduced/oxidized species are calculated by the (U)BHandHLYP method. The γ values of polycyclic p-quinodimethanes and their corresponding one-electron reduced species are all positive and significantly different. The large differences of the γ values are due to a change in the transition energy and are related to the different delocalization of the spin density, which demonstrates that the NLO switching is more effective on one-electron reduction reactions. Therefore, the study on these polycyclic p-quinodimethanes provides a guideline for a molecular design of highly efficient NLO switching. Figure The NLO switching is more effective on one-electron reduction reaction Electronic supplementary material The online version of this article (doi:10.1007/s00894-013-2035-1) contains supplementary material, which is available to authorized users.