ALBERT

Notes: Abstract This paper is a sequel to a paper by the author entitled “Restricted Transition Probabilities and Their Applications to Some Problems in the Dynamics of Biological Populations” (Bull. Math. Biophysics, 1966,28, 315–331). The paper is divided into two parts. In part one some aspects of the maximum size attained by the population during a finite time interval are studied for the case the stochastic process underlying the evolution of the population is a birth process. Two interesting by-products emerge from the study presented in part one; namely a combinatorial method of finding solutions to the Kolmogorov differential equations in special cases, and secondly, a set of criteria for the optimum allocation of genotypes in the host population of a host-pathogen system. The optimum allocation of genotypes in the host population is a problem of practical importance in controlling plant pathogens. In part two the theory of restricted transition probabilities developed in the companion paper is applied in finding the distribution of the time to the appearance of the first mutation for the case of a two dimensional birth process. The distribution of the time to the appearance of the first mutation is of importance in understanding the role mutation plays in the evolution of a population, particularly in the pathogen population of a host-pathogen system.

Notes: Abstract The complex arrangement of the muscle fibers in the ventricular wall and the nonsymmetric contraction and expansion of the ventricle preclude the writing of a differential equation of motion for the ventricle as a whole. We can, however, describe the motion of the ventricle by describing the motion of the dimensional parameters length and diameter; the radius, circumference, cross-sectional area, and volume following naturally from these. The ventricle is assumed to be an ellipsoid of revolution and the dimensional parameters to be periodic functions of time. Each of the parameters is expressed as a Fourier series.

Notes: Abstract Le modèle de Nelson peut-être considéré comme une approximation du modèle de Hodgkin-Huxley. Moins précis, il est plus maniable. Le modèle de Nelson peut également être considéré comme une généralisation du modèle de Hodgkin-Huxley. En effet, il introduit des liaisons synaptiques localisées ou diffusantes, et un processus de facilitation. Le mécanisme des liaisons synaptiques ne se traduit pas facilement dans le langage de Hodgkin-Huxley. Par contre, le processus de facilitation s'interprète facilement. Nelson's model can be taken as an approximation of Hodgkin-Huxley's model. Its precision is lesser, but it is more usable. Nelson's model can also be taken as a generalization of Hodgkin-Huxley's one; for it introduces localized or diffusing synaptic connexions and a facilitating process. The mechanism of synaptic connexions cannot be easily translated into Hodgkin-Huxley's language. On the contrary, the facilitating process is easily interpreted.

Notes: Abstract A spherical model for the human left ventricle with two different types of aneurysms, circular and rectangular-square, is proposed and meaningful relations are derived between the parameters of the aneurysms and ventricle. Such ventricular parameters as stroke volume, end-diastolic volume, and end-systolic volume are given normal human values to compute values for end-systolic radius and percentage shortening of muscle for various sized circular and rectangular-square aneurysms.

Notes: Abstract The Volterra theory of two competing populations is extended to the contemporary social problem of crime control. Domains of stability for the time dependence of the numbers in the criminal and enforcement groups are exposed by a numerical example. Both augmentation and reduction of enforcement can produce a stable system. Average values of the ratio of members in each group show great sensitivity to the control policies adopted by the remaining sector of the total population.

Notes: Abstract The paper deals with interactions of substances via an enzymatic reaction (Bull. Math. Biophysics,25, 141–154, 1963). The substances are the activators, inhibitors and/or substrates of the reaction. Due to the bimolecularity of the processes in the reaction, the quantitative relation between the steady state amount of complexes and the amounts of the substances assumes a typical form. In multiple enzymatic reactions this form is more complicated, though basically similar. Because the substances may influence the steady state amounts of the complexes in opposite directions, the compensation and blocking effects are the properties of enzymatic reactions. The substances with the same direction of influence may potentiate each other. In the enzymatic reaction here considered, the potentiation is always non-negative.

Notes: Abstract Growth-rate functions in analytic form have been obtained for cell cultures in which the doubling times follow the Gaussian and Poisson distributions. The growth-rate functions are calculated by using Laplace transforms to solve an integral equation previously presented. Oscillatory solutions result if a substantial fraction of the cells in a culture are synchronized to divide at some particular time. The synchrony and, hence, the oscillatory character of the growth-rate function eventually disappear because of the non-zero variance of the doubling-time distribution. If their variances are sufficiently small, the Gaussian and Poisson doubling-time distributions lead to growth-rate functions that become identical in the limit of large time.

Notes: Abstract IfN(t) is the expected number of cells in a culture at timet, $$\dot N(t)$$ the corresponding time derivative, andf(t−τ)dt the probability that a cell of aget−τ at timet will divide in the succeeding time intervaldt, then according to Hirsch and Engelberg (this issue) there obtains the integral equation $$\dot N(t) = 2\int_{ - \infty }^t {f(t - \tau )\dot N(\tau )d\tau }$$ for describing the dynamics of the cell population. It is the purpose of this note to give two alternative derivations of this equation, one based on the age density equation of Von Foerster, and the other based on a generalized form of the Harris-Bellman equation describing the first moment of an age dependent, branching process. In addition, a probability model is posed from which the Von Foerster equation and, hence, the Hirsch-Engelberg equation readily follows.

Notes: Abstract A model of the dissolution kinetics of powdered enamel is developed based on the kinetic rate termq, $$q = K'H - k'C \cdot P_1$$ , whereH=[H +],C=[Ca ++] andP 1=[HPO 4 = ]. The differential equations describing the rate of mineral dissolution (and the linearly related rate of appearance of calcium and phosphate in solution) have been derived and solved for three basic cases: (1) when thepH of the solution and surface area of the enamel are considered constant, (2) when thepH is assumed constant, but the reduction in surface area during dissolution is considered, and (3) when the rise ofpH resulting from the buffering effect of the dissolved enamel is considered along with the change in surface area. Analytical solutions have been obtained for cases (1) and (2), while a numerical solution has been found for case (3). Data from a study on enamel dissolution are presented that agree with the theory of case (3), and it is noted that apH rise as large as 0.5 can occur, as has been shown elsewhere in the literature.

Notes: Abstract On the basis of Landahl's theory of two-choice learning it is shown that application of punishment for wrong responses, without giving award for correct ones, does not lead to complete learning, no matter how many trials are used. If initially a “wrong response” was learned, then an attempt to inhibit it by punishment alone will in a class of cases lead only to a 50% suppression of that wrong response. Possible connection with the problem of effectiveness of punishment as a deterrent for crime is mentioned.

Notes: Abstract A set of characteristic parameters is given for electrophoresis accompanied by diffusion, followed by a method of simplification of the transport equation. The concept of electrophoretic similarity is introduced in connection with the presentation of solutions and the final section contains some dimensional considerations of the potential equation.

