Counterpoise技术和He-H2O配合物的平衡几何

在得到MP4, SCF的关于总能量及结合能的四个势能面的基础上,着重分析了电子相关作用, 消除基函数超位误差的Counterpoise(CP)技术对van der Waals配合物He-H2O的平衡几何的影响并得到了更精确的平衡几何预测。     

Structural changes in the molecular sheets along (hk0) planes derived from cellulose Iβ by molecular dynamics simulations

We studied the stability of molecular sheets with four cellotetraoses in an aqueous environment by molecular dynamics simulation to identify the molecular details of first structure as one of the possibilities in the course of crystallization of cellulose I. After simulation, the molecular sheets formed by van der Waals forces along the (11?0) and (110) crystal plane did not change their structures in an aqueous environment, whereas the other ones formed by hydrogen bonds along the (100) and (200) crystal plane changed into a van der Waals associated molecular sheet, similar to the former. These simulated molecular sheets formed by van der Waals forces were structurally stable in water because of their hydrophilic exterior and hydrophobic interior. Therefore, if the molecular sheet structures are formed in the real system, the sheets formed by van der Waals forces are probably the initial structure of crystallization. A close analysis indicated that these sheets could be classified into two groups in terms of the hydrogen bonding networks, camber angle, and main and side chain conformations. One group was the molecular sheets corresponding to the (110) after simulation. This sheet is probably rigid because intramolecular hydrogen bonds of the chains in the sheet are highly developed. The other group was the molecular sheets corresponding to (200), (100), and (11?0) crystal plane: the chains in these sheets seemed to be rather flexible due to their moderately developed intramolecular hydrogen bonds.     

An ab initio investigation of the O(3P)-H2(1sigma(g)+) van der Waals well

We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well. We predict the van der Waals minimum to occur in perpendicular geometry, and located at a closer distance than a secondary well in colinear geometry. The potentials obtained in the present study confirm the previous calculations of Alexander [M. H. Alexander, J. Chem. Phys., 1998, 108, 4467], but disagree with the earlier work of Harding and co-workers [Z. Li, V. A. Apkarian and L. B. Harding, J. Chem. Phys., 1997, 106, 942] as well as with recently refitted surfaces of Brand?o and coworkers [J. Brand?o, C. Mogo and B. C. Silva, J. Chem. Phys., 2004, 121, 8861]. Inclusion of spin-orbit coupling reduces the depth of the van der Waals minimum without causing a change in its position.     

Rotational spectroscopic and ab initio studies of the Xe-H2O van der Waals dimer

An ab initio potential energy surface of the Xe-H(2)O van der Waals dimer was constructed at the coupled cluster level of theory with single, double, and pertubatively included triple excitations. For the Xe atom, the small-core pseudopotential and augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ-PP) basis set was used. Dunning's augmented correlation-consistent polarized valence triple-zeta (aug-cc-pVTZ) basis set was chosen for O and H atoms. Midbond functions were used to supplement the atom-centered basis sets. Rotational spectra of the Xe-H(2)O van der Waals dimer were recorded with a pulsed-nozzle Fourier transform microwave spectrometer. Rotational transitions within two internal rotor states, namely, the 0(00) and 1(01) states, were measured and assigned. Nuclear quadrupole hyperfine structures due to the (131)Xe (I = (3)/(2)), D (I = 1) and (17)O (I = (5)/(2)) nuclei were also observed and analyzed. Information about the molecular structure and the H(2)O angular motions was extracted from the spectroscopic results with the assistance of the ab initio potential.     

X…H2O（X=Li，Na，K）非线性光学性质的从头算理论

在MP2水平上采用6-311G基组计算了van der Waals复合物X…H2O(X=Li, Na, K)的非线性光学性质(μ, α, β), 讨论了基组效应和电子相关效应对计算结果的贡献. 在MP2/6-311++G(2df, 2pd)水平上计算得到的三个复合物分子X(X=Li, Na, K)•••H2O的非线性光学性质. 结果表明, 三种复合物分子均具有巨大的一阶超极化率, 其中最外层电子的弥散特性对一阶超极化率有很大的影响.     

