COSMO-SMD Initialization (regarding solvent accessible surface area)

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
Clicked A Few Times
Threads 2
Posts 9
Hello,

I am using SMD solvation on a CdTe cluster model (input selections are provided below) and am getting odd results regarding the G(SMD-CDS) energy contribution. On a clean CdTe cluster model, after the cosmo initialization section NWCHEM reports: G(SMD-CDS) = 7.190 kcal/mol and a SMS-CDS SASA = 55.418 (Input and Output #1 below), whereas if I have an adsorbed hydrogen atom on the same CdTe cluster model nwchem reports: G(SMD-CDS) = 0.000 kcal/mol and SMD-CDS SASA = 0.000 ang^2 (Input and Output #2 below). In the second case, I would expect SMD-CDS SASA (and the corresponding CDS energy contribution) to be non-zero given the similarity to the geometry employed in the first case. Could anyone help me understand why I am getting values of 0.000?

Thanks in advance,
Tom


Input #1:

start "CdTe"
memory 4000 Mb

geometry noautoz
H1 -1.79253420 -8.58929717 -1.75764609 charge 1.5
H1 0.62359575 -9.68359177 1.76906492 charge 1.5
H1 2.88179988 -8.28778835 -1.75733531 charge 1.5
H1 8.33542745 2.74270851 -1.75660404 charge 1.5
H1 -8.69679715 4.29916838 1.76771176 charge 1.5
H1 -6.54114935 5.84744645 -1.75820071 charge 1.5
H1 5.73643514 6.63986483 -1.75641179 charge 1.5
H1 8.07237021 5.38170371 1.76951997 charge 1.5
H1 -8.62248184 1.64242566 -1.75420551 charge 1.5
H2 -5.68704512 -6.61020021 -1.50237666 charge 0.5
H2 -4.31353070 -5.68356312 -3.90375747 charge 0.5
H2 5.00989173 -5.08224574 -3.90413763 charge 0.5
H2 0.34817411 -5.38290614 -3.87144755 charge 0.5
H2 6.49018134 -5.82416961 -1.50110954 charge 0.5
H2 -6.90597350 -1.79586283 -3.90331241 charge 0.5
H2 -8.28878582 -2.70706910 -1.50243461 charge 0.5
H2 2.41765728 -1.19440769 -3.90269553 charge 0.5
H2 -0.17278133 2.69162039 -3.90300215 charge 0.5
H2 -4.83449968 2.39095948 -3.87081106 charge 0.5
H2 -2.24308513 -1.49562491 -3.90400427 charge 0.5
H2 8.56880507 -1.62003953 -1.50085640 charge 0.5
H2 4.48892220 2.99227668 -3.86869324 charge 0.5
H2 7.08054637 -0.89416926 -3.90288738 charge 0.5
H2 -0.46201485 7.14916943 1.79274252 charge 0.5
H2 -2.76487856 6.57924938 -3.90256003 charge 0.5
H2 -2.88093598 8.23444474 -1.50121972 charge 0.5
H2 6.42274467 -3.17638343 1.79355763 charge 0.5
H2 1.80038850 8.53308269 -1.50215559 charge 0.5
H2 1.89787214 6.87969593 -3.90399890 charge 0.5
H2 -5.96372715 -3.97543323 1.79249501 charge 0.5
Cd -1.92309071 -6.84994572 -1.14912003
Cd 2.78841226 -6.54538404 -1.14976757
Cd -1.83772534 -4.27067100 1.53420074
Cd 0.50378575 -7.80067193 1.63900756
Cd 2.36990564 -3.99945102 1.53339230
Cd -4.46649284 -2.97935249 -1.17201824
Cd -7.06379121 0.85665246 -1.14859236
Cd -2.37699625 1.17513033 -1.15984818
Cd -7.00498473 3.46359896 1.63755684
Cd 0.17062142 -2.64500936 -1.15968290
Cd -4.64749369 -0.05544230 1.53370218
Cd 4.81410015 -2.37917044 -1.17191515
Cd 2.20562418 1.47085427 -1.15881492
Cd -2.77976745 3.72492868 1.53326461
Cd 6.89469647 1.75980268 -1.14630933
Cd 2.27652012 4.05189514 1.53473938
Cd 4.61650778 0.54238791 1.53422665
Cd 6.50101475 4.33456140 1.64000910
Cd -4.96914219 5.09092506 -1.14956936
Cd -0.34607553 5.35918101 -1.17176903
Cd 4.27424650 5.68864051 -1.14698944
Te -2.13733786 -7.01561718 1.68330058
Te -4.30350670 -5.66021492 -2.12483856
Te 3.02182962 -6.68257723 1.68291100
Te 0.34677775 -5.37006471 -2.06373874
Te 4.99660737 -5.05934947 -2.12489649
Te -6.87861876 -1.79744088 -2.12460033
Te -7.29635485 0.72385828 1.68417000
Te -4.44015729 -2.96112363 1.66745286
Te 2.41519477 -1.19358188 -2.13154061
Te -2.11943569 1.04724002 1.66605338
Te -2.24129196 -1.49407411 -2.13255858
Te 1.96768891 1.31108834 1.66710745
Te -4.82542064 2.38485179 -2.06317783
Te 4.78213244 -2.36585839 1.66765661
Te -0.17368693 2.68871227 -2.13162044
Te 0.15175184 -2.36042211 1.66647278
Te 7.14364266 1.65537162 1.68580106
Te 7.05439855 -0.89654157 -2.12360734
Te 4.47794115 2.98592428 -2.06097300
Te -2.74926151 6.56038189 -2.12347509
Te -0.34427436 5.32325374 1.66778582
Te -5.00545069 5.35883773 1.68299528
Te 1.88406299 6.85679987 -2.12476666
Te 4.27482702 5.95765404 1.68544347
end

constraints
fix atom 1:30
end

charge 0

basis "ao basis"
* library 6-31G** except Cd Te
Cd library "Stuttgart RSC 1997 ECP"
Te library "Stuttgart RLC ECP"
end

basis "cd basis"
* library "Ahlrichs Coulomb Fitting"
end

ecp
Cd library "Stuttgart RSC 1997 ECP"
Te library "Stuttgart RLC ECP"
end

driver
xyz output
maxiter 200
end

dft
iterations 200
xc b3lyp
disp vdw 2
direct
end

cosmo
do_cosmo_smd true
solvent water
do_gasphase false
end

task dft optimize


Output #1:

...... end of -cosmo- initialization ......


G(SMD-CDS) energy (kcal/mol)  =     7.190
SMD-CDS SASA (angstrom**2) = 55.418


Input #2:

start "GaP111"

memory 4000 Mb

geometry noautoz
H1 5.59219903 6.77947913 -1.75129237 charge 1.5
H1 3.95626664 8.87397246 1.77169032 charge 1.5
H1 1.31285496 8.68409788 -1.75773952 charge 1.5
H1 -8.64079336 1.45034625 -1.76470281 charge 1.5
H1 5.71365542 -7.83862760 1.78515175 charge 1.5
H1 3.09008260 -8.21076468 -1.74423068 charge 1.5
H1 -8.15021874 -3.20845549 -1.76045484 charge 1.5
H1 -9.63896044 -1.00645844 1.76193638 charge 1.5
H1 6.88658600 -5.45403075 -1.73703023 charge 1.5
H2 8.12096396 3.21805987 -1.49016892 charge 0.5
H2 6.47742799 3.03371307 -3.89348258 charge 0.5
H2 -2.05826599 6.83313184 -3.90737968 charge 0.5
H2 2.20949171 4.93320190 -3.86793053 charge 0.5
H2 -3.02751139 8.17951730 -1.50691773 charge 0.5
H2 6.96684741 -1.61305009 -3.88924288 charge 0.5
H2 8.61155857 -1.44678510 -1.48592311 charge 0.5
H2 -1.56905250 2.18591135 -3.90213676 charge 0.5
H2 -1.08038413 -2.45887024 -3.89839738 charge 0.5
H2 3.18757870 -4.35834006 -3.85944844 charge 0.5
H2 2.69819367 0.28729591 -3.89669082 charge 0.5
H2 -6.82107449 5.42200839 -1.50961706 charge 0.5
H2 -5.34832238 -0.55937957 -3.87084633 charge 0.5
H2 -5.83771060 4.08625957 -3.90908570 charge 0.5
H2 -2.90245385 -6.53630778 1.79805701 charge 0.5
H2 -0.59096430 -7.10563249 -3.89415667 charge 0.5
H2 -1.26045647 -8.62355251 -1.49243219 charge 0.5
H2 -4.20209965 5.80597744 1.78825846 charge 0.5
H2 -5.54449554 -6.71352699 -1.50037822 charge 0.5
H2 -4.85982687 -5.20574214 -3.90210539 charge 0.5
H2 7.13756881 0.75887109 1.80546430 charge 0.5
Cd 4.89733876 5.18226392 -1.14830336
Cd 0.58165091 7.10659610 -1.15065536
Cd 3.64090676 2.99643277 1.53173528
Cd 3.17570117 7.16051346 1.63664155
Cd -0.30293522 4.76909426 1.48772985
Cd 5.33796135 0.56904956 -1.18254196
Cd 5.86354326 -4.04225410 -1.13469655
Cd 1.55563614 -2.21976660 -1.14207436
Cd 4.54744360 -6.35277994 1.65211884
Cd 1.09839994 2.41997288 -1.17294219
Cd 4.22307288 -2.01371877 1.52870155
Cd -3.15209102 4.34736530 -1.17347712
Cd -2.62262803 -0.27444372 -1.16587719
Cd 0.64839967 -4.76360713 1.14061379
Cd -6.90838626 1.64724286 -1.15791632
Cd -3.92586579 -2.51383145 1.52830758
Cd -4.34230851 1.64990803 1.52765541
Cd -7.77196779 -0.79127710 1.63154216
Cd 2.04259194 -6.86709630 -1.09926017
Cd -2.19586282 -4.90723835 -1.11720413
Cd -6.41919629 -3.05307646 -1.14355447
Te 5.18288399 5.27530875 1.68549337
Te 6.45262035 3.01983235 -2.11533897
Te 0.43816738 7.43047580 1.67511932
Te 2.19957237 4.92519961 -2.05914414
Te -2.06154234 6.80400048 -2.12640295
Te 6.93749031 -1.60084051 -2.10723120
Te 6.07805815 -4.04446562 1.69998628
Te 5.32748999 0.63058749 1.66874990
Te -1.56882830 2.18658856 -2.12950725
Te 1.46163028 -1.76276665 1.69201556
Te 2.67227771 0.24741647 -2.11600672
Te -2.34575058 -0.24423355 1.66371016
Te 3.17965729 -4.39039356 -2.04850011
Te -3.15253066 4.31971263 1.66866531
Te -1.09547280 -2.47204092 -2.11944886
Te 0.97748052 2.29707816 1.67077743
Te -7.08853676 1.87822297 1.67473792
Te -5.81847598 4.07613221 -2.13185906
Te -5.33922059 -0.56224569 -2.06238808
Te -0.60690565 -7.07036227 -2.10614102
Te -2.23752662 -4.84765807 1.72346758
Te 1.96295260 -7.24762679 1.71979865
Te -4.84532947 -5.19045147 -2.12461416
Te -6.57688555 -3.26513838 1.68804520
H 1.49447312 -0.00669412 1.68893084
end

constraints
fix atom 1:30
end

charge 0

basis "ao basis"
* library 6-31G** except Cd Te
Cd library "Stuttgart RSC 1997 ECP"
Te library "Stuttgart RLC ECP"
end

basis "cd basis"
* library "Ahlrichs Coulomb Fitting"
end

ecp
Cd library "Stuttgart RSC 1997 ECP"
Te library "Stuttgart RLC ECP"
end

driver
xyz output
maxiter 200
end

dft
iterations 200
xc b3lyp
disp vdw 2
direct
mult 2
end

cosmo
do_cosmo_smd true
solvent water
do_gasphase false
end

task dft optimize

Output #2:

...... end of -cosmo- initialization ......


G(SMD-CDS) energy (kcal/mol)  =     0.000
SMD-CDS SASA (angstrom**2) = 0.000


Who's here now Members 0 Guests 0 Bots/Crawler 0


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC