detci - Online di Cloud

PROGRAM:

NAMA

detci - Program Interaksi Konfigurasi Determinan

DESKRIPSI

Program detci mendiagonalisasi operator Hamiltonan elektronik nonrelativistik dalam a
dasar determinan Slater. Himpunan determinan yang digunakan (ruang CI) dapat dipilih dalam a
berbagai cara. Program dapat menangani ruang CI apa pun yang dapat diformulasikan sebagai
Ruang Aktif Terbatas CI. Ini termasuk CISD, CISDT, CISDTQ, dll., hingga CI Penuh, sebagai
serta CI multireferensi di mana referensi dipilih sebagai semua penentu dalam
yang hingga n elektron tereksitasi di beberapa ruang aktif MO. Ini termasuk CISD[T],
CISD[TQ], dan CI orde kedua (SOCI).

REFERENSI

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).

2. Melewati Batas Satu Miliar dalam Perhitungan Full Configuration-Interaction (FCI),
J. Olsen, P. Jorgensen, dan J. Simons, Chem. fisik. Lett. 169, 463 (1990).

Subruang virtual tersier (RAS IV):

1. Fungsi Gelombang Variasi Ringkas Menggabungkan Triple dan Quadruple Terbatas
Eksitasi, CD Sherrill dan HF Schaefer, J. Phys. Kimia 100, 6069-6075
(1996).

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.

FILE DIBUTUHKAN

input.dat - File masukan
file71 - Integral satu elektron yang diubah
file72 - Integral dua elektron yang diubah

SEMENTARA FILE BEKAS

file50 - Diagonal Hamiltonian
file51 - vektor CI
file52 - Vektor sigma
file53 - file D (vektor koreksi)

FILE UPDATED

output.dat - File keluaran

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.

-e Opsi ini menyebabkan total energi CI atau energi ditulis ke file yang disebut
detci_energies.dat.

-c nilai
Memberikan konvergensi yang lebih longgar pada vektor CI, berguna dalam perhitungan DETCAS. NS
nilai adalah bilangan real, bukan bilangan bulat seperti pada KONVERGENSI. Konvergensi yang digunakan
akan menjadi yang lebih longgar nilai dan KONVERGENSI.

Masukan tambahan untuk program ini dibaca dari file Kata kunci berikut adalah
sah:

KONVERGENSI = bilangan bulat
Konvergensi yang diinginkan pada vektor CI. Konvergensi dicapai ketika RMS dari
kesalahan dalam vektor CI kurang dari 10**(-n). Standarnya adalah 4 untuk energi dan 7
untuk gradien. Ini bukan kriteria konvergensi vektor CI yang sama seperti yang ditemukan di
GUGACI.

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.

DIAG_METHOD = string
Ini menentukan metode mana yang akan digunakan dalam mendiagonalisasi Hamiltonian. NS
opsi yang valid adalah: RSP, untuk membentuk seluruh matriks H dan mendiagonalisasi menggunakan libciomr
untuk mendapatkan semua nilai eigen (nb membutuhkan memori BESAR); OLSEN, untuk menggunakan Olsen's
metode subruang terbalik prakondisi (1990); MITRUSHENKOV, untuk menggunakan 2x2
metode Olsen/Davidson; dan DAVIDSON (Atau SEM) untuk menggunakan Ekspansi Simultan Liu
Metode, yang identik dengan metode Davidson jika hanya satu akar yang ditemukan.
Ada juga mode debugging SEM, UJI SEMT. itu SEM metode adalah yang paling
kuat, tetapi juga membutuhkan 2(N*M)+1 vektor CI pada disk, di mana N adalah maksimum
jumlah iterasi dan M adalah jumlah akar.

