Reduced-size polarized basis sets for calculations of molecular electric properties. III. Second-row atoms

Reduced-size polarized (ZmPolX) basis sets are developed for the second-row atoms X = Si, P, S, and Cl. The generation of these basis sets follows from a simple physical model of the polarization effect of the external electric field which leads to highly compact polarization functions to be added to the chosen initial basis set. The performance of the ZmPolX sets has been investigated in calculations of molecular dipole moments and polarizabilities. Only a small deterioration of the quality of the calculated molecular electric properties has been found. Simultaneously the size of the present reduced-size ZmPolX basis sets is about one-third smaller than that of the usual polarized (PolX) sets. This reduction considerably widens the range of applications of the ZmPolX sets in calculations of molecular dipole moments, dipole polarizabilities, and related properties.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

Energy levels, wavefunction compositions and electric dipole transitions in neutral Ca

A configuration-interaction approach, based on the use of B-spline basis sets combined with a model potential including monoelectronic and dielectronic core polarization effects, is employed to calculate term energies and wavefunctions for neutral Ca. Results are reported for singlet and triplet bound states, and some quasi-bound states above the lowest ionization limit, with angular momentum up to L = 4. Comparison with experiment and with other theoretical results shows that this method yields the most accurate energy values for neutral Ca obtained to date. Wavefunction compositions, necessary for labelling the levels, and the effects of semi-empirical polarization potentials on the wavefunctions are discussed, as are some recent identifications of doubly-excited states. It is shown that taking into account dielectronic core polarization changes the energies of the lowest terms in Ca significantly, in general by a few hundred cm(-1), the effect decreasing rapidly for the higher bound states. For Rydberg states with n approximate to 7 the accuracy of the results is often better than a few cm(-1). For series members (or perturbers) with a pronounced 3d character the error can reach 150 cm(-1). The wavefunctions are used to calculate oscillator strengths and lifetimes for a number of terms and these are compared with existing measurements. The agreement is good but points to a need for improved measurements.

Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio $^{2}\Sigma^+$ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction (MRCI+Q)is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis sets (aug-pCV$n$Z-PP, n= 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X$^2\Sigma^+$ state (extrapolated to the CBS limit) is 16895.12 cm$^{-1}$ (2.094 eV), which agrees within 0.1$\%$ of a revised experimental value of <16910.6 cm$^{-1}$, while the calculated re is within 0.03 pm of the experimental result.

Many-body theory of positron-atom interactions

A many-body theory approach is developed for the problem of positron-atom scattering and annihilation. Strong electron- positron correlations are included nonperturbatively through the calculation of the electron-positron vertex function. It corresponds to the sum of an infinite series of ladder diagrams, and describes the physical effect of virtual positronium formation. The vertex function is used to calculate the positron-atom correlation potential and nonlocal corrections to the electron-positron annihilation vertex. Numerically, we make use of B-spline basis sets, which ensures rapid convergence of the sums over intermediate states. We have also devised an extrapolation procedure that allows one to achieve convergence with respect to the number of intermediate- state orbital angular momenta included in the calculations. As a test, the present formalism is applied to positron scattering and annihilation on hydrogen, where it is exact. Our results agree with those of accurate variational calculations. We also examine in detail the properties of the large correlation corrections to the annihilation vertex.

On velocity-based local model networks for nonlinear identification

This paper exposes the strengths and weaknesses of the recently proposed velocity-based local model (LM) network. The global dynamics of the velocity-based blended representation are directly related to the dynamics of the underlying local models, an important property in the design of local controller networks. Furthermore, the sub-models are continuous-time and linear providing continuity with established linear theory and methods. This is not true for the conventional LM framework, where the global dynamics are only weakly related to the affine sub-models. In this paper, a velocity-based multiple model network is identified for a highly nonlinear dynamical system. The results show excellent dynamical modelling performances, highlighting the value of the velocity-based approach for the design and analysis of LM based control. Three important practical issues are also addressed. These relate to the blending of the velocity-based local models, the use of normalised Gaussian basis functions and the requirement of an input derivative.

B-spline methods in R-matrix theory for scattering in two-electron systems

The use of B-spline basis sets in R-matrix theory for scattering processes has been investigated. In the present approach a B-spline basis is used for the description of the inner region, which is matched to the physical outgoing wavefunctions by the R-matrix. Using B-splines, continuum basis functions can be determined easily, while pseudostates can be included naturally. The accuracy for low-energy scattering processes is demonstrated by calculating inelastic scattering cross sections for e colliding on H. Very good agreement with other calculations has been obtained. Further extensions of the codes to quasi two-electron systems and general atoms are discussed as well as the application to (multi) photoionization.

INFLUENCE OF CORE POLARIZATION ON THE ELECTRON-AFFINITY OF CA

We have developed a method, based on the use of B-spline basis sets and model potentials, for determining properties of systems with two or three electrons outside a polarizable closed-shell core. It is applied to the calculation of the electron affinity of Ca and the resulting value of 17.7 meV is in excellent agreement with the most recent experiments. It is found that the dielectronic core-valence interaction reduces the electron affinity by 39.5 meV.

STO and GTO field-induced polarization functions for H to Kr

Field-induced polarization (FIP) functions were proposed over two decades ago to improve the accuracy of calculated response properties, and the FIP functions in GTO form for H and C to F were tested on small molecules, with encouraging results. The concept of FIP,is now extended to all atoms up to Kr. New simplifying approximations for the description of asymptotic highest occupied atomic orbitals. (HOAOs) are introduced in this study. They provide the basis for STO and GTO exponents of a complete set of FIP functions from H to Kr, which are both listed for the convenience of the users. Tests on the polarizabilities of a series of atoms and molecules demonstrate that addition of the FIP basis functions to a series' of standard basis sets drastically improves the performance of all these basis sets compared to converged results. Moreover, the byproduct of this study (approximate asymptotic HOAOs) provides information for the construction of accurate basis sets for long-range ground state properties. (C) 2003 Wiley Periodicals, Inc.

A study of velocity fields in the transition region of ε Eri (K2 V)

Analyses of the widths and shifts of optically thin emission lines in the ultraviolet spectrum of the active dwarf e Eri (K2 V) are presented. The spectra were obtained using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer. The linewidths are used to find the non-thermal energy density and its variation with temperature from the chromosphere to the upper transition region. The energy fluxes that could be carried by Alfvén and acoustic waves are investigated, to test their possible roles in coronal heating. Acoustic waves do not appear to be a viable means of coronal heating. There is, in principle, ample flux in Alfvén waves, but detailed calculations of wave propagation are required before definite conclusions can be drawn concerning their viability. The high sensitivity and spectral resolution of the above instruments have allowed two-component Gaussian fits to be made to the profiles of the stronger transition region lines. The broad and narrow components that result share some similarities with those observed in the Sun, but in e Eri the broad component is redshifted relative to the narrow component and contributes more to the total line flux. The possible origins of the two components and the energy fluxes implied are discussed. On balance our results support the conclusion of Wood, Linsky & Ayres, that the narrow component is related to Alfvén waves reaching to the corona, but the origin of the broad component is not clear.

Positron scattering and annihilation in hydrogenlike ions

Diagrammatic many-body theory is used to calculate the scattering phase shifts, normalized annihilation rates Zeff, and annihilation ? spectra for positron collisions with the hydrogenlike ions He+, Li2+, B4+, and F8+. Short-range electron-positron correlations and longer-range positron-ion correlations are accounted for by evaluating nonlocal corrections to the annihilation vertex and the exact positron self-energy. The numerical calculation of the many-body theory diagrams is performed using B-spline basis sets. To elucidate the role of the positron-ion repulsion, the annihilation rate is also estimated analytically in the Coulomb-Born approximation. It is found that the energy dependence and magnitude of Zeff are governed by the Gamow factor that characterizes the suppression of the positron wave function near the ion. For all of the H-like ions, the correlation enhancement of the annihilation rate is found to be predominantly due to corrections to the annihilation vertex, while the corrections to the positron wave function play only a minor role. Results of the calculations for s-, p-, and d-wave incident positrons of energies up to the positronium-formation threshold are presented. Where comparison is possible, our values are in excellent agreement with the results obtained using other, e.g., variational, methods. The annihilation-vertex enhancement factors obtained in the present calculations are found to scale approximately as 1+(1.6+0.46l)/Zi, where Zi is the net charge of the ion and l is the positron orbital angular momentum. Our results for positron annihilation in H-like ions provide insights into the problem of positron annihilation with core electrons in atoms and condensed matter systems, which have similar binding energies.

Molecular rearrangement reactions in the gas phase triggered by electron attachment

Bond formation and rearrangement reactions in gas phase electron attachment were studied through dissociative electron attachment (DEA) to pentafluorotoluene (PFT), pentafluoroaniline (PFA) and pentafluorophenol (PFP) in the energy range 0-14 eV. In the case of PFA and PFP, the dominant processes involve formation of [M - HF](-) through the loss of neutral HF. This fragmentation channel is most efficient at low incident electron energy and for PFP it is accompanied by a substantial conformational change of the anionic fragment. At higher energy, HF loss is also observed as well as a number of other fragmentation processes. Thermochemical threshold energies have been computed for all the observed fragments and classical trajectories of the electron attachment process were calculated to elucidate the fragmentation mechanisms. For the dominant reaction channel leading to the loss of HF from PFP, the minimum energy path was calculated using the nudged elastic band method.

Modeling study of woody biomass:Interactions of cellulose, hemicellulose, and lignin

Lignocellulosic biomass pretreatment and the subsequent thermal conversion processes to produce solid, liquid, and gas biofuels are attractive solutions for today's energy challenges. The structural study of the main components in biomass and their macromolecular complexes is an active and ongoing research topic worldwide. The interactions among the three main components, cellulose, hemicellulose, and lignin, are studied in this paper using electronic structure methods, and the study includes examining the hydrogen bond network of cellulose-hemicellulose systems and the covalent bond linkages of hemicellulose-lignin systems. Several methods (semiempirical, Hartree-Fock, and density functional theory) using different basis sets were evaluated. It was shown that theoretical calculations can be used to simulate small model structures representing wood components. By comparing calculation results with experimental data, it was concluded that B3LYP/6-31G is the most suitable basis set to describe the hydrogen bond system and B3LYP/6-31G(d,p) is the most suitable basis set to describe the covalent system of woody biomass. The choice of unit model has a much larger effect on hydrogen bonding within cellulose-hemicellulose system, whereas the model choice has a minimal effect on the covalent linkage in the hemicellulose-lignin system. © 2011 American Chemical Society.