A basis set comparison in a variational scheme for the Yukawa potential

The variational method is used to obtain solutions to Schrödinger's equation for a particle in the radially screened Yukawa potential. A basis set is presented. While the Laguerre basis set shows considerable improvement over the hydrogenic one, problems are still encountered as the screening parameter approaches its threshold value. Variational calculations are also presented using an Eckart-type basis set which looks promising near the critical screening region.     

An auto-coherent variational scheme for Yukawa potential

The variational principle is used to obtain solutions to Schrödinger's equation for a particle in the Yukawa potential. A Laguerre basis set extended by an extra function is employed in the calculations. A special parameter used in the extra function and its relation with the systems energy results in utilizing an auto-coherent scheme. Considerable improvement seems to be achieved especially in the critical region where the screening parameter approaches its threshold value.     

Convergent perturbation studies in screened coulomb potential systems: a high precision numerical algorithm via Laguerre basis set

In this work, a high precision algorithm is developed to determine the discrete spectrum of screened coulomb potential systems. The algorithm is mainly based upon the use of the perturbation of a hydrogen-like operator by a bounded operator. The Laguerre basis set expansion is employed in the procedure to obtain the operator inversion. Although a functional analytic analysis of errors and proof of convergence theorem are still lacking, it appears, numerically, that the method rapidly converges for bounded screened coulomb potential. Extremely accurate numerical results for the bound-state energies, in the case of Yukawa potential, are presented for illustrative purposes.     

Analytical treatment of the oscillating Yukawa potential

Using a suitable Laguerre basis set that ensures a tridiagonal matrix representation of the reference Hamiltonian, we were able to evaluate in closed form the matrix elements of the generalized Yukawa potential with a complex screening parameter. This enabled us to treat analytically both the cosine and sine-like Yukawa potentials on equal footing and compute their bound states spectrum as the eigenvalues of the associated analytical matrix representing their Hamiltonians. Finally we used a carefully designed complex scaling method to evaluate the resonance energies and compared our results satisfactorily with those obtained in the literature for the cosine-like Yukawa potential.     

超球坐标下原子、分子体系的变分计算：—氦原子和氢负离子基态

应用超球会标表示氦原子和氢负离子的薛定谔方程，将二电子原子在三维空间中的运动转化为单电子原子在六维空间中受广义库仑力作用的运动，我们给出了六维空间广义角动量算符的本征值与本征函数，并以此本征函数微基构造超球波函数，得到超球径微分方程，以广义Laguerre 多项式表示超球径波函数，运用密度矩阵和线性变分法得到非正交基下超球径波函数满足的久期方程，最后求得能量和波函数，计算结果与精确的计算符合良好。     

The multiscale coarse-graining method. VIII. Multiresolution hierarchical basis functions and basis function selection in the construction of coarse-grained force fields

The multiscale coarse-graining (MS-CG) method is a method for determining the effective potential energy function for a coarse-grained (CG) model of a molecular system using data obtained from molecular dynamics simulation of the corresponding atomically detailed model. The coarse-grained potential obtained using the MS-CG method is a variational approximation for the exact many-body potential of mean force for the coarse-grained sites. Here we propose a new numerical algorithm with noise suppression capabilities and enhanced numerical stability for the solution of the MS-CG variational problem. The new method, which is a variant of the elastic net method [Friedman et al., Ann. Appl. Stat. 1, 302 (2007)], allows us to construct a large basis set, and for each value of a so-called "penalty parameter" the method automatically chooses a subset of the basis that is most important for representing the MS-CG potential. The size of the subset increases as the penalty parameter is decreased. The appropriate value to choose for the penalty parameter is the one that gives a basis set that is large enough to fit the data in the simulation data set without fitting the noise. This procedure provides regularization to mitigate potential numerical problems in the associated linear least squares calculation, and it provides a way to avoid fitting statistical error. We also develop new basis functions that are similar to multiresolution Haar functions and that have the differentiability properties that are appropriate for representing CG potentials. We demonstrate the feasibility of the combined use of the elastic net method and the multiresolution basis functions by performing a variational calculation of the CG potential for a relatively simple system. We develop a method to choose the appropriate value of the penalty parameter to give the optimal basis set. The combined effect of the new basis functions and the regularization provided by the elastic net method opens the possibility of using very large basis sets for complicated CG systems with many interaction potentials without encountering numerical problems in the variational calculation.     

