Basis set and electron correlation effects on the polarizability and second hyperpolarizability of model open-shell pi-conjugated systems

The basis set and electron correlation effects on the static polarizability (alpha) and second hyperpolarizability (gamma) are investigated ab initio for two model open-shell pi-conjugated systems, the C(5)H(7) radical and the C(6)H(8) radical cation in their doublet state. Basis set investigations evidence that the linear and nonlinear responses of the radical cation necessitate the use of a less extended basis set than its neutral analog. Indeed, double-zeta-type basis sets supplemented by a set of d polarization functions but no diffuse functions already provide accurate (hyper)polarizabilities for C(6)H(8) whereas diffuse functions are compulsory for C(5)H(7), in particular, p diffuse functions. In addition to the 6-31G(*)+pd basis set, basis sets resulting from removing not necessary diffuse functions from the augmented correlation consistent polarized valence double zeta basis set have been shown to provide (hyper)polarizability values of similar quality as more extended basis sets such as augmented correlation consistent polarized valence triple zeta and doubly augmented correlation consistent polarized valence double zeta. Using the selected atomic basis sets, the (hyper)polarizabilities of these two model compounds are calculated at different levels of approximation in order to assess the impact of including electron correlation. As a function of the method of calculation antiparallel and parallel variations have been demonstrated for alpha and gamma of the two model compounds, respectively. For the polarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset methods bracket the reference value obtained at the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples level whereas the projected unrestricted second-order M?ller-Plesset results are in much closer agreement with the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples values than the projected unrestricted Hartree-Fock results. Moreover, the differences between the restricted open-shell Hartree-Fock and restricted open-shell second-order M?ller-Plesset methods are small. In what concerns the second hyperpolarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset values remain of similar quality while using spin-projected schemes fails for the charged system but performs nicely for the neutral one. The restricted open-shell schemes, and especially the restricted open-shell second-order M?ller-Plesset method, provide for both compounds gamma values close to the results obtained at the unrestricted coupled cluster level including singles and doubles with a perturbative inclusion of the triples. Thus, to obtain well-converged alpha and gamma values at low-order electron correlation levels, the removal of spin contamination is a necessary but not a sufficient condition. Density-functional theory calculations of alpha and gamma have also been carried out using several exchange-correlation functionals. Those employing hybrid exchange-correlation functionals have been shown to reproduce fairly well the reference coupled cluster polarizability and second hyperpolarizability values. In addition, inclusion of Hartree-Fock exchange is of major importance for determining accurate polarizability whereas for the second hyperpolarizability the gradient corrections are large.     

Quantum master equation method based on the broken-symmetry time-dependent density functional theory: application to dynamic polarizability of open-shell molecular systems

A novel method for the calculation of the dynamic polarizability (α) of open-shell molecular systems is developed based on the quantum master equation combined with the broken-symmetry (BS) time-dependent density functional theory within the Tamm-Dancoff approximation, referred to as the BS-DFTQME method. We investigate the dynamic α density distribution obtained from BS-DFTQME calculations in order to analyze the spatial contributions of electrons to the field-induced polarization and clarify the contributions of the frontier orbital pair to α and its density. To demonstrate the performance of this method, we examine the real part of dynamic α of singlet 1,3-dipole systems having a variety of diradical characters (y). The frequency dispersion of α, in particular in the resonant region, is shown to strongly depend on the exchange-correlation functional as well as on the diradical character. Under sufficiently off-resonant condition, the dynamic α is found to decrease with increasing y and/or the fraction of Hartree-Fock exchange in the exchange-correlation functional, which enhances the spin polarization, due to the decrease in the delocalization effects of π-diradical electrons in the frontier orbital pair. The BS-DFTQME method with the BHandHLYP exchange-correlation functional also turns out to semiquantitatively reproduce the α spectra calculated by a strongly correlated ab initio molecular orbital method, i.e., the spin-unrestricted coupled-cluster singles and doubles.     

Theoretical Study on Divergence Problems of Single Reference Perturbation Theories

Divergences of the single reference perturbation theories due to the addition of diffusion basis functions have been investigated for both closed-shell and open-shell molecular systems. It is found that the oscillatory range of perturbation energies of open-shell systems increases as the spin multiplicity of systems changes from 2 to 4. Feenberg transformation is exploited to treat the divergence problems. It is found numerically that within the interval of Feenberg parameter there exists a minimum perturbation order at which the perturbation series become convergent. It is also found for the open-shell systems that the magnitude of the corresponding Feenberg parameter becomes larger as the spin multiplicity of the system of interest changes from 2 to 4.     

