How to mesh up Ewald sums (I): A theoretical and numerical comparison of various particle mesh routines

Standard Ewald sums, which calculate e.g. the electrostatic energy or the force in periodically closed systems of charged particles, can be efficiently speeded up by the use of the Fast Fourier Transformation (FFT). In this article we investigate three algorithms for the FFT-accelerated Ewald sum, which attracted a widespread attention, namely, the so-called particle-particle-particle-mesh (P3M), particle mesh Ewald (PME) and smooth PME method. We present a unified view of the underlying techniques and the various ingredients which comprise those routines. Additionally, we offer detailed accuracy measurements, which shed some light on the influence of several tuning parameters and also show that the existing methods -- although similar in spirit -- exhibit remarkable differences in accuracy. We propose combinations of the individual components, mostly relying on the P3M approach, which we regard as most flexible.Comment: 18 pages, 8 figures included, revtex styl

The optimal P3M algorithm for computing electrostatic energies in periodic systems

We optimize Hockney and Eastwood's Particle-Particle Particle-Mesh (P3M) algorithm to achieve maximal accuracy in the electrostatic energies (instead of forces) in 3D periodic charged systems. To this end we construct an optimal influence function that minimizes the RMS errors in the energies. As a by-product we derive a new real-space cut-off correction term, give a transparent derivation of the systematic errors in terms of Madelung energies, and provide an accurate analytical estimate for the RMS error of the energies. This error estimate is a useful indicator of the accuracy of the computed energies, and allows an easy and precise determination of the optimal values of the various parameters in the algorithm (Ewald splitting parameter, mesh size and charge assignment order).Comment: 31 pages, 3 figure

Selected new developments in computational chemistry.

Molecular dynamics is a general technique for simulating the time-dependent properties of molecules and their environments. Quantum mechanics, as applied to molecules or clusters of molecules, provides a prescription for predicting properties exactly (in principle). It is reasonable to expect that both will have a profound effect on our understanding of environmental chemistry in the future. In this review, we consider several recent advances and applications in computational chemistry

Molecular modeling studies suggest that zinc ions inhibit HIV-1 protease by binding at catalytic aspartates.

Human immunodeficiency virus type 1 protease is inhibited in vitro by zinc ions at neutral pH. The binding site of these ions is not known; however, experimental data suggest that binding may occur in the active site. To examine the possibility of zinc binding in the active site, molecular dynamics simulations in the presence and absence of zinc have been carried out to 200 psec. The results are compared with the 2.8-A crystallographic structures of a synthetic HIV-1 protease, and a zinc binding site at the catalytic aspartate residues (Asp-25, Asp-25') is proposed. Molecular dynamics simulations show that the zinc ion remains stably bound in this region, coordinating the carboxylate side chains of both aspartate residues. Interaction with zinc does not disrupt the dimeric structure of the protein or significantly alter the structure of the active site. These data are consistent with experimental studies of HIV-1 protease inhibition by zinc and give strong evidence that this is the binding site that leads to inactivation

Potentiality in Biology

We take the potentialities that are studied in the biological sciences (e.g., totipotency) to be an important subtype of biological dispositions. The goal of this paper is twofold: first, we want to provide a detailed understanding of what biological dispositions are. We claim that two features are essential for dispositions in biology: the importance of the manifestation process and the diversity of conditions that need to be satisfied for the disposition to be manifest. Second, we demonstrate that the concept of a disposition (or potentiality) is a very useful tool for the analysis of the explanatory practice in the biological sciences. On the one hand it allows an in-depth analysis of the nature and diversity of the conditions under which biological systems display specific behaviors. On the other hand the concept of a disposition may serve a unificatory role in the philosophy of the natural sciences since it captures not only the explanatory practice of biology, but of all natural sciences. Towards the end we will briefly come back to the notion of a potentiality in biology

Volume-energy correlations in the slow degrees of freedom of computer-simulated phospholipid membranes

Constant-pressure molecular-dynamics simulations of phospholipid membranes in the fluid phase reveal strong correlations between equilibrium fluctuations of volume and energy on the nanosecond time-scale. The existence of strong volume-energy correlations was previously deduced indirectly by Heimburg from experiments focusing on the phase transition between the fluid and the ordered gel phases. The correlations, which are reported here for three different membranes (DMPC, DMPS-Na, and DMPSH), have volume-energy correlation coefficients ranging from 0.81 to 0.89. The DMPC membrane was studied at two temperatures showing that the correlation coefficient increases as the phase transition is approached

Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation

In this paper, we investigate the conformational dynamics of alanine dipeptide under an external electric field by nonequilibrium molecular dynamics simulation. We consider the case of a constant and of an oscillatory field. In this context we propose a procedure to implement the temperature control, which removes the irrelevant thermal effects of the field. For the constant field different time-scales are identified in the conformational, dipole moment, and orientational dynamics. Moreover, we prove that the solvent structure only marginally changes when the external field is switched on. In the case of oscillatory field, the conformational changes are shown to be as strong as in the previous case, and non-trivial nonequilibrium circular paths in the conformation space are revealed by calculating the integrated net probability fluxes.Comment: 23 pages, 12 figure

Effects of daily porcine somatotropin administration on tile performance and carcass characteristics of finishing swine

One hundred forty-four pigs (72 barrows and 72 gilts) were utilized in six treatments with six pens each to evaluate four levels of daily porcine somatotropin (pST) injections (0, 1, 3, or 5 mg/d) in combination with diets containing 13 or 16% crude protein (CP). One randomly selected pig from each pen was slaughtered when it reached a weight between 230 and 240 lb. Daily feed intake (ADFI), tenth rib backfat, and estimated percentage of lean pork were reduced in pigs fed the 13% CP diet and injected with 3 mg/d pST compared to pigs fed the same diet and injected daily with a placebo. Feed intake was reduced in pigs fed the 16% CP diet and injected daily with 3 and 5 mg/d pST. Improvement in feed conversion (F/G), tenth rib backfat, and estimated percentage lean pork occurred when pigs fed the 16% CP diet were injected with 1 mg/d, with greater improvements occurring at the 3 and 5 mg/d levels of pST. The improvement in F/G and the greater magnitude of response observed in pigs fed 16% CP diets compared to the response of pigs fed the 13% CP diet indicate that both the performance and carcass characteristics of pigs injected with pST are dependent on the dietary CP content.; Swine Day, Manhattan, KS, November 16, 198