Induction correction model for rotation of two or three dihedral angles

In a previous work we have introduced an intramolecular induction correction model. In this work we have used the model to calculate the total dipol moment of six molecules as a function of two or three dihedral angles that are simultaneously varied in the molecules. It is found that the induction model behaves very well for the systems studied when compared with a regular force field model where fixed charges and dipoles are rotated along with the atoms of the molecules. This suggests that the proposed induction correction model can be used to model systems containing several dihedral angles around which rotations are performed.     

An intramolecular induction correction model of the molecular dipole moment

A model for intramolecular polarization is presented. It is used to describe the changes in the molecular charge distribution occurring as a response to changes of dihedral angles in the molecule. The model is based on a multicenter multipole distribution of the molecular charge distribution. The electric field from this charge distribution induce dipole moments in the same molecule. The model contains atom type parameters to describe the damping of the electric field. A total of four atom types are used. The parameters are fitted to a calibration set with various functional groups, and tested against a validation set. The error obtained for the calibration set is reduced by 92% and by 88% for the validation set, if compared to an accurate state-of-the-art force field. It is shown that rotating the non-polarizable multicenter multipole distribution for the equilibrium geometry gives too large dipole moments for dihedral angles deviating from the equilibrium geometry. This will lead to too large long-range attractions in simulations. This problem is overcome by using the dipole polarizability correction suggested here, which gives dipole moments very close to the Hartree-Fock dipole moments obtained from reference calculations.     

The coupling constant polarizability and hyperpolarizabilty of 1J(NH) in N-methylacetamide, and its application for the multipole spin-spin coupling constant polarizability/reaction field approach to solvation

We present a benchmark study of a combined multipole spin-spin coupling constant (SSCC) polarizability/reaction field (MJP/RF) approach to the calculation of both specific and bulk solvation effects on SSCCs of solvated molecules. The MJP/RF scheme is defined by an expansion of the SSCCs of the solvated molecule in terms of coupling constant dipole and quadrupole polarizabilities and hyperpolarizabilities derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics (MD) simulations thereby accounting for solute-solvent dynamical effects. The MJP/RF method is benchmarked against polarizable QM/MM calculations for the one-bond N-H coupling constant in N-methylacetamide. The best agreement between the MJP/RF and QM/MM approaches is found by truncating the electric field expansion in the MJP/RF approach at the linear electric field level. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the coupling constant (hyper-)polarizabilities and find that they are affected by the basis set in a way similar to the coupling constants themselves.     

分子极化效应与烷烃、醇的空腔表面积

烷烃以及醇在水中形成的空腔表面积CSA可由下式计算:CSA=122.984+8.65736SBL-103.862ΔMPEI。其中SBL是分子中所有化学键长度的总和, ΔMPEI是含有相同碳原子数目的(烷烃或者醇)支链异构体分子极化效应指数MPEI支与直链异构体分子极化效应指数MPEI直的差值。     

Energy decomposition in molecular complexes: implications for the treatment of polarization in molecular simulations

This study examines the contribution of electrostatic and polarization to the interaction energy in a variety of molecular complexes. The results obtained from the Kitaura-Morokuma (KM) energy decomposition analysis at the HF/6-31G(d) level indicate that, for intermolecular distances around the equilibrium geometries, the polarization energy can be determined as the addition of the polarization energies of interacting blocks, as the mixed polarization term is typically negligible. Comparison of KM and QM/MM results shows that the electrostatic energy determined in the KM method is underestimated (in absolute value) by QM/MM methods. The reason of such underestimation can be attributed to the simplified representation of treating the interaction between overlapping charge distribution by the interaction of a QM molecule with a set of point charges. Nevertheless, the polarization energies calculated by KM and QM/MM methods are in close agreement. Finally, a consistent, automated strategy to derive charge distributions that include implicitly polarization effects in pairwise, additive force fields is presented. The strategy relies in the simultaneous fitting of electrostatic and polarization energies computed by placing a suitable perturbing particle at selected points around the molecule. The suitability of these charges to describe molecular interactions is discussed.     

The effect of confinement on the electronic energy and polarizability of a hydrogen molecular ion