Geometry optimizations of open-shell systems with the fragment molecular orbital method

The ability to perform geometry optimizations on large molecular systems is desirable for both closed- and open-shell species. In this work, the restricted open-shell Hartree-Fock (ROHF) gradients for the fragment molecular orbital (FMO) method are presented. The accuracy of the gradients is tested, and the ability of the method to reproduce adiabatic excitation energies is also investigated. Timing comparisons between the FMO method and full ab initio calculations are also performed, demonstrating the efficiency of the FMO method in modeling large open-shell systems.     

Revitalizing Spin Natural Orbital Analysis: Electronic Structures of Mixed-Valence Compounds,Singlet Biradicals,and Antiferromagnetically Coupled Systems

Chemical systems with open-shell electronic structure have been gaining attention these days. Their potential applications in first-row transition metal catalysis, molecular wires, photovoltaics and other potential applications have urged the adoption of a simple analysis tool to better understand their open-shell electronic structures, especially the role played by the unpaired electrons. Despite its lack of popularity, spin natural orbital (SNO) analysis is a tool we found to well-suit this purpose. We have therefore re-examined how the SNO could help us analyze some interesting open-shell systems, including mixed-valence compounds, singlet biradicals, and antiferromagnetically coupled systems. We found that some interesting patterns emerge from SNO analysis, especially those associated with exchange interaction. © 2019 Wiley Periodicals, Inc.     

Spin-adapted open-shell time-dependent density functional theory. III. An even better and simpler formulation

The recently proposed spin-adapted time-dependent density functional theory (S-TD-DFT) [Z. Li and W. Liu, J. Chem. Phys. 133, 064106 (2010)] resolves the spin-contamination problem in describing singly excited states of high spin open-shell systems. It is an extension of the standard restricted open-shell Kohn-Sham-based TD-DFT which can only access those excited states due to singlet-coupled single excitations. It is also far superior over the unrestricted Kohn-Sham-based TD-DFT (U-TD-DFT) which suffers from severe spin contamination for those excited states due to triplet-coupled single excitations. Nonetheless, the accuracy of S-TD-DFT for high spin open-shell systems is still inferior to TD-DFT for well-behaved closed-shell systems. The reason can be traced back to the violation of the spin degeneracy conditions (SDC) by approximate exchange-correlation (XC) functionals. Noticing that spin-adapted random phase approximation (S-RPA) can indeed maintain the SDC by virtue of the Wigner-Eckart theorem, a hybrid ansatz combining the good of S-TD-DFT and S-RPA can immediately be envisaged. The resulting formalism, dubbed as X-TD-DFT, is free of spin contamination and can also be viewed as a S-RPA correction to the XC kernel of U-TD-DFT. Compared with S-TD-DFT, X-TD-DFT leads to much improved results for the low-lying excited states of, e.g., N(2)(+), yet with much reduced computational cost. Therefore, X-TD-DFT can be recommended for routine calculations of excited states of high spin open-shell systems.     

How does it become possible to treat delocalized and/or open-shell systems in fragmentation-based linear-scaling electronic structure calculations? The case of the divide-and-conquer method

The authors have developed a fragmentation-based linear-scaling electronic structure calculation strategy named the divide-and-conquer (DC) method, which has been implemented into the Gamess program package. Although there are many sorts of fragmentation-based linear-scaling schemes, most of them require the charge and spin multiplicity of each fragment a priori. Therefore, their applications to delocalized and/or open-shell systems have been limited. However, the DC method is a notable exception because the distribution of electrons in the entire system is automatically determined by the universal Fermi level. In this perspective, the authors have summarized the performance of the linear-scaling self-consistent field (SCF) and post-SCF calculations of delocalized and/or open-shell systems based on the DC method. Furthermore, some future prospects of the method have been discussed.     

The simulated ab initio molecular orbital technique. VI. Open-shell radicals in the spin restricted formalism

The simulated ab initio molecular orbital (SAMO ) method previously applied to RHF closed-shell and UHF open-shell situations has been applied to open-shell radicals, such as the butyl radical and the pentyl radical, within the RHF open-shell framework. The open-shell Hartree–Fock theory is developed such that a rapidly convergent iterative method for evaluating the SAMO wave-function can be employed. Results closely parallel those for the same systems using the UHF method and are of comparable accuracy to SAMO closed-shell results.     

Frequency-dependent polarizability of open-shell atomic systems

