Variational Monte Carlo study of core-valence separation schemes for first-row atoms and positive ions /

All-electron partitioning of wave functions into products ^core^vai of core and valence parts in orbital space results in the loss of core-valence antisymmetry, uncorrelation of motion of core and valence electrons, and core-valence overlap. These effects are studied with the variational Monte Carlo method using appropriately designed wave functions for the first-row atoms and positive ions. It is shown that the loss of antisymmetry with respect to interchange of core and valence electrons is a dominant effect which increases rapidly through the row, while the effect of core-valence uncorrelation is generally smaller. Orthogonality of the core and valence parts partially substitutes the exclusion principle and is absolutely necessary for meaningful calculations with partitioned wave functions. Core-valence overlap may lead to nonsensical values of the total energy. It has been found that even relatively crude core-valence partitioned wave functions generally can estimate ionization potentials with better accuracy than that of the traditional, non-partitioned ones, provided that they achieve maximum separation (independence) of core and valence shells accompanied by high internal flexibility of ^core and Wvai- Our best core-valence partitioned wave function of that kind estimates the IP's with an accuracy comparable to the most accurate theoretical determinations in the literature.

Diffusion Monte Carlo study of electronic properties for H and Be atoms /

We examined three different algorithms used in diffusion Monte Carlo (DMC) to study their precisions and accuracies in predicting properties of isolated atoms, which are H atom ground state, Be atom ground state and H atom first excited state. All three algorithms — basic DMC, minimal stochastic reconfiguration DMC, and pure DMC, each with future-walking, are successfully impletmented in ground state energy and simple moments calculations with satisfactory results. Pure diffusion Monte Carlo with future-walking algorithm is proven to be the simplest approach with the least variance. Polarizabilities for Be atom ground state and H atom first excited state are not satisfactorily estimated in the infinitesimal differentiation approach. Likewise, an approach using the finite field approximation with an unperturbed wavefunction for the latter system also fails. However, accurate estimations for the a-polarizabilities are obtained by using wavefunctions that come from the time-independent perturbation theory. This suggests the flaw in our approach to polarizability estimation for these difficult cases rests with our having assumed the trial function is unaffected by infinitesimal perturbations in the Hamiltonian.

Studies of some cocondensations of nickel atoms with unsaturated organic ligands and with mixed ligand systems

The cocondensation of nickel with a number of unsaturated ligands was studied, as was the cocondensation with a number of mixed ligand systems. Enamines were found not to react with nickel while acrylonitrile was polymerized. In the mixed ligand syst.ems different products were obtained than when the ligands were cocondensed individually. Cocondensations of benzyl halide/allyl halide mixtures gave unstable products that were not observed when the halides were cocondensed individually. The effect of Kao-Wool insulation on nickel/benzyl halide cocondensations was found to be significant. Kao-Wool caused the bulk of the benzyl halide to be polymeri zed to a number of poly-benzylic species. An alkali metal reactor was designed for the evaporation of sodium and potassium atoms into cold solutions of metal halide and an or ganic substrate. This apparatus was used to synthesize Ni(P¢3 )3' but proved unsuccessful for synthesizing a nickel-enamine compound.

Quantum Monte Carlo study of electrostatic polarizabilities of H and He atoms /

The infinitesimal differential quantum Monte Carlo (QMC) technique is used to estimate electrostatic polarizabilities of the H and He atoms up to the sixth order in the electric field perturbation. All 542 different QMC estimators of the nonzero atomic polarizabilities are derived and used in order to decrease the statistical error and to obtain the maximum efficiency of the simulations. We are confident that the estimates are "exact" (free of systematic error): the two atoms are nodeless systems, hence no fixed-node error is introduced. Furthermore, we develope and use techniques which eliminate systematic error inherent when extrapolating our results to zero time-step and large stack-size. The QMC results are consistent with published accurate values obtained using perturbation methods. The precision is found to be related to the number of perturbations, varying from 2 to 4 significant digits.

Optical Properties of Sb2Te3, and Dilute Magnetic Semiconductors Sb1.97 VO.03 Te3 and Sb1.94 CrO.06 Te3

This thesis reports on the optical properties of the dilute magnetic semiconductors, Sb1.97 V 0.03 Te3 and Sb1.94Cr0.06Te3, along with the parent compound Sb2Te3' These materials develop a ferromagnetic state at low temperature with Curie temperatures of 22 K and 16 K respectively. All three samples were oriented such that the electric field vector of the light was perpendicular to the c-axis. The reflectance profile of these samples in the mid-infrared (500 to 3000 cm-1) shows a pronounced plasma edge which retracts with decreasing temperature. The far-infrared region of these samples exhibits a phonon at ~ 60 cm-1 which softens as temperature decreases. Kramers-Kronig analysis and a Drude-Lorentz model were employed to determine the optical constants of the bulk samples. The real part of the optical conductivity is shown to consist of intraband contributions at frequencies below the energy gap (~0.26 eV) and interband contributions at frequencies above the energy gap. The temperature dependence of the scattering rate show that a mix of phonon and impurity scattering are present, while the signature of traditional spin disorder (magnetic) scattering was difficult to confirm.

Calculation of nondifferential properties for atomic ground states /

A new method for sampling the exact (within the nodal error) ground state distribution and nondiflPerential properties of multielectron systems is developed and applied to firstrow atoms. Calculated properties are the distribution moments and the electronic density at the nucleus (the 6 operator). For this purpose, new simple trial functions are developed and optimized. First, using Hydrogen as a test case, we demonstrate the accuracy of our algorithm and its sensitivity to error in the trial function. Applications to first row atoms are then described. We obtain results which are more satisfactory than the ones obtained previously using Monte Carlo methods, despite the relative crudeness of our trial functions. Also, a comparison is made with results of highly accurate post-Hartree Fock calculations, thereby illuminating the nodal error in our estimates. Taking into account the CPU time spent, our results, particularly for the 8 operator, have a relatively large variance. Several ways of improving the eflSciency together with some extensions of the algorithm are suggested.

The osmoadaptive response of the wine yeast Saccharomyces cerevisiae K1-V1116 during icewine fermentation

The adapted metabolic response of commercial wine yeast under prolonged exposure to concentrated solutes present in Icewine juice is not fully understood. Presently, there is no information regarding the transcriptomic changes in gene expression associated with the adaptive stress response ofwine yeast during Icewine fermentation compared to table wine fermentation. To understand how and why wine yeast respond differently at the genomic level and ultimately at the metabolic level during Icewine fermentation, the focus ofthis project was to identify and compare these differences in the wine yeast Saccharomyces cerevisiae KI-Vll16 using cDNA microarray technology during the first five days of fermentation. Significant differences in yeast gene expression patterns between fermentation conditions were correlated to differences in nutrient utilization and metabolite production. Sugar consumption, nitrogen usage and metabolite levels were measured using enzyme assays and HPLC. Also, a small subset of differentially expressed genes was verified using Northern analysis. The high osmotic stress experienced by wine yeast throughout Icewine fermentation elicited changes in cell growth and metabolism correlating to several fermentation difficulties, including reduced biomass accumulation and fermentation rate. Genes associated with carbohydrate and nitrogen transport and metabolism were expressed at lower levels in Icewine juice fermenting cells compared to dilute juice fermenting cells. Osmotic stress, not nutrient availability during Icewine fermentation appears to impede sugar and nitrogen utilization. Previous studies have established that glycerol and acetic acid production are increased in yeast during Icewine fermentation. A gene encoding for a glycerollW symporter (STL1) was found to be highly expressed up to 25-fold in the i Icewine juice condition using microarray and Northern analysis. Active glycerol transport by yeast under hyperosmotic conditions to increase cytosolic glycerol concentration may contribute to reduced cell growth observed in the Icewine juice condition. Additionally, genes encoding for two acetyl CoA synthetase isoforms (ACSl and ACS2) were found to be highly expressed, 19- and II-fold respectively, in dilute juice fermenting cells relative to the Icewine juice condition. Therefore, decreased conversion of acetate to acetyl-CoA may contribute to increased acetic acid production during Icewine fermentation. These results further help to explain the response of wine yeast as they adapt to Icewine juice fermentation. ii

The structural effects of divalent cations and insulin on phospholipid model membranes /