PREKONDISI = string
Ini menentukan jenis prasyarat yang akan digunakan dalam diagonalisasi yang dipilih
metode. Opsi yang valid adalah: DAVIDSON yang mendekati matriks Hamilton
oleh elemen diagonal; H0BLOCK_INV yang menggunakan Hamiltonian eksak dari
H0_BLOCKSIZE dan secara eksplisit membalikkannya; GEN_DAVIDSON yang melakukan spektral
dekomposisi H0BLOCK; ITER_INV menggunakan pendekatan iteratif untuk mendapatkan
vektor koreksi H0BLOCK. H0BLOCK_INV, GEN_DAVIDSON, dan ITER_INV
semua pendekatan secara formal setara tetapi ITER_INV kurang komputasi
mahal. Standarnya adalah DAVIDSON.

REFERENSI = string
Ini menentukan jenis fungsi referensi. Ini adalah RHF atau ROHF. UHF dan
TWOCON tidak didukung. Untuk ROHF, kelipatan 1 menyiratkan kulit terbuka
kaos. Program akan berjalan untuk singlet kulit terbuka, tetapi belum
disesuaikan dengan benar untuk menggunakan referensi dua penentu yang benar dalam kasus ini, jadi
berjalan dengan referensi singlet shell terbuka tidak disarankan kecuali untuk CI lengkap.

UPDATE = string
DAVIDSON menggunakan pembaruan DAVIDSON standar atau rumus vektor koreksi, sementara
OLSEN menggunakan vektor koreksi OLSEN. Standarnya adalah DAVIDSON.

HD_OTF = boolean
Jika TRUE, elemen-elemen diagonal dari matriks Hamilton dihitung secara on-the-fly,
jika tidak, vektor elemen diagonal ditulis ke file terpisah pada disk. Bawaan
adalah benar.

HD_AVE = string
HD_EXACT menggunakan energi diagonal yang tepat yang menghasilkan vektor ekspansi yang
mematahkan simetri putaran. HD_KAVE rata-rata energi diagonal selama spin-coupling
mengatur menghasilkan spin vektor ekspansi murni. ORB_ENER menggunakan jumlah orbital
pendekatan energi memberikan spin vektor ekspansi murni tetapi biasanya menggandakan
jumlah iterasi davidson. PENGINJILAN menggunakan jumlah dan perbedaan orbital
energi dengan energi referensi SCF untuk menghasilkan vektor ekspansi murni spin.
PEMIMPIN pendekatan yang mengurangi kontribusi satu elektron dari
energi orbital, dikalikan dengan 0.5, dan menambahkan kembali kontribusi satu elektron
dalam, menghasilkan vektor ekspansi murni putaran dan dikembangkan oleh Anda dengan benar dan berfungsi sebagai
juga PENGINJILAN.

NODFILE = boolean
Hanya mungkin jika NUM_ROOTS = 1. Menggunakan ruang vektor terakhir dalam file BVEC untuk
menulis DVEC awal daripada menggunakan file DVEC terpisah.

ENERGI_CONVERGENCE = bilangan bulat
Konvergensi yang diinginkan pada energi CI. Standarnya adalah 6 untuk energi titik tunggal
dan 8 untuk gradien atau CASSCF.

EX_LVL = bilangan bulat
Tingkat eksitasi untuk eksitasi ke orbital virtual (default 2, yaitu CISD).

VAL_EX_LVL = bilangan bulat
Tingkat eksitasi untuk referensi dalam orbital RAS II. Default ke nol.

CFI = boolean
Jika bendera ini disetel ke BENAR, maka penyimpanan string disederhanakan untuk Full
CI dan perhitungannya membutuhkan lebih sedikit overhead. Namun, hasil akhir harus
menjadi identik dengan mereka ketika CFI = SALAH. Dapat menyebabkan hasil yang tidak terduga jika CFI =
BENAR tapi EX_LVL tidak konsisten dengan CI Penuh.

BEKU_DOCC = bilangan bulat_array
Jumlah orbital yang menempati dua kali energi terendah di setiap orbital tak tereduksi
representasi dari mana tidak akan ada eksitasi. Pemesanan Kapas dari
representasi yang tidak dapat diubah digunakan. Standarnya adalah vektor nol.

BEKU_UOCC = bilangan bulat_vektor
Jumlah orbital kosong berenergi tertinggi dalam setiap representasi tak tereduksi
di mana tidak akan ada eksitasi. Standarnya adalah vektor nol.

