lammps npt系综iso和aniso控压对比

article/2025/9/21 20:25:10

大家好,我是小马老师。
本文讲解lammps基础命令fix npt设置。
关于系综的设置,前面已经介绍过多次,但在咨询中发现,对于大部分初学者来说,对于系综选择和设置还是有些困惑。
fix npt命令名称就可以看出,它的主要作用是对体系进行控温和控压操作。
npt的控温原理与nvt类似,通过调整原子的速度来调整整个体系的温度。
体系压力的控制则通过调整box的尺寸进行调节,如果体系压力超过设定压力值,会扩大box尺寸,同理,体系压力小于设定压力值,则会缩小box的尺寸。
至于调节哪个方向的尺寸,可以根据需要进行设置。
如只调节z方向的压力,表示z方向的尺寸可随体系压力变化而变化,而x和y方向尺寸则固定不变,对应的命令为:

fix 1 npt temp 300 300 0.1 z 0 0 1

如果三个方向都需要调节压力,可设置iso或aniso:

fix 1 npt temp 300 300 0.1 iso 0 0 1
fix 1 npt temp 300 300 0.1 ansio 0 0 1

这两个命令都是调节三个方向的压力,不同点在哪里?
iso表示x、y、z三个方向压力耦合到一起,三个方向以相同的比例进行缩放。
aniso表示x、y、z三个方向根据各自方向的压力独立调节尺寸,压力大的方向尺寸变化大,三个方向尺寸变化可能会不相同。
下面以sio2非晶建模为例,对比iso和aniso不同控压方式的区别:
iso控压:
iso控压
aniso控压:
在这里插入图片描述

通过对比,很明显的看出,在同一方向最终的结构尺寸是不同的。
在lammps模拟时,如果需要各方向尺寸维持一定的比例,可以用iso控压,如果允许各方向尺寸自由变化,可以用aniso控压。
本公众号号正式推出一对一咨询辅导服务,根据课题方向不同详细讲解对应的in文件编写方法、模拟关键技术、数据后处理,经一对一咨询辅导后能够独立编写出in文件。

具体内容可查看公众号一对一培训说明:lammps一对一培训

公众号案例代码以及学习交流,请加QQ群:236862503
在这里插入图片描述

扫描关注微信公众号:lammps加油站。
在这里插入图片描述


http://chatgpt.dhexx.cn/article/CokNUUem.shtml

相关文章

lammps案例:npt系综练手,H2压缩模拟

上一篇文章中介绍了npt系综的设置方法,本文给出一个lammps模拟氢气压缩的案例,熟悉npt系综的设置方法。 模拟温度300k,压强100GPa,在该条件下,氢气应该为固态。 在MS中建立H2模型,共包含1000个H原子&…

lammps教程:实例讲解npt、nvt系综的选择

对于lammps新手来说,系综的选择令人困惑,到底是应该用nvt系综还是npt系综? 我在初学lammps时就被这个问题困扰了很久。如果去问别人,可能会告诉你一个答案:根据模拟需要确定。 其实和没说一样。 到底如何确定系综&a…

lammps问题-关于NVT与NVE系综产生的不同

相同的材料构型,在不同的系综下产生如下不同的结果: 1.NVE系综下: 2.NVT系综如下: 产生的不同结果如上,原脚本如下: 1. shell mkdir initial-comp shell mkdir config-comp shell mkdir restart-comp #--…

Ensemble Average(系综平均)(集平均)

什么是系综平均?: 下图(来自知乎:卖艺的小青年)中集平均实际就是系综平均,按照随机信号分析的知识来分析,我们可以把这一系列函数看成一个随机信号,取定一个时间,我们就得…

系综理论(Ensemble Theory)

本文参考[Pathria] Statistical Mechanics (3rd Ed.)里的思路整理了系综理论的一些基本知识。Pathria的书有较多繁复的计算,在此不详细给出,只希望对系综理论能形成一个宏观的认知。 1. 前言 给定一个系统的宏观态描述 ( N , V , E ) (N,V ,E) (N,V,E)…

经典系综理论

系综理论是统计力学的理论基础 本质上来说,统计热力学中只有一个问题,即给定能量$E$,如何分布在$N$个全同系统构成的系综上 --薛定谔 1.系综的定义 在下面的介绍中我们可以逐渐体会这句话的含义,此时我们暂且把这个定位我们解决的…

lammps教程:nve/nvt/npt系综设置方法

系综是指大量性质和结构完全相同的、处于各种运动状态的、各自独立的系统的集合。 在lammps模拟中,常用的系综有nve、nvt、npt三种系综,下面说说这三种系综的使用方法。 1. nve(微正则系综) nve系综保证系统中原子数量n、体系总…

数值分析(2)-多项式插值: 拉格朗日插值法

