Accurate pKa Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pKa values for six carboxylic acids. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol to calculate pKa values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit.

Accurate relative pKa calculations for carboxylic acids using complete basis set and Gaussian-n models combined with continuum solvation methods

The complete basis set methods CBS-4, CBS-QB3, and CBS-APNO, and the Gaussian methods G2 and G3 were used to calculate the gas phase energy differences between six different carboxylic acids and their respective anions. Two different continuum methods, SM5.42R and CPCM, were used to calculate the free energy differences of solvation for the acids and their anions. Relative pKa values were calculated for each acid using one of the acids as a reference point. The CBS-QB3 and CBS-APNO gas phase calculations, combined with the CPCM/HF/6-31+G(d)//HF/6-31G(d) or CPCM/HF/6-31+G(d)//HF/6-31+G(d) continuum solvation calculations on the lowest energy gas phase conformer, and with the conformationally averaged values, give results accurate to ½ pKa unit. © 2001 American Institute of Physics.

Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting

Localized short-echo-time (1)H-MR spectra of human brain contain contributions of many low-molecular-weight metabolites and baseline contributions of macromolecules. Two approaches to model such spectra are compared and the data acquisition sequence, optimized for reproducibility, is presented. Modeling relies on prior knowledge constraints and linear combination of metabolite spectra. Investigated was what can be gained by basis parameterization, i.e., description of basis spectra as sums of parametric lineshapes. Effects of basis composition and addition of experimentally measured macromolecular baselines were investigated also. Both fitting methods yielded quantitatively similar values, model deviations, error estimates, and reproducibility in the evaluation of 64 spectra of human gray and white matter from 40 subjects. Major advantages of parameterized basis functions are the possibilities to evaluate fitting parameters separately, to treat subgroup spectra as independent moieties, and to incorporate deviations from straightforward metabolite models. It was found that most of the 22 basis metabolites used may provide meaningful data when comparing patient cohorts. In individual spectra, sums of closely related metabolites are often more meaningful. Inclusion of a macromolecular basis component leads to relatively small, but significantly different tissue content for most metabolites. It provides a means to quantitate baseline contributions that may contain crucial clinical information.

Fast assembly of fock matrices utilising symmetry properties of the basis set

A method of assembling the elements of the Fock matrix is described which is a modification of that due to Dacre. Lists of symmetry equivalent one-electron integrals are used as pointers to abbreviate the process of collecting two-electron integrals into the Fock matrix.

GCHF basis sets and their application in the electronic structure study of PrMnO3

The approach called generator coordinate Hartree-Fock (GCHF) method is used in the selection of Gaussian basis set [25s18p for O ((3)p), 31s21p14d for Mn (S-6), and 33s22p16d9f for Pr ((4)J)] for atoms. The role of the weight functions in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are contracted to (25s18p/9s5p), (31s21p14d/9s6p4d), and (33s22pl6d9f118sl2p5d3f) by segmented contraction scheme of Dunning and they are utilized in calculations of Restricted-Open-HF (ROHF) Total and Orbital energies of the (MnO+1)-Mn-3 and (PrO+1)-Pr-1 fragments, to evaluate their quality in molecular studies. The addition of one d polarization function in the contracted (9s5p) basis set for O(P-3) atom and their application with the contracted (9s6p4d), (18s21p5d3f) basis sets for Mn (S-6) and Pr-Pr ((4)j) atoms lead to the electronic structure study of PrMnO3. The dipole moment, the total energy, and total atomic charges properties were calculated and were carried out at ROHF level with the [PrMnO3](2) fragment. The calculated values show that PrMnO3 does not present piezoelectric properties. (C) 2004 Elsevier B.V. All rights reserved.

Gaussian basis sets to the theoretical study of the electronic structure of perovskite (LaMnO3)

The Generator Coordinate Hartree-Fock (GCHF) method is applied to generate extended (20s14p), (30s19p13d), and (31s23p18d) Gaussian basis sets for the 0, Mn, and La atoms, respectively. The role of the weight functions (WFs) in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are then contracted to [5s3p] and [11s6p6d] for 0 and Mn atoms, respectively, and [17s11p7d] for La atom by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaan (HFR) method for (MnO1+)-Mn-5 and (LaO1+)-La-1 fragments. The results obtained with the contracted basis sets are compared with values obtained with the extended basis sets. The addition of one d polarization function in the contracted basis set for 0 atom and its utilization with the contracted basis sets for Mn and La atoms leads to the calculations of dipole moment and total atomic charges of perovskite (LaMnO3). The calculations were performed at the HFR level with the crystal [LaMnO3](2) fragment in space group C-2v the values of dipole moment, total energy, and total atomic charges showed that it is reasonable to believe that LaMnO3 presents behaviour of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

The generator coordinate Hartree-Fock method as strategy for building Gaussian basis sets to ab initio study of the process of adsorption of sulfur on platinum (200) surface

The scheme named generator coordinate Hartree-Fock method (GCHF) is used to build (22s14p) and (33s22p16d9f) gaussian basis sets to S ((3)P) and Pt ((3)D) atoms, respectively. Theses basis sets are contracted to [13s10p] and [19s13p9d5f] through of Dunning's segmented contraction scheme and are enriched with d and g polarization functions, [13s10p1d] and [19s13p9d5flg]. Finally, the [19s13p9d5f1g] basis Set to Pt ((3)D) was supplemented with s and d diffuse functions, [20s13p10d5flg], and used in combination with [13s10p1d] to study the effects of adsorption of S ((3)D) atom on a pt ((3)D) atom belonged to infinite Pt (200) surface. Atom-atom overlap population, bond order, and infrared spectrum of [pt(_)S](2 -) were calculated properties and were carried out at Hartree-Fock-Roothaan level. The results indicate that the process of adsorption of S ((3)P) on pt ((3)D) in the infinite Pt (200) surface is mainly caused by a strong contribution of sigma between the 3p(z) orbital of S ((3)P) and the 6s orbital of pt ((3)D). (c) 2004 Elsevier B.V. All rights reserved.

Design of Gaussian basis sets to the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, tetraoxochromium (IV) ion, and tetraoxochromium (VI) ion

The Generator Coordinate Hartree-Fock (GCHF) method is employed to design 16s, 16s10p, 24s17p13d, 25s17p13d, and 26s17p Gaussian basis sets for the H ((2)S), O ((3)P), O(2-) ((1)S), Cr(3+) ((4)F), Cr(4+) ((3)F), and Cr(6+) ((1)S) atomic species. These basis sets are then contracted to (4s) for H ((2)S), (6s4p) for O ((3)P), and O(2-) ((1)S), (986p3d) for Cr(3+) ((4)F), (10s8p3d) for Cr(4+) ((3)F), and (13s7p) for Cr(6+) (1S) by a standard procedure. For evaluation of the quality of those basis sets in molecular calculations, we have accomplished studies of total and orbital (HOMO and HOMO-1) energies at the HF-Roothaan level for the molecular species of our interest. The results obtained with the contracted basis sets are compared to the values obtained with our extended basis sets and to the standard 6-311G basis set from literature. Finally, the contracted basis sets are enriched with polarization function and then utilized in the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, [Cr(H(2)O)(6)](3+), tetraoxochromium (IV) ion, [CrO(4)](4-), and tetraoxochromium (VI) ion, [CrO(4)](2-). The respective theoretical harmonic frequencies and IR-intensities were computed at the density functional theory (DFT) level. In the DFT calculations we employed the Becke's 1988 functional using the LYP correlation functional. The comparison between the results obtained and the corresponding experimental values indicates a very good description of the IR-spectra of the molecular ions studied, and that the GCHF method is still a legitimate alternative for selection of Gaussian basis sets. (C) 2003 Elsevier B.V. All rights reserved.

Designing Gaussian basis sets to the theoretical study of the piezoelectric effect of perovskite (BaTiO3)

Gaussian basis sets (24s14p, 30s19p14d, and 33s21p14d for O (P-3), Ti (S-5), and Ba (S-1) atoms, respectively), are designed with the strategy of the Generator Coordinate Hartree-Fock method. The basis sets are then contracted to [6s4p], [10s5p4d], and [16s9p5d] to O, Ti, and Ba atoms, respectively, and used in calculations of total and orbital energies of (TiO+2)-Ti-1 and (BaO)-Ba-1 fragments for quality evaluation in molecular studies. For O atom, the [6s4p] basis set is enriched with d polarization function and used along with the [10s5p4d] and [16s9p5d] basis sets for the theoretical study of the piezoelectric effect of perovskite (BaTiO3). The results of this work evidence that the piezoelectric properties in BaTiO3 can be caused by electrostatic interactions. (C) 2003 Elsevier B.V. All rights reserved.

Development of basis sets to calculations of the electronic structure of YMnO3

The generator coordinate Hartree-Fock method was used to develop 20s17p, 30s20p14d, and 30s21p16d Gaussian basis sets for the O ((3)p), Mn (S-6), and Y (D-2) atoms, respectively. The Gaussian basis sets were contracted to 20s17p/9s7p, 30s20p14d/11s7p7d, and 30s21p16d/14s7p7d and utilized in calculations of total energy and orbital energies of the (MnO1+)-Mn-5 and (YO1+)-Y-3 fragments to evaluate its quality in molecular studies. Finally, the contracted basis set for O atom was supplemented with one polarization function of d symmetry and used along with the other contracted basis sets (for Mn and Y) to calculate dipole moments, total energy, and total atomic charges in YMnO3 in space group D-6h. The analysis of those properties showed that is reasonable to believe that YMnO3 present behavior of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

Contracted GTF basis sets applied to the theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion

Contracted GTF basis sets designed with aid of the Generator Coordinate Hartree-Fock (GCHF) method for H(2S), O2-(1S), and Cr3+(4F) atomic species are applied to perform theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion. The 16s, 16s 10p, and 24s17p13d GTF basis sets were contracted to [4s] for H atom, [6s4p], and [9s6p3d] for O2- and Cr3+, respectively, by Dunning's scheme. For Cr3+, the [9s6p3d] basis set was enriched with f polarization function and used in combination com [4s] and [6s4p] in the study of our interest. The results obtained in this report show that the contracted GTF basis sets used are a useful alternative for the theoretical interpretation of Raman spectrum of hexaaquachromium(III) ion and that GCHF method is an effective alternative to selection of GTF basis sets for theoretical study of vibrational properties of poliatomic species. © 2003 Elsevier Science B.V. All rights reserved.

The ability of the Gaussian-2, Gaussian-3, Complete Basis Set–QB3, and Complete Basis Set–APNO model chemistries to model the geometries of small water clusters

The Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO methods have been used to calculate geometries of neutral clusters of water, (H2O)n, where n = 2–6. The structures are in excellent agreement with those determined from experiment and those predicted from previous high-level calculations. These methods also provide excellent thermochemical predictions for water clusters, and thus can be used with confidence in evaluating the structures and thermochemistry of water clusters.

Using phase-space localized basis functions to obtain vibrational energies of molecules

For decades scientists have attempted to use ideas of classical mechanics to choose basis functions for calculating spectra. The hope is that a classically-motivated basis set will be small because it covers only the dynamically important part of phase space. One popular idea is to use phase space localized (PSL) basis functions. This thesis improves on previous efforts to use PSL functions and examines the usefulness of these improvements. Because the overlap matrix, in the matrix eigenvalue problem obtained by using PSL functions with the variational method, is not an identity, it is costly to use iterative methods to solve the matrix eigenvalue problem. We show that it is possible to circumvent the orthogonality (overlap) problem and use iterative eigensolvers. We also present an altered method of calculating the matrix elements that improves the performance of the PSL basis functions, and also a new method which more efficiently chooses which PSL functions to include. These improvements are applied to a variety of single well molecules. We conclude that for single minimum molecules, the PSL functions are inferior to other basis functions. However, the ideas developed here can be applied to other types of basis functions, and PSL functions may be useful for multi-well systems.

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.