ALBERT

Notes: Abstract The partial differential equation of the random walk problem with persistence of direction and external bias is derived. By persistence of direction or internal bias we mean that the probability a particle will travel in a given direction need not be the same for all directions, but depends solely upon the particle's previous direction of motion. The external bias arises from an anisotropy of the medium or an external force on the particle. The problem is treated by considering that the net displacement of a particle arises from two factors, namely, that neither the probability of the particle traveling in any direction after turning nor the distance the particle travels in a given direction need be the same for all directions. A modified Fokker-Planck equation is first obtained using the assumptions that the particles have a distribution of travel times and speeds and that the average time of travel between turns need not be zero. The fional equation incopporating the assumption of a persistence of direction and an external bias is then derived. Applications to the study of diffusion and to long-chain polymers are then made.

Notes: Abstract The transmission of some information or behavior pattern is treated as a flow of “particles” which execute random motions over a population of individuals and which may multiply or disappear. Equations are derived for the number density of these “particles” and from this is calculated the number of individuals through which the “particles” have passed. The results are applied to a number of situations such as 1) uniform spatial distribution with multiplication factor decreasing with time because of loss of interest or confusion of the information, 2) multiplication factor constant but the rate of spreal decreasing with multiple hearings, 3) one-dimensional region with a small starting region with or without an absorbing barrier 4) two-dimensional region with absorbing barrier, 5) continous sources of information within a small region in one dimension, 6) uniform spatial distribution in which individuals do not respond to more than one hearing.

Notes: Abstract A situation is considered in which a fluid containing a substance flows through a vessel at a constant rate, the substance being permeable to the vessel wall. In the region outside the vessel there is supposed to be rapid mixing in the direction perpendicular to the axis of the vessel but no mixing longitudinally. The solution for the spatial distribution at any time is given for the case of an arbitrary initial distribution along the vessel length in the absence of an input. The solution is also given for the case of a single impulsive input, the concentration being initially zero everywhere.

Notes: Abstract Organisms orient themselves to a stimulus by two general methods. One method is by directed orientation (taxis); the other is by undirected locomotory reaction (kinesis). An equation, and the methods for finding the necessary parameters of this equation, is derived for the distribution of organisms within a container, with the following limitations: (1) the organisms have no accommodation, (2) they are always active, and (3) the stimulus changes slowly with position. Necessary modifications of the equation are then derived, so that the last two limitations may be eliminated. The equation cannot be solved excatly because of its complexity; hence an approximation method must be used. This method is discussed, an approximate solution is found, and a time constant for equilibrium to be established is derived. Applications tovarious experiments in the literature are then made with fairly satisfactory results. A new interpretation of the theory of klino-kinesis with accommodation is found upon application of the equations developed to experimental work. Further limitations and uses of these equations are then discussed.

Notes: Abstract Certain parameters are defined which roughly characterize the internal structure of networks. A given network structure uniquely determines the values of the parameters, but the reverse is not true. The parameters therefore define certain classes of networks. One of the parameters, thedispersion D(S) gives an indication of the “compactness” of the internal structure. Addition theorems and inequalities are derived relating the dispersions of sub-systems to the dispersion of the complete structure.

Notes: Abstract A mathematical theory is developed which permits the determination of certain parameters of an inhomogenous tissue, such as a nerve trunk without its epineurium. The parameters are the permeability coefficients for entrance into an exit of a substance from the nerve fibers, and the diffusion coefficient of the interstitial material. The experimental data required are the dimensions of the cross-section, the average diameter of the fibers, and the ratio of the cross-sectional are of the fibers to the total cross-section, as well as the time course of the decrease of the fraction of the substance left in the nerve trunk, when the trunk is immersed in a bathing solution containing none of it.

Notes: Abstract A model based on enzyme localization is developed which gives rise to an apparent active transport of a metabolite into or out of cells. The model is applied to three simple situations, using Fick's equation and the Rashevsky approximation. It is shown that the apparent efficiency can be made as large as desired if, for constant reaction, the outer cell region is made sufficiently small, or, for autocatalytic reaction, if the metabolite concentration in the outer region is sufficiently small. The physical limitations imposed by this mechanism are developed for all three situations.

