Determination of accurate semiexperimental equilibrium structure of proline using efficient transformations of anharmonic force fields among the series of isotopologues

The complete semiexperimental (SE) equilibrium structure of proline (45 degrees of freedom) is determined using the mixed estimation method. The cubic force fields for the parent and eight isotopologues of proline molecule are evaluated at the MP2-FC/cc-pVTZ level in Cartesian coordinates. The accuracy of the SE structure is verified by optimising the structure with the CCSD(T) model and a basis set of quadruple-ζ quality. A significantly more accurate equilibrium structure of proline is obtained when compared to the previous one. It is shown that the employed technique is efficient for the determination of SE equilibrium structures of rather large molecules. A simple transformation of anharmonic force fields between normal coordinate and Cartesian coordinate representations is proposed. The suggested technique allows efficient evaluation of the rotation–vibration interaction constants for a number of isotopologues, once the cubic force field of any species is found either in normal or Cartesian coordinates.     

Making More Extensive Use of the Coupled-cluster Wave Function: from the Standard Energy Expression to the Energy Expectation Value

The standard coupled-cluster (CC) scheme with single and double excitations in the cluster operator (CCSD) includes only up to quadruple excitations in the equations. The CCSD exponential expansion generates, however, all possible excitations out of the reference function through products of the cluster operators. Clearly, in all standard approximate CC approaches only a part of the CC wave function is used in the equations. If the standard CCSD wave function is inserted into the energy expectation value expression then the complete CCSD wave function contributes to the energy. Such an energy expectation value expression can be presented as a sum of the standard CCSD energy formula plus correction terms. The correction terms provide an information about the quality of the total CC function. Contributions associated with the presence of higher than double excitations in the bra CCSD wave function supplement the CCSD energy obtained within the standard scheme. These contributions can be generated in a sequential way by considering intermediate excitation levels for the bra CCSD wave function in the expectation value expression before reaching the highest excitation level. In this way the importance of particular components differing in the standard and expectation value CCSD energies can be investigated. Some of the contributions can be recognized as close to or identical with the so-called renormalized noniterative corrections to the CC methods. We try to see to what an extent the nonstandard energy expressions, like the energy expectation value or the asymmetric energy formula, can be used to extend the applicability of the CCSD method illustrating our considerations with some numerical examples.Dedicated to Professor Jean-Paul Malrieu to honor his contribution to quantum chemistry and physics     

The coupling constant averaged exchange–correlation energy density

ABSTRACT

The exchange–correlation energy, central to density-functional theory, may be represented in terms of the coupling constant averaged (CCA) exchange–correlation energy density. We present an approach to calculate the CCA energy density using accurate ab initio methods and its application to simple atomic systems. This function provides a link between intrinsically non-local, many-body electronic structure methods and simple local and semi-local density-functional approximations (DFAs). The CCA energy density is resolved into separate exchange and correlation terms and the features of each compared with those of quantities commonly used to construct DFAs. In particular, the more complex structure of the correlation energy density is found to exhibit features that align well with those present in the Laplacian of the density, suggesting its role as a key variable to be used in the construction of improved semi-local correlation functionals. The accurate results presented in this work are also compared with those provided by the Laplacian-dependent Becke–Roussel model for the exchange energy.     

Thied—Order Approximation of 0^＋＋ Glueball Mass and Wavefunction of （2＋1）—Dimensional SU（3） Lattice Gauge Theory

The random phase approximation is applied to the coupled-cluster expansions of lattice gauge theory (LGT) Using this method,wavefuctions are approximated by linear combination of graplh consisting of only one connected wilson loop.We study the excited state energy and wavelunction in (2 1)-D SU(3) LGT up to the thired order,The glueball mass shows a good scaling behavior.     

An alternative definition of the electron propagator in the superoperator form and its relation to linear response theory in a coupled-cluster framework

In this paper it is shown that (i) there exists an alternative definition of the superoperator resolvent for calculation of difference energy satisfying linked cluster theorem for a coupled-cluster choice of the ground-state function which may even be approximate; (ii) the pole-structure of this propagator-like function in superoperator form is shown to contain information similar to that contained in the conventional propagator. (iii) It is demonstrated that suitable “Killer conditions” and completeness of the “operator manifold”—essential for understanding the pole-structure of the propagator—can be established both for an exact and an approximate ground state function in a coupled-cluster form. (iv) It is also demonstrated that difference energies calculated with these propagator-like functions are identical to those obtained from a linear response theory in a coupled-cluster form put forward recently by Mukherjeeet al and Monkhorst.     

硫族铅化物阳离子低电子态的理论研究

采用我们最近发展的含旋轨耦合的运动方程耦合簇计算电离能(EOMIP-CC)方法,在CCSD级别上计算了硫族铅化物PbS、PbSe、PbTe阳离子低电子态的平衡键长和谐振频率以及绝热和垂直电离能,得到的结果与已有的实验值吻合较好.不考虑旋轨耦合(SOC)的情况下通过与CCSD(T)的计算结果比较,考察了三重激发对计算结果的影响,结果显示考虑三重激发的贡献后得到的键长和频率结果与实验值吻合更好.计算结果表明PbTe+中2Π态的能量分裂明显大于PbS+和PbSe+中2Π态的能量分裂,但是PbTe+中2Π1/2和2Σ1/2态之间的相互耦合则明显弱于PbS+和PbSe+中这两个态之间的耦合.PbTe+中2Π1/2和2Σ1/2态之间耦合很弱,一方面是因为2Σ+态和2Π态的能量差比PbS+和PbSe+中2Σ+态和2Π态的能量差大,另一方面还由于PbTe+中2Π1/2和2Σ1/2态之间的旋轨耦合矩阵元只是PbS+和PbSe+中2Π1/2和2Σ1/2态之间的旋轨耦合矩阵元的一半.这些计算结果为PbS+、PbSe+、PbTe+阳离子的低电子态性质提供了新的理论数据,可以为将来的实验数据提供参考.     

AuX (X=O,S)低电子态的理论研究

通常要用多参考态方法才能合理处理需考虑旋轨耦合（SOC）效应的开壳层分子如AuO和AuS的低电子态. 事实上,通过选取合适的参考态,采用运动方程耦合簇方法（EOM-CC）也能计算这些分子的一些低电子态,而且EOM-CC方法是单参考态方法,使用起来比多参考态方法更加简单. 本文采用最近发展的含旋轨耦合的EOM-CC 计算电离能的方法（EOMIP-CC）,选取对应的负离子为参考态,在CCSD 级别上计算了AuO 和AuS低电子态的性质. 在不考虑旋轨耦合时,通过比较EOMIP-CCSD和EOMIP-CCSDT的结果考察EOMIPCCSD的精度. 此外,与EOMIP-CCSDT结果相比,如果自旋污染较为显著而且T₁的模较大时,UCCSD（T）方法对能量最低的某一特定对称性的电子态的所对应的电离能误差约为0.1-0.15 eV. 在考虑了旋轨耦合效应后,我们的方法得到的键长和振动频率与实验值吻合较好. 另一方面,虽然EOMIP-SOC-CCSD高估了能量较高的²Δ_3/2态、²Σ_1/2⁺态和²Π_1/2态的能量,但是对于其它能量更低的电子态,它们的能量与已有实验值误差在0.2 eV 左右. 这显示我们所用的含SOC的EOMIP-CCSD方法对原本需要用多参考态方法才能处理的AuO和AuS低电子态能给出可靠的结果.     

Molecular structure of 5-fluorouracil from gas-phase electron diffraction data and quantum-chemical calculations

1. Download : Download high-res image (125KB)
2. Download : Download full-size image