Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

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.

Integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons

A general graded reflection equation algebra is proposed and the corresponding boundary quantum inverse scattering method is formulated. The formalism is applicable to all boundary lattice systems where an invertible R-matrix exists. As an application, the integrable open-boundary conditions for the q-deformed supersymmetric U model of strongly correlated electrons are investigated. The diagonal boundary K-matrices are found and a class of integrable boundary terms are determined. The boundary system is solved by means of the coordinate space Bethe ansatz technique and the Bethe ansatz equations are derived. As a sideline, it is shown that all R-matrices associated with a quantum affine superalgebra enjoy the crossing-unitarity property. (C) 1998 Elsevier Science B.V.

Role of generalized and episode specific memories in the word frequency effect in recognition

Frequency, recency, and type of prior exposure to very low-and high-frequency words were manipulated in a 3-phase (i.e., familiarization training, study, and test) design. Increasing the frequency with which a definition for a very low-frequency word was provided during familiarization facilitated the word's recognition in both yes-no (Experiment 1) and forced-choice paradigms (Experiment 2). Recognition of very low-frequency words not accompanied by a definition during familiarization first increased, then decreased as familiarization frequency increased (Experiment I). Reasons for these differences were investigated in Experiment 3 using judgments of recency and frequency. Results suggested that prior familiarization of a very low-frequency word with its definition may allow a more adequate episodic representation of the word to be formed during a subsequent study trial. Theoretical implications of these results for current models of memory are discussed.

Observation of generalized synchronization of chaos in a driven chaotic system

We report on the experimental observation of the generalized synchronization of chaos in a real physical system. We show that under a nonlinear resonant interaction, the chaotic dynamics of a single mode laser can become functionally related to that of a chaotic driving signal and furthermore as the coupling strength is further increased, the chaotic dynamics of the laser approaches that of the driving signal.

New integrable boundary conditions for the q-deformed supersymmetric U model and Bethe ansatz equations

New classes of integrable boundary conditions for the q-deformed (or two-parameter) supersymmetric U model are presented. The boundary systems are solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Febrile seizures and generalized epilepsy associated with a mutation in the Na+-channel beta 1 subunit gene SCN1B

Febrile seizures affect approximately 3% of all children under six years of age and are by far the most common seizure disorder(1). A small proportion of children with febrile seizures later develop ongoing epilepsy with afebrile seizures(2). Segregation analysis suggests the majority of cases have complex inheritance(3) but rare families show apparent autosomal dominant: inheritance. Two putative loci have been mapped (FEB1 and FEB2), but specific genes have not yet been identified(4,5). We recently described a clinical subset, termed generalized epilepsy with febrile seizures plus (GEFS(+)), in which many family members have seizures with fever that may persist beyond six years of age or be associated with afebrile generalized seizures(6). We now report linkage, in another large GEFS(+) family, to chromosome region 19q13.1 and identification of a mutation in the voltage-gated sodium (Na+)-channel beta 1 subunit gene (SCN1B). The mutation changes a conserved cysteine residue disrupting a putative disulfide bridge which normally maintains an extracellular immunoglobulin-like fold. Go-expression of the mutant pr subunit with a brain Na+-channel alpha subunit in Xenopus laevis oocytes demonstrates that the mutation interferes with the ability of the subunit to modulate channel-gating kinetics consistent with a loss-of-function allele. This observation develops the theme that idiopathic epilepsies are a family of channelopathies and raises the possibility of involvement of other Na+-channel subunit genes in febrile seizures and generalized epilepsies with complex inheritance patterns.

Generalized optimal lattice covering of finite-dimensional Euclidean space

A generalization of the classical problem of optimal lattice covering of R-n is considered. Solutions to this generalized problem are found in two specific classes of lattices. The global optimal solution of the generalization is found for R-2. (C) 1998 Elsevier Science Inc. All rights reserved.

Risk premiums and benefit measures for generalized-expected-utility theories

Standard tools for the analysis of economic problems involving uncertainty, including risk premiums, certainty equivalents and the notions of absolute and relative risk aversion, are developed without making specific assumptions on functional form beyond the basic requirements of monotonicity, transitivity, continuity, and the presumption that individuals prefer certainty to risk. Individuals are not required to display probabilistic sophistication. The approach relies on the distance and benefit functions to characterize preferences relative to a given state-contingent vector of outcomes. The distance and benefit functions are used to derive absolute and relative risk premiums and to characterize preferences exhibiting constant absolute risk aversion (CARA) and constant relative risk aversion (CRRA). A generalization of the notion of Schur-concavity is presented. If preferences are generalized Schur concave, the absolute and relative risk premiums are generalized Schur convex, and the certainty equivalents are generalized Schur concave.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

Localizing the relativistic electron

A causally well-behaved solution of the localization problem for the free electron is given, with natural space-time transformation properties, in terms of Dirac's position operator x. It is shown that, although x is not an observable in the usual sense, and has no positive-energy (generalized) eigenstates, the four-vector density (rho(x, t), j(x, t)/c) is observable, and can be localized arbitrarily precisely about any point in space, at any instant of time, using only positive energy states. A suitable spin operator can be diagonalized at the same time.

Bounds on integrals of the Wigner function

The integral of the Wigner function over a subregion of the phase space of a quantum system may be less than zero or greater than one. It is shown that for systems with 1 degree of freedom, the problem of determining the best possible upper and lower bounds on such an integral, over an possible states, reduces to the problem of finding the greatest and least eigenvalues of a Hermitian operator corresponding to the subregion. The problem is solved exactly in the case of an arbitrary elliptical region. These bounds provide checks on experimentally measured quasiprobability distributions.

Level-one highest weight representations of U-q[g1(1 vertical bar 1)] and associated vertex operators

We study the level-one irreducible highest weight representations of U-q[gl(1\1)] and associated q-vertex operators. We obtain the exchange relations satisfied by these vertex operators. The characters and supercharacters associated with these irreducible representations are calculated'. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Level-one highest weight representation of U-q[sl((N)over-cap vertical bar 1)] and Bosonization of the multicomponent Super t-J model

We study the level-one irreducible highest weight representations of the quantum affine superalgebra U-q[sl((N) over cap\1)], and calculate their characters and supercharacters. We obtain bosonized q-vertex operators acting on the irreducible U-q[sl((N) over cap\1)] modules and derive the exchange relations satisfied by the vertex operators. We give the bosonization of the multicomponent super t-J model by using the bosonized vertex operators. (C) 2000 American Institute of Physics. [S0022- 2488(00)00508-9].