Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Variational Calculation of the Hydrogen Molecular Cation Using Maple

The variational calculation of the energy of the Hydrogen Molecular (H2+) Cation's LCAO-MO for the sigma and sigma* states as functions of the AO's screening constant and the internuclear distance is carried out explicitly in great detail.

Variational Calculation of the Hydrogen Molecular Cation using Maple (II). The pi and pi-star Orbitals

The pi and pi-star orbitals of the hydrogen molecular cation are obtained using Maple in the same manner as the sigma and sigma-star orbitals were obtained in paper-36.

CIPWFULL: A FORTRAN Program for Calculation of Comprehensive CIPW Norms of Igneous Rocks

CIPWFULL is a user-friendly, stand-alone FORTRAN software program that is designed to calculate the comprehensive CIPW normative mineral composition of igneous rocks and strictly adheres to the original formulation of the CIPW protocol. This faithful adherence alleviates inaccuracies in normative mineral calculations by programs commonly used by petrologists. Additionally, several of the most important petrological and mineralogical parameters of igneous rocks are calculated by the program. Along with all the regular major oxide elements, all the significant minor elements whose contents can potentially effect the CIPW normative mineral composition are included. CIPWFULL also calculates oxidation ratios for igneous rock samples that have only one oxidation state of iron reported in the specimen analysis. It also provides an option for normalization of analyses to unity on a hydrous-free basis in order to facilitate comparison of norms among rock groups. Other capabilities of the program cater for rare situations, like the presence of cancrinite or exclusion from the norm calculation of rare rocks like carbonatite. Several mineralogical, petrological and discriminatory parameters and indexes are additionally calculated by the CIPWFULL program. The CIPWFULL program is very efficient and flexible and allows for a user-defined free-format input of all the chemical species, and it permits feeding of minor elements as parts per million or oxide percentages. Results of calculations are printed in a formatted ASCII text file and may be optionally casted into a space-delimited text files that are ready to be imported to general spreadsheet programs. CIPWFULL is DOS-based and is implemented on WINDOWS and mainframe platforms.

Tab. 1: Approximate dimension of the pingos

A group of nine pingos occurs in the valley of a glacial meltwater river. The pingos rise from a plain of low-center polygons. Some pingos have a typical cone shape, but others are linear, apparently centered on ice wedges . The occurrence of most pingos at the junction of oversize ice wedge polygon ridges suggests that the injection of water and the segregation of ice occurred along pathways provided by the ice wedges.