Monte Carlo studies of isomers, structures, and properties in benzene-cyclohexane clusters: computation strategy and application to the dimer and trimer, (C6H6)(C6H12)n, n = 1-2

Low-temperature isomeric energies, structures, and properties of benzene-cyclohexane clusters are investigated via Monte Carlo simulations. The Monte Carlo strategy is first documented and then applied to (C(6)H(6))(C(6)H(12)) and (C(6)H(6))(C(6)H(12))(2) using four different potential energy surfaces. Results identify a single parallel-displaced dimer isomer. MP2 optimizations and frequency calculations support the Monte Carlo dimer structure and identify the van der Waals mode observed in vibronic spectra. Caloric simulations identify two temperatures where structural transitions occur and imply an experimental temperature below 10 K for dimers in cold supersonic expansions. The (C(6)H(6))(C(6)H(12))(2) studies identify eight independent trimer isomers: three form parallel-stacked (sandwich) arrangements with the two cyclohexane moieties related through a D(6)(h) transformation. The remaining five trimer isomers are trigonal, with no overall symmetry. Caloric studies indicate that the sandwich and trigonal isomeric classes coexist independently below 60 K, consistent with trimer vibronic spectra that contain two independent van der Waals progressions.     

Intramolecular chalcogen-tin interactions in (o-MeE-C6H4)CH2SnPh3-nCln (E = S, O; n = 0, 1, 2), characterized by X-ray diffraction and 119Sn solution and solid-state NMR

Organotin(IV) compounds of the type (o-MeE-C6H4)CH2SnPh3-nCln were synthesized, E = O, n = 0 (1), n = 1 (2), n = 2 (3) and E = S, n = 0 (4), n = 1 (5), n = 2 (6). The complexes exhibit significant trigonal bipyramidal pentacoordination at tin as a consequence of intramolecular Sn-O (1-3) and Sn-S (4-6) interactions upon substitution of the phenyl groups by chloro groups. The intramolecular Sn-O distances in 1, 2, and 3 are 83%, 75%, and 79% of the sum of the van der Waals radii. The equivalent Sn-S values for 4, 5, and 6 are 90%, 73%, and 71%, respectively. The geometry of compound 3 is complicated by intermolecular dimerization via bridging chlorines creating a distorted octahedral geometry at tin. The related dichloro sulfur compound 6 also exhibits an intermolecular association in the form of Sn-Cl-H hydrogen bonding leading to a polymeric structure in the solid state. CPMAS 119Sn NMR spectroscopy suggests that the intramolecular Sn-E interactions persist in solution and also facilitated the discovery of a new crystalline form of 4, 4', that contains a Sn-S distance which is 95% the sum of the van der Waals radii.     

A spectroscopically effective molecular mechanics model for the intermolecular interactions of the hydrogen-bonded N-methylacetamide dimer

An MP2/6-31+G¹ calculation of the N-methylacetamide dimer shows that it has two minimum energy structures, both hydrogen bonded with peptide planes roughly perpendicular to each other. A complete molecular mechanics optimization of the dimer has been done, using a model for the intermolecular interactions consisting of charges, atomic dipoles, and van der Waals interactions and the methodology of our spectroscopically determined force field for the intramolecular interactions. The two structures are satisfactorily reproduced, as are their interaction energies, their dipole moments, and, from the point of view of our goal of a spectroscopically accurate force field, their six intermolecular normal mode frequencies.     

The Weak Interaction on CO2 …CO van der Waals Complex : A Theoretical Study

The geometries of van der Waals complex CO2…CO were optimized at DFT and second-order Moller-Plesset perturbation(MP2) levels with the large basis set,three stable structures were found.The most stable structure has a T-shape geometry in which the CO lies along the C2 axis of CO2,with the two C atoms direct contact and R(C…C)=0.3227nm.The corresponding energies of the most stable structure were calculated by means of MP2,MP4D,MP4DQ,MP4SDTQ,MP4SDQ,CCSD and CCSD(T) methods,The BSSE (basis set superposition error) wads eliminated by the Boys-Bernardi counterpoise correction(CP) method.According to thermodynamics data.van der Waals complex CO2…CO can found at a low temperature and or a high pressure,There is a little charge transferred between the two interacted subunits.In the most stable structure,CO2 is the acceptor and CO is the donor.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