Notes: Abstract A previously derived iteration formula for a random net was applied to some data on the spread of information through a population. It was found that if the axon density (the only free parameter in the formula) is determined by the first pair of experimental values, the predicted spread is much more rapid than the observed one. If the successive values of the “apparent axon density” are calculated from the successive experimental values, it is noticed that this quantity at first suffers a sharp drop from an initial high value to its lowest value and then gradually “recovers”. An attempt is made to account for this behavior of the apparent axon density in terms of the “assumption of transitivity”, based on a certain socio-structural bias, namely, that the likely contacts of two individuals who themselves have been in contact are expected to be strongly overlapping. The assumption of transitivity leads to a drop in the apparent axon density from an arbitrary initial value to the vicinity of unity (if the actual axon density is not too small). However, the “recovery” is not accounted for, and thus the predicted spread turns out to beslower than the observed.

Notes: Abstract The assumption of transitivity treated in part I is modified in various ways to describe an information-diffusion process, in which a certain amount of randomness of contact does occur. In one model a parameter is introduced which is indicative of a tendency to go beyond one's immediate vicinity to spread the information as the vicinity becomes saturated with knowers. In another model the randomness appears in the assumption that new knowers are uniformly distributed among the knowers. Two of the equations thus derived, each with two free parameters are in good agreement with experimental results.

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 The impressed field, “Window Field” (WF), due to a half-wave action potential on a muscle fiber, has been calculated on the basis of potential theory. It has been shown that in spite of the small intensity of the field, its integrated action can transfer the energy needed to induce, contraction from the membrane to the interior of the fiber. The energy of polarization has been found to be sufficient to exceed the energy of, thermal agitation on that length of fiber, which can be identified as the length of a sarcomere. The changes of ion concentration, caused by the WF, if calculated on the assumption of the semipermeability of theZ membranes, was found to be equal to the changes necessary to induce contraction of actomyosinin vitro.

Notes: Abstract Some general properties of the solution of the diffusion equation are deduced for the steady-state, spherically symmetric system. On the basis of these developments some results of N. Rashevsky (Bull. Math. Biophysics,11, 15, 1949) are discussed and the results of a previous investigation (Hearon,Bull. Math. Biophysics,12, 135, 1950b) are extended to more general conditions. In particular these extensions apply to the flow of a soluteagainst its concentration gradient, the nonzero gradient of an inert metabolite, and theaccumulation or exclusion of an inert metabolite in a metabolic system.

Notes: Abstract The approximation method of N. Rashevsky is discussed and reviewed. It is shown that in addition to theexplicit assumptions and approximations there is involved the assumption that the rate of metabolism is the same at every point in the cell and that theaverage rate of metabolism is different from zero. An expression is given for the error in the approximate method when the rate of metabolism is any function of the concentration. It is also shown that a solution in theform of that obtained by the approximate method is not possible if the generalized laws of diffusion are assumed to apply.

Notes: Abstract Rate equations for the enzymatic oxidation of succinic acid are derived on the assumption that when a single molecule of substrate combines with an enzyme molecule, it can do so with either one or two sites on the enzyme, and that oxidation occurs only in the second case. In addition it is assumed that the product of the reaction, fumaric acid, combines reversibly with the enzyme. With certain enzyme preparations the data fitted such an equation satisfactorily. In others the rate was that of a first-order reaction, but addition of cytochrome changed it to the former type. It was concluded that the transfer of hydrogen to oxygen was a first-order reaction and dominated the whole rate when enzyme preparations were used which had been washed relatively free of cytochrome. When the limiting factor was succino-dehydrogenase the rates followed the new equation. Criteria for recognizing noncompetitive inhibition are given, and inhibition by di-tertiary butyl peroxide was shown to be of this type.

Notes: Abstract This paper deals with the following question: Which distributions of radiosensitivity in a population can lead to an exponential survival curve? The problem is solved exactly, with statistical fluctuations in dose fully accounted for. It is shown that only an exponential distribution of sensitivities can give rise to an exponential survival curve.

Notes: Abstract An approximation method is introduced which enables a number of diffusion-type problems to be solved in an approximate but simple manner. Many cases require only the solution of a simple first-order differential equation. The method is applied to a number of cases in which the exact solutions are available. A comparison shows that the method is quite satisfactory in these cases. The method is applied to diffusion problems with rate of consumption proportional to concentration or to the square of the concentration. In the latter case, the result obtained is essentially the same as that found by H. G. Landau (1950) after elaborate calculations.

Notes: Abstract It is shown that a slight modification of a model of excitatory phenomena in irritable tissues, which has been treated before, exhibits spontaneous oscillations. The frequency of these oscillations and the time-course of the potential across the model membrane have been determined, together with the dependence of some of their characteristics on some important parameters, particularly (Ca++).

