Thesis (M.S.)--Miami University, Dept. of Physics, 2008. / Title from first page of PDF document. Includes bibliographical references (p. 36-37).

The typical form for a molecular mechanics force field consists of a foundation of pair-wise terms to describe bonded and non-bonded atomic interactions, with multi-body correction terms to deal with the limitations of pair-wise terms. I present here the first attempts of a molecular mechanics model that is founded on multi-body terms, which are based on the Bond Valence Model (Brown, 2002) and recent developments in the Vectorial Bond Valence Model (Bickmore et al., 2013a; Harvey et al., 2006). I calibrated these models on pressure vs. energy curves for a set of SiO2 polymorphs. The average deviation for the best-fit iteration, with only six adjustable parameters was ±1.98 kcal/mol.

No description available.

A química computacional apresenta a grande vantagem de prover informações fundamentais para espécies moleculares propostas, antes mesmo de sua síntese em laboratório. A Teoria do Funcional da Densidade é bastante utilizada nesta área, produzindo resultados satisfatórios para um grande número de propriedades e sistemas, mas com uma menor demanda por recursos computacionais que métodos mais avançados. Entretanto, o desenvolvimento de funcionais que incluem efeitos relativísticos ainda se encontra num estágio inicial. Em geral, tais efeitos são importantes em compostos de átomos pesados, embora devam ser considerados também em sistemas com átomos mais leves se a propriedade em estudo for particularmente sensível, como é o caso do gradiente de campo elétrico na posição de núcleos em moléculas. Assim, na primeira etapa desta dissertação foi avaliado o desempenho de funcionais comuns de troca-correlação não relativísticos, quando utilizados em conjunto com o formalismo de quatro componentes (tratamento relativístico), no estudo dos gradientes de campo elétrico em núcleos de átomos (índio, antimônio, iodo, lutécio e háfnio) constituindo moléculas diatômicas. Foram investigados funcionais baseados nas aproximações da densidade local e do gradiente generalizado, funcionais híbridos e que incluem correções em termos da atenuação com a distância. Nossos resultados, que estão em acordo com observações da literatura, ressaltam o melhor desempenho de funcionais híbridos e com correções de atenuação para esta propriedade e demonstram a importância do uso do método indireto. Posteriormente, foi feita uma nova parametrização de alguns dos melhores funcionais não relativísticos selecionados na etapa anterior (B3LYP, PBE0 e CAM-B3LYP), dentro do formalismo de quatro componentes, para uso no cálculo destes mesmos gradientes num grupo teste de átomos (cobre, iodo, lantânio e ouro) em moléculas lineares. Nestes casos, os funcionais modificados propostos tiveram um bom desempenho geral e foram particularmente bem sucedidos para cobre e ouro. Finalmente, é possível destacar o funcional híbrido PBE0 e sua modificação, proposta neste estudo, por conta de seu desempenho excelente, tanto para os metais como para os demais elementos que tiveram seus EFGs investigados aqui. / The computational chemistry has the great advantage of providing fundamental information for proposed molecular species even before their synthesis in laboratory. The Density Functional Theory is widely used in this area, producing satisfactory results for a large number of properties and systems, but with a lower demand for computational resources than that of more advanced methods. However, the development of functionals that include relativistic effects is still at an early stage. In general, these effects are important in compounds containing heavy elements, but they must also be considered in systems of lighter atoms if the studied property was particularly sensitive, as occurs for the electric field gradient at the position of nuclei in molecules. Thus, the first step of this dissertation was to evaluate the performance of common non-relativistic exchange-correlation functionals when used in conjunction with the four component formalism (relativistic treatment) in the study of electric field gradients at the nuclei of atoms (indium, antimony, iodine, lutetium and hafnium) forming diatomic molecules. Functionals based on the local density approximation and generalized gradient approximation, hybrid functionals and the ones that include attenuation corrections were investigated. Our results, which are in agreement with observations in the literature, highlight the best performance of hybrid functionals and attenuation corrections for this property and demonstrate the importance of using the indirect approach. Subsequently, there was a new parameterization of some of the best non-relativistic functionals selected in the previous step (B3LYP, PBE0 and CAM - B3LYP) within the four component formalism for calculations of these same gradients in a trial group of atoms (copper, iodine, lanthanum and gold) into linear molecules. In these cases, the modified functionals proposed had a satisfactory overall performance and were particularly successful for copper and gold. Finally, it is possible to mention the excellent performance of the hybrid functional PBE0 and its modification proposed in this study for both metals and the other elements that had their EFGs investigated here.

