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.

Calculation of the magnetic field penetration depth for high-Tc cuprate superconductors based on the Interlayer Pair Tunneling model /

In this work, the magnetic field penetration depth for high-Tc cuprate superconductors is calculated using a recent Interlayer Pair Tunneling (ILPT) model proposed by Chakravarty, Sudb0, Anderson, and Strong [1] to explain high temperature superconductivity. This model involves a "hopping" of Cooper pairs between layers of the unit cell which acts to amplify the pairing mechanism within the planes themselves. Recent work has shown that this model can account reasonably well for the isotope effect and the dependence of Tc on nonmagnetic in-plane impurities [2] , as well as the Knight shift curves [3] and the presence of a magnetic peak in the neutron scattering intensity [4]. In the latter case, Yin et al. emphasize that the pair tunneling must be the dominant pairing mechanism in the high-Tc cuprates in order to capture the features found in experiments. The goal of this work is to determine whether or not the ILPT model can account for the experimental observations of the magnetic field penetration depth in YBa2Cu307_a7. Calculations are performed in the weak and strong coupling limits, and the efi"ects of both small and large strengths of interlayer pair tunneling are investigated. Furthermore, as a follow up to the penetration depth calculations, both the neutron scattering intensity and the Knight shift are calculated within the ILPT formalism. The aim is to determine if the ILPT model can yield results consistent with experiments performed for these properties. The results for all three thermodynamic properties considered are not consistent with the notion that the interlayer pair tunneling must be the dominate pairing mechanism in these high-Tc cuprate superconductors. Instead, it is found that reasonable agreement with experiments is obtained for small strengths of pair tunneling, and that large pair tunneling yields results which do not resemble those of the experiments.

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.

Calculation of the optical conductivity for the Anderson impurity model using the slave boson technique

The optical conductivity of the Anderson impurity mode l has been calculated by emp l oying the slave boson technique and an expansion in powers of l i N, where N is the d egeneracy o f the f electron level . This method has been used to find the effective mass of the conduction electrons for temperatures above and below the Kondo tempera ture. For low temperatures, the mass enhancement is f ound to be large while a t high t emperatures, the mass enhancement is sma ll. The conductivity i s f ound to be Drude like with frequency dependent effective mass and scattering time for low independent effective mass and temperatures and scattering time f requency for high t emperatures. The behavior of both the effective mass and the conductivity is in qualitative agreement with experimental r esul t s .

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

Approximate cost of completing the railway from Port Dalhousie to St. Catharines and an estimate of the cost of the piers at Port Dalhousie signed by William Hamilton Merritt

Approximate cost of completing the railway from Port Dalhousie to St. Catharines and an estimate of the cost of the piers at Port Dalhousie signed by William Hamilton Merritt (5 pages, handwritten), July 8, 1854.

Approximate estimate of the cost of completing the Port Dalhousie Railway to the Grand Central Railway Station at Lock 12

Approximate estimate of the cost of completing the Port Dalhousie Railway to the Grand Central Railway Station at Lock 12. This document is badly torn and burned but most of the text is legible, July 14, 1854.

Approximate estimate of the cost of constructing and completing the Port Dalhousie and Thorold Railway to St. Catharines

Approximate estimate of the cost of constructing and completing the Port Dalhousie and Thorold Railway to St. Catharines signed by S.D. Woodruff (2 pages, handwritten), Jan. 8, 1855.

Approximate estimate of the cost of extending the Port Dalhousie and Thorold Railway from Geneva Street to the Great Western Railway Station at Lock no. 12

Approximate estimate of the cost of extending the Port Dalhousie and Thorold Railway from Geneva Street to the Great Western Railway Station at Lock no. 12 (2 copies) [one appears to be a rough copy] (2 pages, handwritten), Feb. 2, 1855.

Approximate estimate of the cost of macadamizing, grading, bridging and putting in culverts from Hurst’s Bridge above Thorold to Port Robinson

Approximate estimate of the cost of macadamizing, grading, bridging and putting in culverts from Hurst’s Bridge above Thorold to Port Robinson (2 pages, handwritten), n.d.

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.

Chart of approximate quantity of excavation in slides in the deep cut

Chart of approximate quantity of excavation in slides in the deep cut, July 1, 1848.

On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

The kinetics and induced decomposition on the thermal decomposition of hydroperoxides /

This project is focussed on the thermsLl decomposition of t-butyl hydroperoxide and sec-butyl hydroperoxide at 120°C to 160°C in three alcohol solvents. These are methanol, ethajiol and isopropyl alcohol. The aim of the project was to examine the process of induced decomposition. Thermal decomposition of t-hutyl hydroperoxide and sec-butyl hydroperoxide indicate that these reactions have first-order kinetics with activation energies on the order of 20 to 28 K cal/mole, Styrene was used as a free radical trap to inhibit the induced decomposition. The results permitted calculation of how much induced decomposition occurred in its absence. The experimental resvilts indicate that the induced decomposition is important for t-butyl hydroperoxide in alcohol solvents, as shown by both the reaction rate suid product studies. But sec-butyl hydroperoxide results show that the concerted mechanism for the interaction of two sec-butylperoxy radicals occurs in addition to the induced decomposition. Di-sodium E.D,T.A. was added to reduce possible effects of trace transition metal ion .impurities. The result of this experiment were not as expected. The rate of hydroperoxide decomposition was about the same but was zero-order in hydroperoxide concentration.