Notes: Abstract It is demonstrated that, if the variations of viscosity throughout a cell are considered, swelling stresses may produce elongation and division. To do this it is necessary to generalize Betti's theorem to cover systems containing viscosity gradients and such a generalization is presented. On the basis of two special assumptions it is shown that most of the results of the diffusion drag theory of cell division may be duplicated by the present theory.

Notes: Abstract The theory of imitative behavior, developed previously, is applied to the case of two social groups which are separated spatially. If the information of each group as to the behavior of the other is complete, the case reduces to that of a single group. When any information is lacking at all, the two groups are independent. If we have two mutually exclusive behaviorsA andB, all four combinationsAA, AB, BA, andBB are possible. If the mutual information gradually increases from zero, then for a certain value of it, the group which is more informed about the behavior of the other will change to that behavior if it did not already exhibit it. If for constant information the size of the group increases, then above a certain threshold value, the larger group imposes its behavior on the smaller.

Notes: Abstract Local regulation of blood flow as determined by capillary diameter and the number of open capillaries in a region is considered. The local changes in capillary diameter and in the number of open capillaries are assumed to be due to concentration changes of a diffusible, nonspecified metabolite. This metabolite is produced in the tissue and carried away by the blood stream. Using these assumptions and applying pertinent data on capillaries, deductions are made concerning: (a) the law of blood flow as a function of temperature and capillary radius for the hyperemia of high temperature, (b) high flow as it depends on metabolism during strenuous exercise of muscle, and (c) a first approximation to the time duration of occlusion hyperemia.

Notes: Abstract The necessary and sufficient condition is given forn integers to be the score structure of a society with a dominance relation. A proof is also given for a theorem showing that there are members who dominate every other member either directly or indirectly through a single intermediate member.

Notes: Abstract On the basis of a previous general formulation (Bull. Math. Biophysics,15, 21–29, 1953a) a discussion is given of the error in the approximation method of N. Rashevsky. This error, inherent in the method when the metabolic rate is different at each point in the cell, is discussed in detail and numerical values are presented for two particular cases: the rate proportional to the concentration and the rate a prescribed function of the spatial coordinates. It is shown that the formulation for the first case also applies to several other cases, that the error is negligible provided the rate is sufficiently small, and that the error is fairly sensitive to the cell size. If the rate depends upon the coordinatesalone a small rate is not sufficient to insure a negligible error. The relations between the exact method, the standard approximate method, an earlier approximate method (Physics,7 260, 1936), and a more recent refinement (Bull. Math. Biophysics,10, 201, 1948) of the standard method are discussed.

Notes: Abstract It is shown on the basis of (1) conservation of mass, (2) positive concentrations, and (3) the principle of detail balancing that periodic reactions cannot occur in a closed system described bylinear differential equations. The matrix,A, of the rate equations must be such that |A|=0,a ij〉0 fori≠j,a ii〈0, andVAV −1=B, whereV is diagonal andB is symmetric. These properties ofA imply that the latent roots are real and non-positive and that neither catalysis nor inhibition can be described bylinear equations. It is further shown that periodic reactions cannot occur in anopen system for which the matrix associated with the chemical reactions has the above properties and in which thesimple law of diffusion is obeyed. The relation of these results to Onsager's reciprocal relations and to previous work on periodic and cyclic chemical reactions is discussed. The utility of certain of these results for the treatment of isotope kinetics is indicated.

Notes: Abstract It is demonstrated that an explanation of the small radius effect or the so-called sigma phenomenon may be obtained by noting that one of the effects of the presence of suspended particles in a flowing fluid is to increase the velocity of flow near the wall over that existing in the absence of particles. This effect may be considered equivalent to relaxing the boundary conditions at the wall. An expression for the viscosity is compared with data and fit is found to be good.

Notes: Abstract The solution for the spatial distribution of ions in a Donnan equilibrium has been given by J. H. Bartlett and R. A. Kromhout (1952). The present note gives an explicit solution for the case in which the length of the region containing the membrane is large; in biological situations this requires only that the length considered should be greater than a few hundred Ångstrom units. The Donnan equilibrium may be considered to be a special case of a situation in which forces other than electrical act upon the ions; in particular, it represents the case in which only one ion is acted upon and the energy difference on the two sides of the membrane is infinite. An expression is given for the difference in energy of theith in terms of the electrical potential and of the ion concentrations. As an illustration, the results are applied to nerve membrane potentials.

Notes: Abstract A mechanism is described which accounts for the active transport of Na+ ions through a membrane. It is assumed that at one side of the membrane the ion combines with a carrier ion, the resulting carrier compound then diffuses through the membrane and decomposes at the other side of the membrane. The free diffusion of the ions is also taken into account. The time rate of accumulation of the ion in question at the latter side of the membrane is calculated in terms of the concentrations of the ion at both sides of the membrane.

Notes: Abstract The recent extension of the approximation method is applied to enable us to arrive at the time course of the concentrations at both sides of a membrane. From the differential equations which govern these, the steady-state solution is obtained in terms of the parameters, which are determined by the thickness of the diffusion layers, the chemical composition and reactions, and the diffusion constant of the membrane.

Notes: Abstract An equation is derived from the spread of a “state” by contact through a thoroughly mixed population, in which the probability of transmission depends both on the over-all duration of the process and on the time an individual has been in the “state.” Cases in which this probability is a function of only one or the other of the two “times” are worked out. It is shown that in the case of dependence on “private time” alone the asymptotic value of the fraction of the population effected is the same as that derived by the random net approach.

Notes: Abstract The reflection of pressure waves in a fluid enclosed within a tube with an elastic wall is studied for the case of a localized change in diameter of the tube. The concept of impedance is introduced. The relation of the reflection characteristics of the parts of the tube at either side of the change is derived on the basis of the continuity of pressure and mass flow at the site of the change. This relations is used to derive the expression for the ratio of the pressure oscillations measured in front of, and behind, the constriction in terms of the constants of the system. As a result, a method is indicated to locate the coarctation from measurements of the pressures in front of, and behind it.

Notes: Abstract The question is raised concerning the possible causes of abnormally small standard deviations found in counting samples in which particles are distributed at random (e.g., blood cells, fat globules in milk, etc.). The effect of discarding abnormal samples is discounted inasmuch as small standard deviations occur even when all samples are counted. An approximation method is used to calculate the effect of finite particle size, of known repulsive forces between particles and of convection currents. This calculation shows that neither finite size nor the known repulsive forces are sufficient to account for the observed abnormality of standard deviation, but that convection currents can possibly account for it. The possible presence of long-range repulsive forces cannot, however, be excluded.

Notes: Abstract The effect of finite particle size on the standard deviation in sample counts is computed for the one-dimensional case. To a first order of approximation the correction is found to be identifical with that found by H. de Vries (1953) using a general approximation method.

Notes: Abstract A critical examination of the “classical” theories of photoreception in view of more recent experimental findings yields the result that these theories do not possess the property to describe all the more significant phenomena of photoreception correctly, and to some extent suffer the lack of more general applicability. The basis for a new and presumably more general theory of photoreception based on dynamical aspects is laid out. Emphasis is put on the time course of afferent and efferent excitation in the photoreception model, consisting of a receptor element, an afferent and an efferent neuron of the one-factor Rashevsky-type, and an effector organ.

Notes: Abstract A mathematical model for the development of human society, beginning with the earliest stages of urban cultures, is outlined. In the early stages of history, behavior was characterized largely by adherence to a number of beliefs and prejudices of diffeirent kinds, which were accepted on faith and not subject to critical rational analysis. Due to psychobiological variability a very small number of individuals spontaneously appear at all times who challenge the accepted beliefs and prejudices and do not follow the accepted patterns of social behavior. The effect of these individuals upon the rest of the society, especially upon the younger generation, depends on the facilities with which information spreads in society. In earliest societies, when modern methods of mass communication were unknown, the channels of communication were practically identical with the channels of economic transport. The latter in its turn depended on the nature of the roads, and especially on the presence of waterway, which facilitated transportation. The sizes of the earliest cities and the distances between them were largely determined by relative ease of transporation. Expressions are derived for the average size of the earliest cities and for the average distance between them. The calculated average populations of the earliest cities are of the order of 103; the distance of the order of 102 km. Both are in agreement with some archaeological findings. An expression for the time spaon required for the development from the earliest stages of urban cultures to the present time is derived and shown to depend on the specific shoreline of the country, that is, the length of the shorline divided by the area of the country. It is pointed out that western Europe's specific shoreline, including land bordering both seas and rivers, is ten times as large as the shoreline area of other parts of the world. It is shown that this greater specific shoreline may account quantitatively for the faster social and technological development of western Europe in the last few centuries. The calculated total span of time of development from earliest urban cultures to our days is found to be of the order of magnitude of ten thousand years. It is shown that the model accounts for the existence at the present time of primitive cultures. A number of suggestions is made in regard to other possible applications of mathematics to history.

Notes: Abstract The assumptions latent in the derivation of the integral equation of Branson are rendered explicit and discussed. It is shown that the equation is valid only for systems in which the substance disappears according to a linear rate law.

Notes: Abstract It is shown that the validity of Branson's integral description of metabolizing systems is subject to severe limitations. The validity is insured only in cases where the reaction is of first order, or quasi of first order. In all other cases Branson's equation has to be modified to insure general applicability. The consequences of a different definition of the metabolizing functionF have also been investigated. With the new definitionF describes the pure effect of metabolization. It is found that in this case the integral equation is only capable of describing first-order reactions. With a slight modification of the integral equation it is possible to describe metabolites “with age”, which do not have reactions of definite order, but which satisfy the superposition principle.

Notes: Abstract A problem in probability is stated with included the problem of the distribution of bacterial mutants as a special case. This problem is solved exactly but since the resulting expressions are too complicated for practical use, various approximate expressions for the distribution are considered, especially for the bacterial mutation case.

Notes: Abstract Simple reaction and discrimination reaction, under the influence of imitation, are considered for the situation in which the stimulus or the stimuli vary slowly with time. The result is analogous to hysteresis under certain conditions. The calculations are facilitated by the solution of $$x = \int_{ - \infty }^{a + \beta x} {g\left( \xi \right)d} \xi ,$$ g(ξ) being the normal error function. Values ofx(α, β) are given in a table.

Notes: Abstract On the basis of simple physical considerations the blood flow in a branching circulatory system is studied. The case of two groups of parallel vessels is treated. The vessels of the same group are supposed to be identical. The resistance of each group is determined by the resistance of each vessel in the group and by the number of vessels in the group. From the dependence of the resistance of each vessel on its radius an expression is obtained for the blood flow through each group of vessels in terms of the numbers and sizes of the vessels in each group. The number of open vessels in an organ and the radius of each of those vessels are assumed to depend on the metabolic rate of that organ. The relations so obtained, together with the expression above, are applied to derive the blood flow through an organ as a function of the metabolic rate of that organ. It is indicated that the relations obtained might describe the shifting of blood from one organ to another if the activity of one of them changes. A way is pointed out to treat neural regulation of this phenomenon.

Notes: Abstract A stochastic model of population growth is treated using the Bellman-Harris theory of agedependent stochastic branching processes. The probability distribution for the population size at any time and the expectation are obtained when it is assumed that there is probability (1−σ), 0≤σ〈1, of the organism dividing into two at the end of its lifetime, and probability σ that division will not take place.

Notes: Abstract In a previous paper, in which a possible mathematical approach to history was outlined, it was shown that urbanization plays an important part in the propagation of new ideas. The rate of such propagation influences the rate of historical developments. The present paper deals in more detail with possible mechanisms of formation of earliest cities. Equations are derived which give the limiting size of such cities and their rate of growth. Of particular importance for the spread of new ideas is the spread of information. The latter largely depends on the fraction of individuals who travel between city and country. Expressions for this quantity are derived. An approach is outlined to the mathematical study of the earliest social classes, which may have been formed as a result of military, religious, or economic stratifications.

Notes: Abstract The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube is computed. The mathematical model is more general than those used in previous analyses: the tube wall properties are realistic; the convective part of the axial inertia forces is taken into account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are present. Four parameters influencing the velocity of propagation are obtained and discussed. Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial inertia effects can be important in laboratory set-ups.

Notes: Abstract A general theory of the drying of frozen tissue is developed and applied to the measurement of the drying rate of frozen guinea pig liver. It is shown that for a given temperature of the subliming ice crystals the mininum drying time of a piece of guinea pig liver is greater than the minimum sublimation time of a piece of ice of the same size and shape by a factor of the order of one thousand. This fact has many implications in the design of freeze-dry apparatus which will be developed in a following paper.

Notes: Abstract The equations governing the time course of the exchange of substances between the blood in the capillaries and the extracellular space are solved for the case of substances which do not penetrate the cells. The equations given relate the time course of the exchange process to the various tissue and circulation parameters such as the specific capillary wall area, the pore area, the inter-capillary distance, the size of the extra-vascular, extra-cellular space, the diffusion coefficient in this space, and the velocity of blood in the capillaries. Some experimental work on capillary exchange is discussed in relation to the theory and estimates are made of the relative importance of the various tissue and circulation parameters in the exchange of substances in different tissues.

Notes: Résumé Les premiers étages sensoriels sont étudiés en utilisant notre modèle de neurone et en supposant que les réseaux responsables de la perception sont particulièrement solides, stables, économiques. Nous montrons que les premiers neurones doivent être spontanément périodiques et autorégulés. La nécessité fonctionnelle des premiers étages de la voie visuelle est démontrée. Par analogie, nous étudions la voie auditive.