It is well accepted that structural studies with model membranes are of considerable value in understanding the structure of biological membranes. Many studies with models of pure phospholipids have been done; but the effects of divalent cations and protein on these models would make these studies more applicable to intact membrane. The present study, performed with above view, is a structural analysis of divalent io~cardio1ipin complexes using the technique of x-ray diffraction. Cardiolipin, precipitated from dilute solution by divalent ionscalcium, magnesium and barium, contains little water and the structure formed is similar to the structure of pure cardiolipin with low water content. The calcium-cardiolipin complex forms a pure hexagonal type II phase that exists from 40 to 400 C. The molar ratio of calcium and cardiolipin in the complex is 1 : 1. Cardiolipin, precipitated with magnesium and barium forms two co-existing phases, lamellar and hexagonal, the relative quantity of the two phases being dependent on temperature. The hexagonal phase type II consisting of water filled channels formed by adding calcium to cardiolipin may have a remarkable permeability property in intact membrane. Pure cardiolipin and insulin at pH 3.0 and 4.0 precipitate but form no organised structure. Lecithin/cardiolipin and insulin precipitated at pH 3.0 give a pure lamellar phase. As the lecithin/cardiolipin molar ratio changes from 93/7 to SO/50, (a) the repeat distance of the lamellar changes from 72.8 X to 68.2 A; (b) the amount of protein bound increases in such a way that cardiolipin/insulin molar ratio in the complex reaches a maximum constant value at lecithin/cardiolipin molar ratio 70/30. A structural model based on these data shows that the molecular arrangement of lipid and protein is a lipid bilayer coated with protein molecules. The lipid-protein interaction is chiefly electrostatic and little, if any, hydrophobic bonding occurs in this particular system. So, the proposed model is essentially the same as Davson-Daniellifs model of biological membrane.

The design, synthesis and characterization of new building blocks for the preparation of molecule-based magnetic materials /

