价键理论中的组态相互作用

近 2 0年来 ,从头计算水平的价键 (VB)方法得到了人们的重视 ,并广泛应用于化学反应等问题的研究[1～ 5] ,然而目前价键理论的计算方法仍然很不完善 .用 VBSCF方法进行计算虽然比较简单 ,能正确地描述化学反应的形成机理 ,但数值结果不理想 ;而用 BOVB方法[4 ] 进行计算虽然可以得到较好的计算结果 ,但存在收敛困难等问题 .分子轨道理论中的组态相互作用是一种简单直接的电子相关能计算方法 ,显然这一方法可以应用于价键方法中 .然而与分子轨道理论方法不同 ,在价键方法中 ,无法直接得到空轨道 ,此外如何选取激发价键函数使得计算结果…     

Analysis of vibrational spectra of L-alanylglycine based on density functional theory calculations

FT Raman and IR spectra of the crystallized biologically active molecule, L-alanylglycine (L-Ala-Gly) have been recorded and analyzed. The equilibrium geometry, bonding features and harmonic vibrational frequencies of L-Ala-Gly have been investigated with the help of B3LYP density functional theory (DFT) method. The calculated molecular geometry has been compared with the experimental data. The assignments of the vibrational spectra have been carried out with the help of normal coordinate analysis (NCA) following the scaled quantum mechanical force field methodology (SQMFF). The optimized geometry shows the non-planarity of the peptide group of the molecule. Potential energy surface (PES) scan studies has also been carried out by ab initio calculations with B3LYP/6-311+G** basis set. The red shifting of NH3+ stretching wavenumber indicates the formation of N-H...O hydrogen bonding. The change in electron density (ED) in the sigma* antibonding orbitals and E2 energies have been calculated by natural bond orbital analysis (NBO) using DFT method. The NBO analysis confirms the occurrence of strong intermolecular hydrogen bonding in the molecule.     

Spectroscopic analysis of small organic molecules: a comprehensive near-edge x-ray-absorption fine-structure study of C6-ring-containing molecules

We report high-resolution C 1s near-edge x-ray-absorption fine-structure (NEXAFS) spectra of the C6-ring-containing molecules benzene (C6H6), 1,3- and 1,4-cyclohexadiene (C6H8), cyclohexene (C6H10), cyclohexane (C6H12), styrene (C8H8), and ethylbenzene (C8H10) which allow us to examine the gradual development of delocalization of the corresponding pi electron systems. Due to the high experimental resolution, vibrational progressions can be partly resolved in the spectra. The experimental spectra are compared with theoretical NEXAFS spectra obtained from density-functional theory calculations where electronic final-state relaxation is accounted for. The comparison yields very good agreement between theoretical spectra and experimental results. In all cases, the spectra can be described by excitations to pi*- and sigma*-type final-state orbitals with valence character, while final-state orbitals of Rydberg character make only minor contributions. The lowest C 1s-->1pi* excitation energy is found to agree in the (experimental and theoretical) spectra of all molecules except for 1,3-cyclohexadiene (C6H8) where an energy smaller by about 0.6 eV is obtained. The theoretical analysis can explain this result by different binding properties of this molecule compared to the others.     

Quantifying the symmetry content of the electronic structure of molecules: molecular orbitals and the wave function

A scheme to quantify the symmetry content of the electronic wave function and molecular orbitals for arbitrary molecules is developed within the formalism of Continuous Symmetry Measures (CSMs). After defining the symmetry operation expectation values (SOEVs) as the key quantity to gauge the symmetry content of molecular wavefunctions, we present the working equations to be implemented in order to carry out real calculations using standard quantum chemistry software. The potentialities of a symmetry analysis using this new method are shown by means of some illustrative examples such as the changes induced in the molecular orbitals of a diatomic molecule by an electronegativity perturbation, the breaking of orbital symmetry along the dissociation path of the H(2) molecule, the changes in the molecular orbitals upon a geometrical distortion of the benzene molecule, and the inversion symmetry content in the different spin states of the [Fe(CH(3))(4)](2-) complex.     

