Variational Calculation of the Hydrogen Molecular Cation Using Maple

The variational calculation of the energy of the Hydrogen Molecular (H2+) Cation's LCAO-MO for the sigma and sigma* states as functions of the AO's screening constant and the internuclear distance is carried out explicitly in great detail.

LCAO-MO's are not Eigenfunctions of the H_2^+ Hamiltonian

A simple LCAO-MO for the hydrogen molecule cation is tested for eigenfunctionality and found to be flawed.

Variational Calculation of the Hydrogen Molecular Cation using Maple (II). The pi and pi-star Orbitals

The pi and pi-star orbitals of the hydrogen molecular cation are obtained using Maple in the same manner as the sigma and sigma-star orbitals were obtained in paper-36.

The second dissociation threshold bound levels of hydrogen molecule

The near-threshold highly bound states of all three stable isotopic variants of molecular hydrogen have been studied. Numerous perturbations and unexpected transitions are observed as far as 1cm(-1) just below the second dissociation threshold. This complex structure may arise from a combination of nonadiabatic coupling between B, B', C electronic states, perturbations due to. ne and hyperfine interactions, and strong shape resonances. The perturbed near-threshold states and vibrational continuum exhibit finegrained structure, differing greatly between isotopes because of varying nonadiabatic coupling.

Positronium impact excitation of hydrogen molecule to B-1 Sigma(+)(u) and b(3)Sigma(+)(u) states

Calculation for the electronic excitation of the ground state of H-2 to B (1) Sigma(u)(+) and b(3) Sigma(u)(+) states by positronium- (Ps) atom impact has been carried out using the first Born approximation considering discrete Ps excitations up to n = 6 and Ps ionization in the final state. To include the effect of electron exchange, we propose an alternative approximation scheme in the light of the Rudge approach, which takes into account the composite nature of the Ps-atom projectile.

Positronium scattering by a hydrogen molecule including exchange

We perform a three-positronium (Ps) state [Ps(ls,2s,2p)] coupled-channel calculation of Ps-H-2 scattering including the effect of electron exchange. At medium energies, higher excitations and ionization of Ps are treated within the framework of the first Born approximation. In both cases exchange is included using a recently proposed nonlocal model exchange potential which is free of non-orthogonality problems common in the usual antisymmetrization scheme. The present total cross sections at low and medium energies are in encouraging agreement with experiment.

Measurement of the conductance of a hydrogen molecule

Recent years have shown steady progress towards molecular electronics, in which molecules form basic components such as switches, diodes and electronic mixers. Often, a scanning tunnelling microscope is used to address an individual molecule, although this arrangement does not provide long-term stability. Therefore, metal–molecule–metal links using break-junction devices have also been explored; however, it is difficult to establish unambiguously that a single molecule forms the contact. Here we show that a single hydrogen molecule can form a stable bridge between platinum electrodes. In contrast to results for organic molecules, the bridge has a nearly perfect conductance of one quantum unit, carried by a single channel. The hydrogen bridge represents a simple test system in which to understand fundamental transport properties of single-molecule devices.

Hydrogen-bonding in the structures of the 1:1 salts of isonipecotamide with a set of four polyfunctional monocyclic heteroaromatic carboxylic acids

The crystal structures of the 1:1 proton-transfer compounds of isonipecotamide (piperidine-4-carboxamide) with the monocyclic heteroaromatic carboxylic acids, isonicotinic acid, picolinic acid, dipicolinic acid and pyrazine-2,3-dicarboxylic acid have been determined at 200 K and their hydrogen-bonding patterns examined. The compounds are respectively anhydrous 4-carbamoylpiperidinium pyridine-4-carboxylate (1), the partial hydrate 4-carbamoylpiperidinium pyridine-2-carboxylate 0.25 water (2), the solvate 4-carbamoylpiperidinium 6-carboxypyridine-2-carboxylate methanol monosolvate (3), and anhydrous 4-carbamoylpiperidinium 3-carboxypyrazine-2-carboxylate (4). In compounds 1 and 3, hydrogen-bonding interactions give two-dimensional sheet structures which feature enlarged cyclic ring systems, while in compounds 2 and 4, three-dimensional structures are found. The previously described cyclic R2/2(8) hydrogen-bonded amide-amide dimer is present in 2 and 3. The hydrogen-bonding in 2 involves the partial-occupancy water molecule while the structure of 4 is based on inter-linked homomolecular hydrogen-bonded cation-cation and anion-anion associated chains comprising head-to-tail interactions. This work further demonstrates the utility of the isonipecotamide cation in the generation of chemically stable hydrogen-bonded systems, particularly with aromatic carboxylate anions, providing crystalline solids.