Abstract The studies of the deformed nucleus for the Pb isotopes in the A~190 re- gion has been carried out by using the P.S.M. The comparison between the theoretical calculation of the transition energy, kinetic moment of ine- rtia, dynamic moment of inertia, paring gap, rotational energy, quadrupole moment, and B(M1)/B(E2) and the experiment observations are presented. We found that the kinetic moment of inertia is increased as the spin incre- ase for the even-even nuclei. This rise has been interpreted as being due to the combined alignment of the quasiparticles in the high-N intruder or- bit. And the kinetic moment of inertia is constant function of spin in the odd-even nuclei. It is interpreted as been due to the Pauli blocking of odd -quasineutron in the high-N intruder orbit.

A química computacional apresenta a grande vantagem de prover informações fundamentais para espécies moleculares propostas, antes mesmo de sua síntese em laboratório. A Teoria do Funcional da Densidade é bastante utilizada nesta área, produzindo resultados satisfatórios para um grande número de propriedades e sistemas, mas com uma menor demanda por recursos computacionais que métodos mais avançados. Entretanto, o desenvolvimento de funcionais que incluem efeitos relativísticos ainda se encontra num estágio inicial. Em geral, tais efeitos são importantes em compostos de átomos pesados, embora devam ser considerados também em sistemas com átomos mais leves se a propriedade em estudo for particularmente sensível, como é o caso do gradiente de campo elétrico na posição de núcleos em moléculas. Assim, na primeira etapa desta dissertação foi avaliado o desempenho de funcionais comuns de troca-correlação não relativísticos, quando utilizados em conjunto com o formalismo de quatro componentes (tratamento relativístico), no estudo dos gradientes de campo elétrico em núcleos de átomos (índio, antimônio, iodo, lutécio e háfnio) constituindo moléculas diatômicas. Foram investigados funcionais baseados nas aproximações da densidade local e do gradiente generalizado, funcionais híbridos e que incluem correções em termos da atenuação com a distância. Nossos resultados, que estão em acordo com observações da literatura, ressaltam o melhor desempenho de funcionais híbridos e com correções de atenuação para esta propriedade e demonstram a importância do uso do método indireto. Posteriormente, foi feita uma nova parametrização de alguns dos melhores funcionais não relativísticos selecionados na etapa anterior (B3LYP, PBE0 e CAM-B3LYP), dentro do formalismo de quatro componentes, para uso no cálculo destes mesmos gradientes num grupo teste de átomos (cobre, iodo, lantânio e ouro) em moléculas lineares. Nestes casos, os funcionais modificados propostos tiveram um bom desempenho geral e foram particularmente bem sucedidos para cobre e ouro. Finalmente, é possível destacar o funcional híbrido PBE0 e sua modificação, proposta neste estudo, por conta de seu desempenho excelente, tanto para os metais como para os demais elementos que tiveram seus EFGs investigados aqui. / The computational chemistry has the great advantage of providing fundamental information for proposed molecular species even before their synthesis in laboratory. The Density Functional Theory is widely used in this area, producing satisfactory results for a large number of properties and systems, but with a lower demand for computational resources than that of more advanced methods. However, the development of functionals that include relativistic effects is still at an early stage. In general, these effects are important in compounds containing heavy elements, but they must also be considered in systems of lighter atoms if the studied property was particularly sensitive, as occurs for the electric field gradient at the position of nuclei in molecules. Thus, the first step of this dissertation was to evaluate the performance of common non-relativistic exchange-correlation functionals when used in conjunction with the four component formalism (relativistic treatment) in the study of electric field gradients at the nuclei of atoms (indium, antimony, iodine, lutetium and hafnium) forming diatomic molecules. Functionals based on the local density approximation and generalized gradient approximation, hybrid functionals and the ones that include attenuation corrections were investigated. Our results, which are in agreement with observations in the literature, highlight the best performance of hybrid functionals and attenuation corrections for this property and demonstrate the importance of using the indirect approach. Subsequently, there was a new parameterization of some of the best non-relativistic functionals selected in the previous step (B3LYP, PBE0 and CAM - B3LYP) within the four component formalism for calculations of these same gradients in a trial group of atoms (copper, iodine, lanthanum and gold) into linear molecules. In these cases, the modified functionals proposed had a satisfactory overall performance and were particularly successful for copper and gold. Finally, it is possible to mention the excellent performance of the hybrid functional PBE0 and its modification proposed in this study for both metals and the other elements that had their EFGs investigated here.

