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 The theory of complementary variational principles is used to obtain maximum and minimum principles for diffusion problems with Michaelis-Menten kinetics. In an illustrative calculation we obtain an extremely accurate variational solution in good agreement with the numerical solution of McElwain (1978).

Notes: Abstract Necessary and sufficient conditions for primitivity of a product of two Leslie matrices are given. Such a product could be used in modeling the growth of a population governed alternately by two different sets of fertility and survival parameters.

Notes: Abstract Zadeh's transfer function method for linear time-variable systems is used to apply frequency-domain analysis to a periodically time-varying elastance model of the left ventricle. Left ventricular pressure computed from the system function of the time-varying elastance and the phasors of aortic flow shows a typical waveform of the measured ventricular pressure.

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: Abstract A theory of noise fluctuations is developed which is applicable to systems of any size in which unimolecular or bimolecular reactions are occurring. The main difference between small and large reacting systems is that in the former the probability of finding a particle in a particular state does not obey a Gaussian distribution, but satisfies a distribution which reflects the mechanism of the chemical reaction. This difference is reflected in the main result of the theory: an autocorrelation function that is expressible as a sum of exponentials, the amplitudes of which are explicit functions of the moments of the distribution. Thus, by using small systems, the autocorrelation function,in principle, allows the elucidation of reaction mechanisms. Numerical simulations indicate that for reacting systems having ten or fewer particles, the deviation of the autocorrelation function from a single exponential should be easily detectable, and that estimates of the first four moments of the distribution should be possible. Accurate inference of the distribution, however, will require further mathematical and experimental advances.

Notes: Abstract The recent mathematical formalization of the concepts of matter and extrinsical energy, which are used for the relational representation of biological systems, is employed in the analysis of the important experimental discoveries of Comorosanet al. related to low energy electromagnetic irradiations on enzyme substrates. By means of the present analysis one of the properties inherent to the experimental phenomena is more precisely exposed, and theoretical developments corresponding to “energetical evolutions” in a biological system (Leguizamón, 1976) may now have an experimental basis. Important limitations are introduced for the validity of the commutativity and associativity of cartesian product of sets, when they represent matter and its linked extrinsical energy. In connection with this last aspect, new important knowledge is obtained for the relational mathematical representation of biological systems.

Notes: Abstract The structure of solutions to a simple spatially dependent population model involving growth and death is investigated. Two forms of motility of the population are considered: (1) random motion only modeled by a Fickian law, and (2) a directed component of motion (chemotaxis), included in addition to the random motion. Under certain growth conditions a traveling wave of constant speed is approached. This speed can be increased by the addition of the chemotaxis with a corresponding increase in the asymptotic population. Development of initial conditions into a wave is illustrated numerically.

Notes: Abstract This paper describes mechanisms of intracellular and intercellular adaptation that are due to spatial or temporal factors. The spatial mechanisms support self-regulating pattern formation that is capable of directing self-organization in a large class of systems, including examples of directed intercellular growth, transmitter production, and intracellular conductance changes. A balance between intracellular flows and counterflows causes adaptation. This balance can be shifted by environmental inputs. The decrease in Ca2+-modulated outward K+ conductance in certain molluscan nerve cells is a likely example. Examples wherein Ca2+ acts as a second messenger that shunts receptor sensitivity can also be discussed from this perspective. The systems differ in basic ways from recent diffusion models. Chemical transducers driven by membrane-bound intracellular signals can establish long-range intercellular interactions that compensate for variable intercellular distances and are invariant under developmental size changes; diffusional signals do not. The intracellular adaptational mechanisms are formally analogous to intercellular mechanisms that include cellular properties which are omitted in recent reaction-diffusion models of pattern formation. The cellular models use these properties to compute size-invariant properties despite wide variations in their intercellular signals. Mechanisms of temporal adaptation can be derived from the simplest laws of chemical transduction by using a correspondence principle. These mechanisms lead to such properties of intercellular signals as transient overshoot, antagonistic rebound, and an inverted U in sensitivity as intracellular signals or adaptation levels shift. Such effects are implicated in studies of behavioral, reinforcement, motor control, and cognitive coding.

