Charge calculations in molecular mechanics. Part 8 Partial atomic charges from classical calculations

Summary The CHARGE2 programme, which involves the classical calculation of both the inductive and resonance contributions to the partial atomic charges in molecules is described, and the charges and electrostatic potentials obtained presented for some illustrative examples.In substituted methanes (CH₃X, CF₃X, CCl₃X) the effects of varying the electronegativity of the substituents and the - and -substituent contributions are clearly illustrated for a variety of substituent groups X.The problems involved in the inclusion of silicon into this scheme are detailed, together with the methods of overcoming them. The partial atomic charges ( and contributions) and electrostatic potentials for some silicon oxygen compounds are presented and discussed.The partial atomic charges from CHARGE2 for all the natural amino acids as their N-acetyl, N-methyl-amides are given and compared with those obtained from the AMBER and ECEPP/2 force fields. Considerable differences in these figures are observed, with the AMBER charges consistently much larger than those from the other two methods.The CHARGE2 partial atomic charges and electrostatic potentials for the four common nucleic acids, adenine, cytosine, guanine and thymine, are given and compared with those derived from other calculations. Again there is general similarity but also there are considerable differences, with those from the AMBER force field somewhat larger than the other methods.For previous parts in this series, see Refs. 1-7.     

Evaluation of the limit of performance of an analytical method based on a statistical calculation of its critical concentrations according to ISO standard 11843: Application to routine control of banned veterinary drug residues in food according to European Decision 657/2002/EC

Traceability of the measurement of analytical parameters capable of evaluating the performance of methods is an important concept for the assessment of quality for the routine control, especially for residue monitoring of non-authorized medicinal substances in food from animal origin. The European Decision no. 657/2002/EC recommends to calculate two statistical limits, CCα and CCβ, which allow to evaluate the critical concentrations above which the method reliably distinguish and quantify a substance taking into account the variability of the method and the statistical risk to take a wrong decision. The calculation, which can be derived from the ISO standard no. 11843 is applied on a routine basis. An example displays a very simple way for evaluating the performance of an LC-MSMS method which has been validated a few years ago and is qualified onto a Micromass Quattro LCZ tandem mass spectrometer to monitor and confirm the nitrofuran metabolite residues in food from animal origin. Community Reference Laboratory for Antimicrobial Veterinary Drug Residue Control in Food from Animal Origin     

Theoretical studies on the reaction path and dynamics of the reaction CH_3+H_2O→CH_4+OH

The reaction path, the dynamical properties along the reaction path and CVT rate constants are computed by the ab initio MO method, the reaction path Hamiltonian theory and the variational transition state theory. The results show that the effect of the electron correlation energy on activation barrier is large, the recrossing and tunneling effects exist in the reaction.     

Molecular pair potential of chiral amino acid amphiphile in Langmuir monolayers on the basis of an atomistic model

In the present work, the pair potential of enantiomeric N-palmitoyl aspartic acid amphiphile monolayer at the air/water interface is calculated based on an atomistic model. The molecular structure and partial charges are calculated using two semi empirical (PM3, AM1) and one empirical (Gasteiger and Marcili) methods. A distance-dependent dielectric function is used to represent the interfacial dielectric constant at the aqueous subphase. The present study indicates that a pair of molecules have favorable interaction at specific ranges of mutual orientations. Other orientations are favorable but at larger separations. Favorable electrostatic interaction at a specific combination of orientation and short separations of the head groups significantly contribute to the total energy. The curvature of the domain boundary is suggested to be driven by the favorable arrangement which is dependent on the pair potential of molecules. The use of charges obtained by the PM3 and GM do not lead to a significant variation of the orientation-dependent features, while the AM1 predicts higher partial charges and interactions are stronger than the former two methods. However, orientation-dependent features remain the same. The variations in the LJ parameters and charges indicate that the conclusions made are insensitive to the choice of parameters. The mutual favorable interaction predicted by calculation agree with the handedness of curvature of domains.     

