Comparison of model potentials for molecular dynamics simulation of crystalline silica

Simulationen von SiO2 mit dem von van Beest, Kramer und vanSanten (BKS) entwickelten Paarpotenzial erzeugen vielezufriedenstellende Ergebnisse, aber auch charakteristischeSchwachstellen. In dieser Arbeit wird das BKS-Potenzial mitzwei kürzlich vorgeschlagenen Potenzialen verglichen, dieeffektiv Mehrteilchen-Wechselwirkungen beinhalten. Der ersteAnsatz erlaubt dazu fluktuierende Ladungen, der zweiteinduzierbare Polarisierungen auf den Sauerstoffatomen. Die untersuchten Schwachstellen des BKS Potenzialsbeinhalten das Verhältnis der zwei Gitterkonstanten a und cim Quarzübergang, das von BKS falsch beschrieben wird.Cristobalit und Tridymit erscheinen instabil mit BKS.Weiterhin zeigt die BKS-Zustandsdichte charakteristischeAbweichungen von der wahren Zustandsdichte. DerÜbergangsdruck für den Stishovit I-II Übergang wird deutlichüberschätzt. Das Fluktuierende-Ladungs-Modell verbesserteinige der genannten Punkte, reproduziert aber viele andereEigenschaften schlechter als BKS. DasFluktierende-Dipol-Modell dagegen behebt alle genanntenArtefakte. Zusätzlich wird der druckinduzierte Phasenübergang imalpha-Quarz untersucht. Alle Potentiale finden die selbeStruktur für Quarz II. Bei anschliessender Dekompressionerzeugt BKS eine weitere Phase, während die beiden anderenPotentiale wieder zum alpha-Quarz zurückkehren. Weiterhinwerden zwei Methoden entwickelt, um die piezoelektrischenKonstanten bei konstantem Druck zu bestimmen. Die Ergebnissegeben Hinweise auf eine möglicherweisenicht-elektrostatische Natur der Polarisierungen imFluktuierende-Dipole-Modell. Mit dieser Interpretation scheint das Fluktuierende-DipolPotential alle verfügbaren experimentellen Daten am bestenvon allen drei untersuchten Ansätzen zu reproduzieren.

Effet de la symétrie sur la supraconductivité de LaRhSi3 et la frustration magnétique du SrRE2O4 (RE = Dy ou Ho)

Dans la première partie, nous présentons les résultats de l'étude du supraconducteur sans inversion de symétrie LaRhSi3 par spectroscopie muonique. En champ nul, nous n'avons pas détecté de champ interne. Ceci indique que la fonction d'onde de l'état supraconducteur n'est pas dominée par l'état triplet. Les mesures en champ transverse de 35G présentent une transition en accord avec la transition de phase attendue sous le champ critique Hc1. Nous avons répété ces mesures pour un champ entre Hc1 et Hc2, 150G. Le spectre obtenu pour ces mesures conserve l'asymétrie et relaxe rapidement à basse température tel que prédit pour un supraconducteur dans la phase d'Abrikosov. Néanmoins, les relaxations produites par ce balayage en température présentent une transition à près de 2 fois la température critique attendue. Dans la deuxième partie de ce mémoire, nous donnons l'interprétation des résultats de la diffraction neutronique inélastique par l'étude des champs électriques cristallins. Ces mesures ont été effectuées sur des aimants frustrés SrHo2O4 et SrDy2O4 sous la forme de poudre. L'étude des niveaux produits par les champs cristallins par la méthode des opérateurs de Stevens indique une perte du moment cinétique dans les deux matériaux. Pour le SrDy2O4, l'état fondamental serait constitué de quatre états dégénérés quasi accidentellement qui portent un moment magnétique total non-nul. Toute fois, nos mesures de susceptibilité magnétique ne montrent aucun ordre au-dessus de 50mK. Pour le SrHo2O4, le fondamental est formé d'une paire accidentelle. Nous obtenons un moment magnétique de 6.94(8)$\mu_B$ ce qui s'accorde avec les données expérimentales.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Ionization and Structural Changes of the DMPG Vesicle along Its Anomalous Gel-Fluid Phase Transition: A Study with Different Lipid Concentrations