Notes: Abstract Large radiation doses to the lung can cause early death from cardiopulmonary insufficiency resulting from radiation pneumonitis and pulmonary fibrosis. A model for early death following inhalation of insoluble radioactive particles is propose. The model is based on three assumptions: (1) early death results from damage to a cluster of cells from a large number of cell clusters at risk, (2) the dose that causes early death depends on how the radiation is delivered in time and (3) the cell clusters at risk to damage are equally sensitive ro radiation. Results from asymptotic theory of extreme values, along with biophysical considerations, suggest that the cumultive distribution function for the absorbed radiation dose to the production of pulmonary injury sufficient to cause early death is best estimated by the third asymptotic distribution without a threshold. This distribution function is identical to the Weibull cumulative distribution function. Data for Beagle dogs after inhaling relatively insoluble forms of alpha- or beta-gamma-emitting particles are shown to support the Weibull model.

Notes: Abstract Models of the human respiratory tract were developed based on detailed morphometric measurements of a silicone rubber cast of the human tracheobronchial airways. Emphasis was placed on the “Typical Path Lung Model” which used one typical pathway to represent a portion of the lung, such as a lobe, or to represent the whole lung. The models contain geometrical parameters, including airway segment diameters, lengths, branching angles and angles of inclination to gravity, which are needed for estimating inhaled particle deposition. Aerosol depositions for various breathing patterns and particle sizes were calculated using these lung models and the modified Findeisen-Landahl computational scheme. The results agree reasonably well with recent experimental data. Regional deposition, including lobar deposition fractions, are also calculated and compared with results based on the ICRP lung deposition model.

Notes: Abstract The completely symmetrical system is defined as having identical transfer coefficients between pairs of compartments and the same loss coefficient for each compartment. The eigenvalues and eigenvector are explicitly found along with the inverses of the system matrix and the matrix of eigenvectors. Many properties, special instances of more general theorems, can be seen at once from the explicit analytic solution of the initial value, washout and washin problems. The system serves as a known case for testing estimation procedures, algorithms for solutions of linear systems, eigenvalue-eigenvector and inversion routines and is of considerable tutorial value.

Notes: Abstract The mathematical structures underlying the theories of organismic sets, (M, R)-systems and molecular sets are shown to be transformed naturally within the theory of categories and functors. Their natural transformations allow the comparison of distinct entities, as well as the modelling of dynamics in “organismic” structures.

Notes: Abstract To explain the sodium conductance change using Wei's dipole model (Wei, 1969), we may expect that during depolarization the dipole's population difference, ΔN, is first reduced and then returns more slowly to its resting value. This paper shows that the experimental results of gating currents support this idea. Such time course of ΔN, however, is not a usual relaxation process. To account for the unusual behavior of ΔN, we propose two additional assumptions: (1) there exists a special coupling system (probably the intramolecular vibrations) whose coupling strength with the dipoles is much stronger than with the thermal bath (intermolecular vibrations), and (2) there also exist “traps” for the dipole's excitation energy so that this energy is transformed into other energy forms at a rate increasing with the increase of depolarization. Experiments suggest that the traps are proteins located at the inner membrane surface.

Notes: Abstract Current research into the dynamics of iterative ecological and biological models has lead to a number of theorems concerning the existence of various types of iterative dynamical behavior. In particular, much study has been done on the dynamical behavior of the “simplest dynamical system”f b(x)=bx(1−x), which is just the canonical discrete form of logistic growth equations found in ecology, sociobiology, and population biology. In this paper, we make use of some of the techniques and concepts of topological dynamics to construct a number of generalized conjugacy theorems. These theorems are then used to demonstrate that the mappingf b has a number of conjugacy classes in which the dynamics of the iterates is equivalent to within a change of variables. The concepts of fitness and survival in logistic equations are then shown to be independent, if we follow certain intuitive definitions for these concepts. This conclusion follows from a comparison of the conjugacy classes of the functionf b and the extinction sets off b.

Notes: Abstract For chemical reactions not at equilibrium but proceeding in the forward direction in the steady state, a result found by a method first introduced by H. G. Britton (1963, 1965) is generalized to prove that if $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ is the unidirectional flux ratio, $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ exp (−ΔG/RT). The conditions under which the equality or inequality applies are discussed. If the unidirectional fluxes are not in the steady state, the unidirectional flux ratio is time invariant in certain specific situations. One such important case is for chemical reaction systems with an ordered sequence of reactions. For systems with more than one pathway, $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ is not constant except for special cases. These results also apply to diffusional and active transport systems.

