## Decompose LFP into components

Moderator: wwlytton

breakwave922
Posts: 90
Joined: Mon May 19, 2014 10:39 pm

### Decompose LFP into components

I have a network model consisting of thousands of neurons, and I calculated LFP using 'xtra' . In this way, the calculation of LFP is based on the recording of i_membrane (transmembrane current). But now, I want to decompose the LFP into components reflecting different types of sources. For instance, I want to know the contribution of opto-channel current that I built on my cell to the LFP. Is there any method that I can calculate this?

What I plan to do is, record certain channel current, and replace the transmembrane currents with it in the equation for calculating LFP, and the result will be the contribution to LFP from this current source. Is this the correct way to do?

ted
Posts: 5752
Joined: Wed May 18, 2005 4:50 pm
Location: Yale University School of Medicine
Contact:

### Re: Decompose LFP into components

Assuming that one or more cells are embedded in a conductive medium that is linear, the potential Vi at any point i in the medium is

Code: Select all

``````Vi = SUMMA Zx_j Im_j
j``````
where the index j ranges over all compartments in the model, Im_j is the total (i.e. ionic + capacitive) transmembrane current for compartment j, Zx_j is the transfer impedance between i and j. Note that Vi and Im_j are the Fourier (or Laplace) transforms of the corresponding time domain variables vi_j and im_j.

If the medium is nondispersive (can be treated as purely resistive), the formula can be simplified to

Code: Select all

``````vi = SUMMA Rx_j im_j
j``````
where Rx_j is the transfer resistance between i and j.

For any particular component x of membrane current it is possible to calculate the contribution of that component to vi

Code: Select all

``````vxi = SUMMA Rx_j imx_j
j``````

breakwave922
Posts: 90
Joined: Mon May 19, 2014 10:39 pm

### Re: Decompose LFP into components

Thanks, Ted. I guess it is clearer to me now.
I'm wondering, if the channel current is implemented as point process mechanism, I can simply use the distance between the position of the channel and the electrode to calculate contributed LFP. However, if the channel current is implemented as distributed(density) mechanism, how can I calculate the distance to the electrode to calculate contributed LFP, as the channel is distributed all over the compartment.

ted
Posts: 5752
Joined: Wed May 18, 2005 4:50 pm
Location: Yale University School of Medicine
Contact:

### Re: Decompose LFP into components

At large distances (more than about 10 x segment length) it is reasonable to treat the segment as a point source of current located at the segment center. A more fastidious approach is to treat each segment as a line source with the same length as the segment, centered at the location of the segment center, and oriented parallel to the line that joins the segment ends. Segment orientation and center location can be calculated by treating the section's pt3d data as piecewise linear functions of range. The mathematical basis of the line source approximation has been published more than once, e.g. ModelDB entry 84589.

breakwave922
Posts: 90
Joined: Mon May 19, 2014 10:39 pm

### Re: Decompose LFP into components

Thanks Ted.

Another question I forgot to discuss with you was, if I calculate single channel current's (implemented as point-process) contribution to LFP, according to what you suggested, I can obtain it by using this channel's current weighted it by distance and medium conductivity. I'm wondering in this way, do I only consider the channel current flow into the cell? But how about the current flow out of the cell thru passive channel? If I only consider the current flow in at certain point to calculate it's contribution to LFP, is it correct?

ted