Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

Low frequency Raman scattering in amorphous materials: fused quartz, "pyrex" boro-silicate glass and soda-lime silicate glass

Raman scattering in the region 20 to 100 cm -1 for fused quartz, "pyrex" boro-silicate glass, and soft soda-lime silicate glass was investigated. The Raman spectra for the fused quartz and the pyrex glass were obtained at room temperature using the 488 nm exciting line of a Coherent Radiation argon-ion laser at powers up to 550 mW. For the soft soda-lime glass the 514.5 nm exciting line at powers up to 660 mW was used because of a weak fluorescence which masked the Stokes Raman spectrum. In addition it is demonstrated that the low-frequency Raman coupling constant can be described by a model proposed by Martin and Brenig (MB). By fitting the predicted spectra based on the model with a Gaussian, Poisson, and Lorentzian forms of the correlation function, the structural correlation radius (SCR) was determined for each glass. It was found that to achieve the best possible fit· from each of the three correlation functions a value of the SCR between 0.80 and 0.90 nm was required for both quartz and pyrex glass but for the soft soda-lime silicate glass the required value of the SCR. was between 0.50 and 0.60 nm .. Our results support the claim of Malinovsky and Sokolov (1986) that the MB model based on a Poisson correlation function provides a universal fit to the experimental VH (vertical and horizontal polarizations) spectrum for any glass regardless of its chemical composition. The only deficiency of the MB model is its failure to fit the experimental depolarization spectra.

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.

Pressure dependence of superconductivity in amorphous Ni x Zr 100-x alloys

Pressure variations of the superconducting transition temperature Ic of a series of amorphous NixZr 1 OO-x alloys have been studied under quasmydrostatic pressures upto 8 G Pa. For amorphous samples having Ni-concentration less than 40%, i)Tc/dP is positive in sign and it decreases non linearly with increase in I. whereasdTcldP is negative in sign for Ni concentration of 45%. Comparison with the Hall coefficient (I) and the thermoelectric power (2) results for the same amorphous alloys leads to the conclusion that s-d hybridization nature of the d-band (Nil plays a central role in the sign reversal behaviour. Application of pressures greater than 2 G Pa to Ni20ZrgO led to the formation of a new phase, w-Zr. which retains its form after the pressure is released.