Fetal rights

From Wikipedia, the free encyclopedia

https://en.wikipedia.org/wiki/Fetal_rights 

Fetal rights are the moral rights or legal rights of the human fetus under natural and civil law. The term fetal rights came into wide usage after Roe v. Wade, the 1973 landmark case that legalized abortion in the United States. The concept of fetal rights has evolved to include the issues of maternal substance use disorders, including alcohol use disorder and opioid use disorder. Most international human rights charters "clearly reject claims that human rights should attach from conception or any time before birth." While international human rights instruments lack a universal inclusion of the fetus as a person for the purposes of human rights, the fetus is granted various rights in the constitutions and civil codes of several countries.

History

See also: History of abortion law debate

In antiquity, the fetus was sometimes protected by restrictions on abortion. Some versions of the Hippocratic Oath indirectly protected the fetus by prohibiting abortifacients. Until approximately the mid-19th century, philosophical views on the fetus were influenced in part by Aristotelian concept of delayed hominization. According to it, human fetuses only gradually acquire their souls, and in the early stages of pregnancy the fetus is not fully human. Relying on examinations of miscarried fetuses, Aristotle believed that male fetuses acquire their basic form at around day 40, and female ones at day 90. For Pythagoreans, however, fetal life was co-equal in moral worth with adult human life from the moment of conception; similar views were held by Stoics. Ancient Athenian law did not recognise fetal right to life before the ritual acknowledgement of the child. The law, however, allowed for the postponement of the execution of sentenced pregnant women until a baby was delivered.

Several Hindu texts on ethics and righteousness, such as Dharmaśāstra, give fetus a right to life from conception, although in practice such texts are not always followed.

The property law of the Roman Empire granted fetus inheritance rights. As long as the fetus was conceived before the testator's death (usually, the father) and then born alive, their inheritance rights were equal to those born before the testator's death. Even though under Roman law the fetus was not a legal subject, it was a potential person whose property rights were protected after birth. Roman jurist Ulpian noted that "in the Law of the Twelve Tables he who was in the womb is admitted to the legitimate succession, if he has been born". Another jurist Julius Paulus Prudentissimus similarly noted, that "the ancients provided for the free unborn child in such a way that they preserved for it all legal rights intact until the time of birth". The inheritance rights of the fetus were means of fulfilling the testator's will. The interests of the fetus could be protected by a custodian, usually a male relative, but in some cases a woman herself could be appointed the custodian. The Digest granted the fetus consanguinity rights, vesting the protection of fetal interests in the praetor. The Digest also prohibited the execution of pregnant women until delivery. The Roman law also envisaged that if a slave mother had been free for any period between the time of the conception and childbirth, the child would be regarded as born free. Although the mother might have become slave again before the childbirth, it was considered that the unborn should not be prejudiced by the mother's misfortune. At the same time, Greek and Roman sources do not mention issues of alcohol consumption by pregnant women. On that basis it is believed that Greeks and Romans were not aware of the fetal alcohol syndrome.

A baby holding the "Petition of the Unborn Babes", 18th century illustration.

After the spread of Christianity, an issue emerged: whether it was permissible for a pregnant woman to be baptised before childbirth, due to uncertainty as to whether the fetus would be cobaptised with its mother. The Synod of Neo-Caesarea decided that the baptism of a pregnant woman in any stage of gestation did not include the fetus. In the Middle Ages, fetal rights were closely associated with the concept of ensoulment. In some cases the fetus could also inherit or be in the order of succession. In the Byzantine Empire, a fetus was regarded as a natural person and could inherit alongside blood descendants and slaves. Byzantine Emperor Michael VIII Palaiologos allowed soldiers to transfer their pronoiai to their unborn children. The unborn royals were increasingly granted the right to succession. In 1284, King of Scotland Alexander III designated his future unborn children as heirs presumptive by the act of parliament to avoid potential squabbles among loyal descendants of his lineage. The 1315 entail of Scottish king Robert the Bruce allowed the unborn collateral individuals to be in line for the throne beyond his brother Edward and daughter Marjorie Bruce. After the death of Albert II of Germany in 1439, his then-unborn son Ladislaus the Posthumous inherited his father's sovereign rights. In 1536, the British Parliament gave the unborn children of Henry VIII and Jane Seymour precedence in the line of royal succession. The medieval distinction between the ensouled and the unensouled fetus was removed after Pope Pius IX decreed in 1854 that the ensoulment of Virgin Mary occurred at conception.

In 1751, a pamphlet "The Petition of the Unborn Babes to the Censors of the Royal College of Physicians of London" by physician Frank Nicholls was published, advocating fetal right to life and protection. The pamphlet anticipated many of the arguments of the 21st century's pro-life movement. In 1762, English jurist and judge William Blackstone wrote that an "infant in its mother's womb" could benefit from a legacy and receive an estate as if it were actually bom. The fetus was thus considered a person for purposes of inheritance. Similarly to the Roman law, the Napoleonic Code envisaged that if a woman becomes a widow, a male guardian should be appointed for her unborn child.

In the 20th century and particularly after World War II fetal rights issues continued to develop. In 1948, the Declaration of Geneva was adopted which prior to amendments in 1983 and 2005, advised physicians to "maintain the utmost respect for human life from the time of its conception". In 1967, American Bar Association Journal noted "the modern trend of legal decisions that grant every property and personal right to the unborn child, including the right to life itself, from conception on". In 1975, while interpreting the right to life under the Basic Law of Germany, the Federal Constitutional Court opined that "life in the sense of historical existence of a human individual" exists "at least from the 14th day after conception (nidation, individuation)" and thus everyone's right to life under the Basic Law of Germany includes the unborn as human beings. The 1980s witnessed the reappearance of fetal protection in the workplace, aimed at guarding fetal health in potentially hazardous working conditions. In 1983, Ireland was one of the first countries in the world to constitutionalize a fetal right to life by passing the Eighth Amendment to the Constitution, later repealed in September 2018.

Modern regulations

The only modern international treaty specifically tackling the fetal rights is the American Convention on Human Rights which envisages the fetal right to life from the moment of conception. The convention was ratified by twenty five countries of the Americas (two countries later denounced the convention leading the current number of ratifiers to be twenty three) in 1973–1993. Mexico ratified the convention with the reservation that the expression "in general" concerning the fetal right to life does not constitute an obligation and that this matter falls within the domain of the states. While the convention may be interpreted to permit domestic abortion laws in exceptional circumstances, it effectively declares the fetus a person. However, only a minority of state ratifiers completely prohibit abortion without allowing for an exception when the pregnant woman's life is in danger (Dominican Republic, El Salvador and Nicaragua).

Based on the 1959 Declaration of the Rights of the Child, preambular paragraph 9 of the Convention on the Rights of the Child (CRC) states that "the child... needs... appropriate legal protection before as well as after birth", but due to ambiguity the legal protection of the fetus conflicts with the rights of a pregnant girl under the same Convention. Such conflict is sometimes called maternal-fetal conflict. Under CRC, the rights of a pregnant girl are interpreted as superseding those of her fetus. The states retain the power to decide for themselves what prenatal legal protection they would adopt under CRC. A proposal to grant fetus the right to life from conception was put forward by Belgium, Brazil, El Salvador, Mexico and Morocco during drafting of the International Covenant on Civil and Political Rights (ICCPR), but it was rejected in favor of less stringent wording. At the same time, ICCPR prohibits the execution of pregnant women.

The World Medical Association Declaration on Therapeutic Abortion notes that "circumstances bringing the interests of a mother into conflict with the interests of her unborn child create a dilemma and raise the question as to whether or not the pregnancy should be deliberately terminated". The Dublin Declaration on Maternal Health, signed in 2012, prioritizes fetal right to life by noting that "there is a fundamental difference between abortion, and necessary medical treatments that are carried out to save the life of the mother, even if such treatment results in the loss of life of her unborn child". Several organizations, such as World Health Organization (WHO) and Human Rights Watch prioritize women's reproductive rights over fetal rights.

Under European law, a fetus is generally regarded as an in utero part of the mother and thus its rights are held by the mother. The European Court of Human Rights opined that the right to life does not extend to fetuses under Article 2 of the European Convention on Human Rights (ECHR). In H. v. Norway, the European Commission did not exclude that "in certain circumstances" the fetus may enjoy "a certain protection under Article 2, first sentence". Two European Union member states (Hungary and Slovakia) grant the fetus the constitutional right to life. The Constitution of Norway grants the unborn royal children the right of succession to the throne. In English common law, fetus is granted inheritance rights under the born alive rule.

Islamic law grants the fetus the right to life particularly after ensoulment, which according to various Islamic jurists happens after 40–42 days or four months after conception (some Shiite jurists believe the ensoulment occurs after 11 to 14 days, during the implantation of the fertilized egg in the uterine wall). Both the Sunni and Shiite jurists accord the fetus inheritance rights under two conditions: if a man dies and a pregnant wife survives him, the fetal right to inherit is secure and the inheritance cannot be disposed of before the fetus' share is set aside. Under the second condition, if a woman aborts the fetus at any stage and ignores any vital signs, the fetus is entitled to the inheritance of any legitimate legator who dies after its conception.

The legal debate on fetal rights sometimes invokes the notion of fetal viability. Its primary determinant is fetal lung capacity which typically develops at twenty-three to twenty-four weeks. The twenty-three weeks is usually regarded as the lower bound of fetal viability because technology has been unable to surpass the limit set by lung development. It was nonetheless stated that technology has made it possible to regard the fetus as a patient independent of the mother. In Winnipeg Child and Family Services v. G., the judges argued that "technologies like real-time ultrasound, fetal heart monitors and foetoscopy can clearly show us that the fetus is alive" and thus the born alive rule is "outdated and indefensible".

