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.

Phonon spectra and thermal properties of some fcc metals using embedded-atom potentials /

By employing the embedded-atom potentials of Mei et ai.[l], we have calculated the dynamical matrices and phonon dispersion curves for six fee metals (Cu,Ag,Au,Ni,Pd and Pt). We have also investigated, within the quasiharmonic approximation, some other thermal properties of these metals which depend on the phonon density of states, such as the temperature dependence of lattice constant, coefficient of linear thermal expansion, isothermal and adiabatic bulk moduli, heat capacities at constant volume and constant pressure, Griineisen parameter and Debye temperature. The computed results are compared with the experimental findings wherever possible. The comparison shows a generally good agreement between the theoretical values and experimental data for all properties except the discrepancies of phonon frequencies and Debye temperature for Pd, Pt and Au. Further, we modify the parameters of this model for Pd and Pt and obtain the phonon dispersion curves which is in good agreement with experimental data.

On the anharmonic, multi-phonon, Debye-Waller contributions to the phonon-limited resistivity of metals : applications to Na and K

The anharmonic, multi-phonon (MP), and Oebye-Waller factor (OW) contributions to the phonon limited resistivity (;0) of metals derived by Shukla and Muller (1979) by the doubletime temperature dependent Green function method have been numerically evaluated for Na and K in the high temperature limit. The anharmonic contributions arise from the cubic and quartic shift of phonons (CS, QS), and phonon width (W) and the interference term (1). The QS, MP and OW contributions to I' are also derived by the matrix element method and the results are in agreement with those of Shukla and Muller (1979). In the high temperature limit, the contributions to;O from each of the above mentioned terms are of the type BT2 For numerical calculations suitable expressions are derived for the anharmonic contributions to ~ in terms of the third and fourth rank tensors obtained by the Ewald procedure. The numerical calculation of the contributions to;O from the OW, MP term and the QS have been done exactly and from the CS, Wand I terms only approximately in the partial and total Einstein approximations (PEA, TEA), using a first principle approach (Shukla and Taylor (1976)). The results obtained indicate that there is a strong pairwise cancellation between the: OW and MP terms, the QS and CS and the Wand I terms. The sum total of these contributions to;O for Na and K amounts to 4 to 11% and 2 to 7%, respectively, in the PEA while in the TEA they amount to 3 to 7% and 1 to 4%, respectively, in the temperature range.

Phonon dispersion curves and atomic mean square displacement for several fcc and bcc materials

