Mario Piris

A natural orbital functional for multiconfigurational states

An explicit formulation of the Piris cumulant λΔ,Π matrix is described herein, and used to reconstruct the two-particle reduced density matrix (2-RDM). Then, we have derived a natural orbital functional, the Piris Natural Orbital Functional 5, PNOF5, constrained to fulfill the D, Q, and G positivity necessary conditions of the N-representable 2-RDM. This functional yields a remarkable accurate description of systems bearing substantial (near)degeneracy of one-particle states. The theory is applied to the homolitic dissociation of selected diatomic molecules and to the rotation barrier of ethylene, both paradigmatic cases of near-degeneracy effects. It is found that the method describes correctly the dissociation limit yielding an integer number of electrons on the dissociated atoms. PNOF5 predicts a barrier of 65.6 kcal/mol for the ethylene torsion in an outstanding agreement with Complete Active Space Second-order Perturbation Theory (CASPT2). The obtained occupation numbers and pseudo one-particle energies at the ethylene transition state account for fully degenerate π orbitals. The calculated equilibrium distances, dipole moments, and binding energies of the considered molecules are presented. The values obtained are accurate comparing those obtained by the complete active space self-consistent field method and the experimental data.

Communication: The role of the positivity N-representability conditions in natural orbital functional theory

The positivity conditions for the N-representability of the reduced density matrices are considered to propose a new natural orbital functional. The Piris reconstruction functional, which is based on an explicit form of the two-particle cumulant λ(Δ,Π) is used to reconstruct the two-particle reduced density matrix. A new approach for Π matrix, satisfying rigorously D, Q, and G necessary conditions, leads to Piris Natural Orbital Functional 4 (PNOF4). The theory is applied to the dissociation of selected diatomic molecules. The equilibrium distances, dipole moments, harmonic frequencies, anharmonicity constants, and binding energies of the considered molecules are presented. The values we have obtained are very accurate results comparing with the experimental data.

Communications: Accurate description of atoms and molecules by natural orbital functional theory

The spin-conserving density matrix functional theory is used to propose an improved natural orbital functional. The Piris reconstruction functional, PNOF, which is based on an explicit form of the two-particle cumulant lambda(Delta,Lambda) satisfying necessary positivity conditions for the two-particle reduced density matrix, is used to reconstruct the latter. A new approach Lambda((3)), as well as an extension of the known Delta(alphabeta) to spin-uncompensated systems lead to PNOF3. The theory is applied to the calculation of the total energies of the first- and second-row atoms (H-Ne) and a number of selected small molecules. The energy differences between the ground state and the lowest-lying excited state with different spin for these atoms, and the atomization energies of the considered molecules are also presented. The obtained values agree remarkably well with their corresponding both CCSD(T, full) and experimental values.

Iterative diagonalization for orbital optimization in natural orbital functional theory.

A challenging task in natural orbital functional theory is to find an efficient procedure for doing orbital optimization. Procedures based on diagonalization techniques have confirmed its practical value since the resulting orbitals are automatically orthogonal. In this work, a new procedure is introduced, which yields the natural orbitals by iterative diagonalization of a Hermitian matrix F. The off‐diagonal elements of the latter are determined explicitly from the hermiticity of the matrix of the Lagrange multipliers. An expression for diagonal elements is absent so a generalized Fockian is undefined in the conventional sense, nevertheless, they may be determined from an aufbau principle. Thus, the diagonal elements are obtained iteratively considering as starting values those coming from a single diagonalization of the matrix of the Lagrange multipliers calculated with the Hartree‐Fock orbitals after the occupation numbers have been optimized. The method has been tested on the G2/97 set of molecules for the Piris natural orbital functional. To help the convergence, we have implemented a variable scaling factor which avoids large values of the off‐diagonal elements of F. The elapsed times of the computations required by the proposed procedure are compared with a full sequential quadratic programming optimization, so that the efficiency of the method presented here is demonstrated.

Dispersion interactions within the Piris natural orbital functional theory: the helium dimer.

The authors have investigated the description of the dispersion interaction within the Piris natural orbital functional (PNOF) theory. The PNOF arises from an explicit antisymmetric approach for the two-particle cumulant in terms of two symmetric matrices, Delta and Lambda. The functional forms of these matrices are obtained from the generalization of the two-particle system expressions, except for the off-diagonal elements of Delta. The mean value theorem and the partial sum rule obtained for the off-diagonal elements of Delta provide a prescription for deriving practical functionals. In particular, the previous employed approximation {Jpp/2} for the mean values {Jp*} affords several molecular properties but it is incapable to account for dispersion effects. In this work, the authors analyze a new approach for Jp* obtained by factorization of the matrix Delta within the bounds on its off-diagonal elements imposed by the positivity conditions of the two-particle reduced density matrix. Additional terms for the matrix elements of Lambda proportional to the square root of the holes are again introduced to describe properly the occupation numbers of the lowest occupied levels. The authors have found that the cross products between weakly occupied orbitals must be removed from the functional form of Lambda to obtain a correct long-range asymptotic behavior. The PNOF is used to predict the binding energy as well as the equilibrium distance of the helium dimer. The results are compared with the full configuration-interaction calculations and the corresponding experimental data.

Spin conserving natural orbital functional theory.

The natural orbital functional theory is considered for spin uncompensated systems, i.e., systems that have one or more unpaired electrons. The well-known cumulant expansion is used to reconstruct the two-particle reduced density matrix. A new condition to ensure the conservation of the total spin is obtained for the two-particle cumulant matrix. An extension of the Piris natural orbital functional 1 (PNOF1), based on an explicit form for the cumulant, to spin uncompensated systems is also considered. The theory is applied to the calculation of energy differences between the ground state and the lowest lying excited state with different spins for first-row atoms (Li, Be, B, C, N, O, and F) and diatomic oxygen molecule (O(2)). The values we obtained are very accurate results as compared to the CCSD(T) method and the experimental data.

A GENERALIZED SELF-CONSISTENT-FIELD PROCEDURE IN THE IMPROVED BCS THEORY

AbstractThe basic ideas of the Improved Bardeen–Cooper–Schrieffer (IBCS) approach to the first‐ and second‐order Reduced Density Matrices (1‐ and 2‐RDM) are briefly reviewed. The molecular orbital occupations

Performance of PNOF5 Natural Orbital Functional for Radical Formation Reactions: Hydrogen Atom Abstraction and C-C and O-O Homolytic Bond Cleavage in Selected Molecules.

\left\{ {\rho _q } \right\}