Electronic structure of copper phthalocyanine: an experimental and theoretical study of occupied and unoccupied levels

An experimental and theoretical study of the electronic structure of copper phthalocyanine (CuPc) molecule is presented. We performed x-ray photoemission spectroscopy (XPS) and photoabsorption [x-ray absorption near-edge structure (XANES)] gas phase experiments and we compared the results with self-consistent field, density functional theory (DFT), and static-exchange theoretical calculations. In addition, ultraviolet photoelectron spectra (UPS) allowed disentangling several outer molecular orbitals. A detailed study of the two highest occupied orbitals (having a(1u) and b(1g) symmetries) is presented: the high energy resolution available for UPS measurements allowed resolving an extra feature assigned to vibrational stretching in the pyrrole rings. This observation, together with the computed DFT electron density distributions of the outer valence orbitals, suggests that the a(1u) orbital (the highest occupied molecular orbital) is mainly localized on the carbon atoms of pyrrole rings and it is doubly occupied, while the b(1g) orbital, singly occupied, is mainly localized on the Cu atom. Ab initio calculations of XPS and XANES spectra at carbon K edge of CuPc are also presented. The comparison between experiment and theory revealed that, in spite of being formally not equivalent, carbon atoms of the benzene rings experience a similar electronic environment. Carbon K-edge absorption spectra were interpreted in terms of different contributions coming from chemically shifted C 1s orbitals of the nonequivalent carbon atoms on the inner ring of the molecule formed by the sequence of CN bonds and on the benzene rings, respectively, and also in terms of different electronic distributions of the excited lowest unoccupied molecular orbital (LUMO) and LUMO+1. In particular, the degenerate LUMO appears to be mostly localized on the inner pyrrole ring.     

Vibrational effects on the electron momentum distributions of valence orbitals of formamide

The ionization energy spectra and electron momentum distributions of formamide were investigated using the high-resolution electron momentum spectrometer in combination with high level calculations. The observed ionization energy spectra and electron momentum distributions were interpreted using symmetry adapted cluster-configuration interaction theory, outer valence Green function, and DFT-B3LYP methods. The ordering of 10a(') and 2a(") orbitals of formamide was assigned unambiguously by comparing the experimental electron momentum distributions with the corresponding theoretical results, i.e., 10a(') has a lower binding energy. In addition, it was found that the low-frequency wagging vibration of the amino group at room temperature has noticeable effects on the electron momentum distributions. The equilibrium-nuclear-positions-approximation, which was widely used in electron momentum spectroscopy, is not accurate for formamide molecule. The calculations based on the thermal average can evidently improve the agreement with the experimental momentum distributions.     

Probing the shape and stereochemistry of molecular orbitals in locally flexible aromatic chains: a penning ionization electron spectroscopy and Green's function study of the electronic structure of biphenyl

We report on the results of an exhaustive study of the interplay between the valence electronic structure, the topology and reactivity of orbitals, and the molecular structure of biphenyl by means of Penning ionization electron spectroscopy in the gas phase upon collision with metastable He*(2(3)S) atoms. The measurements are compared with one-particle Green's function calculations of one-electron and shake-up valence ionization spectra employing the third-order algebraic diagrammatic construction scheme [ADC(3)]. Penning ionization intensities are also analyzed by means of the exterior electron-density model and comparison with photoelectron spectra: in contrast with the lines originating from sigma orbitals, ionization lines belonging to the pi-band system have large Penning ionization cross sections due to their greater extent outside the molecular van der Waals surface. The involved chemi-ionization processes are further experimentally investigated using collision-energy-resolved Penning ionization electron spectroscopy. The cross sections of pi-ionization bands exhibit a markedly negative collision-energy dependence and indicate that the interaction potential that prevails between the molecule and the He*(2(3)S) atom is strongly attractive in the pi-orbital region. On the other hand, the partial ionization cross sections pertaining to sigma-ionization channels are characterized by more limited collision-energy dependencies, as a consequence of rather repulsive interactions within the sigma-orbital region. A comparison of ADC(3) simulations with the Penning ionization electron spectra and UV photoelectron spectra measured by Kubota et al. [Chem. Phys. Lett. 1980, 74, 409] on thin films of biphenyl deposited at 170 and 109 K on copper demonstrates that biphenyl molecules lying at the surface of polycrystalline layers adopt predominantly a planar configuration, whereas within an amorphous sample most molecules have twisted structures similar to those prevailing in the gas phase.     

