python分子化学模拟_#分子模拟#MDTraj分子模拟python包(一)

python分子化学模拟_#分子模拟#MDTraj分子模拟python包(一)MDTraj是分子动力学模拟的一个python包,相对于MDAnalysis个人觉得操作性更强,更加Python范一些。其能够进行不同模拟软件的轨迹转换,常规计算,分析等等一体化。今天我们介绍其安装方法和简单使用。官方地址:点击进入安装方法推荐使用conda安装mdtraj.$condainstall-cconda-forgemdtraj同样也可以使用pip来安装mdtraj测试安装运…

大家好,又见面了,我是你们的朋友全栈君。

MDTraj是分子动力学模拟的一个python包,相对于MDAnalysis个人觉得操作性更强,更加Python范一些。其能够进行不同模拟软件的轨迹转换,常规计算,分析等等一体化。今天我们介绍其安装方法和简单使用。

19e7a9833ee1c6d324dc18c39af1aebd.png

官方地址:点击进入

安装方法

推荐使用conda安装mdtraj.

$ conda install -c conda-forge mdtraj

同样也可以使用pip来安装mdtraj

测试安装

运行测试是一个非常好的方法来认证工作,测试套件采用nose,如果还没有安装可以使用pip来进行安装

pip install nose

运行测试命令如下

nosetests mdtraj -v

例子

这段文字提供一系列的例子,资源和代码来帮助使用mdtraj

如果是通过编译安装的mdtraj可以进入path-to-mdtraj/examples查看,如果是一键安装,可以通过Github进行查看,其为Ipython notebook格式。

1. MDTraj介绍

首先从硬盘中加载轨迹,MDTraj会自动的使用最合适的方式加载不同的文件格式。

import mdtraj as md

t=md.load(‘trajectory.xtc’,top=’trajectory.pdb’)

print t

一些轨迹例如Gromacs XTC轨迹文件并不包含拓扑信息,我们需要采用top关键字来加载拓扑文件,例如PDB文件.

如果你只对部分轨迹感兴趣,你可以切割(slice)他们

#查看十帧(frames)

print t[1:10]

#查看最后一帧

print t[-1]

轨迹对象包含许多对象,最多显而易见的是卡迪儿(Cartesuab)坐标.作为numpy array存储在xyz下.轨迹中的距离单位均为纳米(nanometers).时间单位为皮秒(picoseconds).角度存储为度(不是弧度).

print t.xyz.shape

print np.mean(t.xyz)

#第一个十帧时间模拟

print t.time[0:10]

# 最后一帧的晶胞长度

t.unitcell_lengths[-1]

(100, 22, 3)

0.89365752249053032

array([ 0.002, 0.004, 0.006, 0.008, 0.01 , 0.012, 0.014, 0.016,0.018, 0.02 ], dtype=float32)

array([ 2., 2., 2.], dtype=float32)

保存轨迹到文件也非常容易操作

# 保存成hdf5格式,最后一个2表示每2帧保存一次

t[::2].save(‘halftraj.h5’)

#个人推荐的保存成dcd格式

t[0:10].save_dcd(‘first-ten-frames.dcd)

轨迹包含拓扑对象的引用,这可以派上用场。例如,如果你想保存你的轨迹一份只含有α碳原子的轨迹,你可以这样做:

atoms_to_keep=[a.index for a in t.topology.atoms if a.name == ‘CA’]

t.restrict_atoms(atoms_to_keep) #在轨迹中扮演适当的作用

t.save(‘CA-only.h5’)

下期我们将会介绍其基本的原子选择操作。

更多原创精彩内容敬请关注生信杂谈:

108da3abdd430b2f7ab159f41ec5c92e.png

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/141933.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • java中文乱码_Java中文乱码问题的解决方案[通俗易懂]

    java中文乱码_Java中文乱码问题的解决方案[通俗易懂]只要掌握了中文乱码问题产生的原因,然后对症下药,就可以顺利地解决这些问题。下面我们对容易产生乱码问题的场景进行分析,并提出解决方案。1.以POST方法提交的表单数据中有中文字符由于Web容器默认的编码方式是ISO-8859-1,在Servlet/JSP程序中,通过请求对象的getParameter()方法得到的字符串是以ISO-8859-1转换而来,这是导致乱码产生的原因之一。为了避免容器以ISO…

  • 链表经典算法

    链表经典算法

  • vagrant up报错:unable to mount VirtualBox shared folders

    vagrant up报错:unable to mount VirtualBox shared folders

  • SSM-Spring(2)_AOP[通俗易懂]

    SSM-Spring(2)_AOP[通俗易懂]AOP用Spring需要导入包<dependency> <groupId>org.aspectj</groupId> <artifactId>aspectjweaver</artifactId> <version>1.9.4</version> </dependency>方式一:使用Spring接口编写javapackage com.kuang.log;

  • NTP协议详解_ntp是安全协议吗

    NTP协议详解_ntp是安全协议吗前言NTP(NetworkTimeProtocol)网络时间协议基于UDP,用于网络时间同步的协议,使网络中的计算机时钟同步到UTC,再配合各个时区的偏移调整就能实现精准同步对时功能。提供NTP对时的服务器有很多,比如微软的NTP对时服务器,利用NTP服务器提供的对时功能,可以使我们的设备时钟系统能够正确运行。NTP报文格式NTP报文格式如上图所示,它的字段含义参考如下:L…

    2022年10月12日
  • FileInputStream的available方法

    FileInputStream的available方法available():返回与之关联的文件的字节数importjava.io.File;importjava.io.FileInputStream;importjava.io.IOException;publicclassFileInputStreamDemo2{ publicstaticvoidmain(String[]args)throwsIOE…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号