detcas - Trực tuyến trên đám mây

CHƯƠNG TRÌNH:

TÊN

detcas - Tương tác cấu hình quyết định / Hoàn thành hoạt động
Chương trình trường tự nhất quán trong không gian

MÔ TẢ

Chương trình detcas thực hiện trường tự nhất quán hoàn toàn không gian hoạt động (CASSCF)
tối ưu hóa các obitan phân tử thông qua quy trình hai bước kết hợp với
chương trình tương tác cấu hình xác định DETCI. Chương trình khá đơn giản và
hiện đang sử dụng phương pháp tiếp cận Newton-Raphson để cập nhật các quỹ đạo, sử dụng
quỹ đạo gần đúng Hessian. Sự hội tụ được đẩy nhanh bằng cách sử dụng phép đảo ngược trực tiếp của Pulay đối với
thủ tục không gian con lặp lại (DIIS). Mã đã được viết để cho phép tổng quát hơn
các chức năng sóng không nhất thiết phải có xử lý CI đầy đủ của không gian hoạt động. Trong
đặc biệt, bất kỳ chức năng sóng CI không gian hoạt động hạn chế (RAS) nào được hỗ trợ bởi DETCI đều có thể
được sử dụng, cho phép các chức năng sóng MCSCF của loại RASSCF.

THAM KHẢO

Orbital Hessian gần đúng:

1. G. Chaban, MW Schmidt, và MS Gordon, Theor. Chim. Acta 97, 88-95 (1997).

CI không gian hoạt động bị hạn chế:

1. Thuật toán tương tác dựa trên cấu hình xác định cho hoàn chỉnh và hạn chế
Không gian tương tác cấu hình, J. Olsen, BO Roos, P. Jorgensen và HJ Aa.
Jensen, J. Chem. Thể chất. 89, 2185 (năm 1988).

Không gian hoạt động hạn chế SCF:

1. P.-A. Malmqvist, A. Rendell, và BO Roos, J. Phys. Chèm. 94, 5477 (1990).

Chương trình DETCI:

1. CD Sherrill, Thuật toán tính toán cho quy mô lớn đầy đủ và đa tham chiếu
Các chức năng sóng tương tác cấu hình, luận án tiến sĩ, Đại học Georgia, Athens,
GA, 1996.

2. CD Sherrill và HF Schaefer, Phương pháp tương tác cấu hình: Những tiến bộ
trong các Phương pháp Tiếp cận Có Tương quan Cao, Lời khuyên. Chem lượng tử. 34, 143-269 (1999).

CÁC TẬP TIN YÊU CẦU

input.dat - Tập tin đầu vào
file78 - Tích phân một electron MO (toán tử fzc)
file72 - Tích phân hai electron MO
file73 - Ma trận mật độ một hạt MO
file74 - Ma trận mật độ hai hạt MO
file75 - MO Lagrangian

CÁC TẬP TIN CẬP NHẬT

output.dat - Tệp đầu ra
file14.dat - Bản ghi năng lượng và độ dốc quỹ đạo

ĐẦU VÀO FORMAT

Các đối số dòng lệnh sau có sẵn:

-Yên lặng Điều này cho kết quả tương tự như PRINT = 0.

-o biệt danh
Cung cấp tên tệp cho tệp đầu ra. Mặc định là output.dat.

Đầu vào bổ sung cho chương trình này được đọc từ tệp input.dat. Thông dụng hơn
từ khóa là:

WFN = chuỗi
Các loại chức năng sóng được hỗ trợ là CASSCF và RASSCF.

CHUYỂN ĐỔI = số nguyên
Sự hội tụ mong muốn trên gradient quỹ đạo. Sự hội tụ đạt được khi RMS
của sai số trong gradient quỹ đạo nhỏ hơn 10 ** (- n). Giá trị mặc định là 4 cho
tính toán năng lượng và 7 cho gradient.

DOCC = số nguyên_mảng
Vectơ này cung cấp số lượng các obitan bị chiếm gấp đôi trong mỗi không thể. Có
không có mặc định.

SOCC = số nguyên_mảng
Vectơ này cung cấp số lượng các obitan bị chiếm riêng lẻ trong mỗi không thể. Có
không có mặc định.

NĂNG LƯỢNG_CONVERGENCE = số nguyên
Sự hội tụ mong muốn trên tổng năng lượng MCSCF. Giá trị mặc định là 7.

FROZEN_DOCC = số nguyên_mảng
Số lượng các obitan chiếm năng lượng thấp nhất trong mỗi obitan không thể điều chỉnh được
biểu diễn theo nghĩa đen sẽ bị đóng băng (không được cập nhật trong MCSCF). Các
Thứ tự cotton của các đại diện không thể tin được được sử dụng. Mặc định là
vectơ không.

FROZEN_UOCC = số nguyên_vector
Số lượng các obitan chưa sử dụng năng lượng cao nhất trong mỗi biểu diễn bất khả quy
theo đúng nghĩa đen sẽ bị đóng băng (không được cập nhật trong MCSCF). Mặc định là số XNUMX
vectơ.

RESTRICTED_DOCC = số nguyên_mảng
Số lượng các obitan chiếm năng lượng thấp nhất trong mỗi obitan không thể điều chỉnh được
đại diện sẽ được tối ưu hóa nhưng vẫn bị chiếm đóng gấp đôi trong MCSCF.
Các quỹ đạo này đến sau FROZEN_DOCC các quỹ đạo. Mặc định là số XNUMX
vectơ.