The electronic energy and the polarizability of a confined hydrogen molecular ion in the ground state and the first excited state, for cavities of different volumes, are calculated using the variational method. In the treatment adopted an alternative molecular wave function is introduced with only one variational parameter and based on wave functions used for confined atoms. © 2016 Wiley Periodicals, Inc.     

Electronic polarization in dipolar organic molecular semiconductors: The case study of 1,2,3,4-tetrafluoro-6,7-dimethylnaphthalene crystal

Electronic polarization has an important impact on the site energies of charge carriers that play a key role in determining the charge transport in organic semiconductors.Dipolar molecules have strong intermolecular interactions and widespread applications in organic optoelectronics.Howeve r,compared with nonpolar organic semiconductors,electronic polarization for dipolar systems has been rarely studied.Here,taking 1,2,3,4-tetrafluoro-6,7-dimethylnaphthalene as representative,we have calculated the electronic polarization energies of dipolar organic molecular crystals by means of a polarizable forcefield method.Surprisingly,our results point to that the polarization energies for this dipolar system are similar to those of nonpolar systems.In addition,theπ-πstack contributes only about 30%~40%to the total polarization energy,thus the polarization effects along the three dimensions should be treated equally even for the one-dimensional stacking crystals.     

The calculation of electric dipole vibronic intensities in centrosymmetric coordination compounds using the crystal field-closure-ligand polarization model

Equations are derived for the evaluation of the contribution to the electric dipole vibronic intensity of d-d and f-f transitions in octahedral complex ions. Expressions for the evaluation of the crystal-field and ligand polarization terms and the cross term between are presented in symmetry adapted form and tables of the required coupling constants are given.     

Application of the FAKE molecular-orbital method to diatomic molecules XY (X,Y = H,F, Cl,Br, I)

Theoretical Chemistry Accounts - The FAKE (fast, accurate kinetic energy) method of semiempirical molecular orbital calculation is applied to diatomic molecules XY (X, Y= H, F, Cl, Br, I). The...     

The effect of electric field on the electronic structure of DNA. I. Calculation of the polarizability and of the permanent dipole moment for the nucleotide bases and base pairs

Generalizing the method of Hirschfelder, Curtiss and Bird for the calculation of the polarizability of molecules, we have calculated both the diagonal and off-diagonal polarizability tensor elements of the nucleotide bases and base pairs. In the course of calculations also the components of the permanent dipole moments have been computed for these systems. On the basis of the results obtained we have determined also the critical field strengths necessary to break off a base pair and thus to induce the unwinding of the DNA double helix. According to the data obtained, a field strength of the order 10⁸ V/cm can induce the unwinding of DNA if its direction lies in the plane of the base pairs.     

磁场效应对理解水分解反应物的可能意义

电催化水分解由两个基元反应构成,即析氢反应(HER)和析氧反应(OER).开发强大的HER和OER技术需要在分子层面理解反应机理,然而,目前水分解反应的反应物还没有完全确定.本文利用磁场来研究HER中的质子传输和OER中的氢离子传输,对确定HER和OER中真实的反应物具有重要意义.磁场是改变离子等带电物质运动的一种有效...     

The effect of electric field geometry on the performance of electromembrane extraction systems: Footprints of a third driving force along with migration and diffusion

The distribution of electric field vectors was first calculated for electromembrane extraction (EME) systems in classical and cylindrical electrode geometries. The results showed that supported liquid membrane (SLM) has a general field amplifying effect due to its lower dielectric constant in comparison with aqueous donor/acceptor solutions. The calculated norms of the electric field vector showed that a DC voltage of 50 V can create huge electric field strengths up to 64 kV m⁻¹ and 111 kV m⁻¹ in classical and cylindrical geometries respectively. In both cases, the electric field strength reached its peak value on the inner wall of the SLM. In the case of classical geometry, the field strength was a function of the polar position of the SLM whereas the field strength in cylindrical geometry was angularly uniform. In order to investigate the effect of the electrode geometry on the performance of real EME systems, the analysis was carried out in three different geometries including classical, helical and cylindrical arrangements using naproxen and sodium diclofenac as the model analytes. Despite higher field strength and extended cross sectional area, the helical and cylindrical geometries gave lower recoveries with respect to the classical EME. The observed decline of the signal was proved to be against the relations governing migration and diffusion processes, which means that a third driving force is involved in EME. The third driving force is the interaction between the radially inhomogeneous electric field and the analyte in its neutral form.     

