Fast computation of DFT nuclear gradient with multiresolution*This article is dedicated to Dr. Russell J. Boyd for his distinguished career in chemical research.

2011 ◽  
Vol 89 (6) ◽  
pp. 657-662 ◽  
Author(s):  
Nicholas J. Russ ◽  
Chun-min Chang ◽  
Jing Kong

We present an efficient algorithm for evaluating the exchange-correlation contribution to the nuclear gradients of density-functional theory calculation within the local spin-density approximation. The algorithm is an extension of the multiresolution exchange-correlation (mrXC) method, which treats smooth and compact parts of the electron density separately. The nuclear gradient of the smooth density is calculated on the even-spaced grid while the compact part of the density is handled on the normal atom-centered grid (ACG). The overall formulism is still formally based on the ACG, and thus does not change the results of the existing ACG-based algorithms for all-electron density-functional theory (DFT) calculations. The variation of the positions and weights of ACG owing to the nuclear perturbation is also handled rigorously. Benchmark calculations with different basis sets and sizes of ACG show that mrXC reduces the computational cost by several times without loss of accuracy. It also lessens the impact on the CPU time when the size of the ACG is increased.

2001 ◽  
Vol 56 (1) ◽  
pp. 13-24 ◽  
Author(s):  
Hans Bock ◽  
Sven Holl ◽  
Volker Krenzel

Abstract The structures of tri-and tetraiodo-substituted carbon compounds are determined either expe­rimentally by X-Ray Structure Analysis or, because crystallization of tetraiodothiophene could not be achieved, approximated by Density Functional Theory optimization of structural data from a donor/acceptor complex. The structures show noteworthy details such as a second po­lymorph of tetraiodoethene crystallized by sublimation or herringbone crystal packing patterns of tetraiodopyrrole derivatives. All molecular geometries are discussed and compared based on relativistic density functional theory calculations with 6 -31G* basis sets including iodine pseudopotentials. They reproduce even finer structural details due to van der Waals repulsion of the bulky iodo substituents. Natural Bond Orbital (NBO) charge distributions suggest positive partial charges at all iodine centers with the strongest polarization Cδ㊀ → Iδ㊉ in HCI3, which contains well over 97% iodine.


2020 ◽  
Author(s):  
Peter Kraus

Improving results of correlated wavefunction theory calculations by extrapolating from successive basis sets is nowadays a common practice. However, such approaches are uncommon in density functional theory, especially due its faster convergence towards the basis set limit. In this work I present approaches for basis set extrapolation in density functional theory that enable users to obtain results of 4-zeta quality from 3- and 2-zeta calculations, i.e. at a fraction of the computational cost. The extrapolation techniques work well with modern density functionals and common basis sets.<br>


2020 ◽  
Vol 860 ◽  
pp. 282-287
Author(s):  
Wan Nurfadhilah Zaharim ◽  
Shukri Sulaiman ◽  
Saidah Sakinah Mohd Tajudin ◽  
Siti Nuramira Abu Bakar ◽  
Nur Eliana Ismail ◽  
...  

The Density Functional Theory method was employed to investigate the electronic structure and muonium hyperfine interaction of muonium trapped near carbon atom labelled as '5' in cytosine nucleobase. Eighteen different basis sets in combination with B3LYP functional were examined in geometry optimization calculations on the muoniated radical. There are significant quantitative differences in the calculated total energy. The employment of basis set that does not include polarization function produces an optimized structure with high total energy. The 6-311++G(d,p) basis set yielded the lowest total energy as compared to other basis sets. The bond order of muonium trapped at C5 atom is in the range of 0.841 to 0.862. The 6-31G basis set produced the muonium Fermi contact coupling constant that is the closest to the experimental value.


2020 ◽  
Author(s):  
Peter Kraus

Improving results of correlated wavefunction theory calculations by extrapolating from successive basis sets is nowadays a common practice. However, such approaches are uncommon in density functional theory, especially due its faster convergence towards the basis set limit. In this work I present approaches for basis set extrapolation in density functional theory that enable users to obtain results of 4-zeta quality from 3- and 2-zeta calculations, i.e. at a fraction of the computational cost. The extrapolation techniques work well with modern density functionals and common basis sets.<br>


