DSRG-MRPT2 PySCF implementation

.github/workflows/ci_conda.yml

@@ -0,0 +1,25 @@
+name: CI_conda
+
+on:
+  workflow_dispatch:
+
+jobs:
+  build-conda-linux:
+    runs-on: ubuntu-latest
+    strategy:
+      fail-fast: false
+    steps:
+      - uses: actions/checkout@v4
+      - name: Setup conda
+        uses: s-weigand/setup-conda@v1
+        with:
+          update-conda: true
+          conda-channels: conda-forge
+      - run: conda --version
+      - run: which python
+      - name: Build conda package
+        run: |
+          export CMAKE_BUILD_PARALLEL_LEVEL=4
+          conda install -y conda-build
+          conda config --set anaconda_upload False
+          conda build --output-folder . conda

.github/workflows/publish.yml

@@ -39,7 +39,7 @@ jobs:
     strategy:
       fail-fast: false
     env:
-      img: quay.io/pypa/manylinux2014_aarch64:2023-03-12-25fd859
+      img: quay.io/pypa/manylinux2014_aarch64:latest
     steps:
     - name: Checkout
       uses: actions/checkout@v4

CHANGELOG

@@ -1,3 +1,41 @@
+PySCF 2.7.0 (2024-09-23)
+------------------------
+* Added
+  - Superposition of Atomic Potentials (SAP) initial guess for SCF methods.
+  - Supports pickle serialization for all methods.
+  - ADC 1-particle density matrix and dipole moment.
+  - Spin-separated 3-RDMs.
+  - Traceless quadrupole moment for SCF methods.
+  - Supports for fractional coordinates in Cell.
+  - Population analysis for KSCF.
+  - A, B matrices for k-point TDRKS.
+  - AutoAux scheme for generating auxiliary basis sets.
+* Improved
+  - Automatic cleanup for HDF5 temporary files.
+  - Saves CI coefficients for SA-CASSCF to chkfile.
+  - UHF/UKS initial guess with better spin-symmetry breaking code.
+  - New attribute .cycles in SCF methods and CC methods to save iteration counts.
+  - FFT performance.
+  - CPHF convergence in nuclear hessian.
+  - Eigenvalue solver for linear response theory.
+  - Enhance TDDFT code with symmetry adaptation in diagonalization.
+  - More efficient PBC-CDERI loading.
+* Fixes
+  - The missing atomic radius parameters for the Treutler-Ahlrichs scheme in DFT grids
+  - Symmetry adapted GHF with complex orbitals.
+  - Complex-valued FCI matvec operation.
+  - Bug in Smearing-RHF for odd number of electrons and disable Smearing-ROHF.
+  - MP2 FNO bug when pct_occ=1 .
+  - DHF dipole moment.
+  - Bug in PBC ECP integral.
+  - Bug in parser for spin-orbit ECPs.
+  - Wrong transition dipoles of triplet TDDFT states.
+  - GHF-X2C dipole moment.
+  - FCI symmetry validation code for cylindrical symmetry.
+  - Missing complex conjugation for complex orbitals in SymAdaptedGHF.
+  - Custom UHF objects using the Hamiltonian in FCIDUMP.
+
+
 PySCF 2.6.2 (2024-06-19)
 ------------------------
 * Fixes

NOTICE

@@ -110,6 +110,8 @@ Christopher Hillenbrand
 scohenjanes5
 Michal Krompiec
 Victor Yu
+James Serna (The Ohio State University)
+Brian Zhao
 
 
 ---

README.md

