Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Variational Localized Active Space Self-Consistent Field Method
67
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Fragmentation methods applied to multireference wave functions constitute a road towards the application of highly accurate ab initio wave function calculations to large molecules and solids. However, it is important for reproducibility and transferability that a fragmentation scheme be well-defined with minimal dependence on initial orbital guesses or user-designed ad hoc fragmentation schemes. One way to improve this sort of robustness is to ensure the energy obeys a variational principle; i.e., that the active orbitals and active space wave functions minimize the electronic energy in a certain ansatz for the molecular wave function. We extended the theory of the localized active space self-consistent field, LASSCF, method (JCTC 2019, 15, 972) to fully minimize the energy with respect to all orbital rotations, rendering it truly variational. The new method, called vLASSCF, substantially improves the robustness and reproducibility of the LAS wave function compared to LASSCF. We analyze the storage and operation cost scaling of vLASSCF compared to orbital optimization using a standard CASSCF approach and we show results of vLASSCF calculations on some simple test systems. We show that vLASSCF is energetically equivalent to CASSCF in the limit of one active subspace, and that vLASSCF significantly improves upon the reliability of LASSCF energy differences, allowing for more meaningful and subtle analysis of potential energy curves of dissociating molecules. We also show that all forms of LASSCF have a lower operation cost scaling than the orbital-optimization part of CASSCF.
rate research
Read More
We report an efficient algorithm using density fitting for the relativistic complete active space self-consistent field (CASSCF) method, which is significantly more stable than the algorithm previously reported by one of the authors [J. E. Bates and T. Shiozaki, J. Chem. Phys. 142, 044112 (2015)]. Our algorithm is based on the second-order orbital update scheme with an iterative augmented Hessian procedure, in which the density-fitted orbital Hessian is directly contracted to the trial vectors. Using this scheme, each microiteration is made less time consuming than one Dirac-Hartree-Fock iteration, and macroiterations converge quadratically. In addition, we show that the CASSCF calculations with the Gaunt and full Breit interactions can be efficiently performed by means of approximate orbital Hessians computed with the Dirac-Coulomb Hamiltonian. It is demonstrated that our algorithm can also be applied to systems under an external magnetic field, for which all of the molecular integrals are computed using gauge-including atomic orbitals.
We present the numerical implementation of the time-dependent complete-active-space self-consistent-field (TD-CASSCF) method [Phys. Rev. A, 88, 023402 (2013)] for atoms driven by a strong linearly polarized laser pulse. The present implementation treats the problem in its full dimensionality and introduces a gauge-invariant frozen-core approximation, an efficient evaluation of the Coulomb mean field scaling linearly with the number of basis functions, and a split-operator method specifically designed for stable propagation of stiff spatial derivative operators. We apply this method to high-harmonic generation in helium, beryllium, and neon and explore the role of electron correlations.
An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large complete active spaces (CAS). The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the CAS, meanwhile maintaining full spin symmetry; (2) the use of Jacobi rotation for the optimization of active orbitals, in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Just like selected iCI being a very accurate approximation to CASCI, iCISCF(2) is a very accurate approximation to CASSCF. Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the performances of iCISCF and iCISCF(2).
We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wavefunction. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) emph{ansatz}, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.
We present a first principles molecular dynamics approach that is based on time-reversible ex- tended Lagrangian Born-Oppenheimer molecular dynamics [Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) costruction are required in each integration time step. The proposed dy- namics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents an ideal starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