Dispersions of saturated anionic phospholipid dimyristoyl phosphatidylglycerol (DMPG) have been extensively studied regarding their peculiar thermostructural behavior. At low ionic strength, the gel-fluid transition is spread along nearly 17 degrees C, displaying several thermal events in the calorimetric profile that is quite different from the single sharp peak around 23 degrees C found for higher ionic strength DMPG dispersions. To investigate the role of charge in the bilayer transition, we carefully examine the temperature dependence of the electrical conductivity of DMPG dispersions at different concentrations, correlating the data with the corresponding differential scanning calorimetry (DSC) traces. Electrical conductivity together with electrophoretic mobility measurements allowed the calculation of the dependence of the degree of ionization of DMPG vesicles on lipid concentration and temperature. It was shown that there is a decrease in vesicle charge as the lipid concentration increases, which is probably correlated with the increase in the concentration of bulk Na(+). Apart from the known increase in the electrical conductivity along the DMPG temperature transition region, a sharp rise was observed at the bilayer pretransition for all lipid concentrations studied, possibly indicating that the beginning of the chain melting process is associated with an increase in bilayer ionization. It is confirmed here that the gel-fluid transition of DMPG at low ionic strength is accompanied by a huge increase in the dispersion viscosity. However, it is shown that this measured macroviscosity is distinct from the local viscosity felt by either charged ions or DMPG charged aggregates in measurements of electrical conductivity or electrophoretic mobility, Data presented here give support to the idea that DMPG vesicles, at low ionic strength, get more ionized along the temperature transition region and could be perforated and/or deformed vesicle structures.

Synthesis, spectroscopic characterization and X-ray structures of five-coordinate diorganotin(IV) complexes containing 5-hydroxypyrazoline derivatives as ligands

Four new diorganotin(IV) complexes have been prepared from R(2)SnCl(2) (R = Me, Ph) with the ligands 5-hydroxy-3-metyl-5-phenyl-1-(S-benzildithiocarbazate)-pyrazoline (H(2)L(1)) and 5-hydroxy-3-methyl-5-phenyl-1-(2-thiophenecarboxylic)-pyrazoline (H(2)L(2)). The complexes were characterized by elemental analysis, IR. (1)H (13)C, (119)Sn NMR and Mossbauer spectroscopes The complexes [Me(2)SnL(1)], [Ph(2)SnL(1)] and [Me(2)SnL(2)] were also studied by single crystal X-ray diffraction and the results showed that the Sn(IV) central atom of the complexes adopts a distorted trigonal bipyramidal (TBP) geometry with the N atom of the ONX-tridentate (X = O and S) ligand and two organic groups occupying equatorial sites. The C-Sn-C angles for [Me(2)Sn(L(1))] and [Ph(2)Sn(L(1))] were calculated using a correlation between (119)Sn Mossbauer and X-ray crystallographic data based on the point-charge model Theoretical calculations were performed with the B3LYP density functional employing 3-21G(*) and DZVP all electron basis sets showing good agreement with experimental findings General and Sn(IV) specific IR harmonic frequency scale factors for both basis sets were obtained from comparison with selected experimental frequencies (C) 2010 Elsevier B V All rights reserved

A theoretical analysis on electronic structure of the (110) surface of TiO2-SnO2 mixed oxide

Mixed oxide compounds, such as TiO2-SnO2 system are widely used as gas sensors and should also provide varistor properties modifying the TiO2 surface. Therefore, a theoretical investigation has been carried out characterizing the effect of SnO2 on TiO2 addition on the electronic structure by means of ab initio SCF-LCAO calculations using all electrons. In order to take into account the finite size of the cluster, we have used the point charge model for the (TiO2)(15) cluster to study the effect on electronic structure of doping the TiO2 (110) Surface. The contracted basis set for titanium (4322/42/3), oxygen (33/3) and tin (43333/4333/43) atoms were used. The charge distributions, dipole moments, and density of states of doping TiO2 and vacancy formation are reported and analysed. (C) 2003 Elsevier B.V. All rights reserved.

Ab Initio correlated all electron dirac-fock calculations for Eka-Francium Fluoride (E119F)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

O leitor dos contos de Machado de Assis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Ab Initio correlated all electron Dirac-Fock calculations for eka-francium fluoride (E119F)