Notes: Abstract The problem of the viscous flow of an incompressible Newtonian liquid in a converging tapered tube has been solved in spherical polar coordinates. The method of the solution involves the Stokes' stream function and a technique introduced by Stokes in the study of a sphere oscillating in a fluid. The theory for the flow in a rigid tube includes: (1) the pulsatile flow with both radial and angular velocity components; (2) the steady state flow with both radial and angular velocity components and (3) the very slow steady state flow with only a radial velocity component present. For a tapered elastic tube, the velocity of the propagated pulse wave is determined. The solution given is in terms of the elastic constants of the system and the coordinates for this type of geometry. The pulse velocity is then related to the velocity in an elastic cylindrical tube with the necessary correction terms to account for the tapered tube.

Notes: Abstract In this paper a class of branching processes applicable to populations reproducing by some asexual means or by a simple selfing system of mating is studied. The paper is divided into three parts. In part one the mathematical model is introduced, part two is a mathematical analysis of the model, and in part three concrete applications and examples are given. Many of the proofs of the theorems in part two are omitted but will appear in a subsequent issue of theBulletin.

Notes: Abstract The application of the earlier results (Pavlidis, T. 1965. “A New Model for Simple Neural Nets and its Application in the Design of a Neural Oscillator.”Bull. Math. Biophysics,27, No. 2, 215–229) to the design of more complex neural nets is attempted. The following cases are considered: 1. Chains of neurons where it is proven that the frequency of the output pulses does not depend on the value of the input as long as it is above a certain threshold. 2. Groups of neurons with backward inhibition which present an intermittent mode of operation. 3. Neural nets with periodic facilitation which permit time sharing of certain components for different functions. 4. A neural net which can detect the sign of the input even if the main receptor is sensitive only to the absolute value of it, is presented. 5. A velocity estimating neural net which in combination with one of the nets with intermittent response provides a model for the smooth eye tracking movements.

Notes: Abstract By assigning coordinates to the information space comprising all knowledge, rigorous mathematical interpretations can be placed on such terms as academic ability, memory and creativity such that these psychometric concepts can be incorporated into a framework of functional analysis which then permits the optimization of long-term academic learning processes through the location of the teaching trajectories in information space which will maximize the knowledge accumulated in a generalized educational system composed of a complex of subject-pupil-teacher interactions. The concepts of discrete and continuous information spaces are discussed in connection with subject-subject, subjectpupil and pupil-pupil interactions, and the advantages of using variational versus dynamic programming methods of optimizing alternative educational systems are evaluated.

Notes: Abstract IfK is a partition of a setK which is partially ordered by the relationR andR is a collection of pairs of sets ofK such that the sets of each pair are related byR in the sense of Rashevsky, thenR is a relation which partially ordersK. Necessary and sufficient conditions thatK be a chain are obtained, and ifK is a chain under these conditions, it is shown thatK is unique.

Notes: Abstract This paper continues a comparison of the Taylor series and spherical harmonic forms of multipole representations initiated by Yeh (Bull. Math. Biophysics,24, 197–207, 1962). It is shown that while transformations from Taylor series form into spherical harmonic form is always possible, the inverse cannot be accomplished as suggested by Yeh; corrected transformation equations are given. It is also shown that direct measurement of Taylor coefficients, as outlined in Yeh, Martinek, and de Beaumont (Bull. Math. Biophysics,20, 203–216, 1958), is actually not possible. Accordingly, only the spherical harmonic coefficients can be determined by measurement of surface potentials, as in electrocardiography.

Notes: Abstract This paper is a continuation of a paper, “Some Multi-Dimensional Branching Processes as Motivated by a Class of Problems in Mathematical Genetics I,” by C. J. Mode, which appeared in a previous issue of theBulletin. Its purpose is two-fold; namely to discuss the mathematical existence of the model and to supply the mathematical proofs of some theorems in section two of the paper mentioned above. This paper should be read in conjunction with the previous paper.

Notes: Abstract Rosen (Bull. Math Biophysics. 1959) has argued that a self-reproducing automaton of the type originally described by von Neumann is impossible because of a logical paradox inherent in its definition. The paradox is resolved by explicitly allowing errors (mutations) in the system and thus introducing evolution. There is no paradox in an automaton, originating from a slightly different ancestor through mutation. The von Neumann model thus becomes realistic and useful for a discussion of biological phenomena.