Photoelectron spectroscopy of the doubly-charged anions [MIVO(mnt)2]2- (M = Mo, W; mnt = S2C2(CN)2(2-): access to the ground and excited states of the [MVO(mnt)2]- anion

Photodetachment photoelectron spectroscopy was used to investigate the electronic structure of the doubly charged complexes [MIVO(mnt)2]2- (M = Mo, W; mnt = 1,2-dicyanoethenedithiolato). These dianions are stable in the gas phase and are minimal models for the active sites of the dimethyl sulfoxide reductase family of molybdenum enzymes and of related tungsten enzymes. Adiabatic and vertical electron binding energies for both species were measured, providing detailed information about molecular orbital energy levels of the parent dianions as well as the ground and excited states of the product anions [MVO(mnt)2]-. Density functional theory calculations were used to assist assignment of the detachment features. Differences in energy between these features provided the energies of ligand-to-metal charge-transfer transitions from S(pi) and S(sigma) molecular orbitals to the singly occupied metal-based orbital of the products [MVO(mnt)2]-. These unique data for the M(V) species were obtained at the C(2)(v)() geometry of the parent M(IV) dianions. However, theoretical calculations and available condensed phase data suggested that a geometry featuring differentially folded dithiolene ligands (Cs point symmetry) was slightly lower in energy. The driving force for ligand folding is a favorable covalent interaction between the singly occupied metal-based molecular orbital (a1 in C2v) point symmetry; highest occupied molecular orbital (HOMO)) and the least stable of the occupied sulfur-based molecular orbitals (b1 in C2v point symmetry, HOMO-1) that is only possible upon reduction to the lower symmetry. This ligand folding induces a large increase in the intensity predicted for the a' S(pi) --> a' dx2 - y2 charge-transfer transition originating from the HOMO-2 of [MVO(mnt)2](-) under Cs point symmetry. Electronic absorption spectra are available for the related species [MoVO(bdt)2]- (bdt = 1,2-benzenedithiolato) and for the oxidized form of dimethyl sulfoxide reductase. The intense absorptions at approximately 1.7 eV have been assigned previously to S(sigma) --> Mo transitions, assuming C2v geometry. The present work indicates that the alternative a' S(pi) --> a' dx2 - y2 of Cs geometry must be considered. Overall, this study confirms that the electronic structure of the M-dithiolene units are exquisitely sensitive to dithiolene ligand folding, reinforcing the proposal that these units are tunable conduits for electron transfer in enzyme systems.     

Structure and electronic properties of a benzene-water solution

Electronic properties of benzene in water were investigated by a sequential quantum mechanical/molecular dynamics approach. Emphasis was placed on the analysis of the structure, polarization effects, and ionization spectrum. By adopting a polarizable model for both benzene and water the structure of the benzene-water solution is in good agreement with data from first principles molecular dynamics. Further, strong evidence that water molecules acquire enhanced orientational order near the benzene molecule is found. Upon hydration, the quadrupole moment of benzene is not significantly changed in comparison with the gas-phase value. We are also reporting results for the dynamic polarizability of benzene in water. Our results indicate that the low energy behaviour of the dynamic polarizability of gas-phase and hydrated benzene is quite similar. Outer valence Green's function calculations for benzene in liquid water show a splitting of the gas-phase energy levels associated with the 1e(1g)(π), 2e(2g), and 2e(1u) orbitals upon hydration. Lifting of the orbitals degeneracy and redshift of the outer valence bands is related to symmetry breaking of the benzene structure in solution and polarization effects from the surrounding water molecules.     

Sigma aromaticity of the bimetallic Au5Zn+ cluster