On the mathematical nature of Guseinov’s rearranged one-range addition theorems for Slater-type functions

Starting from one-range addition theorems for Slater-type functions, which are expansion in terms of complete and orthonormal functions based on the generalized Laguerre polynomials, Guseinov constructed addition theorems that are expansions in terms of Slater-type functions with a common scaling parameter and integral principal quantum numbers. This was accomplished by expressing the complete and orthonormal Laguerre-type functions as finite linear combinations of Slater-type functions and by rearranging the order of the nested summations. Essentially, this corresponds to the transformation of a Laguerre expansion, which in general only converges in the mean, to a power series, which converges pointwise. Such a transformation is not necessarily legitimate, and this contribution discusses in detail the difference between truncated expansions and the infinite series that result in the absence of truncation.     

Representation of cusps in a hyperspherical basis set

The wave functions of Coulomb systems have cusps at points corresponding to two particle coelescences. In this paper, we derive series representing the cusps in terms of hyperspherical harmonics multiplied by functions of the hyperradius. When the hyperspherical method is applied to Coulomb systems, the harmonics which appear in these series should be included in the hyperangular basis set.     

An investigation into the use of expansion methods in the calculation of chemical shifts and diamagnetic susceptibilities

The chemical shift and the diamagnetic susceptibility of the hydrogenic atom with magnetic dipole and origin of the external magnetic field vector potential noncoincident with the hydrogenic nucleus have been calculated from perturbation theory using a set of expansion functions whose radial parts are single exponent associated Laguerre functions. In contrast to hydrogenic expansion functions these functions give rapid convergence to the exact values of the second-order energy summations when centred at the hydrogenic nucleus. The rate of convergence is fairly insensitive to the choice of expansion function exponent.     

The correlation consistent composite approach (ccCA): an alternative to the Gaussian-n methods

An alternative to the Gaussian-n (G1, G2, and G3) composite methods of computing molecular energies is proposed and is named the "correlation consistent composite approach" (ccCA, ccCA-CBS-1, ccCA-CBS-2). This approach uses the correlation consistent polarized valence (cc-pVXZ) basis sets. The G2-1 test set of 48 enthalpies of formation (DeltaHf), 38 adiabatic ionization potentials (IPs), 25 adiabatic electron affinities (EAs), and 8 adiabatic proton affinities (PAs) are computed using this approach, as well as the DeltaHf values of 30 more systems. Equilibrium molecular geometries and vibrational frequencies are obtained using B3LYP density functional theory. When applying the ccCA-CBS method with the cc-pVXZ series of basis sets augmented with diffuse functions, mean absolute deviations within the G2-1 test set compared to experiment are 1.33 kcal mol(-1) for DeltaHf,0.81 kcal mol(-1) for IPs, 1.02 kcal mol(-1) for EAs, and 1.51 kcal mol(-1) for PAs, without including the "high-level correction" (HLC) contained in the original Gn methods. Whereas the HLC originated in the Gaussian-1 method as an isogyric correction, it evolved into a fitted parameter that minimized the error of the composite methods, eliminating its physical meaning. Recomputing the G1 and G3 enthalpies of formation without the HLC reveals a systematic trend where most DeltaHf values are significantly higher than experimental values. By extrapolating electronic energies to the complete basis set (CBS) limit and adding G3-like corrections for the core-valence and infinite-order electron correlation effects, ccCA-CBS-2 often underestimates the experimental DeltaHf, especially for larger systems. This is desired as inclusion of relativistic and atomic spin-orbit effects subsequently improves theoretical DeltaHf values to give a 0.81 kcal mol(-1) mean absolute deviation with ccCA-CBS-2. The ccCA-CBS method is a viable "black box" method that can be used on systems with at least 10-15 heavy atoms.     