Notes: Abstract Previously proposed formulae for the quantitative estimation of bidirectional shunts across ventricular septal defects require determination of the oxygen contents of mixed venous, pulmonary artery, pulmonary venous, and aortic blood. Because these formulae do not take into account the mixing of oxygenated with unoxygenated blood within the ventricles, their use must result in underestimation of shunt flows in each direction. A mathematical model for a ventricular defect is examined, in which it is assumed that mixing of blood occurs in each of six sites in the venae cavae or right atrium, right ventricle, pulmonary artery, left atrium, left ventricle, and aorta. A total of fourteen streams of blood can flow from one to another of these mixing sites. As long as complete mixing occurs in the six specified mixing sites, any degree of mixing or non-mixing of the various streams is permitted. From the equations characterizing the model, formulae are derived in which the shunt flow in each direction is expressed in terms of the oxygen contents in the six mixing sites and the fractions of blood which enter the shunt from either side without prior mixing in a ventricular mixing site. The previously reported formulae, which apply when no ventricular mixing is allowed to occur, lead to theoretical minimum values for the shunt flows in each direction. At the opposite extreme where all the shunting blood is required to mix in a ventricle before entering the shunt, formulae for maximum possible shunt flows are also obtained. The absolute values for the left-to-right and right-to-left shunt flows, which must lie somewhere between the theoretical maximum and minimum values, cannot be computed from blood gas data alone.

Notes: Abstract Experimental evidence strongly suggests that the contractility of the intact heart in situ, in contrast to that of striated muscle elsewhere in the body, is controlled in a close-cycle system. Thus, the variation of intraventricular pressure during systole follows a complex pattern, whose relative form remains quite constant regardless of the duration of ejection. By use of the single-chambered model of the cardiovascular system, a mathematical representation of a feasible feedback mechanism is developed. The requirement that the feedback system must satisfy mathematical principles eliminates relationships apparently reasonable from a physiological viewpoint. A clinical application which the mathematical development suggests is that early arterial hypertension may arise from an abnormal feedback mechanism with excessively large cardiac output in the initial portion of systole.

Notes: Abstract Due to the lack of direct X-ray evidence for base pairing being the only mechanism for the formation of double helix in a DNA crystal, an alternative explanation is suggested so that the observed DNA loop becomes essential.

Notes: Abstract Validity of group ring expression of selfed population is shown for cases in which there are differences in recombination probabilities between two sexual sides of a plant.

Notes: Abstract An integral equation analysis of generaln compartment steady state systems imbedded in static media of arbitrary complexity has been developed. A set of initial entry functions can be found which serve to determine a corresponding set of partitioned initial entry functions. The partitioned functions, in turn, can be used to predict the probabilities and time courses of various transport histories and to determine all steady state rates of flow between measured compartments. The method is quite general, being completely applicable, for example, to closed systems, to cyclic systems and to systems in which relatively rapid (but finite) exchange between compartments occurs.

Notes: Abstract From the definition of a strong and weakn-ary relation betweenn sets, given in a previous paper (Bulletin of Mathematical Biophysics,27, 477–492), it follows that for a given set ofn sets and givenn-ary relationR between them there can exist only one strong relation, but a large number of weak ones. An expression for the total number of possible weak relations is derived and the notion of the degree of weakness of a relation is introduced and discussed.

Notes: Abstract The problem of determining the sequence of a biopolymer from its fragments is stated in mathematical terms. Using concrete properties of a free monoid, certain general classes of biopolymers are shown to be insolvable from fragment data produced by complete digestion where enzymes specific for any possible combination of chemical bonds are employed.

Notes: Abstract To the extent that all biological phenomena are perceivable only through their physical manifestations, it may be justified to assume that all biological phenomena will be eventually represented in terms of physics; perhaps not of present day physics, but of some “extended” form of it. However, even if this should be correct, it must be kept in mind that representing individual biological phenomena in terms of physics is not the same as deducing from known physical laws the necessity of biological phenomena. Drawing an analogy from pure mathematics, it is possible that while every biological phenomenon may be represented in terms of physics, yet biological statements represent a class of “undecidable” statements within the framework of physics. Such a conjecture is reinforced by the history of physics itself and illustrated on several examples. The 19th century physicists tried in vain todeduce electromagnetic phenomena from mechanical ones. A similar situation may exist in regard to biological and social sciences. Quite generally, the possibility of representing a class B phenomena in terms of class A phenomena does not imply that the phenomena of class B can be deduced from those of class A. The consequences of the above on the relation between physics, biology, and sociology are studied. A tentative postulational formulation of basic biological principles are given and some consequences are discussed. It is pointed out that not only can the study of biological phenomena throw light on some physical phenomena, but that the study of social phenomena may be of value for the understanding of the structures and functions of living organisms. The possibility of a sort of “socionics” is indicated.