Two new families of building blocks have been prepared and fully characterized and their coordination chemistry exploited for the preparation of molecule-based magnetic materials. The first class of compounds were prepared by exploiting the chemistry of 3,3'-diamino-2,2'-bipyridine together with 2-pyridine carbonyl chloride or 2-pyridine aldehyde. Two new ligands, 2,2'-bipyridine-3,3'-[2-pyridinecarboxamide] (Li, 2.3) and N'-6/s(2-pyridylmethyl) [2,2'bipyridine]-3,3'-diimine (L2, 2.7), were prepared and characterized. For ligand L4, two copper(II) coordination compounds were isolated with stoichiometrics [Cu2(Li)(hfac)2] (2.4) and [Cu(Li)Cl2] (2.5). The molecular structures of both complexes were determined by X-ray crystallography. In both complexes the ligand is in the dianionic form and coordinates the divalent Cu(II) ions via one amido and two pyridine nitrogen donor atoms. In (2.4), the coordination geometry around both Cu11 ions is best described as distorted trigonal bipyramidal where the remaining two coordination sites are satisfied by hfac counterions. In (2.5), both Cu(II) ions adopt a (4+1) distorted square pyramidal geometry. One copper forms a longer apical bond to an adjacent carbonyl oxygen atom, whereas the second copper is chelated to a neighboring Cu-Cl chloride ion to afford chloride bridged linear [Cu2(Li)Cl2]2 tetramers that run along the c-axis of the unit cell. The magnetic susceptibility data for (2.4) reveal the occurrence of weak antiferromagnetic interactions between the copper(II) ions. In contrast, variable temperature magnetic susceptibility measurements for (2.5) reveal more complex magnetic properties with the presence of ferromagnetic exchange between the central dimeric pair of copper atoms and weak antiferromagnetic exchange between the outer pairs of copper atoms. The Schiff-base bis-imine ligand (L2, 2.7) was found to be highly reactive; single crystals grown from dry methanol afforded compound (2.14) for which two methanol molecules had added across the imine double bond. The susceptibility of this ligand to nucleophilic attack at its imine functionality assisted via chelation to Lewis acidic metal ions adds an interesting dimension to its coordination chemistry. In this respect, a Co(II) quaterpyridine-type complex was prepared via a one-pot transformation of ligand L2 in the presence of a Lewis acidic metal salt. The rearranged complex was characterized by X-ray crystallography and a reaction mechanism for its formation has been proposed. Three additional rearranged complexes (2.13), (2.17) and (2.19) were also isolated when ligand (L2, 2.7) was reacted with transition metal ions. The molecular structures of all three complexes have been determined by X-ray crystallography. The second class of compounds that are reported in this thesis, are the two diacetyl pyridine derivatives, 4-pyridyl-2,6-diacetylpyridine (5.5) and 2,2'-6,6'-tetraacetyl-4,4'-bipyridine (5.15). Both of these compounds have been designed as intermediates for the metal templated assembly of a Schiff-base N3O2 macrocycle. From compound (5.15), a covalently tethered dimeric Mn(II) macrocyclic compound of general formula {[Mn^C^XJCl-FkO^Cl-lO.SFbO (5.16) was prepared and characterized. The X-ray analysis of (5.16) reveals that the two manganese ions assume a pentagonal-bipyramidal geometry with the macrocycle occupying the pentagonal plane and the axial positions being filled by a halide ion and a H2O molecule. Magnetic susceptibility data reveal the occurrence of antiferromagnetic interactions between covalently tethered Mn(II)-Mn(II) dimeric units. Following this methodology a Co(II) analogue (5.17) has also been prepared which is isostructural with (5.16).

An experimental investigation of dechanneling in copper single crystals

This investigation comprises a comparison of experimental and theoretical dechanneling of MeV protons in copper single crystals. Dechanneling results when an ion's transverse energy increases to the value where the ion can undergo small impact parameter collisions with individual atoms. Depth dependent dechanneling rates were determined as functions of lattice temperature, ion beam energy and crystal axis orientation. Ion beam energies were IMeV and 2MeV,temperatures ranged from 35 K to 280 K and the experiment was carried out along both the (lOa) and <110) axes. Experimental data took the form of aligned and random Rutherford backscattered energy spectra. Dechanneling rates were extracted from these spectra using a single scattering theory that took explicit account of the different stopping powers experienced by channeled and dechanneled ions and also included a correction factor to take into account multiple scattering effects along the ion's trajectory. The assumption of statistical equilibrium and small angle scattering of the channeled ions allows a description of dechanneling in terms of the solution of a diffusion like equation which contains a so called diffusion function. The diffusion function is shown to be related to the increase in average transverse energy. Theoretical treatments of increase in average transverse energy due to collisions of projectiles with channel electrons and thermal perturbations in the lattice potential are reviewed. Using the diffusion equation and the electron density in the channel centre as a fitting parameter dechanneling rates are extracted. Excellent agreement between theory and experiment has been demonstrated. Electron densities determined in the fitting procedure appear to be realistic. The surface parameters show themselves to be good indicators of the quality of the crystal.

Investigation of molecular polarizabilities and derivatives in halomethanes

The research undertaken was to obtain absolute Raman intensities for the symmetric stretching vibrations of the methyl halides, CH3X with (X=F, CI, Br), by experiment and theory. The intensities were experimentally measured using the Ar+ ion gas laser as excitation source, a Spex 14018 double monochromator and a RCA C-31034 photomultiplier tube as detector. These intensities arise from changes in the derivative of the polarizability (8 a'), with respect to vibration along a normal coordinate (8qi). It was intended that these derivatives obtained with respect to normal coordinates would be converted to derivatives with respect to internal coordinates, for a quantitative comparison with theory. Theoretical numerical polarizability derivatives for the stretching vibrations are obtained using the following procedure. A vibration was simulated in the molecule by increasi.ng and decreasing the respective bond by the amount ±o.oosA for the C-H bonds and ±o.oIA for the C-X (X=F, CI, Br) bond. The derivative was obtained by taking the difference in the polarizability for the equilibrium geometry and the geometry when a particular bond is changed. This difference, when divided by the amount of change in each bond and the number of bonds present results in the derivative of the polarizability with respect to internal coordinate i.e., !1u/!1r. These derivatives were obtained by two methods: I} ab initio molecular orbital calculation and 2} theory of atoms in molecules (AIM) analysis. Due to errors in the experimental setup only a qualitative analysis of the results was undertaken relative to the theory. Theoretically it is predicted that the symmetric carbonhalogen stretch vibrations are more intense than the respective carbon-hydrogen stretch, but only for the methyl chloride and bromide. The carbon fluorine stretch is less intense than the carbon-hydrogen stretch, a fact which is attributed to the small size and high electronegativity of the fluorine atom. The experimental observations are seen to agree qualitatively with the theory results. It is hoped that when the experiment is repeated, a quantitative comparison can be made. The analysis by the theory of atoms in molecules, along with providing polarizabilities and polarizability derivatives, gives additional information outlined below. The theory provides a pictorial description of the main factors contributing to the molecular polarizability and polarizability derivative. These contributions are from the charge transfer and atomic dipole terms i.e., transfer of charge from one atom to another and the reorganization of atomic electronic charge distribution due to presence of an electric field. The linear relationship between polarizability and molecular volume was also observed.

Modified statistical methods to calculate rare gas interaction potentials

New density functionals representing the exchange and correlation energies (per electron) are employed, based on the electron gas model, to calculate interaction potentials of noble gas systems X2 and XY, where X (and Y) are He,Ne,Ar and Kr, and of hydrogen atomrare gas systems H-X. The exchange energy density functional is that recommended by Handler and the correlation energy density functional is a rational function involving two parameters which were optimized to reproduce the correlation energy of He atom. Application of the two parameter function to other rare gas atoms shows that it is "universal"; i. e. ,accurate for the systems considered. The potentials obtained in this work compare well with recent experimental results and are a significant improvement over those from competing statistical modelS.

NMR studies of the exchange reactions of CH3CN.BX3 with excess CH3CN /|nJoseph Fogelman. -- 260 St. Catharines, Ont. : [s. n.],

Exchange reactions between molecular complexes and excess acid or base are well known and have been extensively surveyed in the literature(l). Since the exchange mechanism will, in some way involve the breaking of the labile donor-acceptor bond, it follows that a discussion of the factors relating to bonding in molecular complexes will be relevant. In general, a strong Lewis base and a strong Lewis acid form a stable adduct provided that certain stereochemical requirements are met. A strong Lewis base has the following characteristics (1),(2) (i) high electron density at the donor site. (ii) a non-bonded electron pair which has a low ionization potential (iii) electron donating substituents at the donor atom site. (iv) facile approach of the site of the Lewis base to the acceptor site as dictated by the steric hindrance of the substituents. Examples of typical Lewis bases are ethers, nitriles, ketones, alcohols, amines and phosphines. For a strong Lewis acid, the following properties are important:( i) low electron density at the acceptor site. (ii) electron withdrawing substituents. (iii) substituents which do not interfere with the close approach of the Lewis base. (iv) availability of a vacant orbital capable of accepting the lone electron pair of the donor atom. Examples of Lewis acids are the group III and IV halides such (M=B, AI, Ga, In) and MX4 - (M=Si, Ge, Sn, Pb). The relative bond strengths of molecular complexes have been investigated by:- (i) (ii) (iii) (iv) (v] (vi) dipole moment measurements (3). shifts of the carbonyl peaks in the IIIR. (4) ,(5), (6) .. NMR chemical shift data (4),(7),(8),(9). D.V. and visible spectrophotometric shifts (10),(11). equilibrium constant data (12), (13). heats of dissociation and heats of reactions (l~), (16), (17), (18), (19). Many experiments have bben carried out on boron trihalides in order to determine their relative acid strengths. Using pyridine, nitrobenzene, acetonitrile and trimethylamine as reference Lewis bases, it was found that the acid strength varied in order:RBx3 > BC1 3 >BF 3 • For the acetonitrile-boron trihalide and trimethylamine boron trihalide complexes in nitrobenzene, an-NMR study (7) showed that the shift to lower field was. greatest for the BB~3 adduct ~n~ smallest for the BF 3 which is in agreement with the acid strengths. If electronegativities of the substituents were the only important effect, and since c~ Br ,one would expect the electron density at the boron nucleus to vary as BF3atoms has been proposed as the predominating factor, i.e. ~rt- back-bond between a lone electron pair on the halogen and the vacant orbital on the boron site. The degree of back-bonding varies inversely as the bo~on halogen distance and one would therefore expect the B-F bond to exhibit greater back-bonding character than the B-Cl or B-Br bonds. Since back-bonding transfers electron density from substituent to the boron atom site, this process would be expected to weaken the Lewis acid strength. This explains the Lewis acid strength increasing in the order BF 3 BC1 3 BBr 3 . When the acetonitrile boron trihalide complex is formed, the boron atom undergoes ~_cbange of hybridization from sp2 to sp3. From a linear relationship between the heat of formation of ethyl acetate adducts and the shift in the carbonyl I.R. stretch, Drago (22) et al have proposed that the angular di~tortion of the X-B-X bonds from sp2 (12 ) to sp3 (10 hybridization is proportional to the amount of charge transferred, i.e. to the nature of the base, and they have rejected the earlier concept of reorganization energy in explaining the formation of the adduct bond (19).

Palladium and platinum complexes of vitamin Bâ

Palladium and platinum complexes of pyridoxamine, pyridoxine and pyridoxal have been prepared. The structures of the complexes PtCI2PM.H20, trans-PdC12 (PN)2 and [PLH+ ]2[PtC16] 2- ,H20 have been determined by use of single crystal x-ray studies. The compounds PdC12PH, trans-PdC12 (PN) 2 , cis-PdCI2 (PN)2 and cis PdC12 (PL)2 were also studied by use of carbon-13 nmr spectroscopy. All the complexes have also been characterised by use of infrared spectral studies. In the complexes, PtCI2PM.H20 and PdC12PM, the ligand pyridoxamine is chela ted to the metal through the aminomethyl nitrogen and the phenolate oxygen atoms whereas in the complexes, trans-PdCI2 (PN)2' cis-PdCI2 (PN)2 and cis-PdC12 (PL)2 the vitamin B6 ligands are coordinated to the metal through the pyridine ring nitrogen. The compounds [PLH+ ]2[PtCI6] 2- .H20 and [PMH2] 2+ [PdCI4] 2- .H20have no direct metal-ligand bonding, In all the complexes, the metal maintains a square planar coordination except in [PLH +] 2[PtCI6] 2- ,H20 where the metal is octahedrally coordinated. PH pyridoxamine [PMH ] 2+ = diprotonated pyridoxamine 2 PN = pyridoxine PL pyridoxal PLH+ protonated pyridoxal

The application of the structural correlation method to determine the reaction coordinate for a putative ring closure reaction

The cr ystal structure of the compound 2-benzoylethylidene-3-(2,4- dibromophenyl)-2,3-dihydro-5-phenyl-l,3,4-thiadiazole* C23H16Br2NZOS (BRMEO) has been determined by using three dimensiona l x-ray diffraction data. The crys tal form is monoclinic, space group P21/c, a = 17.492(4), o -.t' 0 R 0 b =: 16.979(1), c = 14.962(1) A, "X. =o= 90 ',= 106.46(1) , z = 8, graphite-monochromatized Mo~ rad iation, Jl= 0.710J3~, D = 1.62g/cc and o D = 1.65g/cc. The data were col lected on ~ Nonius CAD-4 c diffractometer. The following atoms were made anisotropic: Br, S, N, 0, C7, and C14-C16 for each i ndependent molecu le ; the rest were left isotropic. For 3112 independent refl ec tions with F > 6G\F), R == 0.057. The compound has two independent molecules within the asymmetric unit. Two different conformers were observed which pack well together. /l The S---O interaction distances of 2.493(6) and 2 . 478(7) A were observed for molecules A and B respectively. These values are consistent with earlier findings for 2-benzoylmethylene-3-(2,4-dibromophenyl)- ~~ 2,3-dihydro-5-phenyl-l,3,4-thiadiazole C22H14Br2N20S (BRPHO) and 2-benzoylpropylidene-3-(2,4-dibromophenyl)-2,3-dihydroiii ,'r 5-phenyl-l,3,4-thiadiazole C24H18Br2N20S (BRPETO ) where S---O distances are l ess than the van der Waals (3.251\) but greater than those expected for () a single bond (1.50A). From the results and the literature it appears obvious that the energy/reaction coordinate pathway has a minimum between the end structures (the mono- and bicyclic compounds). * See reference (21) for nomenclature.