MO visualization of Pi Orbitals

In the spirit of trying to convert people to understanding atomic orbitals centered elsewhere than the origin, we continue the discussion of visualizing molecular orbitals, so called LCAO-MO, using various plotting tricks in Maple.

The role of Ti as a catalyst for the dissociation of hydrogen on a Mg(0001) surface

In this paper, the dissociative chemisorption of hydrogen on both pure and Ti-incorporated Mg(0001) surfaces are studied by ab initio density functional theory (DFT) calculations. The calculated dissociation barrier of hydrogen molecule on a pure Mg(0001) surface (1.05 eV) is in good agreement with comparable theoretical studies. For the Ti-incorporated Mg(0001) surface, the activated barrier decreases to 0.103 eV due to the strong interaction between the molecular orbital of hydrogen and the d metal state of Ti. This could explain the experimentally observed improvement in absorption kinetics of hydrogen when transition metals have been introduced into the magnesium materials.

NMR (1H, 13C) and MO (ab initio, CNDO/2, EHT) studies on basic and N-protonated hydrazinecarbothioamide

The conformational preferences of hydrazinecarbothioamide (HCTA, H2NNHCSNH2) in its basic and N-protonated (PHCTA, H3NNNHCSNH2) forms have been studied by 1H and 13C NMR spectroscopy and by theoretical LCAO-MO methods (ab initio, CNDO/2 and EHT). The hindered rotation around the C---N bond has been investigated by a total line shape analysis for the thioamide protons and by the three MO methods. Changes in the molecular conformation and electronic structure on protonation are briefly discussed.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Continuum structures at the second dissociation limit of hydrogen molecules

We describe near-threshold high-resolution spectra and continuum resonance dynamical behaviour of all three stable hydrogen isotopic variants, and finally obtain improved values for the dissociation energies of hydrogen molecule and its ion. The second dissociation limit is determined by analysing the onset of the vibrational continuum.

Hydrogen adsorption in a copper ZSM5 zeolite: an inelastic neutron scattering study

Currently microporous oxidic materials including zeolites are attracting interest as potential hydrogen storage materials. Understanding how molecular hydrogen interacts with these materials is important in the rational development of hydrogen storage materials and is also challenging theoretically. In this paper, we present an incoherent inelastic neutron scattering (INS) study of the adsorption of molecular hydrogen and hydrogen deuteride (HD) in a copper substituted ZSM5 zeolite varying the hydrogen dosage and temperature. We have demonstrated how inelastic neutron scattering can help us understand the interaction of H-2 molecules with a binding site in a particular microporous material, Cu ZSM5, and by implication of other similar materials. The H-2 molecule is bound as a single species lying parallel with the surface. As H-2 dosing increases, lateral interactions between the adsorbed H-2 molecules become apparent. With rising temperature of measurement up to 70 K (the limit of our experiments), H-2 molecules remain bound to the surface equivalent to a liquid or solid H-2 phase. The implication is that hydrogen is bound rather strongly in Cu ZSM5. Using the simple model for the anisotropic interaction to calculate the energy levels splitting, we found that the measured rotational constant of the hydrogen molecule is reduced as a consequence of adsorption by the Cu ZSM5. From the decrease in total signal intensity with increasing temperature, we were able to observe the conversion of para-hydrogen into ortho-hydrogen at paramagnetic centres and so determine the fraction of paramagnetic sites occupied by hydrogen molecules, ca. 60%. (c) 2006 Elsevier B.V. All rights reserved.