Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебники / Hearing - From Sensory Processing to Perception Kollmeier 2007

.pdf
Скачиваний:
157
Добавлен:
07.06.2016
Размер:
6.36 Mб
Скачать

498

ITD(ms)

320 pps

1

0

−1

0

50

100

 

Time (ms)

 

ITD(ms)

H.S. Colburn et al.

640 pps

1

0

−1

0

50

100

 

Time (ms)

 

Fig. 2 Firing rates to 40-Hz cosine amplitude-modulated pulse trains of 320 and 640 pps

solely due to the adaptive mechanisms in the IC model. When high-rate pulse trains are amplitude modulated, ITD sensitivity to on-going stimulus is observed in a strong response (Fig. 2). This result supports the hypothesis that an adaptation mechanism with amplitude-modulation enhances ITD sensitivity.

4Cell Membrane Factors in Modulation Effects

In contrast with the Jeffress (1948) coincidence model for ITD sensitivity, real coincidence mechanisms involve neurons with refractory properties that prevent neural firings and degrade ITD tuning for high frequency inputs (Cook et al. 2003; Dasika et al. 2005). Refractory behavior may also play a significant role in restoring neural discharges with amplitude modulation. Explored here are effects of synaptic parameters and sinusoidally amplitude modulated (SAM) pulse-train stimuli on the ITD sensitivity of a model neuron with an active membrane.

4.1Model Description

The neural model and stimulus generation are similar to a multi-compartment MSO model used in a previous study (Zhou et al. 2005) with two modifications: 1) the cell model has only a soma compartment with the set of ionic channels, and 2) inputs to the model simulate the peripheral responses to amplitude-modulated electrical-pulse stimulations instead of tones. With a probability proportional to the amplitude of the pulse, a pulse triggers a change in post-synaptic conductance with maximum value Ge. The input

Models of Neural Responses to Bilateral Electrical Stimulation

499

pulse train is described as X(t)=Xamp(t)Xfine(t), where the pulse-train carrier is

Xfine(t)=Σk p(t-kT) and the modulated amplitude is Xamp(t)=A (1+ sin(2πƒmt)). The amplitude A is a Gaussian random variable [N(1,σ)], ƒm is the modula-

tion frequency, and 1/T is the carrier pulse rate. Effects of σ, ƒm, T and Ge on ITD tuning are explored below.

4.2Results

Figure 3 shows that the synaptic strength influences the neural entrainment to pulse trains. The model cell responds to each input pulse at low pulse rates but not at high pulse rates when the synaptic strength is weak (i.e., low Ge, upper traces). Entrainment to the input pulse-train improves when the Ge value is increased (middle trace on right). This result is consistent with in vitro observations, which show that activities of specific potassium channels block neural responses to high-frequency stimuli (Reyes et al. 1996). The model action potential has a short duration and fast recovery time (a few milliseconds) due to the combined activities of fast-sodium and high-threshold and low-threshold potassium channels on the soma (Rothman and Manis 2003). The combination of these temporal characteristics and the strength of the conductance determine the relative refractory period (RRP), which accounts for the reduced neural activity.

ITD computations involve temporal summations of bilateral inputs, and either the in-phase or out-of-phase responses may be in the regime where the RRP has a pronounced effect on neural discharges. Noting that the effective rate of out-of-phase inputs is twice the in-phase rate (ignoring strength), it is not surprising that the synaptic strength may have different effects for these conditions. Figure 4 shows ITD tuning curves for the model cell in response to un-modulated pulse trains for three synaptic strength levels. The sub-threshold, near-threshold and super-thresholds levels of strength are defined in terms of unilateral response probabilities. Specifically, the entrainment probability of model discharges to a low-rate (100 pps) unilateral, pulse-trains are at 0%, 14%, and 100% for the three

A B

T=10msT=2ms

Fig. 3 The effect of the membrane refractoriness on neural entrainment to input pulse trains with constant amplitude: A model responses when conductance Ge is low and the inter-pulse interval (T=10 ms) is larger than the relative refractory period of the model cell – note perfect entrainment; B model responses when the inter-pulse interval (T=2 ms) is comparable to the relative refractory period

500

H.S. Colburn et al.

500 pps

 

300

 

 

300

Spike/sec

200

 

 

200

100

 

 

100

 

0−1

0

1

0−1

 

 

ITD(msec)

 

 

1000 pps

0 1

ITD(msec)

ge= 0.5nS(sub)

ge= 0.9nS(near)

ge= 1.5nS(super)

Fig. 4 ITD tuning for unmodulated pulse trains (500 and 1000 pps) for several values of synaptic conductance Ge. Tuning is very sensitive to synaptic strengths (Ge). For the sub-threshold Ge condition, membrane refractoriness decreases model responses to in-phase inputs (relative to near-threshold Gecase), resulting in no ITD-tuning. For the near-threshold Ge condition, membrane refractoriness decreases model responses to out-of-phase inputs, resulting in better ITD tuning. For the super-threshold Ge condition, summed inputs overcome membrane refractoriness at most ITD values, resulting in saturated ITD tunings

threshold levels. Simulation results indicate that refractory processes lead to reduced discharges for in-phase inputs with the sub-threshold Ge and for out-of-phase inputs with the near-threshold Ge, and has less effect on model discharges with a super-threshold Ge. As a result, the corresponding ITD tuning at these three conditions show different overall activities, dynamic ranges, and tuning widths.

Amplitude modulating the inputs reduces the effective synaptic strength (and therefore neural firing) during half of the modulation cycle, thereby minimizing the effects of refractoriness on subsequent inputs. However, the resultant ITD tuning can be either improved or degraded dependent upon the synaptic strength level. Figure 5 shows that amplitude modulation improves ITD tuning at sub-threshold Ge condition and degrades ITD tuning at near-threshold Ge condition. This is due to neural release from the RRP for either in-phase or out-of-phase inputs. For the sub-threshold Ge condition (A and B), amplitude modulation releases in-phase responses from membrane refractoriness and ITD tuning is improved relative to that to un-modulated pulse-trains. For the near-threshold Ge condition (C and D), amplitude modulation releases out-of-phase responses from membrane refractoriness and ITD tuning is degraded relative to that to un-modulated pulse-trains.

Finally, adding amplitude variation decreases the effective synaptic strength by reducing the coincidence of unilateral inputs, and therefore could be used to improve ITD tuning. Results (not shown) indicate that an increase in jitter (larger σ) decreases the response rate such that ITD dependence may increase or decrease with increases in σ.

Models of Neural Responses to Bilateral Electrical Stimulation

501

 

150

A

500 PPS

 

150

B

1000 PPS

 

 

 

 

 

 

Spike/sec

100

 

 

 

100

 

 

 

50

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

 

0−1

0

1

0

−1

0

1

 

200

C

 

 

150

D

 

 

Spike/sec

 

 

 

100

 

 

 

 

 

 

 

 

 

 

100

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 −1

0

1

0−1

0

1

 

 

 

ITD(msec)

 

 

 

ITD(msec)

 

fm=50Hz

fm=100Hz

unmodulated

Fig. 5 Amplitude modulation (ƒm=50 or 100 Hz) can either improve or degrade ITD tuning to pulse-train stimuli (500 or 1000 pps): A,B the sub-threshold Ge condition; C,D the nearthreshold Ge condition

5Effects of Input Rate and Phase on ITD Sensitivity

The output of real a coincidence detector is also highly dependent on the rate and synchrony of its inputs. Principal MSO neurons are believed to derive ITD sensitivity by coincidence detection of their bilateral inputs, which are phase locked to auditory stimuli. This section explores why higher-rate, more highly synchronized responses in the auditory nerve during electrical stimulation may produce a loss in ITD sensitivity. Simulations of MSO neurons under acoustic and electrical stimulation of the auditory periphery were compared using the single compartment model.

5.1Model Description

The single-compartment model MSO cell here is nearly identical to that in Sect. 4, with identical channel dynamics, conductances, and capacitance (Zhou et al. 2005), and Ih calculated in the manner of Rothman and Manis (2003). The model MSO cell has 20 bilateral inputs (10 per side) from explicitly modeled auditory nerve fibers. A 500-Hz tone was the driving stimulus for both acoustic and electric input models. Each stimulus period for each fiber, it is determined if a spike occurs and when according to a Gaussian distribution with σ set by the nominal synchronization index (SI) (Zhou et al. 2005). Average input spike rate and SI were set according to condition, and phase delays between electric inputs were introduced as described below.

