detcas - Online di Cloud

PROGRAM:

NAMA

detcas - Interaksi Konfigurasi Determinan / Aktif Lengkap
Program Bidang Konsisten Luar Angkasa

DESKRIPSI

Program detca melakukan bidang konsisten-diri-ruang-aktif-lengkap (CASSCF)
optimasi orbital molekul melalui prosedur dua langkah dalam hubungannya dengan
program interaksi konfigurasi penentu DETCI. Program ini cukup sederhana dan
saat ini menggunakan pendekatan Newton-Raphson untuk memperbarui orbital, menggunakan pendekatan sederhana
perkiraan orbital Hessian. Konvergensi dipercepat menggunakan inversi langsung Pulay dari
prosedur iterative subspace (DIIS). Kode telah ditulis untuk memungkinkan lebih umum
fungsi gelombang yang tidak selalu menampilkan perlakuan CI penuh dari ruang aktif. Di dalam
khususnya, setiap fungsi gelombang CI ruang aktif terbatas (RAS) yang didukung oleh DETCI dapat:
digunakan, memungkinkan fungsi gelombang MCSCF dari tipe RASSCF.

REFERENSI

Perkiraan Orbital Hessian:

1. G. Chaban, MW Schmidt, dan MS Gordon, Theor. Chim. Akta 97, 88-95 (1997).

CI Ruang Aktif Terbatas:

1. Algoritma Interaksi Konfigurasi Berbasis Determinan untuk Lengkap dan Terbatas
Konfigurasi Ruang Interaksi, J. Olsen, BO Roos, P. Jorgensen, dan HJ Aa.
Jensen, J.Chem. fisik. 89, 2185 (1988).

SCF Ruang Aktif Terbatas:

1. P.-A. Malmqvist, A. Rendell, dan BO Roos, J. Phys. Kimia 94, 5477 (1990).

Program DETCI:

1. CD Sherrill, Algoritma Komputasi untuk Skala Besar Penuh dan Multi-Referensi
Fungsi Gelombang Interaksi Konfigurasi, tesis PhD, University of Georgia, Athens,
GA, 1996.

2. CD Sherrill dan HF Schaefer, Metode Interaksi Konfigurasi: Lanjutan
dalam Pendekatan Sangat Berkorelasi, Adv. Kimia kuantum. 34, 143-269 (1999).

FILE DIBUTUHKAN

input.dat - File masukan
file78 - integral satu elektron MO (operator fzc)
file72 - integral dua elektron MO
file73 - MO matriks kepadatan satu partikel
file74 - matriks kepadatan dua partikel MO
file75 - MO Lagrangian

FILE UPDATED

output.dat - File keluaran
file14.dat - Catatan energi dan gradien orbital

INPUT FORMAT

Argumen baris perintah berikut tersedia:

-diam Ini memberikan hasil yang sama dengan CETAK = 0.

-o nama-nama
Memberikan nama file untuk file output. Default ke output.dat.

Masukan tambahan untuk program ini dibaca dari file input.dat. Lebih umum digunakan
kata kunci adalah:

WFN = tali
Jenis fungsi gelombang yang didukung adalah CASSCF dan RASSCF.

KONVERGENSI = bilangan bulat
Konvergensi yang diinginkan pada gradien orbital. Konvergensi dicapai ketika RMS
kesalahan dalam gradien orbital kurang dari 10**(-n). Standarnya adalah 4 untuk
perhitungan energi dan 7 untuk gradien.

DOKTER = bilangan bulat_array
Vektor ini memberikan jumlah orbital yang terisi ganda di setiap irep. Ada
tidak ada default.

SOCC = bilangan bulat_array
Vektor ini memberikan jumlah orbital yang terisi tunggal di setiap irep. Ada
tidak ada default.

ENERGI_CONVERGENCE = bilangan bulat
Konvergensi yang diinginkan pada total energi MCSCF. Standarnya adalah 7.

BEKU_DOCC = bilangan bulat_array
Jumlah orbital yang menempati dua kali energi terendah di setiap orbital tak tereduksi
representasi yang secara harfiah akan dibekukan (tidak diperbarui di MCSCF). NS
Urutan kapas dari representasi yang tidak dapat direduksi digunakan. Standarnya adalah
vektor nol.

BEKU_UOCC = bilangan bulat_vektor
Jumlah orbital kosong berenergi tertinggi dalam setiap representasi tak tereduksi
yang secara harfiah akan dibekukan (tidak diperbarui di MCSCF). Standarnya adalah nol
vektor.

DIBATASI_DOCC = bilangan bulat_array
Jumlah orbital yang menempati dua kali energi terendah di setiap orbital tak tereduksi
representasi yang akan dioptimalkan tetapi tetap ditempati ganda di MCSCF.
Orbital ini muncul setelah BEKU_DOCC orbital. Standarnya adalah nol
vektor.

