Low-lying states of heavy nuclei within the nucleon pair approximation

In this article we perform systematic calculations on low-lying states of 33 nuclei with A=202-212, using the nucleon pair approximation of the shell model. We use a phenomenological shell-model Hamiltonian that includes single-particle energies, monopole and quadrupole pairing interactions, and quadrupole-quadrupole interactions. The building blocks of our model space include one J=4 valence neutron pair, and one J=4,6,8 valence proton pair, in addition to the usual S and D pairs. We calculate binding energies, excitation energies, electric quadrupole and magnetic dipole moments of low-lying states, and E2 transition rates between low-lying states. Our calculated results are reasonably consistent with available experimental data. The calculated quadrupole moments and magnetic moments, many of which have not yet been measured for these nuclei, are useful for future experimental measurements.

Quantum-mechanical calculations of vibrationally resolved cross sections for non-dissociative charge transfer of O3+ with H-2

Charge transfer due to collisions of ground state O3+ (2s(2)2p P-2) ions with molecular hydrogen is investigated using the quantum-mechanical molecular-orbital close-coupling (MOCC) method, and electronic and vibrational state-selective cross sections along with the corresponding differential cross sections are calculated for projectile energies of 100, 500, 1000 and 5000 eV/u at the orientation angles of 25 degrees,45 degrees and 89 degrees. The adiabatic potentials and radial coupling matrix elements utilized in the QMOCC calculations were obtained with the spin-coupled valence-bond approach. The infinite order sudden approximation (IOSA) and the vibrational sudden approximation (VSA) are utilized to deal with the rotation of H-2 and the coupling between the electron and the vibration of H-2. It is found that the distribution of vibrationally resolved cross sections with the vibrational quantum number upsilon' of H-2(+) (upsilon') varies with the increment of the projectile energy; and the electronic and vibrational stateselective differential cross sections show similar behaviors: there is a highest platform within a very small scattering angle, beyond which the differential cross sections decrease as the scattering angle increases and lots of oscillating structures appear, where the scattering angle of the first structure decreases as E-P(-1/2) with the increment of the projectile energy E-P; and the structure and amplitude of the differential cross sections are sensitive to the orientation of molecule H-2, which provides a possibility to identify the orientations of molecule H-2 by the vibrational state-selective differential scattering processes.

Time-dependent density-functional calculations for optical spectra of Na-2 and Na-4 clusters

With the frame of the time-dependent local density approximation, an efficient description of the optical response of clusters has been used to study the photo-absorption cross section of Na-2 and Na-4 clusters. It is shown that our calculated results are in good agreement with the experiment. In addition, our calculated spectrum for the Na-4 cluster is in better agreement with experiment than the GW absorption spectrum.

Isospin-dependent pairing interaction from nuclear matter calculations

The isospin dependence of the effective pairing interaction is discussed on the basis of the Bardeen, Cooper, and Schrieffer theory of superfluid asymmetric nuclear matter. It is shown that the energy gap, calculated within the mean field approximation in the range from symmetric nuclear matter to pure neutron matter, is not linearly dependent on the symmetry parameter owing to the nonlinear structure of the gap equation. Moreover, the construction of a zero-range effective pairing interaction compatible with the neutron and proton gaps in homogeneous matter is investigated, along with some recent proposals of isospin dependence tested on the nuclear data table.

band Structure calculations on orthorhombic bulk and nanocrystalline SrY2O4

The electronic structure of SrY2O4 is calculated by using a density functional method, and the exchange and correlation have been treated by using a the generalized gradient approximation (GGA) within the scheme due to Perdew, Burke, and Ernzerhof (PBE). SrY2O4 is predicted to be a direct-gap material because the top of the valence band and the bottom of the conduction band are along the same direction at G. The bond length and the bond covalency are also calculated by using a chemical bond method.

Magnetic structure and orbital ordering in tetragonal and monoclinic KCrF3 from first-principles calculations

KCrF3 has been systematically investigated by using the full-potential linearized augmented plane wave plus local orbital method within the generalized gradient approximation and the local spin density approximation plus the on-site Coulomb repulsion approach. The total energies for ferromagnetic and three different antiferromagnetic configurations are calculated in the high-temperature tetragonal and low-temperature monoclinic phases, respectively.

Application of the Sanchez-Lacombe lattice fluid theory to the system pvme/ps and model calculations

Cloud point curves reported in the literature for five representatives of the system poly(vinyl methyl ether)/polystyrene were evaluated theoretically by means of the Sanchez-Lacombe lattice fluid theory. The measured phase separation behavior can be described within experimental error using only one adjustable parameter (quantifying the interaction between the unlike mers). The Flory-Huggins interaction parameters calculated from this theoretical description depend in good approximation linearly on composition (volume fractions) and on the inverse temperature. An evaluation of these data yields a maximum heat effect which is almost one order of magnitude less (ca. -0.25 J/cm(3)) than obtained via Hess's cycle (dissolution of the components and of the blend) from calorimetric measurements. Model calculations on the basis of the present theory demonstrate that the critical points shift to a different extent upon a certain relative change in the molar mass of the blend components. The sensitivity of the calculated phase diagrams against changes in the scaling parameter decreases in the following order: interaction energies between unlike mers, differences in the scaling temperatures, pressures and densities.