The time-dependent variational principle is used to calculate the frequency-dependent dipole polarizabilities of 2p open-shell atomic systems. A general theory is built up to include frequency-dependent perturbations using the Roothaan–Hartree–Fock (HFR) formalism. The excitation energies corresponding to the low-lying excited states of the system are obtained from the poles of the dynamic polarizability values. A fairly accurate knowledge of transition oscillator strengths is also obtained from a knowledge of the polarizability values near the poles. The excitation energies and oscillator strengths are compared with available data.     

13顶点金属碳硼烷结构和非线性光学性质的DFT研究

采用量子化学密度泛函理论(DFT)B3LYP/6-31G(d)方法,对4,1,6-MC2B10H1213顶点金属碳硼烷几何构型进行优化,结合有限场(FF)方法计算了它们的极化率和二阶非线性光学(NLO)系数.结果表明,十个13顶点金属碳硼烷分子中1a～6a的二阶NLO系数与其构型纵向扩张呈现相同的规律.分子的前线分子轨道能级差越小,其二阶NLO系数越大.对于不同自旋态的同种金属碳硼烷分子,其偶极矩值为高自旋态大于相应的低自旋态,极化率和二阶NLO系数与自旋多重度没有一致的对应规律,自旋多重度对NLO性质影响不大.     

Design of molecular ferromagnets

A large variety of molecular ferromagnets have been synthesized since the discovery of the first organic ferromagnets, including pure organic compounds, organometallic charge-transfer complexes, metal complex-organic radical compounds, and transition metal complexes coupled to organic radicals. Besides, there are many reports on the observation of ferromagnetism in polymers and organic matrix composites. Molecular ferromagnets have great potential in different areas of technology such as low frequency magnetic shielding, magnetic imaging, magneto-optics and information storage. We provide a brief review on the current strategies for the design of molecular (organic) ferromagnets. This includes exploiting the inherent advantages of molecular systems, such as the ability to fine-tune the properties at the molecular level, and to control dimensionality, supramolecular structuring and hierarchy of spin interactions etc. for carrying out structural modifications and chemical functionalisations of stable open-shell molecules in order to generate supramolecular structures in which the natural prediction for antiparallel spin alignment (antiferromagnetism) is avoided.     

A mathematical discussion of Pons Viver's implementation of Löwdin's spin projection operator

Recently, the molecular electronic structure theories for efficiently treating static (or strong) correlation in a black-box manner have attracted much attention. In these theories, a spin projection operator is used to recover the spin symmetry of a broken-symmetry Slater determinant. Very recently, Pons Viver proposed the practical and exact implementation of Löwdin's spin projection operator (Int. J. Quantum Chem. 2019, 119, e25770). In the present study, we attempt to supply mathematical proofs to Pons Viver's proposals and show a condition for establishing Pons Viver's implementation. Moreover, we explicitly derive the (spin projected) extended Hartree-Fock (EHF) equations on the basis of the model of common orbitals (ie, closed-shell orbitals used in the restricted open-shell Hartree-Fock (ROHF) method), which was combined by Pons Viver with the EHF method.     

Doping-enhanced hyperpolarizabilities of silicon clusters: a global ab initio and density functional theory study of Si10 (Li, Na, K)n (n=1, 2) clusters

A global theoretical study of the (hyper)polarizabilities of alkali doped Si(10) is presented and discussed. First, a detailed picture about the low lying isomers of Si(10)Li, Si(10)Na, Si(10)K, Si(10)Li(2), Si(10)Na(2), and Si(10)K(2) has been obtained in a global manner. Then, the microscopic first (hyper)polarizabilities of the most stable configurations have been determined by means of ab initio methods of high predictive capability such as those based on the M?ller-Plesset perturbation and coupled cluster theory, paying extra attention to the (hyper)polarizabilities of the open shell mono-doped systems Si(10)Li, Si(10)Na, Si(10)K, and the influence of spin contamination. These results were used to assess the performance of methods of low computational cost based on density functional theory (DFT) in the reliable computation of these properties in order to proceed with an in-depth study of their evolution as a function of the alkali metal, the cluster composition, and the cluster structure. The most interesting outcomes of the performed (hyper)polarizability study indicate that while alkali doping leaves the per atom polarizability practically unaffected, influences dramatically the hyperpolarizabilities of Si(10). The lowest energy structures of the mono-doped clusters are characterized by significantly enhanced hyperpolarizabilities as compared to the analogue neutral or charged bare silicon clusters Si(10) and Si(11), while, certain patterns governed by the type and the number of the doping agents are followed. The observed hyperpolarizability increase is found to be in close connection with specific cluster to alkali metal charge transfer excited states and to the cluster structures. Moreover, an interesting correlation between the anisotropy of the electron density, and the hyperpolarizabilities of these systems has been observed. Finally, it is important to note that the presented method assessment points out that among the various DFT functionals used in this work, (B3LYP, B3PW91, BhandHLYP, PBE0, CAM-B3LYP, LC-BLYP, LC-BPW91) only B3PW91 and PBE0 out of the seven provided a consistent quantitative performance for both polarizabilities and hyperpolarizabilities with respect to the ab initio methods utilized here. On the other hand, the long range corrected functionals LC-(U)BLYP and LC-(U)BPW91 (μ = 0.47) failed to supply quantitatively accurate hyperpolarizability results in all the studied clusters while the CAM-(U)B3LYP functional performs satisfactory only in the case of the Na and K doped systems.     