The creation of human embryos for all research purposes is prohibited by the Convention for the Protection of Human Rights and Dignity of the Human Being with regard to the Application of Biology and Medicine. However, similarly to the abortion debate, in the normative debate on embryo research two views can be distinguished: a "fetalist" view focusing on the moral value of the embryo, and a "feminist" view advocating the interests of women, particularly candidate oocyte donors.

Fetal rights by country

Country	Constitutional protection of fetal rights		Recognition of personhood
Canada	No		No	223. When child becomes human beingA child becomes a human being within the meaning of this Act when it has completely proceeded, in a living state, from the body of its mother , whether or not: (a) it has breathed; (b) it has an independent circulation; or (c) the navel string is severed.
Chile	Yes	Article 19 The Constitution guarantees all persons: 1.The right to life and to the physical and mental integrity of the person. The law protects the life of the unborn.	Yes
Dominican Republic	Yes	Article 37The right to life is inviolable from conception until death. The death penalty may not be established, pronounced, nor applied in any case.	Yes
Ecuador	Yes	Article 45Children and adolescents shall enjoy the rights that are common to all human beings, in addition to those that are specific to their age. The State shall recognize and guarantee life, including care and protection from the time of conception.	Yes
El Salvador	Yes	Article 1El Salvador recognizes the human person as the origin and the end of the activity of the State, which is organized to attain justice, judicial security, and the common good. In that same manner, it recognizes as a human person every human being since the moment of conception.	Yes
Guatemala	Yes	Article 3The State guarantees and protects the human life from its conception, as well as the integrity and security of the person.	Yes
Hungary	Yes	Article 2Human dignity shall be inviolable. Every human being shall have the right to life and human dignity; the life of the foetus shall be protected from the moment of conception.	Yes
Honduras	Yes	Article 67The unborn shall be considered as born for all rights accorded within the limits established by law.	Yes
Madagascar	Yes	Article 19The State recognizes and organizes for all individuals the right to the protection of health from their conception through the organization of free public health care, which gratuitousness results from the capacity of the national solidarity.	Yes
Peru	Yes	Article 2.To life, his identity, his moral, psychical, and physical integrity, and his free development and well-being. The unborn child is a rights-bearing subject in all cases that benefit him.	Yes
Brazil	No		Yes	Article 2. The civil personality of the person starts in the birth with the life, but the law safeguard, since the conception, the rights of the unborn.
Philippines	Yes	Section 12 The State recognizes the sanctity of family life and shall protect and strengthen the family as a basic autonomous social institution. It shall equally protect the life of the mother and the life of the unborn from conception.	Yes
Slovakia	Yes	Article 151. Everyone has the right to life. Human life is worthy of protection already before birth.	Yes
Solomon Islands	No		Yes

Argentina, Barbados, Bolivia, Brazil, Chile, Colombia, Costa Rica, Dominica, Ecuador, El Salvador, Grenada, Guatemala, Haiti, Honduras, Jamaica, Mexico, Nicaragua, Panama, Paraguay, Peru, Dominican Republic, Suriname, Trinidad and Tobago (denounced), Uruguay and Venezuela (denounced).

 

United States President George W. Bush signs the Unborn Victims of Violence Act of 2004

The Eighth Amendment of the Constitution of Ireland gave "the unborn" a right to life equal to that of "the mother". In 2018, the Supreme Court ruled that the fetus' only inherent constitutionally protected right is the right to be born, overturning a High Court ruling that a fetus additionally possessed the children's rights guaranteed by Article 42A of the Constitution. On 25 May 2018, a referendum was passed which amended the Constitution by the substitution of the former provision recognising the right to life of the unborn, with one permitting the Oireachtas, the Irish Parliament, to legislate for the termination of pregnancies. This amendment took effect when it was signed into law by the President of Ireland on 18 September 2018, and abortion was governed by the Protection of Life During Pregnancy Act 2013 until it was replaced and repealed by the Health (Regulation of Termination of Pregnancy) Act 2018, which took effect on 1 January 2019.

In the United States, as of 2014, thirty-eight states provide certain level of criminal protection for the unborn, and twenty-three of these states have laws that protect the fetus from conception until birth. All US states–by statute, court rule or case law–permit a guardian ad litem to represent the interests of the unborn. In 1999, the Unborn Victims of Violence Act was introduced into United States Congress which defines violent assault committed against pregnant women as being a crime against two victims: the woman and the fetus she carries. This law was passed in 2004 after the murder of Laci Peterson and the fetus she was carrying. In 2002, U.S. President George W. Bush announced a plan to ensure health care coverage for fetuses under the State Children's Health Insurance Program (SCHIP).

The civil codes of several countries, such as China (including Hong Kong and Macau) and Russia, as well as some US states, grant fetus inheritance rights, usually under the born alive rule. In the civil code of Iran, fetus can inherit in case of abortion that took place due to a crime, as long as the fetus was alive even for a second after birth. Under the civil code of Japan, for the purposes of inheritance the fetus is deemed to have already been born. The civil codes of the Philippines and Spain envisage that donations to the unborn children can be made and accepted by "persons who would legally represent them if they were already born". The same is allowed by the Malikis.

Alongside Norway, the Constitution of Bhutan grants the unborn royal children the right to succession, but only if there is no male heir.

Behavioral intervention

Various initiatives, prompted by concern for the ill effects which might be posed to the health or development of a fetus, seek to restrict or discourage women from engaging in certain behaviors while pregnant. Also, in some countries, laws have been passed to restrict the practice of abortion based upon the gender of the fetus.

• Many jurisdictions actively warn against the consumption of alcoholic beverages by pregnant women, recommending a maximum intake or total abstinence, due to its association with fetal alcohol syndrome. Countries that encourage those who are pregnant to avoid alcohol either entirely or partially include Australia, Canada, France, Iceland, Israel, the Netherlands, Norway, New Zealand, Spain, the United Kingdom, and the United States.
• Many national and international agencies recommend dietary guidelines for pregnant women due to the health risks posed by the consumption of fish contaminated with methylmercury through industrial pollution. Studies have linked exposure to various levels of methylmercury in utero to neurological disorders in children.
• The use of tobacco products or exposure to secondhand smoke during pregnancy has been linked to low birth weight. Governor Mike Huckabee of Arkansas, citing studies which attribute 10% of infant deaths to tobacco-smoking mothers, considered adopting a smoking ban for pregnant women in 2006 with the aim of reducing infant mortality.
• No U.S. state has enacted a law which criminalizes specific behavior during pregnancy, but, nonetheless, it has been estimated that at least 200 American women have been criminally prosecuted or arrested under existing child abuse statutes for allegedly bringing about harm in-utero through their conduct during pregnancy. Reasons for pressing charges included use of illicit substances, consumption of alcohol, and failure to comply with a doctor's order of bedrest or caesarean section. Drug addicts have been accused of "supplying drugs to a minor" through unintentional chemical subjection via the umbilical cord. Others have been charged with assault with a deadly weapon with the "deadly weapon" in question being an illegal drug. Minnesota, Wisconsin and South Dakota allow women who continue to use substances while pregnant to be civilly committed. Some states require that medical providers report any infant who is born with a physical dependency, or who tests positive for residual traces of alcohol or drugs, to child welfare authorities.
• Cultural preferences for male children in some parts of Asia, such as Mainland China, India, South Korea, and Taiwan, have sometimes led to sex-selective abortion of female fetuses, leading to the disparity between male-to-female birth rates which is observed in some places. It is a crime in all these jurisdictions to procure an abortion for the purposes of sex selection.

Biological neuron model

From Wikipedia, the free encyclopedia

https://en.wikipedia.org/wiki/Biological_neuron_model 

 

Fig. 1. Neuron and myelinated axon, with signal flow from inputs at dendrites to outputs at axon terminals. The signal is a short electrical pulse called action potential or 'spike'.

 

Fig 2. Time course of neuronal action potential ("spike"). Note that the amplitude and the exact shape of the action potential can vary according to the exact experimental technique used for acquiring the signal.

Biological neuron models, also known as a spiking neuron models, are mathematical descriptions of the properties of certain cells in the nervous system that generate sharp electrical potentials across their cell membrane, roughly one millisecond in duration, called action potentials or spikes (Fig. 2). Since spikes are transmitted along the axon and synapses from the sending neuron to many other neurons, spiking neurons are considered to be a major information processing unit of the nervous system. Spiking neuron models can be divided into different categories: the most detailed mathematical models are biophysical neuron models (also called Hodgkin-Huxley models) that describe the membrane voltage as a function of the input current and the activation of ion channels. Mathematically simpler are integrate-and-fire models that describe the membrane voltage as a function of the input current and predict the spike times without a description of the biophysical processes that shape the time course of an action potential. Even more abstract models only predict output spikes (but not membrane voltage) as a function of the stimulation where the stimulation can occur through sensory input or pharmacologically. This article provides a short overview of different spiking neuron models and links, whenever possible to experimental phenomena. It includes deterministic and probabilistic models.

Introduction: Biological background, classification and aims of neuron models

Non-spiking cells, spiking cells, and their measurement

Not all the cells of the nervous system produce the type of spike that define the scope of the spiking neuron models. For example, cochlear hair cells, retinal receptor cells, and retinal bipolar cells do not spike. Furthermore, many cells in the nervous system are not classified as neurons but instead are classified as glia.

