Elliptical orbital studies of H21 and H2 molecules /|nS. K. Gupta. -- 260 St. Catharines, Ont. : [s. n.],

Calculations are performed on the \S <:Jd ground states of d ' + the H and HC) molecules using a basis set of non-integral ~ ~ I elliptical orbitals. Different variational wavefunctions constructed i- for H~ involved one parameter to three par~~eter variation. In order to l"'educe the ntunber of parameters in most commonly 0- used basis orbitals set, the importance of the term (,+~) Y\ over the term ;u 'Where n is a variational pararneter and the value of cr may be given by boundary condition or cusp condition is outlined in Chapters II and III. It is found that the two parameter -+ calculation for the ground state of the hydrogen molecule for R =1.4 (equilibrium internuclear separation) is perfonaed. The results are excellent.

Exchange energy calculation of He-He with supermolecular one-matrix

Exch~nge energy of the He-He system is calculated using the one-density matrix which has been modified according to the supermolecular density formula quoted by Kolos. The exchange energy integrals are computed analytically and by the Monte Carlo method. The results obtained from both ways compared favourably,with the results obtained from the SCF program HONDO

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.

Molecular dynamics calculation of mean square displacement in alkali metals and rare gas solids and comparison with lattice dynamics

Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.

On the formulation and the calculation of the harmonic contributions to the Debye-Waller factor in metals (sodium)

The algebraic expressions for the anharmonic contributions to the Debye-Waller factor up to 0(A ) and 0 L% ) £ where ^ is the scattering wave-vector] have been derived in a form suitable for cubic metals with small ion cores where the interatomic potential extends to many neighbours. This has been achieved in terms of various wave-vector dependent tensors, following the work of Shukla and Taylor (1974) on the cubic anharmonic Helmholtz free energy. The contribution to the various wave-vector dependent tensors from the coulomb and the electron-ion terms in the interatomic metallic potential has been obtained by the Ewald procedure. All the restricted multiple whole B r i l l o u i n zone (B.Z.) sums are reduced to single whole B.Z. sums by using the plane wave representation of the delta function. These single whole B.Z. sums are further reduced to the •%?? portion of the B.Z. following Shukla and Wilk (1974) and Shukla and Taylor (1974). Numerical calculations have been performed for sodium where the Born-Mayer term in the interatomic potential has been neglected because i t is small £ Vosko (1964)3 • *n o^er to compare our calculated results with the experimental results of Dawton (1937), we have also calculated the r a t io of the intensities at different temperatures for the lowest five reflections (110), (200), (220), (310) and (400) . Our calculated quasi-harmonic results agree reasonably well with the experimental results at temperatures (T) of the order of the Debye temperature ( 0 ). For T » © ^ 9 our calculated anharmonic results are found to be in good agreement with the experimental results.The anomalous terms in the Debye-Waller factor are found not to be negligible for certain reflections even for T ^ ©^ . At temperature T yy Op 9 where the temperature is of the order of the melting temperature (Xm) » "the anomalous terms are found to be important almost for all the f i ve reflections.

The influence of water on the glucose affinity of hexokinase

In the present thesis, the role of hydration during the glucose induced conformational change of hexokinase is investigated. This is accomplished by applying the osmotic stress technique. The osmotic stress technique is founded on varying of the activity of water in a system in order to determine ifs effects. This is accomplished by adding inert solute molecules that are excluded from the system under study. The solute molecules used within the present investigation are Polyethylene glycols (PEGs). PEGs aid in the removal of water from hexokinase by exerting osmotic pressure. The osmotic pressures of the PEG solutions are also measured with both vapour pressure osmometry and secondary osmometry with phospholipids. An interesting discovery is made in that the osmotic pressures of PEG and co-solute solutions are non-additive. This indicates that PEG concentrates co-solutes in solution by making a certain proportion of the water inaccessible. Glucose binding was measured fluorometrically and the glucose equilibrium dissociation constant (GEDC) of hexokinase is measured in solutions containing the different MW PEGs. Changes in the sensitivity of the glucose affinity with osmotic pressure allows the calculation of the change in the numbers of polymer-inaccessible water molecules upon the binding of glucose to hexokinase ~Nw. It was determined the ~Nw decreases with increases in osmotic pressure in the presence of all MW PEGs. ~Nw decreases from values between 45-290 water molecules at low pressure to approximately 15 at high pressure. There is also a molecular weight dependence observed. There are large decreases in ~Nw with osmotic pressure in the presence of PEGs above MW 1000. However, below MW 1500 changes in ~Nw with osmotic pressure are relatively small. These findings are interpreted with respect to two possible mechanisms involving changes in the conformation of hexokinase u~der osmotic pressure and the access of the PEG molecules to water surrounding hexokinase.

