We have conducted an experimental and theoretical study on first-and second-order Raman scattering of zinc blende and wurtzite ZnS. Based on the calculated phonon band structure, phonon density of states, and symmetry selection rules, we have clearly identified for the first time the origins of these vibration modes in the second-order Raman spectra from these two polymorphs. For zinc blende ZnS, it is found that the previously estimated frequency of the LA mode at X point in the Brillouin-zone boundary is much smaller than the value obtained from other experiments and our calculation. Considering the involvement of LA phonon at X point, we reassign the second-order Raman active modes and some other modes which have not yet been understood so far. This work clarifies some of the controversial Raman mode assignments in zinc blende and wurtzite ZnS.

The CASPT2/CASSCF method with the 6-311G basis set and an active space up to (14, 11) was used to explore the ultrafast internal conversion mechanism for excited 9H-adenine. Three minima, two transition states, and seven conical intersections were obtained to build up the two deactivation pathways for the internal conversion mechanism. Special efforts were made to explore the excited-state potential energy surfaces near the Franck-Condon region and determine the various barriers in the processes of deactivation. The barrier required from the 1pipi (1La) state to deactivate nonradiatively is found to be lower than that required from the 1pipi (1Lb) state. On 250 nm excitation, the 1pipi (1La) state is populated, and the transition from 1pipi (1La) to the lowest 1npi state involves very low barriers, which may account for the observed short (<50 fs) lifetime of the 1pipi excited state. The deactivation of the lowest 1npi state is required to overcome a barrier of 3.15 kcal/mol, which should be responsible for the 750 fs lifetime of the npi excited state. On 267 nm excitation, the vibrationally active 1pipi (1Lb) state is populated. Excitation at 277 nm prepares the 1pipi (1Lb) state without much excessive vibrational energy, which may be responsible for the observed >2 ps lifetime.

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.

An efficient implementation of the generalized energy-based fragmentation (GEBF) approach (Li, W.; Li, S.; Jiang, Y. J Phys. Chem. A 2007, 111, 2193) for treating a wide range of large molecules is presented. In this implementation, the fragmentation process can be automatically done for a general molecule, with only some functional groups defined by users. A new and fast scheme is designed for the generation of various subsystems and the derivation of their coefficients. The newly implemented GEBF approach has been applied to several large molecules including proteins, nucleic acids, and supermolecules with fused aromatic rings. Test calculations within the Hartree-Fock (HF) and density functional theory (DFT) framework demonstrate that the GEBF approach can provide reasonably accurate ground-state energies and optimized structures, which are in good agreement with those from conventional HF or DFT calculations. The GEBF approach implemented in this work can now be employed by nonexpert users to compute energies, optimized structures, and some molecular properties at various ab initio levels for a broad range of large molecules on ordinary PC workstations.