Bond critical points in the electronic structures of binary hydrides

The bond critical points of the binary hydrides formed by the elements of the first two rows of the periodic table have been calculated. Particular attention has been paid to the basis-set dependence of the bond critical points at the experimental equilibrium geometries, or where necessary at model geometries. With the exception of H₂S, stepwise extension of the basis set leads to a smooth convergence of the bond critical points to a set of values which appear to converge to the Hartree–Fock limit. For H₂S it is shownb that the position of the bodn critical point is not only more sensitive to the presence of polarization functions in the basis set, but depends strongly on the orbital exponents of the polarization functions. Extensive optimizations of the exponents of the polarization functions have been carried out with the (12s9p/5s) basis set for second-row hydrides. The effects of contracting the Huzinaga basis sets have been examined.     

Use of locally dense basis sets for nuclear magnetic resonance shielding calculations

The locally dense basis set approach to the calculation of nuclear magnetic resonance shieldings is one in which a sufficiently large or dense set of basis functions is used for an atom or molecular fragment containing the resonant nucleus or nuclei of interest and fewer or attenuated sets of basis functions employed elsewhere. Provided the dense set is of sufficient size, this approach is capable of determining chemical shieldings nearly as well as a calculation with a balanced basis set of quality equal to the locally dense set, but with considerable savings of CPU time. Detailed comparisons are provided of locally dense and balanced calculations in the gauge including atomic orbital (GIAO) method for the individual principal values, the isotropic shieldings, and the tensor orientations for hydrogen, carbon, nitrogen, oxygen, fluorine, and phosphorus nuclei. It is seen that chemical functional groups can often define the appropriate molecular fragment to be taken locally dense. While the present test cases are for the most part small molecules, the value of the method is that it will allow calculations on systems that would otherwise presently be computationally expensive or inaccessible. © John Wiley & Sons, Inc.     

Prediction of emission wavelengths of phosphorescent NHC based emitters for OLEDs

Some years ago we developed a general method for the prediction of the emission wavelengths for cyclometalated NHC complexes. We intentionally used a standard DFT functional with a moderate double-zeta basis set (B3LYP/6-31G*) to be able to use this method for fast screening. Since then we have synthetized a wide variety of different ligand motifs, which are far beyond the initial scope. We herein report the results obtained using our method on recently synthesized complexes, including, e.g. bimetallic systems. The predicted wavelengths of the emission are in an excellent agreement with the measured experimental values. We report the results for thirty-four complexes with different electronic and steric properties. For the majority of the systems the differences between experimentally determined and predicted values are within 10 nm.     

Vibrational energies of PH3 calculated variationally at the complete basis set limit

The potential energy surface for the electronic ground state of PH(3) was calculated at the CCSD(T) level using aug-cc-pV(Q+d)Z and aug-cc-pVQZ basis sets for P and H, respectively, with scalar relativistic corrections included. A parametrized function was fitted through these ab initio points, and one parameter of this function was empirically adjusted. This analytical PES was employed in variational calculations of vibrational energies with the newly developed program TROVE. The convergence of the calculated vibrational energies with increasing vibrational basis set size was improved by means of an extrapolation scheme analogous to the complete basis set limit schemes used in ab initio electronic structure calculations. The resulting theoretical energy values are in excellent agreement with the available experimentally derived values.     

Calculations of D‐wave bound states and resonance states of the screened helium atom using correlated exponential wave functions

