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.

Flow injection analysis (FIA) determination of lead by hydride generation in banded coral skeletons /

A flow injection hydride generation direct current plasma atomic emission spectrometric (FI-HG-DCP-AES) method was developed for the determination of lead at ng.ml-l level. Potassium ferricyanide (K3Fe(CN)6) was used along with sodium tetrahydroborate(III) (NaBH4) to produce plumbane (PbH4) in an acid medium. The design of a gas-liquid separator (hydride generator) was tested and the parameters of the flow injection system were optimized to achieve a good detection limit and sample throughput. The technique developed gave a detection limit of 0.7 ng.ml-l(3ob). The precision at 20 ng.ml"* level was 1.6 % RSD with 1 1 measurements (n=l 1). Volume of sample loop was 500 |J.l. A sample throughput of 120 h"^ was achieved. The transition elements, Fe(II), FeOH), Cd(n), Co(II), Mn(n), Ni(II) and Zn(n) do not interfere in this method but 1 mg,l'l Cu(II) will suppress 50 % of the signal from a sample containing 20 ng.ml'l Pb. This method was successfully applied to determine lead in a calcium carbonate (CaC03) matrix of banded coral skeletons from Si-Chang Island in Thailand.

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.

An application of cell-cluster theory to a rare gas crystal /|n[by] K. Westera. - 260 St. Catharines [Ont.] : Dept. of Physics, Brock University,

The Lennard-Jones Devonshire 1 (LJD) single particle theory for liquids is extended and applied to the anharmonic solid in a high temperature limit. The exact free energy for the crystal is expressed as a convergent series of terms involving larger and larger sets of contiguous particles called cell-clusters. The motions of all the particles within cell-clusters are correlated to each other and lead to non-trivial integrals of orders 3, 6, 9, ... 3N. For the first time the six dimensional integral has been calculated to high accuracy using a Lennard-Jones (6-12) pair interaction between nearest neighbours only for the f.c.c. lattice. The thermodynamic properties predicted by this model agree well with experimental results for solid Xenon.

A treatment of atomic mean square displacement in higher order perturbation theory

A general derivation of the anharmonic coefficients for a periodic lattice invoking the special case of the central force interaction is presented. All of the contributions to mean square displacement (MSD) to order 14 perturbation theory are enumerated. A direct correspondance is found between the high temperature limit MSD and high temperature limit free energy contributions up to and including 0(14). This correspondance follows from the detailed derivation of some of the contributions to MSD. Numerical results are obtained for all the MSD contributions to 0(14) using the Lennard-Jones potential for the lattice constants and temperatures for which the Monte Carlo results were calculated by Heiser, Shukla and Cowley. The Peierls approximation is also employed in order to simplify the numerical evaluation of the MSD contributions. The numerical results indicate the convergence of the perturbation expansion up to 75% of the melting temperature of the solid (TM) for the exact calculation; however, a better agreement with the Monte Carlo results is not obtained when the total of all 14 contributions is added to the 12 perturbation theory results. Using Peierls approximation the expansion converges up to 45% of TM• The MSD contributions arising in the Green's function method of Shukla and Hubschle are derived and enumerated up to and including 0(18). The total MSD from these selected contributions is in excellent agreement with their results at all temperatures. Theoretical values of the recoilless fraction for krypton are calculated from the MSD contributions for both the Lennard-Jones and Aziz potentials. The agreement with experimental values is quite good.

Investigation into the determination of arsenic by direct current plasma atomic emission spectrometer

Improvements have been made on the currently available hydride generator system manufactured by SpectraMetrics Incorporated, because the system was found to be unsatisfactory with respect to the following: 1. the drying agent, anhydrous calcium chloride, 2. the special sample tube, 3. the direction of argon flow through the Buchner funnel when it came to dealing with real sample, that is, with reference only to aqueous extracts of soil samples. Changes that were made on the system included the replacement of anhydrous calcium chloride with anhydrous calcium sulphate and the replacement of the special sample tube with a modified one made from silica. Re-directing the flow of argon through the top of the Buchner funnel appeared to make the system compatible with aqueous extracts of soil samples. The interferences from 1000 ~g/mL of nickel(II) , cobalt(II), iron(III), copper(II) have been eliminated with the aid of 1.4 M hydrochloric acid and 1% (weight/volume) L-cystine. Greater than 90% recovery of 0.3 ~g/mL arsenic signal was achieved in each case. Furthermore, 103% of arsenic signal was accomplished in the presence of 1000 ~g/mL cadmium with 5 M Hel. tVhen each of the interferents was present in solution at 1000 ppm, a recovery of 85% was achieved by using 5 M hydrochloric acid and 3% (weight/volume) L-cystine. Without L-cystine and when 1.4 M hydrochloric acid was used, the recoveries were 0% (Ni), 0% (Co), 88% (Fe), 15% (Cu), 18% (Cd). Similarly, a solution containing 1000 ppm of each interferent gave a zero percent recovery of arsenic. The reduction of trivalent and pentavalent arsenic at a pH less than one has also been investigated and shown to be quantitative if peak areas are measured. The reproducibility determination of a 0.3 Vg/mL standard arsenic solution by hydride generation shows a relative standard deviation of 3.4%. The detection limits with and without Porapak Q have been found to be 0.6 ng/mL and 1.0 ng/mL, respectively.

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.