Resultados obtidos por cálculos Dirac-Fock correlacionados de 4 componentes para o fluoreto do elemento E119 (Eka-Frâncio) com base estável e precisa, livre de prolapso variacional, são reportados neste trabalho. No nível CCSD(T), a distância de equilíbrio R_e, frequência harmônica ω_e e energia de dissociação D_e são 2,432 Å, 354,97 cm_-1 e 116,92 kcal mol_-1, respectivamente. Também são reportados base livre de prolapso variacional de 4 componentes para o elemento 119, uma curva analítica de energia potencial precisa e o espectro vibracional a partir dos dados obtidos no nível CCSD(T). Nossos resultados sugerem que a molécula E119F deva ser menos iônica que seus fluoretos alcalinos homólogos mais leves, em contraste com o senso químico comum baseado nas propriedades periódicas - era de se esperar nesta molécula a ligação química mais iônica possível. Também encontramos que a correção do tipo modelo de carga para negligenciar as integrais do tipo SS resulta em erros insignificantes e acelera os cálculos cerca de 3 vezes no nível CCSD(T) e cerca de 4 vezes no nível DFT/B3LYP.

Magnetic susceptibility study of Ce3+ in PbCeA (A=Te, Se, S)

The magnetic susceptibility of Pb(1-x)Ce(x)A (A=S, Se and Te) crystals with Ce3+ concentrations 0.006 <= x <= 0.036 was investigated in the temperature range from 2 K to 300 K. The magnetic susceptibility data was found to be consistent with a E-2(5/2) lowest manifold for Ce3+ ions with a crystal-field splitting Delta=E(Gamma(8))-E(Gamma(7)) of about 340 K, 440 K and 540 K for Pb1-xCexTe, Pb1-xCexSe, and Pb1-xCexS, respectively. For all the three compounds the doublet Gamma(7) lies below the Gamma(8) quadruplet which confirms the substitution of Pb2+ by Ce3+ ions in the host crystals. The observed values for the crystal-field splitting are in good agreement with the calculated ones based on the point-charge model. Moreover, the effective Lande factors were determined by X-band (similar to 9.5 GHz), electron paramagnetic measurements (EPR) to be g=1.333, 1.364, and 1.402 for Ce ions in PbA, A = S. Se and Te, respectively. The small difference with the predicted Lande factor g of 10/7 for the Gamma(7) (J=5/2) ground state was attributed to crystal-field admixture. (C) 2012 Elsevier B.V. All rights reserved.

THE DESIGN OF A POWER SYSTEM FOR A WIRELESS SENSOR NETWORK FOR LONG-TERM OPERATION

Wireless sensor networks (WSNs) have shown wide applicability to many fields including monitoring of environmental, civil, and industrial settings. WSNs however are resource constrained by many competing factors that span their hardware, software, and networking. One of the central resource constrains is the charge consumption of WSN nodes. With finite energy supplies, low charge consumption is needed to ensure long lifetimes and success of WSNs. This thesis details the design of a power system to support long-term operation of WSNs. The power system’s development occurs in parallel with a custom WSN from the Queen’s MEMS Lab (QML-WSN), with the goal of supporting a 1+ year lifetime without sacrificing functionality. The final power system design utilizes a TPS62740 DC-DC converter with AA alkaline batteries to efficiently supply the nodes while providing battery monitoring functionality and an expansion slot for future development. Testing tools for measuring current draw and charge consumption were created along with analysis and processing software. Through their use charge consumption of the power system was drastically lowered and issues in QML-WSN were identified and resolved including the proper shutdown of accelerometers, and incorrect microcontroller unit (MCU) power pin connection. Controlled current profiling revealed unexpected behaviour of nodes and detailed current-voltage relationships. These relationships were utilized with a lifetime projection model to estimate a lifetime between 521-551 days, depending on the mode of operation. The power system and QML-WSN were tested over a long term trial lasting 272+ days in an industrial testbed to monitor an air compressor pump. Environmental factors were found to influence the behaviour of nodes leading to increased charge consumption, while a node in an office setting was still operating at the conclusion of the trail. This agrees with the lifetime projection and gives a strong indication that a 1+ year lifetime is achievable. Additionally, a light-weight charge consumption model was developed which allows charge consumption information of nodes in a distributed WSN to be monitored. This model was tested in a laboratory setting demonstrating +95% accuracy for high packet reception rate WSNs across varying data rates, battery supply capacities, and runtimes up to full battery depletion.

PVT behaviour of fluids with potential models optimised for phase equilibria

NPT and NVT Monte Carlo simulations are applied to models for methane and water to predict the PVT behaviour of these fluids over a wide range of temperatures and pressures. The potential models examined in this paper have previously been presented in the literature with their specific parameters optimised to fit phase coexistence data. The exponential-6 potential for methane gives generally good prediction of PVT behaviour over the full range of temperature and pressures studied with the only significant deviation from experimental data seen at high temperatures and pressures. The NSPCE water model shows very poor prediction of PVT behaviour, particularly at dense conditions. To improve this. the charge separation in the NSPCE model is varied with density. Improvements for vapour and liquid phase PVT predictions are achieved with this variation. No improvement was found in the prediction of the oxygen-oxygen radial distribution by varying charge separation under dense phase conditions. (C) 2004 Elsevier B.V. All rights reserved.

THE DESIGN OF A POWER SYSTEM FOR A WIRELESS SENSOR NETWORK FOR LONG-TERM OPERATION

Wireless sensor networks (WSNs) have shown wide applicability to many fields including monitoring of environmental, civil, and industrial settings. WSNs however are resource constrained by many competing factors that span their hardware, software, and networking. One of the central resource constrains is the charge consumption of WSN nodes. With finite energy supplies, low charge consumption is needed to ensure long lifetimes and success of WSNs. This thesis details the design of a power system to support long-term operation of WSNs. The power system’s development occurs in parallel with a custom WSN from the Queen’s MEMS Lab (QML-WSN), with the goal of supporting a 1+ year lifetime without sacrificing functionality. The final power system design utilizes a TPS62740 DC-DC converter with AA alkaline batteries to efficiently supply the nodes while providing battery monitoring functionality and an expansion slot for future development. Testing tools for measuring current draw and charge consumption were created along with analysis and processing software. Through their use charge consumption of the power system was drastically lowered and issues in QML-WSN were identified and resolved including the proper shutdown of accelerometers, and incorrect microcontroller unit (MCU) power pin connection. Controlled current profiling revealed unexpected behaviour of nodes and detailed current-voltage relationships. These relationships were utilized with a lifetime projection model to estimate a lifetime between 521-551 days, depending on the mode of operation. The power system and QML-WSN were tested over a long term trial lasting 272+ days in an industrial testbed to monitor an air compressor pump. Environmental factors were found to influence the behaviour of nodes leading to increased charge consumption, while a node in an office setting was still operating at the conclusion of the trail. This agrees with the lifetime projection and gives a strong indication that a 1+ year lifetime is achievable. Additionally, a light-weight charge consumption model was developed which allows charge consumption information of nodes in a distributed WSN to be monitored. This model was tested in a laboratory setting demonstrating +95% accuracy for high packet reception rate WSNs across varying data rates, battery supply capacities, and runtimes up to full battery depletion.

Atomic Charge Transfer-counter Polarization Effects Determine Infrared Ch Intensities Of Hydrocarbons: A Quantum Theory Of Atoms In Molecules Model.

Atomic charge transfer-counter polarization effects determine most of the infrared fundamental CH intensities of simple hydrocarbons, methane, ethylene, ethane, propyne, cyclopropane and allene. The quantum theory of atoms in molecules/charge-charge flux-dipole flux model predicted the values of 30 CH intensities ranging from 0 to 123 km mol(-1) with a root mean square (rms) error of only 4.2 km mol(-1) without including a specific equilibrium atomic charge term. Sums of the contributions from terms involving charge flux and/or dipole flux averaged 20.3 km mol(-1), about ten times larger than the average charge contribution of 2.0 km mol(-1). The only notable exceptions are the CH stretching and bending intensities of acetylene and two of the propyne vibrations for hydrogens bound to sp hybridized carbon atoms. Calculations were carried out at four quantum levels, MP2/6-311++G(3d,3p), MP2/cc-pVTZ, QCISD/6-311++G(3d,3p) and QCISD/cc-pVTZ. The results calculated at the QCISD level are the most accurate among the four with root mean square errors of 4.7 and 5.0 km mol(-1) for the 6-311++G(3d,3p) and cc-pVTZ basis sets. These values are close to the estimated aggregate experimental error of the hydrocarbon intensities, 4.0 km mol(-1). The atomic charge transfer-counter polarization effect is much larger than the charge effect for the results of all four quantum levels. Charge transfer-counter polarization effects are expected to also be important in vibrations of more polar molecules for which equilibrium charge contributions can be large.