Extended gamma sources modelling using multipole expansion: Application to the Tunisian gamma source load planning

In this work, we focus on the effect of the source geometry for radiation profile studies using the generalized Laplace's expansion. We develop a new method to express the bare gamma photon flux rate in terms of the standard Cartesian multipole moments. Using the properties of these moments we establish a close relationship between the radiation profile and the geometrical features of the source. As application a detailed study of the arrangement of the unit pencil sources of the Tunisian irradiation facility is performed. Using this method, millions of possible configurations for various load plans were investigated, in a reasonable time and even multisteps scenarios were considered. Furthermore, using this analytical method it was even possible to optimize the activity of each new unit source.     

Torque and atomic forces for Cartesian tensor atomic multipoles with an application to crystal unit cell optimization

New equations for torque and atomic force are derived for use in flexible molecule force fields with atomic multipoles. The expressions are based on Cartesian tensors with arbitrary multipole rank. The standard method for rotating Cartesian tensor multipoles and calculating torque is to first represent the tensor with n indexes and 3ⁿ redundant components. In this work, new expressions for directly rotating the unique (n + 1)(n + 2)/2 Cartesian tensor multipole components Θ_pqr are given by introducing Cartesian tensor rotation matrix elements X( R ). A polynomial expression and a recursion relation for X( R ) are derived. For comparison, the analogous rotation matrix for spherical tensor multipoles are the Wigner functions D( R ). The expressions for X( R ) are used to derive simple equations for torque and atomic force. The torque and atomic force equations are applied to the geometry optimization of small molecule crystal unit cells. In addition, a discussion of computational efficiency as a function of increasing multipole rank is given for Cartesian tensors. © 2016 Wiley Periodicals, Inc.     

Cumulative potential-derived atomic multipoles

A method to derive the atomic multipole moments cumulatively up to quadrupole moments was developed. The multipole moments are obtained by least-square simulating the molecular electrostatic potentials. Only the components of the term of highest order in the atomic multipole expansion are optimized while the lower terms remain fixed. The calculations on HF, H₂O and NH₃ show that the cumulative method can give reasonable qualitative and fairly good quantitative results.     

The effects of atomic multipole moments obtained by the potential-derived method on hydrogen bonding

A potential-derived atomic multipole method called the cumulative potential-derived atomic multipole method is developed, with which electrostatic atomic multipole moments are derived by fitting the molecular electric potential in a cumulative way. It is applied to the hydrides of N , O , F , S , Cl , and methanol and the hydrogen-bonded dimers formed between them. The relationship between atomic multipole moments and molecular charge distributions is found. The structures calculated with Buckingham's electrostatic model are in good agreement with experiments. The phenomena of nonlinear structures of most H -bonded complexes—the deviations of symmetry axes of electron donors from H bonds—and correct distinguishing between two alternative structures are attributed to atomic dipole and quadrupole moments. Compared with other methods, this method has a quantitative and qualitative advantage and simple algorithm. The main conclusion is that the atomic multipole moments play a substantial role, although a potential-derived charge model was deemed sufficient previously. © 1993 John Wiley & Sons, Inc.     

Accuracy of distributed multipoles and polarizabilities: comparison between the LoProp and MpProp models

Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.     

Accuracy and tractability of a kriging model of intramolecular polarizable multipolar electrostatics and its application to histidine

We propose a generic method to model polarization in the context of high‐rank multipolar electrostatics. This method involves the machine learning technique kriging, here used to capture the response of an atomic multipole moment of a given atom to a change in the positions of the atoms surrounding this atom. The atoms are malleable boxes with sharp boundaries, they do not overlap and exhaust space. The method is applied to histidine where it is able to predict atomic multipole moments (up to hexadecapole) for unseen configurations, after training on 600 geometries distorted using normal modes of each of its 24 local energy minima at B3LYP/apc‐1 level. The quality of the predictions is assessed by calculating the Coulomb energy between an atom for which the moments have been predicted and the surrounding atoms (having exact moments). Only interactions between atoms separated by three or more bonds (“1, 4 and higher” interactions) are included in this energy error. This energy is compared with that of a central atom with exact multipole moments interacting with the same environment. The resulting energy discrepancies are summed for 328 atom–atom interactions, for each of the 29 atoms of histidine being a central atom in turn. For 80% of the 539 test configurations (outside the training set), this summed energy deviates by less than 1 kcal mol^?1. © 2013 Wiley Periodicals, Inc.     

Polarisable multipolar electrostatics from the machine learning method Kriging: an application to alanine