BEKU_INTI = boolean
Opsi ini menentukan apakah orbital inti beku akan dimasukkan
implisit (benar) atau eksplisit (salah). Dalam kasus sebelumnya, energik
kontribusi dari orbital inti beku dilipat menjadi satu elektron
integral dan menjadi "energi inti beku" yang dihitung oleh program transformasi.
Defaultnya benar.

EKSPOR_VECTOR = boolean
Ini menentukan apakah akan menyimpan vektor konvergensi di akhir proses. NS
vektor disimpan dalam format transparan sehingga program lain dapat menggunakan
dengan mudah. Formatnya ditentukan dalam src/lib/libqt/slatardset.h. Standarnya adalah
Salah.

NUM_EKSPOR = bilangan bulat
If EKSPOR_VECTOR disetel ke true, this menentukan jumlah vektor yang
perlu diekspor pada akhir proses. Standarnya adalah 1.

GUESS_VECTOR = string
Ini menentukan jenis vektor tebakan yang akan digunakan dalam iterasi CI. Saat ini
hanya digunakan dengan metode iterasi SEM. Nilai yang diterima adalah UNIT untuk vektor satuan
Tebak (NUM_ROOTS dan NUM_INIT_VECS keduanya harus 1); H0_BLOK menggunakan vektor eigen
dari submatriks BLOK H0 (default); DFILE untuk menggunakan NUM_ROOTS konvergen sebelumnya
vektor dalam file D; dan MP2 untuk menggunakan fungsi gelombang MP2 sebagai tebakan (tidak berfungsi
saat ini).

H0_BLOCKSIZE = bilangan bulat
Parameter ini menentukan ukuran blok "H0" dari Hamiltonian yang
diselesaikan dengan tepat. n determinan dengan energi SCF terendah dipilih, dan a
submatriks Hamiltonian dibentuk menggunakan determinan ini. Submatriks ini adalah
digunakan untuk mempercepat konvergensi iterasi CI di OLSEN dan MITRUSHENKOV
skema iterasi, dan juga untuk menemukan tebakan awal yang baik untuk SEM metode jika
GUESS_VECTOR = H0_BLOK. Default ke 40. Perhatikan bahwa program dapat mengubah
diberikan ukuran untuk Ms = 0 kasus (Ms0 = BENAR) jika menentukan bahwa blok H0 termasuk
hanya satu anggota dari sepasang determinan yang terkait dengan simetri pembalikan waktu. Untuk
ukuran blok yang sangat kecil, ini mungkin bisa menghilangkan seluruh blok H0; NS
program harus mencetak peringatan jika ini terjadi.

H0_BLOCK_COUPLING_SIZE = bilangan bulat
Parameter yang menentukan ukuran blok kopling dalam generalisasi
prasyarat davidson. Nilai defaultnya adalah 1000.

MAX_DET = bilangan bulat
Menetapkan jumlah maksimum determinan; jika ruang CI lebih besar dari ini,
program dibatalkan. Opsi ini ada untuk memastikan bahwa perhitungan yang sangat besar tidak
dijalankan secara tidak sengaja. Selama fase pengembangan saat ini, standarnya adalah 10000, tetapi
itu akan dibangkitkan tak lama lagi.

MAKSIMUM = bilangan bulat
Jumlah maksimum iterasi untuk mendiagonalisasi Hamiltonian. Default ke 12.

Ms0 = boolean
If BENAR, gunakan komponen status Ms=0. Default untuk BENAR jika kulit tertutup dan
untuk SALAH sebaliknya. Terkait dengan S parameter.

NPRINT = bilangan bulat
Nilai ini menentukan jumlah determinan yang akan dicetak, bersama dengan
koefisien mereka, dalam daftar penentu paling penting dalam CI akhir
vektor. Standarnya adalah 20.

NUM_ROOTS = bilangan bulat
Nilai ini memberikan jumlah akar yang akan diperoleh dari sekuler
persamaan. Standarnya adalah satu. Jika diperlukan lebih dari satu root, atur DIAG_METHOD
untuk SEM (atau, untuk kasus yang sangat kecil, RSP or UJI SEMT).

