Search

A new population analysis, ADCHα-I, based on the interpolation between the Hirshfeld (H) and the iterative Hirshfeld (H-I) methods through a parameter α and on the atomic dipole moment corrected Hirshfeld (ADCH) methodology is proposed, in combination with the constrained dipole moment density functional theory (CD-DFT) previously developed, to determine the charge distributions of force fields. Following this approach, the electronic density of the isolated molecule is determined for the value of the dipole moment that reproduces the experimental dielectric constant, in order to incorporate through this property the effects of the surrounding molecules in the liquid, and to carry on this information to the molecular simulation, the new population analysis is built to obtain the set of charges that reproduces this dipole moment. By selecting α = 1/2, one is led to charges that are larger than the ones obtained through H and ADCH and smaller than those of H-I and that incorporate, at the local level, information about the response of isolated atoms to donate or to accept charge, which is not considered in ADCH. The results obtained for several liquid properties indicate that the combination of CD-DFT with this population analysis leads to a good description of the charge distributions in force fields used in molecular simulations. [ABSTRACT FROM AUTHOR]

The use of molecular simulation methods to reproduce and predict experimental information has become an excellent option for designing new experiments. The success of these methods depends on the force field used. Their quality depends on the interaction parameters and the target properties used to obtain them. In this work, we review and extend the methodology developed by our group to obtain interaction parameters of polar liquids. In addition, an automated workflow is developed to obtain new force field parameters and the analysis of physical properties. It is applied to small molecules of various functional groups in organic chemistry. The dielectric constant, surface tension, liquid density, self-diffusion coefficient and solubility of a solute in water are used as target properties. The simulations are developed in liquid phase, liquid–vapour and liquid–liquid equilibriums. The procedure can be extended to include other target properties and has the advantage that the parameters can be obtained in shorter times. It also opens the possibility of studying a wide range of problems of scientific and technological interest. [ABSTRACT FROM AUTHOR]

The implementation of the cyclic cluster model (CCM) for molecular mechanics is presented in the framework of the computational chemistry program deMon2k. Because the CCM is particularly well-suited for the description of periodic systems with defects, it can be used for periodic QM/MM approaches where the non-periodic QM part is treated as a defect in a periodic MM surrounding. To this end, we present here the explicit formulae for the evaluation of the Ewald sum and its first- and second-order derivatives as implemented in deMon2k. The outlined implementation was tested in molecular dynamics (MD) simulations and periodic structure optimization calculations. MD simulations of an argon system were carried out using the Nosé-Hoover chain (NHC) thermostat and the Martyna-Tobias-Klein (MTK) barostat to control the temperature and pressure of the system, respectively. For the validation of CCM structure optimization a set of molecular crystals were optimized using the Ewald method for the evaluation of the electrostatic interactions. Two optimization procedures for the determination of the atomic positions and CCM cell parameters were tested. Our results show that the simultaneous optimization of the atomic positions and cell parameters is most efficient. [ABSTRACT FROM AUTHOR]