We present a polarisable multipolar interatomic electrostatic potential energy function for force fields and describe its application to the pilot molecule MeNH-Ala-COMe (AlaD). The total electrostatic energy associated with 1, 4 and higher interactions is partitioned into atomic contributions by application of quantum chemical topology (QCT). The exact atom–atom interaction is expressed in terms of atomic multipole moments. The machine learning method Kriging is used to model the dependence of these multipole moments on the conformation of the entire molecule. The resulting models are able to predict the QCT-partitioned multipole moments for arbitrary chemically relevant molecular geometries. The interaction energies between atoms are predicted for these geometries and compared to their true values. The computational expense of the procedure is compared to that of the point charge formalism.     

A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank

A finite field method for calculating spherical tensor molecular polarizability tensors α_lm;l′m′ = ?Δ_lm/??_l′m′* by numerical derivatives of induced molecular multipole Δ_lm with respect to gradients of electrostatic potential ?_l′m′* is described for arbitrary multipole ranks l and l′. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole–hexadecapole level are calculated for water using the following ab initio methods: Hartree–Fock (HF), Becke three‐parameter Lee‐Yang‐Parr exchange‐correlation functional (B3LYP), Møller–Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Ab Initio Study on the Anti-HIV Activities of Hydroxyflavones

1 INTRODUCTION AIDS, in other words, Acquired Immunodefici-encysyndrome, is a viral contagious disease withhigh death rate. However, effective treatment againstit has not appeared so far. At present, about 100kinds of natural compounds that could be extractedfrom natural products or artificially synthesized havebeen found to have anti-HIV activities to some ex-tent. They belong to flavonoid, terpenoid, steroid, cou-marin, peptide, alkaloid etc., respectively. Previousstudies show that…     

Basis set convergence studies of Hartree-Fock calculations of molecular properties within the resolution of the identity approximation

Within the resolution of the identity (RI) method, the convergence of the Hartree-Fock (HF) total molecular energy and the multipole moments in the course of the combined regular expansion of the molecular and auxiliary (RI) basis sets is studied. Dunning's cc-pVXZ series is used for both the molecular and the RI basis sets. The results show the calculated quantities converge to the HF limit when both the molecular and the RI basis sets are expanded from correlation-consistent polarized valence double zeta to correlation-consistent polarized valence sextuple zeta. Combinations of molecular/RI basis sets sufficient for convergence of the total energy and of the multipole moments at various accuracy levels have been determined. A measure of the RI basis set incompleteness is suggested and discussed. As it is significantly faster than the standard HF algorithm for small and midsize molecules, the RI-HF method, together with appropriate expanding series of both molecular and RI basis sets, provide an efficient tool to estimate and control the error of the Hartree-Fock calculations due to the finite basis set.     

Prediction of conformationally dependent atomic multipole moments in carbohydrates

The conformational flexibility of carbohydrates is challenging within the field of computational chemistry. This flexibility causes the electron density to change, which leads to fluctuating atomic multipole moments. Quantum Chemical Topology (QCT) allows for the partitioning of an “atom in a molecule,” thus localizing electron density to finite atomic domains, which permits the unambiguous evaluation of atomic multipole moments. By selecting an ensemble of physically realistic conformers of a chemical system, one evaluates the various multipole moments at defined points in configuration space. The subsequent implementation of the machine learning method kriging delivers the evaluation of an analytical function, which smoothly interpolates between these points. This allows for the prediction of atomic multipole moments at new points in conformational space, not trained for but within prediction range. In this work, we demonstrate that the carbohydrates erythrose and threose are amenable to the above methodology. We investigate how kriging models respond when the training ensemble incorporating multiple energy minima and their environment in conformational space. Additionally, we evaluate the gains in predictive capacity of our models as the size of the training ensemble increases. We believe this approach to be entirely novel within the field of carbohydrates. For a modest training set size of 600, more than 90% of the external test configurations have an error in the total (predicted) electrostatic energy (relative to ab initio) of maximum 1 kJ mol^?1 for open chains and just over 90% an error of maximum 4 kJ mol^?1 for rings. © 2015 Wiley Periodicals, Inc.     

Calculation of molecular quadrupole moments and a demonstration of the importance of overlap densities in the theory of polyatomic molecules

