计算化学公社

标题: 求助：电势与量子电容 [打印本页]

作者
Author: RSS 时间: 2023-10-25 21:04
标题: 求助：电势与量子电容
各位老师，请问下在量子电 (, 下载次数 Times of downloads: 11) 容与电势的关系图中，为什么我拟合出来的图所有的体系都是一样的趋势呢？能否帮我看看是代码哪里出了问题嘛？以下是我的代码：cd ('F:\桌面')
% 导入态密度数据文件
data = readtable('DOSCAR 3.dat');
energy = data(:, 1).Variables; % 能量数据
density = data(:, 2).Variables; % 态密度数据
phi = linspace(-1, 1, 100); % 生成100个从-1到1的电势值

% 定义量子电容函数
quantum_capacitance = @(x, a, b, phiG, T) a * exp(-b*(x-phiG)/(T+1e-10));

% 定义拟合模型
model = fittype(@(a, b, phiG, T, x) quantum_capacitance(x, a, b, phiG, T));
% 进行拟合
fitresult = fit(energy, density, model, 'StartPoint', [-0.5, 0.5, 0, 300]);
a = fitresult.a; % 拟合参数a
b = fitresult.b; % 拟合参数b
phiG = fitresult.phiG; % 拟合参数phiG
T = 300;

% 计算电势-量子电容
C = quantum_capacitance(phi, a, b, phiG, T);

% 绘制电势-量子电容关系图
figure;
plot(phi, C);
xlabel('phi (V)')
ylabel('C (μF/cm^2)');

作者
Author: 幻想银河 时间: 2023-12-8 13:54
请问楼主问题解决了吗我也在计算这个量子电容也遇到一些问题是否可以交流一下

作者
Author: W1312319875 时间: 2024-5-6 20:19
求楼主赐教，最近也在做量子电容的计算，方便加个联系方式学习一下吗

欢迎光临计算化学公社 (http://bbs.keinsci.com/)