MD models

A series of articles were published by Ennari et al. on MD simulation of transport processes in Poly(Ethylene Oxide) and sulfonic acid-based polymer electrolyte.136,137 The work was started by the determination of the parameters for the ions missing from the PCFF forcefield made by MSI (Molecular Simulations Inc.), to create a new forcefield, NJPCFF. In the models, the proton is represented as a hard ball with a positive charge. Zhou et al. used the similar approach to model Nation.138 The repeating unit of Nafion (Fig. 17) was optimized using ab initio VAMP scheme. The protons were modeled with hydronium ions. Three unit cell or molecular models were used for the MD simulation. The unit cell contains 5000 atoms 20 pendent side chains, and branched Nafion backbone created with the repeating unit. Their water uptakes or water contents were 3, 13, or 22 IEO/SO3, which correspond to the room temperature water uptakes at 50% relative humidity (RH), at 100% RH, and in liquid water respectively.18 The temperature was initially set at a value between 298.15 and 423.15 K under NVE ensemble with constant particle number, constant volume (1 bar), and constant energy. 

However, the final temperature was usually different from the initial temperature. The initial velocity was randomly determined according to Boltzmann statistics. The time-step was set at 1.0 fs and duration was 100 ps. Each run of the molecular dynamic (MD) computation lasts about 72 h using a Dell Pentium IV personal computer. The position vector of each particle, Ri(t), is saved as a function of time. The diffusivity of a specific molecule or ion was evaluated using136,137  

As shown in Fig. 18, the simulation conductivity data were generally consistent with the experimental results. However, there are appreciable differences between the simulation and experimental results. At some points, the differences can be 100%. The other observation is that the simulation must underestimate the activation energies of the conduction. The primary reason for this discrepancy is that these simulation models do not take into account interaction between the membrane itself and its environment. In reality, the water uptake at elevated temperatures may be greater than that at room temperature. In the simulations, it was assumed that both 

The simulation models also correctly predicted the diffusivities of hydronium and methanol in a wide range of temperature (Fig. 19). Methanol is a neutral species and weakly interacts with Nation backbone. It is not surprising that the present MD models that do not consider chemical interaction between the molecules can still correctly evaluate the diffusivity of methanol. Because the present experimental setup is limited for liquid samples, whether or not the permeability of diffusivity is strongly depends on water content has not been examined. In summary, this work provided benchmark for the atomistic simulation of the transport processes in Nation at water content above 3 although at some points, the errors can be 100%. 

---

We now describe a relatively simple MD model of a low-index crystal surface, which was conceived for the purpose of studying the rate of mass transport (8). The effect of temperature on surface transport involves several competing processes. A rough surface structure complicates the trajectories somewhat, and the diffusion of clusters of atoms must be considered. In order to simplify the model as much as possible, but retain the essential dynamics of the mobile atoms, we will consider a model in which the atoms move on a "substrate" represented by an analytic potential energy function that is adjusted to match that of a surface of a (100) face-centered cubic crystal composed of atoms interacting with a Lennard-Jones... 

We present and discuss results for MD modeling of fluid systems. We restrict our discussion to systems which are in a macroscopically steady state, thus eliminating the added complexity of any temporal behavior. We start with a simple fluid system where the hydrodynamic equations are exactly solvable. We conclude with fluid systems for which the hydrodynamic equations are nonlinear. Solutions for these equations can be obtained only through numerical methods. 

Uehing MD. Model patient. Bio-IT World, November 1, 2005, http //www.bio-itworld.com/archive/121503/trials.html... 

Further progress in understanding membrane instability and nonlocality requires development of microscopic theory and modeling. Analysis of membrane thickness fluctuations derived from molecular dynamics simulations can serve such a purpose. A possible difficulty with such analysis must be mentioned. In a natural environment isolated membranes assume a stressless state. However, MD modeling requires imposition of special boundary conditions corresponding to a stressed state of the membrane (see Refs. 84,87,112). This stress can interfere with the fluctuations of membrane shape and thickness, an effect that must be accounted for in analyzing data extracted from computer experiments. 

We also emphasize that the MD model does include the vibrational motions of bond, and torsional angles (in the minima of the respective potentials) but, somehow, these small scale fast motions are rapidly damped out in the melt, and do not affect the motion on the nanometer scale (and for corresponding times) significantly. 

As mentioned earlier, in principle, one can model the dynamics of a simple classical fluid by means of MD simulations. This technique, although straightforward, is relatively time-consuming, and therefore not suitable for observation of large-scale macroscopic phenomena in the fluid. However, one often does not need such a detailed description of the microdynamics as provided by MD. In such cases, it would be more efficient to strip the MD model down to its barest essentials, where the only requirement is that the model behaves like a fluid macroscopically, but is still atomistic in character—i.e., the mechanism underlying the fluid motion is the movement of particles. From the derivation of... 

MDM hydantoin, antimicrobial used in cosmetics, 7 831t, 832 MD modeling, 12 576—577 Mean cell residence time (MCRT), in biological waste treatment, 25 830 Mean centering, 6 35—38 Mean diameters, for statistical properties of droplets, 23 186... 

In (20) the titration of a model Fe3+ ion in solution was simulated with the Fe3+-0-H MD model described above. Because of the way the model is designed, allowing hydrolysis reactions to occur spontaneously, the distribution of the M(OH) species can simply be recorded as a function of the amount of protons and hydroxide ions added to the solution. The simulations are used to find the distribution of hydrolysis species in a neutral solution and to determine how this distribution responds to added protons or hydroxide ions in solution. 

Fig. 7. Fe3(0H)7(H20)62+ trimer used to compare gas-phase acidities of m3-OH, groups between density functional theory and the MD model. The MD model and GGA DFT calculations both predict a gas-phase acidity of 179 kcal/mol.

MD Modeling of Interface Structure on Catalysts in Aqueous Solution... 

MD Modeling of Interface Structure of Polymer Electrolyte Catalysts Interface... 

Goddard et al.18 are currently carrying out exploratory calculations on both OER and proton conduction processes using the ReaXFF MD method. They anticipate that this will establish a very realistic MD model of a tiny fuel cell and will allow determining how the performance changes as the composition, configuration, and other conditions are altered. 

Atomistic MD models can be extended to the coarse-grained level introduced in the previous section, which is determined by the dimension of the backbone chain and branch. For the precise description of water molecular behavior, simple point charge (SPC) model was adopted (Krishnan et al., 2001), which can be used to simulate complex composition systems and quantitatively express vibrational spectra of water molecules in vapor, liquid, and solid states. The six-parameter (Doh, o , fi, Lye, Lyy, and Lee) SPC potential used for the water molecules is shown in Equation (24) ... 

Fiber spinning Finite element modeling — MD modeling... 

To estimate the relaxation time of the adsorbed complex, MD modeling of ZrCl4 adsorption on the hydroxylated Zr02 surface was performed. The... 

A periodic model of the zirconia surface was used in MD modeling. It consisted of 2 x 2 elementary Zr02 cells, and ZrCl4 precursors incident perpendicular to the surface. The Zr02 surface was connected with a thermostat with a temperature of 600 K, and the initial velocities of the ZrCl4 precursor corresponded to a temperature of 600 K. 

Multi-well modified Boyden chambers can be obtained from Neuro Probe Inc. (Gaithersburg, MD). Model AP48 has been widely used for endothelial cell chemotaxis assays. The apparatus consists of top and bottom acrylic plates, a silicon gasket and assembly screws. The bottom plate has 48 wells, each with 25 pL final volume. These correspond to holes on the top plate, and form the upper wells when the chamber is assembled. The filter (polycarbonate, 25 x 80 mm) is placed between the top and bottom plates, and a gasket is placed over the filter to create the seal. The apparatus can be purchased with a selection of accessories, such as curved forceps, filter clamps, and wipers. These are required to process the filters after use. Filters can also be obtained from Nucleopore Inc. (Pleasanton, CA) and Costar (Cambridge, MA). 

Recently, Garth Jones (UNSW) and his co-workers have described a semi-classical molecular dynamics (MD) model, which incorporates the trajectory surface hopping (TSH) method, that shows promise in investigating a number of important problems concerning the effects of molecular vibrations on the dynamics of ET reactions, including formally symmetry-forbidden ones. Essentially,... 

---