Experimental and theoretical evidence for a "sigma aromatic" bimetallic cluster is presented. A mass spectrometric analysis of AuNZn+ (N = 2-44) photofragments shows Au5Zn+ to be very abundant, proving its high stability. Calculations predict that Au5Zn+ has a planar geometry and six valence s electrons occupying delocalized sigma-bonded molecular orbitals in a manner similar to that of aromatic organic compounds, except for their nodal properties in the molecular plane. The delocalized sigma electrons provide diamagnetic ring currents, suggesting Au5Zn+ is aromatic.     

A new fragment-based approach for calculating electronic excitation energies of large systems

We present a new fragment-based scheme to calculate the excited states of large systems without necessity of a Hartree-Fock (HF) solution of the whole system. This method is based on the implementation of the renormalized excitonic method [M. A. Hajj et al., Phys. Rev. B 72, 224412 (2005)] at ab initio level, which assumes that the excitation of the whole system can be expressed by a linear combination of various local excitations. We decomposed the whole system into several blocks and then constructed the effective Hamiltonians for the intra- and inter-block interactions with block canonical molecular orbitals instead of widely used localized molecular orbitals. Accordingly, we avoided the prerequisite HF solution and the localization procedure of the molecular orbitals in the popular local correlation methods. Test calculations were implemented for hydrogen molecule chains at the full configuration interaction, symmetry adapted cluster/symmetry adapted cluster configuration interaction, HF/configuration interaction singles (CIS) levels and more realistic polyene systems at the HF/CIS level. The calculated vertical excitation energies for lowest excited states are in reasonable accordance with those determined by the calculations of the whole systems with traditional methods, showing that our new fragment-based method can give good estimates for low-lying energy spectra of both weak and moderate interaction systems with economic computational costs.     

Ab-initio calculations on large molecules and solids by desirable computational procedures

Desirable computational procedures developed here recently for ab-initio calculations on large molecules are outlined. Effective core model potentials (MODPOT) permit calculations of valence electrons only explicitly, yet accurately; a charge-conserving integral prescreening evaluation to decide whether a block of integrals will be larger than a preset threshold and thus be calculated explicitly is effective for spatially extended systems; an efficient MERGE technique to save and reuse common invariant skeletal integrals is useful for geometry variations and for adding basis fcuntions, substituent groups and molecules; and an effective configuration interaction (CI) Hamiltonian into which are folded the effects of the occupied molecular orbitals from which no excitations are allowed is useful for molecular decompositions and intermolecular reactions. These techniques have been extended for CI calculations on breaking a chemical bond in a molecule in a crystal or solid; atom-class/atomic-class potential functions and dispersion calculations have been added. In a new program, POLY-CRYST, all the integral strategies for large molecules are meshed.     

Nature of coordination bond in silatranes and its formation dynamics according to the ab initio calculations

Quantum-chemical calculations of 1-hydrosilatrane molecule with complete optimization of its geometry and at various fixed Si…N distances (2.0 to 3.7 Å) has been carried out at the MP2/6-31G(d) level. The silatranes coordination bond is formed of different atomic orbitals of Si and N atoms participating in a series of molecular orbitals. With the Si…N distance decreasing, contributions of the atomic orbitals in these molecular orbitals have been changed, number of the molecular orbitals has increased, and total energy of the molecule has decreased. At the coordination centers are getting closer, population of the nitrogen valence s and p _z orbitals have changed due to the corresponding bond angle change; the populations of Si and H orbitals are not significantly changed.     

Experimental and theoretical evidence for the formation of zinc tricarbonyl in solid argon

Zinc carbonyls are extremely rare. Here we report experimental and theoretical evidence of unprecedented zinc tricarbonyl, Zn(CO)3, the next member of the series of 18-electron metal carbonyls Cr(CO)6 --> Fe(CO)5 --> Ni(CO)4, whereas there is no evidence for the formation of the zinc mono- and dicarbonyls Zn(CO)n (n = 1, 2). DFT calculations predict that the Zn(CO)3 molecule has a singlet ground state with D3h symmetry. The formation of Zn(CO)3 involves 4s --> 4p promotion of the Zn atom, which increases the Zn-CO bonding by decreasing the sigma repulsion and significantly increasing the Zn 4sp hybrid orbitals --> CO pi* back-donation.     