Laser‐induced periodic surface structure on a polymer surface: The effect of a prerubbing treatment on the surface structure

The influence of the mechanical rubbing of a polyimide (PI) film on the laser‐induced periodic structure (LIPS) was demonstrated. The periodicity and amplitude of LIPS were greater when the rubbing direction was parallel to the laser polarization direction. The amplitude became small and the periodicity of LIPS did not show an obvious change when the rubbing direction was perpendicular to the laser polarization direction. The effect of the rubbing pretreatment on LIPS was explained on the basis of the wave‐guide effect of rubbing‐induced microgrooves on LIPS formation. The orientation of PI chains induced by mechanical rubbing was relaxed after laser irradiation, and a new orientation of PI chains was formed during the LIPS formation. When the rubbing direction was perpendicular to the laser polarization direction, the orientation of PI chains remained in the rubbing direction. The laser‐irradiated, perpendicularly rubbed PI surface could be used to verify the effects of surface morphologies and intermolecular interactions on liquid‐crystal alignment. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 1273–1280, 2003     

Mechanistic aspects of the Wolff-Kishner reaction. V. The cation effect on the reaction of benzophenone hydrazone in butyl carbitol and decyl alcohol

The rate of the Wolff-Kishner reaction of benzophenone hydrazone in butyl carbitol increases as the cation of the alkoxide base is varied in the order K>Na>Li>Mg. The replacement of butyl carbitol by 1-decanol also accelerates the reaction, and an additional increase is caused by the presence of crown ether. On the basis of changes in the activation parameters, it is concluded that the reactivity of the hydrazone anion increases as the contact ion pair becomes more fully solvated and dissociated.     

外偶极电场作用下H_2O,D_2O和T_2O的可逆分解电压

采用密度泛函B3P86方法和6-311++G(3df,3pf)基组,计算了在-0.05～0.05a.u.外偶极电场作用下,H2O,D2O,T2O,H2,D2,T2,O2的电子能量、核运动能量和熵值,在此基础上通过计算H2O(g)→H2(g)+O2(g)、D2O(g)→D2(g)+O2(g)、T2O(g)→T2(g)+O2(g)的焓变ΔH、熵变ΔS、Gibbs函数变化ΔG,最后得到了H2O,D2O,T2O的可逆分解电压Er.计算结果表明,外偶极电场存在时,H2O,D2O,T2O的Gibbs自由能变ΔG和可逆分解电压Er都有明显的变化,当外偶极电场正方向增加时,其Gibbs自由能变ΔG和可逆分解电压Er均趋于线性增加;当外偶极电场负方向增加时,其Gibbs自由能变ΔG和可逆分解电压Er均趋于线性减小;在相同外偶极电场作用下,Gibbs自由能变ΔG和可逆分解电压Er随H2O,D2O,T2O依次增加.     

The temperature window of minimum flow resistance in melt flow of polyethylene. Further studies on the effect of strain rate and branching

The existence of a narrow temperature window (150–153°C) of smooth extrudability coupled with a minimum in flow resistance (extrusion pressure) in high-molecular weight polyethylene (>4 × 10⁵ g mol^?1) was the subject of a previous article where it was associated with strain-induced formation of the mobile hexagonal mesophase. The new findings of this note show that this minimum in flow resistance only sets in above a critical strain rate; this is interpreted in terms of the requirement of a critical strain rate in order to stretch molecules to their fully extended configuration. Furthermore, this critical strain rate is shown to be higher for lower molecular weight materials, in agreement with a priori considerations. Additionally, the temperature at which the pressure minimum occurs in a polyethylene containing methyl branches shifts to a significantly lower value than that for the linear material. This is interpreted in terms of the ? CH₃ groups raising the crystal free energy, thereby lowering the temperature at which the transition to the hexagonal phase occurs.     

Comment on “The effect of confinement on the electronic energy and polarizability of a hydrogen molecular ion” by Josimar Fernando da Silva,Fabricio Ramos Silva and Elso Drigo Filho,Int. J. Quantum Chem. 2016, 116, 497–503

Various conceptual and numerical errors reported in the paper cited in the title of this contribution are pointed out and clarified. Our main concern is the evident violation of the variational principle, which is used in that paper for the study of the confined hydrogen molecular ion to obtain both energies and properties.