Investigation of numerical atomic orbitals for first-principles calculations of the electronic and transport properties of silicon nanowire structures

This thesis is focused on the application of numerical atomic basis sets in studies of the structural, electronic and transport properties of silicon nanowire structures from first-principles within the framework of Density Functional Theory. First we critically examine the applied methodology and then offer predictions regarding the transport properties and realisation of silicon nanowire devices. The performance of numerical atomic orbitals is benchmarked against calculations performed with plane waves basis sets. After establishing the convergence of total energy and electronic structure calculations with increasing basis size we have shown that their quality greatly improves with the optimisation of the contraction for a fixed basis size. The double zeta polarised basis offers a reasonable approximation to study structural and electronic properties and transferability exists between various nanowire structures. This is most important to reduce the computational cost. The impact of basis sets on transport properties in silicon nanowires with oxygen and dopant impurities have also been studied. It is found that whilst transmission features quantitatively converge with increasing contraction there is a weaker dependence on basis set for the mean free path; the double zeta polarised basis offers a good compromise whereas the single zeta basis set yields qualitatively reasonable results. Studying the transport properties of nanowire-based transistor setups with p+-n-p+ and p+-i-p+ doping profiles it is shown that charge self-consistency affects the I-V characteristics more significantly than the basis set choice. It is predicted that such ultrascaled (3 nm length) transistors would show degraded performance due to relatively high source-drain tunnelling currents. Finally, it is shown the hole mobility of Si nanowires nominally doped with boron decreases monotonically with decreasing width at fixed doping density and increasing dopant concentration. Significant mobility variations are identified which can explain experimental observations.

Parallel bandwidth characteristics calculations for thin avalanche photodiodes on a SGI Origin 2000 supercomputer

An important factor for high-speed optical communication is the availability of ultrafast and low-noise photodetectors. Among the semiconductor photodetectors that are commonly used in today’s long-haul and metro-area fiber-optic systems, avalanche photodiodes (APDs) are often preferred over p-i-n photodiodes due to their internal gain, which significantly improves the receiver sensitivity and alleviates the need for optical pre-amplification. Unfortunately, the random nature of the very process of carrier impact ionization, which generates the gain, is inherently noisy and results in fluctuations not only in the gain but also in the time response. Recently, a theory characterizing the autocorrelation function of APDs has been developed by us which incorporates the dead-space effect, an effect that is very significant in thin, high-performance APDs. The research extends the time-domain analysis of the dead-space multiplication model to compute the autocorrelation function of the APD impulse response. However, the computation requires a large amount of memory space and is very time consuming. In this research, we describe our experiences in parallelizing the code in MPI and OpenMP using CAPTools. Several array partitioning schemes and scheduling policies are implemented and tested. Our results show that the code is scalable up to 64 processors on a SGI Origin 2000 machine and has small average errors.

Multi-state CDW approximation to electron capture during slow p-H collisions

The first complete multi-state CDW close coupling calculations which use a fully normalized basis set are performed. The results obtained at impact energies in the region of 10 keV for total and n = 2 capture cross sections are in reasonably good accord with experiment despite the fact that only the ground states of both species and the n = 2 states of the projectile are incorporated into the model. The theory has significant advantages over other atomic and molecular expansions which may require extensive bases to obtain similar accuracy.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

yambo: An ab initio tool for excited state calculations

yambo is an ab initio code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the GW approximation for the self-energy. Optical properties are evaluated either by solving the Bethe-Salpeter equation or by using the adiabatic local density approximation. yambo is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. yambo relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible 110 procedures and is interfaced to several publicly available density functional ground-state codes.

Exciton-Plasmon States in Nanoscale Materials: Breakdown of the Tamm-Dancoff Approximation

Within the Tamm-Dancoff approximation, ab initio approaches describe excitons as packets of electron-hole pairs propagating only forward in time. However, we show that in nanoscale materials excitons and plasmons hybridize, creating exciton-plasmon states where the electron-hole pairs oscillate back and forth in time. Then, as exemplified by the trans-azobenzene molecule and the carbon nanotubes, the Tamm-Dancoff approximation yields errors larger than the accuracy claimed in ab initio calculations. Instead, we propose a general and efficient approach that avoids the Tamm-Dancoff approximation, correctly describes excitons, plasmons, and exciton-plasmon states, and provides a good agreement with experimental results.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.