Flow regime prediction via froude number calculation in a rock-bedded stream

Mathematical predictions of flow conditions along a steep gradient rock bedded stream are examined. Stream gage discharge data and Manning's Equation are used to calculate alternative velocities, and subsequently Froude Numbers, assuming varying values of velocity coefficient, full depth or depth adjusted for vertical flow separation. Comparison of the results with photos show that Froude Numbers calculated from velocities derived from Manning's Equation, assuming a velocity coefficient of 1.30 and full depth, most accurately predict flow conditions, when supercritical flow is defined as Froude Number values above 0.84. Calculated Froude Number values between 0.8 and 1.1 correlate well with observed transitional flow, defined as the first appearance of small diagonal waves. Transitions from subcritical through transitional to clearly supercritical flow are predictable. Froude Number contour maps reveal a sinuous rise and fall of values reminiscent of pool riffle energy distribution.

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.

Calculation of ditching in the south side of the feeder from Brown’s ditch culvert to lot no.5

Calculation of ditching in the south side of the feeder from Brown’s ditch culvert to lot no.5, n.d.

Statistical Inference for Computable General Equilibrium Models with Application to a Model of the Moroccan Economy

We study the problem of measuring the uncertainty of CGE (or RBC)-type model simulations associated with parameter uncertainty. We describe two approaches for building confidence sets on model endogenous variables. The first one uses a standard Wald-type statistic. The second approach assumes that a confidence set (sampling or Bayesian) is available for the free parameters, from which confidence sets are derived by a projection technique. The latter has two advantages: first, confidence set validity is not affected by model nonlinearities; second, we can easily build simultaneous confidence intervals for an unlimited number of variables. We study conditions under which these confidence sets take the form of intervals and show they can be implemented using standard methods for solving CGE models. We present an application to a CGE model of the Moroccan economy to study the effects of policy-induced increases of transfers from Moroccan expatriates.

Political Influence, Economic Interests and Endogenous Tax Structure in a Computable Equilibrium Framework: with Applicatioon to the United States, 1973 and 1983

A full understanding of public affairs requires the ability to distinguish between the policies that voters would like the government to adopt, and the influence that different voters or group of voters actually exert in the democratic process. We consider the properties of a computable equilibrium model of a competitive political economy in which the economic interests of groups of voters and their effective influence on equilibrium policy outcomes can be explicitly distinguished and computed. The model incorporates an amended version of the GEMTAP tax model, and is calibrated to data for the United States for 1973 and 1983. Emphasis is placed on how the aggregation of GEMTAP households into groups within which economic and political behaviour is assumed homogeneous affects the numerical representation of interests and influence for representative members of each group. Experiments with the model suggest that the changes in both interests and influence are important parts of the story behind the evolution of U.S. tax policy in the decade after 1973.

Confidence Regions for Calibrated Parameters in Computable General Equilibrium Models

We consider the problem of accessing the uncertainty of calibrated parameters in computable general equilibrium (CGE) models through the construction of confidence sets (or intervals) for these parameters. We study two different setups under which this can be done.

Methods to Estimate Dynamic Stochastic General Equilibrium Models

This paper employs the one-sector Real Business Cycle model as a testing ground for four different procedures to estimate Dynamic Stochastic General Equilibrium (DSGE) models. The procedures are: 1 ) Maximum Likelihood, with and without measurement errors and incorporating Bayesian priors, 2) Generalized Method of Moments, 3) Simulated Method of Moments, and 4) Indirect Inference. Monte Carlo analysis indicates that all procedures deliver reasonably good estimates under the null hypothesis. However, there are substantial differences in statistical and computational efficiency in the small samples currently available to estimate DSGE models. GMM and SMM appear to be more robust to misspecification than the alternative procedures. The implications of the stochastic singularity of DSGE models for each estimation method are fully discussed.