Investigations into the extraction and the determination of polycyclic aromatic hydrocarbons (PAHs) from water by solid-phase extraction and capillary gas chromatography/ mass spectrometry

Factors affecting the detennination of PAHs by capillary GC/MS were studied. The effect of the initial column temperature and the injection solvent on the peak areas and heights of sixteen PAHs, considered as priority pollutants, USillg crosslinked methyl silicone (DB!) and 5% diphenyl, 94% dimethyl, 1% vinyl polysiloxane (DBS) columns was examined. The possibility of using high boiling point alcohols especially butanol, pentanol, cyclopentanol, and hexanol as injection solvents was investigated. Studies were carried out to optimize the initial column temperature for each of the alcohols. It was found that the optimum initial column temperature is dependent on the solvent employed. The peak areas and heights of the PAHs are enhanced when the initial column temperature is 10-20 c above the boiling point of the solvent using DB5 column, and the same or 10 C above the boiling point of the solvent using DB1 column. Comparing the peak signals of the PAHs using the alcohols, p-xylene, n-octane, and nonane as injection solvents, hexanol gave the greatest peak areas and heights of the PAHs particularly the late-eluted peaks. The detection limits were at low pg levels, ranging from 6.0 pg for fluorene t9 83.6 pg for benzo(a)pyrene. The effect of the initial column temperature on the peak shape and the separation efficiency of the PARs was also studied using DB1 and DB5 columns. Fronting or splitting of the peaks was obseIVed at very low initial column temperature. When high initial column temperature was used, tailing of the peaks appeared. Great difference between DB! and.DB5 columns in the range of the initial column temperature in which symmetrical.peaks of PAHs can be obtained is observed. Wider ranges were shown using DB5 column. Resolution of the closely-eluted PAHs was also affected by the initial column temperature depending on the stationary phase employed. In the case of DB5, only the earlyeluted PAHs were affected; whereas, with DB1, all PAHs were affected. An analytical procedure utilizing solid phase extraction with bonded phase silica (C8) cartridges combined with GC/MS was developed to analyze PAHs in water as an alternative method to those based on the extraction with organic solvent. This simple procedure involved passing a 50 ml of spiked water sample through C8 bonded phase silica cartridges at 10 ml/min, dried by passing a gentle flow of nitrogen at 20 ml/min for 30 sec, and eluting the trapped PAHs with 500 Jll of p-xylene at 0.3 ml/min. The recoveries of PAHs were greater than 80%, with less than 10% relative standard deviations of nine determinations. No major contaminants were present that could interfere with the recognition of PAHs. It was also found that these bonded phase silica cartridges can be re-used for the extraction of PAHs from water.

Investigations into the determinations of hydride-forming elements

