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.

A small disturbance in the strong vortex flow

A small disturbance in the axisymmetric, bathtub-like flow with strong vorticity is considered and the asymptotic representation of the solution is found. It is shown that if the disturbance is smaller than a certain critical scale, the conventional strong vortex approximation cannot describe the field generated by the disturbance not only in the vicinity of the disturbance but also at the distances much larger than the critical scale. (C) 2001 American Institute of Physics.

Simulation of a complete reflected shock tunnel showing a vortex mechanism for flow contamination

Simulations of a complete reflected shock tunnel facility have been performed with the aim of providing a better understanding of the flow through these facilities. In particular, the analysis is focused on the premature contamination of the test flow with the driver gas. The axisymmetric simulations model the full geometry of the shock tunnel and incorporate an iris-based model of the primary diaphragm rupture mechanics, an ideal secondary diaphragm and account for turbulence in the shock tube boundary layer with the Baldwin-Lomax eddy viscosity model. Two operating conditions were examined: one resulting in an over-tailored mode of operation and the other resulting in approximately tailored operation. The accuracy of the simulations is assessed through comparison with experimental measurements of static pressure, pitot pressure and stagnation temperature. It is shown that the widely-accepted driver gas contamination mechanism in which driver gas 'jets' along the walls through action of the bifurcated foot of the reflected shock, does not directly transport the driver gas to the nozzle at these conditions. Instead, driver gas laden vortices are generated by the bifurcated reflected shock. These vortices prevent jetting of the driver gas along the walls and convect driver gas away from the shock tube wall and downstream into the nozzle. Additional vorticity generated by the interaction of the reflected shock and the contact surface enhances the process in the over-tailored case. However, the basic mechanism appears to operate in a similar way for both the over-tailored and the approximately tailored conditions.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Interference in hyperbolic space

The interference in a phase space algorithm of Schleich and Wheeler [Nature 326, 574 (1987)] is extended to the hyperbolic space underlying the group SU(1,1). The extension involves introducing the notion of weighted areas. Analytic expressions for the asymptotic forms for overlaps between the eigenstates of the generators of su(1,1) thus obtained are found to be in excellent agreement with the numerical results.[S1050-2947(98)08602-8].

Heat flow for the Yang-Mills-Higgs field and the Hermitian Yang-Mills-Higgs metric

For a parameter lambda > 0, we study a type of vortex equations, which generalize the well-known Hermitian-Einstein equation, for a connection A and a section phi of a holomorphic vector bundle E over a Kahler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang-Mills-Higgs field on E. Assuming the lambda -stability of (E, phi), we prove the existence of the Hermitian Yang-Mills-Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.

Casimir invariants and characteristic identities for gl(infinity)

A full set of (higher-order) Casimir invariants for the Lie algebra gl(infinity) is constructed and shown to be well defined in the category O-FS generated by the highest weight (unitarizable) irreducible representations with only a finite number of nonzero weight components. Moreover, the eigenvalues of these Casimir invariants are determined explicitly in terms of the highest weight. Characteristic identities satisfied by certain (infinite) matrices with entries from gl(infinity) are also determined and generalize those previously obtained for gl(n) by Bracken and Green [A. J. Bracken and H. S. Green, J. Math. Phys. 12, 2099 (1971); H. S. Green, ibid. 12, 2106 (1971)]. (C) 1997 American Institute of Physics.

Super-Yangian double and its central extension

We introduce the super-Yangian double DY(h over bar), [gl(m/n)] and its central extension DY(h over bar)[gl(m/n)]. We give their defining relations in terms of current generators and obtain Drinfeld co-multiplication. (C) 1997 Elsevier Science B.V.

Topological charge and angular momentum of light beams carrying optical vortices

We analyze the properties of light beams carrying phase singularities, or optical vortices. The transformations of topological charge during free-space propagation of a light wave, which is a combination of a Gaussian beam and a multiple charged optical vortex within a Gaussian envelope, are studied both in theory and experiment. We revise the existing knowledge about topological charge conservation, and demonstrate possible scenarios where additional vortices appear or annihilate during free propagation of such a combined beam. Coaxial interference of optical vortices is also analyzed, and the general rule for angular-momentum density distribution in a combined beam is established. We show that, in spite of any variation in the number of vortices in a combined beam, the total angular momentum is constant during the propagation. [S1050-2947(97)09910-1].

Comments on the Drinfeld realization of the quantum affine superalgebra U-q[gl(m backslash n)(1)] and its Hopf algebra structure

By generalizing the Reshetikhin and Semenov-Tian-Shansky construction to supersymmetric cases, we obtain the Drinfeld current realization for the quantum affine superalgebra U-q[gl(m\n)((1))]. We find a simple coproduct for the quantum current generators and establish the Hopf algebra structure of this super current algebra.

The use and abuse of population viability analysis

A recent study by Brook ef al. empirically tested the performance of population viability analysis (PVA) using data from 21 populations across a wide range of species. The study concluded that PVAs are good at predicting the future dynamics of populations. We suggest that this conclusion is a result of a bias in the studies that Brook et al, included in their analyses, We present arguments that PVAs can only be accurate at predicting extinction probabilities if data are extensive and reliable, and if the distribution of vital rates between individuals and years can be assumed stationary in the future, or if any changes can be accurately predicted. In particular, we note th at although catastrophes are likely to have precipitated many extinctions, estimates of the probability of catastrophes are unreliable.

Temporal relationship between the auditory brainstem response and focal responses of auditory nerve root and cochlear nucleus during development in the tammar wallaby (Macropus eugenii)

Thirty-two pouch-young tammar wallabies were used to discover the generators of the auditory brainstem response (ABR) during development by the use of simultaneous ABR and focal brainstem recordings. A click response from the auditory nerve root (ANR) in the wallaby was recorded from postnatal day (PND) 101, when no central auditory station was functional, and coincided with the ABR, a simple positive wave. The response of the cochlear nucleus (CN) was detected from PND 110, when the ABR had developed 1 positive and 1 negative peak. The dominant component of the focal ANR response, the N-1 wave, coincided with the first half of the ABR P wave, and that of the focal CN response, the N-1 wave, coincided with the later two thirds. In older animals, the ANR response coincided with the ABR's N-1, wave, while the CN response coincided with the ABR's P-2, N-2 and P-3 waves, with its contribution to the ABR P-2 dominant. The protracted development of the marsupial auditory system which facilitated these correlations makes the tammar wallaby a particularly suitable model. Copyright (C) 2001 S. Karger AG, Basel.

Moderately strong vorticity in a bathtub-type flow

Theoretical and numerical analysis is performed for an inviscid axisymmetric vortical bathtub-type flow. The level of vorticity is kept high so that the image of the flow on the radial-axial plane (r-z plane) is not potential. The most significant findings are: (1) the region of validity of the strong vortex approximation is separated from the drain by a buffer region, (2) the power-law asymptote of the stream function, specified by Delta Psi similar to r(4/3) Deltaz, appears near the axis when vorticity in the flow is sufficiently strong and (3) the local Rossby number in the region of the 4/3 power-law the initial vorticity level in the flow and the global Rossby number.