整理一下数值分析的笔记~ 目录: 1. 误差 2. 多项式插值与样条插值(THIS) 3. 函数逼近 4. 数值积分与数值微分 5. 线性方程组的直接解法 6. 线性方程组的迭代解法 7. 非线性方程求根 8. 特征值和特征向量的计算 9. 常微分方程初值问题的数值解 1. 拉格朗日基函数 定…

拉格朗日插值法求多项式系数 (附代码)

写在前面: 学了拉格朗日插值法之后发现大家都说可以在O(n^2)时间内得到多项式系数,但是没有找到代码,网上找了很多资料又因为我太弱了没能看懂,最后在emofunx学长的帮助下终于搞明白了。 由于太弱没能看懂的文章 引入 我们都知…

拉格朗日插值多项式的龙格现象

利用插值基函数可以很容易得到拉格朗日插值多项式,公式结构紧凑,在理论分析中甚为重要。但存在以下问题: 当插值节点增减时,计算要全部重新进行,甚为不便(当然现在我们有MATLAB,只需改变下参数…

拉格朗日多项式插值法

在数值分析中,拉格朗日常用于多项式插值。假定提供一组数据点[xi,yi],拉格朗日插值多项式就是由这些 数据的线性运算得到的。 其中基本的多项式有以下公式计算得到 注意 1.第一提供的xi应该是没有相同的,否则不能应用此算法 2.对于每一个xi&…

拉格朗日插值和求多项式系数

拉格朗日介绍 先说说拉格朗日是啥吧 首先 拉格朗日插值是给你 n1 个点 (x,y) 然后根据这n个点可以O(n^2)的求出多项式的系数。也就是解出这个多项式的答案。 假设给你一个多项式 ya0a1*xa2*x^2 然后给你3个解 (x1,y1)(x2,y2)(x3,y3)你第一个想法是怎么解?解方程啊…

数学建模准备 插值(拉格朗日多项式插值,牛顿多项式插值,分段线性插值,分段三次样条插值,分段三次Hermite插值)

文章目录 摘要(必看)0 基础概念什么是插值插值用途什么是拟合插值和拟合的相同点插值和拟合的不同点 1 常用的基本插值方法1.1 多项式插值法1.1.1 拉格朗日多项式插值法多项式插值并不是次数越大越好(龙格现象)分段低次线性插值以…

数模--拉格朗日多项式插值、matlab实现

拉格朗日多项式公式: matlab中插入一个值的代码 function yhlagrange(x,y,xh) %定义拉格朗日插入函数 nlength(x); %统计x和xh的长度 mlength(xh); yhzeros(1,m); %构建一行m列的zero矩阵 c1ones(n-1,1); …

计算方法实验(一):拉格朗日插值多项式

拉格朗日插值数学原理 给定平面上 n 1 n 1 n1个不同的数据点 ( x k , f ( x k ) ) (x_{k},f(x_{k})) (xk​,f(xk​)), k 0 , 1 , ⋯ , n k 0,1,\cdots,n k0,1,⋯,n, x i ≠ x j x_{i} \neq x_{j} xi​​xj​, i ≠ j i \neq j i​j&…

拉格朗日插值多项式在MATLAB中的实现

拉格朗日插值多项式在MATLAB中的实现 Hi! 这是我的一个CSDN博客 例 以下给出了针对题(2)的使用方法 输入: 1.节点值 2.需要插值的原函数 对于题(2)需要设置参数为: xx [0,0.25,0.5,1];% The nodes(n…

lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题)

lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题) 关于拉格朗日插值法相关理论知识,在这里小编不在赘述,请不明白的小伙伴自行百度。小编只负责给出matlab源码。 **例题:**看下面例题(如图): matlab代码: %%%% 求拉格朗日多项式及基函…

拉格朗日插值算法分析

这几天一直研究拉格朗日多项式,今天将自己对拉格朗日多项式的理解写在这里,方便大家交流。 在数值分析中,拉格朗日常用于多项式插值。假定提供一组数据点[xi,yi],拉格朗日插值多项式就是由这些 数据的线性运算得到的。 其中基本…

拉格朗日(Lagrange)插值多项式

题目: 基本原理: 拉格朗日(Lagrange)插值多项式python实现: # encoding: utf-8 from symtable import Symbol X[0.4,0.5,0.6,0.7,0.8] Y[-0.9163,-0.6931,-0.5108,-0.3567,-0.2231] print(X,Y) Lfloat(0.0) x0.54 len…

使用拉格朗日多项式(Lagrangian polynomials)的插值法(python,数值积分)

第三十五篇 拉格朗日多项式插值 插值多项式 首先考虑一个函数的推导,该函数精准地通过一系列np离散数据点。虽然有无限多的函数具备这个条件,但我们将专注于最简单的一个,一个n阶多项式,其中n np−1。我们称这个函数为“插值多…