914 research outputs found

    Symplectic quaternion scheme for biophysical molecular dynamics

    Get PDF
    Massively parallel biophysical molecular dynamics simulations, coupled with efficient methods, promise to open biologically significant time scales for study. In order to promote efficient fine-grained parallel algorithms with low communication overhead, the fast degrees of freedom in these complex systems can be divided into sets of rigid bodies. Here, a novel Hamiltonian form of a minimal, nonsingular representation of rigid body rotations, the unit quaternion, is derived, and a corresponding reversible, symplectic integrator is presented. The novel technique performs very well on both model and biophysical problems in accord with a formal theoretical analysis given within, which gives an explicit condition for an integrator to possess a conserved quantity, an explicit expression for the conserved quantity of a symplectic integrator, the latter following and in accord with Calvo and Sanz-Sarna, Numerical Hamiltonian Problems (1994), and extension of the explicit expression to general systems with a flat phase space

    Anisotropy of Mechanical Properties of Pinctada margaritifera Mollusk Shell

    Get PDF
    The research was co-financed by the European Union from the resources of the European Social Fund (Project No.WND-POWR.03.02.00-00-I043/16), the DAAD program. A.C. acknowledges project CGL2017-85118-P of the Spanish Ministerio de Ciencia e InnovaciĂłnfor funding.The mechanical properties such as compressive strength and nanohardness were investigated for Pinctada margaritifera mollusk shells. The compressive strength was evaluated through a uniaxial static compression test performed along the load directions parallel and perpendicular to the shell axis, respectively, while the hardness and Young modulus were measured using nanoindentation. In order to observe the crack propagation, for the first time for such material, the in-situ X-ray microscopy (nano-XCT) imaging (together with 3D reconstruction based on the acquired images) during the indentation tests was performed. The results were compared with these obtained during the micro-indentation test done with the help of conventional Vickers indenter and subsequent scanning electron microscopy observations. The results revealed that the cracks formed during the indentation start to propagate in the calcite prism until they reach a ductile organic matrix where most of them are stopped. The obtained results confirm a strong anisotropy of both crack propagation and the mechanical strength caused by the formation of the prismatic structure in the outer layer of P. margaritifera shell.The research was co-financed by the European Union from resources of the European Social Fund (Project No.WND-POWR.03.02.00-00-I043/16)

    First-principles molecular-dynamics simulations of a hydrous silica melt: Structural properties and hydrogen diffusion mechanism

    Full text link
    We use {\it ab initio} molecular dynamics simulations to study a sample of liquid silica containing 3.84 wt.% H2_2O.We find that, for temperatures of 3000 K and 3500 K,water is almost exclusively dissolved as hydroxyl groups, the silica network is partially broken and static and dynamical properties of the silica network change considerably upon the addition of water.Water molecules or free O-H groups occur only at the highest temperature but are not stable and disintegrate rapidly.Structural properties of this system are compared to those of pure silica and sodium tetrasilicate melts at equivalent temperatures. These comparisons confirm the picture of a partially broken tetrahedral network in the hydrous liquid and suggest that the structure of the matrix is as much changed by the addition of water than it is by the addition of the same amount (in mole %) of sodium oxide. On larger length scales, correlations are qualitatively similar but seem to be more pronounced in the hydrous silica liquid. Finally, we study the diffusion mechanisms of the hydrogen atoms in the melt. It turns out that HOSi2_2 triclusters and SiO dangling bonds play a decisive role as intermediate states for the hydrogen diffusion.Comment: 25 pages, 18 figures. submitte

    Molecular dynamics in shape space and femtosecond vibrational spectroscopy of metal clusters

    Full text link
    We introduce a method of molecular dynamics in shape space aimed at metal clusters. The ionic degrees of freedom are described via a dynamically deformable jellium with inertia parameters derived from an incompressible, irrotational flow. The shell correction method is used to calculate the electronic potential energy surface underlying the dynamics. Our finite temperature simulations of Ag_14 and its ions, following the negative to neutral to positive scheme, demonstrate the potential of pump and probe ultrashort laser pulses as a spectroscopy of cluster shape vibrations.Comment: Latex/Revtex, 4 pages with 3 Postscript figure

    Controlling the energy flow in nonlinear lattices: a model for a thermal rectifier

    Full text link
    We address the problem of heat conduction in 1-D nonlinear chains; we show that, acting on the parameter which controls the strength of the on site potential inside a segment of the chain, we induce a transition from conducting to insulating behavior in the whole system. Quite remarkably, the same transition can be observed by increasing the temperatures of the thermal baths at both ends of the chain by the same amount. The control of heat conduction by nonlinearity opens the possibility to propose new devices such as a thermal rectifier.Comment: 4 pages with figures included. Phys. Rev. Lett., to be published (Ref. [10] corrected
    • …
    corecore