502

H.S. Colburn et al.

5.2Results

Figure 6 shows the rate-ITD characteristics of the model MSO cell for acoustic inputs and three variations of electric inputs. Synaptic strength was set and maintained such that simulated acoustic inputs (120 spikes/s for each input, SI=0.7) produced rate-ITD sensitivity similar to real MSO cells (Goldberg and Brown 1969; Yin and Chan 1990). Electrical stimulation was modeled as 250 spikes/s for each input, with SI=0.99. In simulated electrical stimulation without phase-dispersion (top curve), the high-rate, highly phase-locked inputs produced entrainment and saturation of the rate-ITD curve at good phase, and the resulting “monaural coincidences” produced high discharge rates at all phases, including bad phase. To simulate the possible effect of neural delays that may normally compensate for the cochlear traveling wave delay, phase dispersion between inputs was introduced such that added phases were equally spaced across the stimulus period (bottom curve). With this high degree of phase-dispersion, ITD sensitivity was eliminated, and few output discharges occurred across the range of ITD. To simulate the possible effect of neural plasticity after cochlear implantation, phase dispersion was limited to half the stimulus period. With this moderate phase dispersion, enhanced ITD sensitivity occurred in terms of the amplitude of the rate-ITD curve, however with slightly broader peaks (broader tuning) than for the lower rate acoustic stimulation.

 

700

model MSO responses to acoustic and electric stimulation

 

 

 

 

 

 

 

 

 

 

 

 

 

simulated acoustic inputs: rate 120 spk/s, sync index 0.7

 

 

 

600

 

simulated electric inputs: rate 250 spk/s, sync index 0.99

 

 

 

 

simulated electric inputs: rate 250 spk/s, high phase-dispersion

 

 

(spikes/sec)Rate

 

 

simulated electric inputs: rate 250 spk/s, moderate phase-dispersion

 

500

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

400

 

 

 

 

 

 

 

 

Discharge

300

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

0−2

−1.5

−1

−0.5

0

0.5

1

1.5

2

ITD (ms)

Fig. 6 ITD sensitivity of a model MSO neuron in response to acoustic and electric stimulation

Models of Neural Responses to Bilateral Electrical Stimulation

503

6Summary and Conclusions

Cochlear implant studies of ITD sensitivity in human psychoacoustics and animal electrophysiology lead to a consideration of ITD-dependent neural activities. Our simulation results suggest stimulation strategies to improve ITD sensitivity for CI users. Results show that increased neural firing to SAM stimuli can be explained by either release from neural adaptation or from membrane refractoriness. Model results also indicate that, dependent upon the synaptic strength level, the introduction of amplitude modulation can either improve or degrade ITD tuning to high-frequency pulse trains. These results suggest that, to enhance the ITD sensitivity of individual neurons, the effective parameter range of amplitude-modulated stimuli would be small and vary across ITD-sensitive cells, which have different membrane properties and synaptic input strengths. Further studies are needed to explore input-parameter combinations that could further optimize tuning to carrier ITDs.

Acknowledgments. This work was supported by the US Public Health Service NIH/NIDCD grant R01 DC05775 Bertrand Delgutte, PI.

References

Cai HM, Carney LH, Colburn HS (1998) A model for binaural response properties of inferior colliculus neurons. II. A model with interaural time difference-sensitive excitatory and inhibitory inputs and an adaptation mechanism,” J Acoust Soc Am 103:494–506

Cook DL, Schwindt PC, Grande LA, Spain WJ (2003) Synaptic depression in the localization of sound. Nature 421:29–30

Dasika VK, White JA, Carney LH, Colburn HS (2005) Effects of inhibitory feedback in a network model of avian brain stem. J Neurophysiol 94:400–414

Goldberg JM, Brown PB (1969) Responses of binaural neurons of dog superior olive to dichotic time stimuli: some physiological mechanisms of sound localization. J Neurophysiol 32:516–523

Jeffress LA (1948) A place theory of sound localization. J Comp Physiol Psychol 41:35–39

Kiang NYS, Moxon EC (1972) Physiological considerations in artificial stimulation of the inner ear. Ann Otol Rhinol Laryngol 81:714–730

Lane CC, Delgutte B (2005) Neural correlates and mechanisms of spatial release from masking: single-unit and population responses in the inferior colliculus. J Neurophysiol 94:1180–1198. Epub 2005 Apr 27

Lane CC, Kopco N, Delgutte B, Shinn-Cunningham BG, Colburn HS (2005) A cat’s cocktail party: psychophysical, neurophysiological and computational studies of spatial release from masking. In: Pressnitzer D, de Cheveign A, McAdams S, Collet L (eds) Auditory signal processing: physiology, psychoacoustics, and models. Springer, Berlin Heidelberg New York

Moxon E (1967) Electrical stimulation of the inner ear of cat. Doctoral dissertation, MIT Elec Eng Dept

Poon B (2006) Sound localization and interaural time sensitivity with bilateral cochlear implants. PhD Dissertation, M.I.T. Health Science and Technology, Cambridge, MA

Reyes AD, Rubel EW, Spain WJ (1996) In vitro analysis of optimal stimuli for phase-locking and time-delayed modulation of firing in avian nucleus laminaris neurons. J Neurosci 16:993–1007

504

H.S. Colburn et al.

Rothman J, Manis P (2003) Kinetic analyses of three distinct potassium conductances in ventral cochlear nucleus neurons. J Neurophysiol 89:3083–3096

Smith ZM (2006) Binaural interactions with bilateral electric stimulation of the cochlea. PhD Dissertation, M.I.T. Health Science and Technology, Cambridge, MA

Yin TCT, Chan JCK (1990) Interaural time sensitivity in the medial superior olive of the cat. J Neurophysiol 64:465–488

Zhou Y, Carney LH, Colburn HS (2005) A model for interaural time difference sensitivity in the medial superior olive: interaction of excitatory and inhibitory synaptic inputs, channel dynamics, and cellular morphology. J Neurosci 25:3046–3058

54 Neural and Behavioral Sensitivities to Azimuth Degrade with Distance in Reverberant Environments

SASHA DEVORE1, ANTJE IHLEFELD2, BARBARA G. SHINN-CUNNINGHAM2,

AND BERTRAND DELGUTTE1

1Introduction

Reverberation poses a challenge to sound localization in rooms. In an anechoic space, the only energy reaching a listener’s ears arrives directly from the sound source. In reverberant environments, however, acoustic reflections interfere with the direct sound and distort the ongoing directional cues, leading to fluctuations in interaural time and level differences (ITD and ILD) over the course of the stimulus (Shinn-Cunningham et al. 2005). These effects become more severe as the distance from sound source to listener increases, which causes the ratio of direct to reverberant energy (D/R) to decrease (Hartmann et al. 2005; Shinn-Cunningham et al. 2005).

Few neurophysiological and psychophysical studies have systematically examined sensitivity to sound source azimuth as a function of D/R (Rakerd and Hartmann 2005). Here we report the results of two closely-integrated studies aimed at characterizing the influence of acoustic reflections like those present in typical classrooms on both the directional sensitivity of auditory neurons and the localization performance of human listeners. We used low-frequency stimuli to emphasize ITDs, which are the most important binaural cue for sounds containing low-frequency energy (MacPherson and Middlebrooks 2002; Wightman and Kistler 1992). We find that reverberation reduces the directional sensitivity of low-frequency, ITD-sensitive neurons in the cat inferior colliculus (IC), and that this degradation becomes more severe with decreasing D/R (increasing distance). We show parallel degradations in human sensitivity to the azimuth of low-frequency noise.

1Eaton-Peabody Laboratory, Massachusetts Eye and Ear Infirmary, Boston, MA, USA, sashad@mit.edu, Bertrand_Delgutte@meei.harvard.edu

2Hearing Research Center, Boston University, Boston, MA, USA, ihlefeld@bu.edu, shinn@bu.edu

Hearing – From Sensory Processing to Perception

B. Kollmeier, G. Klump, V. Hohmann, U. Langemann, M. Mauermann, S. Uppenkamp, and J. Verhey (Eds.) © Springer-Verlag Berlin Heidelberg 2007

506

S. Devore et al.