We have investigated the effects of screened Coulomb (Yukawa) potentials on the bound ^1,3D states and the doubly excited ^1,3 D^e resonance states of helium atom using highly correlated exponential basis functions. The Density of resonance states are calculated using stabilization method. Highly correlated exponential basis functions are used to consider the correlation effect between the charged particles. A total of 18 resonances (nine each for ¹ D^e and ³ D^e states) below the n = 2 He ⁺ threshold has been calculated. For each spin states, this includes four members in the 2pnp series, three members in the 2snd series, and two members in 2pnf series. The resonance energies and widths for various screening parameters ranging from infinity to a small value for these ^1,3 D^e resonance states are reported along with the bound‐excited 1s3d ^1,3 D state energies. Overall behavior of the spectral profile of 1s3d ¹D state of helium atom due to electron‐electron and electron‐nucleus screening are also presented. Accurate resonance energies and widths are also reported for He in vacuum. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

On the counterpoise correction for the basis set superposition error in large systems

The appropriateness of the use of the counterpoise correction for the basis set superposition error in SCF calculations of the interaction energies for pairs of aliphatic amino acids is analyzed in this paper. Our results show that for this type of molecule where the magnitude of the basis set superposition error can become quite big, the use of the counterpoise method provides interaction energies in good agreement with near Hartree-Fock values. The inaccuracies associated with the counterpoise method are much less important compared with the basis set superposition error itself. It is shown that the use of a well-balanced minimal basis set together with the counterpoise method is a good compromise (quality versus computational cost) for calculating interaction energies in systems involving molecules of biological interest.     

Optimal grids for generalized finite basis and discrete variable representations: definition and method of calculation

The method of optimal generalized finite basis and discrete variable representations (FBR and DVR) generalizes the standard, Gaussian quadrature grid-classical orthonormal polynomial basis-based FBR/DVR method to general sets of grid points and to general, nondirect product, and/or nonpolynomial bases. Here, it is shown how an optimal set of grid points can be obtained for an optimal generalized FBR/DVR calculation with a given truncated basis. Basis set optimized and potential optimized grids are defined. The optimized grids are shown to minimize a function of grid points derived by relating the optimal generalized FBR of a Hamiltonian operator to a non-Hermitian effective Hamiltonian matrix. Locating the global minimum of this function can be reduced to finding the zeros of a function in the case of one dimensional problems and to solving a system of D nonlinear equations repeatedly in the case of D>1 dimensional problems when there is an equal number of grid points and basis functions. Gaussian quadrature grids are shown to be basis optimized grids. It is demonstrated by a numerical example that an optimal generalized FBR/DVR calculation of the eigenvalues of a Hamiltonian operator with potential optimized grids can have orders of magnitude higher accuracy than a variational calculation employing the same truncated basis. Nevertheless, for numerical integration with the optimal generalized FBR quadrature rule basis optimized grids are the best among grids of the same number of points. The notions of Gaussian quadrature and Gaussian quadrature accuracy are extended to general, multivariable basis functions.     

A self‐guided search for good local minima of the sum‐of‐squared‐error in nonlinear least squares regression

Hard modeling of nonlinear chemical or biological systems is highly relevant as a model function together with values for model parameters provides insights in the systems' functionalities. Deriving values for said model parameters via nonlinear regression, however, is challenging as usually one of the numerous local minima of the sum‐of‐squared errors (SSEs) is determined; furthermore, for different starting points, different minima may be found. Thus, nonlinear regression is prone to low accuracy and low reproducibility. Therefore, there is a need for a generally applicable, automated initialization of nonlinear least squares algorithms, which reaches a good, reproducible solution after spending a reasonable computation time probing the SSE‐hypersurface. For this purpose, a three‐step methodology is presented in this study. First, the SSE‐hypersurface is randomly probed in order to estimate probability density functions of initial model parameter that generally lead to an accurate fit solution. Second, these probability density functions then guide a high‐resolution sampling of the SSE‐hypersurface. This second probing focuses on those model parameter ranges that are likely to produce a low SSE. As the probing continues, the most appropriate initial guess is retained and eventually utilized in a subsequent nonlinear regression. It is shown that this “guided random search” derives considerably better regression solutions than linearization of model functions, which has so far been considered the best‐case scenario. Examples from infrared spectroscopy, cell culture monitoring, reaction kinetics, and image analyses demonstrate the broad and successful applicability of this novel method. Copyright © 2014 John Wiley & Sons, Ltd.     

