Just for the record, except for charge density studies most small molecule structures are refined with neutral atom scattering factors even when ions such as Cl- are present. For example SHELX uses "4 Gaussian plus const: International Tables for Crystallography (1992)". Users rarely input ionic scattering factors rather than using the default neutral atom scattering factors, possibly they are afraid of a Category A Alert from CheckCIF if they have a charged crystal. In practice, the difference between ionic and neutral atom scattering factors is mostly absorbed by the displacement parameters (B-values) and the R-factors are very similar. The main reason why routine small molecule R-values are rarely much less than 2% is the assumption that the scattering factors are spherically symmetrical, i.e. bonding and lone-pair electrons are ignored. This could be addressed by the inclusion of invarioms (precalculated aspherical scattering factors that depend on the local chemical environment of an atom) in Refmac and Phenix Refine.
George On Wed, Nov 02, 2011 at 08:39:51PM -0700, Pavel Afonine wrote: > Continuing on the subject, as far as I know there are at least three flavors > of > form-factors currently used in refinement programs: > > "4 gaussian plus const": > International Tables for Crystallography (1992) > > "5 gaussian plus const": > D. Waasmaier & A. Kirfel. Acta Cryst. (1995). A51, 416-431. "New analytical > scattering-factor functions for free atoms and ions" > > "n-gaussian" (n determined dynamically) > Grosse-Kunstleve RW, Sauter NK & Adams PD. Newsletter of the IUCr Commission > on > Crystallographic Computing 2004, 3:22-31. "cctbx news" > > All three are available in PHENIX (the 3rd is used by default), and I presume > the first one is used in CNS and Refmac, if I remember correctly (the authors > of respective programs please correct me). > > Pavel > > On Wed, Nov 2, 2011 at 7:17 PM, Ivan Shabalin <[email protected]> wrote: > > Hi James! > > Thank you very much for the gnuplot-ish version of ${CLIBD}/atomsf.lib!! > It > works very nice and is very useful for education! > > As I understand, the form factor is the Fourier transform of electron > charge density. It is plotted as f(electrons) vs sin(tetta)/lambda and is > approximated as 5 Gaussian (Cromer and Mann) in REFMAC. And you made > reverse Fourier transform of the approximation and plotted the electron > density distribution in the real space. > > So, can I ask, what unit is x? Is it angstrom? > And what is Y? is it e/A3 (electron density)? > > I found, that at Bf=20, density profiles look almost the same for ions and > atoms (Mg2+ and Mg, Cl- and Cl). Does that means, there is no sense to > specify atomic charge in refmac refinement? It looks a bit strange, > because > the numbers of electrons are different. Or decreasing in number of > electrons is compensated with significant decrease in atom size (that can > have the same effect as Bf lowering)? With Bf=0 the difference in curves > is > significant. > > With best regards, > Ivan Shabalin > > -- Prof. George M. Sheldrick FRS Dept. Structural Chemistry, University of Goettingen, Tammannstr. 4, D37077 Goettingen, Germany Tel. +49-551-39-3021 or -3068 Fax. +49-551-39-22582