A computational method for calculating quadrupole moments from molecular wave functions in a Slater orbital basis set is described. Using both IEHT and CNDO wave functions quadrupole moments for a series of polyatomic molecules are calculated. They are compared with experimental results and the IEHT wave functions are found to give agreement with experiment while CNDO wave functions do not. The importance of bicentric densities (overlap densities) in the calculation of multipole moments is shown. This is followed by a discussion of the usefulness of these wave functions for a quantitative characterization of the electronic structure of large molecules.     

SCF Deformation densities and electrostatic potentials of purines and pyrimidines

4-31G wave functions have been computed for five purines and pyrimidines. The calculated deformation densities have been partitioned into atomic fragments, which were integrated to yield atomic multipole moments. The transferability of atomic fragments between related molecules was verified by constructing model maps for uracil and guanine from appropriate fragments of cytosine and adenine. Model electrostatic potentials calculated from the moments of model atoms are similar to the corresponding 4-31G potentials. Comparison of 4-31G and 4-31G** deformation densities of cytosine provides simple rules for estimating the effects of polarization functions on the atomic multipole moments of most atom types occurring in the purines and pyrimidines. These rules were applied to the other molecules and yielded reasonable approximations for their molecular dipole moments. Substituting CH₃ for H has little effect on the deformation density beyond the substitution center.     

The electrostatic potential generated by topological atoms. II. Inverse multipole moments

Quantum chemical topology defines finite atoms, whose bounded electron density generates a well-defined electrostatic potential. A multipole expansion based on spherical tensors provides a potential that is formally convergent outside the divergence sphere. Part I of this series [P. L. A. Popelier and M. Rafat, Chem. Phys. Lett.376, 148 (2003)] showed that a continuous multipole expansion expands the convergence region, thereby allowing the electrostatic potential to be evaluated at short range. Here, we propose a different method, based on "inverse" multipole moments, enabling an expansion that converges everywhere. These moments are defined by inverse (i.e., negative) powers of the magnitude of the position vector describing the electron density inside the atom. We illustrate this technique on nitrogen in N(2), oxygen in H(2)O, and oxygen in the phenolic group of the amino acid tyrosine. The proposed method constitutes a considerable advance over the method presented in Part I.     

Efficient implementation of the fast multipole method

A number of computational techniques are described that reduce the effort related to the continuous fast multipole method, used for the evaluation of Coulomb matrix elements as needed in Hartree-Fock and density functional theories. A new extent definition for Gaussian charge distributions is proposed, as well as a new way of dividing distributions into branches. Also, a new approach for estimating the error caused by truncation of multipole expansions is presented. It is found that the use of dynamically truncated multipole expansions gives a speedup of a factor of 10 in the work required for multipole interactions, compared to the case when all interactions are computed using a fixed multipole expansion order. Results of benchmark calculations on three-dimensional systems are reported, demonstrating the usefulness of our present implementation of the fast multipole method.     

Unified Treatment for Arbitrary-rank Cartesian Electric and Magnetic Multipole Moment Operators

In this study, the theory of cartesian electric and magnetic multipole moments is extended in a unified way. The general analytical expressions for distinct components of arbitrary rank cartesian electric and magnetic multipole moment operators are derived as linear combination of corresponding spherical operators, which can be used as interconversion between cartesian and spherical electric and magnetic multipole moment tensors. The transformation properties, such as translation and rotation of cartesian electric and magnetic multipole moments are given in a very simple general form. The relationship between distinct and linearly independent components of cartesian multipole moment tensors in system of linear symmetry is also presented. The formulae obtained in this paper can be utilized to calculate the interaction energies between charge distributions.     

Atomic multipoles: Electrostatic potential fit,local reference axis systems,and conformational dependence

Currently, all standard force fields for biomolecular simulations use point charges to model intermolecular electrostatic interactions. This is a fast and simple approach but has deficiencies when the electrostatic potential (ESP) is compared to that from ab initio methods. Here, we show how atomic multipoles can be rigorously implemented into common biomolecular force fields. For this, a comprehensive set of local reference axis systems is introduced, which represents a universal solution for treating atom‐centered multipoles for all small organic molecules and proteins. Furthermore, we introduce a new method for fitting atomic multipole moments to the quantum mechanically derived ESP. This methods yields a 50–90% error reduction compared to both point charges fit to the ESP and multipoles directly calculated from the ab initio electron density. It is shown that it is necessary to directly fit the multipole moments of conformational ensembles to the ESP. Ignoring the conformational dependence or averaging over parameters from different conformations dramatically deteriorates the results obtained with atomic multipole moments, rendering multipoles worse than partial charges. © 2012 Wiley Periodicals, Inc.