Neuronal activity can be measured with different experimental techniques, such as the "Whole cell" measurement technique, which captures the spiking activity of a single neuron and produces full amplitude action potentials.

With extracellular measurement techniques an electrode (or array of several electrodes) is located in the extracellular space. Spikes, often from several spiking sources, depending on the size of the electrode and its proximity to the sources, can be identified with signal processing techniques. Extracellular measurement has several advantages: 1) Is easier to obtain experimentally; 2) Is robust and lasts for a longer time; 3) Can reflect the dominant effect, especially when conducted in an anatomical region with many similar cells.

Overview of neuron models

Neuron models can be divided into two categories according to the physical units of the interface of the model. Each category could be further divided according to the abstraction/detail level:

1. Electrical input–output membrane voltage models – These models produce a prediction for membrane output voltage as a function of electrical stimulation given as current or voltage input. The various models in this category differ in the exact functional relationship between the input current and the output voltage and in the level of details. Some models in this category predict only the moment of occurrence of output spike (also known as "action potential"); other models are more detailed and account for sub-cellular processes. The models in this category can be either deterministic or probabilistic.
2. Natural stimulus or pharmacological input neuron models – The models in this category connect between the input stimulus which can be either pharmacological or natural, to the probability of a spike event. The input stage of these models is not electrical, but rather has either pharmacological (chemical) concentration units, or physical units that characterize an external stimulus such as light, sound or other forms of physical pressure. Furthermore, the output stage represents the probability of a spike event and not an electrical voltage.

Although it is not unusual in science and engineering to have several descriptive models for different abstraction/detail levels, the number of different, sometimes contradicting, biological neuron models is exceptionally high. This situation is partly the result of the many different experimental settings, and the difficulty to separate the intrinsic properties of a single neuron from measurements effects and interactions of many cells (network effects). To accelerate the convergence to a unified theory, we list several models in each category, and where applicable, also references to supporting experiments.

Aims of neuron models

Ultimately, biological neuron models aim to explain the mechanisms underlying the operation of the nervous system. However several approaches can be distinguished from more realistic models (e.g., mechanistic models) to more pragmatic models (e.g., phenomenological models). Modeling helps to analyze experimental data and address questions such as: How are the spikes of a neuron related to sensory stimulation or motor activity such as arm movements? What is the neural code used by the nervous system? Models are also important in the context of restoring lost brain functionality through neuroprosthetic devices.

Electrical input–output membrane voltage models

The models in this category describe the relationship between neuronal membrane currents at the input stage, and membrane voltage at the output stage. This category includes (generalized) integrate-and-fire models and biophysical models inspired by the work of Hodgkin–Huxley in the early 1950s using an experimental setup that punctured the cell membrane and allowed to force a specific membrane voltage/current.

Most modern electrical neural interfaces apply extra-cellular electrical stimulation to avoid membrane puncturing which can lead to cell death and tissue damage. Hence, it is not clear to what extent the electrical neuron models hold for extra-cellular stimulation.

Hodgkin–Huxley

Main article: Hodgkin–Huxley model

The Hodgkin–Huxley model (H&H model) is a model of the relationship between the flow of ionic currents across the neuronal cell membrane and the membrane voltage of the cell. It consists of a set of nonlinear differential equations describing the behaviour of ion channels that permeate the cell membrane of the squid giant axon. Hodgkin and Huxley were awarded the 1963 Nobel Prize in Physiology or Medicine for this work.

We note the voltage-current relationship, with multiple voltage-dependent currents charging the cell membrane of capacity C_m

${\displaystyle C_{\mathrm {m} }{\frac {dV(t)}{dt}}=-\sum _{i}I_{i}(t,V).}$

The above equation is the time derivative of the law of capacitance, Q = CV where the change of the total charge must be explained as the sum over the currents. Each current is given by

${\displaystyle I(t,V)=g(t,V)\cdot (V-V_{\mathrm {eq} })}$

where g(t,V) is the conductance, or inverse resistance, which can be expanded in terms of its maximal conductance ḡ and the activation and inactivation fractions m and h, respectively, that determine how many ions can flow through available membrane channels. This expansion is given by

${\displaystyle g(t,V)={\bar {g}}\cdot m(t,V)^{p}\cdot h(t,V)^{q}}$

and our fractions follow the first-order kinetics

${\displaystyle {\frac {dm(t,V)}{dt}}={\frac {m_{\infty }(V)-m(t,V)}{\tau _{\mathrm {m} }(V)}}=\alpha _{\mathrm {m} }(V)\cdot (1-m)-\beta _{\mathrm {m} }(V)\cdot m}$

with similar dynamics for h, where we can use either τ and m_∞ or α and β to define our gate fractions.

The Hodgkin–Huxley model may be extended to include additional ionic currents. Typically, these include inward Ca²⁺ and Na⁺ input currents, as well as several varieties of K⁺ outward currents, including a "leak" current.

The end result can be at the small end 20 parameters which one must estimate or measure for an accurate model. In a model of a complex systems of neurons, numerical integration of the equations are computationally expensive. Careful simplifications of the Hodgkin–Huxley model are therefore needed.

The model can be reduced to two dimensions thanks to the dynamic relations which can be established between the gating variables. it is also possible to extend it to take into account the evolution of the concentrations (considered fixed in the original model).

Perfect Integrate-and-fire

One of the earliest models of a neuron is the perfect integrate-and-fire model (also called non-leaky integrate-and-fire), first investigated in 1907 by Louis Lapicque. A neuron is represented by its membrane voltage V which evolves in time during stimulation with an input current I(t) according

${\displaystyle I(t)=C{\frac {dV(t)}{dt}}}$

which is just the time derivative of the law of capacitance, Q = CV. When an input current is applied, the membrane voltage increases with time until it reaches a constant threshold V_th, at which point a delta function spike occurs and the voltage is reset to its resting potential, after which the model continues to run. The firing frequency of the model thus increases linearly without bound as input current increases.

The model can be made more accurate by introducing a refractory period t_ref that limits the firing frequency of a neuron by preventing it from firing during that period. For constant input I(t)=I the threshold voltage is reached after an integration time t_int=CV_thr/I after start from zero. After a reset, the refractory period introduces a dead time so that the total time until the next firing is t_ref+t_int . The firing frequency is the inverse of the total inter-spike interval (including dead time). The firing frequency as a function of a constant input current is therefore

${\displaystyle \,\!f(I)={\frac {I}{C_{\mathrm {} }V_{\mathrm {th} }+t_{\mathrm {ref} }I}}.}$

A shortcoming of this model is that it describes neither adaptation nor leakage. If the model receives a below-threshold short current pulse at some time, it will retain that voltage boost forever - until another input later makes it fire. This characteristic is clearly not in line with observed neuronal behavior. The following extensions make the integrate-and-fire model more plausible from a biological point of view.

Leaky integrate-and-fire

The leaky integrate-and-fire model which can be traced back to Louis Lapicque, contains, compared to the non-leaky integrate-and-fire model a "leak" term in the membrane potential equation, reflecting the diffusion of ions through the membrane. The model equation looks like

${\displaystyle C_{\mathrm {m} }{\frac {dV_{\mathrm {m} }(t)}{dt}}=I(t)-{\frac {V_{\mathrm {m} }(t)}{R_{\mathrm {m} }}}}$

A neuron is represented by an RC circuit with a threshold. Each input pulse (e.g. caused by a spike from a different neuron) causes a short current pulse. Voltage decays exponentially. If the threshold is reached an output spike is generated and the voltage is reset.

where V_m is the voltage across the cell membrane and R_m is the membrane resistance. (The non-leaky integrate-and-fire model is retrieved in the limit R_m to infinity, i.e. if the membrane is a perfect insulator). The model equation is valid for arbitrary time-dependent input until a threshold V_th is reached; thereafter the membrane potential is reset.

For constant input, the minimum input to reach the threshold is I_th = V_th / R_m. Assuming a reset to zero, the firing frequency thus looks like

${\displaystyle f(I)={\begin{cases}0,&I\leq I_{\mathrm {th} }\\\left[t_{\mathrm {ref} }-R_{\mathrm {m} }C_{\mathrm {m} }\log \left(1-{\tfrac {V_{\mathrm {th} }}{IR_{\mathrm {m} }}}\right)\right]^{-1},&I>I_{\mathrm {th} }\end{cases}}}$

which converges for large input currents to the previous leak-free model with refractory period. The model can also be used for inhibitory neurons.

The biggest disadvantage of the Leaky integrate-and-fire neuron is that it does not contain neuronal adaptation so that it cannot describe an experimentally measured spike train in response to constant input current. This disadvantage is removed in generalized integrate-and-fire models that also contain one or several adaptation-variables and are able to predict spike times of cortical neurons under current injection to a high degree of accuracy.

Adaptive integrate-and-fire

Neuronal adaptation refers to the fact that even in the presence of a constant current injection into the soma, the intervals between output spikes increase. An adaptive integrate-and-fire neuron model combines the leaky integration of voltage V with one or several adaptation variables w_k

${\displaystyle \tau _{\mathrm {m} }{\frac {dV_{\mathrm {m} }(t)}{dt}}=RI(t)-[V_{\mathrm {m} }(t)-E_{\mathrm {m} }]-R\sum _{k}w_{k}}$

${\displaystyle \tau _{k}{\frac {dw_{k}(t)}{dt}}=-a_{k}[V_{\mathrm {m} }(t)-E_{\mathrm {m} }]-w_{k}+b_{k}\tau _{k}\sum _{f}\delta (t-t^{f})}$

