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.

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 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.

The crystal and molecular structure of bis (pyridoxamine) copper (II) dinitrate monohydrate

The crystal structure of Cu(PM)2(N03hoH20 (where PM is pyridoxamine, CSHI2N202) has been determined from three dimensional x-ray diffraction data. The crystals are triclinic, space group pI, a = 14.248 (2), b = 8.568 (1), c = 9.319 (1) 1, a = 94.08 (1), e = 89.73 (1), y~~ 99.18 (1)°, z = 2, jl(MoK) = 10.90 em-I, Po = 1.61 g/cm3 and Pc = 1.61 g/em3• The structure a was solved by Patterson techniques from data collected on a Picker 4-circle diffractometer to 26max = 45°. All atoms, including hydrogens, have been located. Anisotropic thermal parameters have been refined for all nonhydrogen atoms. For the 2390 independent reflections with F ? 3cr(F) , R = 0.0408. The results presented here provide the first detailed structural information of a metal complex with PM itself. The copper atoms are located on centres of symmetry and each is chela ted by two PM zwitterions through the amino groups and phenolate oxygen atoms. The zwitterionic form found in this structure involves the loss of a proton from the phenolate group and protonation of the pyridine ring nitrogen atoms. The two independent Cu(PM)2 moieties are symmetrically bridged by a single oxygen atom from one of the nitrate groups. The second nitrate group is not coordinated to the copper atoms but is central to an extensive hydrogen bonding network involving the water molecule and uncoordinated functional groups of PM.

The crystal and molecular structure of thiamine hydroiodide /|nWilliam Edward Lee. -- 260 St. Catharines [Ont.] : Dept. of Chemistry, Brock University,

The x-ray crystal structure of thiamine hydroiodide,C1ZH18N40S12' has been determined. The unit cell parameters are a = 13.84 ± 0.03, o b = 7.44 ± 0.01, c = 20.24 ± 0.02 A, 8 = 120.52 ± 0.07°, space group P2/c, z = 4. A total of 1445 reflections having ,2 > 2o(F2), 26 < 40° were collected on a Picker four-circle diffractometer with MoKa radiation by the 26 scan technique. The structure was solved by the heavy atom method. The iodine and sulphur atoms were refined anisotropically; only the positional parameters were refined for the hydrogen atoms. Successive least squares cycles yielded an unweighted R factor of 0.054. The site of protonation of the pyrimidine ring is the nitrogen opposite the amino group. The overall structure conforms very closely to the structures of other related thiamine compounds. The bonding surrounding the iodine atoms is distorted tetrahedral. The iodine atoms make several contacts with surrounding atoms most of them at or near the van der Waal's distances A thiaminium tetrachlorocobaltate salt was produced whose molecular and crystal structure was j~dged to be isomorphous to thiaminium tetrachlorocadmate.

The crystal and molecular structure of 18-Crown-6 HgC12 and 18-Crown-6 Cdc12

Hg(18-Crown-6)C12 and Cd(18-Crown-6)C12 are isostructura1, space group Cl~ Z = 2. For the mercury compound, a = 10.444(2) A° , b = 11. 468(1) A° , c = 7.754(1) A° , a = 90.06(1)°, B = 82.20(1)°, Y = 90.07(1)°, Dobs = 1.87, Dca1c = 1.93, V = 920.05 13, R = 4.66%. For the cadmium compound, 000 a = 10.374(1) A, b = 11.419(2) A, c = 7.729(1) A, a = 89.95(1)°, B = 81.86(2)°, Y = 89.99(1)°, Dobs = 1.61, Dcalc = 1.64, V = 906.4613, R = 3.95%. The mercury and cadmium ions exhibit hexagonal bipyramidal coordination, with the metal ion located on a centre of symmetry in the plane of the oxygen atoms. The main differences between the two structures are an increase in the metal-oxygen distance and a reduction in the metalchloride distance when the central ion changes from Cd2+ to Hg2+. These differences may be explained in terms of the differences in hardness or softness of the metal ions and the donor atoms.