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.

The normal co-ordinate analysis, vibrational spectra and theoretical infrared intensities of some thiocarbonyl halide molecules /

TITLE: The normal co-ordinate analysis, vibrational spectra and theoretical infrared intensities of some thiocarbonyl halides. AUTHOR: J. L. Brema SUPERVISOR: Dr. D. C. Moule NUMBER OF PAGES: 89 ABSTRACT: The vibrational assignment of the five-in-plane fundamental modes of CSClBr has been made on the basis of infrared gas phase and liquid Raman spectral analyses to supplement our earlier vibrational studies. Even though the one out-of-plane fundamental was not observed spectroscopically an attempt has been made to predict its frequency. The vibrational spectra contained impurity bands and the CSClBr assignment was made only after a thorough analysis of the impurities themselves. A normal co-ordinate analysis calculation was performed assuming a Urey-Bradley force field. This calculation yielded the fundamental frequencies in good agreement with those observed after refinement of the originally transferred force constants. The theoretical frequencies are the eigenvalues of the secular equation and the calculation also gave the corresponding eigenvectors in the form of the very important LLj matrix. The [l] matrix is the transfoirmation between internal co-ordinates and normal co-ordinates and it is essential for Franck-Condon calculations on electronically excited molecules and for infrared Integrated band intensity studies. Using a self-consistent molecular orbital calculation termed "complete neglect of differential overlap" (CNDO/2) , theoretical values of equilibrium bond lengths and angleswere calcuted for a series of carbonyl and thlocarbonyl molecules. From these calculations valence force field force constants were also determined but with limited success. With the CNIX)/2 method theoretical dipole moment derivatives with respect to symmetrized internal co-ordinates were calculated and the results should be useful in a correlation with experimentally determined values.

Interaction between vortices and impurities in a d-wave superconductor /

High temperature superconductors were discovered in 1986, but despite considerable research efforts, both experimental and theoretical, these materials remain poorly understood. Because their electronic structure is both inhomogeneous and highly correlated, a full understanding will require knowledge of quasiparticle properties both in real space and momentum space. In this thesis, we will present a theoretical analysis of the scanning tunneling microscopy (STM) data in BSCCO. We introduce the Bogoliubov-De Gennes Hamiltonian and solve it numerically on a two-dimensional 20 x 20 lattice under a magnetic field perpendicular to the surface. We consider a vortex at the center of our model. We introduce a Zn impurity in our lattice as a microscopic probe of the physical properties of BSCCO. By direct numerical diagonalization of the lattice BogoliubovDe Gennes Hamiltonian for different positions of the impurity, we can calculate the interaction between the vortex and the impurity in a d-wave superconductor.

Quantum Monte Carlo : some theoretical and numerical studies

In Part I, theoretical derivations for Variational Monte Carlo calculations are compared with results from a numerical calculation of He; both indicate that minimization of the ratio estimate of Evar , denoted EMC ' provides different optimal variational parameters than does minimization of the variance of E MC • Similar derivations for Diffusion Monte Carlo calculations provide a theoretical justification for empirical observations made by other workers. In Part II, Importance sampling in prolate spheroidal coordinates allows Monte Carlo calculations to be made of E for the vdW molecule var He2' using a simplifying partitioning of the Hamiltonian and both an HF-SCF and an explicitly correlated wavefunction. Improvements are suggested which would permit the extension of the computational precision to the point where an estimate of the interaction energy could be made~

Rotation-electronic interaction in the 3p-complex Rydberg state of water

An energy theory is formulated for the rotational energy levels in a p-complex Rydberg state of an asymmetric top molecule of symmetry C2v. The effective Hamiltonian used consists of the usual rigid rotor Hamiltonian augmented with terms representing electronic spin and orbital angular momentum effects. Criteria for assigning symmetry species to the rotational energy levels, following Houganfs scheme that uses the full molecular group,are established and given in the form of a table. This is particularly suitable when eigenvectors are calculated on a digital computer. Also, an intensity theory for transitions to the Rydberg p-complex singlet states is presented and selection rules in terms of symmetry species of energy states are established. Finally, applications to HpO and DpO are given.

Exact property estimation from diffusion Monte Carlo with minimal stochastic reconfiguration /

