intro intermediate state

references.bib

@@ -675,6 +675,24 @@ pages = {811--824}
   VOLUME = {104},
   PAGES = {2806-2820} }
 
+@incollection{Benda2013,
+	author = {Benda, Jan and Grewe, Jan and Krahe, R{\~A}ƒ{\^A}ƒ{\~A}‚{\^A}ƒ{\~A}ƒ{\^A}‚{\~A}‚{\^A}¼diger},
+	booktitle = {Animal Communication and Noise},
+	doi = {10.1007/978-3-642-41494-7_12},
+	editor = {Brumm, Henrik},
+	isbn = {978-3-642-41493-0},
+	language = {English},
+	pages = {331-372},
+	publisher = {Springer Berlin Heidelberg},
+	refid = {876},
+	series = {Animal Signals and Communication},
+	timestamp = {2014.02.03},
+	title = {Neural Noise in Electrocommunication: From Burden to Benefits},
+	url = {http://dx.doi.org/10.1007/978-3-642-41494-7_12},
+	volume = {2},
+	year = {2013},
+	bdsk-url-1 = {http://dx.doi.org/10.1007/978-3-642-41494-7_12}}
+
 @incollection{Benda2020,
   title={The physics of electrosensory worlds.},
   author={Benda, Jan},
@@ -6455,6 +6473,18 @@ groups and electrogenic mechanisms.},
   PAGES = {341--354}
 }
 
+@article{Wiesenfeld1995,
+	author = {Wiesenfeld, K. and Moss, F.},
+	journal = {Nature},
+	keywords = {noise, statistics},
+	pages = {33--36},
+	refid = {440},
+	timestamp = {2008.09.26},
+	title = {Stochastic resonance and the benefits of noise: from ice ages to crayfish and squids},
+	volume = {373},
+	year = {1995}
+}
+
 @Article{Wilkens2002,
   Title                    = {The electric sense of the paddlefish: a passive system for the detection and capture of zooplankton prey.},
   Author                   = {Lon A. Wilkens and Michael H. Hofmann and Winfried Wojtenek},

susceptibility1.tex

@@ -406,21 +406,20 @@ In this work, the influence of nonlinearities on stimulus encoding in the primar
 
 %with nonlinearities being observed in all sensory modalities
 
-Neuronal systems are inherently nonlinear \citealp{Adelson1985, Brincat2004, Chacron2000, Chacron2001, Nelson1997, Gussin2007, Middleton2007, Longtin2008}. A prominent example of a nonlinearity is rectification that is assumed to occur through the transduction machinery of inner hair cells \citealp{Peterson2019}, signal rectification in receptor cells \citealp{Chacron2000, Chacron2001} or the rheobase of action-potential generation \citealp{Middleton2007, Longtin2008}. Nonlinearity can be necessary to explain the behavior of complex cells in the visual system \citealp{Adelson1985}, to extract information about the stimulus \citealp{Barayeu2023} and to encode stimulus features as up- and down-strokes \citealp{Gabbiani1996}. 
+Nonlinear processes are fundamental in neuronal information processing. On the systemic level: deciding to take one or another action is a nonlinear process. On a finer scale, individual neurons are inherently nonlinear units: whether an action potential is elicited depends on the membrane potential to exceed a certain threshold\cite{Adelson1985, Brincat2004, Chacron2000, Chacron2001, Nelson1997, Gussin2007, Middleton2007, Longtin2008}. We can find nonlinearities in many sensory systems such as is rectification in the transduction machinery of inner hair cells \citealp{Peterson2019}, signal rectification in electroreceptor cells \citealp{Chacron2000, Chacron2001} or in complex cells in the visual system \citealp{Adelson1985}. In the auditory or the active electric sense, for example, nonlinear processes are needed to extract envelopes\citealp{Joris2004, Barayeu2023}. In conjunction with neuronal noise nonlinear mechanism facilitate the encoding of weak stimuli via stochastic resonance\citealp{Wiesenfeld1995, Stocks2000,Neiman2011fish}.
 
+While the sensory periphery can be often well described by linear models, upstream neurons no longer encode linearly. Rather, nonlinear processes are implemented to extract special stimulus features\citealp{Gabbiani1996}. In active electrosensation, the self-generated electric field (electric organ discharge, EOD)\citealp{Salazar2013} that is quasi sinusoidal in wavetype electric fish and acts as the carrier signal that is amplitude modulated in the context of communication\citealp{Benda2013, Fotowat2013, Walz2014, Henninger2018} as well as object detection and navigation\citealp{Fotowat2013, Nelson1999}. In social contexts the interference of the EODs of two interacting animals result in a new signal with a characteristic periodic amplitude modulation, the so-called beat. The beat amplitude is defined by the smaller EOD amplitude. The beat frequency is defined as the difference between the two EOD frequencies ($\Delta f = f-\feod{}$). Beats are common stimuli in different sensory modalities enabling rhythm and pitch perception in human hearing \citealp{Roeber1834, Plomp1967, Joris2004, Grahn2012} and providing context for electrocommunication in weakly electric fish \citealp{Engler2001, Hupe2008, Henninger2018, Benda2020}. Cutaneous electroreceptor organs in the skin are distributed over their body\citealp{Carr1982} and are tuned to the own field\citealp{Hopkins1976,Viancour1979}. P-type electroreceptor afferents (P-units) innervate these organs via ribbon synapses\citealp{Szabo1965, Wachtel1966} and project via the lateral line nerve to the hindbrain where they trifurcate and project onto pyramidal cells in the electrosensory lateral line lobe (EL)\citealp{Krahe2014}. The P-units ot the gymnotiform electric fish \lepto{} encode such amplitude modulations (AMs) by modulation of their firing rate\citealp{Gabbiani1996}. 
 
-The time-resolved firing rate of a neuron can be described by the Volterra series where the first-order term describes the linear contribution and all higher-order terms the nonlinear contributions \citealp{Voronenko2017}. In previous work, the second-order response of the Volterra series, the second-order susceptibility, was analytically retrieved based on leaky integrate-and-fire (LIF) models, where the input were two sine waves \citealp{Voronenko2017}. There second-order susceptibility was demonstrated to be very pronounced at specific input frequencies. A triangular nonlinear shape, with nonlinearities appearing if one of the input frequencies, the sum or the difference of the frequencies was equal to the baseline firing rate \fbase{} of the cell, was predicted. Since these nonlinear effects were especially pronounced if one signal was weak, they might influence the faint signal detection as it was observed in the field in the framework of the electrosensory cocktail party \citealp{Henninger2018}. 
+Even though the extraction of the AM itself requires a nonlinearity\citealp{Middleton2006,Stamper2012Envelope,Savard2011} the firing rate encodes its time-course faithfully\citealp{Gussin2007,Grewe2017,Savard2011}. In the context of social signalling among three fish we observe an AM of the AM, also referred to as second-order envelope or just social envelope\citealp{Middleton2006, Savard2011, Stamper2012Envelope}. Encoding this again requires nonlinearities and it was shown that a subpopulation of P-units is sensitive to envelopes\citealp{Savard2011}. Nonlinearities can thus be required to increase sensitivity for certain stimulus features.
+Theoretical work on leaky integrate-and-fire (LIF) model neurons revealed specific pattern of nonlinear interactions when the model was driven with two sinewaves\citealp{Voronenko2017}. This situation is reminiscent of the real situation when three electric fish interact. The active electrosensory system is then exposed to the two AMs arising from the interactions with either foreign fish's field. Previous recordings in the natural habitat showed interactions of three animals in which two animals interacted closely (strong signals) and were interrupted by an intruding animal. The intruder was detected at spatial distances that led to an extremely faint intruder signal while the strong signal of the other animal is present\citep{Henninger2018}. How can the intruder be detected despite the low amplitude and the presence of the dominating signal of the close-by animal? Could the nonlinearities described theoretically and in modelling studies facilitate the detection? As mentioned above, the firing rate of the P-units can be partially explained by a simple linear filter. The linear relation of the cellular response and the stimulus can be captured by the first-order susceptibility (or transfer function). As in Volterra series, higher-order terms describe the nonlinear interactions. We here apply the same approach as Voronenko and colleagues\citealp{Voronenko2017} and quantify the nonlinearity by estimating the second-order susceptibility. 
 
 
 %. The EOD
-In this work, nonlinearities are studied in the gymnotiform weakly electric fish of the species \lepto. These fish are equipped with an electric organ \citealp{Salazar2013}, that is constantly active and produces a quasi-sinusoidal electric organ discharge (EOD), with a long-term stable and fish-specific frequency (\feod, \citealp{Knudsen1975, Henninger2020}), that is used for electrolocation \citealp{Fotowat2013, Nelson1999} and communication \citealp{Fotowat2013, Walz2014, Henninger2018}. If two fish meet, their EODs interfere and result in a new signal with a characteristic periodic amplitude modulation, the so-called beat. The beat amplitude is defined by the smaller EOD of the encountered fish and is expressed in relation to the receiver EOD as a contrast. The beat frequency is defined as the difference between the two EOD frequencies ($\Delta f_{1} = f_{1}-\feod{}$). Beats are common stimuli in different sensory modalities enabling rhythm and pitch perception in human hearing \citealp{Roeber1834, Plomp1967, Joris2004, Grahn2012} and providing the context for electrocommunication in weakly electric fish \citealp{Engler2001, Hupe2008, Henninger2018, Benda2020}. Whereas the predictions of the nonlinear response theory \citealp{Voronenko2017} are applicable if the main driving force for a neuron are beats and not the whole signal, will be addressed in this work. 
+If two fish meet, their EODs interfere   Whereas the predictions of the nonlinear response theory \citealp{Voronenko2017} are applicable if the main driving force for a neuron are beats and not the whole signal, will be addressed in this work. 
 
 %apply to P-units,
 
-Cutaneous tuberous organs, that are distributed all over the body of these fish\citealp{Carr1982}, sense the actively generated electric field and its modulations.
+The EOD is encoded by the p-type electroreceptor afferents, the P-units. These fire phase-locked to the EOD of the fish with a spike at a similar phase in the EOD cycle, with the phase depending on the receptor location \citealp{Hladnik2023}. P-units fire in a probabilistic manner, resulting in an interspike interval (ISI) distribution with maxima at integer multiples of the EOD period.  In previous works, P-units have mainly been considered to be linear encoders \citealp{Xu1996, Benda2005, Gussin2007, Grewe2017}. Still, P-units exhibit nonlinear effects, especially for strong beat amplitudes \citealp{Nelson1997}. Although some P-units can be described by a linear decoder, other P-units require a nonlinear decoder \citealp{Chacron2004}.
-Within a single tuberous electroreceptor organ about 30 primary electroreceptors form ribbon synapses onto the dendrites of the same afferent fiber
-\citealp{Szabo1965, Wachtel1966} that projects via the lateral line nerve to the hindbrain. There, the afferent fiber synapses onto pyramidal neurons in the
-electrosensory lateral line lobe (ELL, \citealp{Maler2009a}). The EOD is encoded by the p-type electroreceptor afferents, the P-units. These fire phase-locked to the EOD of the fish with a spike at a similar phase in the EOD cycle, with the phase depending on the receptor location \citealp{Hladnik2023}. P-units fire in a probabilistic manner, resulting in an interspike interval (ISI) distribution with maxima at integer multiples of the EOD period.  In previous works, P-units have mainly been considered to be linear encoders \citealp{Xu1996, Benda2005, Gussin2007, Grewe2017}. Still, P-units exhibit nonlinear effects, especially for strong beat amplitudes \citealp{Nelson1997}. Although some P-units can be described by a linear decoder, other P-units require a nonlinear decoder \citealp{Chacron2004}.
 
 When the receiver fish with EOD frequency \feod{} is alone, a peak at \fbase{} is present in the power spectrum of the firing rate of a recorded P-unit (\subfigrefb{motivation}{A}, bottom). P-units also represent \feod{} in their firing rate, but this peak is beyond the range of frequencies addressed in this figure. The beat frequency in a two fish setting is represented in the spike trains, firing rate and the corresponding power spectrum of P-units (\subfigrefb{motivation}{B, C}). When three fish encounter, all their waveforms interfere with both beat frequencies being present in the superimposed signal (\subfigrefb{motivation}{D}, top). The power spectrum of the P-unit's firing rate contains both beat frequencies, but also nonlinear peaks at the sum frequency \bsum{} and the difference frequency \bdiff{} (\subfigrefb{motivation}{D}, bottom). This difference of the two beat frequencies is known as the social envelope \citealp{Stamper2012Envelope, Savard2011} that often emerges as the modulation of two beats in the superimposed signal. The question whereas P-units encode envelopes has been the subject of controversy, with some works not considering P-units as envelope encoders \citealp{Middleton2006}, while others identify some P-unit populations as successful envelope encoders \citealp{Savard2011}. In this work the second-order susceptibility at the sum and difference frequency will be characterized for the three fish setting, addressing whereas P-units encode envelopes and how this encoding is influenced by their baseline firing properties.