NUM_INIT_VECS = bilangan bulat
Jumlah vektor awal yang akan digunakan dalam prosedur iterasi CI. Default untuk
jumlah akar.

OPDM = boolean
If BENAR hitung matriks kepadatan satu partikel dan buat OPDM_MENULIS default untuk
BENAR. Nilai default dari OPDM is SALAH.

OPDM_FILE = bilangan bulat
File (nomor unit) untuk menulis matriks kepadatan satu partikel jika OPDM_MENULIS =
BENAR. Nilai default saat ini adalah 73.

OPDM_MENULIS = boolean
Tandai apakah akan menulis matriks kepadatan satu partikel ke disk atau tidak.

OPDM_PRINT = boolean
Tandai apakah matriks kepadatan satu partikel akan dicetak atau tidak.

OPDM_DIAG = boolean
Tandai apakah matriks kerapatan satu partikel akan didiagonalisasi atau tidak.

WRTNO = boolean
Tandai apakah akan menulis orbital alami CI ke PSIF_CHKPT atau tidak.

FILE ORBS = bilangan bulat
File (nomor unit) untuk menulis berbagai orbital alami CI. Nilai defaultnya adalah
76.

OPDM_AVE = boolean
Tandai apakah akan merata-ratakan OPDM pada beberapa akar untuk mendapatkan a
matriks kepadatan satu partikel rata-rata negara bagian. Matriks densitas ini dapat didiagonalisasi
untuk mendapatkan orbital alami CI.

ORBS_ROOT = bilangan bulat
Bendera untuk menetapkan nomor akar tempat penulisan orbital alami CI
PSIF_CHKPT. Nilai default adalah 1 (akar terendah).

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).

ROOT = bilangan bulat
Akar untuk menuliskan matriks kerapatan dua partikel untuk (kerapatan satu partikel
matriks ditulis untuk semua akar). Berguna untuk CASSCF atau CI . khusus negara bagian
optimasi pada keadaan tereksitasi.

S = bilangan bulat
Nilai bilangan kuantum spin S diberikan oleh opsi ini. Standarnya adalah 0
(kaos). Satu-satunya hal yang sebenarnya digunakan adalah menentukan fase
setengah redundan dari vektor CI ketika komponen Ms=0 digunakan (yaitu, Ms0 =
BENAR). Untuk kasus di mana S bukan bilangan bulat, parameter ini tidak perlu dimasukkan
karena keadaan seperti itu tidak dapat memiliki komponen Ms=0.

TPDM = boolean
If BENAR hitung matriks kerapatan dua partikel dan buat TPDM_MENULIS default untuk
BENAR. Nilai default dari TPDM is SALAH.

TPDM_FILE = bilangan bulat
File (nomor satuan) untuk menulis matriks kerapatan dua partikel jika TPDM_MENULIS =
BENAR. Nilai default saat ini adalah 74.

TPDM_MENULIS = boolean
Tandai apakah akan menulis matriks kepadatan dua partikel ke disk atau tidak.

TPDM_PRINT = boolean
Tandai apakah matriks densitas dua partikel akan dicetak atau tidak. Biasanya sangat
ide buruk kecuali untuk men-debug kasus kecil.

Ada juga beberapa input yang kurang umum digunakan, yang tidak perlu dilakukan oleh pemula yang menggunakan PSI
digunakan.

BENDAZZOLI = boolean
Gunakan beberapa rutinitas untuk menghitung sigma berdasarkan makalah Bendazzoli et al.
Tampaknya lebih lambat dan tidak bermanfaat; mungkin menghilang pada akhirnya. Bekerja hanya untuk
CI penuh dan saya tidak ingat apakah saya bisa melihat bagaimana skema pintar mereka
diperluas ke RAS pada umumnya.

CALC_SSQ = boolean
Jika BENAR, hitung nilai harapan operator S^2 untuk CI akhir
fungsi gelombang untuk setiap akar. Pada prinsipnya, DETCI harus menghasilkan fungsi eigen S^2.
Nilai defaultnya salah.

