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.

Crystal structure and ab initio calculations of a cyano-carbamimidic acid ethyl ester

The structure of one tautomer (amine form) of cyano-carbamimidic acid ethyl ester or (amino-ethoxy-methylidene)aminoformonitrile (CAS: 13947-84-7) was determined by single crystal X-ray diffraction. Ab initio quantum chemical calculations at the B3LYP, MP2 and G3 levels were performed to investigate the stability and the formation of the different tautomers and conformers. The calculations indicate that the amine form is the more stable tautomer, showing a high degree of election conjugation. The most stable amine conformer located by the calculations corresponds to the crystallized structure. On the contrary, in the less stable imine form, the conjugation is separated by a N2-C2 single bond. (C) 2007 Elsevier B.V. All rights reserved.

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Quasiparticle calculations of the electronic properties of ZrO2 and HfO2 polymorphs and their interface with Si

Quasiparticle calculations are performed to investigate the electronic band structures of various polymorphs of Hf and Zr oxides. The corrections with respect to density-functional-theory results are found to depend only weakly on the crystal structure. Based on these bulk calculations as well as those for bulk Si, the effect of quasiparticle corrections is also investigated for the band offsets at the interface between these oxides and Si assuming that the lineup of the potential at the interface is reproduced correctly within density-functional theory. On the one hand, the valence-band offsets are practically unchanged with a correction of a few tenths of electron volts. On the other hand, conduction-band offsets are raised by 1.3-1.5 eV. When applied to existing calculations for the offsets at the density-functional-theory level, our quasiparticle corrections provide results in good agreement with the experiment.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Monte Carlo calculations of quantum yield in inhomogeneous PtSi/p-Si Schottky barriers

Monte Carlo calculations of quantum yield in PtSi/p-Si infrared detectors are carried out taking into account the presence of a spatially distributed barrier potential. In the 1-4 mu m wavelength range it is found that the spatial inhomogeneity of the barrier has no significant effect on the overall device photoresponse. However, above lambda = 4.0 mu m and particularly as the cut-off wavelength (lambda approximate to 5.5 mu m) is approached, these calculations reveal a difference between the homogeneous and inhomogeneous barrier photoresponse which becomes increasingly significant and exceeds 50% at lambda = 5.3 mu m. It is, in fact, the inhomogeneous barrier which displays an increased photoyield, a feature that is confirmed by approximate analytical calculations assuming a symmetric Gaussian spatial distribution of the barrier. Furthermore, the importance of the silicide layer thickness in optimizing device efficiency is underlined as a trade-off between maximizing light absorption in the silicide layer and optimizing the internal yield. The results presented here address important features which determine the photoyield of PtSi/Si Schottky diodes at energies below the Si absorption edge and just above the Schottky barrier height in particular.

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.

Relativistic radiatively damped R -matrix calculations of the electron-impact excitation of helium-like iron and krypton

Electron-impact excitation data for He-like ions are of significant importance for diagnostic applications to both laboratory and astrophysical plasmas. Here we report on the first fully relativistic R -matrix calculations with radiation damping for the He-like ions Fe 24+ and Kr 34+ . Effective collision strengths for these two ions have been determined with and without damping over a wide temperature range for all transitions between the 49 levels through n = 5. We find that damping has a pronounced effect on the effective collision strengths for excitation to some of the low-lying levels, but its effect on excitation to the vast majority of levels is small. At the energy of a resonance peak, we also investigate the effect of radiation damping on the angular distribution of scattered electrons. Finally, we compare our results for Fe 24+ with an earlier intermediate coupling frame transformation R -matrix calculation with radiation damping by Whiteford et al ( J. Phys. B: At. Mol. Opt. Phys. 34 3179) and find good agreement, especially for excitation to the lower levels.

Intermediate-coupling R -matrix calculations of electron-impact excitation of Fe 5+

For applications to laboratory and astrophysical plasmas, there is a great need for accurate electron-impact excitation data between individual levels in the lower charge-state ions of iron. Recently, we have reported on the first intermediate-coupling R -matrix calculation of electron-impact excitation in Fe 4+ , in which the close-coupling expansion of the target included levels from both ground and excited configurations (Ballance et al 2007 J. Phys. B: At. Mol. Opt. Phys. [/0953-4075/40/23/f01] 40 F327 , 2008 Europhys. News 39 14). In this paper, we present the results of two large intermediate-coupling Dirac R -matrix calculations of electron-impact excitation of Fe 5+ . The results from the two calculations, which differ only in the configuration–interaction expansions of the target, are compared. These comparisons provide some indication of the accuracy of the calculations and the resulting data should be useful in modelling plasmas containing iron.