On the nature of the Møller-Plesset critical point

It has been suggested [F. H. Stillinger, J. Chem. Phys. 112, 9711 (2000)] that the convergence or divergence of M?ller-Plesset perturbation theory is determined by a critical point at a negative value of the perturbation parameter z at which an electron cluster dissociates from the nuclei. This conjecture is examined using configuration-interaction computations as a function of z and using a quadratic approximant analysis of the high-order perturbation series. Results are presented for the He, Ne, and Ar atoms and the hydrogen fluoride molecule. The original theoretical analysis used the true Hamiltonian without the approximation of a finite basis set. In practice, the singularity structure depends strongly on the choice of basis set. Standard basis sets cannot model dissociation to an electron cluster, but if the basis includes diffuse functions then it can model another critical point corresponding to complete dissociation of all the valence electrons. This point is farther from the origin of the z plane than is the critical point for the electron cluster, but it is still close enough to cause divergence of the perturbation series. For the hydrogen fluoride molecule a critical point is present even without diffuse functions. The basis functions centered on the H atom are far enough from the F atom to model the escape of electrons away from the fluorine end of the molecule. For the Ar atom a critical point for a one-electron ionization, which was not previously predicted, seems to be present at a positive value of the perturbation parameter. Implications of the existence of critical points for quantum-chemical applications are discussed.     

Cluster models and ab initio calculations of (19)F NMR isotropic chemical shifts for inorganic fluorides

(19)F NMR isotropic chemical shift (delta(iso)) calculations are performed in crystallized compounds using the GIAO method with the B3LYP hybrid functional at DFT level. Clusters centered on the studied fluorine atoms mimic the crystalline structures. The 6-311+G(d) basis set is chosen for the central fluorine atom, and the LanL2DZ basis set for the others. The metal atoms are described by the 3-21G(2d) basis set or, when not available, by the CRENBL basis set with the corresponding ECP, and augmented with 2d polarization functions when existing. First, for high-symmetry systems (MF, MF(2), and MF(3) compounds), a systematization of the cluster building up from coordination spheres is proposed, generalized to fluoroperovskites and fluoroaluminates KAlF(4) and RbAlF(4). When applied to rather low symmetry systems such as barium fluorometalates BaMgF(4), BaZnF(4), and Ba(2)ZnF(6), the definition of the coordination spheres is far from easy. Then, for structures built up from a MF(6) octahedron network, we may define different "starting clusters": [FM(2)F(8)] for the shared fluorine atoms, [FMF(4)] for the unshared ones, and [FBa(4)](7+) for the "free" ones. Analogous "starting clusters" are then tested on compounds from the NaF-AlF(3), BaF(2)-AlF(3), and CaF(2)-AlF(3) binary systems and for alpha-BaCaAlF(7) that are also built up from a MF(6) octahedron network. For each of these corresponding fluorine sites, delta(iso) values are calculated with the "starting clusters" and several larger clusters and compared to the experimental delta(iso) values. For the barium-containing clusters, the RMS deviation is equal to 51 ppm. It is suggested that this result may be related to the poor quality of the barium basis sets for which no polarization functions are available for the moment. In total, chemical shifts were calculated for 122 fluorine sites, in a various range of compounds. For the clusters without barium, the ab initio method leads to a RMS equal to 22 ppm, which is a quite nice result keeping in mind that the (19)F chemical shift range is larger than 200 ppm.