[molpro-user] Linear dependent basis
Kirk Peterson
kipeters at wsu.edu
Mon Jul 12 16:45:43 BST 2010
John,
could you send your basis set block to me?
best wishes,
Kirk
On Jul 12, 2010, at 3:11 AM, John Keith wrote:
> Hi all,
>
> I'm trying to run single reference ab initio calculations on
> transition metal complexes and build up enough calculations to
> extrapolate to the CBS limit.
>
> My molecule is the PdCl4 dianion, and I'm using the Stuttgart Pd
> pseudopotential and the (aug)-cc-pvNZ basis sets on all other valence
> electrons. All are taken from the EMSL basis set library. I get SCF
> convergence with cc-pvDZ, aug-cc-pvDZ, cc-pvTZ, and aug-cc-pvTZ basis
> sets, but at cc-pvQZ I run into an error saying:
>
> ------------------------
> Eigenvalues of metric
>
> 1 0.212E-08 0.406E-05 0.681E-05 0.681E-05 0.835E-05 0.102E-04
> 0.103E-04 0.106E-04
>
> ERROR: BASIS LINEARLY DEPENDENT OR WRONG S
>
> The current threshold for the smallest tolerated eigenvalue of S is 1.0D-08
> ------------------------
>
> I spent some time searching the user-list and google, and I couldn't
> find a solution I needed. What I gather is that since my geometries
> are reasonable, the jobs run fine with smaller basis sets, and since I
> only have one eigenvalue below the threshold, this error probably just
> comes from the fact my QZ basis is (not surprisingly) linearly
> dependent and should be scaled down one or a few basis functions.
>
> Pardon the naive question, but does someone know the "right" way to
> modify this basis set so it is not linearly dependent, but that it is
> also of desired quality for a meaningful CBS extrapolation? Is it as
> straightforward as starting at the QZ basis and removing individual
> diffuse functions until the basis set is not linearly dependent, or is
> there a higher art to it?
>
> My input file is:
> ***Pd_4Cl
>
> cartesian
>
> geomtyp=xyz
> orient,mass
> geometry={
> 5
> Pd_4Cl calc
> Pd -0.213992 0.027380 -0.117721
> Cl -1.943705 -1.515957 0.437439
> Cl 1.304488 -1.784426 -0.423439
> Cl 1.515859 1.570547 -0.672946
> Cl -1.732472 1.839247 0.188003
> }
> set,charge=-2
>
> <BASIS>
>
> hf
>
> In my <BASIS> block, I've been using a script to pasting the explicit
> basis sets from the EMSL website into the input file. I don't want to
> paste all those messy numbers here unless that's actually requested.
>
> Thanks in advance!
>
> --john
>
> #####################################################################
> Ph D. John A. Keith
> Institut für Elektrochemie Email: john.keith at uni-ulm.de
> Universität Ulm Telephone: +49-731-50-26061
> Albert-Einstein-Allee 47 URL: http://www.wag.caltech.edu/home/johnk
> Ulm D-89069 URL:
> http://www.uni-ulm.de/theo_echem/index.shtml
> #####################################################################
> There is as much difference between us and ourselves as between us and others.
>
> -Michel de Montaigne, essayist (1533-1592)
> #####################################################################
> _______________________________________________
> Molpro-user mailing list
> Molpro-user at molpro.net
> http://www.molpro.net/mailman/listinfo/molpro-user
More information about the Molpro-user
mailing list