Biomedical subjects

Donald L Thompson

Publications and source records attributed to Donald L Thompson.

23 records · Page 2Linked to original sources

PubMed2004

Interpolating moving least-squares methods for fitting potential energy surfaces: applications to classical dynamics calculations.

As a continuation of our efforts to develop efficient and accurate interpolating moving least-squares (IMLS) methods for generating potential energy surfaces, we carry out classical trajectories and compute kinetics properties on higher degree IMLS surfaces. In this study, we have investigated the choice of coordinate system, the range of points (i.e., the cutoff radius) used in fitting, and strategies for selections of data points and basis elements. We illustrate and test the method by applying it to hydrogen peroxide (HOOH). In particular, reaction rates for the O-O bond breaking in HOOH are calculated on fitted surfaces using the classical trajectory approach to test the accuracy of the IMLS method for providing potentials for dynamics calculations.

Journal Article↗

PubMed2004

A molecular-dynamics study of structural and physical properties of nitromethane nanoparticles.

The structural and physical properties of nanoparticles of nitromethane are studied by using molecular dynamics methods with a previously developed force field. [Agrawal et al., J. Chem. Phys. 119, 9617 (2003).] This force field accurately predicts solid- and liquid-state properties as well as melting of bulk nitromethane. Molecular dynamics simulations of nanoparticles with 480, 240, 144, 96, 48, and 32 nitromethane molecules have been carried out at various temperatures. The carbon-carbon radial distribution function, dipole-dipole correlation function, core density, internal enthalpy, and atomic diffusion coefficients of the nanoparticles were calculated at each temperature. These properties were used to characterize the physical phases and thus determine the melting transitions of the nanoparticles. The melting temperatures predicted by the various properties are consistent with one another and show that the melting temperature increases with particle size, approaching the bulk limit for the largest particle. A size dependence of melting points has been observed in experimental and theoretical studies of atomic nanoparticles, and this is a further demonstration of the effect for large nanoparticles of complex molecular materials.

Journal Article↗

PubMed2004

Molecular dynamics studies of melting and solid-state transitions of ammonium nitrate.

Molecular dynamics simulations are used to calculate the melting point and some aspects of high-temperature solid-state phase transitions of ammonium nitrate (AN). The force field used in the simulations is that developed by Sorescu and Thompson [J. Phys. Chem. A 105, 720 (2001)] to describe the solid-state properties of the low-temperature phase-V AN. Simulations at various temperatures were performed with this force field for a 4 x 4 x 5 supercell of phase-II AN. The melting point of AN was determined from calculations on this supercell with voids introduced in the solid structure to eliminate superheating effects. The melting temperature was determined by calculating the density and the nitrogen-nitrogen radial distribution functions as functions of temperature. The melting point was predicted to be in the range 445 +/- 10 K, in excellent agreement with the experimental value of 442 K. The computed temperature dependences of the density, diffusion, and viscosity coefficient for the liquid are in good agreement with experiment. Structural changes in the perfect crystal at various temperatures were also investigated. The ammonium ions in the phase-II structure are rotationally disordered at 400 K. At higher temperatures, beginning at 530 K, the nitrate ions are essentially rotationally unhindered. The density and radial distribution functions in this temperature range show that the AN solid is superheated. The rotational disorder is qualitatively similar to that observed in the experimental phase-II to phase-I solid-state transition.

Journal Article↗

PubMed2004

Improving the accuracy of interpolated potential energy surfaces by using an analytical zeroth-order potential function.

We present a method for improving the accuracy and efficiency of interpolation methods, in which an analytical zeroth-order potential-energy surface is employed as a reference surface. To investigate and test the method, we apply it to hydrogen peroxide where there exists an accurate analytical surface which we take as the "exact" surface for obtaining the energies and derivatives for fitting and assessing the accuracy. Examples are given for four-dimensional and six-dimensional surfaces interpolated by using either the modified Shepard or second-degree interpolating moving least-squares approach, with comparisons for cases with and without using the zeroth-order potential.

Journal Article↗

PubMed2004

On combining molecular dynamics and stochastic dynamics simulations to compute reaction rates in liquids.

An approach that combines molecular dynamics and stochastic dynamics calculations for obtaining reaction rates in liquids is investigated by studying the cis-->trans isomerization of HONO in liquid krypton. The isomerization rates are computed for several liquid densities by employing full-dimensional molecular-dynamics simulations. The rates are also computed by employing the stochastic dynamics method for a wide range of collision frequencies. Comparisons of the two sets of the computed rates show that for a wide range of liquid densities there is a simple linear relation between the liquid density rho and the collision frequency alpha, that is, alpha=crho. This suggests that once the constant c is determined from a molecular-dynamics calculation at a single density, the reaction rates can be obtained from stochastic dynamics calculations for the entire range of liquid densities where alpha=crho holds. The applicability of the combined molecular dynamics and stochastic dynamics approach provides a practical means for obtaining rate constants at considerable savings of computer time compared to that required by using full-dimensional molecular-dynamics simulations alone.

Journal Article↗