@@ -7,13 +7,14 @@ Python-based Simulations of Chemistry Framework
 [![Build Status](https://github.com/pyscf/pyscf/workflows/CI/badge.svg)](https://github.com/pyscf/pyscf/actions?query=workflow%3ACI)
 [![codecov](https://codecov.io/gh/pyscf/pyscf/branch/master/graph/badge.svg)](https://codecov.io/gh/pyscf/pyscf)
 
-2024-06-19
+2024-09-23
 
-* [Stable release 2.6.2](https://github.com/pyscf/pyscf/releases/tag/v2.6.2)
+* [Stable release 2.7.0](https://github.com/pyscf/pyscf/releases/tag/v2.7.0)
 * [Changelog](../master/CHANGELOG)
 * [Documentation](http://www.pyscf.org)
 * [Installation](#installation)
 * [Features](../master/FEATURES)
+* [News](https://pyscf.org/news.html): **2nd PySCF Developers Meeting!**
 
 
 # Installation

examples/df/01-auxbasis.py

@@ -65,3 +65,8 @@
 mf.with_df.auxbasis = df.aug_etb(mol, beta=1.7)
 mf.kernel()
 
+# Generating auxiliary basis using the AutoAux algorithms proposed by Stoychev
+# (JCTC, 13, 554)
+mf = scf.RHF(mol).density_fit()
+mf.with_df.auxbasis = df.auto_aux(mol)
+mf.kernel()

examples/gto/04-input_basis.py

@@ -65,7 +65,7 @@
 mol = gto.M(
     atom = '''O 0 0 0; H1 0 1 0; H2 0 0 1''',
     basis = {'O': gto.parse('''
-# Parse NWChem format basis string (see https://bse.pnl.gov/bse/portal).
+# Parse NWChem format basis string (see https://www.basissetexchange.org/).
 # Comment lines are ignored
 #BASIS SET: (6s,3p) -> [2s,1p]
 O    S

examples/gto/31-basis_set_exchange.py

@@ -14,6 +14,22 @@
 a local copy of the Basis Set Exchange.
 '''
 
+# PySCF has a native interface to the Basis Set Exchange, which is
+# used as a fall-through if the basis set was not found within PySCF
+# itself
+mol = gto.M(
+    atom = '''
+N          0.6683566134        0.2004327755        0.0000000000
+H          0.9668193796       -0.3441960976        0.8071193402
+H          0.9668193796       -0.3441960976       -0.8071193402
+F         -0.7347916126       -0.0467759204        0.0000000000
+''',
+    basis = 'HGBSP1-7',
+    verbose = 4
+)
+
+
+# One can also use the Basis Set Exchange in a more direct fashion
 mol = gto.M(
     atom = '''
 N          0.6683566134        0.2004327755        0.0000000000

examples/pbc/01-input_output_geometry.py

@@ -0,0 +1,48 @@
+#!/usr/bin/env python
+
+'''
+This example demonstrates some ways to input the Cell geometry, including
+lattice vectors, and to write the Cell geometry to file.
+
+Example solid is wurtzite BN.
+'''
+
+from pyscf.pbc import gto
+
+# Input Cartesian coordinates for the lattice vectors a and the atomic
+# positions. Coordinates are in Angstrom by default
+
+cell = gto.Cell()
+cell.a = [[  2.5539395809,  0.0000000000,  0.0000000000],
+          [ -1.2769697905,  2.2117765568,  0.0000000000],
+          [  0.0000000000,  0.0000000000,  4.2268548012]]
+cell.atom = '''
+            B     1.276969829         0.737258874         4.225688066
+            N     1.276969829         0.737258874         2.642950986
+            B     0.000000000         1.474517748         2.112260792
+            N     0.000000000         1.474517748         0.529523459
+            '''
+cell.pseudo = 'gth-pade'
+cell.basis = 'gth-szv'
+cell.build()
+
+# Write cell geometry to file
+# - Format is guessed from filename
+# - These can be read by VESTA, Avogadro, etc.
+# - XYZ is Extended XYZ file format, which includes lattice vectors in the
+#   comment line
+cell.tofile('bn.vasp')
+cell.tofile('POSCAR')
+cell.tofile('bn.xyz')
+
+# Read a and atom from file
+from pyscf.pbc.gto.cell import fromfile
+a, atom = fromfile('bn.vasp')
+a, atom = fromfile('bn.xyz')
+
+# Read a and atom from file directly into Cell
+cell = gto.Cell()
+cell.fromfile('bn.vasp')
+cell.pseudo = 'gth-pade'
+cell.basis = 'gth-szv'
+cell.build()

examples/pbc/09-band_ase.py

@@ -4,7 +4,7 @@
 
 import matplotlib.pyplot as plt
 
-from ase.lattice import bulk
+from ase.build import bulk
 from ase.dft.kpoints import sc_special_points as special_points, get_bandpath
 
 c = bulk('C', 'diamond', a=3.5668)

examples/solvent/05-pcm.py

@@ -63,4 +63,13 @@
 td.kernel()
 
 mf = mol.RHF().PCM().run()
-td = mf.TDA().run()
+td = mf.TDA().run()
+
+# infinite epsilon with any PCM computations
+cm = pcm.PCM(mol)
+cm.eps = float('inf') # ideal conductor
+mf = dft.RKS(mol, xc='b3lyp')
+mf = mf.PCM(cm)
+mf.kernel()
+
+

examples/solvent/07-cosmors.py

@@ -0,0 +1,115 @@
+#!/usr/bin/env python
+'''
+An example of using COSMO-RS functionality.
+'''
+
+#%% Imports
+
+import io
+
+import numpy as np
+import matplotlib.pyplot as plt
+
+from pyscf import gto, dft
+from pyscf.solvent import pcm
+from pyscf.solvent.cosmors import get_sas_volume
+from pyscf.solvent.cosmors import write_cosmo_file
+from pyscf.solvent.cosmors import get_cosmors_parameters
+
+
+#%% Set parameters
+
+# get molecule
+coords = '''
+C        0.000000    0.000000             -0.542500
+O        0.000000    0.000000              0.677500
+H        0.000000    0.9353074360871938   -1.082500
+H        0.000000   -0.9353074360871938   -1.082500
+'''
+mol = gto.M(atom=coords, basis='6-31G*', verbose=1)
+
+# configure PCM
+cm = pcm.PCM(mol)
+cm.eps = float('inf') # f_epsilon = 1 is required for COSMO-RS
+cm.method = 'C-PCM' # or COSMO, IEF-PCM, SS(V)PE, see https://manual.q-chem.com/5.4/topic_pcm-em.html
+cm.lebedev_order = 29 # lebedev grids on the cavity surface, lebedev_order=29  <--> # of grids = 302
+
+
+#%% COSMO-files
+
+# run DFT SCF (any level of theory is OK, though DFT is optimal)
+mf = dft.RKS(mol, xc='b3lyp')
+mf = mf.PCM(cm)
+mf.kernel()
+
+# generate COSMO-file
+with io.StringIO() as outp: # with open('formaldehyde.cosmo', 'w') as outf:
+    write_cosmo_file(outp, mf)
+    print(outp.getvalue())
+
+
+# if PCM DFT were computed with fepsilon < 1 <=> eps != inf the ValueError will be raised
+# use ignore_low_feps=True to overrule it, but please be sure you know what you're doing
+
+# run DFT SCF
+cm = pcm.PCM(mol)
+cm.eps = 32.613 # methanol
+mf = dft.RKS(mol, xc='b3lyp')
+mf = mf.PCM(cm)
+mf.kernel()
+
+# try to get COSMO-file
+with io.StringIO() as outp:
+    try:
+        write_cosmo_file(outp, mf)
+        print(outp.getvalue())
+    except ValueError as e:
+        print(e)
+
+# overruling
+with io.StringIO() as outp:
+    write_cosmo_file(outp, mf, ignore_low_feps=True)
+    print(outp.getvalue())
+
+
+# The molecular volume is computed for the solvent-accessible surface generated
+# within pcm.PCM object. Please note that r_sol is assumed to be equal to zero,
+# so that SAS is a union of vdW spheres.
+# Increasing integration step increases accuracy of integration, but for most COSMO-RS
+# modelling and ML step=0.2 should be enough:
+
+# compute volumes with different step values
+steps = [1.0, 0.5, 0.2, 0.1, 0.05, 0.02, 0.01]
+Vs = [get_sas_volume(mf.with_solvent.surface, step) for step in steps]
+# plot
+ax = plt.gca()
+ax.plot(Vs)
+ax.set_xticks(range(len(steps)))
+_ = ax.set_xticklabels(steps)
+plt.show()
+
+
+#%% Sigma-profiles
+
+# compute SCF PCM
+cm = pcm.PCM(mol)
+cm.eps = float('inf')
+mf = dft.RKS(mol, xc='b3lyp')
+mf = mf.PCM(cm)
+mf.kernel()
+
+# compute sigma-profile and related parameters
+params = get_cosmors_parameters(mf)
+print(params)
+plt.plot(params['Screening charge, e/A**2'],
+         params['Screening charge density, A**2'])
+plt.show()
+
+# custom sigma grid
+sigmas_grid = np.linspace(-0.025, 0.025, 101)
+params = get_cosmors_parameters(mf, sigmas_grid)
+plt.plot(params['Screening charge, e/A**2'],
+         params['Screening charge density, A**2'])
+plt.show()
+
+

examples/tools/01-fcidump.py

@@ -104,3 +104,4 @@
 mf.mol.verbose = 4
 mf.run()
 mf.MP2().run()
+mf.to_uhf().newton().run()

pyproject.toml

@@ -57,8 +57,9 @@ cppe = ["cppe"]
 pyqmc = ["pyqmc"]
 bse = ["basis-set-exchange"]
 dispersion = ["pyscf-dispersion"]
+ccpy = ["coupled-cluster-py"]
 
-all = ["pyscf[forge,geomopt,doci,properties,semiempirical,cppe,pyqmc,bse,dispersion]"]
+all = ["pyscf[forge,geomopt,doci,properties,semiempirical,cppe,pyqmc,bse,dispersion,ccpy]"]
 
 # extras which should not be installed by "all" components
 cornell_shci = ["pyscf-cornell-shci"]

pyscf/__init__.py

@@ -35,7 +35,7 @@
 
 '''
 
-__version__ = '2.6.2'
+__version__ = '2.7.0'
 
 import os
 import sys