DIBATASI_UOCC = bilangan bulat_array
Jumlah orbital kosong berenergi tertinggi dalam setiap representasi tak tereduksi
yang akan dioptimalkan tetapi tetap kosong di MCSCF. Orbital ini datang
sebelum BEKU_UOCC orbital. Standarnya adalah vektor nol.

NCASITER = bilangan bulat
Jumlah maksimum iterasi untuk mengoptimalkan orbital. Pilihan ini harus be
ditentukan in itu DEFAULT bagian of memasukkan, karena harus terlihat oleh
program pengendalian PSI. Default ke 1.

PRINT = bilangan bulat
Opsi ini menentukan verbositas output. Nilai 1 atau 2 menentukan
pencetakan minimal, nilai 3 menentukan pencetakan verbose. Nilai 4 atau 5 adalah
digunakan untuk debugging. Jangan gunakan level 5 kecuali kasus uji sangat kecil (mis
STO H2O CISD).

Kata kunci yang kurang umum digunakan adalah:

DIIS_FREQ = bilangan bulat
Prosedur ekstrapolasi DIIS akan dicoba setiap n iterasi. Standarnya
adalah 1.

DIIS_MAX_VECS = bilangan bulat
Jumlah maksimum vektor subruang untuk prosedur DIIS. Setelah nomor ini
vektor tercapai, vektor yang lebih tua akan dijatuhkan dari subruang sebagai
diperlukan. Standarnya adalah 8.

DIIS_MIN_VECS = bilangan bulat
Jumlah minimum vektor subruang sebelum interpolasi DIIS adalah
dilakukan. Standarnya adalah 2. Tidak masuk akal untuk memiliki nilai kurang dari 2.

DIIS_START = bilangan bulat
Nomor iterasi saat Pulay's Direct Inversion of the Iterative Subspace (DIIS)
prosedur untuk percepatan konvergensi harus diaktifkan. Standarnya adalah 3.
Sebelum iterasi ini, vektor tidak ditambahkan ke subruang DIIS. Pertama
Langkah DIIS tidak akan diambil sampai DIIS_MIN_VECS vektor berada di subruang DIIS.

LAG_FILE = bilangan bulat
File (nomor unit) untuk membaca matriks lagrangian. Nilai defaultnya adalah
saat ini 75.

OPDM_FILE = bilangan bulat
File (nomor unit) untuk membaca matriks kepadatan satu partikel. Nilai default
saat ini 73.

SKALA_GRAD = boolean
Memberitahu apakah akan menskalakan gradien orbital dengan perkiraan (diagonal) orbital
Goni. Standarnya adalah BENAR.

SKALA_LANGKAH = nyata
Faktor skala untuk langkah rotasi orbital. Standarnya adalah 1.0.

TPDM_FILE = bilangan bulat
File (nomor unit) untuk membaca matriks kepadatan dua partikel. Nilai default
saat ini 74.

TINGKAT_SHIFT = boolean
Memberitahu apakah akan memungkinkan untuk pergeseran level dari matriks hessian. Ini bisa digunakan
untuk memastikan bahwa hessian adalah pasti positif untuk iterasi awal. Jika
hessian tidak pasti positif mendekati konvergensi, maka ini mungkin merupakan tanda
ketidakstabilan numerik di MCSCF. Standarnya adalah BENAR.

SHIFT = nyata
Memberitahu berapa banyak elemen diagonal goni harus digeser untuk level
bergeser. Standarnya adalah 0.01.

DETERM_MIN = nyata
Nilai minimum yang diizinkan untuk determinan matriks hessian jika level shifting
aktif. Ini tidak digunakan jika LEVEL_SHIFT=FALSE. Standarnya adalah 0.00001.

LANGKAH_MAX = nyata
Ini adalah rotasi orbital tunggal maksimum yang diizinkan. Standarnya adalah 0.30.

KEKUATAN_LANGKAH = boolean
Ini memungkinkan pengguna untuk mengabaikan langkah yang dihitung dan memaksakan langkah dalam a
arah tertentu. Ini bisa berguna jika mencoba memaksakan perhitungan
dari titik pelana. Nilai defaultnya salah.

FORCE_PAIR = bilangan bulat
Ini adalah indeks untuk pasangan independen yang akan diputar jika
FORCE_STEP=BENAR. Ini diabaikan sebaliknya. Tidak ada default.

FORCE_VALUE = nyata
Ini adalah nilai rotasi orbital untuk memutar orbital yang ditentukan oleh
FORCE_PAIR. Ini hanya digunakan jika FORCE_STEP disetel ke TRUE. Nilai defaultnya adalah
0.0.

CHECK_HESSIAN = boolean
Hitung nilai eigen dari orbital hessian. Ini bagus untuk memeriksa
solusi titik pelana. Harus digunakan tanpa pergeseran level. Standarnya adalah
SALAH.

EIGEN_VECTORS = boolean
Apakah Anda ingin menghitung vektor eigen dari orbital goni juga. Ini adalah
diabaikan jika CHECK_HESSIAN=FALSE. Standarnya adalah SALAH.

8 May, 1998 detca(1)

Gunakan detcas online menggunakan layanan onworks.net