2Single-Unit Neurophysiology

2.1Methods

Methods for recording from low-frequency, ITD-sensitive neurons in the IC of anesthetized cats were as described by Hancock and Delgutte (2004). We focused on measuring neural responses as a function of source azimuth in simulated rooms.

Binaural room impulse responses (BRIRs) were simulated using the roomimage method (Allen and Berkeley 1979) for a pair of receivers corresponding to the left and right ears of a cat in the center of a simulated reverberant room (11 × 13 × 3 m). We did not include a model of the head in the simulations, so that the resulting BRIRs contained ITD but essentially no ILD cues. BRIRs were calculated for azimuths spanning the frontal hemifield (−90° to 90°) at distances of 1 m and 3 m with respect to the midpoint of the receivers. Anechoic impulse responses were created by time-windowing the direct wavefront from the 1-m reverberant BRIRs. The direct-to-reverberant energy ratio (D/R) was calculated as the ratio of the energy in the direct sound (timewindowed from the BRIR) to the energy of the remaining impulse response. An overall D/R was determined for each source distance by averaging across all azimuths. Virtual room stimuli were created by convolving the BRIRs with a reproducible 400-ms broadband noise burst. The first 400 ms of the resulting signals were presented to dial-in-urethane anesthetized cats over calibrated, closed acoustic systems.

Neural responses were measured as a function of source azimuth for each virtual room condition (anechoic, 1 m, and 3 m). We typically used 11 azimuths (15° steps) or, occasionally, 7 azimuths (30° steps). The noise stimulus was repeated 16 times at each azimuth, with random order across azimuths. We computed the average firing rate by counting the number of action potentials over the stimulus duration.

2.2Results

We measured neural responses as a function of azimuth for 25 IC units from 7 cats. Rate-azimuth curves for two typical units are shown in Fig. 1A. Neural rate responses depend strongly on source azimuth in the anechoic condition (D/R = ∞ dB), with a preference for contralateral azimuths. Reverberation reduces the range of firing rates over all source azimuths (a “demodulation”), although rates still vary systematically with azimuth. To quantify these observations, we define the relative response range as the range of firing rates for a given room condition normalized by the range of firing rates in the anechoic condition. The relative range for the anechoic condition is 1, by definition. For most neurons, the relative ranges in the reverberant conditions are less

Neural and Behavioral Sensitivities

507

(spikes/sec)RateFiring

A

 

 

 

 

 

 

 

 

RangeRelative

B

 

 

 

20

 

 

 

 

 

 

 

 

1.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.0

 

 

 

 

10

D/R:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∞ dB

 

 

 

 

 

 

 

 

0.5

 

 

 

 

 

0 dB

 

 

 

 

 

 

 

 

 

 

 

 

 

−10 dB

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

unit sd049b

 

unit sd054b

 

0

 

 

 

 

−90 −60 −30

0

30

60 90−90 −60 −30

0

30

60

90

 

0

−10

 

 

Azimuth (º)

 

Azimuth (º)

 

 

 

Direct-to-Reverberant Energy Ratio (dB)

 

 

 

 

 

 

 

 

 

 

Fig. 1 A Mean neural response rate (±s.e.m.) as a function of source azimuth for two neurons in the cat IC. B Mean (large filled symbols) relative response range as a function of D/R. Small symbols indicate relative response range for individual units

Relative Information Transfer (%)

50

40

30

20

10

0

0

−10

 

Direct-to-Reverberant Energy Ratio (dB)

Fig. 2 Mean (large filled symbols) rIT as a function of D/R. Small symbols show rIT for individual units

than 1, indicating a compression of the rate-azimuths curves; furthermore, relative range decreases with decreasing D/R (Fig. 1B).

Neural sensitivity to azimuth depends not only on the response range, but also on the variability in responses at each azimuth. To quantify sensitivity, we computed the mutual information (MI) between the source azimuth and the neural spike counts for each neuron and room condition. MI characterizes the precision with which the source azimuth can be estimated from the neural spike counts without making additional assumptions about the neural code. To compare across experiments using different numbers of stimuli, MI was normalized by the stimulus entropy to get the relative information transfer (rIT). Figure 2 shows that the rIT systematically decreases with decreasing D/R.