Redox-based spin diversity in a 6-oxophenalenoxyl system: generation,ESR/ENDOR/TRIPLE,and theoretical studies of 2,5,8-tri-tert-butylphenalenyl- 1,6-bis(olate) salts

[reaction: see text] Novel open-shell molecular salts, 2,5,8-tri-tert-butylphenalenyl-1,6-bis(olate) salts, were designed on the basis of the 6-oxophenalenoxyl system and generated by the chemical reduction of 6-hydroxyphenalenone derivatives. ESR/ENDOR/TRIPLE measurements and DFT calculations provide unequivocal determination of the structure and spin density distribution, which demonstrate redox-based spin diversity of the 6-oxophenalenoxyl system.     

Analysis of wave functions for open-shell molecules

During the past decade we have looked at several ways to track the distribution of unpaired electrons during chemical reactions and in different spin states. These methods were inspired by our previous work on singlet di-radicals where the spin density is zero yet there are clearly singly occupied orbitals. More recently we have been concerned with analysis of wave functions for single molecule magnets. This review discusses the mathematical framework by which open-shell systems can be described, in addition to methods that extract the effectively unpaired electron density, the spin state of atoms in a molecule, and other useful properties from a molecular wave function. Some of the difficulties associated with using broken spin Slater determinants to evaluate the exchange coupling parameters in the Heisenberg Hamiltonian are also mentioned.     

Theoretical investigation of the (hyper)polarizabilities of pyrrole homologues C4H4XH (X = N, P, As, Sb, Bi). A coupled-cluster and density functional theory study

Geometries, inversion barriers, static and dynamic electronic and vibrational dipole polarizability (alpha), and first (beta) and second (gamma) hyperpolarizability of the pyrrole homologues C(4)H(4)XH (X = N, P, As, Sb, Bi) have been calculated by Hartree-Fock, M?ller-Plesset second-order perturbation theory, coupled-cluster theory accounting for singles, doubles, and noniterative triple excitations methods, as well as density functional theory using B3LYP and PBE1PBE functionals and Sadlej's Pol and 6-311G basis sets. Relativistic effects on the heavier homologues stibole and bismole have been taken into account within effective core potential approximation. The results show that the electronic (hyper)polarizabilities monotonically increase with the atomic number of the heteroatom, consistent with the decrease in the molecular hardness. Ring planarization reduces the carbon-carbon bond length alternation of the cis-butadienic unit, enhancing the electronic polarizability values (alpha(e)) by 4-12% and the (hyper)polarizability values (and gamma(e)) by 30-90%. Pure vibrational and zero-point vibrational average contributions to the (hyper)polarizabilities have been determined within the clamped nucleus approach. In the static limit, the pure vibrational hyperpolarizabilities have a major contribution. Anharmonic corrections dominate the pure vibrational hyperpolarizabilities of pyrrole, while they are less important for the heavier homologues. Static and dynamic electronic response properties of the pyrrole homologues are comparable to or larger than the corresponding properties of the furan and cyclopentadiene homologue series.     

Performance of the M06 family of exchange-correlation functionals for predicting magnetic coupling in organic and inorganic molecules

The performance of the M06 family of exchange-correlation potentials for describing the electronic structure and the Heisenberg magnetic coupling constant (J) is investigated using a set of representative open-shell systems involving two unpaired electrons. The set of molecular systems studied has well defined structures, and their magnetic coupling values are known experimentally. As a general trend, the M06 functional is about equally as accurate as B3LYP or PBE0. The performance of local functionals is important because of their economy and convenience for large-scale calculations; we find that M06-L local functional of the M06 family largely improves over the local spin density approximation and the generalized gradient approximation.