SEARCH · PubMed Health

Results for “Computer Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35Linked to original sources

PubMed1991

Different types of rectification at electrical synapses made by a single crayfish neurone investigated experimentally and by computer simulation.

The rectification properties of electrical synapses made by the segmental giant (SG) neurone of crayfish (Pacifastacus leniusculus) were investigated. The SG acts as an interneurone, transmitting information from the giant command fibres (GFs) to the abdominal fast flexor (FF) motoneurones. The GF-SG (input) synapses are inwardly-rectifying electrical synapses, while the SG-FF (output) synapses are outwardly rectifying electrical synapses. This implies that a single neurone can make gap junction hemichannels with different rectification properties. The coupling coefficient of these synapses is dependent upon transjunctional potential. There is a standing gradient in resting potential between the GFs, SG and FFs, with the GFs the most hyperpolarized, and the FFs the most depolarized. The gradient thus biases each synapse into the low-conductance state under resting conditions. There is functional double rectification between the bilateral pairs of SGs within a single segment, such that depolarizing membrane potential changes of either SG pass to the other SG with less attenuation than do hyperpolarizing potential changes. Computer simulation suggests that this may result from coupling through the intermediary FF neurones.

Animals↗

PubMed1985

Computer simulation for hormones related to primary thyropathy.

We propose a mathematical model of the human hypothalamus-anterior pituitary-thyroid system regulating basal metabolism, and practice computer simulation concerning primary thyropathy such as Graves' disease, hypothyroidism, T4-toxicosis and T3-toxicosis by use of this model. In order to throw light on properties of the system, indicial responses of the hormones, T4, T3, rT3, and TSH, and the function of the thyroid gland are computed. Medical treatments for Graves' disease and for hypothyroidism are simulated with a view to enhancing clinical significance. Performance of the simulation leads to an interesting result that when the convertion rate of blood T4 to blood T3 increases, explicit T3-toxicosis occurs, although the function of the thyroid gland is normal.

Computers↗

PubMed1996

Loci of movement of selected points on the femoral head during normal gait. Three-dimensional computer simulation.

Wear of ultrahigh-molecular-weight polyethylene and the subsequent lytic response to the particulate wear debris are the dominant problems in total joint arthroplasty surgery. Wear testing apparatus can play a vital role in the in vitro evaluation of the many factors involved in wear, such as head size, surface roughness, materials for the head, and new materials for the socket. Wear of ultrahigh-molecular-weight polyethylene may be influenced by the wear path. For the related polymer, high-density polyethylene, the wear path is critical to wear magnitude. What is the actual path taken by a single point (or by multiple representative points) on the femoral head of a total hip arthroplasty as it passes through the gait cycle? The goal of this computer simulation study was to trace the paths of specific points on the femoral head as they moved against the polyethylene cup during a single cycle of normal gait to illustrate the motions occurring at the intraarticular surface of the hip joint. This study also yielded unusual data on the "distance traversed" by these points during a single gait cycle. It was found that there was not one path, but rather there were many, and the paths varied widely in both shape and length depending on the location on the femoral head. Moreover, the differences in excursion and direction at different sites during the loaded phase were great. In addition, distances traveled by different points on the femoral head of any given size varied by a factor greater than 2. Most of the points traced quasielliptical paths. This automatically means that the paths of neighboring points cross each other, creating multidirectional shear forces on the acetabular cup surface which may be important in the localization and extent of wear. The plots of traces of the points derived from this study can serve as benchmarks for the ability of hip simulators to reproduce the actual distances and paths of travel of individual points on the femoral head.

Computer Simulation↗

PubMed1994

Computer simulation of brain cooling during cardiopulmonary bypass.

A mathematical model of heat transport was used to analyze the effects of convection, metabolism, and conduction on the rate of brain cooling and the final brain temperature during cardiopulmonary bypass. Convection, a function of cerebral blood flow and arterial blood temperature, is by far the most important process to determine the rate of brain cooling. Arterial blood temperature almost entirely determines the final brain temperature. Although conduction (head surface cooling) has little effect on the rate of brain cooling or final brain temperature in adults, it may have moderate effects in infants. Brain metabolic heat production has insignificant direct effects on the rate of brain cooling and final brain temperature in both adults or infants. Computer simulation of convective cooling of the adult brain to 27 degrees C shows that, with routine perfusion techniques, brain temperature equilibration is rapid (16 minutes) and small brain-blood temperature gradients are achieved. Simulation of infant brain cooling to 17 degrees C shows that, to avoid excessive brain-blood temperature gradients, 22 to 26 minutes may be required to achieve brain temperature equilibration.

Adult↗

PubMed1998

Posterior tilting of the tibial component decreases femoral rollback in posterior-substituting knee replacement: a computer simulation study.

Posterior tilting of the tibial component is thought to increase the range of motion in posterior cruciate-retaining total knee replacement, but its effect on implant motion in posterior cruciate-substituting total knee replacement is unknown. This issue has become of interest recently because manufacturers have introduced instrumentation that produces a posteriorly tilted tibial cut for both implant types. The purpose of this study was to investigate how motion of posterior cruciate-substituting total knee replacement is affected when the tibial component is installed with posterior tilt. Sagittal plane implant motions were predicted from prosthesis geometry with use of a computer simulation in which the femoral condyles were assumed to sit in the bottoms of the tibial condylar wells when the knee was in extension. Rollback of the femoral component was produced by a cam-spine mechanism at higher angles of flexion. The simulations revealed that even small degrees of posterior tilt reduced rollback by limiting the interaction between the cam and spine. Tilting the component posteriorly by 5 degrees caused the cam to contact the spine at a knee flexion angle that was 18 degrees higher than with the untilted component. The results suggest that posterior tilting of the tibial component in posterior cruciate-substituting knee replacement may not produce the same beneficial effects that have been reported for the tilting of tibial components in posterior cruciate-retaining knee replacement.

Arthroplasty, Replacement, Knee↗

PubMed1992

Is the experience with CHART compatible with experimental data? A new model of repair kinetics and computer simulations.

A new incomplete repair model is introduced that differs from previous models of this type by not assuming that repair is complete during long intervals, e.g. "overnight" intervals of 12-24 h. The model was used to assess the risk of myelopathy resulting from continuous hyperfractionated accelerated radiotherapy treatment (CHART) in light of recent experimental data on the rat spinal cord. Model calculations employing biexponential repair kinetics showed that CHART treatments might result in a higher myelopathy risk than an equal dose given in conventional 2-Gy fractions if the parameters obtained from the animal data hold. The probability of observing what has been reported for CHART was determined in computer simulations for different variance scenarios. The chance to observe four myelopathies in the 74 cervical cord patients was estimated to range between 25 and 62%, while the probability to see 0 in 68 thoracic cord patients ranged from 48 to 27%. These numbers were derived from reasonable assumptions about the repair kinetics (e.g. 60% of damage repaired with a half-time of 8 h) so that the over-all probability to observe 4/74 and 0/68 was maximized, and depending on the scenario fell in the range 12-17%. Finally, from these simulations a myelopathy risk of approximately 0.3-1.2% is predicted for the currently employed maximal CHART dose to the spinal cord, i.e. 42 Gy. We conclude that the CHART experience is not compatible with the new experimental data (p < 5%). Incomplete repair is unlikely to be the sole reason for the unexpected toxicity of CHART (p < or = 17%).(ABSTRACT TRUNCATED AT 250 WORDS)

Animals↗

PubMed1996

The effects of nested primer binding sites on the reproducibility of PCR: mathematical modeling and computer simulation studies.

The polymerase chain reaction (PCR) has become an indispensable tool in modern biological research. Although the application of PCR is a standard routine, we widely lack a theoretical understanding of the dynamic processes involved, especially with respect to the amplification of nonreproducible and/or unexpected amplification products. For one potential source of uncertainty, the presence of nested primer binding sites within an amplifyable DNA locus, we consider a simple stochastic model for the dynamics of PCR amplification of competing products. For commonly used thermostable DNA polymerases lacking a 5'-3'-exonuclease activity, we predict the relative amplification frequencies of competing PCR products dependent on the primer binding probability, the number of PCR cycles, and the number of initial DNA template molecules. At low primer binding probabilities and low numbers of initial DNA template molecules and PCR cycles, we expect the amplification of two products. At high primer binding probabilities and/or high copy numbers of initial template molecules only one main amplification product is predicted at increasing cycle numbers. Furthermore, by means of computer simulation studies we quantify the stochastic variation for the amplification frequencies of competing products.

Bacteriophage lambda↗

PubMed1997

Computer simulation of a breast cancer metastasis model.