Notes: Abstract We show that when we represent (ℓ, ℛ)-systems with fixed genome as automata (sequential machines), we get automata with output-dependent states. This yields a short proof that ((ℓ, ℛ)-systems from a subcategory of automata—and with more homomorphisms than previously exhibited. We show how ((ℓ, ℛ)-systems with variable genetic structure may be represented as automata and use this embedding to set up a larger subcategory of the category of automata. An analogy with dynamical systems is briefly discussed. This paper presents a formal exploration and extension of some of the ideas presented by Rosen (Bull. Math. Biophyss,26, 103–111, 1964;28, 141–148;28 149–151). We refer the reader to these papers, and references cited therein, for a discussion of the relevance of this material to relational biology.

Notes: Abstract A two-dimensional nonlinear integro-differential equation with time-varying coefficients describing the behavior of the fluttering wing-body systems typical of natural flight mechanisms has been deduced from the Navier-Stokes equation which generalizes local pressure and velocity distributions in the externally oscillating air field. The resulting equation for the wing forces is combined with an analogous expression for the forces of gravitation and acceleration associated with the body. The air acceleration force, not previously considered in bio-physical models of insect and bird flight, is shown to arise from a formal analysis of unsteady or time-varying contributions to the velocity field, while the square form of the conventional steady state aerodynamic forces is derived from the intertial terms in the Navier-Stokes equation with the aid of the approximations of Newtonian impact theory. Previous calculations (Houghton, 1964) have indicated that the contribution to gravitational stability of air acceleration and aerodynamic life are roughly in the ratio of 3:1.

Notes: Abstract Certain types of cortical electrical events are non-propagated so that the associated electric fields must have standing wave characteristics. However, cortical electric events typically are generated by neurone populations which cannot be activated simultaneously on impulse driving. Hence the sum of the standing wave fields due to asynchronous activation of adjoining regions of cortical neurones must give the appearance of a traveling wave. Analysis of cortical waveforms is further complicated by curvature in cortical surfaces. A model is presented that shows the effects of curvature and time lag in activation on the form of the potential at points in space around a laminar array of elements simulating a population of cortical neurones. The results are compared with waveforms evoked by single-shock stimulation of the prepyriform cortex in cats.

Notes: Abstract A continuity equation for cell-number density in a population of cells is derived, and a system of equations for eliminating parameters between the general solution and the initial distribution obtained.

Notes: Abstract Use of an electrical model of the left ventricle of the heart and the arterial system permits analysis of the changes which take place as the capacity of the myocardium for generation of force decreases. The model is simple in structure, and its construction and practical testing would not be difficult. It demonstrates that, as the heart muscle weakens, the peak of intracardiac force occurs later in systole, and the difference between the intracardiac pressure and the aortic pressure in the second half of systole is much greater than for the normal heart. The feedback mechanisms which are proposed to affect myocardial contractility would affect this compensation for cardiac weakening. Indices to categorize the behavior of the normal, compensated though weakened, and decompensated myocardium are proposed.

Notes: Abstract The frequency distribution in a population of cells of the quantityCD (defined as the amount of some chromosomal substance in a cell which containsC chromosomes) is calculated using the distribution in the population of the amount per chromosome,D, and the distribution of chromosome number,C.

Notes: Abstract The rate of chromosomal DNA synthesis in an exponentially growing population of cells having chromosome-number dispersion is calculated using DNA histogram data, chromosome-number distribution data, and the assumptions that the synthesis rate is constant and DNA double exactly.

Notes: Abstract An estimate is made of the rate coefficient for linear DNA synthesis with exact doubling in an exponentially growing population of Ehrlich ascites tumor cells having chromosome-number dispersion. Comparison of calculated and experimental results suggest that the assumptions used in the calculation are tenable, but further experimental evidence is needed to prove this.

Notes: Abstract The paper develops further some suggestions made previously (Bulletin of Mathematical Biophysics,28, 283–308, 1966) that certain biological phenomena may be more easily interpreted from a “sociological” point of view by considering the organism as a social aggregate of cells and a cell as a social aggregate of genes. In this light the problems of origin of life on earth, of aging, and of parasitism and symbiosis are discussed. The notion of social aggregates of different orders is introduced.

Notes: Abstract A theoretical and experimental study was made of the mechanical behavior of the cornea. The theoretical analysis included an analytical solution for the symmetrical constraint of a thin, shallow, spherical shell by a rigid indenter. The experimental study investigated the rheology of the cornea with particular emphasis on its compliance with the requirements of the Boltzmann Superposition Principle. Representative results of tests on twenty enucleated hog eyes and two human eyes have been reported. The corneas of the human and hog eyes behaved as linear viscoelastic solids; the human eyes differed from the hog eyes in having a long term creep component. Several eyes were tested at the site of procurement, six to seven minutes after the animal's death, and it was established that creep is not an artifact due to aging or enucleation. The analytical and experimental results were combined to study some instruments used to detect the level of pressure in the eye. The theoretical analysis predicted that a type of elastic instability occurs during the process of flattening a small portion of the cornea; this is discussed with reference to the Goldmann and Mackay-Marg tonometers. The role of corneal creep was considered with reference to the response of the Schiøtz indentation tonometer during the time dependent process known as tonography.

Notes: Abstract Following previous studies, differential equations are established which determine the variation of the stimulus towards a corrective turn of the steering wheel and its effect on the excitation of the centers in the brain which results in the production of the corrective turn. The equations are derived under the highly oversimplified assumption that all excitation thresholds are so small that they can be neglected. Under these assumptions it is found that the tracking curve of a car is a sinusoid with negative damping, that is, with an ever increasing amplitude. Driving under these assumptions is imposible since the car will always eventually jump off the road. The possible effects of the threshold as well as stimuli towards corrective turns other than the distance from the edge of the lane are very briefly discussed. In spite of the negative results of the paper, its interest lies in the circumstance that with the complication of the model, we find that driving depends not only on the reaction times as the only “purely biological” parameter, but on three other neurobiophysical constants. In a subsequent paper (Rashevsky, 1967) it is shown how the introduction of one or more purely biological parameters of the driver makes a stable driving regime possible.

Notes: Abstract A vast number of biologically important processes are based upon bimolecular systems. In these systems intermediate complexes are formed. Bimolecular systems in which no complex-complex interactions occur are called linear systems of complexes. A definition and some characteristic properties of these systems are given here. There may exist a contradiction of Onsager's principle of detailed balancing in these systems; however, no principal differences are found between the steady state behavior of an open system and that of a closed system. It is shown that the steady state behavior of a linear system of complexes of arbitrary complexity has some similarities with the steady state behavior of a simple bimolecular system, e.g., Michaelis-Menten enzymatic reaction. Multiplicity of action of the substances participating in biomolecular processes may produce some qualitative differences in the steady state behavior of the system.

Notes: Abstract A time-dependent DNA histogram is calculated for an irradiated population of cells under the limiting assumption that the cells cannot pass through prophase due to the effects of the radiation. The population is assumed to increase exponentially prior to irradiation, but after irradiation to neither gain nor lose cells. Chromosome-number dispersion is taken into account in the calculation. The qualitative behavior of the calculated and experimental histograms are in reasonable agreement. The quantitative agreement between the two is relatively good at short post-irradiation times but is poor at long post-irradiation times (say, greater than half the doubling time). This suggests that recovery phenomena cannot be neglected at long post-irradiation times.

Notes: Abstract It is pointed out that the three different stimuli for a corrective turn, namely the distance from the edge of the lane, the rate of approach to the edge, and the angle between the direction of the car and the direction of the lane (Bull. Math. Biophysics,28, 645–654, 1966,29, 181–186, 1967) may act all three simultaneously. It is found that in that case the tracking curve of the car is stable below a critical speed and becomes unstable above it.

Notes: Abstract Continuing a previous study (Bull. Math. Biophysics, 28, 645–654, 1966), the biophysical mechanism of a corrective turn is investigated for the case where the stimulus for the corrective turn is produced not only by the perception of the nearness of an edge of the lane, but also by the rate of approach of the car towards the edge. In that case it is found that the tracking curve of the car may consist of a series of damped sinusoids and safe driving would be possible at any speed if it were not for the endogenous fluctuation in the driver's central nervous system. If the effect of the rate of approach increases sufficiently rapidly as the distance to the edge of the lane decreases, then a stable undamped oscillating tracking curve is possible. The case is also studied where the driver makes a corrective turn in response to a direct perception of the angle between the direction of the lane and the longitudinal axis of the car.

Notes: Abstract The principle of minimal work requires that the conducting airways of the human lung should have a maximum radius for minimal resistance to gas flow. At the same time there is a requirement that the airways should have a minimal volume for economy of space. These two opposing requirements have been investigated mathematically, and a method for calculating the angle of branching which produces minimal volume has been derived. The relationship of the radii of the parent and daughter branches to produce minimal resistance has been similarly defined. By measurement of a bronchial cast from a human lung the extent to which the predicted optimum structure is realized in practice has been shown. The change in structure associated with change of function at the transition from conducting airway to diffusion zone has been demonstrated.

Notes: Abstract This paper considers a class of set-theoretical entities, calledn-rank Linnaean structures, which are intended as abstract models of the taxonomic classificatory systems of biology. In the first part, devoted to formalism, finite Linnaean structures are discussed in complete generality; but, in addition, eight distinct subclasses are noted and some of the properties of their elements are explored. In the second part, concerned with applications, it is shown that taxonomic systems may be recast in the form of finite Linnaean structures, and an effort is made to show that some undesirable features of earlier models are avoided without artificiality and without abandoning extensional mathematics.

Notes: Abstract Generalizing an idea of M. Richardson (Fundamentals of Mathematics, New York: Macmillan Co., 1958), an APS on a given populationP is a non-empty collection of non-empty subsets ofP such that ifA is in the collection andA⊆B, thenB is in the collection. From a structure of this kind a partial ordering ofP, called therelated bumping order, is derived. The question is raised as to what kinds of partial orderings can be so obtained. For structures determined by voting weights of the members of the population, a complete characterization of all possible bumping orders is obtained.

Notes: Abstract The “second method” of Liapunov is used to perform a stability analysis of a mathematical model of the neuron. This analysis is based on the hypothesis that the firing of the neuron coincides with a temporary state of instability of the system, and that the initiation of all-or-none process depends on the magnitude of membrane depolarization and its first time derivative. It is found that the stability (and hence the possibility of a second firing) is restored approximately when the rate of membrane repolarization is at a maximum. This result predicts that the duration of the period of absolute refractoriness in neurons would be about 75 per cent of the spike duration, and thus shorter than the value usually obtained from experimental measurements.

Notes: Abstract Some aspects of masking phenomena are considered in terms of the simplest possible model of two-factor neural elements. The effect of a number of variables can be accounted for, but the introduction of an internuncial element results in a masking function which need not be symmetric about zero delay interval. As an illustration, the results for a special case are compared with available data. In general, such a model results in a masking function which depends on the intensity, area, and duration of the stimuli, as well as on the temporal and spatial separation between them.

Notes: Abstract The general equations are discussed describing two species in competition or in symbiosis or feeding one on the other.

Notes: Abstract It is shown that the principle of biological epimorphism (Rashevsky,Mathematical Principles in Biology and Their Applications, Springfield, Ill.: Charles Thomas, 1960) is contained in the theory of organismic sets (Bull. Math. Biophysics,29, 139–152, 1967) if an additional postulate not directly connected to mappings is made.

Notes: Abstract This paper deals with bimolecular systems in which also complex-complex interactions occur. Because of the complexity of the problem, an approximation in a form of coupled linear systems of complexes (Bull. Math. Biophysics,29, 1–16, 1967) is considered. Two types of couplings, serial and parallel, are studied. In the serial coupling the nonlinear system of complexes has the same behavior as its subsystems. An entity, initial sensitivity, has interesting properties: in serial coupling it is at most equal to the product and in parallel coupling, at most equal to the sum of partial initial sensitivities.

Notes: Abstract Amplification effect in the catalytic bimolecular systems is a consequence of the kinetic characteristic of the catalyst. Two types of the coefficient of amplification are defined. The applicability of these definitions is given by the type of the bimolecular system. In a simple example it is shown that the concept of amplification is meaningful in these systems. Furthermore, two rules, analogous to those for a coupling of amplifiers, are derived for the two basic modes of coupling of catalytic systems. Thus, in biological systems the catalytic reactions may be regarded as biologically effective amplifiers.

Notes: Abstract It is postulated that cell hydration is governed by adsorption of water on cell proteins in accord with the Bradley adsorption isotherm, and that the action of a solute in the surrounding solution is to lower the vapor pressure of the solution so that cell water adsorption is decreased by moving down the Bradley isotherm. From these concepts, it is derived that cell volume (V) should be related to solute concentration (x) by the equationV=−E log10 x+F whereE andF are constants which are independent of type of solute. For a non-adsorbed solute this agrees well with experimental data. For solutes which are adsorbed by cell proteins, a correction in the above equation may be necessary at higher solute concentrations, which is shown to be compatible with various experimental data. The types of experiments which are generally used to support the osmotic pressure theory of cell hydration agree equally well with the adsorption theory. The virtue of the adsorption theory is that, unlike the osmotic pressure theory of cell swelling, it is compatible with permeability of the cell membrane to solutes, which has been experimentally observed for various solutes.

Notes: Abstract Conditions under which a time varying electromagnetic field problem (such as arises in electrophysiology, electrocardiography, etc.) can be reduced to the conventional quasistatic problem are summarized. These conditions are discussed for typical physiological parameters.

Notes: Abstract A compartmental lung model with any number of synchronously filling and emptying functional chambers and a common dead space or conducting region is considered. It is shown that the model gives rise to an output, in an open circuit washout determination, which is a weighted sum of exponentials. From estimates of these weights and exponential components, estimates of the model parameters can be recovered. Relations giving the unique correspondence between the output parameters and the model parameters are derived and the existence and uniqueness of solutions established.

Notes: Abstract A physical model that incorporates all the experimental information on the formation of the visual pigment rhodopsin is presented. The visual pigments consist of a chromophore bound to an appropriate protein. Thus rhodopsin (λm 505 mμ) is formed by a Schiff’s base linkage C19H27CH=NH+-opsin (λm 440 mμ) between 11-cis retinal (λm 380 mμ) and the protein opsin (λm 280 mμ). It is found that there exists a red shift in the spectrum of rhodopsin from the Schiff’s base. The model brings an explanation for this red shift. It is shown that such a shift may be due to a charge transfer process (R. S. Mulliken,J. Am. Chem. Soc.,74, 811–824, 1952) between an electron at the double bond of carbons C11−C12 and an atomic orbital of the sulphur present in cysteine. This provides an explanation of the presence of SH-groups in the protein after the absorption of light. A one-electron approximation is used and the dipole momentμ NV ; hence, the oscillator strengthf of the transitionNV is estimated and compared with the experimentally determined extinction coefficient ∈m by mixing 3.5×10−3 M of 11-cis retinal with 8.3×10−5 M of cysteine at pH ranges 6 through 8. Reasonable agreement is found. Solvent, concentration and temperature dependence are shown also.

Notes: Abstract By assigning coordinates to the environmental function space comprising all physical and mental stimuli, mathematical interpretations can be based on such terms as adaptability, and reactivity which relate to individuals interacting with their environment within a society. These psychometric concepts are incorporated into a framework of functional analysis, which permits the optimization of social change by maximizing the satisfaction integral through the use of variational or dynamic programming methods in conjunction with some optimal social policy. The approach provides a mathematical connection between psychology and sociology, and further demonstrates that existing forms of government are simulated by differential equations belonging to the same general class. The synthesis of new classes of functional equations describing social progress is visualized as a legitimate objective for abstract mathematical sociology.

Notes: Abstract The theory of imitative behavior as developed hitherto by the author was based on the assumption that each individual has a natural preference for one of the two mutually exclusive behaviors. The endogenous fluctuations in the central nervous system then result in the individual’s exhibiting the two behaviors alternately with a relative frequency determined by the natural preference. Imitation shifts the natural preference towards one or the other of the two mutually exclusive behaviors. In the present approach it is suggested that the relative frequency of the two mutually exclusive behaviors exhibited alternately is determined by maximizing the “satisfaction function” of the individual, that is by hedonistic factors rather than by purely random fluctuations. Corresponding equations are developed. It is shown that in certain cases, even when the imitation effect is absent, a sort of “pseudoimitation” may occur. Another situation leads, in the case of two individuals only, to a complete “division of labor” between them, with respect to the two behaviors. Each one exhibits only one behavior. After that imitation is introduced explicitly by assuming that imitation by one individual or another increases the satisfaction function of the imitating individual. Results thus obtained show similarities to the results of the old theory.

Notes: Abstract More than 20 years after its proposal, Keller and Segel's model (1971,J. theor. Biol.,30, 235–248) remains by far the most popular model for chemical control of cell movement. However, before the Keller-Segel equations can be applied to a particular system, appropriate functional forms must be specified for the dependence on chemical concentration of the cell transport coefficients and the chemical degradation rate. In the vast majority of applications, these functional forms have been chosen using simple intuitive criteria. We focus on the particular case of eukaryotic cell movement, and derive an approximation to the detailed model of Sherrattet al. (1993,J. theor. Biol.,162, 23–40). The approximation consists of the Keller-Segel equations, with specific forms predicted for the cell transport coefficients and chemical degradation rate. Moreover, the parameter values in these functional forms can be directly measured experimentally. In the case of the much studied neutrophil-peptide system, we test our approximation using both the Boyden chamber and under-agarose assays. Finally, we show that for other cell-chemical interactions, a simple comparison of time scales provides a rapid check on the validity of our Keller-Segel approximation.

Notes: Abstract The formal structure of evolutionary theory is based upon the dynamics of alleles, individuals and populations. As such, the theory must assume the prior existence of these entities. This existence problem was recognized nearly a century ago, when DeVries (1904,Species and Varieties: Their Origin by Mutation) stated. “Natural selection may explain the survival of the fittest, but it cannot explain the arrival of the fittest.” At the heart of the existence problem is determining how biological organizations arise in ontogeny and in phylogeny. We develop a minimal theory of biological organization based on two abstractions from chemistry. The theory is formulated using λ-calculus, which provides a natural framework capturing (i) the constructive feature of chemistry, that the collision of molecules generates specific new molecules, and (ii) chemistry's diversity of equivalence classes, that many different reactants can yield the same stable product. We employ a well-stirred and constrained stochastic flow reactor to explore the generic behavior of large numbers of applicatively interacting λ-expressions. This constructive dynamical system generates fixed systems of transformation characterized by syntactical and functional invariances. Organizations are recognized and defined by these syntactical and functional regularities. Objects retained within an organization realize and algebraic structure and possess a grammar which is invariant under the interaction between objects. An organization is self-maintaining, and is characterized by (i) boundaries established by the invariances, (ii) strong self-repair capabilities responsible for a robustness to perturbation, and (iii) a center, defined as the smallest kinetically persistent and self-maintaining generator set of the algebra. Imposition of different boundary conditions on the stochastic flow reactor generates different levels of organization, and a diversity of organizations within each level. Level 0 is defined by selfcopying objects or simple ensembles of copying objects. Level 1 denotes a new object class, whose objects are self-maintaining organizations made of Level 0 objects, and Level 2 is defined by self-maintaining metaorganizations composed of Level 1 organizations. These results invite analogy to the history of life, that is, to the progression from self-replication to self-maintaining procaryotic organizations to ultimately yield self-maintaining eucaryotic organizations. In our system self-maintaining organizations arise as a generic consequence of two features of chemistry, without appeal to natural selection. We hold these findings as calling for increased attention to the structural basis of biological order.

Notes: Abstract A model is developed to describe neuronal elongation as a result of the polymerization of microtubules and elastic stretching of the neurites by force produced by the growth cone. The model for a single segment with a single growth cone revealed a constant elongation rate, while the concentration of tubulin in the soma rises, and the concentration of tubulin becomes constant in the growth cone. Extending the model to a neurite with a single branch point and two growth cones revealed the same results. When the assembly or the disassembly rate of microtubules is unequal in both growth cones, transient retraction of one of the terminal segments occurs, which results in complete retraction of the segment when the difference in (dis)assembly rate between the two growth cones is large enough. When the model is applied to large trees, a maximal sustainable number of terminal segments as a function of the production rate of tubulin appears. Mechanisms to stop outgrowth are discussed in relation to the establishment of synaptical contacts between cells.

Notes: Abstract We have developed a new model describing the relationship between plasma and red cell tracers flowing through the lung. The model is the result of an analysis of the transport of radiolabeled plasma albumin between two flowing phases and shows that differences between red cell and plasma tracer curves are related to microvascular hematocrit. The model was tested in an isolated, blood-perfused dog lung preparation in which we injected51Cr-labeled red cells and125I-labeled plasma albumin into the pulmonary artery. From the tracer concentration-time curves at the venous outflow, we calculatedh r, the ratio of microvascular hematocrit to large-vessel hematocrit. In 18 baseline experiments,h r=0.92±0.01 (mn±sem) at a blood flow rate of 10.7±0.3 ml s−1. We determined the effects of (a) glass bead embolization, (b) alloxan, and (c) lobe ligation onh r. Embolization attenuated the separation between plasma and red cells (increasedh r), probably as a consequence of passive vasodilation. Alloxan enhanced separation of plasma and red cells (decreasedh r), possibly as a result of arteriolar vasoconstriction. Ligation of a fraction of the perfused tissue at constant flow did not cause significant change inh r in the remaining perfused tissue. The model assumes that large-vessel transit times are uniform and that all dispersion occurs in the microvasculature. A theoretical analysis apportioning dispersion between large and small vessels disclosed that the error associated with these assumptions is likely to be less than 15% of the measuredh r. We conclude from this study that the microvascular hematocrit model describes experimental plasma and red cell curves. The results imply thath r can be readily deduced from tagged red cells and plasma and can be accounted for in calculating permeability-surface area in diffusing tracer experiments.

Notes: Abstract We present a mathematical model of the cytotoxic T lymphocyte response to the growth of an immunogenic tumor. The model exhibits a number of phenomena that are seenin vivo, including immunostimulation of tumor growth, “sneaking through” of the tumor, and formation of a tumor “dormant state”. The model is used to describe the kinetics of growth and regression of the B-lymphoma BCL1 in the spleen of mice. By comparing the model with experimental data, numerical estimates of parameters describing processes that cannot be measuredin vivo are derived. Local and global bifurcations are calculated for realistic values of the parameters. For a large set of parameters we predict that the course of tumor growth and its clinical manifestation have a recurrent profile with a 3- to 4-month cycle, similar to patterns seen in certain leukemias.

Notes: Abstract We present a new, practical algorithm to resolve the experimental data in restriction site analysis, which is a common technique for mapping DNA. Specifically, we assert that multiple digestions with a single restriction enzyme can provide sufficient information to identify the positions of the restriction sites with high probability. The motivation for the new approach comes from combinatorial results on the number of mutually homeometric sets in one dimension, where two sets ofn points are homeometric if the multiset ofn(n−1)/2 distances they determine are the same. Since experimental data contain errors, we propose algorithms for reconstructing sets from noisy interpoint distances, including the possibility of missing fragments. We analyse the performance of these algorithms under a reasonable probability distribution, establishing a relative error limit ofr=Θ(1/n 2) beyond which our technique becomes infeasible. Through simulations, we establish that our technique is robust enough to reconstruct data with relative errors of up to 7.0% in the measured fragment lengths for typical problems, which appears sufficient for certain biological applications.

Notes: Abstract A simple chemical model of the idiotypic network of immune systems, namely the AB model, has been developed by De Boeret al. The complexity of the system, such as the steady states, periodic oscillations and chaotic motions, has been examined by the authors mentioned above. In the present paper, the periodic motions and chaotic behaviours exhibited by the system are intuitively described. To clarify in which parameter domains concerned the system exhibits periodic oscillations and in which parameter domains the system demonstrates chaotic behaviours the Lyapounov exponent is explored. To characterize the strangeness of the attractors, the fractal dimension problem is worked out.

Notes: Abstract We consider a stochastic mechanism of the loss of resistance of cancer cells to cytotoxic agents, in terms of unstable gene amplification. Two models being different versions of a time-continuous branching random walk are presented. Both models assume strong dependence in replication and segregation of the extrachromosomal elements. The mathematical part of the paper includes the expression for the expected number of cells with a given number of gene copies in terms of modified Bessel functions. This adds to the collection of rare explicit solutions to branching process models. Original asymptotic expansions are also demonstrated. Fitting the model to experimental data yields estimates of the probabilities of gene amplification and deamplification. The thesis of the paper is that purely stochastic mechanisms may explain the dynamics of reversible drug resistance of cancer cells. Various stochastic approaches and their limitations are discussed.

Notes: Abstract Dextran has been the most commonly employed test molecule for probing the selectivity of glomerular filtration to macromolecules of varying size. The usual theories for hindered transport of solid spheres through pores have limited utility in interpreting clearance data for dextran or other linear polymers because such polymers in solution more closely resemble random, solvent-filled coils than solid spheres. To provide a model for glomerular filtration of random-coil macromolecules, the equilibrium partitioning of random coils between cylindrical pores and bulk solution was simulated using Monte Carlo calculations, and those results were combined with a hydrodynamic theory for restricted motion of solvent-filled polymer coils in pores. The rates of transport predicted for either neutral random coils or for solid spheres of the same Stokes-Einstein radius were significantly lower than observed transport rates of dextran through the glomerular capillary wall or across synthetic porous membranes. This facilitation of dextran transport was modeled by postulating weak, attractive interactions between dextran monomers and the pore wall. The random-coil model with attractive interactions, modeled using a short-range, square-well potential, was found to adequately represent dextran sieving data in normal rats. Various limitations of this approach are discussed.

Notes: Abstract Method-dependent mechanisms that may affect dynamic numerical solutions of a hyperbolic partial differential equation that models concentration profiles in renal tubules are described. Some numerical methods that have been applied to the equation are summarized, and ways by which the methods may misrepresent true solutions are analysed. Comparison of these methods demonstrates the need for thoughtful application of computational mathematics when simulating complicated time-dependent phenomena.

Notes: Abstract The regulation of the interactions between the actin binding proteins and the actin filaments are known to affect the cytoskeletal structure of F-actin. We develop a model depicting the formation of actin cytoskeleton, bundles and orthogonal networks, via activation or inactivation of different types of actin binding proteins. It is found that as the actin filament density increases in the cell, a spontaneous tendency to organize into bundles or networks occurs depending on the active actin binding protein concentration. Also, a minute change in the relative binding affinity of the actin binding proteins in the cell may lead to a major change in the actin cytoskeleton. Both the linear stability analysis and the numerical results indicate that the structures formed are highly sensitive to changes in the parameters, in particular to changes in the parameter ϕ, denoting the relative binding affinity and concentration of the actin binding proteins.

Notes: Abstract To investigate morphogenesis and in particular circularization mechanisms in young mycelia, we observe cultures of the zygomyceteMucor spinosus and develop discrete models of two-dimensional filamental branching growth. The models are based on the hypothesis that the fungus secretes a regulatory substance that diffuses into the surrounding medium and is detected by the growing hyphae. We also present a simple Markovian growth model without such a feedback, but yielding to analytical computations.

Notes: Abstract In vivo volume growth of two murine tumor cell lines was compared by mathematical modeling to their volume growth as multicellular spheroids. Fourteen deterministic mathematical models were studied. For one cell line, spheroid growth could be described by a model simpler than needed for description of growthin vivo. A model that explicitly included the stimulatory role for cell-cell interactions in regulation of growth was always superior to a model that did not include such a role. The von Bertalanffy model and the logistic model could not fit the data; this result contradicted some previous literature and was found to depend on the applied least squares fitting method. By the use of a particularly designed mathematical method, qualitative differences were discriminated from quantitative differences in growth dynamics of the same cells cultivated in two different three-dimensional systems.

Notes: Abstract This paper develops and applies a dynamic mathematical model for optimal scheduling of nitrogen fertilization and irrigation that minimizes nitrogen leaching subject to a target level of yield. The analysis assumes a single crop grown during a single growing season of a given length. It is shown that substitution of water for nitrogen along a given plant growth path decreases nitrogen leaching and, therefore, groundwater contamination. It is proved that a minimum leaching solution to the optimization problem is obtained with a single nitrogen application at the beginning of the season and irrigation scheduling that maintains a wet soil throughout the growing period. A numerical example utilizing experimental data for an irrigated summer corn in Israel confirms and quantifies the analytical findings.

Notes: Abstract We analyse the stochastic properties of dynamical systems with finite populations of a few differentreplicator species. Our main interest is to evaluate the typicallifetime, i.e. the time for the extinction of the first species in the network, for different catalytic structures, as a function of the population size.

Notes: Abstract The capacity of a model immune network in terms of the number of different antigens that can be vaccinated against without any memory lost is computed and tested by numerical simulations. We also investigate memory loss and failure to vaccinate due to overcrowding the network with too many antigens. The computations are done for two different strategies for proliferation, one implying all the antigen specific clones and the second one being more thrifty.

Notes: Abstract The technique of model-building a protein of known sequence but unknown tertiary structure from the structures of homologous proteins is probably so far the most reliable means of mapping from primary to tertiary structure. A key step towards the realization of the aim is to develop ways of aligning three-dimensional structures of homologus proteins, thereby deriving the rules useful for protein modelling. We have developed a generalized differential-geometric representation of protein local conformation for use in a protein comparison program which aligns protein sequences on the basis of their sequence and conformational knowledge. Because the differetial-geometric distance measure between local conformations is independent of the coordinate frame and remains chirality information, the comparison program is easily implemented, relatively rational and reasonably fast. The utility of this program for aligning closely and distantly related homologous proteins is demonstrated by multiple alignment of globins, serine proteinases and aspartic proteinase domains. Particularly, the method has reached the rational alignment between the mammalian and microbial serine proteinases as compared with many published alignment programs.

Notes: Abstract New formulas for deriving the sensitivities of stable stage structures and reproductive values to changes in vital rates are presented. They enable comparison of the sensities to changes of different elements in the projection matrix; in other words, comparison of partial derivatives of the eigenvectors. These kinds of sensitivities can be used in applied problems such as an analysis of the effect of harvesting on the population structure. However, in this paper, we examine the application of the sensitivities in a more general ecological context. We investigate why the stable stage structure of the mustard aphid,Lipaphis erysimi, changes very little in the temperature interval 10–30°C. The sensitivities of the stable stage structure at 15°C and 25°C were derived. The character of the sensitivites were the same in both temperatures although the stage structure was more sensitive to changes at 15°C than at 25°C. The sensitivity analysis also revealed that the temperature variation results in changes in fecundity and developmental rate that have a counteractive effect on the population structure.

Notes: Abstract Plankton populations undergo dramatic surges. Rapid increases and decreases by a factor of 10 or more are observed, often separated by relatively stable interludes. We propose a description of plankton communities as excitable systems. In particular, we present a model for the evolution of phytoplankton and zooplankton populations which resembles models for the behaviour of excitable media. The parameter dependency of the various “excitable” phenomena, trigger mechanism, threshold, and slow recovery, is clear, and permits ready investigation of the influence of properties of the physical environment, including variations in nutrient fluxes, temperature or pollution levels.

Notes: Abstract Analysis schemes for the classification of synergism and antagonism for mixed agents operate on the discrepancies between observed and calculated results. As such they cannot be confirmed by experiments and therefore have to be tested in terms of mathematical and logical self-consistency. The concept of independent action is close to the literal meaning of the term “non-interaction”. Since this concept does not depend on the mechanisms of actions nor on the type of effect scale used, it is suitable as one of the basic criterion for the definition of synergism and antagonism. A general mathematical framework of independent action is presented in this paper based on the concept of “relative effect” as used in the literature. The, different equations for independent action currently used in various areas are shown to be manifestations, of a general formula under different sets of boundary conditions, which are the natural limiting values of the effects of the corresponding system observed at low and at high doses of the agents. The framework can, be generalized to the combined action ofn-agents as well as to the interaction of an agent with itself. In addition, the differential form of the formula for independent action is derived. This framework of systematic definitions and derived equations enable a more in-depth study of the implications of the concept of independent action and its relation to other concepts of non-interaction.

Notes: Abstract We model how auto-reactiveB cells are kept under control by an idiotypic network. Autoimmunity occurs when the control is broken by an infection or not achieved through an abnormal ontogenetic evolution. We describe the idiotypic network, viz., the central immune system, by idiotype-anti-idiotype pairs which are coupled to a set of highly connected clones, which interact with each clone of the network. Some clones of the central immune system recognize self-antigen. We find a huge variety of fixed points which can be classified as tolerant, autoimmune, and neutral states according to the concentration of the auto-reactive antibody. Most significant are auto-reactive clones which are a member of an idiotype-anti-idiotype pair. In a healthy individual, an autoimmune disease is induced by an antigen infection which triggers a transition from a tolerant to an autoimmune state. Autoimmunity is induced more readily by an antigen coupling to theanti-idiotype than by one interacting with the auto-reactive clone itself. We indicate a possible therapy which simply reverses the processes that have lead to the autoimmune disease. In the early development of the central immune system its highly connected, core part serves to draw the more specific clones of idiotype-anti-idiotype pairs into the network. In order to avoid autoimmunity in ontogenetic evolution the anti-idiotype of an auto-reactive clone must be formed in advance by a sufficiently long period of time. Thus, a well ordered succession of the appearance of the more specific clones is required.

Notes: Abstract A method of dimensionless time-scaling based on extrinsic expectation of life at birth but intrinsic to a system generating a survival distribution is introduced. Such scaling allows the survival fraction function and its associated mortality function to serve as Green's functions for their generalized equivalents. i.e. a “population” function and a “death” function. The analytical mechanics of utilizing these concepts are formulated, applied to the classical Gompertz and Weibull survival models, and discussed with respect to biological relevance.

Notes: Abstract It is now widely accepted that localized high concentrations of Ca2+ (Ca2+ domains) play a major role in controlling the time course of neurotransmitter release. In the present work we calculate the magnitude and the time course of Ca2+ domains that evolve in the vicinity of a Ca2+ channel and an adjacent release site. In the calculations we consider a accurately dimensioned Ca2+ channel. Moreover, the Ca2+ current is continuously adjusted with regard to the accumulated intracellular Ca2+ and, in addition, endogenous buffers are considered. The calculations, carried out by the software FIDAP, based on finite element method, show that the Ca2+ concentrations achieved near the release sites are significantly lower than claimed by other investigators. Furthermore, we present arguments indicating that the Ca2+ domains, regardless of their magnitude, do not play a role in controlling the time course of release of neurotransmitter.

Notes: Abstract Mammalian white blood cells are known to bias the direction of their movement along concentration gradients of specific chemical stimuli, a phenomenon called chemotaxis. Chemotaxis of leukocyte cells is central to the acute inflammatory response in living organisms and other critical physiological functions. On a molecular level, these cells sense the stimuli termed chemotactic factor (CF) through specific cell surface receptors that bind CF molecules. This triggers a complex signal transduction process involving intracellular biochemical pathways and biophysical events, eventually leading to the observable chemotactic response. Several investigators have shown theoretically that statistical fluctuations in receptor binding lead to “noisy” intracellular signals, which may explain the observed imperfect chemotactic response to a CF gradient. The most recent dynamic model (Tranquillo and Lauffenburger,J. Math. Biol. 25, 229–262. 1987) couples a scheme for intracellular signal transduction and cell motility response with fluctuations in receptor binding. However, this model employs several assumptions regarding receptor dynamics that are now known to be oversimplifications. We extend the earlier model by accounting for several known and speculated chemotactic receptor dynamics, namely, transient G-protein signaling, cytoskeletal association, and receptor internalization and recycling, including statistical fluctuations in the numbers of receptors among the various states. Published studies are used to estimate associated constants and ensure the predicted receptor distribution is accurate. Model analysis indicates that directional persistence in uniform CF concentrations is enhanced by increasing rate constants for receptor cytoskeletal inactivation, ternary complex dissociation, and binary complex dissociation, and by decreasing rate constants for receptor internalization and recycling. For most rate constants, we have detected an optimal range that maximizes orientation bias in CF gradients. We have also examined different desensitization and receptor recycling mechanisms that yield experimentally documented orientation behavior. These yield novel insights into the relationship between receptor dynamics and leukocyte chemosensory movement behavior.

Notes: Abstract Given two independent sequences of letters, we seek the probability distribution of the length of the longest matching word. This word can be in different positions in the two sequences and we consider both perfect and nearly perfect matching. We derive bounds and approximations for the probability and compare them with other bounds and approximations. The results can be applied to DNA sequences in molecular biology and generalized matching between two independent random sequences.

Notes: Abstract A basic control engineer's adage-the poles of a feedback compensator become zeros of the closed-loop system-admits difficulties of interpretation even in the most simple of cases; that of single-input, single-output. An earlier investigation has provided an analysis of this adage in a module-theoretic context for systems in state space form while avoiding restrictive assumptions on system minimality or squareness. The main result is expressed concisely in terms of an exact sequence of modules which include Ω-zero modules corresponding to the feedback system and the plant. Extended zero modules of Ω-type incorporate both finite invariant zero structure, and generic zero information which occurs when a system fails to be right-invertible. In the case of compensation in the feedback path, this main exact sequence reduces to a mathematically clear expression of the aforementioned adage: the Ω-zero module of the feedback system is precisely the direct sum of the Ω-zero module of the plant and the system pole module of the feedback compensator. This paper extends the previous work in order to avoid assumptions on causality in the plant. Implicit dynamical systems are employed, in lieu of systems in state space form. Once again, it is not assumed that the system is one-to-one or onto; and so the concepts of generic zeros and their modules are brought into the arena of implicit systems. The implicit system itself is assumed in this work to be regular; however, decoupling zeros are permitted. Moreover, input-decoupling zeros and system pole feedback relationships are considered.

Notes: Abstract When the problem is considered of obtaining a periodic description in state-space form of a linear process which can be modelled by linear difference equations with periodic coefficients, it is natural to ask whether it is possible to preliminarily derive a polynomial equivalent form of such equations, which in the periodic case plays a role similar to the Rosenbrock's polynomial matrix description of a linear time-invariant process. In this paper a polynomial time-invariant description of a linear periodic process is introduced. It is shown that such a polynomial description gives a simple characterization of the dimension of the space of the solutions corresponding to the null input function, i.e., of the order of the periodic model under consideration. In addition, it allows us to introduce a transfer matrix for the computation of the output responses corresponding to null initial conditions, and to deduce conditions for the periodic model to be causal. These results, as well as the possibility of defining strict system equivalence between two periodic models through their time-invariant polynomial descriptions, in a similar sense as in the time-invariant case, show the relevance of such a polynomial time-invariant description for the problem under consideration.

Notes: Abstract An actual sampling process can be modeled as a random process, which consists of the regular (uniform) deterministic sampling process plus an error in the sampling times which constitutes a zero-mean noise (the jitter). In this paper we discuss the problem of estimating the jitter process. By assuming that the jitter process is an i.i.d. one, with standard deviation that is small compared to the regular sampling time, we show that the variance of the jitter process can be estimated from thenth order spectrum of the sampled data,n=2, 3, i.e., the jitter variance can be extracted from the 2nd-order spectrum or the 3rd-order spectrum (the bispectrum) of the sampled data, provided the continuous signal spectrum is known. However when the signal skewness exceeds a certain level, the potential performance of the bispectrum-based estimation is better than that of the spectrum-based estimation. Moreover, the former can also provide jitter variance estimates when the continuous signal spectrum is unknown while the latter cannot. This suggests that the bispectrum of the sampled data is potentially better for estimating any parameter of the sampling jitter process, once the signal skewness is sufficiently large.

Notes: Abstract Functions belonging to various Paley-Wiener spaces have representations in sampling series. When a function does not belong to such a space, the sampling series may converge, not to the object function but to an "alias" of it, and an aliasing error is said to occur. Aliasing error bounds are derived for one- and two-channel sampling series analogous to the Whittaker-Kotel’nikov-Shannon series, and for the multi-band sampling series, and a "derivative" extension of it, due to Dodson, Beaty, et al. The Poisson summation formula is a basic tool throughout. Aliasing in the one-channel case is shown to arise from a transformation with similarities to a projection. Where possible, the sharpness of the error bounds is discussed.

Notes: Abstract This paper is devoted to a study of the Hausdorff-Young theorems from a historical perspective, beginning with the F. Riesz-Fischer theorem. Introduced by W. H. Young (1912), these theorems were considered and extended by F. Hausdorff (1923), F. Riesz (1923), E.C. Titchmarsh (1924), G. H. Hardy and J.E. Littlewood (1926), M. Riesz (1927), and O. Thorin (1939/48). Special emphasis is placed upon the development of the proofs of the two Hausdorff-Young inequalities and their impact upon Fourier analysis as a whole, in particular on the M. Riesz-Thorin convexity theoremand on the interpolation of operators. The golden thread connecting the various extensions and generalizations is the concept of logarithmic convexity, one that goes back to the work of J. Hadamard (1896), A. Liapounoff (1901), J.L.W.V. Jensen (1906), and O. Blumenthal (1907).

Notes: Abstract In this paper we give a further investigation of the method introduced by the author in [1, Frequency-domain bounds for nonnegative unsharply band-limited functions] for proving bounds for functions with nonnegative Fourier transforms. We also dealt with the question of how large the supremum KS of all numbers |f(u)| is with f the Fourier transform of a nonnegative integrable function F and f(0) = 1, |f(ku)| ≤ ε for k ∈ S. Here u 〉 0 and S ⊂ {2, 3, . . .}. This problem was related in [1] to finding the infimum MS of all numbers Mh = maxϑ [(1−h(ϑ))/(1− cos ϑ)] over all 2π-periodic even, smooth functions h whose Fourier cosine coefficients ak vanish for k ∉ S, and it was proved and announced for several cases that MS (1−KS ) = 1. In this paper we prove the results announced in [1]. To that end we generalize the method given in [1] to include Fourier transforms f of probability measures on R and a certain generalized function h, and we show that the numbers KS, MS are assumed as |f(u)|, Mh for certain allowed f,h. Moreover, we establish a fundamental relation between finding the numbers KS, MS and the numbers KT, MT where T = {2, 3, . . .}\S. In particular, we show that MT = 2KS (2KS − 1)−1,KT = 1/2 MS(MS − 1)−1 and that MT (1 − KT) = 1,KSKT = 1/2 , whenever MS (1 − KS) = 1.

Notes: Abstract Finite energy band-limited functions are reconstructed iteratively from nonuniform sample values of the functions and its derivatives. It is shown that the maximum gap allowed between the sampling points increases linearly with the number of derivatives considered. Moreover, a more precise result is presented for the first derivative case and another reconstruction of the functions using the frame algorithm is deduced.

Notes: Abstract In the early 1960s research into radar signal synthesis produced important formulas describing the action of the two-dimensional Fourier transform on auto- and crossambiguity surfaces. When coupled with the Poisson Summation formula, these results become applicable to the theory of Weyl-Heisenberg systems, in the form of lattice sum formulas that relate the energy of the discrete crossambiguity function of two signals f and g over a lattice with the inner product of the discrete autoambiguity functions of f and g over a "complementary" lattice. These lattice sum formulas provide a framework for a new proof of a result of N.J. Munch characterizing tight frames and for establishing an important relationship between l1-summability (condition A) of the discrete ambiguity function of g over a lattice and properties of the Weyl-Heisenberg system of g over the complementary lattice. This condition leads to formulas for upper frame bounds that appear simpler than those previously published and provide guidance in choosing lattice parameters that yield the most snug frame at a stipulated density of basis functions.

Notes: Abstract Let $a〉0, b〉0, ab〈1;$ and let $g\in L^2({\Bbb R}).$ In this paper we investigate the relation between the frame operator $S:f\in L^2({\Bbb R})\rightarrow \sum_{n,m}\,(f,g_{na,mb})\,g_{na,mb}$ and the matrix $H$ whose entries $H_{k,l\,;\,k',l'}$ are given by $(g_{k'/b,l'/a},g_{k/b,l/a})$ for $k,l,k',l'\in{\Bbb Z}.$ Here $f_{x,y}(t)={\rm exp}(2\pi iyt)\,f(t-x),$ $t\in{\Bbb R}$ , for any $f\in L^2({\Bbb R}).$ We show that $S$ is bounded as a mapping of $L^2({\Bbb R})$ into $L^2({\Bbb R})$ if and only if $H$ is bounded as a mapping of $l^2({\Bbb Z}^2)$ into $l^2({\Bbb Z}^2).$ Also we show that $AI\leq S\leq BI$ if and only if $AI\leq\frac{1}{ab}\,H\leq BI,$ where $I$ denotes the identity operator of $L^2({\Bbb R})$ and $l^2({\Bbb Z}^2),$ respectively, and $A\geq 0,$ $B〈\infty.$ Next, when $g$ generates a frame, we have that $(g_{k/b,l/a})_{k,l}$ has an upper frame bound, and the minimal dual function $^{\circ}\gamma$ can be computed as $ab\,\sum_{k,l}\,(H^{-1})_{k,l\,;\,o,o}\,g_{k/b,l/a}.$ The results of this paper extend, generalize, and rigourize results of Wexler and Raz and of Qian, D. Chen, K. Chen, and Li on the computation of dual functions for finite, discrete-time Gabor expansions to the infinite, continuous-time case. Furthermore, we present a framework in which one can show that certain smoothness and decay properties of a $g$ generating a frame are inherited by $^{\circ}\gamma.$ In particular, we show that $^{\circ}\gamma\in{\cal S}$ when $g\in{\cal S}$ generates a frame $({\cal S}$ Schwartz space). The proofs of the main results of this paper rely heavily on a technique introduced by Tolimieri and Orr for relating frame bound questions on complementary lattices by means of the Poisson summation formula.

Notes: Abstract For any ε 〉 0, we construct an orthonormal Schauder basis of C(K) consisting of trigonometric polynomials Tn n = 1, 2, . . . , such that deg(Tn) ≤ (1/2)(1 + ε)n. This is best possible with regard to the degree. The construction uses wavelet techniques.

Notes: Abstract We study the general question of the existence of self-similar lattice tilings of Euclidean space. A necessary and sufficient geometric condition on the growth of the boundary of approximate tiles is reduced to a problem in Fourier analysis that is shown to have an elegant simple solution in dimension one. In dimension two we further prove the existence of connected self-similar lattice tilings for parabolic and elliptic dilations. These results apply to produce Haar wavelet bases and certain canonical number systems.

Notes: Abstract In the spirit of work of Kerman and Sawyer, a condition is given that is necessary and sufficient for the Fourier transform norm inequality $\Big(\int_{{\Bbb R}_d} \vert\hat{f}\vert^q d\mu\Big)^{1/q} \leq C\Big(\int_{{\Bbb R}_d} \vert f\vert^p v\Big)^{1/p}$ provided v is a radial weight for which v−1/p is convexly decreasing and μ is a suitable measure. We also establish alternative conditions for such inequalities by proving corresponding trace type inequalities and maximal function inequalities that underlie the Fourier transform estimates. Our conditions are relatively simple to compute. Among applications we give extensions of a Sobolev restriction theorem.

Notes: Abstract We present two-sided singular value estimates for a class of convolution-product operators related to time-frequency localization.

Notes: Abstract A linear-quadratic (LQ) control problem subject to a standard continuous-time system is called regular if the input weighting matrix is invertible, and singular if this is not the case. Consequently, optimal inputs for regular LQ problems are ordinary functions (state feedbacks), whereas optical controls for singular problems are in general distributions, e.g., impulses. We will show that regularity and singularity in LQ problems subject to ageneral (implicit) system depends not so much on the input weighting matrix, as on the property that the integrand of the cost criterion is a function only if inputs and state trajectories are, as is the case for LQ problems, subject to astandard system. In particular, we will provide a simple criterion for distinguishing between regularity and singularity in LQ problems subject to a general system. Our criterion is expressed in the system coefficients only and reduces to the classical one if the underlying system is standard.

Notes: Abstract Completions of linear time varying singular systems of the formE(t)x′(t)+F(t)x′(t)=f(t) are explicitly computed using recent results on rational matrix functions. The algorithm and the theory behind it are carefully described. Computational issues are discussed.

Notes: Abstract In this paper a spectral method using orthogonal periodic basis functions for the analysis of linear time invariant descriptor systems is discussed, and the case of the trigonometric Fourier functions is investigated in detail. The method is shown to be convergent, in the distributional sense. However, for any finite number of basis functions, the periodicity induced by the chosen basis can give rise to spurious impulsive components in the computed system response, even in the case of correct initial conditions.

Notes: Abstract In this paper we will study topological properties of the class of proper and improperp×m transfer functions of a fixed McMillan degreen. A natural generalization of this class is all autoregressive systems of degreen under external system equivalence. The subset of irreducible systems has in a natural way the structure of a manifold and we show how to extend this topology to the set of all autoregressive systems of degree at mostn. We will describe the subset of systems with fixed Kronecker indicesv=(v1,...,v p ) as an orbit space, which will enable us to calculate the topological dimension for each collection of indicesv. Finally, we will describe the topological closure of those sets in the space of all autoregressive systems.

Notes: Abstract It has been shown in [B.M.90] that non-square implicit differential equations allow for the description of variable structure systems (variable order, variable sign, variable parameters). We combine here the possible control strategy developed in [L.91] for rectangular systems (insuring a unique output behavior for the system compensated with a proportional or proportional and derivative state feedback) with the detector and the observer introduced in [B.M.90] in order to obtain a closed-loop system where the initial structure variation disappears on the output. We also give necessary and sufficient conditions for the free assignment of the associated output dynamics.

Notes: Abstract An exposition of joint cumulants and cumulant spectra is presented. A distinction is emphasized in this paper between the cumulant spectrum of a time series and its stationary version, here called apolyspectrum. The variance and covariance of the sample bispectrum is then derived using a relationship between cumulant spectra of the finite Fourier transform for the 2nd and 4th cumulant function, and the bispectrum and trispectrum of the time series.

Notes: Abstract The detection of a general class of transient (i.e., finite energy) signals in additive stationary interference using the spectral correlation function (second order cumulant spectrum) is presented. Observable features in the two-dimensional spectral correlation function due to properties of signals in the assumed class of transients are exploited to derive a detection statistic. The performance of the proposed detection statistic relative to a conventional power spectral detector is presented.