MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Theoretical and Experimental Studies in Nuclear Magnetic Resonance

Nuclear magnetic resonance (NMR) is a tool used to probe the physical and chemical environments of specific atoms in molecules. This research explored small molecule analogues to biological materials to determine NMR parameters using ab initio computations, comparing the results with solid-state NMR measurements. Models, such as dimethyl phosphate (DMP) for oligonucleotides or CuCl for the active site of the protein azurin, represented computationally unwieldy macromolecules. 31P chemical shielding tensors were calculated for DMP as a function of torsion angles, as well as for the phosphate salts, ammonium dihydrogen phosphate (ADHP), diammonium hydrogen phosphate, and magnesium dihydrogen phosphate. The computational DMP work indicated a problem with the current standard 31P reference of 85% H3PO4(aq.). Comparison of the calculations and experimental spectra for the phosphate salts indicated ADHP might be a preferable alternative as a solid state NMR reference for 31P. Experimental work included magic angle spinning experiments on powder samples using the UNL chemistry department’s Bruker Avance 600 MHz NMR to collect data to determine chemical shielding anisotropies. For the quadrupolar nuclei of copper and scandium, the electric field gradient was calculated in diatomic univalent metal halides, allowing determination of the minimal level of theory necessary to compute NMR parameters for these nuclei.