Matrix elements of U(2n) generators in a multishell spin-orbit basis. I. The del-operator MEs in a two-shell composite Gelfand-Paldus basis

This is the first in a series of three articles which aimed to derive the matrix elements of the U(2n) generators in a multishell spin-orbit basis. This is a basis appropriate to many-electron systems which have a natural partitioning of the orbital space and where also spin-dependent terms are included in the Hamiltonian. The method is based on a new spin-dependent unitary group approach to the many-electron correlation problem due to Gould and Paldus [M. D. Gould and J. Paldus, J. Chem. Phys. 92, 7394, (1990)]. In this approach, the matrix elements of the U(2n) generators in the U(n) x U(2)-adapted electronic Gelfand basis are determined by the matrix elements of a single Ll(n) adjoint tensor operator called the del-operator, denoted by Delta(j)(i) (1 less than or equal to i, j less than or equal to n). Delta or del is a polynomial of degree two in the U(n) matrix E = [E-j(i)]. The approach of Gould and Paldus is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. Hence, to generalize this approach, we need to obtain formulas for the complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The nonzero shift coefficients are uniquely determined and may he evaluated by the methods of Gould et al. [see the above reference]. In this article, we define zero-shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis which are appropriate to the many-electron problem. By definition, these are proportional to the corresponding two-shell del-operator matrix elements, and it is shown that the Racah factorization lemma applies. Formulas for these coefficients are then obtained by application of the Racah factorization lemma. The zero-shift adjoint reduced Wigner coefficients required for this procedure are evaluated first. All these coefficients are needed later for the multishell case, which leads directly to the two-shell del-operator matrix elements. Finally, we discuss an application to charge and spin densities in a two-shell molecular system. (C) 1998 John Wiley & Sons.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. III. General formulas

This is the third and final article in a series directed toward the evaluation of the U(2n) generator matrix elements (MEs) in a multishell spin/orbit basis. Such a basis is required for many-electron systems possessing a partitioned orbital space and where spin-dependence is important. The approach taken is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. A complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis (which is appropriate to the many-electron problem) were obtained in the first and second articles of this series. Ln the first article we defined zero-shift coupling coefficients. These are proportional to the corresponding two-shell del-operator matrix elements. See P. J. Burton and and M. D. Gould, J. Chem. Phys., 104, 5112 (1996), for a discussion of the del-operator and its properties. Ln the second article of the series, the nonzero shift coupling coefficients were derived. Having obtained all the necessary coefficients, we now apply the formalism developed above to obtain the U(2n) generator MEs in a multishell spin-orbit basis. The methods used are based on the work of Gould et al. (see the above reference). (C) 1998 John Wiley & Sons, Inc.

Violation of Kohler's rule by the magnetoresistance of a quasi-two-dimensional organic metal

The interlayer magnetoresistance of the quasi-two-dimensional metal alpha-(BEDT-TTF)(2)KHg(SCN)(4) is considered. In the temperature range from 0.5 to 10 K and for fields up to 10 T the magnetoresistance has a stronger temperature dependence than the zero-field resistance. Consequently Kohler's rule is not obeyed for any range of temperatures or fields. This means that the magnetoresistance cannot be described in terms of semiclassical transport on a single Fermi surface with a single scattering time. Possible explanations for the violations of Kohler's rule are considered, both within the framework of semiclassical transport theory and involving incoherent interlayer transport. The issues considered are similar to those raised by the magnetotransport of the cuprate superconductors. [S0163-1829(98)13219-8].

Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

N-mesityl-C-acylketenimines: 1,5-sigmatropic shifts and electrocyclization to quinolines

Flash vacuum thermolysis (FVT) of triazoles 6a-c generates alpha-oxoketenimines 10, the ester 10a being isolable. FVT of pyrroledione 8 generates the isomeric imidoylketene 9a. Ketenes 9 and ketenimines 10 undergo thermal interconversion by 1,3-shifts of methoxy and dimethylamino groups under mild FVT conditions (ca. 350-400 degrees C). Both 9 and 10 are directly observable by IR spectroscopy at either 77 K or on Ar matrix isolation at 12 K. On FVT at temperatures above ca. 400 degrees C, the ketenimines 10 undergo a 1,5-H shift to o-quinoid imines 12/13, followed by electrocyclization to dihydroquinolines 14 (unobserved) and 15 (observed by NMR). The latter are easily oxidized to alkylquinoline-3-carboxylates or quinoline-3-carboxamides 16 by atmospheric oxygen. Ab initio calculations on model compounds 18-23 predict an energy barrier of ca. 38 kcal mol(-1) (161 kJ mol(-1)) for the 1,5-H shift in N-(o-methylphenyl)ketenimines via the transition state TS19 followed by an electrocyclization barrier to dihydroquinoline 23a via TS22a of ca. 16 kcal mol(-1).

Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Comparison of coherent and weakly incoherent transport models for the interlayer magnetoresistance of layered Fermi liquids

The interlayer magnetoresistance of layered metals in a tilted magnetic field is calculated for two distinct models for the interlayer transport. The first model involves coherent interlayer transport, and makes use of results of semiclassical or Bloch-Boltzmann transport theory. The second model involves weakly incoherent interlayer transport where the electron is scattered many times within a layer before tunneling into the next layer. The results are relevant to the interpretation of experiments on angular-dependent magnetoresistance oscillations (AMRO) in quasi-one- and quasi-two-dimensional organic metals. We find that the dependence of the magnetoresistance on the direction of the magnetic field is identical for both models except when the field is almost parallel to the layers. An important implication of this result is that a three-dimensional Fermi surface is not necessary for the observation of the Yamaji and Danner oscillations seen in quasi-two- and quasi-one-dimensional metals, respectively. A universal expression is given for the dependence of the resistance at AMRO maxima and minima on the magnetic field and scattering time (and thus the temperature). We point out three distinctive features of coherent interlayer transport: (i) a beat frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angular-dependent magnetoresistance when the field is sufficiently large and parallel to the layers, and (iii) a crossover from a linear to a quadratic field dependence for the magnetoresistance when the field is parallel to the layers. Properties (i) and (ii) are compared with published experimental data for a range of quasi-two-dimensional organic metals. [S0163-1829(99)02236-5].

