NVT和minimize问题
NVT和minimize
1、NVT弛豫之前需不需要minimize
当然需要!! minimization的主要作用是消除原子间的 overlap, 不做其他的事情.
2、用NVT弛豫是不是就是fix NVT然后run一段时间,直到势能和总能在某一个值上下波动,即平衡?
对的.
只需要看总能量就可以了, 因为N和T是定值, 那么kinetic energy是定值, 所以势能和总能是一个概念了.
3、给定初始速度应该在哪里给定,minimize之前,还是minimize与fix NVT之间,还是别的?还有关于velocity命令的迷惑,velocity create XX seed(自己随便写还是用别的方法产生)temp temp_compute(这个计算的温度和前面的XX问题有什么区别,它对于velocity这个命令到底有什么作用)
参见 1. 回答. minimization 不需要你初速度的信息, 所以放哪里都无所谓, 你可以尝试的. 便于理解的话, 你就这样放在minize后面好了.
velocity的问题, 手册讲得很清楚了, 再去读一下吧.
4、为什么我fix NVT之后,跑了很长时间还是无法达到我指定的温度,比如
fix NVT 293.0 293.0 5 drag 0.2
显然你需要找到一个优化值, 而不是就用0.2, 参见这个帖子:
http://muchong.com/bbs/viewthread.php?tid=2435901&fpage=1
虽然fix NVT 用的是Nose-Hoover 和 berendsen thermostat 实现方法上有些不同, 但是基本的道理还是一样的.
NVT和minimize问题相关推荐
- 关于minimize跟驰豫的问题
minimize跟驰豫 不要把minimize跟驰豫同等概念而言. 大家在用LAMMPS时minimize的命令都会用到的, 一般是用来最优化初始构型, 保证初始构型的势能最优化. LAMMPS中mi ...
- lammps问题-关于NVT与NVE系综产生的不同
相同的材料构型,在不同的系综下产生如下不同的结果: 1.NVE系综下: 2.NVT系综如下: 产生的不同结果如上,原脚本如下: 1. shell mkdir initial-comp shell mk ...
- GROMACS运行参数之nvt.mdp文件详解
GROMACS(5.1.4)教程:蛋白质配体复合物 官网:点击打开链接 李老师博客:点击打开链接 蛋白质配体复合物模拟nvt平衡过程中需要用到输入文件nvt.mdp,现对里面的各种编辑项目做简单注释. ...
- scipy minimize当目标函数需要参数、当约束条件特别特别多时
from scipy.optimize import minimize 求解 官方说明文档 简单开一下开头 官方下面有例子会让人容易明白 注意我们的目的是:为了求出让fun函数最小的最优解x 当遇见目 ...
- python网页优化公司_使用python优化scipy.optimize.minimize公司
我将逐行检查您的代码,并强调一些问题:from scipy.optimize import minimize import numpy as np prices=np.array([[1.5,50,3 ...
- python实例 优化目标函数_Scipy优化算法--scipy.optimize.fmin_tnc()/minimize()
scipy中的optimize子包中提供了常用的最优化算法函数实现,我们可以直接调用这些函数完成我们的优化问题. scipy.optimize包提供了几种常用的优化算法. 该模块包含以下几个方面 使用 ...
- lammps计算聚合物例子_LAMMPS中的系综(NPT/NVT)命令
官方手册的翻译及补充,关于系综命令可见: fix NVT/NPT/NPHlammps.sandia.gov 补充部分在段落结尾以"添注"形式补充. 语法: fix ID grou ...
- tf.train.GradientDescentOptimizer().minimize()细致探索
1.tf.train.GradientDescentOptimizer().minimize() 其中的minimize可以拆为以下两个步骤: ① 梯度计算 ② 将计算出来的梯度应用到变量的更新中 拆 ...
- tensorflow中optimizer minimize自动训练简介和选择训练variable的方法
本文主要介绍tensorflow的自动训练的相关细节,并把自动训练和基础公式结合起来.如有不足,还请指教. 写这个的初衷:有些教程说的比较模糊,没体现出用意和特性或应用场景. 面向对象:稍微了解点代码 ...
最新文章
- 基于 SpringBoot 的仿豆瓣平台【源码分享】
- 伏威谈淘宝网的高并发处理与压力测试(转)
- 送一个全新的AirPods Pro
- VTK:PolyData之ExtractSelectedIds
- 一步一步实现自己的模拟控件(6)——控件树及控件区域
- Grafana plugins zabbix 安装
- 深度解析数据分析、大数据工程师和数据科学家的区别
- Socket编程 涵盖代码和函数参数介绍
- MyBatis 的执行流程,学废了!
- 2020 大会员惊喜上线!大特权限时领取!
- Network 20Q--Q2 How does Google sell ad spaces?
- 接口做的好怎么形容_大连SIEMENS西门子M174接口plc模块
- 系统学习深度学习(二十三)--SqueezeNet
- 心音与心电信号分析之一--6.26--心音信号检测系统软件设计
- POJ - 1384 Piggy-Bank(完全背包)
- 火山PC编辑框组件详解2
- listmanager htdocs
- 论文阅读-Detecting Deepfakes with Self-Blended Images (利用自混合图像检测深度伪造)
- matlab怎么改变程序字体大小,matlab怎么把程序字体变大
- 解决U盘弹出问题:无法停用 通用卷 设备,因为某一程序仍在访问它
热门文章
- 图的遍历之BSF广度优先算法6.2.2(网络整理)
- 2012 成都网络赛小记
- Objective-C之Category(分类)
- 东原地产首获“蓝筹地产” 千亿货值业绩增速成亮点
- Linux的基本指令(一):常用基础指令
- 1(基于SpringBoot)疫情下的在线办公签到系统-从零构建后端项目基础篇
- win10神州网信+UOS双系统安装步骤
- 基于Scrapy框架的当当网编程开发图书定向爬虫
- 基于LinkedList实现HashMap代码
- Apache Kafka实战读书笔记(推荐指数:☆☆☆☆☆)