2013 ◽  
Vol 16 ◽  
pp. 117-122
Author(s):  
Emerson Rengifo ◽  
Gustavo Murillo

Calculations of the electronic structure of a stacked dimmer sequence from the D(GCAAACGTTTGC)2 B-DNA dodecamer resolved in a PDB file 1HQ7 are performed within density functional theory. Seeking to understand the minimum level of theory that yields a reliable description for these systems, the basis sets 6-31g*, 6-31g*+BSSE, 6-311g*, 6-311g**, 6-311++g** along with the B3LYP and PBE0 exchange-correlation functionals were employed. These results are then used to implement a one dimensional model of long stacked systems to obtain a new semiempirical method that can be employed at low computational cost.


2007 ◽  
Vol 06 (04) ◽  
pp. 675-685 ◽  
Author(s):  
XIAO-HONG LI ◽  
RUI-ZHOU ZHANG ◽  
XIAN-ZHOU ZHANG ◽  
XIN-LU CHENG ◽  
XIANG-DONG YANG

The heats of formation (HOFs) for 15 thiol compounds are calculated by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86) methods with 6-311G** basis set and the complete basis set (CBS-Q) method. It is demonstrated that the B3P86 and CBS-Q methods are accurate to compute the reliable HOFs for thiol compounds. In order to test whether the B3P86 method has a low basis set sensitivity, the HOFs for six thiol compounds are also calculated by using the B3P86 method with 6-31+G*, 6-31+G**, and 6-311+G** basis sets for comparison. We also extend our study by employing the nonlocal BLYP method together with 6-31+G* basis set to calculate the HOFs for thiol compounds. The obtained results are compared with the experimental results. It is noted that the B3P86 method is not sensitive to the basis sets. Considering the inevitably computational cost of CBS-Q method and the reliability of the B3P86 calculation, the B3P86 method with a moderate or a larger basis set such as 6-311G** and 6-311+G** may be more suitable to calculate the HOFs of thiol compounds. In addition, we believe that the maximum error associated with the calculated HOFs is less than 6 kcal/mol for the B3P86/6-311G** method and it is expected that the error bar is more likely 1–5 kcal/mol for the HOFs of thiol compounds.


2019 ◽  
Author(s):  
S. Giarrusso ◽  
Paola Gori-Giorgi

We analyze in depth two widely used definitions (from the theory of conditional probablity amplitudes and from the adiabatic connection formalism) of the exchange-correlation energy density and of the response potential of Kohn-Sham density functional theory. We introduce a local form of the coupling-constant-dependent Hohenberg-Kohn functional, showing that the difference between the two definitions is due to a corresponding local first-order term in the coupling constant, which disappears globally (when integrated over all space), but not locally. We also design an analytic representation for the response potential in the strong-coupling limit of density functional theory for a model single stretched bond.<br>


2019 ◽  
Author(s):  
Mark Iron ◽  
Trevor Janes

A new database of transition metal reaction barrier heights – MOBH35 – is presented. Benchmark energies (forward and reverse barriers and reaction energy) are calculated using DLPNO-CCSD(T) extrapolated to the complete basis set limit using a Weizmann1-like scheme. Using these benchmark energies, the performance of a wide selection of density functional theory (DFT) exchange–correlation functionals, including the latest from the Truhlar and Head-Gordon groups, is evaluated. It was found, using the def2-TZVPP basis set, that the ωB97M-V (MAD 1.8 kcal/mol), ωB97X-V (MAD 2.1 kcal/mol) and SCAN0 (MAD 2.1 kcal/mol) hybrid functionals are recommended. The double-hybrid functionals PWPB95 (MAD 1.6 kcal/mol) and B2K-PLYP (MAD 1.8 kcal/mol) did perform slightly better but this has to be balanced by their increased computational cost.


Sign in / Sign up

Export Citation Format

Share Document