Our objective is to develop a diffusion Monte Carlo (DMC) algorithm to estimate the exact expectation values, ($o|^|^o), of multiplicative operators, such as polarizabilities and high-order hyperpolarizabilities, for isolated atoms and molecules. The existing forward-walking pure diffusion Monte Carlo (FW-PDMC) algorithm which attempts this has a serious bias. On the other hand, the DMC algorithm with minimal stochastic reconfiguration provides unbiased estimates of the energies, but the expectation values ($o|^|^) are contaminated by ^, an user specified, approximate wave function, when A does not commute with the Hamiltonian. We modified the latter algorithm to obtain the exact expectation values for these operators, while at the same time eliminating the bias. To compare the efficiency of FW-PDMC and the modified DMC algorithms we calculated simple properties of the H atom, such as various functions of coordinates and polarizabilities. Using three non-exact wave functions, one of moderate quality and the others very crude, in each case the results are within statistical error of the exact values.

An analysis of the methyl rotation dynamics in the So (x' A) and T (a A) states of thioacetaldehyde from Pyrolysis jet spectra /

Jet-cooled, laser-induced phosphorescence excitation spectra (LIP) of thioacetaldehyde CH3CHS, CH3CDS, CD3CHS and CD3CDS have been observed over the region 15800 - 17300 cm"^ in a continuous pyrolysis jet. The vibronic band structure of the singlet-triplet n -* n* transition were attributed to the strong coupling of the methyl torsion and aldehydic hydrogen wagging modes . The vibronic peaks have been assigned in terms of two upper electronic state (T^) vibrations; the methyl torsion mode v^g, and the aldehydic hydrogen wagging mode v^^. The electronic origin O^a^ is unequivocally assigned as follows: CH3CHS (16294.9 cm"'' ), CH3CDS (16360.9 cm"'' ), CD3CHS (16299.7 cm"^ ), and CD3CDS (16367.2 cm"'' ). To obtain structural and dynamical information about the two electronic states, potential surfaces V(e,a) for the 6 (methyl torsion) and a (hydrogen wagging) motions were generated by ab initio quantum mechanical calculations with a 6-3 IG* basis in which the structural parameters were fully relaxed. The kinetic energy coefficients BQ(a,e) , B^(a,G) , and the cross coupling term B^(a,e) , were accurately represented as functions of the two active coordinates, a and 9. The calculations reveal that the molecule adopts an eclipsed conformation for the lower Sq electronic state (a=0°,e=0"') with a barrier height to internal rotation of 541.5 cm"^ which is to be compared to 549.8 cm"^ obtained from the microwave experiment. The conformation of the upper T^ electronic state was found to be staggered (a=24 . 68° ,e=-45. 66° ) . The saddle point in the path traced out by the aldehyde wagging motion was calculated to be 175 cm"^ above the equilibrium configuration. The corresponding maxima in the path taken by methyl torsion was found to be 322 cm'\ The small amplitude normal vibrational modes were also calculated to aid in the assignment of the spectra. Torsional-wagging energy manifolds for the two states were derived from the Hamiltonian H(a,e) which was solved variationally using an extended two dimensional Fourier expansion as a basis set. A torsionalinversion band spectrum was derived from the calculated energy levels and Franck-Condon factors, and was compared with the experimental supersonic-jet spectra. Most of the anomalies which were associated with the interpretation of the observed spectrum could be accounted for by the band profiles derived from ab initio SCF calculations. A model describing the jet spectra was derived by scaling the ab initio potential functions. The global least squares fitting generates a triplet state potential which has a minimum at (a=22.38° ,e=-41.08°) . The flatter potential in the scaled model yielded excellent agreement between the observed and calculated frequency intervals.

Application of Reptation Quantum Monte Carlo and Related Methods to LiH

This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.

Letter - Mrs. Olive Diefenbaker to Sean O'Sullivan, 6 December 1966

A handwritten letter from Olive Diefenbaker to Sean O'Sullivan, 6 December 1966. Mrs. Diefenbaker remarks about a phone call she received from a Hamiltonian who had met O'Sullivan at "the Convention" and had "'tried to help him over his disillusion"'. Mrs. Diefenbaker took offense to the comments and responded to the individual "I wouldn't worry about Sean. He is mature and wise, and will think his way through like the rest of us".