Notes: Abstract The nonlinear second-order difference equationx n+1=axn(1-xn−1), where 0≦x nX≦1 anda ≧1, is examined from varying points of view, analytical, numerical and geometrical. An analytic expression is obtained for an invariant attracting curveC ∞ (a) in phase space, which becomes the central object of study. This basic curve, which replaces the simple parabolic shape typical of many analogous first-order models, may have a complicated geometrical structure. As the parametera increases,C ∞(a) undergoes transformations characterized by the dynamical descriptions: stable node→stable focus→stable limit cycle →chaotic attractor. Although the limited characterization ofchaos by the appearance of nonperiodic solutions and solutions of arbitrarily large period is relied upon, this appears to be only a simplified approximation of the real behavior of solutions. Trajectories (x n, xn+1),n=0,1,…, are calculated using the related nonlinear planar mapT a(x,y)=(y,ay(1−x)), and regions of persistence and escape are described for characteristic values ofa. The study of persistence, of even more fundamental interest than the associated problems of periodicity and stability, receives special attention. We introduce a geometrical model, similar in many respects to that for the well-known analoguex n+1=axn(1−x n), but having several new and important features. It appears that as the parametera increases in the chaotic regime there are infinitely many intermittent bursts of increase in the probability that any initial point (x 0, x1) will persist in the unit square under successive iterations of the mappingT a, an unexpected property that should be of interest for applications. A discussion of the applicability of these results to population dynamics theory is given, and it is suggested that such equations might find useful application to problems in developmental biology as well.

Notes: Abstract Based on the principle of minimum power, a mathematical model of the functional state of the oxygen transport system is presented. The optimization model minimizes the power expenditure of the heart, bone marrow, lung and other tissues. The model is used to determine the functional parameters of the oxygen transport system in man under both normal and varying barometric pressures. Theoretical results are compared with experimental data.

Notes: Abstract A quantitative model of ion binding and molecular interactions in the lipid bilayer membrane is proposed and found to be useful in examining the factors underlying such membrane characteristics as shape, sidedness, stability and vesicle size at various cation concentrations. The lipid membrane behaves as a bilayer couple whose preferential radius of curvature depends on the expansion or contraction of one monolayer relative to the other. It is proposed that molecular packing may be altered by electrostatic repulsion of adjacent like-charged phospholipid headgroups, or by bringing two headgroups closer together by divalent cation crossbridging. The surface concentrations of each type of cation-phospholipid complex can be described by simple binding equilibria and the Gouy-Chapman-Stern formulation for the surface potential in a diffuse double layer. The asymmetric distribution of acidic phospholipids in most biological membranes can account for the differential effects of identical ionic environments on either side of the bilayer. The fraction of vesicle material which tends to have a right-side-out orientation may be approximated by a normal distribution about the mean curvature. The theory generates vesicle sidedness distributions that, when fitted to experimental results from human erythrocyte membranes, provide an alternative method of estimating intrinsic cationphospholipid dissociation constants and other molecular parameters of the bilayer. The results also corroborate earlier suggestions that the Gouy-Chapman theory tends to overestimate free counter-ion concentrations at the surface under large surface potentials.

Notes: Abstract The “yellow strips” on the cuticle of the Oriental Hornet (Vespa orientalis, Hymenoptera, Vespinae), present photoelectric properties. A mathematical model for the relative changes in resistance as a photoconductive process conforms to the general model for a semiconductor with traps.

Notes: Abstract Damped nonlinear oscillations in biological and biochemical systems are investigated by the extended Krylov-Bogoliubov-Mitropolskii (KBM) method. A review on the extension made by Popov to the KBM method is given and also further improvements are presented. Applications are made to models of oscillating chemical reactions (Lefever and Nicolis, 1971), FitzHugh (1961) equations, and population dynamics (Gatto and Rinaldi, 1977). Comparison to damped oscillating physical and engineering systems is made.

Notes: Abstract The conditions that will allow the lumping together of several age classes in the Leslie model are investigated. We show that if the lumping is to be valid for all population distributions, then the parameters of the model must be periodic. Lumping is valid when the population is in equilibrium, but equilibrium should be tested before the model is lumped.