Recent analysis of relapse data from 1173 untreated early stage breast cancer patients with 16-20 year follow-up shows that frequency of relapse has a double peaked distribution. There is a sharp peak at 18 months, a nadir at 50 months and a broad peak at 60 months. Patients with larger tumors more frequently relapse in the first peak while those with smaller tumors relapse equally in both peaks. No existing theory of tumor growth predicts this effect. To help understand this phenomenon, a model of metastatic growth has been proposed consisting of three distinct phases: a single cell, an avascular growth, and a vascularized lesion. Computer simulation of this model shows that the second relapse peak can be explained by a steady stochastic progression from one phase to the next phase. However, to account for the first relapse peak, a sudden perturbation of the development at the time of surgery is necessary. Model simulations predict that patients who relapse in the second peak would have micrometastases in states of relatively low chemosensitivity when adjuvant therapy is normally administered. The simulation predicts that 15% of T1, 39% of T2, and 51% of T3 staged patients benefit from adjuvant chemotherapy, partially offsetting the advantage of early detection. This suggests that early detection and adjuvant chemotherapy may not be symbiotic strategies. New therapies are needed to benefit patients who would relapse in the second peak.

Breast Neoplasms↗

PubMed1994

A computer simulation study on the input function sampling schedules in tracer kinetic modeling with positron emission tomography (PET).

Tracer kinetic modeling with positron emission tomography (PET) requires measurements of the time-activity curves in both plasma (PTAC) and tissue (TTAC) to estimate physiological parameters, i.e. to fit the parameters of certain compartmental models using PTAC and TTAC as the model input and output functions, respectively. In this paper, we first explored the optimal blood sampling schedule (OBSS) for the input function, based on the tracer [18F]2-fluoro-2-deoxy-D-glucose (FDG) blood sample experimental data. Then using a 5-parameter FDG model we investigated the effects of the plasma sampling schedule, as well as PTAC measurement noise, on the estimation accuracy and reliability of FDG model macro- and micro-parameters and the physiological parameter local cerebral metabolic rates of glucose (LCMRGlc), using computer simulation. Three different methods were used: (a) estimation of the FDG model parameters ignoring PTAC noise using the traditional PTAC schedule (non-OBSS); (b) estimation of the PTAC model parameters and FDG model parameters simultaneously using both non-OBSS and OBSS; (c) estimation of the PTAC model parameters first, then the FDG model parameters using both non-OBSS and OBSS. The results show that OBSS can provide more reliable estimates and largely simplifies the experiment operations.

Blood Specimen Collection↗

PubMed1991

Computer simulation of DNA supercoiling.

Major goals of this research are to comprehend and visualize the detailed three-dimensional arrangements of supercoiled DNA. Attention has been focused in the initial stages on mathematical procedures to generate the spatial coordinates of the B-DNA double helix constrained to specific spatial pathways and on simple energy models of chain conformation. The new treatment of superhelical DNA in terms of parametric curves is an important first step in being able to generate and examine tertiary structure systematically. The location of every residue is implicitly determined by the equation of the closed curve, with the number of computational variables sharply reduced compared to the number required for explicit specification of all chain units. Furthermore, the constraints of ring closure in cyclic chains and/or the end-to-end limitations on constrained open chains are automatically satisfied by the formulations (cubic B-splines and finite Fourier series) chosen in this work. The predicted conformations of elastic DNA do not appear to be tied to either the form of chain representation or the computer simulation method. Significantly, two very different minimization and modeling approaches come to the same structural conclusions. The most stable configurations of the closed circular elastic DNA model are found to be interwound superhelices that are critically dependent on the specified linking number difference. The total elastic energy is proportional to the imposed linking number difference, and beyond the critical linking number difference separating the circular and figure-eight forms, the writhing number of the DNA superhelices is directly proportional to delta Lk. The measured proportionality constant between Wr and delta Lk, however, is somewhat greater than that deduced from experimental observations of plectonemically interwound DNA chains and an assumed structural model. Furthermore, at large delta Lk, the interwound structures appear to curve. The treatment of the DNA double helix as an ideal elastic rod is clearly incorrect. The chain cannot bend with the same ease in all directions. The degree of bending observed in atomic level models is also tied to the angular twist so that the presumed partitioning of bending and twisting components is in error. Furthermore, the local chain bending and twisting are base sequence dependent, with certain residues able to flex more symmetrically than others. The polyelectrolyte character of the DNA is additionally expected to govern the overall folding of the chain and to influence the local secondary structure. The next step in this work is to compare the properties of such "real" DNA with conventional elastic models.(ABSTRACT TRUNCATED AT 400 WORDS)

Computer Simulation↗

PubMed1980

The computer leucocyte. Analysis of the random movement of leucocytes in a visual field by means of computer simulation.

In vital preparations moving polymorphonuclear leucocytes (PMNs) disappear from the visual field. The present study provides evidence by means of geometrical analysis and computer simulation, that this disappearance is due to a heterogeneous distribution of the PMNs at the beginning of observation. The path of random-moving computer leucocytes is more tortuous than the path of comparable vital PMNs without obvious attractant, i.e. vital PMNs supposed to be random-moving. This finding is in agreement with the Allan-Wilkinson observation of 'persistent random walk' of vital cells.

Cell Movement↗

PubMed1998

Solution conformation of the Pseudomonas syringae MSU 16H phytotoxic lipodepsipeptide Pseudomycin A determined by computer simulations using distance geometry and molecular dynamics from NMR data.

Pseudomycin A is a cyclic lipodepsinonapeptide phytotoxin produced by a strain of the plant pathogenic bacterium Pseudomonas syringae. Like other members of this family of bacterial metabolites, it is characterised by a fatty acylated cyclic peptide with mixed chirality and lactonic closure. Several biological activities of Pseudomycin A are lower than those found for some of its congeners, a difference which might depend on the diverse number and distribution of charged residues in the peptide moiety. Hence, it was of interest to investigate its conformation in solution. After the complete interpretation of the two-dimensional NMR spectra, NOE data were obtained and the structure was determined by computer simulations, applying distance geometry and molecular dynamics procedures. The conformation of the large ring of Pseudomycin A in solution includes three rigid structural regions interrupted by three short flexible regions that act as hinges. The overall three-dimensional structure of the cyclic moiety is similar to that of previously studied bioactive lipodepsinonapeptides produced by other pseudomonads.

Computer Simulation↗

PubMed2000

Scheduling surgical cases into overflow block time- computer simulation of the effects of scheduling strategies on operating room labor costs.

UNLABELLED: "Overflow" block time is operating room (OR) time for a surgical group's cases that cannot be completed in the regular block time allocated to each surgeon in the surgical group. Having such overflow block time increases OR utilization. The optimal way to schedule patients into a surgical group's overflow block time is unknown. In this study, we developed a scheduling strategy that balances the OR manager's need to reduce staffing costs and the needs of patients and surgeons for flexibility in choosing the dates and times of cases. We used computer simulation to evaluate our scheduling strategy. Surgeons and patients (i) can schedule the case into any overflow block within 2 wk; (ii) can only schedule the case into a "first case of the day" start time more than 2 wk in the future if there is not enough open time for the case within 2 wk; (iii) must schedule the case to be done within 4 wk; and (iv) are encouraged to perform the case on the earliest possible date. Staffing costs were lowest when the OR manager did not incorporate surgeon and patient preferences when scheduling cases into overflow block time. The strategy we developed provides surgeons and patients with some flexibility in scheduling, while only increasing OR staffing costs slightly over the minimum achieved when the OR manager controls scheduling. IMPLICATIONS: The strategy we developed provides surgeons and patients with some flexibility in scheduling, while increasing OR staffing costs only slightly over the minimum achieved when the OR manager controls scheduling. Staffing costs were lowest when the operating room (OR) manager did not incorporate surgeon and patient preferences when scheduling cases into overflow block time.

Appointments and Schedules↗

PubMed1993

Computer simulations of NMDA and non-NMDA receptor-mediated synaptic drive: sensory and supraspinal modulation of neurons and small networks.

1. The segmental locomotor network in lamprey can generate the rhythmic burst pattern underlying locomotion when it is driven via synaptic glutamate receptors. Lower rates of activity can be evoked by activation of N-methyl-D-aspartate (NMDA) receptors, whereas a rapid activity can only be induced by non-NMDA receptors [kainate/alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)]. The reticulospinal and sensory inputs are known to act via both NMDA and non-NMDA receptors, but it is unclear how these inputs can provide an appropriate control of the locomotor rate. We have examined the effects of different types of excitatory synaptic input to neurons of the locomotor network with the use of a computer-simulated electrical neuron model, with Na+, K+, Ca(2+)-dependent K+ channels, and with inherent oscillatory properties linked to the NMDA conductance. Synapses were modeled as a modulated ionic conductance in the membrane of the postsynaptic cell comprising a voltage-dependent NMDA component (Na+, K+, Ca2+ conductances) of long duration, and/or a non-NMDA component (Na+, K+ conductance) of short duration. 2. By using two neurons to drive a postsynaptic cell with non-NMDA-type synapses, a continuous range of firing frequencies could be evoked in the postsynaptic cell, by altering the firing rate of the presynaptic cells. If a single presynaptic neuron was used, there was a tendency toward spike synchronization between the pre- and postsynaptic cells. 3. When a postsynaptic neuron was driven via NMDA synapses, an oscillatory burst activity could be evoked. The rate of the oscillations was, however, little affected by the presynaptic firing rate. When a drive neuron with mixed (NMDA and non-NMDA) synapses was used, the rate of the oscillations could be changed within a limited frequency range by altering the presynaptic firing rate. By adding another larger drive neuron, having a larger rheobase current and mixed synapses with smaller relative NMDA components, the frequency range of the postsynaptic oscillations could be markedly increased. The frequency range depended on the parameters selected for each of the two types of mixed synapses. 4. A small rhythm-generating neuronal network, comprising six cells connected as the principal interneurons of the lamprey spinal locomotor network, was used to test the role of a tonic NMDA and non-NMDA receptor activation to drive the network and produce bursting.(ABSTRACT TRUNCATED AT 400 WORDS)

Animals↗

PubMed1998

Cryptococcus neoformans chemotyping by quantitative analysis of 1H nuclear magnetic resonance spectra of glucuronoxylomannans with a computer-simulated artificial neural network.

The complete assignment of the proton chemical shifts obtained by nuclear magnetic resonance (NMR) spectroscopy of de-O-acetylated glucuronoxylomannans (GXMs) from Cryptococcus neoformans permitted the high-resolution determination of the total structure of any GXM. Six structural motifs based on an alpha-(1-->3)-mannotriose substituted with variable quantities of 2-O-beta- and 4-O-beta-xylopyranosyl and 2-O-beta-glucopyranosyluronic acid were identified. The chemical shifts of only the anomeric protons of the mannosyl residues served as structure reporter groups (SRG) for the identification and quantitation of the six triads present in any GXM. The assigned protons for the mannosyl residues resonated at clearly distinguishable positions in the spectrum and supplied all the information essential for the assignment of the complete GXM structure. This technique for assigning structure is referred to as the SRG concept. The SRG concept was used to analyze the distribution of the six mannosyl triads of GXMs obtained from 106 isolates of C. neoformans. The six mannosyl triads occurred singularly or in combination with one or more of the other triads. The identification and quantitation of the SRG were simplified by using a computer-simulated artificial neural network (ANN) to automatically analyze the SRG region of the one-dimensional proton NMR spectra. The occurrence and relative distribution of the six mannosyl triads were used to chemotype C. neoformans on the basis of subtle variations in GXM structure determined by analysis of the SRG region of the proton NMR spectrum by the ANN. The data for the distribution of the six SRGs from GXMs of 106 isolates of C. neoformans yielded eight chemotypes, Chem1 through Chem8.

Bacterial Typing Techniques↗

PubMed2000

How does GAP catalyze the GTPase reaction of Ras? A computer simulation study.

