Approximate VB-MO schemes for excitation gaps in Hubbard models

The excitation gaps in the singlet and triplet manifolds for finite Hubbard models in one, two and three dimensions have been obtained using different approximate configuration interaction (CI) schemes, as a function of the correlation strength, by using valence bond (VB) functions constructed over the molecular orbital (MO) basis. These are compared with numerically exact results and it is found that the scheme in which all particle hole excitations below a given threshold are included is the method of choice. The excitation energies are well reproduced, in trend as well as magnitude, particularly when the threshold equals the bandwidth of the corresponding noninteracting system.

Electron energy loss spectroscopic study of Cr(CO)6, Mo(CO)6 and W(CO)6 in the vapour phase

Electron energy loss spectra (EELS) of Cr, Mo and W hexacarbonyls in the vapour phase are reported. Most of the bands observed are similar to those in optical spectra, but the two high energy transitions in the 9·8–11·2 eV region are reported here for the first time. Based on the orbital energies from the ultraviolet photoelectron spectra and the electronic transition energies from EELS and earlier optical studies, the molecular energy level schemes of these molecules are constructed.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

Molecular structure of electron donor-acceptor complexes of metallotetraphenylporphyrins with trinitrobenzene

The metallotetraphenylporphyrins, MTPPs, where M=Co(II), Cu(II) and Ag(II) form one to one molecular complexes in solution with 1,3,5-trinitrobenzene (TNB). The crystal structure of CoTPP.TNB.2CH3OH revealed that the mean separation between the porphyrin and TNB planes is 3·27 Å and the centre of the aryl ring of TNB is displaced by 0·90 Å from the centre of the porphyrin plane. Extended Huckel molecular orbital calculations suggest that acceptor orbitals are predominantly nitro-group based rather than an aryl π framework. The contribution of the metal orbitals in the donation is also seen as predicted by the structural data. Electrochemical redox measurements in solution provide evidence for charge-transfer stabilisation in these complexes.

NDDO-based CI methods for the prediction of electronic spectra and sum-over-states molecular hyperpolarization

NDDO-based (AM1) configuration interaction (CI) calculations have been used to calculate the wavelength and oscillator strengths of electronic absorptions in organic molecules and the results used in a sum-over-states treatment to calculate second-order-hyperpolarizabilities. The results for both spectra and hyperpolarizabilities are of acceptable quality as long as a suitable CI-expansion is used. We have found that using an active space of eight electrons in eight orbitals and including all single and pair-double excitations in the CI leads to results that agree well with experiment and that do not change significantly with increasing active space for most organic molecules. Calculated second-order hyperpolarizabilities using this type of CI within a sum-over-states calculation appear to be of useful accuracy.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.

Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Ab initio GW calculations are a standard method for computing the spectroscopic properties of many materials. The most computationally expensive part in conventional implementations of the method is the generation and summation over the large number of empty orbitals required to converge the electron self-energy. We propose a scheme to reduce the summation over empty states by the use of a modified static remainder approximation, which is simple to implement and yields accurate self-energies for both bulk and molecular systems requiring a small fraction of the typical number of empty orbitals.

Dislocations emission and crack extension at the atomistic crack tip in body-centered-cubic metal Mo

The behaviors of a crack in body-centered-cubic metal Mo under different loading modes were studied using the molecular dynamics method. Dislocation emission was observed near the crack tip in response to mode II loading with theta = 0 degrees in which theta is the inclination angle of the slip plane with respect to the crack plane, and two full dislocations were observed at the stress level of K-II = 1.17 MPa m(1/2) without any evidence of crack extension. Within the range of 0 degrees less than or equal to theta less than or equal to 45 degrees, crack extension was observed in response to mode I loading, and the effect of crystal orientation on the crack propagation was studied, The crack propagated along the [111] slip direction without any evidence of dislocations emission.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

Variable-angle electron spectroscopic studies of various polyatomic molecular systems

The complementary techniques of low-energy, variable-angle electron-impact spectroscopy and ultraviolet variable-angle photoelectron spectroscopy have been used to study the electronic spectroscopy and structure of several series of molecules. Electron-impact studies were performed at incident beam energies between 25 eV and 100 eV and at scattering angles ranging from 0° to 90°. The energy-loss regions from 0 eV to greater than 15 eV were studied. Photoelectron spectroscopic studies were conducted using a HeI radiation source and spectra were measured at scattering angles from 45° to 90°. The molecules studied were chosen because of their spectroscopic, chemical, and structural interest. The operation of a new electron-impact spectrometer with multiple-mode target source capability is described. This spectrometer has been used to investigate the spin-forbidden transitions in a number of molecular systems.

