An all‐electron numerical method for solving the local density functional for polyatomic molecules
TL;DR: In this paper, a method for accurate and efficient local density functional calculations (LDF) on molecules is described and presented with results using fast convergent threedimensional numerical integrations to calculate the matrix elements occurring in the Ritz variation method.
Abstract: A method for accurate and efficient local density functional calculations (LDF) on molecules is described and presented with results The method, Dmol for short, uses fast convergent three‐dimensional numerical integrations to calculate the matrix elements occurring in the Ritz variation method The flexibility of the integration technique opens the way to use the most efficient variational basis sets A practical choice of numerical basis sets is shown with a built‐in capability to reach the LDF dissociation limit exactly Dmol includes also an efficient, exact approach for calculating the electrostatic potential Results on small molecules illustrate present accuracy and error properties of the method Computational effort for this method grows to leading order with the cube of the molecule size Except for the solution of an algebraic eigenvalue problem the method can be refined to quadratic growth for large molecules
Citations
More filters
TL;DR: In this paper, the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids is described and the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets.
Abstract: Recent extensions of the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids are described. Furthermore the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets. The method is applied to calculations of the enthalpy of formation of a set of molecules and solids. We find that the present numerical localized basis sets yield improved results as compared to previous results for the same functionals. Enthalpies for the formation of H, N, O, F, Cl, and C, Si, S atoms from the thermodynamic reference states are calculated at the same level of theory. It is found that the performance in predicting molecular enthalpies of formation is markedly improved for the Perdew–Burke–Ernzerhof [Phys. Rev. Lett. 77, 3865 (1996)] functional.
8,496 citations
TL;DR: In this paper, a general all-atom force field for atomistic simulation of common organic molecules, inorganic small molecules, and polymers was developed using state-of-the-art ab initio and empirical parametrization techniques.
Abstract: A general all-atom force field for atomistic simulation of common organic molecules, inorganic small molecules, and polymers was developed using state-of-the-art ab initio and empirical parametrization techniques. The valence parameters and atomic partial charges were derived by fitting to ab initio data, and the van der Waals (vdW) parameters were derived by conducting MD simulations of molecular liquids and fitting the simulated cohesive energies and equilibrium densities to experimental data. The combined parametrization procedure significantly improves the quality of a general force field. Validation studies based on large number of isolated molecules, molecular liquids and molecular crystals, representing 28 molecular classes, show that the present force field enables accurate and simultaneous prediction of structural, conformational, vibrational, and thermophysical properties for a broad range of molecules in isolation and in condensed phases. Detailed results of the parametrization and validation f...
4,722 citations
TL;DR: It is shown how derivatives of the GPW energy functional, namely ionic forces and the Kohn–Sham matrix, can be computed in a consistent way and the computational cost is scaling linearly with the system size, even for condensed phase systems of just a few tens of atoms.
4,047 citations
TL;DR: The emphasis of this review is on the origin of the electrocatalytic activity of nanostructured catalysts toward a series of key clean energy conversion reactions by correlating the apparent electrode performance with their intrinsic electrochemical properties.
Abstract: A fundamental change has been achieved in understanding surface electrochemistry due to the profound knowledge of the nature of electrocatalytic processes accumulated over the past several decades and to the recent technological advances in spectroscopy and high resolution imaging. Nowadays one can preferably design electrocatalysts based on the deep theoretical knowledge of electronic structures, via computer-guided engineering of the surface and (electro)chemical properties of materials, followed by the synthesis of practical materials with high performance for specific reactions. This review provides insights into both theoretical and experimental electrochemistry toward a better understanding of a series of key clean energy conversion reactions including oxygen reduction reaction (ORR), oxygen evolution reaction (OER), and hydrogen evolution reaction (HER). The emphasis of this review is on the origin of the electrocatalytic activity of nanostructured catalysts toward the aforementioned reactions by correlating the apparent electrode performance with their intrinsic electrochemical properties. Also, a rational design of electrocatalysts is proposed starting from the most fundamental aspects of the electronic structure engineering to a more practical level of nanotechnological fabrication.
3,918 citations
TL;DR: This work synthesized uniform anatase TiO2 single crystals with a high percentage (47 per cent) of {001} facets using hydrofluoric acid as a morphology controlling agent and demonstrates that for fluorine-terminated surfaces this relative stability is reversed.
Abstract: [Yang, Hua Gui; Sun, Cheng Hua; Qiao, Shi Zhang; Liu, Gang; Smith, Sean Campbell; Lu, Gao Qing] Univ Queensland, ARC Ctr Excellence Funct Nanomat, Sch Engn, Brisbane, Qld 4072, Australia. [Yang, Hua Gui; Sun, Cheng Hua; Qiao, Shi Zhang; Liu, Gang; Smith, Sean Campbell; Lu, Gao Qing] Univ Queensland, Australian Inst Bioengn & Nanotechnol, Ctr Computat Mol Sci, Brisbane, Qld 4072, Australia. [Zou, Jin] Univ Queensland, Ctr Microscopy & Microanal, Brisbane, Qld 4072, Australia. [Zou, Jin] Univ Queensland, Sch Engn, Brisbane, Qld 4072, Australia. [Liu, Gang; Cheng, Hui Ming] Chinese Acad Sci, Met Res Inst, Shenyang Natl Lab Mat sci, Shenyang 110016, Peoples R China.;Lu, GQ (reprint author), Univ Queensland, ARC Ctr Excellence Funct Nanomat, Sch Engn, Brisbane, Qld 4072, Australia;s.qiao@uq.edu.au maxlu@uq.edu.au
3,656 citations
Cites background from "An all‐electron numerical method fo..."
...pseudopotential for core electrons: Effective core potential [15] basis set: Double-numeric quality basis set [16,17]...
[...]
...[16] Delley, B....
[...]
References
More filters
TL;DR: In this paper, the Hartree and Hartree-Fock equations are applied to a uniform electron gas, where the exchange and correlation portions of the chemical potential of the gas are used as additional effective potentials.
Abstract: From a theory of Hohenberg and Kohn, approximation methods for treating an inhomogeneous system of interacting electrons are developed. These methods are exact for systems of slowly varying or high density. For the ground state, they lead to self-consistent equations analogous to the Hartree and Hartree-Fock equations, respectively. In these equations the exchange and correlation portions of the chemical potential of a uniform electron gas appear as additional effective potentials. (The exchange portion of our effective potential differs from that due to Slater by a factor of $\frac{2}{3}$.) Electronic systems at finite temperatures and in magnetic fields are also treated by similar methods. An appendix deals with a further correction for systems with short-wavelength density oscillations.
47,477 citations
TL;DR: In this article, the ground state of an interacting electron gas in an external potential was investigated and it was proved that there exists a universal functional of the density, called F[n(mathrm{r})], independent of the potential of the electron gas.
Abstract: This paper deals with the ground state of an interacting electron gas in an external potential $v(\mathrm{r})$. It is proved that there exists a universal functional of the density, $F[n(\mathrm{r})]$, independent of $v(\mathrm{r})$, such that the expression $E\ensuremath{\equiv}\ensuremath{\int}v(\mathrm{r})n(\mathrm{r})d\mathrm{r}+F[n(\mathrm{r})]$ has as its minimum value the correct ground-state energy associated with $v(\mathrm{r})$. The functional $F[n(\mathrm{r})]$ is then discussed for two situations: (1) $n(\mathrm{r})={n}_{0}+\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{n}(\mathrm{r})$, $\frac{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{n}}{{n}_{0}}\ensuremath{\ll}1$, and (2) $n(\mathrm{r})=\ensuremath{\phi}(\frac{\mathrm{r}}{{r}_{0}})$ with $\ensuremath{\phi}$ arbitrary and ${r}_{0}\ensuremath{\rightarrow}\ensuremath{\infty}$. In both cases $F$ can be expressed entirely in terms of the correlation energy and linear and higher order electronic polarizabilities of a uniform electron gas. This approach also sheds some light on generalized Thomas-Fermi methods and their limitations. Some new extensions of these methods are presented.
38,160 citations
TL;DR: In this article, a general and natural choice is to share the charge density at each point among the several atoms in proportion to their free-atom densities at the corresponding distances from the nuclei.
Abstract: For quantitative description of a molecular charge distribution it is convenient to dissect the molecule into well-defined atomic fragments. A general and natural choice is to share the charge density at each point among the several atoms in proportion to their free-atom densities at the corresponding distances from the nuclei. This prescription yields well-localized bonded-atom distributions each of which closely resembles the molecular density in its vicinity. Integration of the atomic deformation densities — bonded minus free atoms — defines net atomic charges and multipole moments which concisely summarize the molecular charge reorganization. They permit calculation of the external electrostatic potential and of the interaction energy between molecules or between parts of the same molecule. Sample results for several molecules indicate a high transferability of net atomic charges and moments.
5,234 citations
4,691 citations