Ab initio calculations of H2-H2 potential surfaces near the van der Waals minimum

The interaction energy between two hydrogen molecules near the van der Waals minimum is computed, for four relative orientations, and for intermolecular distances ranging from 5.5 to 16 au. A partially optimized basis set limited to 52 independent gaussian functions was used throughout the energy calculations. A new method based on the polarized atomic orbital technique has been used to reduce subsequently the size of the CI calculations which makes this method tractable for heavier molecular systems than H₂-H₂.     

Ab initio calculation of the NH(3sigma-)-NH(3sigma-) interaction potentials in the quintet, triplet, and singlet states

We present the ab initio potential-energy surfaces of the NH-NH complex that correlate with two NH molecules in their 3sigma- electronic ground state. Three distinct potential-energy surfaces, split by exchange interactions, correspond to the coupling of the S(A) = 1 and S(B) = 1 electronic spins of the monomers to dimer states with S = 0, 1, and 2. Exploratory calculations on the quintet (S = 2), triplet (S = 1), and singlet (S = 0) states and their exchange splittings were performed with the valence bond self-consistent-field method that explicitly accounts for the nonorthogonality of the orbitals on different monomers. The potential surface of the quintet state, which can be described by a single Slater determinant reference function, was calculated at the coupled cluster level with single and double excitations and noniterative treatment of the triples. The triplet and singlet states require multiconfiguration reference wave functions and the exchange splittings between the three potential surfaces were calculated with the complete active space self-consistent-field method supplemented with perturbative configuration interaction calculations of second and third orders. Full potential-energy surfaces were computed as a function of the four intermolecular Jacobi coordinates, with an aug-cc-pVTZ basis on the N and H atoms and bond functions at the midpoint of the intermolecular vector R. An analytical representation of these potentials was given by expanding their dependence on the molecular orientations in coupled spherical harmonics, and representing the dependence of the expansion coefficients on the intermolecular distance R by the reproducing kernel Hilbert space method. The quintet surface has a van der Waals minimum of depth D(e) = 675 cm(-1) at R(e) = 6.6a0 for a linear geometry with the two NH electric dipoles aligned. The singlet and triplet surfaces show similar, slightly deeper, van der Waals wells, but when R is decreased the weakly bound NH dimer with S = 0 and S = 1 converts into the chemically bound N2H2 diimide (also called diazene) molecule with only a small energy barrier to overcome.     

The mechanism of H-bond rupture: the vibrational pre-dissociation of C2H2-HCl and C2H2-DCl

Pair correlated fragment rovibrational distributions are presented following vibrational predissociation of the C2H2-DCl van der Waals dimer initiated by excitation of the asymmetric (asym) C-H stretch. The only observed fragmentation pathways are DCl (v= 0; j= 6-9)+ C2H2(nu2= 1; j= 1-5). These and previously reported data on the related C2H2-HCl species are analysed using the angular momentum (AM) method. Calculations accurately reproduce fragment rovibrational distributions following dissociation of the C2H2-HCl dimer initiated either by excitation of the asym C-H stretch or via the HCl stretch, and those from C2H2-DCl initiated via asym C-H stretch excitation. The calculations demonstrate that the dimer is bent at the moment of dissociation. Several geometries are found that lead to H-bond breakage via a clearly identified set of fragment quantum states. The results suggest a hierarchy in the disposal of excess energy and angular momentum between fragment vibration, rotation and recoil. Deposition of the largest portion of energy into a C2H2 vibrational state sets an upper limit on HCl rotation, which then determines the energy and AM remaining for C2H2 rotation and fragment recoil. Acceptor C2H2 vibrational modes follow a previously noted propensity, implying that the dissociating impulse must be able to induce appropriate nuclear motions both in the acceptor vibration and in rotation of the C2H2 fragment.     

Electronic structure and chemical bond in naphthalene and anthracene

We investigated the electronic structure of crystalline naphthalene and anthracene within the framework of density functional theory including van der Waals interactions (DFT-D). It is established that for better agreement with experimental values it is necessary to use the increased values of the van der Waals radii, which is caused by an overestimated value of the van der Waals interactions in crystalline linear oligoacenes. Utilization of the DFT-D leads to a correct account of the dispersion forces, which results in a high precision of the computed lattice parameters and cohesive energy. Based on the relaxed crystal structures, we have computed the total and deformation electron density and determined the mechanism of chemical bonds formation in crystals of naphthalene and anthracene. It has been established that the chemical bond in molecular crystals is formed under the influence of not only intramolecular but also intermolecular interactions. On the basis of the Mulliken population analysis it was revealed that two C(3) atoms in naphthalene (or C(3) and C(4) in anthracene) have a positive charge and the population of the rest of the carbon atoms increased, as compared with isolated molecule.     

MBPT studies of van der Waals molecules. III. The reliability of apparently accurate calculations for the magnesium dimer

The interaction potential for the magnesium dimer is calculated by using the diagrammatic many-body perturbation theory within the framework of the supermolecule approach. Different approximations for the perturbation treatment of the electron correlation effects are analysed with particular emphasis on the results of the complete fourth-order many-body perturbation theory. The influence of the composition of the basis set on the calculated interaction potentials and the basis set superposition effects are investigated. It is concluded that the interaction potentials for van der Waals systems calculated by the supermolecule technique might be highly uncertain because of the basis set superposition effects at the correlated level. The advantages and disadvantages of different methods for the calculation of the energy of weakly interacting systems are discussed in the light of their ability to cope with both the electron correlation problem and the basis set superposition effect.     

邻二甲苯…Ar, N2, NH3(ND3)范德华复合物的共振双光子电离光谱

在超声分子束中，使用双光子共振电离光谱技术和飞行时间质谱技术研究了复合物邻二甲苯…Ar．N2，NH3(ND3)．通过理论计算及同位素光谱效应．合理地归属了这些复合物的光谱．并由此获得这些复合物分子问各种模式的振动频率．     

Calculations on the vertical and adiabatic ionization energies of (H2S)2

After determining reliable procedures for calculating the ionization energy of H₂S using many-body perturbation theory and electron propagator theory, the same procedures are used to calculate the vertical ionization energy of the van der Waals molecule (H₂S)₂. The adiabatic ionization energy calculated for a dimer cation with a H₃S⁺...SH structure is in satisfactory agreement with photoionization threshold experiments. An alternative dimer cation structure with a three-electron S-S bond is also discussed.     

Interaction-induced dipole moment of the Ar-H2 dimer: dependence on the H2 bond length

We present ab initio calculations of the interaction-induced dipole moment of the Ar-H2 van der Waals dimer. The primary focus of our calculations is on the H2 bond length dependence of the dipole moment, which determines the intensities of both the collision-induced H2 upsilon = 1 <-- 0 fundamental band in gaseous Ar-H2 mixtures and the dopant-induced H2 upsilon = 1 <-- 0 absorption feature in Ar-doped solid H2 matrices. Our calculations employ large atom-centered basis sets, diffuse bond functions positioned between the two monomers, and a coupled cluster treatment of valence electron correlation; core-valence correlation effects appear to make negligible contributions to the interaction-induced dipole moment for the Ar-H2 configurations considered here.     

Van der Waals surface graphs and molecular shape

A van der Waals surface graph is the graph defined on a van der Waals surface by the intersections of the atomic van der Waals spheres. A van der Waals shape graph has a vertex for each atom with a visible face on the van der Waals surface, and edges between vertices representing atoms with adjacent faces on the van der Waals surface. These are discrete invariants of three‐dimensional molecular shape. Some basic properties of van der Waals surface graphs are studied, including their relationship with the Voronoi diagram of the atom centres, and a class of molecular embeddings is identified for which the dual of the van der Waals surface graph coincides with the van der Waals shape graph. This revised version was published online in July 2006 with corrections to the Cover Date.