大体系多电子相关研究中应用群对称定域轨道的构想

大体系多电子相关研究中应用群对称定域轨道的构想周泰锦，刘爱民（厦门大学化学系，厦门３６１００５）关键词：组态相关，多构型自治叠代，多中心积分，群对称定域轨道，对称约化有关原子簇化合物及化学吸附、过渡态、激发态、催化反应等大体系的量子化学研究，对于探讨...     

Attachment of different donor groups to a cryptand for modulation of two-photon absorption cross-section

Two-photon absorption (TPA) properties of a laterally nonsymmetric aza cryptand with attached side arms have been investigated. This series of Schiff base derivatives supports the mechanistic approach for enhancing the TPA process, which is usually dictated by molecular geometry, pi-bridging, delocalization length, and corresponding charge-transfer possibilities. The results described here suggest that on increasing the branching units, the TPA cross-section, sigma((2)), can be tuned to a larger value. The TPA activity is "switched on" when a metal atom enters the cavity and serves as a conduit of electronic delocalization. The sigma((2)) value increases as the donor strength increases. The maximum value is obtained on moving from the single-branched system to the nearly threefold symmetry. This serves as a useful synthetic strategy for designing novel octupolar molecules with high sigma((2)) values. Theoretical calculations at the B3LYP functional with the 6-31G* basis set under DFT formalism provide supporting evidence that the communication between the side arms through the metal d orbital and more ordered geometry of chromophores leads to a smaller HOMO-LUMO gap, which has a great influence upon the electronic properties of the molecules.     

The time scale for electronic reorganization upon sudden ionization of the water and water-methanol hydrogen bonded dimers and of the weakly bound NO dimer

When the valence molecular orbital is localized sudden ionization can cause the nascent hole to move rapidly even before any relaxation of the geometry occurs. Hydrogen bonded clusters offer suitable test systems where the hole is initially localized on one moiety. Computational studies are reported for the water dimer and water-methanol bimer. The local ionization potential of water is different in the methanol-water and water-methanol conformers and this difference is very clearly reflected in the dynamics of charge migration. For the NO dimer the results are that its structure is symmetric so that the two NO molecules are equivalent and do not exhibit the required localization. The role of symmetry is also evident in the charge propagation for holes created in different orbitals. Localization of the initial hole distribution even if absent in the bare molecule can still be induced by the intense electric field of a sudden photoionization. This effect is computationally studied for the NO dimer in the presence of a static electric field.     

Core-shell photoabsorption and photoelectron spectra of gas-phase pentacene: experiment and theory

The C K-edge photoabsorption and 1s core-level photoemission of pentacene (C22H14) free molecules are experimentally measured, and calculated by self-consistent-field and static-exchange approximation ab initio methods. Six nonequivalent C atoms present in the molecule contribute to the C 1s photoemission spectrum. The complex near-edge structures of the carbon K-edge absorption spectrum present two main groups of discrete transitions between 283 and 288 eV photon energy, due to absorption to pi* virtual orbitals, and broader structures at higher energy, involving sigma* virtual orbitals. The sharp absorption structures to the pi* empty orbitals lay well below the thresholds for the C 1s ionizations, caused by strong excitonic and localization effects. We can definitely explain the C K-edge absorption spectrum as due to both final (virtual) and initial (core) orbital effects, mainly involving excitations to the two lowest-unoccupied molecular orbitals of pi* symmetry, from the six chemically shifted C 1s core orbitals.     

The structure of linear SiCC: An ab initio SCF CI study including vibrational effects

Large-scale SCF CI calculations have been performed for the ground ¹Σ⁺ state of linear SiCC. The calculation includes up to quadruple excitations in the valence space plus all single and double excitations from the valence localized orbitals of the single HF configuration. Vibrational wavefunctions have been derived from the CI potential surface. Vibrational frequencies including anharmonicity corrections are calculated together with the zero-point vibrational correction to the rotational constant. The large dipole moment, 4.62 D, should make this molecule suitable for radioastronomic searches.