Monte Carlo study of the XY-model on quasi-periodic lattices

Monte Carlo Simulations were carried out using a nearest neighbour ferromagnetic XYmodel, on both 2-D and 3-D quasi-periodic lattices. In the case of 2-D, both the unfrustrated and frustrated XV-model were studied. For the unfrustrated 2-D XV-model, we have examined the magnetization, specific heat, linear susceptibility, helicity modulus and the derivative of the helicity modulus with respect to inverse temperature. The behaviour of all these quatities point to a Kosterlitz-Thouless transition occuring in temperature range Te == (1.0 -1.05) JlkB and with critical exponents that are consistent with previous results (obtained for crystalline lattices) . However, in the frustrated case, analysis of the spin glass susceptibility and EdwardsAnderson order parameter, in addition to the magnetization, specific heat and linear susceptibility, support a spin glass transition. In the case where the 'thin' rhombus is fully frustrated, a freezing transition occurs at Tf == 0.137 JlkB , which contradicts previous work suggesting the critical dimension of spin glasses to be de > 2 . In the 3-D systems, examination of the magnetization, specific heat and linear susceptibility reveal a conventional second order phase transition. Through a cumulant analysis and finite size scaling, a critical temperature of Te == (2.292 ± 0.003) JI kB and critical exponents of 0:' == 0.03 ± 0.03, f3 == 0.30 ± 0.01 and I == 1.31 ± 0.02 have been obtained.

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.

Polytypism and Silicon carbide : a solid state nuclear magnetic resonance study

A survey of predominantly industrial silicon carbide has been carried out using Magic Angle Spinning nuclear magnetic resonance (MAS nmr); a solid state technique. Three silicon carbide polytypes were studied; 3C, 6H, and 15R. The 13C and 29 Si MAS nmr spectra of the bulk SiC sample was identified on the basis of silicon (carbon) site type in the d iff ere n t pol Y t Y pes • Out to 5.00 A fro mac en t r a lsi 1 i con (0 r carbon) atom four types of sites were characterized using symmetry based calculations. This method of polytype analysis was also considered, in the prelminary stages, for applications with other polytypic material; CdBr 2 , CdI 2 , and PbI 2 " In an attempt to understand the minor components of silicon carbide, such as its surface, some samples were hydrofluoric acid washed and heated to extreme temperatures. Basically, an HF removable species which absorbs at -110 ppm (Si0 2 ) in the 29 Si MAS nmr spectrum is found in silicon carbide after heating. Other unidentified peaks observed at short recycle delays in some 29 Si MAS nmr spectra are considered to be impurities that may be within the lattice. These components comprise less than 5% of the observable silicon. A Tl study was carried out for 29 Si nuclei in a 3C ii polytype sample, using the Driven Equilibrium Single-Pulse Observation of T1 (DESPOT) technique. It appears as though there are a number of nuclei that have the same chemical shift but different T1 relaxation times. The T1 values range from 30 seconds to 11 minutes. Caution has to be kept when interpreting these results because this is the first time that DESPOT has been used for solid samples and it is not likely in full working order. MAS nmr indicates that the 13C and 29 Si ~sotropic chemical shifts of silicon carbide appear to have a reciprocal type of relationship_ Single crystal nmr analysis of a 6H sample is accordance with this finding when only the resultant isotropic shift is considered. However, single crystal nmr also shows that the actual response of the silicon and carbon nuclear environment to the applied magnetic field at various angles is not at all reciprocal. Such results show that much more single crystal nmr work is required to determine the actual behavior of the local magnetic environment of the SiC nuclei.

Solid state NMR chemical shifts as an alternative to diffraction data in the determination of SIC polytypic structures