The electron-impact spectroscopy of the six chloro-substituted ethylenes has been studied over the energy-loss region from 0-15 eV. Spin-forbidden excitations corresponding to the π → π*, N → T transition have been observed at excitation energies ranging from 4.13 eV in vinyl chloride to 3.54 eV in tetrachloroethylene. Symmetry-forbidden transitions of the type π → np have been oberved in trans-dichloroethyene and tetrachlor oethylene. In addition, transitions to many states lying above the first ionization potential were observed for the first time. Many of these bands have been assigned to Rydberg series converging to higher ionization potentials. The trends observed in the measured transition energies for the π → π*, N → T, and N → V as well as the π → 3s excitation are discussed and compared to those observed in the methyl- and fluoro- substituted ethylenes.

The electron energy-loss spectra of the group VIb transition metal hexacarbonyls have been studied in the 0 eV to 15 eV region. The differential cross sections were obtained for several features in the 3-7 eV energy-loss region. The symmetry-forbidden nature of the ¹A_1g → ¹A_1g, 2t_2g(π) → 3t_2g(π*) transition in these compounds was confirmed by the high-energy, low-angle behavior of their relative intensities. Several low lying transitions have been assigned to ligand field transitions on the basis of the energy and angular behavior of the differential cross sections for these transitions. No transitions which could clearly be assigned to singlet → triplet excitations involving metal orbitals were located. A number of states lying above the first ionization potential have been observed for the first time. A number of features in the 6-14 eV energy-loss region of the spectra of these compounds correspond quite well to those observed in free CO.

A number of exploratory studies have been performed. The π → π*, N → T, singlet → triplet excitation has been located in vinyl bromide at 4.05 eV. We have also observed this transition at approximately 3.8 eV in a cis-/trans- mixture of the 1,2-dibromoethylenes. The low-angle spectrum of iron pentacarbonyl was measured over the energy-loss region extending from 2-12 eV. A number of transitions of 8 eV or greater excitation energy were observed for the first time. Cyclopropane was also studied at both high and low angles but no clear evidence for any spin- forbidden transitions was found. The electron-impact spectrum of the methyl radical resulting from the pyrolysis of tetramethyl tin was obtained at 100 eV incident energy and at 0° scattering angle. Transitions observed at 5.70 eV and 8.30 eV agree well with the previous optical results. In addition, a number of bands were observed in the 8-14 eV region which are most likely due to Rydberg transitions converging to the higher ionization potentials of this molecule. This is the first reported electron-impact spectrum of a polyatomic free radical.

Variable-angle photoelectron spectroscopic studies were performed on a series of three-membered-ring heterocyclic compounds. These compounds are of great interest due to their highly unusual structure. Photoelectron angular distributions using HeI radiation have been measured for the first time for ethylene oxide and ethyleneimine. The measured anisotropy parameters, β, along with those measured for cyclopropane were used to confirm the orbital correlations and photoelectron band assignments. No high values of β similar to those expected for alkene π orbitals were observed for the Walsh or Forster-Coulson-Moffit type orbitals.

I. Electronic states of heavily-doped molecular crystals--naphthalene. I. Theoretical. II. Experimental. II. Lowest singlet and triplet excited states of pyrazine

PART I

The energy spectrum of heavily-doped molecular crystals was treated in the Green’s function formulation. The mixed crystal Green’s function was obtained by averaging over all possible impurity distributions. The resulting Green’s function, which takes the form of an infinite perturbation expansion, was further approximated by a closed form suitable for numerical calculations. The density-of-states functions and optical spectra for binary mixtures of normal naphthalene and deuterated naphthalene were calculated using the pure crystal density-of-state functions. The results showed that when the trap depth is large, two separate energy bands persist, but when the trap depth is small only a single band exists. Furthermore, in the former case it was found that the intensities of the outer Davydov bands are enhanced whereas the inner bands are weakened. Comparisons with previous theoretical calculations and experimental results are also made.

PART II

The energy states and optical spectra of heavily-doped mixed crystals are investigated. Studies are made for the following binary systems: (1) naphthalene-h₈ and d₈, (2) naphthalene--h₈ and αd₄, and (3) naphthalene--h₈ and βd₁, corresponding to strong, medium and weak perturbations. In addition to ordinary absorption spectra at 4˚K, band-to-band transitions at both 4˚K and 77˚K are also analyzed with emphasis on their relations to cooperative excitation and overall density-of-states functions for mixed crystals. It is found that the theoretical calculations presented in a previous paper agree generally with experiments except for cluster states observed in system (1) at lower guest concentrations. These features are discussed semi-quantitatively. As to the intermolecular interaction parameters, it is found that experimental results compare favorably with calculations based on experimental density-of-states functions but not with those based on octopole interactions or charge-transfer interactions. Previous experimental results of Sheka and the theoretical model of Broude and Rashba are also compared with present investigations.

