at present a major difficulty in understanding drug action in the human brain is the ability to relate the multiple physical scales of investigation . how do the findings of molecular biology relate to cellular neurophysiology and in turn to neuroimaging and cognitive testing ?
it is often only via inference and synthesis that investigators can relate complementary findings across research fields in order to construct a theory . to support such a synthesis ,
has sought to do this in regards to cognition and functional neuroanatomy ( i.e. by quantifying the relationship between behaviour and measures of functional brain activity ) . at the microscopic scale
there is also success in quantifying the relationship between molecular and cellular levels ( e.g. by quantifying the relationship between gene expression and electrophysiology ( hanna 2006 ) ) .
however , integrating findings from the microscopic ( single neuron ) level to the macroscopic ( neural population ) level is underdeveloped ( freeman 1992 ; breakspear and stam 2005 ) .
one potential approach is the growing area of mean field population or neural mass modelling ( freeman 1975 ; lopes da silva et  al .
1975 ; nunez 1981 ; jirsa and haken 1996 ; robinson et  al . 1997
david and friston 2003 ) . by spatially averaging the behaviour of neurons which constitute cytoarchitechtonically defined populations ( macrocolumn ) ,
what follows is an introduction to the dynamical mean field modelling of cortical drug response , using recent findings from anaesthesia as an example .
focus is placed on the physiological relevance and utility of this approach with less emphasis on the details of mathematical formulations , which are often quite involved . those readers interested in the mathematical details
are directed to the cited materials below , particularly ( steyn ross et  al .
the remainder of the current section will cover the applicability of mean field modelling to anaesthesia and in doing so begin to highlight the general principles of this approach .
mean field modelling of neural activity  will provide more details of mean field models and their biological basis . in section 
 the neurophysiological effects of anaesthetics will be introduced along with their pharmacology ( targets of anaesthetic action ) , followed by a review of network models of anaesthetic action in section  cellular network models of anaesthetic action  and a more detailed review of mean field models in section 
the clinical relevance of mean field models will be covered in section  clinical implications of mean field population models  followed by summary remarks in section 
the modelling of anaesthetics is an expedient test case of population models for several reasons .
firstly , the clinical utility of anaesthetics has made these drugs a popular research domain , with particular focus on the sites and mechanisms of action ( rudolph and antkowiak 2004 ) .
as will be discussed below this focus on the development of cellular descriptions of action greatly assist the fidelity of model parameterization .
secondly , volatile anaesthetic pharmacology shows a predominate affinity for ligand - gated gamma - amino - butyric - acid ( gaba ) neurotransmission ( franks and lieb 1994 ; rudolph and antkowiak 2004 ) .
it is therefore notable that population models incorporate fast neurotransmitter kinetics ( ligand - gated ionic channels ) and constitute cell populations not by diverse morphology , but by the chief functional properties of inhibition ( i ) and excitation ( e ) . making the dominant cortical neurotransmitter systems ( gaba ( i ) and glutamate ( e ) , petroff 2002 ) primary interests in modelling cortical pharmacology .
thirdly , as the potent influence of anaesthetics is mediated by a dominant neurotransmitter system , its resulting effects on neurophysiology and behaviour are manifestly global .
for example , neocortical populations and the related cognitive faculties are progressively suppressed by anaesthetic agents ( heinke and koelsch 2005 ) .
this global influence also fits with the spatial scale of mean field models . in summary
, anaesthetic pharmacology is fortuitously accommodated in the biological content of mean field models of cerebral cortex .
clinical benefits also arise from the development of theoretical descriptions of anaesthesia . with the growing utilization of electroencephalogram ( eeg ) based depth of anaesthesia ( doa ) monitors , the development of real - time drug monitoring in the central nervous system ( cns ) is taking shape . as a majority of mean field models aim to simulate the eeg , they are particularly germane to doa monitoring . for example , recently the effects of alprazolam , a benzodiazepine , were sensitively monitored ( off - line ) using signal processing methods derived from a population theory of eeg ( liley et  al .
furthermore , these findings provided a sensitive and mechanistic description of how the benzodiazepine class of sedative anxiolytics paradoxically increase eeg beta ( 1525  hz ) power in moderate to high concentrations ( liley et  al .
this effort at theoretically driven analysis reflects biologically motivated methods , which are hoped to provide more sensitive and efficient measures than data driven approaches .
in recent decades both theoretical neuroscience and cognitive science have shown a growing interest in the use of dynamical methods as a research tool ( kelso 1995 ; van gelder 1998 ) .
such research utilizes the tools and vocabulary of the mathematics of dynamical systems for the modelling and empirical analysis of electrical brain activity ( freeman 1992 ; stam 2005 ) .
this emerging field of neurodynamics ( freeman 2000 ) has particularly strong ties with electroencephalography ( eeg ) ( freeman 1975 ; basar 1980 ; nunez 1995 ) and magnetoencephalography ( meg ) ( jirsa and haken 1996 ; david and friston 2003 ) , as the time scales of these measurements ( milliseconds ) are more amenable to a dynamical description than bold ( blood oxygen level dependent ) fmri ( functional magnetic resonance imaging ) and pet ( positron emission tomography ) in addition to being more closely matched to the time scales of cognition and behaviour .
oscillatory dynamics of the eeg has been a main area of research focus in neurodynamics . in general terms
the theoretical methods employed to understand the dynamical genesis of the eeg can be divided into two essentially separate , though complementary , approaches .
the first approach models the activity of large networks of individual model neurons imbued with a range of voltage and ligand dependent ionic conductances ( koch and segev 1998 ) .
while this approach enables an arbitrary level of physiological detail to be included they are limited in as much as the eeg is a macroscopic or bulk property of populations of neurons .
biophysically the eeg arises from the ionic current flowing in response to synchronised synaptic activity .
a typical eeg electrode is recording the synaptic currents generated by well over 100,000 pyramidal neurons ( nunez 1981 ) ( fig .  1 ) .
thus it is simply computationally inefficient to model the eeg using this approach . a preferable approach , in which the spatial scale of the model better reflects the spatial scale of the eeg , is the continuum or mean field model ( nunez 1981 ; jirsa and haken 1996 ; robinson et  al . 1997 ; liley et  al .
it is the synaptic activity generated by spatially circumscribed populations of cortical neurons that are modelled rather than the synaptic activity of individual neurons . in general the columnar organization of cortex
( e.g. the macrocolumn ) is used to define this spatial scale and to subsequently define the spatial scales of neuronal interactions.fig . 
1schematic representation of a cortical macrocolumn and the two major subtypes of neurons in mammalian cerebral cortex  excitatory ( e ) and inhibitory ( i ) . within a macrocolumn excitatory and
inhibitory neurons are able to form all possible combinations of feed forward and feedback synaptic connectivity , whereas synaptic connections between cortical macrocolumns are exclusively excitatory . in human neocortex a cortical macrocolumn is estimated to contain of the order of 100,000300,000  neurons ( nunez 1981 ) .
electric field lines , drawn directed along the axis of the apical dendritic trees of the two excitatory neurons on the right , are the result of inward current flow due to excitatory synaptic activity in the distal portions of the apical dendritic tree , and make the greatest contribution to the surface recordable electrocorticogram ( ecog ) and eeg .
figure from liley and bojak ( 2005 ) schematic representation of a cortical macrocolumn and the two major subtypes of neurons in mammalian cerebral cortex  excitatory ( e ) and inhibitory ( i ) . within a macrocolumn excitatory and inhibitory neurons
are able to form all possible combinations of feed forward and feedback synaptic connectivity , whereas synaptic connections between cortical macrocolumns are exclusively excitatory . in human neocortex a cortical macrocolumn is estimated to contain of the order of 100,000300,000  neurons ( nunez 1981 ) .
electric field lines , drawn directed along the axis of the apical dendritic trees of the two excitatory neurons on the right , are the result of inward current flow due to excitatory synaptic activity in the distal portions of the apical dendritic tree , and make the greatest contribution to the surface recordable electrocorticogram ( ecog ) and eeg . diagram not draw to scale .
figure from liley and bojak ( 2005 ) typical of these mean field theories of electrocortical activity is that of liley et  al .
this theory of eeg is capable of producing the main features of spontaneous human eeg and is able to account quantitatively and qualitatively for the effects that gaba enhancing anaesthetics and sedatives have on the eeg .
it assumes that all the various types of cortical neurons can be divided into synaptically interacting excitatory or inhibitory subpopulations .
each of these populations has the spatial extent of a barrel - shaped region of approximately 0.53  mm in diameter penetrating the entire thickness of cortex .
locally ( i.e. within a macrocolumn ) , excitatory and inhibitory neuronal populations interact with each other and themselves .
only the excitatory populations are capable of forming connections ( long range or cortico  cortical ) with excitatory and inhibitory populations within other nearby macrocolumns .
this theory , like most other mean field theories , considers the eeg as being proportional to variations of the mean excitatory soma membrane potential he .
the soma membrane is driven from its resting potential he by a weighted sum of synaptic inputs ile from excitatory ( e ) and inhibitory ( i ) populations : 
 1\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ \tau_e\frac{\partial}{\partial t } { \bf h}_e = h_{\rm e}^{r}-{\bf h}_e+\sum_{l = e , i } \psi_{le}\left({\bf h}_e\right){\bf i}_{le}\ ; , $ $ \end{document}with a similar equation describing the inhibitory soma membrane potential hi . in the following a bold quantity f stands for mean fields depending on position on cortex and time .
we will provide here only an outline of the mathematical description , the interested reader can consult ( bojak and liley 2005 ; liley and bojak 2005 ) for necessary details .
this theory incorporates a number of key physiological properties which are parameterised in such a way that the electroencephalographic effect of changes in these properties by disease or pharmacology can be systematically studied and investigated : postsynaptic potentials and neurotransmitter kinetics : neurons are considered as passive electrical compartments into which all synaptically induced ionic current flows .
the excitatory and inhibitory synaptic currents , which give rise to excitatory ( epsp ) and inhibitory ( ipsp ) post - synaptic potentials respectively , determine deviations of the mean soma membrane potential from its resting value .
in particular the time course and the corresponding reversal potentials of fast ionotropic excitatory ( ampa / kainate ) and inhibitory ( gabaa ) neurotransmitter kinetics are explicitly included .
the ipsp response iik to presynaptic pulses pik is 
 2\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ \left[\frac{\partial}{\partial t}+\gamma_{ik}\right]\left[\frac{\partial}{\partial t}+\tilde{\gamma}_{ik}\right]{\bf i}_{ik}=\upgamma_{ik}\gamma_{ik}e^{\tilde{\gamma}_{ik}\delta_{ik}}\times{\bf p}_{ik } , $ $ \end{document}and a similar equation applies to the epsps iek .
mean field models of anaesthetic action  below discusses the connection of  and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\tilde{\gamma}$$\end{document } to the characteristic times  and  , and their particular significance to the modelling of anaesthesia.fig
.  2a schematic illustration of the essential parameters used to characterise the shape and amplitude of an inhibitory - postsynaptic potential ( ipsp ) or inhibitory postsynaptic current for both theory and experiment .
the decay constant of the ipsp , ik is defined as the value at which the ipsp has decayed to approximately 1/e of its peak value .
further details regarding this parameterisation can be found in liley and bojak ( 2005 ) .
figure from liley and bojak ( 2005)short - range intracortical connectivity : in addition to the reciprocal shortrange excitatory 
the former because they make up approximately 60% of all synaptic connections in cortex , the latter because they have been shown to be important determinants of the stability and frequency of alpha band oscillatory activity .
because these connections are short - range ( generally less than 3  mm ) conduction delays are considered negligible .
3schematic illustration of the topology of short - range ( intracortical ) , njj ) , and long - range ( cortico  cortical ) , nej , fibre systems that serve to connect excitatory ( e ) and inhibitory ( i ) neuronal populations in liley et  al.s dynamical theory of eeg generation ( liley et  al .
pjj denotes extracortical input if type j to neural population of type j. figure from liley and bojak ( 2005)long - range cortico  cortical connectivity : as already mentioned the long range cortico 
cortical fibres that serve to connect distant areas of cortex with each other are , on the basis of extensive anatomic and physiological evidence , modelled as being purely excitatory .
however because they traverse much larger distances than the intracortical axonal fibres their conduction delays can no longer be considered negligible .
the existence of conduction delays is incorporated by ascribing a finite conduction velocity to cortico  cortical connectivity .
the result is a two dimensional telegraph equation3\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ \left[\left(\frac{1}{v_{ek}}\frac{\partial}{\partial t}+\uplambda_{ek}\right)^2-\frac{3}{2}\nabla^2\right]{\bf \upphi}_{ek}=n^\alpha_{ek}\uplambda_{ek}^2{\bf s}_e , $ $ \end{document}for activity \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$${\upphi}_{ek}$$\end{document } conducted from local excitatory firing \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$${\text{s}}_e$$\end{document } to distant excitatory k  =  e orinhibitory k  =  i populations.extra-cortical input : excitatory or inhibitory thalamic input onto cortical excitatory or inhibitory is incorporated however the theory currently excludes the reciprocal thalamo - cortical connections that a number of other dynamical theories of the eeg include ( robinson et  al .
1997 ) because on the basis of current data their cortical effects can not be meaningfully parameterised .
postsynaptic potentials and neurotransmitter kinetics : neurons are considered as passive electrical compartments into which all synaptically induced ionic current flows .
the excitatory and inhibitory synaptic currents , which give rise to excitatory ( epsp ) and inhibitory ( ipsp ) post - synaptic potentials respectively , determine deviations of the mean soma membrane potential from its resting value .
in particular the time course and the corresponding reversal potentials of fast ionotropic excitatory ( ampa / kainate ) and inhibitory ( gabaa ) neurotransmitter kinetics are explicitly included .
the ipsp response iik to presynaptic pulses pik is 
 2\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ \left[\frac{\partial}{\partial t}+\gamma_{ik}\right]\left[\frac{\partial}{\partial t}+\tilde{\gamma}_{ik}\right]{\bf i}_{ik}=\upgamma_{ik}\gamma_{ik}e^{\tilde{\gamma}_{ik}\delta_{ik}}\times{\bf p}_{ik } , $ $ \end{document}and a similar equation applies to the epsps iek .
mean field models of anaesthetic action  below discusses the connection of  and \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\tilde{\gamma}$$\end{document } to the characteristic times  and  , and their particular significance to the modelling of anaesthesia.fig
.  2a schematic illustration of the essential parameters used to characterise the shape and amplitude of an inhibitory - postsynaptic potential ( ipsp ) or inhibitory postsynaptic current for both theory and experiment .
the decay constant of the ipsp , ik is defined as the value at which the ipsp has decayed to approximately 1/e of its peak value .
further details regarding this parameterisation can be found in liley and bojak ( 2005 ) .
figure from liley and bojak ( 2005 ) a schematic illustration of the essential parameters used to characterise the shape and amplitude of an inhibitory - postsynaptic potential ( ipsp ) or inhibitory postsynaptic current for both theory and experiment .
the decay constant of the ipsp , ik is defined as the value at which the ipsp has decayed to approximately 1/e of its peak value .
further details regarding this parameterisation can be found in liley and bojak ( 2005 ) .
figure from liley and bojak ( 2005 ) short - range intracortical connectivity : in addition to the reciprocal shortrange excitatory  inhibitory connections intracortical excitatory  excitatory and local inhibitory 
inhibitory connections are included  the former because they make up approximately 60% of all synaptic connections in cortex , the latter because they have been shown to be important determinants of the stability and frequency of alpha band oscillatory activity . because these connections are short - range ( generally less than 3  mm ) conduction delays are considered negligible .
3schematic illustration of the topology of short - range ( intracortical ) , njj ) , and long - range ( cortico  cortical ) , nej , fibre systems that serve to connect excitatory ( e ) and inhibitory ( i ) neuronal populations in liley et  al.s dynamical theory of eeg generation ( liley et  al . 2002 ) .
pjj denotes extracortical input if type j to neural population of type j. figure from liley and bojak ( 2005 ) schematic illustration of the topology of short - range ( intracortical ) , njj ) , and long - range ( cortico  cortical ) , nej , fibre systems that serve to connect excitatory ( e ) and inhibitory ( i ) neuronal populations in liley et  al.s dynamical theory of eeg generation ( liley et  al . 2002 ) .
pjj denotes extracortical input if type j to neural population of type j. figure from liley and bojak ( 2005 ) long - range cortico  cortical connectivity : as already mentioned the long range cortico 
cortical fibres that serve to connect distant areas of cortex with each other are , on the basis of extensive anatomic and physiological evidence , modelled as being purely excitatory .
however because they traverse much larger distances than the intracortical axonal fibres their conduction delays can no longer be considered negligible .
the existence of conduction delays is incorporated by ascribing a finite conduction velocity to cortico  cortical connectivity .
the result is a two dimensional telegraph equation3\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ \left[\left(\frac{1}{v_{ek}}\frac{\partial}{\partial t}+\uplambda_{ek}\right)^2-\frac{3}{2}\nabla^2\right]{\bf \upphi}_{ek}=n^\alpha_{ek}\uplambda_{ek}^2{\bf s}_e , $ $ \end{document}for activity \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$${\upphi}_{ek}$$\end{document } conducted from local excitatory firing \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$${\text{s}}_e$$\end{document } to distant excitatory k  =  e orinhibitory k  =  i populations .
extra - cortical input : excitatory or inhibitory thalamic input onto cortical excitatory or inhibitory is incorporated however the theory currently excludes the reciprocal thalamo - cortical connections that a number of other dynamical theories of the eeg include ( robinson et  al . 1997 ) because on the basis of current data their cortical effects can not be meaningfully parameterised . like the majority of its counterparts ( freeman 1975 ; nunez 1981 ; jirsa and haken 1996 ; robinson et  al .
1997 ) this theory is cast mathematically as a set of coupled non - linear partial differential equations .
this general mathematical form applies to a wide variety of phenomena in nature , such as the initiation and propagation of the action potential along an axon ( the celebrated hodgkin 
huxley equations ) and the flow of compressible and incompressible fluids ( the navier  stokes equations ) ( keener and sneyd 2001 ) .
because these equations are continuous in space it implies that the macroscopic dynamics of cortex are that of a spatially continuous excitable two dimensional sheet .
this theory , and other continuum mean field theories like it , are capable of producing a range of dynamical activity bearing striking similarity to that actually seen in real eeg and ecog recordings . for example
( 2002 ) is able to produce a physiologically plausible alpha rhythm when configured by physiologically admissible parameter sets , as well as producing a rich repertoire of dynamical activity that is expected to accord with the complexity of cognitive activity ( dafilis et  al .
while these mean field theories provide powerful methods for understanding a range of eeg phenomena their real physiological relevance lies in establishing a mesoscopic link between microscopic physiology and macroscopic behaviour ( fig .  4 ) .
the advantage of mean field population models of cortical tissue is the ability to explore how perturbations in a range of microscopic parameters influences neural activity macroscopically .
for example , exploring how impairing inhibitory neurotransmission microscopically may lead to hyperexcitability at the macroscopic level ( e.g. epileptic bursting in eeg ) ( kiss and erdi 2002 ; lopes da silva et  al .
4mean field theories of neural populations are best defined as mesoscopic and thus intervene between microscopic and macroscopic scales .
based on the microscopic study of neurons and their interactions , experimental values are used to parameterise inhibitory and excitatory neurons ( a ) and their connectivity at intracortical ( b ) , cortico  cortical ( c ) and extracortical ( d ) scales .
subsequently these mean values of interaction define summed physiological properties which reflect macroscopic state variables mean field theories of neural populations are best defined as mesoscopic and thus intervene between microscopic and macroscopic scales .
based on the microscopic study of neurons and their interactions , experimental values are used to parameterise inhibitory and excitatory neurons ( a ) and their connectivity at intracortical ( b ) , cortico  cortical ( c ) and extracortical ( d ) scales .
although a variety of phenomena have been explored with population models of the eeg ( e.g. epilepsy , sleep , chaos ) , only recently have researchers highlighted the dynamical modelling of drug effects ( aradi and erdi 2006 ) . in recent times , a number of researchers have developed detailed mean field descriptions of general anaesthetic drug action ( steyn - ross et  al .
these models have not only shown to be capable of simulating unexpected experimental phenomena ( e.g. the paradoxical rise and fall in eeg power that occurs during anaesthetic induction  know as the biphasic response ( kuizenga et  al .
1998 , 2001 ) ) , but have also shown the capacity to explain a number of clinical phenomenon ( e.g. anaesthetic induced seizure ) ( liley and bojak 2005 ; wilson et  al .
2006a ) . an elementary spectral comparison of the eeg between normal resting states and deep anaesthesia shows a strong suppression of alpha ( 813  hz ) and beta ( 1325  hz ) power bands , and a dominance of slow wave delta / theta ( 0.58  hz ) power ( john et  al .
additionally , during deep anaesthesia the eeg signal shows periods of spindle ( burst ) like oscillations which are interleaved with isoelectric burst suppression ( rampil 1998 ) .
these eeg phenomena follow intuitively as a psychophysiological analogue to the well know eeg characterisation of sleep stages ( hobson and pace schott 2002 ) .
however , the progression from resting state to deep anaesthesia does not appear to be a simple linear depression .
a number of researchers have noted that the eeg increases in alpha and/or beta frequencies with low doses of general anaesthetics , particularly propofol ( kishimoto et  al .
2004 ) , high doses of benzodiazepines ( liley et  al . 2003 ; jensen et  al . 2005 ) and high doses of dissociative analgesics ( rampil et  al . 1998 ; tsuda et  al .
1998 , 2001 ) , where band limited amplitude increases and decays with anaesthetic induction , and reverses with emergence .
thus , a set of clear quantitative eeg ( qeeg ) phenomena exist for empirical comparison with theoretical simulation of anaesthesia ( john et  al .
furthermore , these empirical phenomenon highlight the dependence of qeeg doa monitoring upon the capacity to detect the nuance of anaesthetic drug effect , in order to claim sensitivity .
anaesthesia is generally defined by loss of consciousness , immobility , amnesia and lack of response to noxious stimuli ( e.g. incision ) ( ishizawa 2007 ) , brought about by a diverse range of neuropharmacological action .
however , the predominant pharmacology of general anaesthetics in cerebral cortex is thought to be agonism of gabaa ionotropic receptor function ( krasowski and harrison 1999 ) .
specifically , anaesthetics varyingly potentiate the amplitude and time course of inhibitory post - synaptic potentials ( ipsp ) via positive allosteric modulation of the gabaa receptor ( hemmings et  al . 2005 ) .
this positive influence , in a simplified sense , imposes a net increase in inhibitory neurotransmission ( increased cl conductance ) and thus a depression of the cns ( rudolph and antkowiak 2004 ) . having a heteropentameric structure the gabaa receptor contains a large subunit / subtype diversity ( fig . 
recent findings in molecular biology have suggested that particular subunit compositions of the gabaa receptor mediate general anaesthesia and its sedative / hypnotic influence ( nelson et  al .
. in particular benzodiazepines may be sensitive to  and  subtypes , whilst anaesthetics differ in affinity relative to  subtype ( rudolph and antkowiak 2004).fig . 
5gaba ionotropic receptors ( gabaa and gabac ) can be formed via combinations of (16 ) , (14 ) , (13 ) , (1 ) , (1 ) , (1 ) , (1 ) and (13 ) subunits ( sieghart and sperk 2002 ) . predominately gabaa receptors are of the 212212 isoform , as shown above ( mckernan and whiting 1996 ; johnston et  al .
resulting in a cl specific ion channel that is gated by gaba , shown above as g. benzodiazepines modulate gabaa receptors via a different binding location , represented as b. anaesthetics show preference for the 2 subunit ( campagna - slater and weaver 2007 ) gaba ionotropic receptors ( gabaa and gabac ) can be formed via combinations of (16 ) , (14 ) , (13 ) , (1 ) , (1 ) , (1 ) , (1 ) and (13 ) subunits ( sieghart and sperk 2002 ) . predominately gabaa receptors are of the 212212 isoform , as shown above ( mckernan and whiting 1996 ; johnston et  al .
resulting in a cl specific ion channel that is gated by gaba , shown above as g. benzodiazepines modulate gabaa receptors via a different binding location , represented as b. anaesthetics show preference for the 2 subunit ( campagna - slater and weaver 2007 ) in modelling anaesthetic action , investigators need to pay attention to these differences in agent affinity , receptor composition and postsynaptic physiological influence . in the case of current mean field models
it is mainly the shape and amplitude of the ipsp which is parametrically altered to model anaesthetic effects ( see fig .  2 ) . due to the considerable diversity of channel composition and lack of empirical data
, individual receptor kinetics can not be comprehensively modelled at present . in the future more detailed experimental investigations
however , differences between anaesthetic agents and their gabaa receptor subtype affinity may well be accounted for sufficiently by their precise impact on the ipsps alone . as will be highlighted in section  clinical implications of mean field population models 
, minor differences in ipsp shape can lead to significant macroscopic effects . before exploring some of the details of mean field models of anaesthesia
, a brief summary of the more discrete cellular and network models will provide some of the principles necessary to parameterise drug action .
the chief aim of cellular approaches is to successfully create discrete models of spiking neurons which can be subject to analytical or numerical exploration ( arhem et  al . 2003 ) .
it is through the coupling of neurons ( nodes ) that larger cell assemblies / networks can be explored in a similar fashion .
a key example of the cellular / network approach to modelling anaesthesia is the work of traub et  al .
( 1999 ) applied their work on hippocampal gamma oscillations to anaesthetics , benzodiazepines and pro - convulsants ( thiopentone , diazapam and bicuculline respectively ) .
simulations of highly coupled inhibitory networks ( 96 all - to - all connected interneurons ) , simulated empirical findings obtained from hippocampus .
specifically , the authors found that the interneuron network peak field potential frequency ( predominantly gamma ) was sensitive to the parameter gaba(a ) , the decay time constant of hyperpolarising gabaa receptor - mediated conductance .
increases in gaba(a ) modelled benzodiazepine and anaesthetic action , consistent with postsynaptic hyperpolarisation and lowered network local field potential frequency .
conversely , the reduction of gaba(a ) reflected those agents which block or otherwise antagonise the gabaa receptor ( pro - convulsants ) , leading to hyperexcitability and increased network frequency .
one interesting development of these findings relates to capacity of high dose benzodiazepines to increase eeg beta power .
( 1999 ) speculated that increased beta activity after the ingestion of benzodiazepines may be due to the frequency of background gamma power being lowered , more synchronous and larger in amplitude thus making it more easily detectable at the scalp .
these simulation studies may prove to be helpful in understanding similar phenomena induced by low dose intravenous anaesthetics and high dose dissociative anaesthetics ( e.g. propofol and ketamine respectively ) . in their studies faulkner et  al .
explored the influence of anaesthetics upon rat hippocampal slices , both empirically and computationally . in the first of two studies faulkner et  al . explored the influence of anaesthetic and analgesic agents ( thiopental , diazepam / temazepam , morphine and ketamine ) on elicited synchronous gamma ( 40  hz ) oscillations in rat hippocampus ( ca1 )
complementary computational simulations used a neural network of individual excitatory and inhibitory layers ( each being 30   20 nodes ) .
nodes were then coupled with 50 excitatory and 49 inhibitory inputs , and a subsection ( slice ) of the array selected for simulation ( 6    20 ; 120 pyramidal ; 120 interneurons ) .
therefore from a model topology which approximated the ca1 hippocampal morphology , a subsection of this model ( slice ) could then be compared with the results obtained from experimental hippocampal slice preparations .
baseline oscillatory activity was created though tonic excitatory driving of both cell groups , with inhibitory drug effects modelled by modifying gabaa conductance .
. showed strong similarities between slice recordings and computational simulations of thiopental s suppressive influence on gamma synchrony .
the authors postulated that decreased gamma synchrony lowered distributed microscopic communication in the hippocampus , providing a putative amnestic mechanism .
extrapolations of this mechanism therefore suggests that anaesthesia results from wide spread loss of high frequency synchronized network activity .
are sensitive to changes in the specific subunit composition of the gabaa receptor . in a recent investigation ujfalussy et  al .
( 2007 ) combined experimental and computational methods to explore gabaa receptor  subunit sensitivity of two positive allosteric modulators ( zolpidem and l-838417 ) .
computational models were constructed to simulate temporal pacemaker properties of the septo - hippocampal system .
( 2007 ) assumed that the maximal synaptic conductance of all gabaa ionotropic receptor channels in hippocampus ( ca1 ) and medial septum would increase .
therefore , either the inhibitory time ( decay ) constant or amplitude of the ipsp was modified to increase the net charge transfer ( synaptic current ) .
in contrast , the other agent modelled , l-838417 , had specific influence on hippocampal pyramidal neurons and was therefore modelled by lowering the driving phasic input of septo - hippocampal projecting neurons .
results showed that the simulated effect of zolpidem and l-838417 produced differing patterns of spike and oscillatory behaviour , however both reduced the commonly observed theta oscillation of the septo - hippocampal circuit ( ujfalussy et  al .
as zolpidem binds to 1gabaa receptors and l-838417 to 2/3/5 gabaa receptors , these methods can assist in drug development . in this particular case
, the authors suggest that some of the side effects of the broadly binding benzodiazepines ( e.g. diazepam ) may be avoided by compounds with more specific anxiolytic gabaergic targets ( ujfalussy et  al .
this study highlighted the capacity to not only simulate drug action by specific pharmacological parameterisation , but also through altering specific cell group connectivity ( i.e. the tonic driving strength of excitatory pyramidal neurons ) .
over the past 510  years a small number of research groups have focused on developing mean field population models of anaesthetic action , which sceniak and maciver ( 2006 ) referred to as anaesthesia in  silico in a recent editorial .
these mesoscopic models are complimentary to single neuronal network based methods , as they allow collective properties of neural tissue dynamics to be simulated on much larger scales .
thus it can be said that while discrete approaches to macroscopic cortical modelling are impractical , so too are the application of population field models to the modelling of microscopic circuits .
the first reported exploration of anaesthetic action in a mean field population model was that of steyn - ross et  al .
( 1998 ) to a set of equations which facilitated analytical exploration using a range of established stochastic methods . with this new set of equations steyn - ross et  al .
( 1999 ) aimed to simulate the effects of anaesthesia by systematically modulating parameters that correspond to the known pharmacological targets associated with anaesthetic action . in particular , of the many system parameters steyn - ross et  al .
( 1999 ) specifically modulated the inhibitory neurotransmitter rate constant ( i ) , which quantifies the shape of the ipsp . as highlighted above
, general anaesthetics typically exert their influence by positive allosteric modulation of the gabaa receptor , increasing negative charge transfer and prolonging the tail of the unitary ipsp . thus reducing i , which increases the time scale of inhibitory neurotransmission , models the increasing effects of anaesthetics
. therefore , anaesthesia was parameterised by scaling the inhibitory neurotransmitter rate constant relative to empirical findings of ipsp prolongation by anaesthetics ( franks and lieb 1994 ) . by observing how the state variables equivalent to the eeg ( excitatory mean soma membrane potential ; he ) changed with modelled increases in drug concentration , the effect of anaesthesia on the eeg was simulated .
( 1999 ) found that through the modulation of the inhibitory neurotransmitter rate constant , the model showed three stages of analytical equilibrium which they termed : coma ( i ) , wakefulness ( ii ) and epilepsy ( iii ) ( fig . 
6d ) . each of these three system states reflected a progressive reduction of anaesthetic concentration .
high levels of modelled anaesthesia ( large prolonging of ipsp ) produced a state neurophysiologically equivalent to coma or very deep anaesthesia .
subsequently , approximate zero levels of simulated anaesthetic effect produced wakefulness , which in turn could be perturbed to an epileptic like state when anaesthetic influence became negative ( equivalent to simulating a pro - convulsant compound which antagonises the gabaa ipsp , e.g. bicuculline ) .
( 1999 ) also were able to model a number of important empirical eeg phenomena associated with anaesthetic action , including the biphasic effect ( fig . 
as indicated above when comparing normal resting levels with deep anaesthesia the eeg shows strong slow wave ( delta / theta ) dominance .
however , the transition between these two states shows an unexpected rise and fall of band limited ( 1115  hz ) amplitude , hence the term biphasic response ( kuizenga et  al .
( 1999 ) found that modelled eeg spectra also showed a clear biphasic response during simulated anaesthesia.fig . 
( 2004 ) modelled the abrupt changes in eeg spectral power that have been commonly observed to occur during anaesthetic induction and emergence .
( a ) the abrupt biphasic response in total modelled eeg spectral power during simulated anaesthetic induction .
anaesthetic action was modelled as a prolongation of the duration of the gabaergic ipsp (  ) .
( c ) actual eeg recordings of anaesthetic induction and emergence with the intravenous anaesthetic agent propofol which illustrate clearly the biphasic variations in total eeg power as anaesthetic concentration is increased and decreased .
propose that anaesthetic induced variations in eeg spectral power can be explained theoretically by a stochastically driven system switching between two stable states via a saddle node bifurcation caused by anaesthetic induced variations in ipsp duration .
( 2004 ) modelled the abrupt changes in eeg spectral power that have been commonly observed to occur during anaesthetic induction and emergence .
( a ) the abrupt biphasic response in total modelled eeg spectral power during simulated anaesthetic induction .
anaesthetic action was modelled as a prolongation of the duration of the gabaergic ipsp (  ) .
( c ) actual eeg recordings of anaesthetic induction and emergence with the intravenous anaesthetic agent propofol which illustrate clearly the biphasic variations in total eeg power as anaesthetic concentration is increased and decreased .
propose that anaesthetic induced variations in eeg spectral power can be explained theoretically by a stochastically driven system switching between two stable states via a saddle node bifurcation caused by anaesthetic induced variations in ipsp duration .
2004 ) , with permission the primary conclusion of steyn - ross et  al .
( 1999 ) was that the empirical biphasic response could be understood as a parametrically induced phase transition between two dynamically distinct states .
( 1999 ) went on to further develop their phase transition description of anaesthesia , drawing analogy to the phase transitions commonly seen in a range of thermodynamic systems , such as the sudden transition of a liquid to a gas during heating ( steyn - ross et  al .
( 2001b ) developed an excitation parameter which reflected the inverse of anaesthetic effect , and therefore constituted the analogue of heat .
it was through a thorough analytical analysis of changes in excitation ( heat ) that steyn - ross et  al .
the authors justify the physiological relevance of decreases in entropy after the phase transition to unconsciousness ( low excitation and strong anaesthetic effect ) , by pointing to empirical findings of lowered metabolic state ( stullken et  al . 1977 ) and lowered eeg spectral entropy during anaesthesia .
( 1999 ) had used a simplified and spatially lumped version of the liley et  al .
( 1998 ) mean field continuum theory , some of the qualitative characteristics of eeg dynamics were not preserved .
therefore in order to begin exploration of the spatial character of the dynamic changes induced by anaesthesia , steyn - ross et  al .
( 2003 ) constructed a one - dimensional cortex of linked neural masses ( cortical rod ) .
this theoretical construction allowed for rudimentary comparison of spatial inhomogeneities during anaesthetic induction ( i.e. the theoretical comparison of two separated monopolar electrodes , with a common reference ) ( steyn - ross et  al .
( 2003 ) related the findings of their analytical work to recent discoveries in anaesthesia research , highlighting their theoretical findings of increased spatial covariance as being supported by those empirical results showing increased inter - electrode coherence at the points of loss and recovery of consciousness ( i.e. the points of proposed phase transition ) ( john et  al .
( 2004 ) posited that the transition from consciousness to unconsciousness occurs as a consequence of a phase transition in cortical neural population dynamics brought about due to the increasing inhibitory effects of modelled anaesthetic concentration . during this transition
there is a surge in band limited eeg power which is of greater magnitude during emergence than during induction , with this hysteresis being observed empirically .
furthermore , the phase transition is also characterised by increased spatial and temporal correlation , for both induction and emergence , with cortical entropy being lowered during unconsciousness . on this basis
the authors argue that anaesthetics induce a sharp rise in eeg power that produces a rapid shift into a state of slowed and more ordered eeg rhythmicity , which is reversed , though not symmetrically , with emergence .
while the results of steyn - ross et  al . were able to account for a number of important macroscopic eeg properties associated with anaesthetic action , they did so in a simplified model that was unable to produce normal resting oscillatory behaviour , such as alpha . in order to remedy this bojak et  al .
( 2002 ) equations to model the effects of anaesthetics on the eeg . unlike the original investigations of steyn - ross et  al .
( 2002 ) equations , rather the full set of equations where numerically solved on a spatio - temporal scale that would more accurately approximate observed eeg . in doing so
, the influence of microscopically defined parameter modification could be observed on a macroscopic scale that is comparable to empirical eeg .
( 2004 ) were also able to qualitatively reproduce the biphasic response through only the modulation of the inhibitory decay constant ( ipsp prolongation ) .
a further development in bojak and liley ( 2005 ) was the reparametrisation of postsynaptic potentials to more accurately fit experimental data .
alpha form had been used , i.e. , for an incoming presynaptic pulse at time t  =  0 the assumed psp response at t    0 was:4\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ r_\alpha(t)=\upgamma\gamma\exp(1)\times t \exp\left(-\gamma t\right ) , $ $ \end{document } as illustrated in fig .  2 for an ipsp
, this response has a maximum amplitude  at rise time 1/. unfortunately the psp decay time  for which the response falls again to /exp(1 ) is for the 
alpha form fixed at   3.1462    . anaesthetic agents often potentiate gabaa induced currents and/or attenuate glutamate mediated ionic currents ( rudolph and antkowiak 2004 ) , thereby selectively prolonging the inhibitory ik without affecting the inhibitory rise time ik . hence bojak and liley ( 2005 ) introduced eq .  2 yielding a response of the form5\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 

\usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ r(t)=\upgamma\gamma \exp\left(\tilde{\gamma}\delta\right)\times \frac{\exp\left(-\gamma t\right)-\exp\left(-\tilde{\gamma}t\right)}{\tilde{\gamma}-\gamma } , $ $ \end{document}with \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\delta\equiv\hbox{ln}(\tilde{\gamma}/\gamma)/(\tilde{\gamma}-\gamma).$$\end{document } this bi - exponential has the same maximum amplitude r(t  =   )  =  . however , the additional parameter \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\tilde{\gamma}$$\end{document } can be used to adjust the decay time freely .
this allowed bojak and liley ( 2005 ) to directly incorporate into their model experimental data on the effect of isoflurane on psps ( maciver et  al . 1996
; banks and pearce 1999 ; nishikawa and maciver 2000 ) , and to predict the eeg under anaesthesia by calculating the corresponding changes to the excitatory mean soma membrane potential he , see eq .  1 .
bojak and liley ( 2005 ) were able to show systematic changes in modelled eeg as a function of isoflurane concentration .
an additional selection criteria was then imposed upon these sets of data , which aimed at extracting sets producing a clear biphasic response ( bojak and liley 2005 ) .
the outcomes of bojak and liley s ( 2005 ) explorations showed some difference to that of steyn - ross et  al .
it is particularly worth highlighting that while both theories describe changes in macroscopic dynamics during anaesthesia , steyn - ross et  al .
( 2004 ) suggest a phase transition ( e.g. boiling water to steam ) , whereas bojak and liley ( 2005 ) posit a continuum transition .
bojak and liley ( 2005 ) highlight that the sharp rise in neural firing rate required for steyn - ross et  al.s ( 2004 ) phase transition is not supported by empirical findings of smooth suppression of action potential frequency with progressive anaesthesia .
( 2007 ) sought to further develop the work of steyn - ross et  al .
( 2004 ) and bojak and liley ( 2005 ) , by implementing a model parameterisation of slow voltage dependent mechanisms on excitatory firing rate dynamics .
the common sigmoidal relationship between mean soma membrane potential and mean firing rate was modified to included the influence of slow shifts in membrane potential .
this voltage dependant slow ionic leak , was motivated by the empirically observed switching ( i.e. up and down states ) in firing rates of neurons during anaesthesia and non - rem sleep , thought to be controlled by cation conductance .
( 2007 ) explored progressive increases in simulated desflurane anaesthesia , showing up state dominance during waking levels , up / down state switching ( delta ) at moderate anaesthesia and down state predominance during strong anaesthetic effect ( fig . 
this switching in firing rate between bursting up states and quiescent down states relates biophysically to the synchronised delta power local field potential predominately seen during deep anaesthesia and stage 3/4 sleep .
these findings further promote the ability to theoretically relate cellular experimental observations with bulk macroscopic cortical dynamics via mean field theoretics.fig . 
( 2007 ) simulated eeg ( mean soma membrane potential he ) at various concentrations of simulated anaesthesia by assuming that anaesthetics , in addition to increasing the duration of the ipsp , also cause slow adaptive changes in the shape of the mesoscopic firing rate function which relates mean soma membrane potential to mean population firing rate . ( 2 ) actual eeg recorded in children undergoing anaesthesia with the volatile halogenated agent desflurane .
( a ) 0 mac ( minimum alveolar concentration ) , ( b ) 1 mac , ( c ) >  1 mac and ( d ) 2 mac .
both real and simulated data show that increasing anaesthetic effect is associated with enhanced total eeg power and reductions in mean frequency .
( 2007 ) , with permission ( 1 ) molaee - ardekani et  al .
( 2007 ) simulated eeg ( mean soma membrane potential he ) at various concentrations of simulated anaesthesia by assuming that anaesthetics , in addition to increasing the duration of the ipsp , also cause slow adaptive changes in the shape of the mesoscopic firing rate function which relates mean soma membrane potential to mean population firing rate .
( 2 ) actual eeg recorded in children undergoing anaesthesia with the volatile halogenated agent desflurane . ( a )
0 mac ( minimum alveolar concentration ) , ( b ) 1 mac , ( c ) >  1 mac and ( d ) 2 mac .
both real and simulated data show that increasing anaesthetic effect is associated with enhanced total eeg power and reductions in mean frequency .
anaesthesia is generally defined by loss of consciousness , immobility , amnesia and lack of response to noxious stimuli ( e.g. incision ) ( ishizawa 2007 ) , brought about by a diverse range of neuropharmacological action .
however , the predominant pharmacology of general anaesthetics in cerebral cortex is thought to be agonism of gabaa ionotropic receptor function ( krasowski and harrison 1999 ) .
specifically , anaesthetics varyingly potentiate the amplitude and time course of inhibitory post - synaptic potentials ( ipsp ) via positive allosteric modulation of the gabaa receptor ( hemmings et  al . 2005 ) .
this positive influence , in a simplified sense , imposes a net increase in inhibitory neurotransmission ( increased cl conductance ) and thus a depression of the cns ( rudolph and antkowiak 2004 ) . having a heteropentameric structure the gabaa receptor contains a large subunit / subtype diversity ( fig . 
recent findings in molecular biology have suggested that particular subunit compositions of the gabaa receptor mediate general anaesthesia and its sedative / hypnotic influence ( nelson et  al .
in particular benzodiazepines may be sensitive to  and  subtypes , whilst anaesthetics differ in affinity relative to  subtype ( rudolph and antkowiak 2004).fig . 
5gaba ionotropic receptors ( gabaa and gabac ) can be formed via combinations of (16 ) , (14 ) , (13 ) , (1 ) , (1 ) , (1 ) , (1 ) and (13 ) subunits ( sieghart and sperk 2002 ) . predominately gabaa receptors are of the 212212 isoform , as shown above ( mckernan and whiting 1996 ; johnston et  al .
resulting in a cl specific ion channel that is gated by gaba , shown above as g. benzodiazepines modulate gabaa receptors via a different binding location , represented as b. anaesthetics show preference for the 2 subunit ( campagna - slater and weaver 2007 ) gaba ionotropic receptors ( gabaa and gabac ) can be formed via combinations of (16 ) , (14 ) , (13 ) , (1 ) , (1 ) , (1 ) , (1 ) and (13 ) subunits ( sieghart and sperk 2002 ) . predominately gabaa receptors are of the 212212 isoform , as shown above ( mckernan and whiting 1996 ; johnston et  al .
resulting in a cl specific ion channel that is gated by gaba , shown above as g. benzodiazepines modulate gabaa receptors via a different binding location , represented as b. anaesthetics show preference for the 2 subunit ( campagna - slater and weaver 2007 ) in modelling anaesthetic action , investigators need to pay attention to these differences in agent affinity , receptor composition and postsynaptic physiological influence . in the case of current mean field models
it is mainly the shape and amplitude of the ipsp which is parametrically altered to model anaesthetic effects ( see fig .  2 ) . due to the considerable diversity of channel composition and lack of empirical data
, individual receptor kinetics can not be comprehensively modelled at present . in the future more detailed experimental investigations
however , differences between anaesthetic agents and their gabaa receptor subtype affinity may well be accounted for sufficiently by their precise impact on the ipsps alone . as will be highlighted in section  clinical implications of mean field population models 
before exploring some of the details of mean field models of anaesthesia , a brief summary of the more discrete cellular and network models will provide some of the principles necessary to parameterise drug action .
the chief aim of cellular approaches is to successfully create discrete models of spiking neurons which can be subject to analytical or numerical exploration ( arhem et  al .
it is through the coupling of neurons ( nodes ) that larger cell assemblies / networks can be explored in a similar fashion .
a key example of the cellular / network approach to modelling anaesthesia is the work of traub et  al .
( 1999 ) applied their work on hippocampal gamma oscillations to anaesthetics , benzodiazepines and pro - convulsants ( thiopentone , diazapam and bicuculline respectively ) .
simulations of highly coupled inhibitory networks ( 96 all - to - all connected interneurons ) , simulated empirical findings obtained from hippocampus .
specifically , the authors found that the interneuron network peak field potential frequency ( predominantly gamma ) was sensitive to the parameter gaba(a ) , the decay time constant of hyperpolarising gabaa receptor - mediated conductance .
increases in gaba(a ) modelled benzodiazepine and anaesthetic action , consistent with postsynaptic hyperpolarisation and lowered network local field potential frequency .
conversely , the reduction of gaba(a ) reflected those agents which block or otherwise antagonise the gabaa receptor ( pro - convulsants ) , leading to hyperexcitability and increased network frequency .
one interesting development of these findings relates to capacity of high dose benzodiazepines to increase eeg beta power .
( 1999 ) speculated that increased beta activity after the ingestion of benzodiazepines may be due to the frequency of background gamma power being lowered , more synchronous and larger in amplitude thus making it more easily detectable at the scalp .
these simulation studies may prove to be helpful in understanding similar phenomena induced by low dose intravenous anaesthetics and high dose dissociative anaesthetics ( e.g. propofol and ketamine respectively ) . in their studies faulkner et  al . explored the influence of anaesthetics upon rat hippocampal slices , both empirically and computationally . in the first of two studies faulkner et  al . explored the influence of anaesthetic and analgesic agents ( thiopental , diazepam / temazepam , morphine and ketamine ) on elicited synchronous gamma ( 40  hz ) oscillations in rat hippocampus ( ca1 ) .
complementary computational simulations used a neural network of individual excitatory and inhibitory layers ( each being 30   20 nodes ) .
nodes were then coupled with 50 excitatory and 49 inhibitory inputs , and a subsection ( slice ) of the array selected for simulation ( 6    20 ; 120 pyramidal ; 120 interneurons ) .
therefore from a model topology which approximated the ca1 hippocampal morphology , a subsection of this model ( slice ) could then be compared with the results obtained from experimental hippocampal slice preparations .
baseline oscillatory activity was created though tonic excitatory driving of both cell groups , with inhibitory drug effects modelled by modifying gabaa conductance .
. showed strong similarities between slice recordings and computational simulations of thiopental s suppressive influence on gamma synchrony .
the authors postulated that decreased gamma synchrony lowered distributed microscopic communication in the hippocampus , providing a putative amnestic mechanism .
extrapolations of this mechanism therefore suggests that anaesthesia results from wide spread loss of high frequency synchronized network activity .
are sensitive to changes in the specific subunit composition of the gabaa receptor . in a recent investigation
( 2007 ) combined experimental and computational methods to explore gabaa receptor  subunit sensitivity of two positive allosteric modulators ( zolpidem and l-838417 ) .
computational models were constructed to simulate temporal pacemaker properties of the septo - hippocampal system .
( 2007 ) assumed that the maximal synaptic conductance of all gabaa ionotropic receptor channels in hippocampus ( ca1 ) and medial septum would increase .
therefore , either the inhibitory time ( decay ) constant or amplitude of the ipsp was modified to increase the net charge transfer ( synaptic current ) .
in contrast , the other agent modelled , l-838417 , had specific influence on hippocampal pyramidal neurons and was therefore modelled by lowering the driving phasic input of septo - hippocampal projecting neurons .
results showed that the simulated effect of zolpidem and l-838417 produced differing patterns of spike and oscillatory behaviour , however both reduced the commonly observed theta oscillation of the septo - hippocampal circuit ( ujfalussy et  al .
as zolpidem binds to 1gabaa receptors and l-838417 to 2/3/5 gabaa receptors , these methods can assist in drug development . in this particular case ,
the authors suggest that some of the side effects of the broadly binding benzodiazepines ( e.g. diazepam ) may be avoided by compounds with more specific anxiolytic gabaergic targets ( ujfalussy et  al .
this study highlighted the capacity to not only simulate drug action by specific pharmacological parameterisation , but also through altering specific cell group connectivity ( i.e. the tonic driving strength of excitatory pyramidal neurons ) .
over the past 510  years a small number of research groups have focused on developing mean field population models of anaesthetic action , which sceniak and maciver ( 2006 ) referred to as anaesthesia in  silico in a recent editorial .
these mesoscopic models are complimentary to single neuronal network based methods , as they allow collective properties of neural tissue dynamics to be simulated on much larger scales .
thus it can be said that while discrete approaches to macroscopic cortical modelling are impractical , so too are the application of population field models to the modelling of microscopic circuits .
the first reported exploration of anaesthetic action in a mean field population model was that of steyn - ross et  al .
( 1998 ) to a set of equations which facilitated analytical exploration using a range of established stochastic methods . with this new set of equations steyn - ross et  al .
( 1999 ) aimed to simulate the effects of anaesthesia by systematically modulating parameters that correspond to the known pharmacological targets associated with anaesthetic action .
( 1999 ) specifically modulated the inhibitory neurotransmitter rate constant ( i ) , which quantifies the shape of the ipsp . as highlighted above
, general anaesthetics typically exert their influence by positive allosteric modulation of the gabaa receptor , increasing negative charge transfer and prolonging the tail of the unitary ipsp . thus reducing i , which increases the time scale of inhibitory neurotransmission , models the increasing effects of anaesthetics
. therefore , anaesthesia was parameterised by scaling the inhibitory neurotransmitter rate constant relative to empirical findings of ipsp prolongation by anaesthetics ( franks and lieb 1994 ) . by observing how the state variables equivalent to the eeg ( excitatory mean soma membrane potential ; he ) changed with modelled increases in drug concentration , the effect of anaesthesia on the eeg was simulated .
( 1999 ) found that through the modulation of the inhibitory neurotransmitter rate constant , the model showed three stages of analytical equilibrium which they termed : coma ( i ) , wakefulness ( ii ) and epilepsy ( iii ) ( fig . 
6d ) . each of these three system states reflected a progressive reduction of anaesthetic concentration .
high levels of modelled anaesthesia ( large prolonging of ipsp ) produced a state neurophysiologically equivalent to coma or very deep anaesthesia .
subsequently , approximate zero levels of simulated anaesthetic effect produced wakefulness , which in turn could be perturbed to an epileptic like state when anaesthetic influence became negative ( equivalent to simulating a pro - convulsant compound which antagonises the gabaa ipsp , e.g. bicuculline ) .
( 1999 ) also were able to model a number of important empirical eeg phenomena associated with anaesthetic action , including the biphasic effect ( fig . 
as indicated above when comparing normal resting levels with deep anaesthesia the eeg shows strong slow wave ( delta / theta ) dominance .
however , the transition between these two states shows an unexpected rise and fall of band limited ( 1115  hz ) amplitude , hence the term biphasic response ( kuizenga et  al .
( 1999 ) found that modelled eeg spectra also showed a clear biphasic response during simulated anaesthesia.fig . 
( 2004 ) modelled the abrupt changes in eeg spectral power that have been commonly observed to occur during anaesthetic induction and emergence .
( a ) the abrupt biphasic response in total modelled eeg spectral power during simulated anaesthetic induction .
anaesthetic action was modelled as a prolongation of the duration of the gabaergic ipsp (  ) .
( c ) actual eeg recordings of anaesthetic induction and emergence with the intravenous anaesthetic agent propofol which illustrate clearly the biphasic variations in total eeg power as anaesthetic concentration is increased and decreased .
propose that anaesthetic induced variations in eeg spectral power can be explained theoretically by a stochastically driven system switching between two stable states via a saddle node bifurcation caused by anaesthetic induced variations in ipsp duration .
( 2004 ) modelled the abrupt changes in eeg spectral power that have been commonly observed to occur during anaesthetic induction and emergence .
( a ) the abrupt biphasic response in total modelled eeg spectral power during simulated anaesthetic induction .
anaesthetic action was modelled as a prolongation of the duration of the gabaergic ipsp (  ) .
( c ) actual eeg recordings of anaesthetic induction and emergence with the intravenous anaesthetic agent propofol which illustrate clearly the biphasic variations in total eeg power as anaesthetic concentration is increased and decreased .
propose that anaesthetic induced variations in eeg spectral power can be explained theoretically by a stochastically driven system switching between two stable states via a saddle node bifurcation caused by anaesthetic induced variations in ipsp duration .
2004 ) , with permission the primary conclusion of steyn - ross et  al .
( 1999 ) was that the empirical biphasic response could be understood as a parametrically induced phase transition between two dynamically distinct states .
( 1999 ) went on to further develop their phase transition description of anaesthesia , drawing analogy to the phase transitions commonly seen in a range of thermodynamic systems , such as the sudden transition of a liquid to a gas during heating ( steyn - ross et  al .
( 2001b ) developed an excitation parameter which reflected the inverse of anaesthetic effect , and therefore constituted the analogue of heat .
it was through a thorough analytical analysis of changes in excitation ( heat ) that steyn - ross et  al .
( 2001b ) constructed a dynamical mean field description of anaesthesia . as with thermodynamic phase transitions , steyn - ross et  al .
the authors justify the physiological relevance of decreases in entropy after the phase transition to unconsciousness ( low excitation and strong anaesthetic effect ) , by pointing to empirical findings of lowered metabolic state ( stullken et  al .
( 1999 ) had used a simplified and spatially lumped version of the liley et  al .
( 1998 ) mean field continuum theory , some of the qualitative characteristics of eeg dynamics were not preserved .
therefore in order to begin exploration of the spatial character of the dynamic changes induced by anaesthesia , steyn - ross et  al .
( 2003 ) constructed a one - dimensional cortex of linked neural masses ( cortical rod ) .
this theoretical construction allowed for rudimentary comparison of spatial inhomogeneities during anaesthetic induction ( i.e. the theoretical comparison of two separated monopolar electrodes , with a common reference ) ( steyn - ross et  al .
( 2003 ) related the findings of their analytical work to recent discoveries in anaesthesia research , highlighting their theoretical findings of increased spatial covariance as being supported by those empirical results showing increased inter - electrode coherence at the points of loss and recovery of consciousness ( i.e. the points of proposed phase transition ) ( john et  al .
( 2004 ) posited that the transition from consciousness to unconsciousness occurs as a consequence of a phase transition in cortical neural population dynamics brought about due to the increasing inhibitory effects of modelled anaesthetic concentration . during this transition
there is a surge in band limited eeg power which is of greater magnitude during emergence than during induction , with this hysteresis being observed empirically .
furthermore , the phase transition is also characterised by increased spatial and temporal correlation , for both induction and emergence , with cortical entropy being lowered during unconsciousness . on this basis
the authors argue that anaesthetics induce a sharp rise in eeg power that produces a rapid shift into a state of slowed and more ordered eeg rhythmicity , which is reversed , though not symmetrically , with emergence .
while the results of steyn - ross et  al . were able to account for a number of important macroscopic eeg properties associated with anaesthetic action , they did so in a simplified model that was unable to produce normal resting oscillatory behaviour , such as alpha . in order to remedy this bojak et  al .
( 2002 ) equations to model the effects of anaesthetics on the eeg . unlike the original investigations of steyn - ross et  al .
( 2002 ) equations , rather the full set of equations where numerically solved on a spatio - temporal scale that would more accurately approximate observed eeg . in doing so
, the influence of microscopically defined parameter modification could be observed on a macroscopic scale that is comparable to empirical eeg .
( 2004 ) were also able to qualitatively reproduce the biphasic response through only the modulation of the inhibitory decay constant ( ipsp prolongation ) .
a further development in bojak and liley ( 2005 ) was the reparametrisation of postsynaptic potentials to more accurately fit experimental data .
alpha form had been used , i.e. , for an incoming presynaptic pulse at time t  =  0 the assumed psp response at t    0 was:4\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ r_\alpha(t)=\upgamma\gamma\exp(1)\times t \exp\left(-\gamma t\right ) , $ $ \end{document } as illustrated in fig .  2 for an ipsp
, this response has a maximum amplitude  at rise time 1/. unfortunately the psp decay time  for which the response falls again to /exp(1 ) is for the 
anaesthetic agents often potentiate gabaa induced currents and/or attenuate glutamate mediated ionic currents ( rudolph and antkowiak 2004 ) , thereby selectively prolonging the inhibitory ik without affecting the inhibitory rise time ik .
hence bojak and liley ( 2005 ) introduced eq .  2 yielding a response of the form5\documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 

\usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$ r(t)=\upgamma\gamma \exp\left(\tilde{\gamma}\delta\right)\times \frac{\exp\left(-\gamma t\right)-\exp\left(-\tilde{\gamma}t\right)}{\tilde{\gamma}-\gamma } , $ $ \end{document}with \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\delta\equiv\hbox{ln}(\tilde{\gamma}/\gamma)/(\tilde{\gamma}-\gamma).$$\end{document } this bi - exponential has the same maximum amplitude r(t  =   )  =  . however , the additional parameter \documentclass[12pt]{minimal } 
 \usepackage{amsmath } 
 \usepackage{wasysym } 
 \usepackage{amsfonts } 
 \usepackage{amssymb } 
 \usepackage{amsbsy } 
 \usepackage{mathrsfs } 
 \usepackage{upgreek } 
 \setlength{\oddsidemargin}{-69pt } 
 \begin{document}$$\tilde{\gamma}$$\end{document } can be used to adjust the decay time freely .
this allowed bojak and liley ( 2005 ) to directly incorporate into their model experimental data on the effect of isoflurane on psps ( maciver et  al .
1996 ; banks and pearce 1999 ; nishikawa and maciver 2000 ) , and to predict the eeg under anaesthesia by calculating the corresponding changes to the excitatory mean soma membrane potential he , see eq .  1 .
bojak and liley ( 2005 ) were able to show systematic changes in modelled eeg as a function of isoflurane concentration .
an additional selection criteria was then imposed upon these sets of data , which aimed at extracting sets producing a clear biphasic response ( bojak and liley 2005 ) .
the outcomes of bojak and liley s ( 2005 ) explorations showed some difference to that of steyn - ross et  al .
it is particularly worth highlighting that while both theories describe changes in macroscopic dynamics during anaesthesia , steyn - ross et  al .
( 2004 ) suggest a phase transition ( e.g. boiling water to steam ) , whereas bojak and liley ( 2005 ) posit a continuum transition .
bojak and liley ( 2005 ) highlight that the sharp rise in neural firing rate required for steyn - ross et  al.s ( 2004 ) phase transition is not supported by empirical findings of smooth suppression of action potential frequency with progressive anaesthesia .
( 2007 ) sought to further develop the work of steyn - ross et  al .
( 2004 ) and bojak and liley ( 2005 ) , by implementing a model parameterisation of slow voltage dependent mechanisms on excitatory firing rate dynamics .
the common sigmoidal relationship between mean soma membrane potential and mean firing rate was modified to included the influence of slow shifts in membrane potential .
this voltage dependant slow ionic leak , was motivated by the empirically observed switching ( i.e. up and down states ) in firing rates of neurons during anaesthesia and non - rem sleep , thought to be controlled by cation conductance .
( 2007 ) explored progressive increases in simulated desflurane anaesthesia , showing up state dominance during waking levels , up / down state switching ( delta ) at moderate anaesthesia and down state predominance during strong anaesthetic effect ( fig . 
this switching in firing rate between bursting up states and quiescent down states relates biophysically to the synchronised delta power local field potential predominately seen during deep anaesthesia and stage 3/4 sleep .
these findings further promote the ability to theoretically relate cellular experimental observations with bulk macroscopic cortical dynamics via mean field theoretics.fig . 
( 2007 ) simulated eeg ( mean soma membrane potential he ) at various concentrations of simulated anaesthesia by assuming that anaesthetics , in addition to increasing the duration of the ipsp , also cause slow adaptive changes in the shape of the mesoscopic firing rate function which relates mean soma membrane potential to mean population firing rate . ( 2 ) actual eeg recorded in children undergoing anaesthesia with the volatile halogenated agent desflurane .
( a ) 0 mac ( minimum alveolar concentration ) , ( b ) 1 mac , ( c ) >  1 mac and ( d ) 2 mac .
both real and simulated data show that increasing anaesthetic effect is associated with enhanced total eeg power and reductions in mean frequency .
( 2007 ) , with permission ( 1 ) molaee - ardekani et  al .
( 2007 ) simulated eeg ( mean soma membrane potential he ) at various concentrations of simulated anaesthesia by assuming that anaesthetics , in addition to increasing the duration of the ipsp , also cause slow adaptive changes in the shape of the mesoscopic firing rate function which relates mean soma membrane potential to mean population firing rate . ( 2 ) actual eeg recorded in children undergoing anaesthesia with the volatile halogenated agent desflurane .
( a ) 0 mac ( minimum alveolar concentration ) , ( b ) 1 mac , ( c ) >  1 mac and ( d ) 2 mac .
both real and simulated data show that increasing anaesthetic effect is associated with enhanced total eeg power and reductions in mean frequency .
the volatile general anaesthetic enflurane , among others , has been shown to paradoxically promote epileptiform activity ( modica et  al .
this pro - convulsant phenomenon is clearly at odds with the suppressive nature of anaesthetic compounds and lacks any clear physiological explanation .
liley and bojak ( 2005 ) explored the pro - convulsant properties of enflurane by parameterising the ipsp with data from in  vivo studies .
these findings suggested that a distinguishing factor between enflurane and other volatile anaesthetics ( e.g. isoflurane ) was a significant reduction in ipsp amplitude ( banks and pearce 1999 ) even though both agents increase net inhibition . because of the early developments in ipsp parameterisation liley and bojak ( 2005 ) were able to systematically explore how independent modulation of the ipsp amplitude influences system dynamics .
it was found that anaesthetic agents , such as enflurane , that reduce ipsp amplitude more than other anaesthetics , were more likely to produce eeg activity characteristic of epilepsy ( liley and bojak 2005 ) .
again liley and bojak ( 2005 ) used biologically constrained parameter sets which displayed suitable resting eeg spectra .
it was found using a bifurcation analysis of a spatially homogeneous model , that isofluranes minimal effect on ipsp amplitude did not produce any progression to system instabilities .
however , in over a third of the parameter sets the lower amplitude enflurane ipsps were shown to induce a strong shift from low firing , low amplitude resting state to an epileptic - like regime ( fig . 
8liley and bojak ( 2005 ) showed that systematically modulating ipsp shape and amplitude could induce high firing large amplitude activity in a mean field model parameterised to produce electroencephalographically plausible activity (  1/f low frequency activity and strong alpha resonance ) .
( b ) it is predicted that anaesthetics that more rapidly reduce the amplitude of the ipsp ( ik ) in addition to prolonging it ( ik ) , are more likely to produce high firing large amplitude activity electroencephalographic characteristic of epileptiform activity .
5 and 11 of liley and bojak ( 2005 ) liley and bojak ( 2005 ) showed that systematically modulating ipsp shape and amplitude could induce high firing large amplitude activity in a mean field model parameterised to produce electroencephalographically plausible activity (  1/f low frequency activity and strong alpha resonance ) .
( b ) it is predicted that anaesthetics that more rapidly reduce the amplitude of the ipsp ( ik ) in addition to prolonging it ( ik ) , are more likely to produce high firing large amplitude activity electroencephalographic characteristic of epileptiform activity .
figure adapted from figs .  5 and 11 of liley and bojak ( 2005 ) at approximately the same time wilson et  al .
similar to the work of liley and bojak ( 2005 ) , wilson et  al .
( 2006a ) showed too that it was the sensitivity towards ipsp amplitude which made enflurane unique in its effects when compared with isoflurane .
an additionally interesting suggestion put forward on the basis of these findings was that the anticonvulsant thiopentone , which has been shown to be proepileptic during administration with low doses of enflurane , may enhance enflurane s influence because of its negligible effect on ipsp amplitude ( wilson et  al .
( 2006a ) was the suggestion that the apparent sensitivity of macroscopic dynamics upon the magnitude of the drug induced peak ipsp amplitude may reflect a general phenomenon in which epileptic - like behaviour may be produced by interference with the integrative properties of dendrites .
more generally , mean field models have an immediate potential to assist in the development of physiologically meaningful methods for monitoring anaesthetic depth using eeg .
the use of biophysical models of cerebral cortical dynamics in parameterising or otherwise motivating eeg analysis methods for doa monitoring have the potential to more sensitively detect clinically relevant changes in the eeg given that they are derived from quantitative physiological descriptions of the anaesthetic process . at present
most monitoring techniques are based on the construction of discriminating statistical functions from the quantitative analysis of eeg , in large populations of normal and clinical groups undergoing a range of anaesthetic procedures , in order to capture the generalised quantitative changes in the eeg that show systematic variations with anaesthetic action . despite the early promise of such black - box approaches
, recent evidence suggests that they do not provide the clinician with any additional information regarding hypnotic state . in a recent study involving 2,000 patients avidan et  al .
( 2008 ) found that the use of bispectral index ( bis ) monitoring was not associated with a reduction in the incidence of awareness during anaesthesia when compared to other more standard techniques for monitoring anaesthetic depth .
this is significant as the bis is clinically the most widely used black - box eeg approach to monitoring anaesthetic depth . furthermore , the bis and other heuristic eeg monitoring approaches have been reported to be blind to the effects of a number of anaesthetic agents which are in routine clinical use , including nitrous oxide and the opioids ( e.g. rampil et  al . 1998 ; anderson and kakobsson 1992 ) . by using methods of analysis derived from physiological descriptions of the anaesthetic process , quantitative indexing of the eeg for
doa monitoring can perhaps be more robust to the variety of macroscopic physiological consequences of anaesthetic action as observed in the eeg .
in an attempt to understand the bulk dynamical properties of cortical tissue , researchers have found mean field theories of particular utility due to their ability to incorporate a large range of important physiological phenomena in a relatively tractable mathematical formulation . through developing biologically plausible models of cortical tissue that span a number of traditionally relevant scales
, a quantitative link can be made theoretically between a range of empirical measurements . theoretically spanning this explanatory gap between molecular / cellular drug action and macroscopic population dynamics assists in dealing with the often inscrutable complexities of brain dynamics .
essential to the success of this approach is the biological plausibility of not only the models but also the parameterisation of drug action .
as has been discussed above , some of the first attempts at using such an approach have already been fruitful in the study of anaesthesia .
these encouraging results add to the growing utility of dynamical systems theory in the analysis and conceptualisation of neural activity .
central to these limitations is the trade off between model utility and model complexity , where maintaining tractability is a serious technical challenge .
firstly , mean field models provide a tractable way of simulating the eeg on large spatial scales , but due to limited physiological and anatomical data often do so by assuming spatially homogeneous parameterisation . given the well know heterogeneity of neocortical architectonics , improving the spatial plausibility ( topology ) of these models , and comparing their differences , is a central issue ( e.g. qubbaj and jirsa 2007 ) .
additionally , the inclusion of more biophysical details to these models often incurs a disproportionate increase in model complexity .
therefore , as part of the improvements to neural field models more generally ( biophysical / topological accuracy ) , attempts should be made to simplify current models where possible in order to limit redundancy and improve parsimony ( e.g. kim and robinson 2007 ) .
more immediate future developments in mean field or neural mass modelling of cerebral cortex include the extension of the current literature of simulated anaesthesia , especially the exploration of dissociative anaesthetic agents such as ketamine , nitrous oxide and xenon whose pharmacological action is the antagonism of the excitatory glutamatergic nmda receptor ( hirota 2006 ) .
such developments naturally lead into other domains of brain function , such as learning and memory in the general case of nmda receptor functioning .
indeed , mean field models of anaesthesia have been extended into sleep ( e.g. wilson et  al .
2006b ) , where the similarities of the eeg spectrum and behavioural state suggest some shared putative mechanisms .
development of these models requires some judicial extension as their already large degrees of freedom must be constrained by empirical data .
some constraints can be provided by utilising mean field models for eeg doa monitoring as mentioned above , where clinical requirements insist on rapid online assessment of eeg dynamics and therefore limit computational complexity . together
these theoretical and practical improvements may produce a better understanding not only of the quantitative changes in the eeg during anaesthesia but also the physiological processes underlying conscious states .
 a great explanatory gap lies between the molecular pharmacology of psychoactive agents and the neurophysiological changes they induce , as recorded by neuroimaging modalities . 
 causally relating the cellular actions of psychoactive compounds to their influence on population activity is experimentally challenging . 
 recent developments in the dynamical modelling of neural tissue have attempted to span this explanatory gap between microscopic targets and their macroscopic neurophysiological effects via a range of biologically plausible dynamical models of cortical tissue . 
 such theoretical models allow exploration of neural dynamics , in particular their modification by drug action . 
 the ability to theoretically bridge scales is due to a biologically plausible averaging of cortical tissue properties . in the resulting macroscopic neural field , 
 individual neurons need not be explicitly represented ( as in neural networks ) . 
 the following paper aims to provide a non - technical introduction to the mean field population modelling of drug action and its recent successes in modelling anaesthesia . 
Introduction
Mean field modelling of neural activity
Modelling anaesthesia
Targets of anaesthetic action
Cellular network models of anaesthetic action
Mean field models of anaesthetic action
Clinical implications of mean field population models
Conclusion