PubMed Health⌕ Search

Biomedical subjects

Anatoliy Volkov

Publications and source records attributed to Anatoliy Volkov.

5 recordsLinked to original sources
PubMed2006

On the calculation of the electrostatic potential, electric field and electric field gradient from the aspherical pseudoatom model.

Accurate, yet simple and efficient, formulae are presented for calculation of the electrostatic potential (ESP), electric field (EF) and electric field gradient (EFG) from the aspherical Hansen-Coppens pseudoatom model of electron density [Hansen & Coppens (1978). Acta Cryst. A34, 909-921]. They are based on the expansion of |r' - r|(-1) in spherical harmonics and the incomplete gamma function for a Slater-type function of the form R(l)(r) = r(n) exp(-alpha ). The formulae are valid for 0 < or = r < or = infinity and are easily extended to higher values of l. Special treatment of integrals is needed only for functions with n = l and n = l + 1 at r = 0. The method is tested using theoretical pseudoatom parameters of the formamide molecule obtained via reciprocal-space fitting of PBE/6-31G** densities and experimental X-ray data of Fe(CO)(5). The ESP, EF and EFG values at the nuclear positions in formamide are in very good agreement with those directly evaluated from density-functional PBE calculations with 6-31G**, aug-cc-pVDZ and aug-cc-pVTZ basis sets. The small observed discrepancies are attributed to the different behavior of Gaussian- and Slater-type functions near the nuclei and to imperfections of the reciprocal-space fit. An EF map is displayed which allows useful visualization of the lattice EF effects in the crystal structure of formamide. Analysis of experimental 100 K X-ray data of Fe(CO)(5) yields the value of the nuclear quadrupole moment Q((57)Fe(m)) = 0.12 x 10(-28) m(2) after taking into account Sternheimer shielding/antishielding effects of the core. This value is in excellent agreement with that reported by Su & Coppens [Acta Cryst. (1996), A52, 748-756] but slightly smaller than the generally accepted value of 0.16 +/- 5% x 10(-28) m(2) obtained from combined theoretical/spectroscopic studies [Dufek, Blaha & Schwarz (1995). Phys. Rev. Lett. 25, 3545-3548].

Journal Article↗
PubMed2004

The interplay between experiment and theory in charge-density analysis.

The comparison of theory and experiment remains a cornerstone of scientific inquiry. Various levels of such comparison applicable to charge-density analysis are discussed, including static and dynamic electron densities, topological properties, d-orbital occupancies and electrostatic moments. The advantages and drawbacks of the pseudoatom multipole are discussed, as are the experimentally constrained wavefunctions introduced by Jayatilaka and co-workers, which combine energy minimization with the requirement to provide a reasonable fit to the X-ray structure factors. The transferability of atomic densities can be exploited through construction of a pseudoatom databank, which may be based on analysis of ab initio molecular electron densities, and can be used to evaluate a host of physical properties. Partitioning of theoretical energies with the Morokuma-Ziegler energy decomposition scheme allows direct comparison with electrostatic interaction energies obtained from electron densities represented by the pseudoatom formalism. Compared with the Buckingham expression for the interaction between non-overlapping densities, the agreement with theory is much improved when a newly developed hybrid EP/MM (exact potential/multipole model) method is employed.

Models, Chemical↗
PubMed2004

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning.

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.

Algorithms↗
PubMed2002

Aspherical-atom scattering factors from molecular wave functions. 1. Transferability and conformation dependence of atomic electron densities of peptides within the multipole formalism.

In this study, the feasibility of building a database of theoretical atomic deformation density parameters applicable to the construction of the densities of biomacromolecules and to the interpretation of their X-ray diffraction data is discussed. The procedure described involves generation of valence-only structure factors of tripeptides calculated from theoretical densities at the B3LYP level and the refinement of multipole parameters against these simulated data. Our results so far indicate that the backbone pseudoatoms extracted in such a way are highly transferable and fairly invariant with respect to rotations around single bonds in the peptide framework. The ultimate goal is to use the aspherical-atom database for improved macromolecular refinements that are based on high-resolution data and for prediction of electrostatic properties of larger molecules.

Algorithms↗