Dirac R -matrix calculations of electron-impact excitation of neon-like krypton

We have employed the Dirac R -matrix method to determine electron-impact excitation cross sections and effective collision strengths in Ne-like Kr 26+ . Both the configuration-interaction expansion of the target and the close-coupling expansion employed in the scattering calculation included 139 levels up through n = 5. Many of the cross sections are found to exhibit very strong resonances, yet the effects of radiation damping on the resonance contributions are relatively small. Using these collisional data along with multi-configuration Dirac–Fock radiative rates, we have performed collisional-radiative modeling calculations to determine line-intensity ratios for various radiative transitions that have been employed for diagnostics of other Ne-like ions.

Collisional-radiative calculations of He line emission in low-temperature plasmas

We present spectral modeling results for neutral helium. Our underlying atomic data contains radiative transition rates that are generated from atomic structure calculations and electron-impact excitation rates, that are determined from both the standard R-matrix method and the R-matrix with pseudostates RMPS method. In this paper, we focus on transitions of particular importance to diagnostic line ratios. For example, our calculated rate coefficient for the electron-impact transition 1s3s 1S→1s3p 1P, which has a pronounced effect on the 728.1 nm diagnostic spectral line, is found to be in good agreement with previous experimental mea- surements. We also consider transitions from the 1s2 1S ground and 1s2s 3S terms to terms of the n=4 shell. They are found to be affected significantly by coupling of the bound states to the target continuum continuum coupling, which is included in our RMPS calculation, but not in our standard R-matrix calculation. We perform collisional-radiative calculations to determine spectral line intensity ratios for three ratios of particular interest, namely the 504.8 nm/471.3 nm, 492.2 nm/471.3 nm, and 492.2 nm/504.8 nm line ratios. Comparing our results determined from the RMPS excitation rates with those from the standard R-matrix excitation rates, we find that continuum coupling affects the rate coefficients significantly, leading to different values for all three line ratios. We also compare our modeling results with spectral measurements taken recently on the Auburn Helicon plasma device, finding that the ground and metastable populations are not in equilibrium, and that the experimental measurements are more consistent with the 1s2s 3S metastable term populations being short lived in the plasma.

Molecular mechanics calculations of tin-bis (triphenylphosphine oxide) complexes

The capability of molecular mechanics for modeling the wide distribution of bond angles and bond lengths characteristic of coordination complexes was investigatecl. This was the preliminary step for future modeling of solvent extraction. Several tin-phosphine oxide COrnI)le:){es were selected as the test groUl) for t.he d,esired range of geometry they eX!libi ted as \-vell as the ligands they cOD.tained r Wllich were c\f interest in connection with solvation. A variety of adjustments were made to Allinger's M:M2 force·-field ill order to inl.prove its performance in the treatment of these systems. A set of u,nique force constants was introduced for' those terms representing the metal ligand bond lengths, bond angles, and, torsion angles. These were significantly smaller than trad.itionallY used. with organic compounds. The ~1orse poteIlt.ial energ'Y function was incorporated for the M-X l')ond lE~ngths and the cosine harmonic potential erlerg-y function was invoked for the MOP bond angle. These functions were found to accomodate the wide distribution of observed values better than the traditional harmonic approximations~ Crystal packing influences on the MOP angle were explored thr"ollgh ttle inclusion of the isolated molecule withil1 a shell cc)ntaini11g tl1e nearest neigl1'bors duri.rlg energy rninimization experiments~ This was found to further improve the fit of the MOP angle.

Systematic trends in self-consistent calculations of linear quantum wires

Systematic trends in the properties of a linear split-gate heterojunction are studied by solving iteratively the Poisson and Schrödinger equations for different gate potentials and temperatures. A two-dimensional approximation is presented that is much simpler in the numerical implementation and that accurately reproduces all significant trends. In deriving this approximation, we provide a rigorous and quantitative basis for the formulation of models that assumes a two-dimensional character for the electron gas at the junction.