Amino-acid-dependent main-chain torsion-energy terms for protein systems

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters of main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for $\psi$ and $\psi'$ angles for each amino acid by using our optimization method, which is one of the knowledge-based approach. In order to test the validity of the new force-field parameters, we then performed folding simulations of $\alpha$-helical and $\beta$-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff03 force field.Comment: 10 pages, (Revtex4.1), 3 tables, 6 figure

Mapping atomistic to coarse-grained polymer models using automatic simplex optimization to fit structural properties

We develop coarse-grained force fields for poly (vinyl alcohol) and poly (acrylic acid) oligomers. In both cases, one monomer is mapped onto a coarse-grained bead. The new force fields are designed to match structural properties such as radial distribution functions of various kinds derived by atomistic simulations of these polymers. The mapping is therefore constructed in a way to take into account as much atomistic information as possible. On the technical side, our approach consists of a simplex algorithm which is used to optimize automatically non-bonded parameters as well as bonded parameters. Besides their similar conformation (only the functional side group differs), poly (acrylic acid) was chosen to be in aqueous solution in contrast to a poly (vinyl alcohol) melt. For poly (vinyl alcohol) a non-optimized bond angle potential turns out to be sufficient in connection with a special, optimized non-bonded potential. No torsional potential has to be applied here. For poly (acrylic acid), we show that each peak of the radial distribution function is usually dominated by some specific model parameter(s). Optimization of the bond angle parameters is essential. The coarse-grained forcefield reproduces the radius of gyration of the atomistic model. As a first application, we use the force field to simulate longer chains and compare the hydrodynamic radius with experimental data.Comment: 34 pages, 3 tables, 16 figure

Optimization of classical nonpolarizable force fields for OH− and H3O+

We optimize force fields for H3O+ and OH− that reproduce the experimental solvation free energies and the activities of H3O+ Cl− and Na+ OH− solutions up to concentrations of 1.5 mol/l. The force fields are optimized with respect to the partial charge on the hydrogen atoms and the Lennard-Jones parameters of the oxygen atoms. Remarkably, the partial charge on the hydrogen atom of the optimized H3O+ force field is 0.8 ± 0.1|e|—significantly higher than the value typically used for nonpolarizable water models and H3O+ force fields. In contrast, the optimal partial charge on the hydrogen atom of OH− turns out to be zero. Standard combination rules can be used for H3O+ Cl− solutions, while for Na+ OH− solutions, we need to significantly increase the effective anion- cation Lennard-Jones radius. While highlighting the importance of intramolecular electrostatics, our results show that it is possible to generate thermodynamically consistent force fields without using atomic polarizability

The dynamics of copper intercalated molybdenum ditelluride

Layered transition metal dichalcogenides are emerging as key materials in nanoelectronics and energy applications. Predictive models to understand their growth, thermomechanical properties and interactions with metals are needed in order to accelerate their incorporation into commercial products. Interatomic potentials enable large-scale atomistic simulations at the device level, beyond the range of applications of first principle methods. We present a ReaxFF reactive force field to describe molybdenum ditelluride and its interactions with copper. We optimized the force field parameters to describe the properties of layered MoTe2 in various phases, the intercalation of Cu atoms and clusters within its van der Waals gap, including a proper description of energetics, charges and mechanical properties. The training set consists of an extensive set of first principle calculations computed from density functional theory. We use the force field to study the adhesion of a single layer MoTe2 on a Cu(111) surface and the results are in good agreement with density functional theory, even though such structures were not part of the training set. We characterized the mobility of the Cu ions intercalated into MoTe2 under the presence of an external electric fields via molecular dynamics simulations. The results show a significant increase in drift velocity for electric fields of approximately 0.4 V/A and that mobility increases with Cu ion concentration.Comment: 21 pages, 9 Figure

Mapping atomistic to coarse-grained polymer models using automatic simplex optimization to fit structural properties

We develop coarse-grained force fields for poly (vinyl alcohol) and poly (acrylic acid) oligomers. In both cases, one monomer is mapped onto a coarse-grained bead. The new force fields are designed to match structural properties such as radial distribution functions of various kinds derived by atomistic simulations of these polymers. The mapping is therefore constructed in a way to take into account as much atomistic information as possible. On the technical side, our approach consists of a simplex algorithm which is used to optimize automatically non-bonded parameters as well as bonded parameters. Besides their similar conformation (only the functional side group differs), poly (acrylic acid) was chosen to be in aqueous solution in contrast to a poly (vinyl alcohol) melt. For poly (vinyl alcohol) a non-optimized bond angle potential turns out to be sufficient in connection with a special, optimized non-bonded potential. No torsional potential has to be applied here. For poly (acrylic acid), we show that each peak of the radial distribution function is usually dominated by some specific model parameter(s). Optimization of the bond angle parameters is essential. The coarse-grained forcefield reproduces the radius of gyration of the atomistic model. As a first application, we use the force field to simulate longer chains and compare the hydrodynamic radius with experimental data.Comment: 34 pages, 3 tables, 16 figure

Strong scaling of general-purpose molecular dynamics simulations on GPUs

We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, arXiv:1308.5587). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, J. Comp. Phys. 117, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., J. Comp. Phys. 227, 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3,375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x.Comment: 30 pages, 14 figure

An optimized intermolecular force field for hydrogen bonded organic molecular crystals using atomic multipole electrostatics

We present a re-parameterization of the a popular intermolecular force field for describing intermolecular interactions in the organic solid state. Specifically, we optimize the performance of the exp-6 force field when used in conjunction with atomic multipole electrostatics. We also parameterize force fields that are optimized for use with multipoles derived from polarized molecular electron densities, to account for induction effects in molecular crystals. Parameterization is performed against a set of 186 experimentally determined, low temperature crystal structures and 53 measured sublimation enthalpies of hydrogen bonding organic molecules. The resulting force fields are tested on a validation set of 129 crystal structures and show improved reproduction of the structures and lattice energies of a range of organic molecular crystals compared to the original force field with atomic partial charge electrostatics. Unit cell dimensions of the validation set are typically reproduced to within 3% with the re-parameterized force fields. Lattice energies, which were all included during parameterisation, are systematically underestimated when compared to measured sublimation enthalpies, with mean absolute errors of between 7.4 and 9.0%

Impact of Gd doping on morphology and superconductivity of NbN sputtered thin films

We report effect of Gd inclusion in the NbN superconductor thin films. The films are deposited on single crystalline Silicon (100) by DC reactive sputtering technique i.e., deposition of Nb and Gd in presence of reactive N2 gas. The fabricated relatively thick films (400 nm) are crystallized in cubic structure. These films are characterized for their morphology, elemental analysis and roughness by Scanning Electron Microscopy (SEM), Energy Dispersive X-ray spectroscopy (EDAX) and Atomic Force Microscopy (AFM) respectively. The optimized film (maximum Tc) is achieved with gas ratio of Ar:N2 (80:20) for both pristine and Gd doped films. The optimized NbN film possesses Tc (R=0) in zero and 140kOe fields are at 14.8K and 8.8K respectively. The Gd doped NbN film showed Tc (R=0) in zero and 130kOe fields at 11.2K and 6.8 K respectively. The upper critical field Hc2(0) of the studied superconducting films is calculated from the magneto-transport [R(T)H] measurements using GL equations. It is found that Gd doping deteriorated the superconducting performance of NbN.Comment: 14 pages Text+Figs: comments/suggestions ([email protected])/www.freewebs.com/vpsawana

A sobering assessment of small-molecule force field methods for low energy conformer predictions

We have carried out a large scale computational investigation to assess the utility of common small-molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn–Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening