libcint 开源项目使用教程

羿妍玫Ivan

于 2024-10-10 07:58:27 发布

阅读量649

点赞数 8

本文链接：https://blog.csdn.net/gitblog_00255/article/details/142806242

版权

libcint 开源项目使用教程

libcint general GTO integrals for quantum chemistry 项目地址: https://gitcode.com/gh_mirrors/li/libcint

1. 项目介绍

libcint 是一个开源的库，用于计算量子化学中的高斯积分。它提供了 C/Fortran API，用于评估一电子和二电子积分，支持笛卡尔、实球谐和旋量高斯型函数。libcint 的主要特点包括：

支持多种高斯型函数：包括笛卡尔高斯型函数、实球谐高斯型函数和旋量高斯型函数。
一电子和二电子积分：支持常规动能积分、核吸引积分（支持高斯核模型）和二电子积分（忽略小于 1e-15 的值）。
线程安全：提供统一的 API 用于所有类型的积分。
最小化外部依赖：仅依赖于 BLAS 库，性能差异通常小于 1%。
低内存使用：仅存储必要的中间数据，内存占用通常小于 1MB。

2. 项目快速启动

2.1 安装依赖

在开始使用 libcint 之前，需要确保系统中安装了以下依赖：

BLAS 库
Python 2.5 或更高版本（可选，用于测试）
Numpy（可选，用于测试）
clisp / SBCL（可选，用于生成新积分）

2.2 下载和编译 libcint

首先，从 GitHub 仓库下载 libcint：

git clone https://github.com/sunqm/libcint.git
cd libcint

接下来，创建并进入构建目录，然后使用 CMake 进行配置和编译：

mkdir build
cd build
cmake ..
make
sudo make install

2.3 使用示例

以下是一个简单的示例，展示如何使用 libcint 计算一电子积分：

#include <stdio.h>
#include "cint.h"

int main() {
    double buf[100];
    int atm[100];
    int bas[100];
    double env[100];

    // 初始化 atm, bas, env 数组
    // ...

    int natm = 2; // 原子数
    int nbas = 2; // 基函数数

    // 计算一电子积分
    cint1e_kin_sph(buf, atm, natm, bas, nbas, env);

    // 输出结果
    for (int i = 0; i < 100; i++) {
        printf("%f\n", buf[i]);
    }

    return 0;
}