RESTRICTED_UOCC = số nguyên_mảng
Số lượng các obitan chưa sử dụng năng lượng cao nhất trong mỗi biểu diễn bất khả quy
sẽ được tối ưu hóa nhưng không bị trống trong MCSCF. Các quỹ đạo này đến
trước FROZEN_UOCC các quỹ đạo. Mặc định là vectơ không.

NCASITER = số nguyên
Số lần lặp lại tối đa để tối ưu hóa các obitan. Tùy chọn này nên be
quy định in các MẶC ĐỊNH phần of đầu vào, bởi vì nó cần được hiển thị với
chương trình điều khiển PSI. Mặc định là 1.

IN = số nguyên
Tùy chọn này xác định độ chi tiết của đầu ra. Giá trị 1 hoặc 2 chỉ định
in tối thiểu, giá trị 3 chỉ định in chi tiết. Giá trị của 4 hoặc 5 là
được sử dụng để gỡ lỗi. Không sử dụng mức 5 trừ khi trường hợp thử nghiệm rất nhỏ (ví dụ:
STO H2O CISD).

Các từ khóa ít được sử dụng hơn là:

DIIS_FREQ = số nguyên
Quy trình ngoại suy DIIS sẽ được thử mọi n các lần lặp lại. Mặc định
là 1.

DIIS_MAX_VECS = số nguyên
Số vectơ không gian con tối đa cho thủ tục DIIS. Sau số này
trong tổng số vectơ được đạt đến, các vectơ cũ hơn sẽ bị loại bỏ khỏi không gian con dưới dạng
cần thiết. Giá trị mặc định là 8.

DIIS_MIN_VECS = số nguyên
Số vectơ không gian con tối thiểu trước khi nội suy DIIS có thể là
đã thực hiện. Giá trị mặc định là 2. Sẽ không hợp lý nếu có các giá trị nhỏ hơn 2.

DIIS_START = số nguyên
Số lần lặp khi Đảo ngược trực tiếp của Pulay của không gian con lặp lại (DIIS)
nên bật quy trình tăng tốc hội tụ. Giá trị mặc định là 3.
Trước lần lặp này, vectơ không được thêm vào không gian con DIIS. Người đầu tiên
Bước DIIS sẽ không được thực hiện cho đến khi DIIS_MIN_VECS vectơ nằm trong không gian con DIIS.

LAG_FILE = số nguyên
Tệp (số đơn vị) để đọc ma trận lagrang. Giá trị mặc định là
hiện 75.

OPDM_FILE = số nguyên
Tệp (số đơn vị) để đọc ma trận mật độ một hạt. Giá trị mặc định
hiện là 73.

SCALE_GRAD = boolean
Cho biết có nên chia tỷ lệ gradient quỹ đạo theo quỹ đạo gần đúng (đường chéo) hay không
Hessian. Mặc định là đúng.

SCALE_STEP = thực
Hệ số tỷ lệ cho bước quay của quỹ đạo. Mặc định là 1.0.

TPDM_FILE = số nguyên
Tệp (số đơn vị) để đọc ma trận mật độ hai hạt. Giá trị mặc định
hiện là 74.

LEVEL_SHIFT = boolean
Cho biết có cho phép chuyển cấp của ma trận hessian hay không. Cái này có thể được sử dụng
để đảm bảo rằng hessian là xác định dương cho các lần lặp đầu tiên. Nếu như
hessian không phải là dương xác định gần hội tụ, thì đây có thể là dấu hiệu của
sự bất ổn định số trong MCSCF. Mặc định là đúng.

SHIFT = thực
Cho biết các phần tử đường chéo của hessian nên được dịch chuyển như thế nào đối với cấp độ
sự dịch chuyển. Mặc định là 0.01.

DETERM_MIN = thực
Giá trị tối thiểu cho phép đối với yếu tố quyết định của ma trận hessian nếu chuyển cấp
đang bật. Điều này không được sử dụng nếu LEVEL_SHIFT = FALSE. Mặc định là 0.00001.

BƯỚC_MAX = thực
Đây là vòng quay quỹ đạo đơn tối đa được phép. Mặc định là 0.30.

FORCE_STEP = boolean
Điều này cho phép người dùng ghi đè bước được tính toán và buộc một bước trong
hướng cụ thể. Điều này có thể hữu ích nếu bạn cố gắng loại bỏ phép tính
từ một điểm yên ngựa. Mặc định này sai.

FORCE_PAIR = số nguyên
Đây là chỉ số cho cặp độc lập sẽ được xoay vòng nếu
FORCE_STEP = TRUE. Điều này được bỏ qua nếu không. Không có mặc định.

FORCE_VALUE = thực
Đây là giá trị quay của quỹ đạo để quay các quỹ đạo được chỉ định bởi
FORCE_PAIR. Điều này chỉ được sử dụng nếu FORCE_STEP được đặt thành TRUE. Giá trị mặc định là
0.0.

CHECK_HESSIAN = boolean
Tính các giá trị riêng của hessian quỹ đạo. Điều này rất tốt cho việc kiểm tra
giải pháp điểm yên. Nên được sử dụng mà không cần chuyển cấp. Mặc định là
SAI.

EIGEN_VECTORS = boolean
Bạn có muốn tính toán eigenvectors của hessian quỹ đạo không. Đây là
bị bỏ qua nếu CHECK_HESSIAN = FALSE. Mặc định là FALSE.

8 Tháng Năm, 1998 detcas(1)

Sử dụng detcas trực tuyến bằng các dịch vụ onworks.net