Notes: Abstract Catastrophe theory is a mathematical theory which, allied with a new and controversial methodology, has claimed wide application, particularly in the biological and the social sciences. These claims have recently been heatedly opposed. This article describes the debate and assesses the merits of the different arguments advanced.

Notes: Abstract Estimates of capillary tracer permeability calculated using multiple indicator data depend upon the particular model adopted to describe blood tissue exchange. The model proposed by Crone (1963) is appropriate when some of the injected tracer diffuses into the tissue but does not return appreciably to the bloodstream before data collection is terminated. Under these conditions extraction of tracer by the tissue depends on a single dimensionless parameter, αcap, defined as the ratio of capillary permeability surface area to water flow. The effects of finite red cell tracer permeability on the Crone model estimate of capillary permeability are examined in the present study. The results indicate that even when back diffusion from the extravascular space is negligible, significant errors in the Crone model estimate can be expected when capillary permeability is relatively high and the ratio of red cell to capillary permeability is less than unity. However, when an aliquot of blood is equilibrated with tracer prior to injection and the dimensionless capillary permeability is relatively low (i.e. αcap ≦ 0.25 for a haematocrit≦50%), the whole blood Crone model estimate of αcap will be within 10% of the actual value, irrespective of red cell permeability. Red cell-plasma exchange for commonly used tracer-organ combinations should not significantly affect Crone estimates of capillary permeability under normal physiological conditions, but may be important in low flow situations.

Notes: Abstract Assuming truncated ellipsoidal geometry for the right and left ventricles, a model is developed for the myocardium enabling biventricular mechanical behavior to be studied. Employing pressure-volume data taken from normal dog hearts and from hearts in which the pulmonary artery has been banded over periods of 2–40 weeks, it is shown that: (a) right ventricular wall stresses are higher than left ventricular stresses; (b) right ventricular wall stress increases initially to a maximum after 3–4 weeks followed by a decline to normal and even subnormal levels, attaining a minimum value at 32–33 weeks; (c) left ventricular stresses behave in a similar manner, attaining their maximum and minimum levels after 7–8 weeks and 32–33 weeks respectively. These results suggest that surgical or medical therapy in patients with hypertrophied ventricles might be more appropriate during the period of wall stress reduction.

Notes: Abstract In this paper we describe a mathematical model of the oscillations of the diaphragm which limits the vitreous body from the anterior segment of the human eye after the lens has been removed in a cataract operation. We study the motion of this diaphragm driven by movements of the eye. Firstly, a mathematical statement of the problem is given and then we solve the problem exactly for a given class of eye movements. From the analysis we deduce that significant oscillations of the membrane are driven by saccades and that it is the angular acceleration of the eye which causes these types of oscillations. A numerical example is given.

Notes: Abstract The Lotka-Volterra system of prey-predator equations is considered with a special type of continuous time delay. In the case of equal diffusion coefficients Hopf’s bifurcation technique is used to show the existence of travelling wave train solutions for the prey-predator system.

Notes: Abstract A mathematical model of prothrombin activation is being proposed which includes the feedback mechanism of thrombin and the alteration of factor V by thrombin. This model is in good agreement with experimental data for the dependence of the rate of thrombin formation on the concentrations of factors V and X a . In particular, it correctly predicts the existence and location of a maximum in both of these cases.

Notes: Abstract A new model of the upper tracheobronchial tree is proposed to account for the three-dimensional nature of the airway system. In addition to the tube length, the tube diameter, and the branching angle, the model includes information on the orientation angle of each tube relative to its parent tube. The orientation angle, defined as the angle between two successive bifurcations, is useful for calculating the gravitational inclination of each tube. The information on orientation angle is further used to construct a binary coding system for identifying individual tubes in the airway tree. The proposed model is asymmetrical, but the same principles can be readily used to construct a symmetrical one.

Notes: Abstract In any control system for which the number of independent controls is smaller than the number of degrees of freedom to be controlled, our choice of control in any state is restricted to a submanifold of smaller dimension than the tangent space. This simple fact has a number of important consequences for questions of biological import; we consider its implications for adaptation, for senescent phenomena and for the determination of tertiary structures of polypeptides through control of certain average properties. We also formulate the Pontryagin Maximum Principle of Optimal control theory in such a way as to inquire whether specific biodynamic systems can be regarded as optimal with respect to rate of accumulation of particular quantities of the system. We find that if this is possible, the quantity in question must play the role of a clock.

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.