Silicon carbide, which has many polytypic modifications of a very simple and very symmetric structure, is an excellent model system for exploring, the relationship between chemical shift, long-range dipolar shielding, and crystal structure in network solids. A simple McConnell equation treatment of bond anisotropy effects in a poly type predicts chemical shifts for silicon and carbon sites which agree well with the experiment, provided that contributions from bonds up to 100 A are included in the calculation. The calculated chemical shifts depend on three factors: the layer stacking sequence, electrical centre of gravity, and the spacings between silicon and carbon layers. The assignment of peaks to lattice sites is proved possible for three polytypes (6H, 15R, and 3C). The fact that the calculated chemical shifts are very sensitive to layer spacings provides us a potential way to detennine and refine a crystal structure. In this work, the layer spacings of 6H SiC have been calculated and are within X-ray standard deviations. Under this premise, the layer spacings of 15R have been detennined. 29Si and 13C single crystal nmr studies of 6H SiC polytype indicate that all silicons and carbons are magnetically anisotropic. The relationship between a magnetic shielding tensor component and layer spacings has been derived. The comparisons between experimental and semi-empirical chemical shielding tensor components indicate that the paramagnetic shielding of silicon should be included in the single crystal chemical shift calculation.

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.

The effects of gramicidin on the structure of phospholipid assemblies /

Gramicidin is an antibiotic peptide that can be incorporated into the monolayers of cell membranes. Dimerization through hydrogen bonding between gramicidin monomers in opposing leaflets of the membrane results in the formation of an iontophoretic channel. Surrounding phospholipids, with various associated mechanical properties, have been shown to influence the gating properties of this channel. Conversely, gramicidin incorporation has been shown to affect the structure of spontaneously formed lipid assemblies. Using small-angle x-ray diffraction and model systems composed of phospholipids and gramicidin, the physical effects incurred by gramicidin incorporation were measured. The reverse hexagonal (H^) phase composed of dioleoylphosphatidylethanolamine (DOPE) monolayers decreased in lattice dimension with increasing incorporation of gramicidin. This indicated that gramicidin was adding negative curvature to the monolayers. In this system, gramicidin was measured to have an apparent intrinsic radius of curvature (Rop*™") of -7. 1 A. The addition of up to 4 mol% gramicidin in mixtures with DOPE did not result in the monolayers becoming stiffer, as indicated by unaltered bending moduli for each composition. Dioleoylphosphatidylcholine (DOPC) alone forms the lamellar (LJ phase when hydrated, but undergoes a transition into the H^ phase when mixed with gramicidin. The lattice repeat dimension decreases systematically with increased gramicidin content. Again, this indicated that gramicidin was adding negative curvature to the monolayers. At 12 mol% gramicidin in mixtures with DOPC, the apparent radius of intrinsic curvature of gramicidin (Rop*"^) was measured to be -7.4 A. This mixture formed monolayers that were very resistant to bending under osmotic pressure, with a measured bending modulus of 1 15 kT. The measurements made in this study demonstrate that peptides are able to modulate the spontaneous curvature and other mechanical properties of phospholipid assemblies.

An anharmonic contribution to the Helmholtz free energy O(lambda 6)

The anharmonic contributions of order A6 to the Helmholtz free energy for a crystal in which every atom is on a site of inversion symmetry, have been evaluated The cor~esponding diagrams in the various orders of the perturbation theory have been presented The validity of the expressions given is for high temperatures. Numerical calculations for the diagrams which contribute to the free energy have been worked out for a nearest-n~ighbour central-force model of a facecentered cubic lattice in the high-temperature limit and in the leading term and the Ludwig approximations. The accuracy of the Ludwig approximation in evaluating the Brillouin-zone sums has been investigated. Expansion for all diagrams in the high-temperature limit has been carried out The contribution to the specific heat involves a linear as well as cubic term~ We have applied Lennard-Jones, Morse and Exponential 6 types of potentials. A comparison between the contribution to the free energy of order A6 to that of order A4 has been made.

Elucidation of a mechanism of cell lysis by chorhexidine : a biophysical approach