Analytical methods for the determination of trace amounts of germanium, tin and arsenic were established using hydride generation coupled with direct current plasma atomic emission spectrometry. A continuous gas flowing batch system for the hydride generation was investigated and was applied to the determination of germanium(Ge), tin(Sn), antimony(Sb) and lead(Pb) (Preliminary results suggest that it is also applicable to arsenic)As) ). With this system, the reproducibility of signals was improved and the determination was speeded up, compared with the conventional batch type hydride generation system. Each determination was complete within one minute. Interferences from a number of transition metal ions, especially from Pd(II), Pt(IV), Ni(II), Cu(II), Co(II), and Fe(II, III), have proven to be very serious under normal conditions, in the determination of germanium, tin, and arsenic. These interference effects were eliminated or significantly reduced in the presence of L-cystine or L-cysteine. Thus, a 10-1000 fold excess of Ni(II), Cu(II), Co(II), Fe(II), Pt(IV), Pd(II), etc. can be tolerated without interference, In the presence of L-cystine or L-cysteine, compared with absence of interference reducing agent. The methods for the determination of Ge, Sn, and As were examined by the analyses of standard reference materials. Interference effects from the sample matrix, for example, in transition metal-rich samples, copper, iron and steel, were eliminated by L-cystine (for As and Sn) and by LI cysteine (for Ge). The analysis of a number of standard reference materials gave excellent results of As and Sn contents in agreement with the certified values, showing there was no systematic interference. The detection limits for both germanium and tin were 20 pg ml- I . Preliminary studies were carried out for the determination of antimony and lead. Antimony was found to react with NaBH4, remaInIng from the previous determinations, giving an analytical signal. A reversed injection manner, i.e., injection of the NaBH4 solution prior to the analyte solution was used to avoid uncertainty caused by residual NaBH4 present and to ensure that an excess of NaB H4 was available. A solution of 0.4% L-cysteine was found to reduce the interference from selected transition metal ions, Co(II), Cu(II), Ni(II) and Pt(IV). Hydrochloric acid - hydrogen peroxide, nitric acid - ammonium persulphate, and potassium dichromate malic acid reaction systems for lead hydride generation were compared. The potassium dichromate - malic acid reaction medium proved to be the best with respect to reproducibility and minimal interference. Cu(II), Ni(II), and Fe(II) caused strong interference In lead determinations, which was not reduced by L-cysteine or Lcystine. Sodium citrate, ascorbic acid, dithizone, thiosemicarbazide and penicillamine reduced interferences to some extent. Further interference reduction studies were carried out uSIng a number of amino acids, glycine, alanine, valine, leucine and histidine, as possible interference reducing agents in the determination of germanium. From glycine, alanine, valine to leucine, the interference reduction effect in germanium determinations decreased. Histidine II was found to be very promising In the reduction of interference. In fact, histidine proved more efficient than L-cystine and L-cysteine In the reduction of interference from Ni(II) in the determination of germanium. Signal enhancement by easily ionized elements (EIEs), usually regarded as an interference effect in analysis by DCP-AES, was studied and successfully applied to advantage in improving the sensitivity and detection limit in the determination of As, Ge, Sn, Sb, and Pb. The effect of alkali and alkaline-earth elements on the determination of the above five hydride forming elements was studied. With the appropriate EIE, a signal enhancement of 40-115% was achieved. Linear calibration and good reproducibility were also obtained in the presence of EIEs. III

Modified statistical methods to calculate rare gas interaction potentials

New density functionals representing the exchange and correlation energies (per electron) are employed, based on the electron gas model, to calculate interaction potentials of noble gas systems X2 and XY, where X (and Y) are He,Ne,Ar and Kr, and of hydrogen atomrare gas systems H-X. The exchange energy density functional is that recommended by Handler and the correlation energy density functional is a rational function involving two parameters which were optimized to reproduce the correlation energy of He atom. Application of the two parameter function to other rare gas atoms shows that it is "universal"; i. e. ,accurate for the systems considered. The potentials obtained in this work compare well with recent experimental results and are a significant improvement over those from competing statistical modelS.

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 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.

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.

Response rainbow trout (Salmon gairdneri) blood gas transport system to temperature, oxygen availability and photoperiod

Hematological status in rainbow trout, Salmo gairdneri, was examined in relation to eight combinations of three environmental fa ctors; temperature (5°, 20°C), oxygen availability «35%, >70% saturation) and photoperiod (16L:8D, 8L:16D) and evaluated by 3-factor analysis of variance. Hemog l obin and hematocrit , indicators of oxygenc arrying capacity increased significantly at the higher temperature, following exposure to hypoxia and in relation to reduced light period. Significant variations in mean corpuscular hemoglobin concentration were not detected. The effects of temperature and oxygen availability were more pronounced than that of photoperiod which was generally masked. Although oxygen availability and photoperiod did not interact with temperature, the interaction of the former fac tors was significant. Elec trophoresis revealed twelve hemoglobin isomorphs. Relative concentration changes were found in re lation to the factors c onsidered with temperature>hypoxia>photoperiod. Howeve r , in terms of absolute concentration, effects were hypoxia>temperature>photoperiod. Photoperiod effects were again masked by temperature and (or) hypoxia. Red cell +2 l eve ls of [CI ] and [Mg ], critical elements in the hemoglobin-oxygen affinity regulating system, were also significantly altered. Red cell CI +2 was influenced only by temperature ; Mg by temper ature and oxygen. No photoperiod influence on either ions was observed. Under nominal 'summer' conditions, these changes point to the likelihood of increases in oxygen-c arrying c apac ity coupled with low Hb-02 affinity adjustments which would be expected to increase oxygen delivery rates to their more rapidly metabolising tissues.