Periodic orbit resonances in layered metals in tilted magnetic fields

The frequency dependence of the interlayer conductivity of a layered Fermi liquid in a magnetic field that is tilted away from the normal to the layers is considered. For both quasi-one- and quasi-two-dimensional systems resonances occur when the frequency is a harmonic of the frequency at which the magnetic field causes the electrons to oscillate on the Fermi surface within the layers. The intensity of the different harmonic resonances varies significantly with the direction of the field. The resonances occur for both coherent and weakly incoherent interlayer transport and so their observation does not imply the existence of a three-dimensional Fermi surface. [S0163-1829(99)51240-X].

Novel vascular graft grown within recipient's own peritoneal cavity

A method by which to overcome the clinical symptoms of atherosclerosis is the insertion of a graft to bypass an artery blocked or impeded by plaque. However, there may be insufficient autologous mammary artery for multiple or repeat bypass, saphenous vein may have varicose degenerative alterations that can lead to aneurysm in high-pressure sites, and small-caliber synthetic grafts are prone to thrombus induction and occlusion. Therefore, the aim of the present study was to develop an artificial blood conduit of any required length and diameter from the cells of the host for autologous transplantation. Silastic tubing, of variable length and diameter, was inserted into the peritoneal cavity of rats or rabbits. By 2 weeks, it had become covered by several layers of myofibroblasts, collagen matrix, and a single layer of mesothelium. The Silastic tubing was removed from the harvested implants, and the tube of living tissue was everted such that it now resembled a blood vessel with an inner lining of nonthrombotic mesothelial cells (the intima), with a media of smooth muscle-like cells (myofibroblasts), collagen, and elastin, and with an outer collagenous adventitia. The tube of tissue (10 to 20 mm long) was successfully grafted by end-to-end anastomoses into the severed carotid artery or abdominal aorta of the same animal in which they were grown. The transplant remained patent for at least 4 months and developed structures resembling elastic lamellae. The myofibroblasts gained a higher volume fraction of myofilaments and became responsive to contractile agonists, similar to the vessel into which they had been grafted. It is suggested that these nonthrombogenic tubes of living tissue, grown in the peritoneal cavity of the host, may be developed as autologous coronary artery bypass grafts or as arteriovenous access fistulae for hemodialysis patients.

Dynamic spatiotemporal expression patterns of neurocan and phosphacan indicate diverse roles in the developing and adult mouse olfactory system

The chondroitin sulfate proteoglycans neurocan and phosphacan are believed to modulate neurite outgrowth by binding to cell adhesion molecules, tenascin, and the differentiation factors heparin-binding growth-associated molecule and amphoterin. To assess the role of these chondroitin sulfate proteoglycans in the olfactory system, we describe here their expression patterns during both embryonic and postnatal development in the mouse. Immunoreactivity for neurocan was first detected in primary olfactory neurons at embryonic day 11.5 (E11.5). Neurocan was expressed by primary olfactory axons as they extended toward the rostral pole of the telencephalon as well as by their arbors in glomeruli after they contacted the olfactory bulb. The role of neurocan was examined by growing olfactory neurons on an extracellular matrix substrate containing neurocan or on extracellular matrix in the presence of soluble neurocan. In both cases, neurocan strongly promoted neurite outgrowth. These results suggest that neurocan supports the growth of primary olfactory axons through the extracellular matrix as they project to the olfactory bulb during development. Phosphacan, unlike neurocan, was present within the mesenchyme surrounding the E11.5 and E12.5 nasal cavity. This expression decreased at E13.5, concomitant with a transient appearance of phosphacan in nerve fascicles. Within the embryonic olfactory bulb, phosphacan was localised to the external and internal plexiform layers. However, during early postnatal development phosphacan was concentrated in the glomerular layer. These results suggest that phosphacan may play a role in delineating the pathway of growing olfactory axons as well as defining the laminar organization of the bulb. Together, the spatiotemporal expression patterns of neurocan and phosphacan indicate that these chondroitin sulfate proteoglycans have diverse in situ roles, which are dependent on context-specific interactions with extracellular and cell adhesion molecules within the developing olfactory nerve pathway. (C) 2000 Wiley-Liss, Inc.

Blood vessels from bone marrow

Lengths of silastic tubing were inserted into the peritoneal cavity of rats or rabbits. By two weeks the free-floating implants had become covered by a capsule consisting of several layers of macrophage-derived myofibroblasts and collagen matrix overlaid by a single layer of mesothelial cells. The tubing was removed from the harvested implant and the tissue everted. This now resembled an artery with an inner lining of mesothelial cells (the intima), a media of myofibroblasts, and an outer collagenous adventitia. The tube of living tissue was grafted by end-to-end anastomoses into the transected carotid artery or abdominal aorta of the same animal in which the tissue had been grown, where it remained parent for four months and developed structures resembling elastic lamellae, The myofibroblasts developed a high volume fraction of myofilaments and became responsive to contractile and relaxing agents similar to smooth muscle cells of the adjacent artery wall.