Chlorhexidine is an effective antiseptic used widely in disinfecting products (hand soap), oral products (mouthwash), and is known to have potential applications in the textile industry. Chlorhexidine has been studied extensively through a biological and biochemical lens, showing evidence that it attacks the semipermeable membrane in bacterial cells. Although extremely lethal to bacterial cells, the present understanding of the exact mode of action of chlorhexidine is incomplete. A biophysical approach has been taken to investigate the potential location of chlorhexidine in the lipid bilayer. Deuterium nuclear magnetic resonance was used to characterize the molecular arrangement of mixed phospholipid/drug formulations. Powder spectra were analyzed using the de-Pake-ing technique, a method capable of extracting both the orientation distribution and the anisotropy distribution functions simultaneously. The results from samples of protonated phospholipids mixed with deuterium-labelled chlorhexidine are compared to those from samples of deuterated phospholipids and protonated chlorhexidine to determine its location in the lipid bilayer. A series of neutron scattering experiments were also conducted to study the biophysical interaction of chlorhexidine with a model phospholipid membrane of DMPC, a common saturated lipid found in bacterial cell membranes. The results found the hexamethylene linker to be located at the depth of the glycerol/phosphate region of the lipid bilayer. As drug concentration was increased in samples, a dramatic decrease in bilayer thickness was observed. Differential scanning calorimetry experiments have revealed a depression of the DMPC bilayer gel-to-lamellar phase transition temperature with an increasing drug concentration. The enthalpy of the transition remained the same for all drug concentrations, indicating a strictly drug/headgroup interaction, thus supporting the proposed location of chlorhexidine. In combination, these results lead to the hypothesis that the drug is folded approximately in half on its hexamethylene linker, with the hydrophobic linker at the depth of the glycerol/phosphate region of the lipid bilayer and the hydrophilic chlorophenyl groups located at the lipid headgroup. This arrangement seems to suggest that the drug molecule acts as a wedge to disrupt the bilayer. In vivo, this should make the cell membrane leaky, which is in agreement with a wide range of bacteriological observations.

The Synthesis and Coordination Chemistry of Two Families of Polydentate Ligands - Exploring Their Potential for the Preparation of Molecule-Based Magnets

The synthesis and studies of two classes of poly dentate ligands are presented as two projects. In project 1, four new carboxamide ligands have been synthesised via the condensation of 2,2',6,6'-tetrachloroformyl-4,4'-bipyridine or 2,6-dichloroformyl pyridine together with heterocyclic amines containing pyridine or pyrazole substituents. The coordination chemistry of these ligands has been investigated and studies have shown that with a Cu(II) salt, two carboxamide ligands LJ and L2 afford large clusters with stoichiometries [Cu8(L1)4Cl16].CHCl3.5H2O.7CH3OH (I) and [Cu9(L2)6Cl6].CH3OH.5H2O.(C2H5)3N (II) respectively. [molecular diagram availabel in pdf]. X-ray diffraction studies of cluster (I) reveal that it has approximate S4 symmetry and is comprised of four ligands and eight copper (II) centers. Here, coordination takes place via amide 0 atoms, and pyrazole nitrogens. This complex is the first reported example of an octanuclear copper cluster with a saddle-shaped structure. The second cluster comprises nine copper ions that are arranged in a cyclic array. Each ligand coordinates three copper centers and each copper ion shares two ligands to connect six ligands with nine copper ions. The amide nitrogens are completely deprotonated and both amide Nand 0 atoms coordinate the metal centres. The cluster has three-fold symmetry. There are six chloride ions, three of which are bridging two neighbouring Cu(II) centres. Magnetic studies of (I) and (II) reveal that both clusters display weak antiferromagnetic interactions between neighbouring Cu(II) centers at low temperature. In the second project, three complexes with stoichiometries [Fe[N302](SCN)2]2 (III), R,R-[Fe[N3O2](SCN)2 (IV) and R,R-]Fe[N3O2](CN)2] (V) were prepared and characterized, where [N302] is a pentadentate macrocycle. Complex (III) was prepared via the metal templated Schiff-base condensation of 2,2',6,6'-tetraacetyl-4,4'-bipyridine together with 3,6-dioxaoctane-I,8-diamine and comprises of a dimeric macro cycle where the two Fe(II) centres are in a pentagonal-bipyramidal environment with the [N302] ligands occupying the equatorial plane and two axial NCS ligands. Complexes (IV) and (V) were prepared via the condensation of 2,6-diacetylpyridine together with a chiral diamine in the presence of FeCh. The synthetic strategy for the preparation of the chiral diamine (4R,5R)-4,5-diphenyl-3,6-dioxa-I,8-octane-diamine was elucidated. The chirality of both macrocycles (IV) and (V) was probed by circular dichroism spectroscopy. The crystal structure of (IV) at 200 K contains two independent molecules in the unit cell, both of which contain a hepta-coordinated Fe(II) and axial NCS ligands. Variable temperature magnetic susceptibility and structural studies are consistent with a high spin Fe(II) complex and show no evidence of any spin crossover behaviour. In contrast, the bis cyanide derivative (V) crystallizes with two independent molecules in the unit cell, both of which have different coordination geometries consistent with different spin states for the two Fe(II) centres. At 250 K, the molecular structure of (V) shows the presence of both 7- and a 6-coordinate Fe(II) complexes in the crystal lattice. As the temperature is lowered, the molecules undergo a structural change and at 100 K the structural data is consistent with a 6- and 5-coordinate Fe(II) complex in the unit cell. Magnetic studies confirm that this complex undergoes a gradual, thermal, spin crossover transition in the solid state. Photomagnetic measurements indicate this is the first chiral Fe (II) sea complex to exhibit a LIESST.

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.

