Studies on the conformation of amino acids XII. Energy calculations on prolyl residue

The favoured conformations of the prolyl residue have been obtained by calculating their potential energies arising from bond-angle strain, torsion-angle strain, non-bonded and electrostatic interatomic energies. In addition to the five membered ring, the peptide unit at the amino end (with ω = 180°) and the C′ atom at the carboxyl end have been taken into account. It is found that there are two local minima in the configurational space of the parameters defining the conformation, as is actually observed-one (denoted by B) with Cγ displaced on the same side as C′, which is lower in energy than the other (denoted by A) with Cγ displaced on the opposite side of C′. The other four atoms Cδ, N, Cα, Cβ are nearly in a plane. The conformations of minimum energy (for both A and B) have bond angles very close to the mean observed values while the torsion angles are well within the range observed in various structures for each type. Taking into account the fact that the influence of neighbouring molecules in a crystal structure may make the conformation of a molecule different from the minimal one, the ranges of the conformational parameters for which the energy is within 0.6 kcal/mole above the minimum value (called the "most probable range") and within 1.2 kcal/mole (called the "probable range") have been determined. The ranges thus obtained, agree well with observation, and most of the observed data lie within the most probable ranges, although differing appreciably from the conformation of minimum energy. The study has been extended, in a limited way, to the conformation of the ring in the amino acid proline. Since the nitrogen is tetrahedral in this (as contrasted with being planar in the prolyl residue), it is found that any one of the five atoms can be out of plane (either way), with the other four lying nearly in a plane. These correspond to low energy conformations (up to 1.2 kcal/mole above the minimum). One such example, in which the Cα atom is out of plane is known for dl-proline · HCl. It is also shown that in these calculations energies due to bond length distortions can be neglected to a good degree of approximation, provided the 'best' values of the bond lengths for the particular compound are used in the theoretical calculations.

Electronic structures of electron donor-acceptor complexes: results from ultraviolet photoelectron spectroscopy and molecular orbital calculations

HeI photoelectron spectra of 1:1 electron donor-acceptor complexes are discussed in the light of molecular orbital calculations. The complexes discussed include those formed by BH3, BF3 and SO2. Some systematics have been found in the ionization energy shifts of the complexes compared to the free components and these are related to the strength of the donor-acceptor bond. Hel spectra of hydrogen bonded complexes are discussed in comparison with results from MO calculations. Limitations of such studies as well as scope for further investigations are indicated.

Late Svecofennian leucogranites of southern Finland : Chronicles of an orogenic collapse

Leucogranite magmatism occurred in southern Finland during the later stages of the Paleoproterozoic Svecofennian orogeny. The leucogranites are considered to have formed from pre-existing crustal rocks that have undergone anatexis in the extensional stage of the orogeny, following continental collision and resultant crustal thickening. The leucogranites have been studied in the field using petrographic and mineralogical methods, elemental and isotope geochemistry on whole rocks and minerals, and U-Pb geochronology. On outcrop scale, these granites typically form heterogeneous, layered, sheet-like bodies that migmatize their country rocks. All of the leucogranites are peraluminous and rich in SiO2, but otherwise display significant geochemical variation. Their Nd isotope composition ranges from fairly juvenile to very unradiogenic, and the Hf isotope composition of their zircon shows a varying degree of mixing in the source, the zircon populations becoming more heterogeneous and generally less radiogenic towards the east. The leucogranites have been dated using U-Pb isotopic analyses, utilizing thermal ionization mass spectrometry, secondary ion mass spectrometry, and laser ablation multicollector ICP mass spectrometry on zircon and monazite. The results show that the granites were emplaced between 1.85 Ga and 1.79 Ga, which is a considerably longer period than has traditionally been perceived for these rocks. The rocks tend to become younger towards the east. Single crystal data also display a wide array of inherited zircons, especially in the eastern part of the leucogranite belt. The most common inherited age groups are ~2.8 2.5 Ga, ~2.1 2.1 Ga, and ~1.9 Ga. Magmatic zircon and monazite usually record similar ages for any one sample.Thermobarometric calculations indicate that the leucogranites in the Veikkola area of southcentral Finland were formed from relatively low-temperature melts, and emplaced at 17-25 km depth, i.e. at mid-crustal level. It is likely that these conditions apply to the Svecofennian leucogranites in general. Large differences in the Hf and Nd isotope compositions, emplacement ages, and distributions of inherited zircon ages show that these granites were formed from different types of source rocks, which probably included both sedimentary and igneous rocks.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

Activation barriers for d.c. conductivity in ionic glasses: Classical calculations using the small-cluster approximation

A small-cluster approximation has been used to calculate the activation barriers for the d.c. conductivity in ionic glasses. The main emphasis of this approach is on the importance of the hitherto ignored polarization energy contribution to the total activation energy. For the first time it has been demonstrated that the d.c. conductivity activation energy can be calculated by considering ionic migration to a neighbouring vacancy in a smali cluster of ions consisting of face-sharing anion polyhedra. The activation energies from the model calculations have been compared with the experimental values in the case of highly modified lithium thioborate glasses.

Chiral polyhydroxylated tetrahydrothiophene derivatives: novel synthesis and structural elucidation by X-ray crystallography, NMR spectroscopy and molecular mechanics calculations

The dideoxygenation reaction of 1,3;4,6-di-O-alkylidene-2,5-di-S-methylthiocarbonyl-D-mannitol derivatives under Barton-McCombie reaction conditions gave the hexahydrodipyranothiophenes 4 and 7 instead of the expected 2,5-dideoxy products. Structural and conformational information on these novel derivatives has been obtained by NMR spectroscopy, single-crystal X-ray crystallography and molecular mechanics calculations.

Size dependence of the bulk modulus of semiconductor nanocrystals from first-principles calculations

The variation in the bulk modulus of semiconductor nanoparticles has been studied within first-principles electronic-structure calculations using the local density approximation (LDA) for the exchange correlation. Quantum Monte Carlo calculations carried out for a silicon nanocrystal Si87H76 provided reasonable agreement with the LDA results. An enhancement was observed in the bulk modulus as the size of the nanoparticle was decreased, with modest enhancements being predicted for the largest nanoparticles studied here, a size just accessible in experiments. To access larger sizes, we fit our calculated bulk moduli to the same empirical law for all materials, the asymptote of which is the bulk value of the modulus. This was found to be within 2-10% of the independently calculated value. The origin of the enhancement has been discussed in terms of Cohen's empirical law M.L. Cohen, Phys. Rev. B 32, 7988 (1985)] as well as other possible scenarios.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.