1
|
On point perforating defects in bilayer structures. Phys Chem Chem Phys 2023; 25:30477-30487. [PMID: 37921407 DOI: 10.1039/d3cp03719c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
This article deals with the issue of perforating point defects (pores) in a bilayer heterostructure composed of striped borophene and graphene. Three types of non-equivalent vacancies of the minimum size are considered. These include a single vacancy and two double vacancies. The study of the properties and stability of the perforating defects in borophene-graphene heterostructures is important given the increasing role of such structures in membranes for water purification, renewable energy generation, and other osmotic applications. Using the DFT method, the atomic configurations and main energy characteristics of the proposed defects are obtained. The results show that the formation of a single boron vacancy on the borophene side of borophene-graphene requires less energy than the formation of a carbon vacancy in graphene. Comparisons between double vacancies in nanoscale materials are unreliable because different reference systems produce the different chemical potentials. The problem of choosing the reference system for reliable calculation of the vacancy formation energies is posed and discussed. Using borophene-graphene as an example, it is shown that the reference system strongly affects the magnitude and sign of the vacancy formation energy. Hydrogenation is tested to stabilize the proposed defects.
Collapse
|
2
|
Semiempirical Approach to the Fukui Function Analysis of Uric Acid under Different pH Conditions. J Phys Chem A 2023; 127:8228-8237. [PMID: 37751600 DOI: 10.1021/acs.jpca.3c04695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
Analytic Fukui functions calculated at a first-principles level are combined with experimental pKa values and the calculation of tautomerization energies to obtain the effective regioselectivity of uric acid toward electron-transfer reactions under different pH conditions. Second-order electron binding energies are also computed to determine which of the tautomers is more likely to participate in the electron transfer. A comparison of vertical and adiabatic proton detachment energies allows us to conclude that tautomerization is not mediating deprotonation and that two monoanionic species are of comparable relevance. The main difference between these monoanionic species is the ring that has been deprotonated. Both monoanionic species are produced from a single neutral tautomer and mainly produce a single dianionic tautomer. As a method for the analysis of systems affected by pH such as uric acid, we propose to plot condensed Fukui functions versus pH, allowing us to draw the effect of pH on the regioselectivity of electron transfer in a single image.
Collapse
|
3
|
On the energetic and magnetic stability of neutral and charged lithium clusters doped with one and two yttrium atoms. Phys Chem Chem Phys 2023; 25:9656-9668. [PMID: 36943209 DOI: 10.1039/d3cp00128h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2023]
Abstract
DFT calculations were performed to study the effect on energetic and magnetic stability when clusters with up to 24 lithium atoms were doped with one and two atoms of yttrium. In this, the effect of the charge was considered. As a result, some stable structures were identified as possible magnetic superatoms, among them, the YLi12+ cluster with an icosahedron geometry with a spin magnetic moment of 4 bohr magnetons. The participation of yttrium in the electron density of the unpaired electrons providing magnetism in clusters was corroborated at the level of a density of states (DOS) calculation and a spin density calculation. In particular, in the Y2Li12+ superatom, it was found that the encapsulated yttrium atom participates with 35.02% and the second yttrium atom with 15.04%. These percentages, with a contribution from p orbitals, but to a greater extent by d orbitals. The complementation to these percentages is due to the participation of the s and p orbitals of the lithium atoms. In general, doping with a second yttrium atom allowed to obtain a greater amount of high magnetic moments, and considering charged clusters allowed to obtain also high magnetic moments.
Collapse
|
4
|
Fitting Multiplet Simulations to L-Edge XAS Spectra of Transition-Metal Complexes Using an Adaptive Grid Algorithm. Inorg Chem 2023; 62:3738-3760. [PMID: 36808900 DOI: 10.1021/acs.inorgchem.2c02830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2023]
Abstract
A new methodology based on an adaptive grid algorithm followed by an analysis of the ground state from the fit parameters is presented to analyze and interpret experimental XAS L2,3-edge data. The fitting method is tested first in a series of multiplet calculations for d0-d7 systems and for which the solution is known. In most cases, the algorithm is able to find the solution, except for a mixed-spin Co2+ Oh complex, where it instead revealed a correlation between the crystal field and the electron repulsion parameters near spin-crossover transition points. Furthermore, the results for the fitting of previously published experimental data sets on CaO, CaF2, MnO, LiMnO2, and Mn2O3 are presented and their solution discussed. The presented methodology has allowed the evaluation of the Jahn-Teller distortion in LiMnO2, which is consistent with the observed implications in the development of batteries, which use this material. Moreover, a follow-up analysis of the ground state in Mn2O3 has demonstrated an unusual ground state for the highly distorted site which would be impossible to optimize in a perfect octahedral environment. Ultimately, the presented methodology can be used in the analysis of X-ray absorption spectroscopy data measured at the L2,3-edge for a large number of materials and molecular complexes of first-row transition metals and can be expanded to the analysis of other X-ray spectroscopic data in future studies.
Collapse
|
5
|
A Fukui Analysis of an Arginine-Modified Carbon Surface for the Electrochemical Sensing of Dopamine. MATERIALS (BASEL, SWITZERLAND) 2022; 15:6337. [PMID: 36143660 PMCID: PMC9506051 DOI: 10.3390/ma15186337] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 09/05/2022] [Accepted: 09/08/2022] [Indexed: 06/16/2023]
Abstract
Amino acid-modified carbon interfaces have huge applications in developing electrochemical sensing applications. Earlier reports suggested that the amine group of amino acids acted as an oxidation center at the amino acid-modified electrode interface. It was interesting to locate the oxidation centers of amino acids in the presence of guanidine. In the present work, we modeled the arginine-modified carbon interface and utilized frontier molecular orbitals and analytical Fukui functions based on the first principle study computations to analyze arginine-modified CPE (AMCPE) at a molecular level. The frontier molecular orbital and analytical Fukui results suggest that the guanidine (oxidation) and carboxylic acid (reduction) groups of arginine act as additional electron transfer sites on the AMCPE surface. To support the theoretical observations, we prepared the arginine-modified CPE (AMCPE) for the cyclic voltammetric sensing of dopamine (DA). The AMCPE showed excellent performance in detecting DA in blood serum samples.
Collapse
|
6
|
Pre-post electron transfer regioselectivity at glycine modified graphene electrode interface for voltammetric sensing applications. JOURNAL OF ELECTROCHEMICAL SCIENCE AND ENGINEERING 2022. [DOI: 10.5599/jese.1438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
In the last few years, glycine (GL) showed good experimental evidence as an electron transfer (ET) mediator at the carbon (in particular graphene (GR)) interface. However, ET properties of GL modified GR interface are still not known completely. These can be achieved using density functional theory-based models. Modelling of modified carbon electrode interfaces is essential in electroanalytical chemistry to get insights into their electronic and redox properties. Here we have modelled glycine modified graphene interface to find out its interfacial redox ET properties. Conceptual density functional theory concepts like frontier molecular orbital (FMO) theory and analytical Fukui functions were utilized to predict the ET sites on the modified graphene surface. It is shown that at the glycine-modified graphene interface, amine groups act as additional oxidation sites and carboxylic acid groups as additional reduction sites. Therefore, glycine acts as an ET mediator at the graphene-based electrode interface. The obtained results are well supported by previously published experimental reports.
Collapse
|
7
|
Studies of Monoamine Neurotransmitters at Nanomolar Levels Using Carbon Material Electrodes: A Review. MATERIALS (BASEL, SWITZERLAND) 2022; 15:5782. [PMID: 36013918 PMCID: PMC9415512 DOI: 10.3390/ma15165782] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 08/14/2022] [Accepted: 08/16/2022] [Indexed: 06/15/2023]
Abstract
Neurotransmitters (NTs) with hydroxyl groups can now be identified electrochemically, utilizing a variety of electrodes and voltammetric techniques. In particular, in monoamine, the position of the hydroxyl groups might alter the sensing properties of a certain neurotransmitter. Numerous research studies using electrodes modified on their surfaces to better detect specific neurotransmitters when other interfering factors are present are reviewed to improve the precision of these measures. An investigation of the monoamine neurotransmitters at nanoscale using electrochemical methods is the primary goal of this review article. It will be used to determine which sort of electrode is ideal for this purpose. The use of carbon materials, such as graphite carbon fiber, carbon fiber micro-electrodes, glassy carbon, and 3D printed electrodes are only some of the electrodes with surface modifications that can be utilized for this purpose. Electrochemical methods for real-time detection and quantification of monoamine neurotransmitters in real samples at the nanomolar level are summarized in this paper.
Collapse
|
8
|
Electronic structure and reactivity indexes of cobalt clusters, both pure and mixed with NO and [Formula: see text] ([Formula: see text], [Formula: see text] and [Formula: see text]). J Mol Model 2022; 28:197. [PMID: 35729282 DOI: 10.1007/s00894-022-05165-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Accepted: 05/19/2022] [Indexed: 11/26/2022]
Abstract
Among the most popular motivations for environmental scientists is improving materials that could be useful to fight or avoid pollution. This work shows a study of neutral and cationic cobalt clusters from 4 to 9 atoms ([Formula: see text], q = 0,1 and n = 4-9) to model their separate interaction with contaminant nitric and nitrous oxides. This study is within the framework of the density functional theory in the Kohn-Sham scheme by using BPW91 functional and 6-311G and 6-31G* basis sets to calculate global and local reactivity indexes. The effect of spin multiplicity is also determined. Results on the geometries of pure cobalt clusters agree with previously reported structures. Global minimum energy structures showed a marked preference towards the interaction of nitric and nitrous oxide molecules with cobalt clusters through chemisorptive dissociation, with the dissociation of the corresponding nitrogen oxide. Reactivity indexes reveal an even-odd alternate, which is related to electron counts. Moreover, the chemical potential is lowering after interaction with nitrogen oxides. The Fukui function illustrates the reactive zones with a high probability of chemisorption of more nitrogen oxide molecules.
Collapse
|
9
|
Cis and trans platinum(II) N-heterocyclic carbene isomers: synthesis, characterization and biological activity. NEW J CHEM 2022. [DOI: 10.1039/d2nj02508f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The synthesis of cis and trans geometrical isomers of platinum(II) complexes with a symmetric N-heterocyclic carbene ligand (MeNHC) is reported. These complexes were obtained from 1,3-dimethylimidazolium-2-carboxylate, a masked NHC precursor....
Collapse
|
10
|
Analytic Alchemical Derivatives for the Analysis of Differential Acidity Assisted by the h Function. J Phys Chem A 2021; 125:10463-10474. [PMID: 34812636 DOI: 10.1021/acs.jpca.1c07364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Analytical calculation of alchemical derivatives based on auxiliary density perturbation theory is described, coded, and validated. For the case where the nucleus is a hydrogen atom and the nuclear charge is changed from 1 to 0, it turns out that a good estimate of the proton binding energies can be obtained very efficiently. First-order results correspond exactly to the molecular electrostatic potential evaluated at the hydrogen nucleus location (removing self-repulsion), in agreement with previously reported extensive studies. Therefore, the second-order results reported here are refinements in accuracy that finally allow a quantitative exploration of differential acidity. Furthermore, the recently reported h function is produced in its analytical form as a byproduct and local descriptor associated with the proton binding energy values found with this approach. In an example application, proton binding energies are computed for a family of imidazolium derivatives to demonstrate the capabilities and the stability of the method with respect to changes in basis set or exchange-correlation functional.
Collapse
|
11
|
Asymmetric Density Fitting with Modified Cholesky Decomposition Applied to Second-Order Electron Propagator. J Chem Theory Comput 2020; 16:1597-1605. [PMID: 31967819 DOI: 10.1021/acs.jctc.9b01215] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computation of molecular orbital electron repulsion integrals (MO-ERIs) as a transformation from atomic orbital ERIs (AO-ERIs) is the bottleneck of second-order electron propagator calculations when a single orbital is studied. In this contribution, asymmetric density fitting is combined with modified Cholesky decomposition to generate efficiently the required MO-ERIs. The key point of the presented algorithms is to keep track of integrals through partial contractions performed on three-center AO-ERIs; these contractions are stored in RAM instead of the AO-ERIs. Two implementations are provided, an in-core, which reduces the arithmetic and memory scaling factors as compared to the four-center AO-ERIs contraction method, and a semidirect, which overcomes memory limitations by evaluating antisymmetrized MO-ERIs in batches. On the numerical side, the proposed approach is fast and stable. The bad effects due to ill conditioning, namely, several negative and close to zero eigenvalues due to machine round off errors of the matrix associated with the density fitting process, are effectively controlled by means of a modified Cholesky factorization that avoids the matrix inversion needed to perform the asymmetrical density fitting implementation. The numerical experience presented shows that the in-core implementation is highly competitive to perform calculations on medium and large basis sets, while the semidirect implementation has small variations in time by changes in the available memory. The general applicability is illustrated on a set of selected relatively large-size molecules.
Collapse
|
12
|
Calculation of the EPR g-tensor from auxiliary density functional theory. J Chem Phys 2020; 152:014105. [DOI: 10.1063/1.5130174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Exploring the Structure, Energetic, and Magnetic Properties of Neutral Small Lithium Clusters Doped with Yttrium: Supermagnetic Atom Research. ACS OMEGA 2018; 3:11252-11261. [PMID: 31459233 PMCID: PMC6645110 DOI: 10.1021/acsomega.8b01463] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/26/2018] [Accepted: 09/04/2018] [Indexed: 06/10/2023]
Abstract
Density functional theory calculations based on magnetic and energetic stability criteria were performed to study a series of yttrium-doped lithium neutral clusters. A relativistic approximation was employed to properly describe the energy and multiplicity of the given clusters' fundamental states. The interaction of the 4d-Y atomic orbitals with the sp-Li states had an important role in the magnetic and energetic behavior of the selected systems. The spin density was concentrated over the yttrium atom regardless of the size of the cluster. Li7Y is a new stable superatom due to its enhanced magnetic properties.
Collapse
|
14
|
Quantum chemical and electrochemical studies of lysine modified carbon paste electrode surfaces for sensing dopamine. NEW J CHEM 2018. [DOI: 10.1039/c7nj04998f] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have improved the sensitivity of a carbon paste electrode from lysine for the sensitive detection of dopamine.
Collapse
|
15
|
Abstract
The active electron transfer (ET) sites on the graphene surface can be controlled by hexagonal boron nitride (h-BN) doping.
Collapse
|
16
|
Analytical Fukui and cyclic voltammetric studies on ferrocene modified carbon electrodes and effect of Triton X-100 by immobilization method. Electrochim Acta 2017. [DOI: 10.1016/j.electacta.2017.11.154] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
17
|
Abstract
The recently developed approach to auxiliary density perturbation theory (<em>J. Chem. Phys.</em> <strong>2008</strong>, <em>128</em>, 134105) for the purpose of calculating molecular properties is here extended to include open-shell systems. Both unrestricted and restricted formalisms are considered. A linear equation system, twice as large as the auxiliary function set is obtained in both cases. For the first time, the formulation for auxiliary density perturbation theory for restricted open-shell formalism is derived.
Collapse
|
18
|
|
19
|
Theoretical and cyclic voltammetric studies on electrocatalysis of benzethonium chloride at carbon paste electrode for detection of dopamine in presence of ascorbic acid. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.05.093] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
20
|
Theoretical calculation of polarizability isotope effects. J Mol Model 2017; 23:90. [DOI: 10.1007/s00894-017-3236-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2016] [Accepted: 01/13/2017] [Indexed: 11/29/2022]
|
21
|
Initial stage of the degradation of three common neonicotinoids: theoretical prediction of charge transfer sites. NEW J CHEM 2017. [DOI: 10.1039/c6nj02655a] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Tautomerization of acetamiprid gives alternatives to search new pathways for its degradation in water.
Collapse
|
22
|
Analytic second derivatives from auxiliary density perturbation theory. J Chem Phys 2016; 145:224103. [PMID: 27984884 DOI: 10.1063/1.4971292] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
|
24
|
Extended source model for diffusive coupling. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2016; 39:4. [PMID: 26802012 DOI: 10.1140/epje/i2016-16004-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Accepted: 12/23/2015] [Indexed: 06/05/2023]
Abstract
Motivated by the prevailing approach to diffusion coupling phenomena which considers point-like diffusing sources, we derived an analogous expression for the concentration rate of change of diffusively coupled extended containers. The proposed equation, together with expressions based on solutions to the diffusion equation, is intended to be applied to the numerical solution of systems exclusively composed of ordinary differential equations, however is able to account for effects due the finite size of the coupled sources.
Collapse
|
25
|
Calculation of positron binding energies using the generalized any particle propagator theory. J Chem Phys 2014; 141:114103. [DOI: 10.1063/1.4895043] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Symbolic Algebra Development for Higher-Order Electron Propagator Formulation and Implementation. J Chem Theory Comput 2014; 10:2363-70. [DOI: 10.1021/ct500204v] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
27
|
13C NMR and Fukui function analysis on C82hydroxylated fullerene through density functional theory. Mol Phys 2014. [DOI: 10.1080/00268976.2013.877168] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
28
|
Antiradical capacity of a series of organotin(IV) compounds: A chemical reactivity study in the Density Functional Theory framework. Inorganica Chim Acta 2014. [DOI: 10.1016/j.ica.2014.01.008] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
29
|
A generalized any-particle propagator theory: Prediction of proton affinities and acidity properties with the proton propagator. J Chem Phys 2013; 138:194108. [DOI: 10.1063/1.4805030] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
A generalized any particle propagator theory: Assessment of nuclear quantum effects on electron propagator calculations. J Chem Phys 2012; 137:074105. [DOI: 10.1063/1.4745076] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
|
32
|
|
33
|
|
34
|
Integral approximations in ab initio, electron propagator calculations. J Chem Phys 2009; 131:124110. [DOI: 10.1063/1.3238243] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Influence of Endohedral Confinement on the Electronic Interaction between He atoms: A He2@C20H20Case Study. Chemistry 2009; 15:1985-90. [DOI: 10.1002/chem.200801399] [Citation(s) in RCA: 101] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
36
|
|
37
|
|
38
|
|
39
|
Assessment of transition operator reference states in electron propagator calculations. J Chem Phys 2007; 127:134106. [DOI: 10.1063/1.2784638] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Abstract
A half-numeric algorithm for the evaluation of effective core potential integrals over Cartesian Gaussian functions is described. Local and semilocal integrals are separated into two-dimensional angular and one-dimensional radial integrals. The angular integrals are evaluated analytically using a general approach that has no limitation for the l-quantum number. The radial integrals are calculated by an adaptive one-dimensional numerical quadrature. For the semilocal radial part a pretabulation scheme is used. This pretabulation simplifies the handling of radial integrals, makes their calculation much faster, and allows their easy reuse for different integrals within a given shell combination. The implementation of this new algorithm is described and its performance is analyzed.
Collapse
|
41
|
Abstract
Density functional calculations of neutral and anionic niobium trimer monoxides are presented. The calculations were performed employing scalar quasirelativistic effective core potentials. In order to test the accuracy of the used effective core potentials in the framework of density functional theory the pulsed field ionization-zero electron kinetic energy photoelectron spectrum of Nb(3)O was simulated and compared to experiment. Different isomers of Nb(3)O and Nb(3)O(-) were studied in order to determine the ground state structures. For both neutral and anionic systems a planar C(2v) structure with an edge-bound oxygen atom was found as a ground state. Equilibrium structure parameters, harmonic frequencies, and adiabatic electron affinity are reported. The calculated electron affinity and frequencies are in good agreement with the available experimental data obtained recently from vibrationally resolved negative ion photoelectron spectroscopy.
Collapse
|
42
|
Efficient and reliable numerical integration of exchange-correlation energies and potentials. J Chem Phys 2004; 121:681-90. [PMID: 15260594 DOI: 10.1063/1.1759323] [Citation(s) in RCA: 133] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An adaptive numerical integrator for the exchange-correlation energy and potential is presented. It uses the diagonal elements of the exchange-correlation potential matrix as a grid generating function. The only input parameter is the requested grid tolerance. In combination with a defined cell function the adaptive grid generation scales almost linear with the number of basis functions in a system. With the adaptive numerical integrator the self-consistent field energy error, which is due to the numerical integration of the exchange-correlation energy, converges with increasing adaptive grid size to a reference value. The performance of the adaptive numerical integration is analyzed using molecules with first, second, and third row elements. Especially for transition metal systems the adaptive numerical integrator shows considerably improved performance and reliability.
Collapse
|