An Abstract Algebraic Theory of L-Fuzzy Relations for Relational Databases

Classical relational databases lack proper ways to manage certain real-world situations including imprecise or uncertain data. Fuzzy databases overcome this limitation by allowing each entry in the table to be a fuzzy set where each element of the corresponding domain is assigned a membership degree from the real interval [0…1]. But this fuzzy mechanism becomes inappropriate in modelling scenarios where data might be incomparable. Therefore, we become interested in further generalization of fuzzy database into L-fuzzy database. In such a database, the characteristic function for a fuzzy set maps to an arbitrary complete Brouwerian lattice L. From the query language perspectives, the language of fuzzy database, FSQL extends the regular Structured Query Language (SQL) by adding fuzzy specific constructions. In addition to that, L-fuzzy query language LFSQL introduces appropriate linguistic operations to define and manipulate inexact data in an L-fuzzy database. This research mainly focuses on defining the semantics of LFSQL. However, it requires an abstract algebraic theory which can be used to prove all the properties of, and operations on, L-fuzzy relations. In our study, we show that the theory of arrow categories forms a suitable framework for that. Therefore, we define the semantics of LFSQL in the abstract notion of an arrow category. In addition, we implement the operations of L-fuzzy relations in Haskell and develop a parser that translates algebraic expressions into our implementation.

L-Fuzzy Structured Query Language

Lattice valued fuzziness is more general than crispness or fuzziness based on the unit interval. In this work, we present a query language for a lattice based fuzzy database. We define a Lattice Fuzzy Structured Query Language (LFSQL) taking its membership values from an arbitrary lattice L. LFSQL can handle, manage and represent crisp values, linear ordered membership degrees and also allows membership degrees from lattices with non-comparable values. This gives richer membership degrees, and hence makes LFSQL more flexible than FSQL or SQL. In order to handle vagueness or imprecise information, every entry into an L-fuzzy database is an L-fuzzy set instead of crisp values. All of this makes LFSQL an ideal query language to handle imprecise data where some factors are non-comparable. After defining the syntax of the language formally, we provide its semantics using L-fuzzy sets and relations. The semantics can be used in future work to investigate concepts such as functional dependencies. Last but not least, we present a parser for LFSQL implemented in Haskell.