Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Electron localization function at the correlated level

The electron localization function (ELF) has been proven so far a valuable tool to determine the location of electron pairs. Because of that, the ELF has been widely used to understand the nature of the chemical bonding and to discuss the mechanism of chemical reactions. Up to now, most applications of the ELF have been performed with monodeterminantal methods and only few attempts to calculate this function for correlated wave functions have been carried out. Here, a formulation of ELF valid for mono- and multiconfigurational wave functions is given and compared with previous recently reported approaches. The method described does not require the use of the homogeneous electron gas to define the ELF, at variance with the ELF definition given by Becke. The effect of the electron correlation in the ELF, introduced by means of configuration interaction with singles and doubles calculations, is discussed in the light of the results derived from a set of atomic and molecular systems

Aplicació de la semblança molecular quàntica en la reducció de l'espai configuracional per a l'estat fonamental i primers excitats de l'àtom d'Heli

Es presenta un mètode de selecció d'orbitals atòmics relacionat amb la teoria de la Semblança Molecular Quàntica, que permet reduir l'espai actiu quan es vol dur a terme un càlcul a nivell d'Interacció de Configuracions per a l'àtom d'heli

The mapping of the local contributions of Fermi and Coulomb correlation into intracule and extracule density distributions

The contributions of the correlated and uncorrelated components of the electron-pair density to atomic and molecular intracule I(r) and extracule E(R) densities and its Laplacian functions ∇2I(r) and ∇2E(R) are analyzed at the Hartree-Fock (HF) and configuration interaction (CI) levels of theory. The topologies of the uncorrelated components of these functions can be rationalized in terms of the corresponding one-electron densities. In contrast, by analyzing the correlated components of I(r) and E(R), namely, IC(r) and EC(R), the effect of electron Fermi and Coulomb correlation can be assessed at the HF and CI levels of theory. Moreover, the contribution of Coulomb correlation can be isolated by means of difference maps between IC(r) and EC(R) distributions calculated at the two levels of theory. As application examples, the He, Ne, and Ar atomic series, the C2-2, N2, O2+2 molecular series, and the C2H4 molecule have been investigated. For these atoms and molecules, it is found that Fermi correlation accounts for the main characteristics of IC(r) and EC(R), with Coulomb correlation increasing slightly the locality of these functions at the CI level of theory. Furthermore, IC(r), EC(R), and the associated Laplacian functions, reveal the short-ranged nature and high isotropy of Fermi and Coulomb correlation in atoms and molecules

A variational approach for calculating Franck-Condon factors including mode-mode anharmonic coupling

We have implemented our new procedure for computing Franck-Condon factors utilizing vibrational configuration interaction based on a vibrational self-consistent field reference. Both Duschinsky rotations and anharmonic three-mode coupling are taken into account. Simulations of the first ionization band of Cl O2 and C4 H4 O (furan) using up to quadruple excitations in treating anharmonicity are reported and analyzed. A developer version of the MIDASCPP code was employed to obtain the required anharmonic vibrational integrals and transition frequencies

Correcting Measurement Error Bias in Interaction Models with Small Samples

Several methods have been suggested to estimate non-linear models with interaction terms in the presence of measurement error. Structural equation models eliminate measurement error bias, but require large samples. Ordinary least squares regression on summated scales, regression on factor scores and partial least squares are appropriate for small samples but do not correct measurement error bias. Two stage least squares regression does correct measurement error bias but the results strongly depend on the instrumental variable choice. This article discusses the old disattenuated regression method as an alternative for correcting measurement error in small samples. The method is extended to the case of interaction terms and is illustrated on a model that examines the interaction effect of innovation and style of use of budgets on business performance. Alternative reliability estimates that can be used to disattenuate the estimates are discussed. A comparison is made with the alternative methods. Methods that do not correct for measurement error bias perform very similarly and considerably worse than disattenuated regression

Desenvolupament de mètodes de preconcentració emprant membranes líquides suportades i extracció en fase sòlida per a la determinació de l'herbicida glifosat i el seu metabòlit AMPA en aigües naturals

El glifosat, N-(fosfonometil) glicina, és un dels herbicides més utilitzats arreu del món a causa de la seva baixa toxicitat i al seu ampli espectre d'aplicació. A conseqüència del gran ús que se'n fa, és necessari monitoritzar aquest compost i el seu principal metabòlit, l'àcid aminometilfosfònic (AMPA), en el medi ambient. S'han descrit diversos mètodes instrumentals basats en cromatografia de gasos (GC) i de líquids (HPLC), sent aquesta darrera l'opció més favorable a causa del caràcter polar dels anàlits. Per assolir nivells de concentració baixos cal, però, la preconcentració dels anàlits. En aquest treball s'estudien diferents alternatives amb aquest objectiu. S'ha avaluat la tècnica de membrana líquida suportada (SLM) on la membrana consisteix en una dissolució orgànica, que conté un transportador (en el nostre cas, un bescanviador d'anions comercial, Aliquat 336), que impregna un suport polimèric microporós que se situa entre dues solucions aquoses: la de càrrega, que conté els anàlits inicialment, i la receptora, on es retenen els anàlits després del seu transport a través de la membrana. Les condicions d'extracció més adequades s'obtenen treballant en medi bàsic amb NaOH on els anàlits estan en forma aniònica i les majors recuperacions s'obtenen amb HCl 0,1 M o NaCl 0,5 M, la qual cosa indica que l'ió clorur és la força impulsora del transport. Un cop dissenyat el sistema, es duen a terme experiments de preconcentració amb dues geometries diferents: un sistema de membrana laminar (LSLM) on recircula la fase receptora i un sistema de fibra buida (HFSLM). Els millors resultats s'obtenen amb el mòdul de fibra buida, amb factors de concentració de 25 i 3 per a glifosat i AMPA, respectivament, fent recircular durant 24 hores 100 ml de solució de càrrega i 4 ml de solució receptora. També s'aplica una tècnica més selectiva, la cromatografia d'afinitat amb ió metàl·lic immobilitzat (IMAC), basada en la interacció entre els anàlits i un metall immobilitzat en una resina a través d'un grup funcional d'aquesta. En aquest estudi s'immobilitza pal·ladi al grup funcional 8-hidroxiquinoleïna de la resina amb matriu acrílica Spheron Oxine 1000 i s'avalua per a l'extracció i preconcentració de glifosat i AMPA. Per a ambdós anàlits l'adsorció és del 100 % i les recuperacions són superiors al 80 % i al 60 % per a glifosat i AMPA, respectivament, utilitzant HCl 0,1 M + NaCl 1 M com a eluent. Aquests resultats es comparen amb els obtinguts amb dues resines més, també carregades amb pal·ladi: Iontosorb Oxin 100, que té el mateix grup funcional però matriu de cel·lulosa, i Spheron Thiol 1000, on el grup funcional és un tiol i la matriu també és acrílica. Per al glifosat els resultats són similars amb totes les resines, però per a l'AMPA la resina Spheron Thiol és la única que proporciona recuperacions superiors al 93 %. Finalment, una altra opció estudiada és l'acoblament de dues columnes de cromatografia líquida (LC-LC). En l'estudi l'objectiu és millorar el mètode existent per a glifosat i AMPA en aigües naturals on el LOD era de 0,25 ug/l. El mètode consisteix en la derivatització precolumna amb el reactiu fluorescent FMOC i l'anàlisi amb l'acoblament LC-LC-fluorescència. Variant lleugerament les condicions de derivatització s'aconsegueix quantificar 0,1 ug/l de glifosat i AMPA. Es fortifiquen aigües naturals amb 0,1, 1 i 10 ug/l dels anàlits per validar el mètode. S'obtenen recuperacions d'entre el 85 % i el 100 %, amb desviacions estàndard relatives inferiors al 8 %. Aplicant una tècnica de preconcentració prèvia a la derivatització i anàlisi utilitzant una resina de bescanvi aniònic, Amberlite IRA-900, es millora la sensibilitat del mètode i s'assoleix un LOD per al glifosat de 0,02 ug/l.