In this thesis, we study the phenomenon of thermally induced polarisation using a combination of theory and computer simulation. Molecules of sufficiently low symmetry exhibit thermo-molecular orientation when subjected to a temperature gradient, leading to considerable electrostatic fields in polar liquids. Here, we first use non-equilibrium molecular dynamics simulations to study this interesting effect numerically. To this end, we propose an integration algorithm to impose a constant heat flux in simulations and show that it greatly improves energy conservation compared to a previous algorithm. We next investigate the thermal polarisation of water and find that truncation of electrostatic interactions can lead to severe artefacts, such as the wrong sign of polarisation and an overestimation of the electric field. We further show that the quadrupole-moment contribution to the electric field is significant and responsible for an inversion of its sign. To facilitate the theoretical description of electrostatic interactions, we propose a new dipolar model fluid as a perturbation of a Stockmayer fluid. Using this modified Stockmayer model, we provide numerical evidence for the recently proposed phenomenon of thermally induced monopoles. We show that the electrostatic field generated by a pair of heated/cooled colloidal particles immersed in such a solvent can be trivially described by two Coulomb charges. Finally, we propose a mean-field theory to predict the thermo-polarisation effect exhibited by our model fluid theoretically, and demonstrate near quantitative agreement with simulation results.

In the first part of the thesis, the molecular structure and electronic properties of phospholipids at the single molecule level and also for a monolayer structure are investigated via ab initio calculations under different degrees of hydration. The focus of the study is on phosphatidylcholines, in particular dipalmitoylphosphatidylcholine (DPPC), which are the most abundant phospholipids in biological membranes. Upon hydration, the phospholipid shape into a sickle-like structure. The hydration dramatically alters the surface potential, dipole and quadrupole moments of the lipids, and probably guides the interactions of the lipids with other molecules and the communication between cells. The vibrational spectrum of DPPC and DPPC-water complexes are completely assigned and it is shown that water hydrating the lipid head groups enables efficient energy transfer across membrane leaflets on sub-picosecond time scales. Moreover, the vibrational modes and lifetimes of pure and hydrated DPPC lipids, at human body temperature, are estimated by performing ab initio molecular dynamics simulations. The vibrational modes of the water molecules close to the head group of DPPC are active in the frequency range between 0.5 - 55 THz, with a peak at 2.80 THz in the energy spectrum. The computed lifetimes for the high-frequency modes agree well with recent data measured at room temperature, where high-order phonon scattering is not negligible. The structure and auto-ionization of water at the water-phospholipid interface are investigated by ab initio molecular dynamics and ab initio Monte Carlo simulations using local density approximation and generalized gradient approximation for the exchange-correlation energy functional. Depending on the lipid head group, strongly enhanced ionization is observed, leading to dissociation of several water molecules into H+ and OH- per lipid. The results can shed light on the phenomena of the high proton conductivity along membranes that has been reported experimentally. In the second part of the thesis, Monte Carlo simulations of the lipid bilayer, on the basis of a coarse grained model, are performed to gain insight into the mechanical properties of planar lipid bilayers. By using a rescaling method, the Poisson's ratio is calculated for different phases. Additional information on the bending rigidity, determined from height fluctuations on the basis of the Helfrich Hamiltonian, allows for calculation of the Young's modulus for each phase. In addition, the free energy barrier for lipid flip-flop process in the fluid and gel phases are estimated. The main rate-limiting step to complete a flip-flop process is related to a free energy barrier that has to be crossed in order to reach the center of the bilayer. The free energy cost for performing a lipid flip-flop in the gel phase is found to be five times greater than in the fluid phase, demonstrating the rarity of such events in the gel phase. Moreover, an energy barrier is estimated for formation of transient water pores that often precedes lipid translocation events and accounts for the rate-limiting step of these pore-associated lipid translocation processes.