The formation of a complex between p21(ras) and GAP accelerates the GTPase reaction of p21(ras) and terminates the signal for cell proliferation. The understanding of this rate acceleration is important for the elucidation of the role of Ras mutants in tumor formation. In principle there are two main options for the origin of the effect of GAP. One is a direct electrostatic interaction between the residues of GAP and the transition state of the Ras-GAP complex and the other is a GAP-induced shift of the structure of Ras to a configuration that increases the stabilization of the transition state. This work examines the relative importance of these options by computer simulations of the catalytic effect of Ras. The simulations use the empirical valence bond (EVB) method to study the GTPase reaction along the alternative associative and dissociative paths. This approach reproduces the trend in the overall experimentally observed catalytic effect of GAP: the calculated effect is 7 +/- 3 kcal/mol as compared to the observed effect of approximately 6.6 kcal/mol. Furthermore, the calculated effect of mutating Arg789 to a nonpolar residue is 3-4 kcal/mol as compared to the observed effect of 4.5 kcal/mol for the Arg789Ala mutation. It is concluded, in agreement with previous proposals, that the effect of Arg789 is associated with its direct interaction with the transition state charge distribution. However, calculations that use the coordinates of Ras from the Ras-GAP complex (referred to here as Ras') reproduce a significant catalytic effect relative to the Ras coordinates. This indicates that part of the effect of GAP involves a stabilization of a catalytic configuration of Ras. This configuration increases the positive electrostatic potential on the beta-phosphate (relative to the corresponding situation in the free Ras). In other words, GAP stabilizes the GDP bound configuration of Ras relative to that of the GTP-bound conformation. The elusive oncogenic effect of mutating Gln61 is also explored. The calculated effect of such mutations in the Ras-GAP complex are found to be small, while the observed effect is very large (8.7 kcal/mol). Since the Ras is locked in its Ras-GAP configuration in our simulations, we conclude that the oncogenic effect of mutation of Gln61 is indirect and is associated most probably with the structural changes of Ras upon forming the Ras-GAP complex. In view of these and the results for the Ras' we conclude that GAP activates Ras by both direct electrostatic stabilization of the transition state and an indirect allosteric effect that stabilizes the GDP-bound form. The present study also explored the feasibility of the associative and dissociative mechanism in the GTPase reaction of Ras. It is concluded that the reaction is most likely to involve an associative mechanism.

Computer Simulation↗

PubMed1983

Babesia bovis: computer simulation of the relationship between the tick vector, parasite, and bovine host.

Components contributing to the survival of Babesia bovis in its tick vector (Boophilus microplus) and bovine host were identified and their relationship simulated with an interactive BASIC computer program. Estimates for major independent variables (host-finding success, feeding success, filial infection rate, recovery rate) were derived from published laboratory and field data. Values for dependent variables (tick burden, inoculation rate, tick and bovine infection rates) were calculated during successive tick generations until equilibrium conditions were reached. The simulation was used to predict the effects of alternate disease control strategies. Babesiosis outbreaks were associated with simulated daily infestations of approximately two to eight engorged ticks. Below this range, B. bovis parasites disappeared, while above it all cattle were infected during the period of calfhood resistance (through 9 months of age). The babesial inoculation rate was most unstable over the above range due to instability in tick populations and tick infection rates. The risk of babesiosis outbreaks among older, less resistant cattle was also greatest over this range. The optimum daily infestation was eight to nine engorged ticks. This level would maintain herd immunity to babesiosis without causing significant physiologic stress to the bovine host. The difficulty in maintaining ecological and immunological stability at low infestation rates discourages "strategic dipping" as a means of maintaining herd immunity to babesiosis. Introduction of tick-resistant cattle provided a more effective means of reducing the risk of babesiosis outbreaks. The economic implications of computer-simulated babesiosis control alternatives should be analyzed. The validity of simulation parameters should be confirmed by field studies on the relationship between tick burdens and the babesial inoculation rate. These studies would be greatly facilitated by development of an assay system for the detection of babesiae in field-collected ticks.

Aging↗

PubMed1993

Mathematical description and computer simulation of retinal cometlike afterimages: a modified neural equation with stability analysis.

A mathematical model for the spatiotemporal description of a well-known psychophysical phenomenon, the cometlike afterimage effect (CLAIE), is presented. The CLAIE occurs when a bright circular light spot moves slowly in the peripheral human retina. Under these conditions, the leading edge of the dot looks circular, but the trailing edge becomes elongated like a comet's tail whose length increases with speed and luminance, and the illusion is more prominent for photopic backgrounds. This cometlike motion smear is described on the basis of the temporal responsiveness and adaptation of rods. The model is an extension of an existing neural model of M. N. Oğuztöreli et al., with an additional term that allows prolonged saturation and long decay time following exposure to intense stimuli, and these effects are held responsible for the cometlike smear. The model predicts the response of photoreceptors through a nonlinear ordinary integrodifferential equation, which includes known biophysical terms for response dynamics, adaptation, saturation, and kinetics of intermediate components of the phototransduction process. The introduction of a saturation coefficient into the neural equation makes it possible to distinguish the different saturation thresholds of the rod-and-cone system. Numerical determination of the stationary solutions and complete linear stability analysis of the improved neural equation are given for a neuron of second order, and some computational results are presented for phase flows around different singular points in the phase field. A computer simulation based on the improved neural equation is presented for modeling the development and features of the CLAIE as a function of the speed and luminance of the stimulus and the background intensity. The computational results agree well with the psychophysical findings relating to the CLAIE.

Afterimage↗