COLLAPSE_SIZE bilangan bulat
Memberikan jumlah vektor untuk dipertahankan ketika subruang Davidson diciutkan (lihat
MAXNVECT di bawah). Jika lebih besar dari satu, subruang yang diciutkan mempertahankan yang terbaik
estimasi vektor CI untuk n iterasi sebelumnya. Default ke 1.

PERTAMA_TMP_UNIT = bilangan bulat
Memberikan nomor file (unit) yang terkait dengan file awal pertama yang digunakan oleh DETCI.
File awal lainnya diberi nomor secara berurutan dari titik ini, sesuai urutannya
H(diag), C, S, D. Masing-masing file logis ini membutuhkan sejumlah file fisik
ditentukan oleh parameter input yang lebih tidak jelas NUM_HD_TMP_UNIT,
NUM_C_TMP_UNITS, NUM_S_TMP_UNIT, NUM_D_TMP_UNIT. Pengguna juga dapat menentukan
titik awal yang berbeda untuk masing-masing set ini menggunakan parameter
PERTAMA_HD_TMP_UNIT Dan seterusnya. Memisahkan file di beberapa unit dapat membantu
menghindari masalah ukuran bilangan bulat dalam menangani file besar yang ada di
DETCI dan di perpustakaan I/O PSI; tapi sekali lagi, saya belum mengujinya untuk melihat apa
terjadi. Unit pertama dari setiap bagian dicetak di bawah judul FILES in
output parameter memulai proses DETCI.

FZC = boolean
Menentukan apakah orbital inti beku diperlakukan sebagai benar-benar beku (yaitu,
absen sama sekali dari perhitungan, FZC = BENAR) atau apakah mereka hadir tetapi
dibatasi untuk ditempati ganda (FZC = SALAH). Dalam program GUGA CI, ini adalah
perbedaan antara apa yang disebut orbital FZC dan COR. Umumnya,
integral untuk orbital inti beku tidak diperlukan oleh DETCI tetapi mungkin diperlukan
untuk MCSCF atau gradien.

ICORE = bilangan bulat
Menentukan cara menangani buffering vektor CI. Nilai 0 membuat program
melakukan I/O satu subblok RAS pada satu waktu; 1 menggunakan seluruh vektor CI sekaligus; dan 2
menggunakan satu blok irep pada suatu waktu. Nilai 0 atau 2 menyebabkan beberapa inefisiensi dalam
I/O (memerlukan beberapa pembacaan vektor C saat membangun H dalam iteratif
subruang jika DIAG_METHOD = SEM), tetapi membutuhkan lebih sedikit memori inti.

SAYA BERHENTI = boolean
If BENAR maka DETCI akan berhenti setelah informasi string terbentuk dan sebelum
integral dibaca. Mungkin akhirnya berubah menjadi bilangan bulat sehingga pengguna dapat
pilih dari beberapa titik pemberhentian.

MAXNVECT = bilangan bulat
Memberikan jumlah maksimum vektor subruang Davidson yang dapat disimpan pada disk untuk
koefisien CI dan vektor sigma. (Ada satu vektor H(diag) dan bilangan
vektor D sama dengan jumlah akar). Ketika jumlah vektor pada disk
mencapai nilai MAXNVECT, subruang Davidson akan diciutkan menjadi
COLLAPSE_SIZE vektor untuk setiap akar. Ini sangat membantu untuk menghemat ruang disk.
Default untuk MAKSIMUM * NUM_ROOTS + NUM_INIT_VECS.

CAMPURAN = boolean
Ini menentukan apakah eksitasi RAS II/RAS III "campuran" diizinkan masuk ke CI
ruang angkasa. Ini berguna untuk menempatkan batasan tambahan pada RAS CI.

CAMPURAN4 = boolean
Ini mirip dengan kata kunci MIXED, tetapi mengacu pada eksitasi ke dalam RAS IV.

