Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Moller-Plesset perturbation theory

Bozkaya, Ugur
Turney, Justin M.
Yamaguchi, Yukio
Schaefer, Henry F.
Sherrill, C. David
Using a Lagrangian-based approach, we present a more elegant derivation of the equations necessary for the variational optimization of the molecular orbitals (MOs) for the coupled-cluster doubles (CCD) method and second-order Moller-Plesset perturbation theory (MP2). These orbital-optimized theories are referred to as OO-CCD and OO-MP2 (or simply "OD" and "OMP2" for short), respectively. We also present an improved algorithm for orbital optimization in these methods. Explicit equations for response density matrices, the MO gradient, and the MO Hessian are reported both in spin-orbital and closed-shell spin-adapted forms. The Newton-Raphson algorithm is used for the optimization procedure using the MO gradient and Hessian. Further, orbital stability analyses are also carried out at correlated levels. The OD and OMP2 approaches are compared with the standard MP2, CCD, CCSD, and CCSD(T) methods. All these methods are applied to H(2)O, three diatomics, and the O(4)(+) molecule. Results demonstrate that the CCSD and OD methods give nearly identical results for H(2)O and diatomics; however, in symmetry-breaking problems as exemplified by O(4)(+), the OD method provides better results for vibrational frequencies. The OD method has further advantages over CCSD: its analytic gradients are easier to compute since there is no need to solve the coupled-perturbed equations for the orbital response, the computation of one-electron properties are easier because there is no response contribution to the particle density matrices, the variational optimized orbitals can be readily extended to allow inactive orbitals, it avoids spurious second-order poles in its response function, and its transition dipole moments are gauge invariant. The OMP2 has these same advantages over canonical MP2, making it promising for excited state properties via linear response theory. The quadratically convergent orbital-optimization procedure converges quickly for OMP2, and provides molecular properties that are somewhat different than those of MP2 for most of the test cases considered (although they are similar for H(2)O). Bond lengths are somewhat longer, and vibrational frequencies somewhat smaller, for OMP2 compared to MP2. In the difficult case of O(4)(+), results for several vibrational frequencies are significantly improved in going from MP2 to OMP2. (C) 2011 American Institute of Physics. [doi:10.1063/1.3631129]


Linear-linear basis functions for MLFMA solutions of magnetic-field and combined-field integral equations
Ergül, Özgür Salih (2007-04-01)
We present the linear-linear (LL) basis functions to improve the accuracy of the magnetic-field integral equation (MFIE) and the combined-field integral equation (CFIE) for three-dimensional electromagnetic scattering problems involving closed conductors. We consider the solutions of relatively large scattering problems by employing the multilevel fast multipole algorithm. Accuracy problems of MFIE and CFIE arising from their implementations with the conventional Rao-Wilton-Glisson (RWG) basis functions can...
Numerical solution of nonlinear reaction-diffusion and wave equations
Meral, Gülnihal; Tezer, Münevver; Department of Mathematics (2009)
In this thesis, the two-dimensional initial and boundary value problems (IBVPs) and the one-dimensional Cauchy problems defined by the nonlinear reaction- diffusion and wave equations are numerically solved. The dual reciprocity boundary element method (DRBEM) is used to discretize the IBVPs defined by single and system of nonlinear reaction-diffusion equations and nonlinear wave equation, spatially. The advantage of DRBEM for the exterior regions is made use of for the latter problem. The differential quad...
Exact Pseudospin Symmetric Solution of the Dirac Equation for Pseudoharmonic Potential in the Presence of Tensor Potential
AYDOĞDU, OKTAY; Sever, Ramazan (Springer Science and Business Media LLC, 2010-04-01)
Under the pseudospin symmetry, we obtain exact solution of the Dirac equation for the pseudoharmonic potential in the presence of the tensor potential with arbitrary spin-orbit coupling quantum number kappa. The energy eigenvalue equation of the Dirac particles is found and the corresponding radial wave functions are presented in terms of confluent hypergeometric functions. We investigate the tensor potential dependence of the energy of the each state in the pseudospin doublet. It is shown that degeneracy b...
IKHDAİR, SAMEER; Sever, Ramazan (World Scientific Pub Co Pte Lt, 2009-03-01)
The point canonical transformation (PCT) approach is used to solve the Schrodinger equation for an arbitrary dimension D with a power-law position-dependent effective mass (PDEM) distribution function for the pseudoharmonic and modified Kratzer (Mie-type) diatomic molecular potentials. In mapping the transformed exactly solvable D-dimensional (D >= 2) Schrodinger equation with constant mass into the effective mass equation by using a proper transformation, the exact bound state solutions including the energ...
Relativistic Two-Dimensional Harmonic Oscillator Plus Cornell Potentials in External Magnetic and AB Fields
Ikhdair, Sameer M.; Sever, Ramazan (Hindawi Limited, 2013)
<jats:p>The Klein-Gordon (KG) equation for the two-dimensional scalar-vector harmonic oscillator plus Cornell potentials in the presence of external magnetic and Aharonov-Bohm (AB) flux fields is solved using the wave function ansatz method. The exact energy eigenvalues and the wave functions are obtained in terms of potential parameters, magnetic field strength, AB flux field, and magnetic quantum number. The results obtained by using different Larmor frequencies are compared with the results in the absenc...
Citation Formats
U. Bozkaya, J. M. Turney, Y. Yamaguchi, H. F. Schaefer, and C. D. Sherrill, “Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Moller-Plesset perturbation theory,” JOURNAL OF CHEMICAL PHYSICS, pp. 0–0, 2011, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/68080.