Notes: Abstract A computational model to predict deposition of a wide variety of particulate pollutants in several species of mammals is presented. The model incorporates breathing pattern and detailed anatomical models of the respiratory tract based on extensive morphometric measurements of individual airways. The predicted deposition from this general model is in close agreement with observed deposition of monodisperse aerosols in rats. Particle size and density and respiratory breathing patterns are the critical parameters affecting regional deposition.

Notes: Abstract A theory of antigen-antibody induced particulate aggregation is developed by investigating the stability of model systems of particles. Conditions for the formation of large aggregates are derived by imposing the requirement that at equilibrium a statistically significant number of redundant bonds would occur in a reduced monomer-dimer model system. A relationship is obtained which predicts the fractional agglutination in the reduced dimer system as a function of the antigen, antibody and particulate concentrations: $$\frac{g}{{2f c_0 (1 - g)^{2^ - } }} = \frac{{s_1 }}{r} + \frac{{s_1 s_2 }}{{2!r^2 }} + ... + \frac{{s_1 s_2 ...s_j }}{{j!r^j }},$$ wherec 0 is the initial concentration of monomer,f is a proximity factor,g is the fractional agglutination,s i is the average rate of formation of theith bond from an (i−1)th bound dimer, andr is the average rate of dissociation of a single antibody-antigen bond.

Notes: Abstract The roles of the concentrations of the three interacting constituents in the aggregation process (antibodies, antigens and particulates) are analyzed in detail. It is shown that the basic equation derived in Part I is consistent over a broad range of conditions with experimental findings previously reported.

Notes: Abstract A general mathematical model describing the biochemical interactions of the hormones luteinizing hormone releasing hormone (LHRH), luteinizing hormone (LH) and testosterone (T) in the male is presented. The model structure consists of a negative feedback system of three ordinary differential equations, in which the qualitative behavior is either a stable constant equilibrium solution or oscillatory solutions. A specific realization of the model is used to describe the experimental observations of pulsatile hormone release, its experimental suppression, the onset of puberty, the effects of castration, and several other qualitative and quantitative results. This model is presented as a first step in understanding the physicochemical interactions of the hypothalamic-pituitary-gonadal axis.

Notes: Abstract Based upon the transition rate equation of dipoles in the membrane, we deal with two important aspects of interaction of nerve signals: (1) conditions for nerve excitation and (2) frequency spectrum analysis of nerve impulse. Interrelations between signal amplitudes and frequencies are formulated in detail. There are several important conclusions which can be drawn from our calculations. First, toexcite the nerve, low frequencies are generally more effective than high frequencies. Second, tosedate the nerve (i.e. to suppress undesired activities), high frequencies would suit better. Third, harmonics produced through interactions of nerve signals are not necessarily weaker than the fundamental frequencies. The great significance of our theory is that it indicates in principle the feasibility to alter or rewrite the information contents of a nerve message in our body by applying stimulations of appropriate strengths and frequencies. Thus, the theory provides a physical basis and hence some understanding for a new branch of medicine—neuro therapy such as Nogier's auriculotherapy, Lamy's phonophoresis, Voll's electroacupuncture and the fast rising TENS (transcutaneous electro-neuro stimulation).

Notes: Abstract A new physical property, called resonance of the B-type is hypothetically attached to the λ =546 nm irradiated crystalline (small) molecules. In this respect an up or down configuration is assumed for those states obtained through irradiation times that are multiples of 5 sec. With these assumptions, the cellular receptors that may detect these states appear to belong to three classes: the up, down and alternatively mixed up-down. Using the classic formalism of eigenvectors and eigenvalues, a simple spectroscopic type of formula is derived, through which all the possible states of the above characteristic may be obtained.

Notes: Abstract A model of ecosystems with migration is proposed from the viewpoint of flow. This model explains the following two points: (1) How the density-dependent terms in population dynamics arise as a consequence of migration. (2) How the ecosystem exhibits a hierarchy in energy per unit biomass.