NUNIT = bilangan bulat
Jumlah file awal yang akan digunakan dalam menyimpan vektor C (dan juga untuk sigma
vektor).

OEI_ERASE = boolean
Ini menentukan apakah program akan menghapus file integral satu elektron setelahnya
telah dibaca. Defaultnya pada akhirnya akan menjadi true, tetapi selama pengembangan
defaultnya salah.

OEI_FILE = bilangan bulat
Kata kunci ini memungkinkan pengguna untuk menentukan file integral satu elektron yang diubah.
Standarnya adalah 71.

PRINT_CIBLKS = boolean
Menentukan apakah program harus mencetak ringkasan semua blok di
Vektor CI (yang dapat dicetak ke dalam bentuk matriks, lihat referensi.)

R4S = boolean
Membatasi string RAS IV ke set minimal, menghemat memori. Jika Anda
khawatir tentang opsi ini, Anda harus menulis kepada David untuk meminta nasihat kecuali jika Anda
ahli DETCI.

REF_SYM = bilangan bulat
Opsi ini memungkinkan pengguna untuk mencari vektor CI dari irep yang berbeda dari
referensi. Ini mungkin hanya masuk akal untuk CI Penuh, dan mungkin tidak
bekerja dengan tebakan vektor satuan. Penomoran dimulai dari nol untuk total-
simetris

REPL_OTF = boolean
Memberi tahu DETCI apakah akan melakukan penggantian string dengan cepat atau tidak. Dapat menyimpan
jumlah memori yang sangat besar (terutama untuk CI yang terpotong) tetapi agak terkelupas dan
belum diuji untuk sementara waktu. Seingat saya, itu hanya berfungsi untuk kelas-kelas tertentu
perhitungan RAS. Hubungi David untuk bantuan. Akhirnya, on-the-fly
barang pengganti harus dikerjakan ulang dengan cara yang jauh lebih cerdas sehingga tidak memakan waktu
ribuan tahun waktu CPU. Pekerjaan di sepanjang garis ini pernah dimulai dan dapat diselesaikan
akhirnya.

MENGULANG KEMBALI = boolean
Opsi ini memungkinkan pengguna untuk melanjutkan iterasi DETCI yang dihentikan
sebelum waktunya. Diasumsikan bahwa vektor CI dan sigma berada pada disk; jumlah
vektor ditentukan oleh RESTART_VECS diciutkan menjadi satu vektor per akar.

RESTART_VECS = bilangan bulat
If MENGULANG KEMBALI = BENAR ini menentukan jumlah vektor CI (dan sigma) untuk dibaca dari
disk. Biasanya ini adalah jumlah iterasi yang berhasil diselesaikan dari a
run sebelumnya dikalikan dengan jumlah root untuk run tersebut.

TEI_ERASE = boolean
Ini menentukan apakah program akan menghapus file integral dua elektron setelahnya
telah dibaca. Defaultnya pada akhirnya akan menjadi true, tetapi selama pengembangan
defaultnya salah.

TEI_FILE = bilangan bulat
Kata kunci ini memungkinkan pengguna untuk menentukan file integral dua elektron yang diubah.
Standarnya adalah 72.

MPN = boolean
Jika opsi ini BENAR, DETCI akan menghitung deret MPn ke urutan ke-k di mana k
ditentukan oleh maxnvect. Untuk sistem kulit terbuka (REF=ROHF, WFN = ZAPTN), DETCI
akan menghitung seri ZAPTn. GUESS_VECTOR harus disetel ke TRUE. HD_OTF harus
disetel ke BENAR. HD_AVE harus disetel ke orb_ener.

SIMPAN_MPN2 = bilangan bulat
Ketika MPN adalah TRUE dan WIGNER adalah TRUE maka opsi ini menjadi valid. Jika diatur ke 1
maka energi MP(2n-1) disimpan. Jika diatur ke 2 maka energi MP(2n-2) akan dihemat. Jika ada
nilai lain energi MPn disimpan. Standarnya adalah 0.

9 Feb, 1996 detci(1)

Gunakan detci online menggunakan layanan onworks.net