ab Initio Calculation on Tautomerism of Pyrazolin—5—ones

ａｂＩｎｉｔｉｏＣａｌｃｕｌａｔｉｏｎｏｎＴａｕｔｏｍｅｒｉｓｍｏｆＰｙｒａｚｏｌｉｎ－５－ｏｎｅｓ￥ＳｕｎＲｅｎ－Ａｎ；ＷａｎｇＣｈａｎｇ－Ｓｈｅｎｇ；ＹａｎｇＺｈｏｎｇ－Ｚｈｉ（ＤｅｐａｒｔｍｅｎｔｏｆＣｈｅｍｉｓｔｒｙ，ＬｉａｏｎｉｎｇＮｏｒｍ...     

AN INVESTIGATION ON THE ab initio CALCULATION OF INTERMOLECULAR INTERACTION——NON-BONDED INTERACTION IN THE HYDROGEN BOND COMPLEX (HF)_2

An ab initio calculation of the hydrogen bond complex (HF)_2 is given with the 6-311 G~(**) basis set, according to which the potential surface around the balance point of the distancès and the orientations between two HF molecules is obtained. The atomic charges in the system are calculated with the PD/LSF method (potential-derived/least-square-fitting method) and then an analysis of the hydrogen bond interaction between two HF molecules is given with the (exp-6-1) potential function, by means of which it is shown that the main interaction between them is not an electro-static but a charge transfer one. The potential curve between two HF molecules is like a Morse function.     

Crystal,Molecular, and Electronic Structure of 1-Acetyl-indoline and Derivatives

The crystal and molecular structures of the following molecules have been determined: 1-acetyl-indoline, 1-acetyl-5-nitro-indoline, l-acetyl-5-nitro-7-bromo-indoline, 1-acetyl-5-bromo-7-nitroindoline, and l-acetyl-5-bromo-7-nitro-indol. Molecular orbital calculations are performed for these compounds and two related species.     

Theoretical study on interaction of different coordination modes of BH4 ligand with transition metal in [TM(BH4)(CO)4] (TM = Cr, Mo)

Density functional calculations were performed on bonding and structural features of [(ηⁿ-BH₄)TM(CO)₄]⁻ (n = 1, 2, 3; TM = Cr, Mo) complexes. Calculations show that the ground state is bidentate which is in good agreement with experimental results. It has been found that the bridge and terminal hydrogen atoms will interchange by two pathways: (i) twist of BH₄ about one of the bridge B-H and (ii) twist of BH₄ about one of the terminal B-H. The molecular orbital calculations and natural bond orbital methodologies for different isomers of these complexes have been evaluated. The final results indicate that case (i) is more preferable relative to another case.     

Study on the effects of fluorine and oxygen deficiency on YBa2Cu307 byab initio method

The calculations of clusters modeling the fluorine-doping and oxygen deficiency of YBa₂Cu₃0₇ have been performed by the method of all-electronab initio Hartree-Fock with self-consistent crystal field. Results show that in CuO planes electric charge significantly increases, the chemical valence of Cu decreases and the covalent bonding of Cu-O greatly weakens owing to oxygen deficiency, while the effect of F restores the local electronic structure of YBa₂Cu₃O₇. The reported opinion that F occupied the oxygen vacancy in Cu-O chains seems disputable according to the calculated bonding characteristics Project supported by the National Natural Science Foundation of China.     

由NO分子从头计算结果和紫外光电子能谱确定NO的分子结构

对于NO分子轨道的能级次序现有两种不同的说法,不少人根据NO与O_2~+是等电子体,由O_2~+分子的能级次序确定NO分子组态。本文用NO和O_2~(+)从头计算结果以及NO的紫外光电子能谱相结合的方法说明NO的5σ轨道是弱成键轨道,5σ的轨道能稍高于1π轨道能,也就是说NO的能级次序是与N_2分子相同的;NO与O_2~+虽是等电子体,但是能级次序并不相同,因此由O_2~+的能级次序确定NO分子的电子组态是不妥的。