PART III

The phosphorescence, fluorescence and absorption spectra of pyrazine-h₄ and d₄ have been obtained at 4˚K in a benzene matrix. For comparison, those of the isotopically mixed crystal pyrazine-h₄ in d₄ were also taken. All these spectra show extremely sharp and well-resolved lines and reveal detailed vibronic structure.

The analysis of the weak fluorescence spectrum resolves the long-disputed question of whether one or two transitions are involved in the near-ultraviolet absorption of pyrazine. The “mirror-image relationship” between absorption and emission shows that the lowest singlet state is an allowed transition, properly designated as ¹B_3u ← ¹A_1g. The forbidden component ¹B_2g, predicted by both “exciton” and MO theories to be below the allowed component, must lie higher. Its exact location still remains uncertain.

The phosphorescence spectrum when compared with the excitation phosphorescence spectra, indicates that the lowest triplet state is also symmetry allowed, showing a strong 0-0 band and a “mirror-image relationship” between absorption and emission. In accordance with previous work, the triplet state is designated as ³B_3u.

The vibronic structure of the phosphorescence spectrum is very complicated. Previous work on the analysis of this spectrum all concluded that a long progression of v_6a exists. Under the high resolution attainable in our work, the supposed v_6a progression proves to have a composite triplet structure, starting from the second member of the progression. Not only is the v_9a hydrogen-bending mode present as shown by the appearance of the C-D bending mode in the d₄ spectrum, but a band of 1207 cm^-1 in the pyrazine in benzene system and 1231 cm^-1 in the mixed crystal system is also observed. This band is assigned as 2v_6b and of a_1g symmetry. Its anonymously strong intensity in the phosphorescence spectrum is interpreted as due to the Fermi resonance with the 2v_6a and v_9a band.

To help resolve the present controversy over the crystal phosphorescence spectrum of pyrazine, detailed vibrational analyses of the emission spectra were made. The fluorescence spectrum has essentially the same vibronic structure as the phosphorescence spectrum.

ELECTRON-STATES OF A STACKING-FAULT RIBBON IN SILICON

The electronic structure of a bounded intrinsic stacking fault in silicon is calculated. The method used is an LCAO-scheme (Linear Combinations of Atomic Orbitals) taking ten atomic orbitals of s-, p-, and d-type into account. The levels in the band gap are extracted using Lanczos' algorithm and a continued fraction representation of the local density of states. We find occupied states located up to 0.3 eV above the valence band maximum (E(v)). This significantly differs from the result obtained for the ideal infinite fault for which the interface state is located at E(v)+ 0.1 eV.

INFLUENCE OF THE 90-DEGREES PARTIAL DISLOCATION CORE STRUCTURE IN SILICON ON ENERGY-LEVELS

An LCAO-scheme taking into account 10 atomic orbitals (s-, p-, and d-type) is used to calculate the electronic structure of the reconstructed 90-degrees partial dislocation in Si. Two different valence force fields producing deviating results are used for modelling the core structure. Geometrical data published by another group is also used. The aim is to explore the influence of geometry on energy levels. We find that the band structure depends sensitively on bond angles. Using data determined by the Tersoff potential we obtain two bands of which the upper one penetrates deeply into the indirect band gap while the geometry minimizing the simple Keating potential leaves the gap completely clear of dislocation states. Thus, from a theoretical point of view, the chief difficulty in calculating the electronic structure of the reconstructed 90-degrees partial is the lack of accurate structural information.

ENERGY-LEVEL CALCULATIONS OF THE RECONSTRUCTED 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme taking into account 10 atomic orbitals (s-, p-, and d-type) applied to a supercell containing 256 atoms is used to calculate the bound states of the reconstructed 90-degrees partial dislocation in Si. The results differ significantly from our earlier calculations on the unreconstructed 90-degrees partial using the same method. We find two bands separate from each other in the entire Brillouin zone and the upper band penetrates deep into the indirect band gap which is in contradiction with the general opinion that core reconstruction clears the band gap of dislocation states.

ELECTRON-STATES OF A VACANCY IN THE CORE OF THE 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme (linear combination of atomic orbitals) taking into account ten atomic orbitals (s-, p-, and d-type) is used to calculate the electronic structure of a vacancy present in the core of the reconstructed 90 degrees partial dislocation in silicon. The levels in the band gap are extracted using Lanczos' algorithm and a continued fraction representation of the local density of states. The three-fold degenerate stale of the ideal vacancy is split into three levels with energies 0.26, 1.1, and 1.9 eV measured from the valence band edge.