where ${\displaystyle \tau _{m}}$ is the membrane time constant , w_k is the adaptation current number, with index k, ${\displaystyle \tau _{k}}$ is the time constant of adaptation current w_k, E_m is the resting potential and t^f is the firing time of the neuron and the Greek delta denotes the Dirac delta function. Whenever the voltage reaches the firing threshold the voltage is reset to a value V_r below the firing threshold. The reset value is one of the important parameters of the model. The simplest model of adaptation has only a single adaptation variable w and the sum over k is removed.

Spike times and subthreshold voltage of cortical neuron models can be predicted by generalized integrate-and-fire models such as the adaptive integrate-and-fire model, the adaptive exponential integrate-and-fire model, or the spike response model. In the example here, adaptation is implemented by a dynamic threshold which increases after each spike.

Integrate-and-fire neurons with one or several adaptation variables can account for a variety of neuronal firing patterns in response to constant stimulation, including adaptation, bursting and initial bursting. Moreover, adaptive integrate-and-fire neurons with several adaptation variables are able to predict spike times of cortical neurons under time-dependent current injection into the soma.

Fractional-order leaky integrate-and-fire

Recent advances in computational and theoretical fractional calculus lead to a new form of model, called Fractional-order leaky integrate-and-fire. An advantage of this model is that it can capture adaptation effects with a single variable. The model has the following form

${\displaystyle I(t)-{\frac {V_{\mathrm {m} }(t)}{R_{\mathrm {m} }}}=C_{\mathrm {m} }{\frac {d^{\alpha }V_{\mathrm {m} }(t)}{d^{\alpha }t}}}$

Once the voltage hits the threshold it is reset. Fractional integration has been used to account for neuronal adaptation in experimental data.

'Exponential integrate-and-fire' and 'adaptive exponential integrate-and-fire'

Main article: Exponential integrate-and-fire

In the exponential integrate-and-fire model, spike generation is exponential, following the equation:

${\displaystyle {\frac {dV}{dt}}-{\frac {R}{\tau _{m}}}I(t)={\frac {1}{\tau _{m}}}\left[E_{m}-V+\Delta _{T}\exp \left({\frac {V-V_{T}}{\Delta _{T}}}\right)\right].}$

where ${\displaystyle V}$ is the membrane potential, ${\displaystyle V_{T}}$ is the intrinsic membrane potential threshold, ${\displaystyle \tau _{m}}$ is the membrane time constant, ${\displaystyle E_{m}}$is the resting potential, and ${\displaystyle \Delta _{T}}$ is the sharpness of action potential initiation, usually around 1 mV for cortical pyramidal neurons. Once the membrane potential crosses ${\displaystyle V_{T}}$, it diverges to infinity in finite time. In numerical simulation the integration is stopped if the membrane potential hits an arbitrary threshold (much larger than ${\displaystyle V_{T}}$) at which the membrane potential is reset to a value V_r . The voltage reset value V_r is one of the important parameters of the model. Importantly, the right-hand side of the above equation contains a nonlinearity that can be directly extracted from experimental data. In this sense the exponential nonlinearity is strongly supported by experimental evidence.

In the adaptive exponential integrate-and-fire neuron  the above exponential nonlinearity of the voltage equation is combined with an adaptation variabe w

${\displaystyle \tau _{m}{\frac {dV}{dt}}=RI(t)+\left[E_{m}-V+\Delta _{T}\exp \left({\frac {V-V_{T}}{\Delta _{T}}}\right)\right]-Rw}$

${\displaystyle \tau {\frac {dw(t)}{dt}}=-a[V_{\mathrm {m} }(t)-E_{\mathrm {m} }]-w+b\tau \delta (t-t^{f})}$

Firing pattern of initial bursting in response to a step current input generated with the Adaptive exponential integrate-and-fire model. Other Firing patterns can also be generated.

where w denotes the adaptation current with time scale ${\displaystyle \tau }$. Important model parameters are the voltage reset value V_r, the intrinsic threshold ${\displaystyle V_{T}}$, the time constants ${\displaystyle \tau }$ and ${\displaystyle \tau _{m}}$ as well as the coupling parameters a and b. The adaptive exponential integrate-and-fire model inherits the experimentally derived voltage nonlinearity of the exponential integrate-and-fire model. But going beyond this model, it can also account for a variety of neuronal firing patterns in response to constant stimulation, including adaptation, bursting and initial bursting. However, since the adaptation is in the form of a current, aberrant hyperpolarization may appear. This problem was solved by expressing it as a conductance.

Stochastic models of membrane voltage and spike timing

The models in this category are generalized integrate-and-fire models that include a certain level of stochasticity. Cortical neurons in experiments are found to respond reliably to time-dependent input, albeit with a small degree of variations between one trial and the next if the same stimulus is repeated.  Stochasticity in neurons has two important sources. First, even in a very controlled experiment where input current is injected directly into the soma, ion channels open and close stochastically and this channel noise leads to a small amount of variability in the exact value of the membrane potential and the exact timing of output spikes. Second, for a neuron embedded in a cortical network, it is hard to control the exact input because most inputs come from unobserved neurons somewhere else in the brain.

Stochasticity has been introduces into spiking neuron models in two fundamentally different forms: either (i) a noisy input current is added to the differential equation of the neuron model; or (ii) the process of spike generation is noisy. In both cases, the mathematical theory can be developed for continuous time, which is then, if desired for the use in computer simulations, transformed into a discrete-time model.

The relation of noise in neuron models to variability of spike trains and neural codes is discussed in Neural Coding and in Chapter 7 of the textbook Neuronal Dynamics.

Noisy input model (diffusive noise)

A neuron embedded in a network receives spike input from other neurons. Since the spike arrival times are not controlled by an experimentalist they can be considered as stochastic. Thus a (potentially nonlinear) integrate-and-fire model with nonlinearity f(v) receives two inputs: an input ${\displaystyle I(t)}$ controlled by the experimentalists and a noisy input current ${\displaystyle I^{\rm {noise}}(t)}$ that describes the uncontrolled background input.

${\displaystyle \tau _{m}{\frac {dV}{dt}}=f(V)+RI(t)+RI^{\text{noise}}(t)}$

Stein's model is the special case of a leaky integrate-and-fire neuron and a stationary white noise current ${\displaystyle I^{\rm {noise}}(t)=\xi (t)}$ with mean zero and unit variance. In the subthreshold regime, these assumptions yield the equation of the Ornstein–Uhlenbeck process

${\displaystyle \tau _{m}{\frac {dV}{dt}}=[E_{m}-V]+RI(t)+R\xi (t)}$

However, in contrast to the standard Ornstein–Uhlenbeck process, the membrane voltage is reset whenever V hits the firing threshold V_th . Calculating the interval distribution of the Ornstein–Uhlenbeck model for constant input with threshold leads to a first-passage time problem. Stein's neuron model and variants thereof have been used to fit interspike interval distributions of spike trains from real neurons under constant input current.

In the mathematical literature, the above equation of the Ornstein–Uhlenbeck process is written in the form

${\displaystyle dV=[E_{m}-V+RI(t)]{\frac {dt}{\tau _{m}}}+\sigma \,dW}$

where ${\displaystyle \sigma }$ is the amplitude of the noise input and dW are increments of a Wiener process. For discrete-time implementations with time step dt the voltage updates are

${\displaystyle \Delta V=[E_{m}-V+RI(t)]{\frac {\Delta t}{\tau _{m}}}+\sigma {\sqrt {\tau _{m}}}y}$

where y is drawn from a Gaussian distribution with zero mean unit variance. The voltage is reset when it hits the firing threshold V_th .

The noisy input model can also be used in generalized integrate-and-fire models. For example, the exponential integrate-and-fire model with noisy input reads

${\displaystyle \tau _{m}{\frac {dV}{dt}}=E_{m}-V+\Delta _{T}\exp \left({\frac {V-V_{T}}{\Delta _{T}}}\right)+RI(t)+R\xi (t)}$

For constant deterministic input ${\displaystyle I(t)=I_{0}}$ it is possible to calculate the mean firing rate as a function of ${\displaystyle I_{0}}$. This is important because the frequency-current relation (f-I-curve) is often used by experimentalists to characterize a neuron. It is also the transfer function in

The leaky integrate-and-fire with noisy input has been widely used in the analysis of networks of spiking neurons. Noisy input is also called 'diffusive noise' because it leads to a diffusion of the subthreshold membrane potential around the noise-free trajectory (Johannesma, The theory of spiking neurons with noisy input is reviewed in Chapter 8.2 of the textbook Neuronal Dynamics.

Noisy output model (escape noise)

In deterministic integrate-and-fire models, a spike is generated if the membrane potential V(t) hits the threshold ${\displaystyle V_{th}}$. In noisy output models the strict threshold is replaced by a noisy one as follows. At each moment in time t, a spike is generated stochastically with instantaneous stochastic intensity or 'escape rate' 

${\displaystyle \rho (t)=f(V(t)-V_{th})}$

that depends on the momentary difference between the membrane voltage V(t) and the threshold ${\displaystyle V_{th}}$. A common choice for the 'escape rate' ${\displaystyle f}$ (that is consistent with biological data) is

${\displaystyle f(V-V_{th})={\frac {1}{\tau _{0}}}\exp[\beta (V-V_{th}]}$

Stochastic spike generation (noisy output) depends on the momentary difference between the membrane potential V(t) and the threshold. The membrane potential V of the spike response model (SRM) has two contributions. First, input current I is filtered by a first filter k. Second the sequence of output spikes S(t) is filtered by a second filter η and fed back. The resulting membrane V(t) potential is used to generate output spikes by a stochastic process ρ(t) with an intensity that depends on the distance between membrane potential and threshold. The spike response model (SRM) is closely related to the Generalized Linear Model (GLM).

where ${\displaystyle \tau _{0}}$is a time constant that describes how quickly a spike is fired once the membrane potential reaches the threshold and ${\displaystyle \beta }$ is a sharpness parameter. For ${\displaystyle \beta \to \infty }$ the threshold becomes sharp and spike firing occurs deterministically at the moment when the membrane potential hits the threshold from below. The sharpness value found in experiments is ${\displaystyle 1/\beta \approx 4mV}$ which means that neuronal firing becomes non-negligible as soon the membrane potential is a few mV below the formal firing threshold.

The escape rate process via a soft threshold is reviewed in Chapter 9 of the textbook Neuronal Dynamics.

For models in discrete time, a spike is generated with probability

${\displaystyle P_{F}(t_{n})=F[V(t_{n})-V_{th}]}$

that depends on the momentary difference between the membrane voltage V at time ${\displaystyle t_{n}}$ and the threshold ${\displaystyle V_{th}}$. The function F is often taken as a standard sigmoidal ${\displaystyle F(x)=0.5[1+\tanh(\gamma x)]}$ with steepness parameter ${\displaystyle \gamma }$, similar to the update dynamics in artificial neural networks. But the functional form of F can also be derived from the stochastic intensity ${\displaystyle f}$ in continuous time introduced above as ${\displaystyle F(y_{n})\approx 1-\exp[y_{n}\Delta t]}$ where ${\displaystyle y_{n}=V(t_{n})-V_{th}}$ is the distance to threshold.

Integrate-and-fire models with output noise can be used to predict the PSTH of real neurons under arbitrary time-dependent input. For non-adaptive integrate-and-fire neurons, the interval distribution under constant stimulation can be calculated from stationary renewal theory. 

Spike response model (SRM)

main article: Spike response model

The spike response model (SRM) is a general linear model for the subthreshold membrane voltage combined with a nonlinear output noise process for spike generation. The membrane voltage V(t) at time t is

${\displaystyle V(t)=\sum _{f}\eta (t-t^{f})+\int _{0}^{\infty }\kappa (s)I(t-s)\,ds+V_{\mathrm {rest} }}$

where t^f is the firing time of spike number f of the neuron, V_rest is the resting voltage in the absence of input, I(t-s) is the input current at time t-s and ${\displaystyle \kappa (s)}$ is a linear filter (also called kernel) that describes the contribution of an input current pulse at time t-s to the voltage at time t. The contributions to the voltage caused by a spike at time ${\displaystyle t^{f}}$ are described by the refractory kernel ${\displaystyle \eta (t-t^{f})}$. In particular, ${\displaystyle \eta (t-t^{f})}$ describes the reset after the spike and the time course of the spike-afterpotential following a spike. It therefore expresses the consequences of refractoriness and adaptation. The voltage V(t) can be interpreted as the result of an integration of the differential equation of a leaky integrate-and-fire model coupled to an arbitrary number of spike-triggered adaptation variables.

Spike firing is stochastic and happens with a time-dependent stochastic intensity (instantaneous rate)

${\displaystyle f(V-\vartheta (t))={\frac {1}{\tau _{0}}}\exp[\beta (V-\vartheta (t)]}$

with parameters ${\displaystyle \tau _{0}}$ and ${\displaystyle \beta }$ and a dynamic threshold ${\displaystyle \vartheta (t)}$ given by

${\displaystyle \vartheta (t)=\vartheta _{0}+\sum _{f}\theta _{1}(t-t^{f})}$

Here ${\displaystyle \vartheta _{0}}$ is the firing threshold of an inactive neuron and ${\displaystyle \theta _{1}(t-t^{f})}$ describes the increase of the threshold after a spike at time ${\displaystyle t^{f}}$. In case of a fixed threshold, one sets ${\displaystyle \theta _{1}(t-t^{f})}$=0. For ${\displaystyle \beta \to \infty }$ the threshold process is deterministic.

The time course of the filters ${\displaystyle \eta ,\kappa ,\theta _{1}}$ that characterize the spike response model can be directly extracted from experimental data. With optimized parameters the SRM describes the time course of the subthreshold membrane voltage for time-dependent input with a precision of 2mV and can predict the timing of most output spikes with a precision of 4ms. The SRM is closely related to linear-nonlinear-Poisson cascade models (also called Generalized Linear Model). The estimation of parameters of probabilistic neuron models such as the SRM using methods developed for Generalized Linear Models is discussed in Chapter 10 of the textbook Neuronal Dynamics.

Spike arrival causes postsynaptic potentials (red lines) which are summed. If the total voltage V reaches a threshold (dashed blue line) a spike is initiated (green) which also includes a spike-afterpotential. The threshold increases after each spike. Postsynaptic potentials are the response to incoming spikes while the spike-afterpotential is the response to outgoing spikes.

The name spike response model arises because in a network, the input current for neuron i is generated by the spikes of other neurons so that in the case of a network the voltage equation becomes

${\displaystyle V_{i}(t)=\sum _{f}\eta _{i}(t-t_{i}^{f})+\sum _{j=1}^{N}w_{ij}\sum _{f'}\varepsilon _{ij}(t-t_{j}^{f'})+V_{\mathrm {rest} }}$

where ${\displaystyle t_{j}^{f'}}$ are the firing times of neuron j (i.e., its spike train) , and ${\displaystyle \eta _{i}(t-t_{i}^{f})}$ describes the time course of the spike and the spike after-potential for neuron i, ${\displaystyle w_{ij}}$ and ${\displaystyle \varepsilon _{ij}(t-t_{j}^{f'})}$ describe the amplitude and time course of an excitatory or inhibitory postsynaptic potential (PSP) caused by the spike ${\displaystyle t_{j}^{f'}}$ of the presynaptic neuron j. The time course ${\displaystyle \varepsilon _{ij}(s)}$ of the PSP results from the convolution of the postsynaptic current ${\displaystyle I(t)}$ caused by the arrival of a presynaptic spike from neuron j with the membrane filter ${\displaystyle \kappa (s)}$.

SRM0

The SRM0 is a stochastic neuron model related to time-dependent nonlinear renewal theory and a simplification of the Spike Renose Model (SRM). The main difference to the voltage equation of the SRM introduced above is that in the term containing the refractory kernel ${\displaystyle \eta (s)}$ there is no summation sign over past spikes: only the most recent spike (denoted as the time ${\displaystyle {\hat {t}}}$) matters. Another difference is that the threshold is constant. The model SRM0 can be formulated in discrete or continuous time. For example, in continuous time, the single-neuron equation is

${\displaystyle V(t)=\eta (t-{\hat {t}})+\int _{0}^{\infty }\kappa (s)I(t-s)\,ds+V_{\mathrm {rest} }}$

and the network equations of the SRM0 are

${\displaystyle V_{i}(t\mid {\hat {t}}_{i})=\eta _{i}(t-{\hat {t}}_{i})+\sum _{j}w_{ij}\sum _{f}\varepsilon _{ij}(t-{\hat {t}}_{i},t-t^{f})+V_{\mathrm {rest} }}$

where ${\displaystyle {\hat {t}}_{i}}$ is the last firing time neuron i. Note that the time course of the postsynaptic potential ${\displaystyle \varepsilon _{ij}}$ is also allowed to depend on the time since the last spike of neuron i so as to describe a change in membrane conductance during refractoriness. The instantaneous firing rate (stochastic intensity) is

${\displaystyle f(V-\vartheta )={\frac {1}{\tau _{0}}}\exp[\beta (V-V_{th})]}$

where ${\displaystyle V_{th}}$ is a fixed firing threshold. Thus spike firing of neuron i depends only on its input and the time since neuron i has fired its last spike.

With the SRM0, the interspike-interval distribution for constant input can be mathematically linked to the shape of the refractory kernel ${\displaystyle \eta }$ . Moreover the stationary frequency-current relation can be calculated from the escape rate in combination with the refractory kernel ${\displaystyle \eta }$. With an appropriate choice of the kernels, the SRM0 approximates the dynamics of the Hodgkin-Huxley model to a high degree of accuracy. Moreover, the PSTH response to arbitrary time-dependent input can be predicted.

Galves–Löcherbach model

3D visualization of the Galves–Löcherbach model for biological neural nets. This visualization is set for 4,000 neurons (4 layers with one population of inhibitory neurons and one population of excitatory neurons each) at 180 intervals of time.

 

Main article: Galves–Löcherbach model

The Galves–Löcherbach model is a stochastic neuron model closely related to the spike response model SRM0  and to the leaky integrate-and-fire model. It is inherently stochastic and, just like the SRM0 linked to time-dependent nonlinear renewal theory. Given the model specifications, the probability that a given neuron ${\displaystyle i}$ spikes in a time period ${\displaystyle t}$ may be described by

${\displaystyle \mathop {\mathrm {Prob} } (X_{t}(i)=1\mid {\mathcal {F}}_{t-1})=\varphi _{i}{\Biggl (}\sum _{j\in I}W_{j\rightarrow i}\sum _{s=L_{t}^{i}}^{t-1}g_{j}(t-s)X_{s}(j),~~~t-L_{t}^{i}{\Biggl )},}$

where ${\displaystyle W_{j\rightarrow i}}$ is a synaptic weight, describing the influence of neuron ${\displaystyle j}$ on neuron ${\displaystyle i}$, ${\displaystyle g_{j}}$ expresses the leak, and ${\displaystyle L_{t}^{i}}$ provides the spiking history of neuron ${\displaystyle i}$ before ${\displaystyle t}$, according to

${\displaystyle L_{t}^{i}=\sup\{s<t:X_{s}(i)=1\}.}$

Importantly, the spike probability of neuron i depends only on its spike input (filtered with a kernel ${\displaystyle g_{j}}$ and weighted with a factor ${\displaystyle W_{j\to i}}$) and the timing of its most recent output spike (summarized by ${\displaystyle t-L_{t}^{i}}$).

Didactic toy models of membrane voltage

The models in this category are highly simplified toy models that qualitatively describe the membrane voltage as a function of input. They are mainly used for didactic reasons in teaching but are not considered valid neuron models for large-scale simulations or data fitting.

FitzHugh–Nagumo

Main article: FitzHugh–Nagumo model

Sweeping simplifications to Hodgkin–Huxley were introduced by FitzHugh and Nagumo in 1961 and 1962. Seeking to describe "regenerative self-excitation" by a nonlinear positive-feedback membrane voltage and recovery by a linear negative-feedback gate voltage, they developed the model described by

${\displaystyle {\begin{array}{rcl}{\dfrac {dV}{dt}}&=&V-V^{3}/3-w+I_{\mathrm {ext} }\\\tau {\dfrac {dw}{dt}}&=&V-a-bw\end{array}}}$

where we again have a membrane-like voltage and input current with a slower general gate voltage w and experimentally-determined parameters a = -0.7, b = 0.8, τ = 1/0.08. Although not clearly derivable from biology, the model allows for a simplified, immediately available dynamic, without being a trivial simplification. The experimental support is weak, but the model is useful as a didactic tool to introduce dynamics of spike generation through phase plane analysis. See Chapter 7 in the textbook Methods of Neuronal Modeling.

Morris–Lecar

Main article: Morris–Lecar model

In 1981 Morris and Lecar combined the Hodgkin–Huxley and FitzHugh–Nagumo models into a voltage-gated calcium channel model with a delayed-rectifier potassium channel, represented by

${\displaystyle {\begin{aligned}&C{\frac {dV}{dt}}&=&-I_{\mathrm {ion} }(V,w)+I\\[6pt]&{\frac {dw}{dt}}&=&\varphi \cdot {\frac {w_{\infty }-w}{\tau _{w}}}\end{aligned}}}$

where ${\displaystyle I_{\mathrm {ion} }(V,w)={\bar {g}}_{\mathrm {Ca} }m_{\infty }\cdot (V-V_{\mathrm {Ca} })+{\bar {g}}_{\mathrm {K} }w\cdot (V-V_{\mathrm {K} })+{\bar {g}}_{\mathrm {L} }\cdot (V-V_{\mathrm {L} })}$. The experimental support of the model is weak, but the model is useful as a didactic tool to introduce dynamics of spike generation through phase plane analysis. See Chapter 7 in the textbook Methods of Neuronal Modeling.

A two-dimensional neuron model very similar to the Morris-Lecar model can be derived step-by-step starting from the Hodgkin-Huxley model. See Chapter 4.2 in the textbook Neuronal Dynamics.

Hindmarsh–Rose

Main article: Hindmarsh–Rose model

Building upon the FitzHugh–Nagumo model, Hindmarsh and Rose proposed in 1984 a model of neuronal activity described by three coupled first-order differential equations:

${\displaystyle {\begin{aligned}&{\frac {dx}{dt}}&=&y+3x^{2}-x^{3}-z+I\\[6pt]&{\frac {dy}{dt}}&=&1-5x^{2}-y\\[6pt]&{\frac {dz}{dt}}&=&r\cdot (4(x+{\tfrac {8}{5}})-z)\end{aligned}}}$

with r² = x² + y² + z², and r ≈ 10⁻² so that the z variable only changes very slowly. This extra mathematical complexity allows a great variety of dynamic behaviors for the membrane potential, described by the x variable of the model, which include chaotic dynamics. This makes the Hindmarsh–Rose neuron model very useful, because being still simple, allows a good qualitative description of the many different firing patterns of the action potential, in particular bursting, observed in experiments. Nevertheless, it remains a toy model and has not been fitted to experimental data. It is widely used as a reference model for bursting dynamics.

Theta model and quadratic integrate-and-fire.

The theta model, or Ermentrout–Kopell canonical Type I model, is mathematically equivalent to the quadratic integrate-and-fire model which in turn is an approximation to the exponential integrate-and-fire model and the Hodgkin-Huxley model. It is called a canonical model because it is one of the generic models for constant input close to the bifurcation point, which means close to the transition from silent to repetitive firing.

The standard formulation of the theta model is

${\displaystyle {\frac {d\theta (t)}{dt}}=(I-I_{0})[1+\cos(\theta )]+[1-\cos(\theta )]}$

The equation for the quadratic integrate-and-fire model is (see Chapter 5.3 in the textbook Neuronal Dynamics)

${\displaystyle \tau _{\mathrm {m} }{\frac {dV_{\mathrm {m} }(t)}{dt}}=(I-I_{0})R+[V_{\mathrm {m} }(t)-E_{\mathrm {m} }][V_{\mathrm {m} }(t)-V_{\mathrm {T} }]}$

The equivalence of theta model and quadratic integrate-and-fire is for example reviewed in Chapter 4.1.2.2 of spiking neuron models.

For input I(t) that changes over time or is far away from the bifurcation point, it is preferable to work with the exponential integrate-and-fire model (if one wants the stay in the class of one-dimensional neuron models), because real neurons exhibit the nonlinearity of the exponential integrate-and-fire model.

Sensory input-stimulus encoding neuron models

The models in this category were derived following experiments involving natural stimulation such as light, sound, touch, or odor. In these experiments, the spike pattern resulting from each stimulus presentation varies from trial to trial, but the averaged response from several trials often converges to a clear pattern. Consequently, the models in this category generate a probabilistic relationship between the input stimulus to spike occurrences. Importantly, the recorded neurons are often located several processing steps after the sensory neurons, so that these models summarize the effects of the sequence of processing steps in a compact form

The non-homogeneous Poisson process model (Siebert)

Siebert modeled the neuron spike firing pattern using a non-homogeneous Poisson process model, following experiments involving the auditory system. According to Siebert, the probability of a spiking event at the time interval ${\displaystyle [t,t+\Delta _{t}]}$ is proportional to a non negative function ${\displaystyle g[s(t)]}$, where ${\displaystyle s(t)}$ is the raw stimulus.:

${\displaystyle P_{\text{spike}}(t\in [t',t'+\Delta _{t}])=\Delta _{t}\cdot g[s(t)]}$

Siebert considered several functions as ${\displaystyle g[s(t)]}$, including ${\displaystyle g[s(t)]\propto s^{2}(t)}$ for low stimulus intensities.

The main advantage of Siebert's model is its simplicity. The shortcomings of the model is its inability to reflect properly the following phenomena:

• The transient enhancement of the neuronal firing activity in response to a step stimulus.
• The saturation of the firing rate.
• The values of inter-spike-interval-histogram at short intervals values (close to zero).

These shortcoming are addressed by the age-dependent point process model and the two-state Markov Model.

Refractoriness and age-dependent point process model

Berry and Meister studied neuronal refractoriness using a stochastic model that predicts spikes as a product of two terms, a function f(s(t)) that depends on the time-dependent stimulus s(t) and one a recovery function ${\displaystyle w(t-{\hat {t}})}$ that depends on the time since the last spike

${\displaystyle \rho (t)=f(s(t))w(t-{\hat {t}})}$

The model is also called an inhomogeneous Markov interval (IMI) process. Similar models have been used for many years in auditory neuroscience. Since the model keeps memory of the last spike time it is non-Poisson and falls in the class of time-dependent renewal models. It is closely related to the model SRM0 with exponential escape rate. Importantly, it is possible to fit parameters of the age-dependent point process model so as to describe not just the PSTH response, but also the interspike-interval statistics.

Linear-nonlinear Poisson cascade model and GLM

Main article: Linear-nonlinear-Poisson cascade model

The linear-nonlinear-Poisson cascade model is a cascade of a linear filtering process followed by a nonlinear spike generation step. In the case that output spikes feed back, via a linear filtering process, we arrive at a model that is known in the neurosciences as Generalized Linear Model (GLM). The GLM is mathematically equivalent to the spike response model SRM) with escape noise; but whereas in the SRM the internal variables are interpreted as the membrane potential and the firing threshold, in the GLM the internal variables are abstract quantities that summarizes the net effect of input (and recent output spikes) before spikes are generated in the final step.

The two-state Markov model (Nossenson & Messer)

The spiking neuron model by Nossenson & Messer produces the probability of the neuron to fire a spike as a function of either an external or pharmacological stimulus. The model consists of a cascade of a receptor layer model and a spiking neuron model, as shown in Fig 4. The connection between the external stimulus to the spiking probability is made in two steps: First, a receptor cell model translates the raw external stimulus to neurotransmitter concentration, then, a spiking neuron model connects between neurotransmitter concentration to the firing rate (spiking probability). Thus, the spiking neuron model by itself depends on neurotransmitter concentration at the input stage.

Fig 4: High level block diagram of the receptor layer and neuron model by Nossenson & Messer.

 

Fig 5. The prediction for the firing rate in response to a pulse stimulus as given by the model by Nossenson & Messer.

An important feature of this model is the prediction for neurons firing rate pattern which captures, using a low number of free parameters, the characteristic edge emphasized response of neurons to a stimulus pulse, as shown in Fig. 5. The firing rate is identified both as a normalized probability for neural spike firing, and as a quantity proportional to the current of neurotransmitters released by the cell. The expression for the firing rate takes the following form:

${\displaystyle R_{\text{fire}}(t)={\frac {P_{\text{spike}}(t;\Delta _{t})}{\Delta _{t}}}=[y(t)+R_{0}]\cdot P_{0}(t)}$

where,

• P0 is the probability of the neuron to be "armed" and ready to fire. It is given by the following differential equation:

${\displaystyle {\dot {P}}_{0}=-[y(t)+R_{0}+R_{1}]\cdot P_{0}(t)+R_{1}}$

P0 could be generally calculated recursively using Euler method, but in the case of a pulse of stimulus it yields a simple closed form expression.

• y(t) is the input of the model and is interpreted as the neurotransmitter concentration on the cell surrounding (in most cases glutamate). For an external stimulus it can be estimated through the receptor layer model:

${\displaystyle y(t)\simeq g_{\text{gain}}\cdot \langle s^{2}(t)\rangle ,}$

with ${\displaystyle \langle s^{2}(t)\rangle }$ being short temporal average of stimulus power (given in Watt or other energy per time unit).

• R₀ corresponds to the intrinsic spontaneous firing rate of the neuron.
• R₁ is the recovery rate of the neuron from the refractory state.

Other predictions by this model include:

1) The averaged evoked response potential (ERP) due to the population of many neurons in unfiltered measurements resembles the firing rate.

2) The voltage variance of activity due to multiple neuron activity resembles the firing rate (also known as Multi-Unit-Activity power or MUA).

3) The inter-spike-interval probability distribution takes the form a gamma-distribution like function.

Experimental evidence supporting the model by Nossenson & Messer

Property of the Model by Nossenson & Messer	Description of experimental evidence
The shape of the firing rate in response to an auditory stimulus pulse	The Firing Rate has the same shape of Fig 5.
The shape of the firing rate in response to a visual stimulus pulse	The Firing Rate has the same shape of Fig 5.
The shape of the firing rate in response to an olfactory stimulus pulse	The Firing Rate has the same shape of Fig 5.
The shape of the firing rate in response to a somato-sensory stimulus	The Firing Rate has the same shape of Fig 5.
The change in firing rate in response to neurotransmitter application (mostly glutamate)	Firing Rate change in response to neurotransmitter application (Glutamate)
Square dependence between an auditory stimulus pressure and the firing rate	Square Dependence between Auditory Stimulus pressure and the Firing Rate (- Linear dependence in pressure square (power)).
Square dependence between visual stimulus electric field (volts) and the firing rate	Square dependence between visual stimulus electric field (volts) - Linear Dependence between Visual Stimulus Power and the Firing Rate.
The shape of the Inter-Spike-Interval Statistics (ISI)	ISI shape resembles the gamma-function-like
The ERP resembles the firing rate in unfiltered measurements	The shape of the averaged evoked response potential in response to stimulus resembles the firing rate (Fig. 5).
MUA power resembles the firing rate	The shape of the empirical variance of extra-cellular measurements in response to stimulus pulse resembles the firing rate (Fig. 5).

Pharmacological input stimulus neuron models

The models in this category produce predictions for experiments involving pharmacological stimulation.

Synaptic transmission (Koch & Segev)

See also: Neurotransmission

According to the model by Koch and Segev, the response of a neuron to individual neurotransmitters can be modeled as an extension of the classical Hodgkin–Huxley model with both standard and nonstandard kinetic currents. Four neurotransmitters primarily have influence in the CNS. AMPA/kainate receptors are fast excitatory mediators while NMDA receptors mediate considerably slower currents. Fast inhibitory currents go through GABA_A receptors, while GABA_B receptors mediate by secondary G-protein-activated potassium channels. This range of mediation produces the following current dynamics:

• ${\displaystyle I_{\mathrm {AMPA} }(t,V)={\bar {g}}_{\mathrm {AMPA} }\cdot [O]\cdot (V(t)-E_{\mathrm {AMPA} })}$
• ${\displaystyle I_{\mathrm {NMDA} }(t,V)={\bar {g}}_{\mathrm {NMDA} }\cdot B(V)\cdot [O]\cdot (V(t)-E_{\mathrm {NMDA} })}$
• ${\displaystyle I_{\mathrm {GABA_{A}} }(t,V)={\bar {g}}_{\mathrm {GABA_{A}} }\cdot ([O_{1}]+[O_{2}])\cdot (V(t)-E_{\mathrm {Cl} })}$
• ${\displaystyle I_{\mathrm {GABA_{B}} }(t,V)={\bar {g}}_{\mathrm {GABA_{B}} }\cdot {\tfrac {[G]^{n}}{[G]^{n}+K_{\mathrm {d} }}}\cdot (V(t)-E_{\mathrm {K} })}$

where ḡ is the maximal conductance (around 1S) and E is the equilibrium potential of the given ion or transmitter (AMDA, NMDA, Cl, or K), while [O] describes the fraction of receptors that are open. For NMDA, there is a significant effect of magnesium block that depends sigmoidally on the concentration of intracellular magnesium by B(V). For GABA_B, [G] is the concentration of the G-protein, and K_d describes the dissociation of G in binding to the potassium gates.

The dynamics of this more complicated model have been well-studied experimentally and produce important results in terms of very quick synaptic potentiation and depression, that is, fast, short-term learning.

The stochastic model by Nossenson and Messer translates neurotransmitter concentration at the input stage to the probability of releasing neurotransmitter at the output stage. For a more detailed description of this model, see the Two state Markov model section above.

HTM neuron model

The HTM neuron model was developed by Jeff Hawkins and researchers at Numenta and is based on a theory called Hierarchical Temporal Memory, originally described in the book On Intelligence. It is based on neuroscience and the physiology and interaction of pyramidal neurons in the neocortex of the human brain.

Comparing the artificial neural network (A), the biological neuron (B), and the HTM neuron (C).

 

Artificial Neural Network (ANN)	Neocortical Pyramidal Neuron (Biological Neuron)	HTM Model Neuron
- Few synapses - No dendrites - Sum input x weights - Learns by modifying weights of synapses	- Thousands of synapses on the dendrites - Active dendrites: cell recognizes hundreds of unique patterns - Co-activation of a set of synapses on a dendritic segment causes an NMDA spike and depolarization at the soma - Sources of input to the cell: Feedforward inputs which form synapses proximal to the soma and directly lead to action potentials NMDA spikes generated in the more distal basal Apical dendrites that depolarize the soma (usually not sufficient enough to generate a somatic action potential) - Learns by growing new synapses	- Inspired by the pyramidal cells in neocortex layers 2/3 and 5 - Thousands of synapses - Active dendrites: cell recognizes hundreds of unique patterns - Models dendrites and NMDA spikes with each array of coincident detectors having a set of synapses - Learns by modeling growth of new synapses

Applications

Main article: Brain–computer interface

Spiking Neuron Models are used in a variety of applications that need encoding into or decoding from neuronal spike trains in the context of neuroprosthesis and brain-computer interfaces such as retinal prosthesis: or artificial limb control and sensation. Applications are not part of this article; for more information on this topic please refer to the main article.

Relation between artificial and biological neuron models

The most basic model of a neuron consists of an input with some synaptic weight vector and an activation function or transfer function inside the neuron determining output. This is the basic structure used for artificial neurons, which in a neural network often looks like

${\displaystyle y_{i}=\varphi \left(\sum _{j}w_{ij}x_{j}\right)}$

where y_i is the output of the i th neuron, x_j is the jth input neuron signal, w_ij is the synaptic weight (or strength of connection) between the neurons i and j, and φ is the activation function. While this model has seen success in machine-learning applications, it is a poor model for real (biological) neurons, because it lacks time-dependence in input and output.

When an input is switched on at a time t and kept constant thereafter, biological neurons emit a spike train. Importantly this spike train is not regular but exhibits a temporal structure characterized by adaptation, bursting, or initial bursting followed by regular spiking. Generalized integrate-and-fire model such as the Adaptive Exponential Integrate-and-Fire model, the spike response model, or the (linear) adaptive integrate-and-fire model are able to capture these neuronal firing patterns.

Moreover, neuronal input in the brain is time-dependent. Time-dependent input is transformed by complex linear and nonlinear filters into a spike train in the output. Again, the spike response model or the adaptive integrate-and-fire model enable to predict the spike train in the output for arbitrary time-dependent input, whereas an artificial neuron or a simple leaky integrate-and-fire does not.

If we take the Hodkgin-Huxley model as a starting point, generalized integrate-and-fire models can be derived systematically in a step-by-step simplification procedure. This has been shown explicitly for the exponential integrate-and-fire model and the spike response model.

In the case of modelling a biological neuron, physical analogues are used in place of abstractions such as "weight" and "transfer function". A neuron is filled and surrounded with water containing ions, which carry electric charge. The neuron is bound by an insulating cell membrane and can maintain a concentration of charged ions on either side that determines a capacitance C_m. The firing of a neuron involves the movement of ions into the cell that occurs when neurotransmitters cause ion channels on the cell membrane to open. We describe this by a physical time-dependent current I(t). With this comes a change in voltage, or the electrical potential energy difference between the cell and its surroundings, which is observed to sometimes result in a voltage spike called an action potential which travels the length of the cell and triggers the release of further neurotransmitters. The voltage, then, is the quantity of interest and is given by V_m(t).

If the input current is constant, most neurons emit after some time of adaptation or initial bursting a regular spike train. The frequency of regular firing in response to a constant current I is described by the frequency-current relation which corresponds to the transfer function ${\displaystyle \varphi }$ of artificial neural networks. Similarly, for all spiking neuron models the transfer function ${\displaystyle \varphi }$ can be calculated numerically (or analytically).

Cable theory and compartmental models

See also: Cable theory

All of the above deterministic models are point-neuron models because they do not consider the spatial structure of a neuron. However, the dendrite contributes to transforming input into output. Point neuron models are valid description in three cases. (i) If input current is directly injected into the soma. (ii) If synaptic input arrives predominantly at or close to the soma (closeness is defined by a length scale ${\displaystyle \lambda }$ introduced below. (iii) If synapse arrive anywhere on the dendrite, but the dendrite is completely linear. In the last case the cable acts as a linear filter; these linear filter properties can be included in the formulation of generalized integrate-and-fire models such as the spike response model.

The filter properties can be calculate from a cable equation.

Let us consider a cell membrane in the form a cylindrical cable. The position on the cable is denoted by x and the voltage across the cell membrane by V. The cable is characterized by a longitudinal resistance ${\displaystyle r_{l}}$ per unit length and a membrane resistance ${\displaystyle r_{m}}$ . If everything is linear, the voltage changes as a function of time

${\displaystyle {\frac {r_{m}}{r_{l}}}{\frac {\partial ^{2}V}{\partial x^{2}}}=c_{m}r_{m}{\frac {\partial V}{\partial t}}+V}$

(19)

We introduce a length scale ${\displaystyle \lambda ^{2}={r_{m}}/{r_{l}}}$ on the left side and time constant ${\displaystyle \tau =c_{m}r_{m}}$ on the right side. The cable equation can now be written in its perhaps best known form:

${\displaystyle \lambda ^{2}{\frac {\partial ^{2}V}{\partial x^{2}}}=\tau {\frac {\partial V}{\partial t}}+V}$

(20)

The above cable equation is valid for a single cylindrical cable.

Linear cable theory describes the dendritic arbor of a neuron as a cylindrical structure undergoing a regular pattern of bifurcation, like branches in a tree. For a single cylinder or an entire tree, the static input conductance at the base (where the tree meets the cell body, or any such boundary) is defined as

${\displaystyle G_{in}={\frac {G_{\infty }\tanh(L)+G_{L}}{1+(G_{L}/G_{\infty })\tanh(L)}}}$,

where L is the electrotonic length of the cylinder which depends on its length, diameter, and resistance. A simple recursive algorithm scales linearly with the number of branches and can be used to calculate the effective conductance of the tree. This is given by

${\displaystyle \,\!G_{D}=G_{m}A_{D}\tanh(L_{D})/L_{D}}$

where A_D = πld is the total surface area of the tree of total length l, and L_D is its total electrotonic length. For an entire neuron in which the cell body conductance is G_S and the membrane conductance per unit area is G_md = G_m / A, we find the total neuron conductance G_N for n dendrite trees by adding up all tree and soma conductances, given by

${\displaystyle G_{N}=G_{S}+\sum _{j=1}^{n}A_{D_{j}}F_{dga_{j}},}$

where we can find the general correction factor F_dga experimentally by noting G_D = G_mdA_DF_dga.

The linear cable model makes a number of simplifications to give closed analytic results, namely that the dendritic arbor must branch in diminishing pairs in a fixed pattern and that dendrites are linear. A compartmental model allows for any desired tree topology with arbitrary branches and lengths, as well as arbitrary nonlinearities. It is essentially a discretized computational implementation of nonlinear dendrites.

Each individual piece, or compartment, of a dendrite is modeled by a straight cylinder of arbitrary length l and diameter d which connects with fixed resistance to any number of branching cylinders. We define the conductance ratio of the ith cylinder as B_i = G_i / G_∞, where ${\displaystyle G_{\infty }={\tfrac {\pi d^{3/2}}{2{\sqrt {R_{i}R_{m}}}}}}$ and R_i is the resistance between the current compartment and the next. We obtain a series of equations for conductance ratios in and out of a compartment by making corrections to the normal dynamic B_out,i = B_in,i+1, as

• ${\displaystyle B_{\mathrm {out} ,i}={\frac {B_{\mathrm {in} ,i+1}(d_{i+1}/d_{i})^{3/2}}{\sqrt {R_{\mathrm {m} ,i+1}/R_{\mathrm {m} ,i}}}}}$
• ${\displaystyle B_{\mathrm {in} ,i}={\frac {B_{\mathrm {out} ,i}+\tanh X_{i}}{1+B_{\mathrm {out} ,i}\tanh X_{i}}}}$
• ${\displaystyle B_{\mathrm {out,par} }={\frac {B_{\mathrm {in,dau1} }(d_{\mathrm {dau1} }/d_{\mathrm {par} })^{3/2}}{\sqrt {R_{\mathrm {m,dau1} }/R_{\mathrm {m,par} }}}}+{\frac {B_{\mathrm {in,dau2} }(d_{\mathrm {dau2} }/d_{\mathrm {par} })^{3/2}}{\sqrt {R_{\mathrm {m,dau2} }/R_{\mathrm {m,par} }}}}+\ldots }$

where the last equation deals with parents and daughters at branches, and ${\displaystyle X_{i}={\tfrac {l_{i}{\sqrt {4R_{i}}}}{\sqrt {d_{i}R_{m}}}}}$. We can iterate these equations through the tree until we get the point where the dendrites connect to the cell body (soma), where the conductance ratio is B_in,stem. Then our total neuron conductance for static input is given by

${\displaystyle G_{N}={\frac {A_{\mathrm {soma} }}{R_{\mathrm {m,soma} }}}+\sum _{j}B_{\mathrm {in,stem} ,j}G_{\infty ,j}.}$

Importantly, static input is a very special case. In biology inputs are time dependent. Moreover, dendrites are not always linear.

Compartmental models enable to include nonlinearities via ion channels positioned at arbitrary locations along the dendrites. For static inputs, it is sometimes possible to reduce the number of compartments (increase the computational speed) and yet retain the salient electrical characteristics.

See also: Multi-compartment model

Conjectures regarding the role of the neuron in the wider context of the brain principle of operation

The neurotransmitter-based energy detection scheme

The neurotransmitter-based energy detection scheme suggests that the neural tissue chemically executes a Radar-like detection procedure.

Fig. 6 The biological neural detection scheme as suggested by Nossenson et al.

As shown in Fig. 6, the key idea of the conjecture is to account neurotransmitter concentration, neurotransmitter generation and neurotransmitter removal rates as the important quantities in executing the detection task, while referring to the measured electrical potentials as a side effect that only in certain conditions coincide with the functional purpose of each step. The detection scheme is similar to a radar-like "energy detection" because it includes signal squaring, temporal summation and a threshold switch mechanism, just like the energy detector, but it also includes a unit that emphasizes stimulus edges and a variable memory length (variable memory). According to this conjecture, the physiological equivalent of the energy test statistics is neurotransmitter concentration, and the firing rate corresponds to neurotransmitter current. The advantage of this interpretation is that it leads to a unit consistent explanation which allows to bridge between electrophysiological measurements, biochemical measurements and psychophysical results.

The evidence reviewed in suggest the following association between functionality to histological classification:

1. Stimulus squaring is likely to be performed by receptor cells.
2. Stimulus edge emphasizing and signal transduction is performed by neurons.
3. Temporal accumulation of neurotransmitters is performed by glial cells. Short term neurotransmitter accumulation is likely to occur also in some types of neurons.
4. Logical switching is executed by glial cells, and it results from exceeding a threshold level of neurotransmitter concentration. This threshold crossing is also accompanied by a change in neurotransmitter leak rate.
5. Physical all-or-non movement switching is due to muscle cells and results from exceeding a certain neurotransmitter concentration threshold on muscle surroundings.

Note that although the electrophysiological signals in Fig.6 are often similar to the functional signal (signal power / neurotransmitter concentration / muscle force), there are some stages in which the electrical observation is different from the functional purpose of the corresponding step. In particular, Nossenson et al. suggested that glia threshold crossing has a completely different functional operation compared to the radiated electrophysiological signal, and that the latter might only be a side effect of glia break.

General comments regarding the modern perspective of scientific and engineering models

• The models above are still idealizations. Corrections must be made for the increased membrane surface area given by numerous dendritic spines, temperatures significantly hotter than room-temperature experimental data, and nonuniformity in the cell's internal structure. Certain observed effects do not fit into some of these models. For instance, the temperature cycling (with minimal net temperature increase) of the cell membrane during action potential propagation not compatible with models which rely on modeling the membrane as a resistance which must dissipate energy when current flows through it. The transient thickening of the cell membrane during action potential propagation is also not predicted by these models, nor is the changing capacitance and voltage spike that results from this thickening incorporated into these models. The action of some anesthetics such as inert gases is problematic for these models as well. New models, such as the soliton model attempt to explain these phenomena, but are less developed than older models and have yet to be widely applied.
• Modern views regarding of the role of the scientific model suggest that "All models are wrong but some are useful" (Box and Draper, 1987, Gribbin, 2009; Paninski et al., 2009).
• Recent conjecture suggests that each neuron might function as a collection of independent threshold units. It is suggested that a neuron could be anisotropically activated following the origin of its arriving signals to the membrane, via its dendritic trees. The spike waveform was also proposed to be dependent on the origin of the stimulus.