Notes: Abstract For precise boundary conditions of biological relevance, it is proved that the steadily propagating plane-wave solution to the Fisher equation requires the unique (eigenvalue) velocity of advance 2(Df)1/2, whereD is the diffusivity of the mutant species andf is the frequency of selection in favor of the mutant. This rigorous result shows that a so-called “wrong equation”, i.e. one which differs from Fisher's by a term that is seemingly inconsequential for certain initial conditions, cannot be employed readily to obtain approximate solutions to Fisher's, for the two equations will often have qualitatively different manifolds of exact solutions. It is noted that the Fisher equation itself may be inappropriate in certain biological contexts owing to the manifest instability of the lowerconcentration uniform equilibrium state (UES). Depicting the persistence of a mutantdeficient spatial pocket, an exact steady-state solution to the Fisher equation is presented. As an alternative and perhaps more faithful model equation for the propagation of certain species properties through a homogeneous population, we consider a reaction-diffusion equation that features a cubic-polynomial rate expression in the species concentration, with two stable UES and one intermediate unstable UES. This equation admits a remarkably simple exact analytical solution to the steadily propagating plane-wave eigenvalue problem. In the latter solution, the sign of the eigenvelocity is such that the wave propagates to yield the “preferred” stable UES (namely, the one further removed from the unstable intermediate UES) at all spatial points ast→∞. The cubic-polynomial equation also admits an exact steady-state solution for a mutant-deficient or mutant-isolated spatial pocket. Finally, the perpetuating growth of a mutant population from an arbitrary localized initial distribution, a mathematical problem analogous to that for ignition in laminar flame theory, is studied by applying differential inequality analysis, and rigorous sufficient conditions for extinction are derived here.

Notes: Abstract The binding of mono-, di- and trivalent cations to negatively charged surfaces is studied within the framework of a modified Gouy-Chapman equation. For any given combination of ions of the above valences, the existence and uniqueness of the solution for the surface potential is shown. The treatment provides the surface potential and charge density. For a system containing only monovalent and divalent ions, analytical solutions are given. When trivalent ions are also present, a procedure based on numerical integration is described. The distance dependence of the electrostatic potential for planar surfaces is given. The calculations provide the amount of cations tightly bound and the amount trapped in the double layer region. The competition between cations for binding to surfaces is elucidated.

Notes: Abstract This paper deals with a model describing the behavior of barium-treatedApalysia neurons. The model is represented by a dynamical system, so-called “complete system”, defined in R4 and depending on a small parameter. The study of this system under zero membrane current conditions was performed with the use of the qualitative theory of singular perturbations. We show that this system has a stable periodic solution of the discontinuous type when the small parameter tends to 0+. A reduced system defined in R3, associated to the complete system was also studied: it corresponds to a constant activation of the inward current. We demonstrate that the corresponding hypothetical cell remains silent under zero current conditions.

Notes: Abstract In the study of chemical modification of enzymes and other biologically active proteins, plots of fractional residual activity as a function of number of groups modified per enzyme molecule are often used to establish a correlation between the chemical modification and enzyme inactivation reactions and to determine the stoichiometry of the modification reaction. This paper presents a critical examination of the underlying theoretical framework of such graphs. Whereas these plots are usually presented as linear functions, it is shown here that the general equation describing the relationship between inactivation and modification contains an exponential term; therefore, in the general case, the plot is actually a curve. It is suggested that caution be exercised in the interpretation of such plots and that equations such as those derived in the text be used to fit theoretical curves to the data, in order to maximize the information gained from chemical modification experiments.

Notes: Abstract This communiction argues that so-called “hermaphroditic” tracer systems, which are neither open nor closed, do not exist physically. The argument is based on the assumption that any observable (possibly nonhomogeneous) macroscopic compartment can be approximated by a compartmentC with a finite number of entry points for the tracer, each associated with an abstract subcompartment ofC. It is shown that the “hermaphroditic” property requires that the mean waiting time be infinite in at least one of the subcompartments, or in a subcompartment elsewhere in the system. A subcompartment with infinite mean waiting time must have some sort of memory, of infinite duration, which knows how long a given particle has been retained, however long that is, and thereby determines its probability of departure. Assuming, as seems likely, that no physical basis exists for such an infinite memory, it follows that “hermaphroditic” systems do not exist.