Notes: Abstract It is shown that any (ℳ ℛ) has some component which cannot be re-established after it has been inhibited. If there is only one such component, it must be central, that is, its inhibition stops the whole system. These results hold even when it is not assumed that ℳ is connected.

Notes: Abstract In this paper a theory of a class of restricted transition probabilities is developed and applied to a problem in the dynamics of biological populations under the assumption that the underlying stochastic process is a continuous time parameter Markov chain with stationary transition probabilities. The paper is divided into three parts. Part one contains sufficient background from the theory of Markov processes to define restricted transition probabilities in a rigorous manner. In addition, some basic concepts in the theory of stochastic processes are interpreted from the biological point of view. Part two is concerned with the problem of finding representations for restricted transition probabilities. Finally, in part three the theory of restricted transition probabilities is applied to the problem of finding and analyzing some properties of the distribution function of the maximum size attained by the population in a finite time interval for a rather wide class of Markov processes. Some other applications of restricted transition probabilities to other problems in the dynamics of biological populations are also suggested. These applications will be discussed more fully in a companion paper.

Notes: Abstract Equilibrium solubility considerations are presented based on the assumption that equating the kinetic expressionq, developed in part I, to zero can describe the equilibrium or steady state between hydroxyapatite and salt solutions. From this expression is derived Hodge's empirical equilibrium equation,C=KH. Further, a lograithmic transformation of this equation results in an expression that accounts for the equilibrium calcium, phosphorus andpH relation found by Levinskas and Neuman. Finally, it also shows the relation between log (C·P) andpH necessary for typical artificial carious lesions as found by Coolidge, Besic and Jacobs. A discussion of a recent theory of hydroxyapatite solubility of LaMer reveals calculation errors that vitiate his results. It is shown that logK 1 (K 1 is the ratio of the rate constants inq and can serve as a solubility equilibrium constant for hydroxyapatite) varies by only 1.2 units when calculated from three diverse sets of data. This variation is less than that reported by LaMer (when the errors of calculation in that work are corrected) and considerably less than the range of 11 among attempts to calculate a conventionalpK sp , as summarized by Hodge.

Notes: Abstract In imitative behavior, as studied previously by N. Rashevsky (Mathematical Biology of Sociol Behavior, Chapter XIII, The University of Chicago Press, 1950), the reason for the majority of a society to accept a particular behavior is based on purely voluntary action (band-wagon effect). In the present paper effects of coercion of the majority by a small minority group which poses the means for coercion, are studied. More general types of equations are thus obtained and threshold effects found, which bear a resemblance to some such effects studied previously.

Notes: Abstract Part III attempts to develop a diffusion controlled model of caries in the intact enamel employing the kinetic results of the previous two parts. A model of the enamel as a granular bed with a diffusible organic matrix filling the interstices is considered. The basic equations of diffusion and simultaneous reaction are developed under the assumption that all the reactions are so rapid as compared with the diffusion rate, that they are in a quasi-equilibrium state. The resultant system of seven coupled, non-linear parabolic partial differential equations is of such complexity that only numerical solutions could be attempted. Stability restrictions inherent in the problem dictated the use of the DuFort-Frankel numerical solution for parabolic boundary problems. Numerical solutions giving the concentration of all reactants, the rate of mineral loss, and the enamel porosity were obtained for a variety of boundary conditions. It is found that departure from the equilibrium condition expressed in part II is necessary for the occurrence of an attack on the enamel. The rate and pattern of penetration is then determined primarily by the concentrations of undissociated buffer, and salts, together with the rate of diffusion in the surrounding medium. The possibility of a relatively intact surface layer persisting over a demineralized subsurface region due solely to the composition of the demineralizing medium is noted. Remineralization behavior in portions of the carious lesion occurs in the model under certain boundary conditions.

Notes: Abstract Previous work (Macey, 1952) in the application of the one-factor theory to the heart is extended. The rate of production of the excitatory state is assumed to be linear. Two possible mechanisms are indicated whereby such a situation might arise. Assumptions are made regarding the mode of action of the chemical mediators on the heart, and an equation is derived relating the heart rate to the frequency of nerve impulses traveling along the cardiac nerves. This result compares favorably with the experimental findings of A. Rosenblueth and F. A. Simeone (1934). Other experimental results are interpreted in terms of the theory.