Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations

In electronic structure theory, when an approximate wavefunction tends to artifactually break the symmetry of the exact Hamiltonian, the corresponding method is referred to as having a "symmetry dilemma" problem. Such types of artifacts were often reported when Hartree-Fock (HF) and the low-level post-HF methods were used, while the traditional Kohn-Sham density functional theory (KS-DFT) methods were usually found to be more resistant to this breakdown. In this work, we present a systematic study on the reliability of the doubly hybrid (DH) DFT methods for several violable cases. Almost all the commonly used B2PLYP-type (bDH) functionals are shown to have a severe "symmetry dilemma" problem and yield dramatically unreliable molecular properties, such as dipole moment, vibrational frequency, and static polarizability at the equilibrium geometry. A one-parameter bDH functional model study demonstrates that such a problem is a combined effect of the inappropriate portion of the HF exchange (over 50%) for the self-consistent field (SCF) calculation and the augmentation of the second-order perturbative contribution. It is remarkable that the XYG3-type (xDH) functionals show a good capability to resist the artifactual symmetry breaking and yield reliable molecular properties when the same critical cases are calculated. In the xDH, there are two functionals of different purposes, namely, the SCF functional and the energy functional, which have different amounts of the HF exchange and different portions of the correlation contributions. The success of the xDHs can be attributed to this flexibility in xDH construction to avoid using an improperly large portion of the HF exchange in the SCF functional. The insights gained in this work are of significance for the development of an improved DH functional.

Dynamic Symmetry Breaking Hidden in Fano Resonance of a Molecule: S1 State of Diazirine Using Quantum Wave Packet Propagation

Fano resonance in the predissociation of the S1 state of diazirine was studied by applying a time-dependent wave packet propagation method, and dynamic symmetry breaking (DSB) around the stationary structure of S1 was disclosed in a detailed analysis of this theoretical result. The DSB was found to originate in coupling between the asymmetric C-N2 stretching and CH2 wagging modes, suggesting that there is a slight time gap between ring opening and the concurrent dragging of two H atoms of the CH2 moiety. Although the depth of the double well due to DSB is just 0.011 eV, its presence noticeably affects the early time dynamics and observed spectrum.

Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory

Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory (OMP3) [U. Bozkaya, J. Chem. Phys. 135, 224103 (2011)] are presented. The OMP3 method is applied to problematic chemical systems with challenging electronic structures. The performance of the OMP3 method is compared with those of canonical second-order Møller-Plesset perturbation theory (MP2), third-order Møller-Plesset perturbation theory (MP3), coupled-cluster singles and doubles (CCSD), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)] for investigating equilibrium geometries, vibrational frequencies, and open-shell reaction energies. For bond lengths, the performance of OMP3 is in between those of MP3 and CCSD. For harmonic vibrational frequencies, the OMP3 method significantly eliminates the singularities arising from the abnormal response contributions observed for MP3 in case of symmetry-breaking problems, and provides noticeably improved vibrational frequencies for open-shell molecules. For open-shell reaction energies, OMP3 exhibits a better performance than MP3 and CCSD as in case of barrier heights and radical stabilization energies. As discussed in previous studies, the OMP3 method is several times faster than CCSD in energy computations. Further, in analytic gradient computations for the CCSD method one needs to solve λ-amplitude equations, however for OMP3 one does not since λ(ab)(ij(1))=t(ij)(ab(1)) and λ(ab)(ij(2))=t(ij)(ab(2)). Additionally, one needs to solve orbital Z-vector equations for CCSD, but for OMP3 orbital response contributions are zero owing to the stationary property of OMP3. Overall, for analytic gradient computations the OMP3 method is several times less expensive than CCSD (roughly ~4-6 times). Considering the balance of computational cost and accuracy we conclude that the OMP3 method emerges as a very useful tool for the study of electronically challenging chemical systems.

Polyradicals of polycyclic aromatic hydrocarbons as finite size models of graphene: highly open-shell nature, symmetry breaking, and enhanced-edge electron density

Properties of polyradicals (all CH bonds dissociated) of benzene and certain polycyclic aromatic hydrocarbons (PAHs) were studied. The occurrence of symmetry breaking is revealed in going from benzene and the PAHs to their polyradicals. Polyradicals would serve as finite size models of graphene with unpassivated edges in a more realistic way than the PAHs. Monoradicals (one CH bond dissociated) of benzene and all of the PAHs and higher radicals of benzene and one PAH (two to all CH bonds successively dissociated) were also investigated. Reliability of the methodology employed was ascertained by a comparison of our calculated single CH bond dissociation energy of benzene with the available previous experimental and theoretical results. Besides ground-state geometries, the aspects studied include single and successive CH bond dissociation energies, and electron density, molecular electrostatic potential (MEP), and spin density distributions. All of the monoradicals studied were found to have doublet spin multiplicity, while polyradicals with 4 to 16 rings and zigzag or mixed-type edges were found to have spin multiplicities varying from triplet to 11et. Bond lengths and bond angles of rings located at the edges are appreciably modified in going from PAHs to polyradicals. Electron density and spin density are found to be enhanced at the edges of monoradicals and polyradicals of PAHs, as found previously for PAHs. However, MEP maps of polyradicals have significantly different features from those of monoradicals and PAHs, which has a significant implication.