The atomic mean square displacement (MSD) and the phonon dispersion curves (PDC's) of a number of face-centred cubic (fcc) and body-centred cubic (bcc) materials have been calclllated from the quasiharmonic (QH) theory, the lowest order (A2 ) perturbation theory (PT) and a recently proposed Green's function (GF) method by Shukla and Hiibschle. The latter method includes certain anharmonic effects to all orders of anharmonicity. In order to determine the effect of the range of the interatomic interaction upon the anharmonic contributions to the MSD we have carried out our calculations for a Lennard-Jones (L-J) solid in the nearest-neighbour (NN) and next-nearest neighbour (NNN) approximations. These results can be presented in dimensionless units but if the NN and NNN results are to be compared with each other they must be converted to that of a real solid. When this is done for Xe, the QH MSD for the NN and NNN approximations are found to differ from each other by about 2%. For the A2 and GF results this difference amounts to 8% and 7% respectively. For the NN case we have also compared our PT results, which have been calculated exactly, with PT results calculated using a frequency-shift approximation. We conclude that this frequency-shift approximation is a poor approximation. We have calculated the MSD of five alkali metals, five bcc transition metals and seven fcc transition metals. The model potentials we have used include the Morse, modified Morse, and Rydberg potentials. In general the results obtained from the Green's function method are in the best agreement with experiment. However, this improvement is mostly qualitative and the values of MSD calculated from the Green's function method are not in much better agreement with the experimental data than those calculated from the QH theory. We have calculated the phonon dispersion curves (PDC's) of Na and Cu, using the 4 parameter modified Morse potential. In the case of Na, our results for the PDC's are in poor agreement with experiment. In the case of eu, the agreement between the tlleory and experiment is much better and in addition the results for the PDC's calclliated from the GF method are in better agreement with experiment that those obtained from the QH theory.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.

Far infrared reflectance and optical properties of organic superconductor (TMTSF)2ClO4 /

Polarized reflectance measurements of the quasi I-D charge-transfer salt (TMTSFh CI04 were carried out using a Martin-Puplett-type polarizing interferometer and a 3He refrigerator cryostat, at several temperatures between 0.45 K and 26 K, in the far infrared, in the 10 to 70 cm- 1 frequency range. Bis-tetramethyl-tetraselena-fulvalene perchlorate crystals, grown electrochemically and supplied by K. Behnia, of dimensions 2 to 4 by 0.4 by 0.2 mm, were assembled on a flat surface to form a mosaic of 1.5 by 3 mm. The needle shaped crystals were positioned parallel to each other along their long axis, which is the stacking direction of the planar TMTSF cations, exposing the ab plane face (parallel to which the sheets of CI04 anions are positioned). Reflectance measurements were performed with radiation polarized along the stacking direction in the sample. Measurements were carried out following either a fast (15-20 K per minute) or slow (0.1 K per minute) cooling of the sample. Slow cooling permits the anions to order near 24 K, and the sample is expected to be superconducting below 1.2 K, while fast cooling yields an insulating state at low temperatures. Upon the slow cooling the reflectance shows dependence with temperature and exhibits the 28 cm- 1 feature reported previously [1]. Thermoreflectance for both the 'slow' and 'fast' cooling of the sample calculated relative to the 26 K reflectance data indicates that the reflectance is temperature dependent, for the slow cooling case only. A low frequency edge in the absolute reflectance is assigned an electronic origin given its strong temperature dependence in the relaxed state. We attribute the peak in the absolute reflectance near 30 cm-1 to a phonon coupled to the electronic background. Both the low frequency edge and the 30 cm-1 feature are noted te shift towards higher frequcncy, upon cntering the superconducting state, by an amount of the order of the expected superconducting energy gap. Kramers-Kronig analysis was carried out to determine the optical conductivity for the slowly cooled sample from the measured reflectance. In order to do so the low frequency data was extrapolated to zero frequency using a Hagen-Ru bens behaviour, and the high frequency data was extended with the data of Cao et al. [2], and Kikuchi et al. [3]. The real part of the optical conductivity exhibits an asymmetric peak at 35 cm-1, and its background at lower frequencies seems to be losing spectral weight with lowering of the temperature, leading us to presume that a narrow peak is forming at even lower frequencies.

The infared optical properties of Sr2RuO4 and SmTiO3 including an object-oriented resistivity interface /

The perovskite crystal structure is host to many different materials from insulating to superconducting providing a diverse range of intrinsic character and complexity. A better fundamental description of these materials in terms of their electronic, optical and magnetic properties undoubtedly precedes an effective realization of their application potential. SmTiOa, a distorted perovskite has a strongly localized electronic structure and undergoes an antiferromagnetic transition at 50 K in its nominally stoichiometric form. Sr2Ru04 is a layered perovskite superconductor (ie. Tc % 1 K) bearing the same structure as the high-tem|>erature superconductor La2_xSrrCu04. Polarized reflectance measurements were carried out on both of these materials revealing several interesting features in the far-infrared range of the spectrum. In the case of SmTiOa, although insulating, evidence indicates the presence of a finite background optical conductivity. As the temperature is lowered through the ordering temperature a resonance feature appears to narrow and strengthen near 120 cm~^ A nearby phonon mode appears to also couple to this magnetic transition as revealed by a growing asymmetry in the optica] conductivity. Experiments on a doped sample with a greater itinerant character and lower Neel temperature = 40 K also indicate the presence of this strongly temperature dependent mode even at twice the ordering temperature. Although the mode appears to be sensitive to the magnetic transition it is unclear whether a magnon assignment is appropriate. At very least, evidence suggests an interesting interaction between magnetic and electronic excitations. Although Sr2Ru04 is highly anisotropic it is metallic in three-dimensions at low temperatures and reveals its coherent transport in an inter-plane Drude-like component to the highest temperatures measured (ie. 90 K). An extended Drude analysis is used to probe the frequency dependent scattering character revealing a peak in both the mass enhancement and scattering rate near 80 cm~* and 100 cm~* respectively. All of these experimental observations appear relatively consistent with a Fermi-liquid picture of charge transport. To supplement the optical measurements a resistivity station was set up with an event driven object oriented user interface. The program controls a Keithley Current Source, HP Nano-Voltmeter and Switching Unit as well as a LakeShore Temperature Controller in order to obtain a plot of the Resistivity as a function of temperature. The system allows for resistivity measurements ranging from 4 K to 290 K using an external probe or between 0.4 K to 295 K using a Helium - 3 Cryostat. Several materials of known resistivity have confirmed the system to be robust and capable of measuring metallic samples distinguishing features of several fiQ-cm.

Far infrared optical properties of V and Cr doped Sb2Te3 /

The far infrared reflectance of Sb2Te3 , Sbi.97Vo.o3Te3 and Sbi.94Cr .o6Te3 was measured near normal incidence at different temperatures (between 45K and 300K). The direct current resistivities of the above samples were also measured between the temperatures of 4K and 300K. Also Kramers Kronig (KK) analyses were performed on the reflectance spectra to obtain the optical conductivities. In the doped samples, it was observed that a phonon at 62cm-1 softens to about 55cm-1 on decreasing the temperature from 295K to 45K. Also, it was observed that the plasma frequency of the doped samples is independent of doping. The scattering rate for the vanadium doped sample was seen to be greater than that for the chromium doped sample despite the fact that vanadium impurity density is less than that of chromium. The Drude-Lorentz model fits to the KK optical conductivity show that the samples used in this work are conventional metals. Definitive measurements of the temperature dependence of the scattering rate across the ferromagnetic transition await equipment changes allowing measurements at low temperature using the mercury cadmium telluride (MCT) detector.

Crystal growth, Raman scattering and optical properties of the superconductor Cd2Re2O7(16O, 18O)/

Using the Physical Vapor Transport method, single crystals of Cd2Re207 have been grown, and crystals of dimensions up to 8x6x2 mm have been achieved. X-ray diffraction from a single crystal of Cd2Re207 has showed the crystal growth in the (111) plane. Powder X-ray diffraction measurements were performed on ^^O and ^^O samples, however no difference was observed. Assigning the space group Fd3m to Cd2Re207 at room temperature and using structure factor analysis, the powder X-ray diffraction pattern of the sample was explained through systematic reflection absences. The temperatiure dependence of the resistivity measurement of ^^O has revealed two structural phase transitions at 120 and 200 K, and the superconducting transition at 1.0 K. Using Factor Group Analysis on three different structiures of Cd2Re207, the number of IR and Raman active phonon modes close to the Brillouin zone centre have been determined and the results have been compared to the temperature-dependence of the Raman shifts of ^^O and ^*0 samples. After scaling (via removing Bose-Einstein and Rayleigh scattering factors from the scattered light) all spectra, each spectrum was fitted with a number of Lorentzian peaks. The temperature-dependence of the FWHM and Raman shift of mode Eg, shows the effects of the two structurjil phase transitions above Tc. The absolute reflectance of Cd2Re207 - '^O single crystals in the far-infrared spectral region (7-700 cm~^) has been measured in the superconducting state (0.5 K), right above the superconducting state (1.5 K), and in the normal state (4.2 K). Thermal reflectance of the sample at 0.5 K and 1.5 K indicates a strong absorption feature close to 10 cm~^ in the superconducting state with a reference temperature of 4.2 K. By means of Kramers-Kronig analysis, the absolute reflectance was used to calculate the optical conductivity and dielectric function. The real part of optical conductivity shows five distinct active phonon modes at 44, 200, 300, 375, and 575 cm~' at all temperatures including a Drude-like behavior at low frequencies. The imaginary part of the calculated dielectric function indicates a mode softening of the mode 44 cm~' below Tc.

An analysis of decoupling procedures in generating thermal Green's functions of O([lambda]â ´) by the Zubarev equation of motion method

The Zubarev equation of motion method has been applied to an anharmonic crystal of O( ,,4). All possible decoupling schemes have been interpreted in order to determine finite temperature expressions for the one phonon Green's function (and self energy) to 0()\4) for a crystal in which every atom is on a site of inversion symmetry. In order to provide a check of these results, the Helmholtz free energy expressions derived from the self energy expressions, have been shown to agree in the high temperature limit with the results obtained from the diagrammatic method. Expressions for the correlation functions that are related to the mean square displacement have been derived to 0(1\4) in the high temperature limit.

Magnetic and high pressure studies in the YPd5B3C3 system

The macroscopic properties of the superconducting phase in the multiphase compound YPd5B3 C.3 have been investigated. The onset of superconductivity was observed at 22.6 K, zero resistance at 21.2 K, the lower critical field Hel at 5 K was determined to be Hel (5) rv 310 Gauss and the compound was found to be an extreme type-II superconductor with the upper critical field in excess of 55000 Gauss at 15 K. From the upper and lower critical field values obtained, several important parameters of the superconducting state were determined at T = 15 K. The Ginzburg-Landau paramater was determined to be ~ > 9 corresponding to a coherence length ~ rv 80A and magnetic penetration depth of 800A. In addition measurements of the superconducting transition temperature Te(P) under purely hydrostatically applied pressure have been carried out. Te(P) of YPd5B3 C.3 decreases linearly with dTe/dP rv -8.814 X 10-5 Jphonon spectrum and the density of states near the Fermi level.

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.

Pressure dependence of spin fluctuations in metallic glasses Zrâ â Niâ â ,Zrâ â (Nix/Feâ â â -x)â â and /Zrx/feâ â â -x

The superconducting transition temperature Tc of metallic glasses ZrxFelOO-x (x=80, 75), Zr75(NixFelOO-x)25 (x=75, 50, 25), and CU2SZr75 were measured under quasi-hydrostatic pressure up to 8 OPa (80kbar). The volume (pressure) dependence of the electron-phonon coupling parameters Aep for CU25Zr75 was calculated using the McMillan equatio11. Using this volume dependence of Aep and the modified McMillan equation which incorporates spin-fluctuations, the volume dependence of the spin fluctuation parameter, Asf, was determined in Zr75Ni25, ZrxFelOO-x , a11d Zr75(NixFelOO-x)25. It was found that with increasing pressure, spinfluctuations are suppressed at a faster rate in ZrxFe lOO-x and Zr75(NixFelOO-x)25, as Fe concentration is increased. The rate of suppression of spin-fluctuations with pressure was also found to be higher in Fe-Zr glasses than in Ni-Zr glasses of similar composition.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Pressure dependence of superconducting properties of MgB2 /

We prepared samples of MgB2 and ran sets of experiments aimed for investigation of superconducting properties under pressure. We found the value of pressure derivative of the transition temperature -1.2 ± 0.05 K/GPa. Then, using McMillan formula, we found that the main contribution to the change of the transition temperature under the pressure is due to the change in phonon frequencies. Griineisen parameter was calculated to be 7g = 2.4. Our results suggest that MgB2 is a conventional superconductor.