Symmetry-adapted-cluster/symmetry-adapted-cluster configuration interaction methodology extended to giant molecular systems: ring molecular crystals

The symmetry adapted cluster (SAC)/symmetry adapted cluster configuration interaction (SAC-CI) methodology for the ground, excited, ionized, and electron-attached states of molecules was extended to giant molecular systems. The size extensivity of energy and the size intensivity of excitation energy are very important for doing quantitative chemical studies of giant molecular systems and are designed to be satisfied in the present giant SAC/SAC-CI method. The first extension was made to giant molecular crystals composed of the same molecular species. The reference wave function was defined by introducing monomer-localized canonical molecular orbitals (ml-CMO's), which were obtained from the Hartree-Fock orbitals of a tetramer or a larger oligomer within the electrostatic field of the other part of the crystal. In the SAC/SAC-CI calculations, all the necessary integrals were obtained after the integral transformation with the ml-CMO's of the neighboring dimer. Only singles and doubles excitations within each neighboring dimer were considered as linked operators, and perturbation selection was done to choose only important operators. Almost all the important unlinked terms generated from the selected linked operators were included: the unlinked terms are important for keeping size extensivity and size intensivity. Some test calculations were carried out for the ring crystals of up to 10 000-mer, confirming the size extensivity and size intensivity of the calculated results and the efficiency of the giant method in comparison with the standard method available in GAUSSIAN 03. Then, the method was applied to the ring crystals of ethylene and water 50-mers, and formaldehyde 50-, 100-, and 500-mers. The potential energy curves of the ground state and the polarization and electron-transfer-type excited states were calculated for the intermonomer distances of 2.8-100 A. Several interesting behaviors were reported, showing the potentiality of the present giant SAC/SAC-CI method for molecular engineering.     

Singularity-free analytical energy gradients for the SAC/SAC-CI method: coupled perturbed minimum orbital-deformation (CPMOD) approach

A new procedure for evaluating energy gradients in a singularity-free manner is presented for use in the SAC/SAC-CI program in which computational dimensions are reduced by the perturbation selection method. The singularity in the energy gradients stemming from a breakdown of the unitary invariance is effectively removed by the minimum orbital-deformation (MOD) method proposed in the previous study. All calculations can be done analytically via new two sets of linear equations combined with the coupled-perturbed Hartree–Fock method. Geometry optimizations for malonaldehyde in the ground and lowest singlet excited states are performed by the new method.     

BH分子X 1Σ+、A 1Π和B 1Σ+ 态的势能函数

利用SAC/SAC-CI方法，使用D95++、6-311++g及cc-PVTZ等基组，对BH分子的基态(X ¹Σ⁺)、第一简并激发态(A ¹Π)及第二激发态(B ¹Σ⁺)的平衡结构和谐振频率进行了优化计算． 通过对三个基组计算结果的比较，得出了cc-PVTZ基组为三个基组中的最优基组的结论；使用cc-PVTZ基组，利用SAC的GSUM(group sum of operators)方法对基态(X ¹Σ⁺)， SAC-CI的GSUM方法对激发态(A ¹Π 和B ¹Σ⁺)进行单点能扫描计算, 用正规方程组拟合Murrell-Sorbie函数，得到了相应电子态的完整势能函数；从得到的势能函数计算了与基态(X ¹Σ⁺)、第一简并的激发态(A ¹Π)和第二激发态(X ¹Σ⁺)相对应的光谱常数(Be、αe、ωe 和ωeχe)，结果与实验数据较为一致． 其中基态、第一激发态与实验数据吻合得较好．     

BF分子X1Σ＋, A1P和B1Σ＋态的势能函数

使用SAC/SAC-CI和D95＋＋, 6-311＋＋g, 6-311＋＋g**及D95(d)基组, 分别对BF分子的基态X¹Σ^＋、第一简并激发态A¹Π和第二激发态B¹Σ^＋的平衡结构和谐振频率进行优化计算. 对所有计算结果进行比较, 得出6-311＋＋g**基组为最优基组. 运用6-311＋＋g**基组和SAC方法对基态X¹Σ^＋, SAC-CI方法对激发态A¹Π和B¹Σ^＋进行单点能扫描计算, 并用正规方程组拟合Murrell-Sorbie函数, 得到相应电子态的势能函数解析式, 由得到的势能函数计算了与X¹Σ^＋, A¹Π和B¹Σ^＋态相对应的光谱常数, 结果与实验数据较为一致.     

LiH分子X 1Σ+ 、A 1Σ+和B 1Π态的势能函数

The energies, equilibrium geometries and harmonic frequencies of the three electronic states (the ground state X 1Σ+, the first excitation state A 1Σ+ and the second excitation degenerate state B 1Π) of LiH molecule have been calculated by using the GSUM (Group Sum of Operators) method of SAC/ SAC-CI with the basis sets D95(d), 6-311G**, and cc-PVTZ. Comparing with the above-mentioned three basis sets, the conclusion is gained that the basis set D95(d) is the most suitable for the energy calculation of LiH molecule. The whole potential curves for these three electronic states are further scanned, using SAC/D95(d) method for the ground state and SAC-CI/D95(d) methods for the excited states. Murrell-Sorbie function were fitted using a least square and then the spectroscopy constants are calculated, which are in good agreement with the experimental data.     

Cluster expansion of the wave function. Electron correlations in singlet and triplet excited states,ionized states,and electron attached states by SAC and SAC–CI theories

The symmetry-adapted-cluster (SAC ) and SAC –CI theories reported previously have been applied to the study of electron correlations in ground state, singlet and triplet excited states, ionized state, and electron attached state. Formulas for calculations of one-electron properties and transition properties from the SAC and SAC –CI wave functions are given. Calculations were carried out for the ground and Rydberg excited states of water and its positive and negative ions, with the use of the simpler computational scheme than the previous one. The results compare well with experiments.     

Inner-shell ionizations and satellites studied by the open-shell reference symmetry-adapted cluster/symmetry-adapted cluster configuration-interaction method

Open-shell reference version of the symmetry-adapted cluster (SAC) and SAC-configuration-interaction (CI) methods, termed open-shell reference (OR)-SAC and OR-SAC-CI methods, are developed and applied to inner-shell ionizations of CH4, NH3, H2O, and HF. The inner-shell ionization potentials and spectra calculated by the OR-SAC and OR-SAC-CI methods are in excellent agreement with the experimental data. Including both of the electron correlation and orbital relaxation is important for quantitative agreements. Timing comparisons with the SAC-CI general-R calculations that give similar high accuracies show an efficiency of the present OR-SAC and OR-SAC-CI methods.     

Singly and doubly excited states of butadiene, acrolein, and glyoxal: Geometries and electronic spectra

Excited-state geometries and electronic spectra of butadiene, acrolein, and glyoxal have been investigated by the symmetry adapted cluster configuration interaction (SAC-CI) method in their s-trans conformation. Valence and Rydberg states below the ionization threshold have been precisely calculated with sufficiently flexible basis sets. Vertical and adiabatic excitation energies were well reproduced and the detailed assignments were given taking account of the second moments. The deviations of the vertical excitation energies from the experiment were less than 0.3 eV for all cases. The SAC-CI geometry optimization has been applied to some valence and Rydberg excited states of these molecules in the planar structure. The optimized ground- and excited-state geometries agree well with the available experimental values; deviations lie within 0.03 A and 0.7 degrees for the bond lengths and angles, respectively. The force acting on the nuclei caused by the excitations has been discussed in detail by calculating the SAC-CI electron density difference between the ground and excited states; the geometry relaxation was well interpreted with the electrostatic force theory. In Rydberg excitations, geometry changes were also noticed. Doubly excited states (so-called 2 (1)A(g) states) were investigated by the SAC-CI general-R method considering up to quadruple excitations. The characteristic geometrical changes and large energetic relaxations were predicted for these states.     

Electronic structure and optical properties of chelating heteroatomic conjugated molecules: a SAC-CI study

The electronic structure and optical properties of 13 chelating heteroatomic conjugated molecules such as pyridine, benzoxazole, and benzothiazole derivatives, which are used as C–N ligands in organometallic compounds, have been investigated. The geometries of the ground and first excited states were obtained by the DFT and CIS methods, respectively, followed by the SAC-CI calculations of the transition energies for absorption and emission. For six compounds whose experimental data are available, the SAC-CI calculations reproduced the experimental values satisfactorily with deviations of less than 0.3 eV for absorption and 0.1 eV for emission except for benzoxazoles. For other molecules, the theoretical absorption and emission spectra were predicted. The lowest ππ* excited-state geometries was calculated to be planar for most of the molecules with two or three conjugated rings connected by single bond. The geometry change due to the ππ* excitation was qualitatively interpreted by electrostatic force theory based on SAC/SAC-CI electron density difference. The excitations are relatively localized in the central region and in the lowest ππ* excited state, the inter-ring single bond shows large change, with a contraction of 0.05–0.09 Å. The present calculations provide reliable information regarding the energy levels of these chelating heteroatomic conjugated compounds.     

Investigation of the electronic spectra and excited-state geometries of poly(para-phenylene vinylene) (PPV) and poly(para-phenylene) (PP) by the symmetry-adapted cluster configuration interaction (SAC-CI) method

The symmetry-adapted cluster-configuration interaction (SAC-CI) method has been used to investigate the optical and geometric properties of the oligomers of poly(para-phenylene vinylene) (PPV) and poly(para-phenylene) (PP). Vertical singlet and triplet absorption spectra and emission spectra have been calculated accurately; the mean average deviation from available experimental results lies within 0.2 eV. The chain length dependence of the transition energies has been improved in comparison to earlier TDDFT and MRSDCI calculations. The present analysis suggests that conventional TDDFT with the B3LYP functional should be used carefully, as it can provide inaccurate estimates of the chain length dependence of the excitation energies of these molecules with long pi conjugation. The T1 state was predicted to be at a lower energy, by 1.0-1.5 eV for PPV and by 0.9-1.7 eV for PP, than the S1 state, which indicates a localized T1 state with large exchange energy. By calculating the SAC-CI electron density difference between the ground and excited states, the geometry relaxations due to excitations can be analyzed in detail using electrostatic force theory. For trans-stilbene, the doubly excited 21Ag state was studied, and the calculated transition energy of 4.99 eV agrees very well with the experimental value of 4.84 eV. In contrast to previous ab initio calculations, we predict this doubly excited 21Ag state to lie above the 11Bu state.     

He2+和He2++分子离子基态和激发态的特性研究

采用SAC/SAC-CI方法在CC-PV5Z基组下, 计算研究了He2+、He2++的基态及低激发态的分子特性, 给出了其基态和一些激发态的势能函数和光谱数据(Be、αe、ωe和ωeχe). 从群论出发推导了相应状态的离解极限；与已有实验结果的He2+(X2Σu+)相比, 计算结果令人满意. 还计算了激发态2Πu、4Σu+和4Πg的结构与光谱数据. 对于He2++, 计算的九个电子态中只有三个态(X1Σg+、1Σg+和1Σu+)属束缚态, 并得到了其光谱常数. 用价键理论模型的不相交规则对He2++基态的势能曲线极大点产生的原因做了较好的分析.     

SAC-CI study of the excited states of free base tetrazaporphin

The SAC (symmetry adapted cluster)/SAC-CI method is applied to calculations of the ground and excited states of free base tetrazaporphin (FBTAP). The electronic spectrum is reproduced in fairly good agreement with experiments, and unknown absorption bands in the energy region higher than 4 eV are predicted. The effect of meso-tetraaza-substitution on the excited states of porphin macrocycle explains why the compound is colored, as in the closely related compounds, phthalocyanines, used as pigments.     

Active-space symmetry-adapted-cluster configuration-interaction and equation-of-motion coupled-cluster methods for high accuracy calculations of potential energy surfaces of radicals

The electron-attached (EA) and ionized (IP) symmetry-adapted-cluster configuration-interaction (SAC-CI) methods and their equation-of-motion coupled-cluster (EOMCC) analogs provide an elegant framework for studying open-shell systems. As shown in this study, these schemes require the presence of higher-order excitations, such as the four-particle-three-hole (4p-3h) or four-hole-three-particle (4h-3p) terms, in the electron attaching or ionizing operator R in order to produce accurate ground- and excited-state potential energy surfaces of radicals along bond breaking coordinates. The full inclusion of the 4p-3h/4h-3p excitations in the EA/IP SAC-CI and EOMCC methods leads to schemes which are far too expensive for calculations involving larger radicals and realistic basis sets. In order to reduce the large costs of such schemes without sacrificing accuracy, the active-space EA/IP EOMCC methodology [J. R. Gour et al., J. Chem. Phys. 123, 134113 (2005)] is extended to the EA/IP SAC-CI approaches with 4p-3h/4h-3p excitations. The resulting methods, which use a physically motivated set of active orbitals to pick out the most important 3p-2h/3h-2p and 4p-3h/4h-3p excitations, represent practical computational approaches for high-accuracy calculations of potential energy surfaces of radicals. To illustrate the potential offered by the active-space EA/IP SAC-CI approaches with up to 4p-3h/4h-3p excitations, the results of benchmark calculations for the potential energy surfaces of the low-lying doublet states of CH and OH are presented and compared with other SAC-CI and EOMCC methods, and full CI results.     

Ab initio study of molecular structures and excited states in anthocyanidins

The structural and electronic characters of four types of hydroxyl group-substituted anthocyanidins (pelargonidin, cyanidin, delphinidin, and aurantinidin) were examined using quantum chemical calculations. For these cationic molecules, both the planar and non-planar structures in the electronic ground state were determined at the B3LYP/D95 level of theory. We revealed that the planar structure is slightly more stable than the non-planar structure for each molecule. For the optimized planar structures, single excitation-configuration interaction (SE-CI) based on the restricted Hartree-Fock (RHF) wave function was evaluated and the electronic character in the low-excited states was discussed in terms of the MO theory. Symmetry adapted cluster (SAC)/SAC-CI calculations were also carried out to estimate the excitation energies precisely. The results showed that hydroxylation of the phenyl group causes a change in the excitation energies without taking the solvent effects into account. The results are in agreement with spectral experiments and previous MO calculations.     

Nonequilibrium solvation for vertical photoemission and photoabsorption processes using the symmetry-adapted cluster-configuration interaction method in the polarizable continuum model

In this paper, we present the theory and implementation of a nonequilibrium solvation model for the symmetry-adapted cluster (SAC) and symmetry-adapted cluster-configuration interaction (SAC-CI) method in the polarizable continuum model. For nonequilibrium solvation, we adopted the Pekar partition scheme in which solvent charges are divided into dynamical and inertial components. With this nonequilibrium solvation scheme, a vertical transition from an initial state to a final state may be described as follows: the initial state is described by equilibrium solvation, while in the final state, the inertial component remains in the solvation for the initial state; the dynamical component will be calculated self-consistently for the final state. The present method was applied to the vertical photoemission and absorption of s-trans acrolein and methylenecyclopropene. The effect of nonequilibrium solvation was significant for a polar solvent.     

Near-infrared absorption in symmetric squarylium and croconate dyes: a comparative study using symmetry-adapted cluster-configuration interaction methods

Symmetric croconate (CR) and squarylium dyes (SQ) are well-known near-infrared (NIR) dyes and, in general, are considered to be donor-acceptor-donor type molecules. It is established in the literature that CR dyes absorb in a longer wavelength region than the corresponding SQ dyes. This has been attributed to the CR ring being a better acceptor than the SQ ring. Thus increasing the donor capacity should lead to a bathochromic shift in both SQ and CR. On the other hand, some experiments reported in the literature have revealed that increasing the conjugation in the donor part of the SQ molecule leads first to red shift, which upon a further increase of the conjugation changes to a blue shift. Hence, to understand the role of the central ring and the substitutions in the absorption of these dyes, we carried out high-level symmetry-adapted cluster-configuration interaction (SAC-CI) calculations of some substituted SQ and CR dyes and compare the absorption energy with the existing experimental data. We found that there is very good agreement. We also carried out SAC-CI calculations of some smaller model molecules, which contain the main oxyallyl substructure. We varied the geometry (angle) of the oxyallyl subgroup and the substitution in these model molecules to establish a correlation with the bathochromic shift. We found that the charge transfer is very small and does not play the key role in the red shift, but on the other hand, the perturbation of the HOMO-LUMO gap (HLG) from both the geometry and substitution seems to be responsible for this shift. We suggest as a design principle that increasing the donor capacity of the groups may not help in the red shift, but introducing groups which perturb the HLG and decrease it without changing the MO character should lead to a larger bathochromic shift.     

Theoretical study of valance photoelectron spectra of hypoxanthine, xanthine, and caffeine using direct symmetry-adapted cluster/configuration interaction methodology

UV photoelectron spectra of hypoxanthine, xanthine, and caffeine, up to 20 eV, were calculated and compared with the experimental spectra reported in literature. The calculations were performed using a novel version of the quantum mechanical symmetry-adapted cluster/configuration interaction (SAC-CI) method termed, direct SAC-CI. The Duning/Huzinaga valance double-zeta D95+(d,p) Gaussian basis set was also employed with this method. The ionization energies and intensities were calculated, and the corresponding spectral bands were assigned. Natural bonding orbital (NBO) calculations were employed for better spectral band assignment. The calculated ionization energies and intensities reasonably produced the experimental photoelectron spectra.     

C4Cl: bent or linear?

The ground state structure for the CCCCCl radical was computed by using symmetry-adapted cluster configuration-interaction (SAC-CI) theory along with density functional theory to overcome the differences raised in the recently published paper [Y. Sumiyoshi et al., Chem. Phys. Lett. 414, 82 (2005)] between the theory and the experiment. SAC-CI results clearly support the earlier experimental conclusion that the radical has the bent ground state structure corresponding to 2Pi symmetry. Contrarily, probably due to spin contamination, mixing of a bent doublet ground state with the quartet components of a linear structure, coupled-cluster singles and doubles (CCSD) calculations were unable to provide reliable results. Results obtained using density functional theory also show that the radical has a bent structure. Some low-lying doublet excited states were also studied using the SAC-CI theory. The energy difference between the ground Pi state and the nearby Sigma state is around 0.2 eV. The excitation energy for the transition with the largest oscillator strength agrees with the strongest absorption peak.