Exact solutions to cable equations in branching neurons with tapering dendrites

Neurons are biological cells with uniquely complex dendritic morphologies that are not present in other cell types. Electrical signals in a neuron with branching dendrites can be studied by cable theory which provides a general mathematical modelling framework of spatio-temporal voltage dynamics. Typically such models need to be solved numerically unless the cell membrane is modelled either by passive or quasi-active dynamics, in which cases analytical solutions can be reduced to calculation of the Green’s function describing the fundamental input-output relationship in a given morphology. Such analytically tractable models often assume individual dendritic segments to be cylinders. However, it is known that dendritic segments in many types of neurons taper, i.e. their radii decline from proximal to distal ends. Here we consider a generalised form of cable theory which takes into account both branching and tapering structures of dendritic trees. We demonstrate that analytical solutions can be found in compact algebraic forms in an arbitrary branching neuron with a class of tapering dendrites studied earlier in the context of single neuronal cables by Poznanski (Bull. Math. Biol. 53(3):457–467, 1991). We apply this extended framework to a number of simplified neuronal models and contrast their output dynamics in the presence of tapering versus cylindrical segments.


Introduction
Most neurons share a common structure consisting of a soma, an axon and dendrites. Dendrites are typically the most extended parts. Their distinct morphologies started to be appreciated by many scientists from the 1890s through the exemplary drawings of Ramón y Cajal [32], and gradually the significance of dendritic functions in single cell neuronal computation became apparent (informative overviews from both an experimental and a theoretical perspective can be found in the series of books [10,33,37,41]). A successful application of cable theory in modelling spatio-temporal voltage dynamics in dendritic arborisations can be attributed to Rall [31,35]. The idea behind cable theory is to build models of dendritic voltage dynamics using the analogy of electrical circuits. The electrical circuits can be either passive (linear) or active (nonlinear), mimicking the absence or presence of voltage-gated ion channels in dendritic membranes. While the combination of realistic dendritic morphologies and active voltage-gated ion channels restricts such models to be solved by only numerical methods using a compartmental approach, mathematical analysis is possible when dendritic voltage dynamics can be described by a linearised (quasi-active) model. Although the linearised models ignore any active dynamics for spike generation, they provide the fundamental foundation for a better understanding of neuronal signal filtration and integration in single dendritic/axonal cables (e.g. [15,26,34]) and in more complex morphologies (e.g. [17,18]).
Most models based on cable theory assume cylindrical dendritic segments, despite the fact that they are known to taper in many types of neurons, demonstrating initially rapid and then moderate decreases in dendritic radius from proximal to distal ends [3,6,20,42]. Earlier mathematical modelling of dendritic voltage dynamics on continuous tapering structures was also attributed to Rall and co-authors with the focus on single tapering segments or dendritic trees that can be reduced to equivalent tapering cables [16,31]. Later Poznanski [27] followed up the theoretical investigation of such tapering cables and identified the geometric types for continuous dendritic tapers that permit analytical solutions. In parallel, multiple studies focussed on finding analytical solutions in arbitrary branching neurons with cylindrical dendritic segments [2,19,[22][23][24][25]). More recently, Glenn and Knisley [11] developed a method of obtaining analytical solutions in dendritic trees with continuous tapering structures and linearised voltage dynamics by solving a recursive transcendental equation, based on the original method designed by Major et al. [24,25] for passive branching neurons with cylindrical dendritic segments.
In this work we use an alternative approach, the sum-over-trips framework, which was initially developed by Abbott et al. [2] from the path integral formulation. The original framework provided a method for calculating analytical Green's functions (voltage response functions given a Dirac delta current injected at some discrete location) in arbitrary branching dendrites with cylindrical segments and passive membranes. Later this method was generalised to support quasi-active cell membranes [8] and electrically coupled neuronal networks [39]. Here we combine the results of [27] and [8], and introduce an extended theoretical framework for calculating Green's functions in an arbitrary branching neuron with cylindrical or tapering segments and passive or quasi-active membranes. Moreover, we demonstrate that the set of geometric types identified in [27] makes it possible to describe the voltage dynamics in the entire branching structure by a unique Helmholtz equation, and thus allows us to apply the method of local point matching [43] for finding the Green's functions in compact algebraic forms.
The paper is structured as follows. In Sect. 2 we introduce and study a model of a single tapering dendritic cable of infinite length. We demonstrate how this model can be reduced to a Helmholtz equation and can then be solved by calculating the corresponding Green's function. Next, in Sect. 3 we consider an arbitrary branching neuron with tapering segments and quasi-active membranes, and extend the sum-over-trips framework for calculating Green's functions in such neurons. In Sect. 4 we consider a number of illustrative examples and demonstrate the application of the extended sum-over-trips framework. In particular, we justify a consideration of a parabolic taper as one of the biologically realistic dendritic geometries, and then study (i) a soma and dendrite model and (ii) a 'Y'-shaped dendritic tree model. Finally in Sect. 5 we provide a discussion of the further potential of this work. Some detailed mathematical derivations are collected in appendices for the interested reader.

A mathematical model of a single tapering dendritic cable 2.1 The generalised cable equation
We consider a single one-dimensional dendritic cable of an infinite length with its radius described by a smooth function r(x), where x ∈ R denotes the spatial location along the cable. The cell membrane is modelled by an electrical circuit consisting of a passive part and an active part which mimics the dynamics of ion channels. Passive dynamics of the cell membrane are governed by the membrane capacitance C m , the leak conductance g l and the resting membrane potential E l . The ion current generated by an arbitrary type of ion channel can be modelled as I ion = g ion (V )(V -E ion ), where g ion (V ) is the voltage-dependent conductance and E ion is the ion's reversal potential. According to the Hodgkin-Huxley formalism [14,15] the form of g ion (V ) depends on the dynamics of the gating variables for a particular type of ion channel. The generalised cable equation with a single type of voltage-gated ion channel can then be written down as where the axial resistivity of the dendritic cytosol R a is a constant and The term models the input in the form of a current with a time course I inj (t) injected at the location x 0 . A detailed derivation of Eq. (1) is provided in Appendix 1. It has been demonstrated earlier that the fully active (nonlinear) I ion dynamics can be well approximated by a quasi-active (linearised) current I ion [8,18]. Equation (1) can thus be reduced to the following linearised system (from now on all membrane potentials are measured from rest): where r ion and L ion are the effective resistance and inductance of the linearised ion channels respectively. Next we will demonstrate how a set of transformations can convert the system of Eqs. (4a) and (4b) into an equivalent form that is ideal for the subsequent mathematical analysis. By defining to be the characteristic length parameter [27,30,31], we can introduce a spatial scaling mapping where μ : R → R is bijective. Applying μ to Eqs. (4a) and (4b) leads to the following system of equations: where τ = C m /g l is the membrane time constant, and formally x = μ -1 (X) for μ -1 the inverse mapping of μ. Similarly to the approach in [27] for the passive tapering cable, we define to be the geometric ratio, and reduce Eq. (7a) to the following form: where Introducing a voltage transformation for we can further reduce Eq. (9) to where Note S φ : R → R is a bijective mapping, because φ(X) is always a positive constant given any X except the case when r(x) = 0. Applying the voltage transformation S φ to Eq. (7b) gives us Therefore the above steps convert the model (4a) and (4b) defined in the (x, t; V )coordinate into the following model in the (X, t; V * )-coordinate: Since both the spatial and voltage mappings μ and S φ are bijective, the two models (4a) and (4b) and (16a) and (16b) are equivalent.

The analytical solution
Applying the Laplace transform (L : f (t) → f (ω)) to Eqs. (16a) and (16b) and assuming zero initial conditions we obtain which can be reduced to the following compact form by simple substitution and rearrangement: where If γ (X, ω) is a constant in X, i.e. γ (X, ω) = γ (ω) for all X, Eq. (18) can be rewritten as where ∇ 2 is the Laplacian of V * on X ∈ R. Note Eq. (20) has the form of a one-dimensional inhomogeneous Helmholtz equation with a complex-valued wavenumber k = γ (ω) √ -1, which is analytically solvable. Under the assumption that the parameters τ , g l , r ion and L ion are constants, the special case when γ (X, ω) = γ (ω) is equivalent to the case of β(X) being simply a constant. One trivial scenario is when r(x) is a constant and therefore β(X) = 1, representing the case of a cylindrical cable. More generally, β(X) being a constant reduces Table 1 The six types of dendritic tapers that permit analytical solutions to a quasi-active cable equation. κ, L are positive constants and n ∈ Z. Modified from [27]. Note that outside of the taper T a p e r r a n g e [L + π κ n, L + π κ (n + 1 2 )] Trigonometric The six geometric types that permit analytical solutions in a single dendritic cable. The geometric ratios F(X) (in blue) were computed as listed in Table 1 with L = 1.5 and κ = π /3, except the Sinusoidal case with L = 0.15 and κ = π /2.7. The dendritic radii (in red) were obtained numerically by solving Eq. (8) for r(x).
For plotting purposes, we assume F(0) = r(0) = 1, and scale the x-coordinate by a factor of 10, where all coordinates and parameters are in arbitrary units Eq. (14) to a Riccati equation, solutions to which provide six geometric types [27]. These six types are listed in Table 1, and their geometric ratios F(X) and the corresponding dendritic radii r(x) are illustrated in Fig. 1. Our modelling framework allows one to have an arbitrary choice for a spatial coordinate system and therefore to properly control the taper range associated with each type of dendritic taper.
Assuming the dendritic cable belongs to one of the six geometric types, we can introduce an additional spatial scaling γ : X →x defined bȳ which normalises Eq. (20) into where The Green's function of Eq. (22) can be found as as Eq. (22) is a special case of the inhomogeneous Helmholtz equation without any boundary conditions, whose Green's function is well known. The derivation of Eq. (24) can be found in Appendix 2. It allows us to find the general solution to Eq. (22) as where C ⊂ C is a curve obtained by γ : R → C given any fixed ω ∈ C. Hence, in the (x, ω; V )-coordinate we have where Here we define z(x 0 , ω) to be the characteristic admittance at location x 0 , while z -1 (x 0 , ω) and r a (x 0 ) are more commonly referred as the characteristic impedance and the axial resistance respectively. If I inj (t) = δ(t), i.e. I inj (ω) = 1, Eq. (26) is reduced to the so-called transfer function of the model: which is also known as the transfer impedance [18].

Figure 2
The Green's function snapshots at three different time points in two (tapering and cylindrical) passive dendritic cables demonstrating the asymmetry in the spread of voltage in the tapering case (a slice of the tapering cable is illustrated in red and that of the cylindrical cable is in dashed black). All the coordinates and parameters are in arbitrary units: 2R a g l = 1, R a = 4π 3/2 , τ = 4, radius of the cylindrical dendrite r c = 1, radius of the tapering cable at the origin r(x = 0) = 1. The tapering dendritic cable belongs to the exponential type with κ = π /3. The Dirac-delta input is placed at x = 0 Applying the inverse Laplace transform (L -1 : f (ω) → f (t)) to Eq. (26), we obtain the solution in the time domain as where given in (30). In the limit r ion → ∞, the system is reduced to a purely passive model with the Green's function The system can be further reduced to the cylindrical case by considering r(x) to be some constant r c . Examples of the Green's function profiles in the passive tapering and cylindrical dendritic cables are illustrated in Fig. 2, clearly demonstrating the asymmetry in the dispersion of the voltage in the tapering case.

The neuronal morphology
Here we consider an arbitrary branching neuron as illustrated in Fig. 3. The global morphology is modelled as a graph Γ = (N, S), where N is the set of somatic, branching and terminal nodes (vertices) and S is the set of dendritic segments (edges). In particular, Γ is a metric graph whose weighted edges are associated with the physical lengths l i for all dendritic segments i ∈ S. Similar to the spatial transformations (6) and (21) for the case of a single tapering cable, we define the mappings μ i : x → X and γ i : X →x locally on each segment i as assuming that each tapering segment is either cylindrical or belongs to one of the six types listed in Table 1. For convenience we simplify the notations by defining the global spatial Although Γ * and Γ share the same graph structure (see Fig. 3), their edge weights are different;

The voltage dynamics
The voltage dynamics in the entire branching structure Γ with quasi-active membranes can be fully characterised by the dynamics V i (x, t) for x ∈ (0, l i ) on each dendritic segment i ∈ S and the boundary conditions at all the nodes χ ∈ N . In particular, V i (x, t) is governed by the following system of equations similar to Eqs. (4a) and (4b) for a single cable: where the input which, comparing to Eq. (3), has an extra term δ ij , a Kronecker delta (i.e. δ ij = 1 if i = j and 0 otherwise) denoting the input location on segment j. Although each segment can have different individual parameters indicated by index i, all the equations can be normalised if we investigate the equivalent dynamics on Γ * instead of Γ . In particular, by introducing the spatial transformations γ · μ in (33a) and (33b), a voltage transformation S φ similar to Eq. (11) and the Laplace transformation L, we can reduce Eqs. (35a) and (35b) into the following equation defined on each segment i: where Note and all the other variables are in the same forms as in Sect. 2, but with index i. Since Eq. (37) is in the same form as the normalised Helmholtz equation (22) for all i ∈ S, (∇ 2 -1) is a linear operator acting on all the segments and thus on the entire graph Γ * .

The boundary conditions
All the boundary conditions are assumed to be governed by two physical laws: continuity of potentials and conservation of currents. Here we consider three types of boundary conditions: terminal, branching and somatic nodes (see Fig. 3). Working on Γ * instead of Γ requires to transform these boundary conditions by γ · μ, S φ and L, details of which are provided in Appendix 3.

Terminal node
A terminal node χ T is either killed or sealed, and it can be described by one of the following boundary conditions: or

Branching node
For a branching node χ B with K attached individual segments, the following two conditions are required: where r a,i (χ B ) = R a /(πr 2 i (χ B )).

Somatic node
A lumped soma χ S can be treated as a special node connecting K segments, and the somatic membrane potential V S (t) can be described by a quasi-active model as where the constants C S , g S , r S and L S determine the linearised dynamics of the somatic membrane.

The general solution
Let H ij (x,ȳ) be the Green's function for the operator (∇ 2 -1) on Γ * . We can find the analytical solution to the system (37) as where C j ⊂ C is a curve obtained by γ · μ : (0, l j ) → C j given any fixed ω ∈ C. When Γ * is an infinite single cable, Eq. (44) is reduced to Eq. (25), as H ij (x,ȳ) is replaced by H ∞ (x,ȳ). In general cases, we can use the sum-over-trips framework introduced earlier in [2,8], which demonstrate that the Green's function H ij (x,ȳ) is closely linked to H ∞ (x,ȳ). In particular, it can be written down as where a trip on Γ * starting at x on segment i and terminating y on segment j is similar to a random walk, but with the restrictions of changing direction only at the nodes. Since the physical trip length l trip (x, y) ≥ 0, Eq. (24) for each trip is reduced to wherel trip (x, y) = γ · μ(l trip (x, y)) is the normalised trip length. In addition, each trip weight (the summand in Eq. (45)) has to satisfy the boundary conditions at all the nodes along the trip, which is encoded in the trip coefficients A trip (ω). We will introduce in Sect. 3.5 the new rules for constructing A trip (ω). More details on the explanation and construction of trips and the sum-over-trips framework can be found in [2,8].
Finally H ij (x,ȳ) in Eq. (45) allows us to construct the Green's function for V i (x, y, ω) on Γ as and the voltage responses for an arbitrary applied current as where From now on we assume for simplicity that the parameters describing the electrical properties of the cell membrane are identical for all dendritic segments and thus drop the subscript i for them, because our focus is on dendritic tapers of individual segments only. The above expressions are similar to those in the case of a single cable introduced in Sect. 2.2, but they are defined locally on each segment. From now on we use the term Green's function (in the frequency domain) when referring to G ij (x, y, ω) rather than H ij (x,ȳ), unless stated otherwise.

Rules for constructing trip coefficients
A trip coefficient A trip (ω) for any trip is calculated by firstly initialising its value to be 1, and then multiplying it by a node factor α nm (ω) every time when the trip travels from segment n to m. Node factors are essentially derived from, and thus encode complete information of, boundary conditions at individual nodes. A detailed derivation is given in Appendix 3. This derivation serves as a constructive proof that the function H ij (x,ȳ) in Eq. (45) constructed using the appropriate node factors satisfies the required boundary conditions in Γ . Similar derivations for cylindrical dendrites can be found in [2,8]. Below we list the node factors for three types of boundary conditions discussed in Sect. 3.3, omitting ω for compactness.

Terminal node
A trip has to reflect at a terminal node χ T . For a killed terminal, the node factors is For a sealed terminal, where The expression for ξ m has a similar form to Eq. (10) for an infinite single cable, but it is defined locally for each segment. Since each node factor is associated with the direction of the trip travelling away from the corresponding node, this direction is assumed to be positive orientation of the local spatial coordinate.

Branching node
At a branching node χ B a trip can either reflect or pass through. If the trip reflects at χ B , the node factor is If the trip passes through χ B , where are defined for segment m and n connected at χ B .

Somatic node
At a somatic node χ S , the node factors have the same expressions as the branching node factors, that is, for a reflective trip, and, for a transitive trip, whereas and is the somatic admittance.

The method of local point matching
In order to find the analytical Green's function of the system (35a) and (35b) in Γ , we apply the extended sum-over-trips framework with the new node factors described above.
Although the theoretical convergence of the infinite sum in the Green's function (47) can be proved using a similar argument to [1], truncation of terms is always required for any numerical computation of Green's functions in an arbitrary branching neuron [5]. Instead, the method of local point matching [43] can be used to find Green's functions in compact algebraic forms. We refer the reader to [43] for a summary of the algorithmic steps of this method, which requires modifications for the case of dendritic tapers on the following three entities: • the normalisation of the spatial coordinates; • the voltage transformation; • the node factors. As a result, we can obtain the Green's function (47) in the following form: where J y can be found by the method of local point matching [43]. We can then use it in Eq. (48) to calculate V i (x, y, ω) and take the inverse Laplace transform to get V i (x, y, t).

Some examples with parabolic dendritic taper 4.1 Basic geometry of a parabolic dendritic segment
In the upcoming examples of simplified neuronal models, we assume that all dendritic segments are either cylindrical or parabolic in geometry. The choice of the parabolic taper is justified by a number of experimental observations [3,6,20,42], indicating that the tapering slope is steeper at the start of the segment and is getting flatter towards the end. In particular, the radius of a parabolic dendritic segment satisfies for x ∈ [0, l 0 ], as we choose the local x-coordinate for the segment so that r 0 = r(0), r 1 = r(l 0 ) are the initial and terminal radii respectively, where l 0 is the length of the segment, and a = (1 -√ r 1 /r 0 )/l 0 defines the slope of the parabola. Under the assumption [r (x)] 2 1, this parabolic taper can be well approximated by the exponential type (see [12] for the proof ), whose geometric ratio can be found as where λ(x) and X are defined in Eqs. (5) and (6), and κ = 3aλ(0)/2 is a positive constant characterising the slope of F(X). Figure 4 justifies the fact that there are negligible difference between the two geometries. Therefore, we model all parabolic dendritic segments by the exponential type for its mathematical tractability demonstrated in Sect. 2. Combining Eqs. (8) and (67) we obtain which defines the spatial mapping μ locally for each parabolic taper. In addition we note that both hyperbolic sine and hyperbolic cosine types asymptotically approach the exponential type in the limit l 0 → ∞ (cf. Table 1), while the other three types have biologically unrealistic geometries since their r(x) functions are concave (see Fig. 1). It has also been computationally validated in [9] and later mathematically proved in [4] that the parabolic taper is optimal in the transfer of current signals along a single dendritic cable.

A soma and dendrite model: parabolic versus cylindrical cable
Here we consider a simple model of a single dendritic cable with one end (x = 0) attached to a lumped soma and the other end (x = l 0 ) being a sealed terminal node (see Fig. 5). Our earlier study of this model with a cylindrical cable of radius r c led to the the following Green's function for the somatic response [43]: where and z S is given in Eq. (64). The subscript c in the above expressions denotes the parameters in the cylindrical cable. The characteristic length parameter λ c and the axial resistance r a,c are simply constants. In the parabolic dendritic cable with radius described by Eq. (66), the geometry can be rewritten in the form of Eq. (67) with the spatial scaling μ given by Eq. (68), and then normalised by The subscript p here and below denotes the parameters in the parabolic cable. There are two node factors associated with this model: the reflective somatic node factor (2p S p -1) and the reflective terminal node factor (2p T p -1), where We can then apply the method of local point matching to construct the following linear system of equations (with the arrows in Fig. 5): where h(x) = exp(-x). Solving the system (74a) and (74b), we obtain which gives The Green's function in the frequency domain can then be calculated using Eq. (65). In particular, the Green's function for the somatic response can be found as For the comparison between the cylindrical and parabolic models, we focus only on their dendritic geometries and consider all other parameters to be the same. For the parabolic dendritic cable we assume r 0 = r c and r 1 = 0. We can then rewrite the somatic Green's functions (69) and (77) as The geometric parameters: r c = 1 μm for the cylindrical model, r 0 = 1 μm and r 1 = 0 μm for the parabolic model, r S = 12.5 μm for both models. The electrical parameters of the dendritic and somatic membranes are the same, and identical in both models: C m = 1 μF · cm -2 , g -1 l = 2000 Ω · cm 2 , R a = 100 Ω · cm, r ion = 1000 Ω · cm 2 , L ion = 5 H · cm 2 In the limiting case of a semi-infinite dendritic cable, and in the limiting case of an infinitesimal dendrite, Both limits can be derived analytically from Eqs. (78) and (79), but can also be heuristically obtained from a geometric perspective. In the first limit (80) the parabola asymptotically becomes a cylinder as κ → 0, while in the second limit (81) both models reduce to a single lumped soma with the somatic impedance z -1 S . To compare the two models in the presence of the quasi-active membranes, the preferred frequency Ω * can be introduced as the frequency at which the magnitude of the Green's function is maximised. In Fig. 6 we plot the preferred frequencies for both models using the Green's functions (78) and (79) with x = y = 0. It is clear that the two curves are close to each other when the dendritic length l 0 is either extremely small or large, which means the two models behave similarly. This can be inferred from the two limits (80) and (81). We note that the curve for the parabolic model is monotonic while the curve for the cylindrical model is not. Within the dendritic length range investigated, the difference is maximised around l 0 = 150 μ m, with which the two models should behave the most differently; for this dendritic length, the time profiles of the two models' somatic responses to a chirp current are illustrated in Fig. 6.

A soma and dendrite model: parabolic versus compartmental cable
Here we still study the same parabolic model as in the previous section, but compare it to a compartmental model, consisting of an array of N cylindrical dendritic segments. We assume that these cylindrical segments share the same length l 0 /N while their radii are successively decreasing. Such compartmental models are commonly used in computational work (e.g. [9,21,42]) to approximate continuous tapering structures. In order to investigate the geometric effect on voltage amplitudes, we consider the two models to be purely passive, and set the dendritic geometries in such a way that the total membrane areas of the two models are equal, while assuming all the other parameters to be the same. In particular, we choose the dendritic radius of segment i ∈ {1, 2, . . . , N} in the compartmental model to be where and the function r(x) is defined by Eq. (66).
Since the two models are both passive, we only study the somatic responses at steady state (as t → ∞) to a step current of strength I step switched on at time t 0 . Equivalently in the frequency domain, we have Therefore, once we obtain the Green's functions in the frequency domain, we can write down the voltage at steady state simply as by applying the final value theorem for the Laplace transform. We can clearly see from Fig. 7 that the curve of V ss (0, y) for the parabolic model is concave, while the curve for the compartmental model is convex on all individual segments. We checked that this property is always valid regardless of parameter choices. One could take Eqs. (78) and (79) as an example: when ω = 0, G c (0, y, 0) is almost exponential in y while G p (0, y, 0) follows a power law in y and the power 3/2 + 3γ p /2κ > 3. This results in a large range of y that permits V ss,p (0, y) > V ss,c (0, y). We also verified that when N is very large (e.g. N = 1000), the V ss,c (0, y) curve became indistinguishable from the V ss,p (0, y) curve in Fig. 7 (not illustrated here). However, considering computational expenses, a dendritic taper is usually approximated by a compartmental model with only a few segments (N < 10 in [9,21,42]). Errors of such approximations are not negligible on segments whose radii are tiny (typically the segment near terminal), even when N is relatively large (e.g. N = 100), because the input (characteristic) impedances are extremely different in the two models.

A 'Y'-shaped dendritic tree: parabolic versus cylindrical segments
Here we consider a parabolic model of a passive neuron with a simplified 'Y'-shaped dendritic tree. The dendritic tree consists of one cylindrical primary dendritic segment and  Fig. 8. We also consider a cylindrical model which differs from the parabolic model by only reducing the geometry of the secondary dendritic segments into cylinders using Eq. (82) with N = 1.
We investigated how the somatic responses at steady state V ss (-l 0 , y) vary with the length of the primary dendritic segments l 0 (see Fig. 9a), using the Green's functions obtained by the method of local point matching (the detailed calculations can be found in Appendix 4).
A noticeable difference can be found locally on the secondary segments. The scale of the difference is modulated by l 0 , which is mainly due to the signal loss along the primary segments. Introducing the normalised difference Δ norm = V ss,p (-l 0 , y)/V ss,c (-l 0 , y) -1, we can see from Fig. 9b that the overall trend of Δ norm remains unchanged; Δ norm > 0 for most of y ∈ (0, l 1 ]. In contrast, if the primary dendritic segment tapers but the secondary segments not, we observed Δ norm > 0 for most of y ∈ [-l 0 , 0) instead. When both primary and secondary segments are tapering, Δ norm > 0 can be observed for most of y ∈ [-l 0 , l 1 ].  Figure 9 (a) Somatic responses at steady state of the cylindrical (red) and parabolic (blue) models to different input locations on the 'Y'-shaped dendritic tree. Both models have a soma and a cylindrical primary dendrite of a radius r c = r 0 = 1 μm whose parameters are the same as in Fig. 7. The cylindrical model has two cylindrical secondary dendrites, while the parabolic model has two parabolic secondary dendrites. Note that for each case the origin of the spatial coordinate is chosen to be in the branching node, placing the somatic node at the coordinate x = -l 0 . (b) The normalised difference Δ norm between the steady-state somatic responses of the two models

Discussion
In this paper we present a unified approach for calculating the Green's functions in a quasiactive (or passive) neuron with an arbitrary branching dendritic structures and tapering segments. It extends the sum-over-trips framework introduced for passive cylindrical dendrites [2] and quasi-active cylindrical dendrites [8], and generalises the work in [27] for single tapering cables. We also demonstrate how the solutions can be found in compact algebraic forms (instead of being represented as infinite sums) using the previously developed method of local point matching [43]. The obtained compact Green's function solutions allow one to conduct mathematical analysis and efficient numerical simulations for better understanding the role of dendritic morphology on neuronal signal modulations. These solutions can be naturally reduced to the models with cylindrical dendritic structures to recover the results in [2,8,43] by considering r(x) = r c . Note that, since the spatial scaling factor of the injected current has been incorporated in the form of the Green's function in this paper (instead of being present in the function I 0 (x, t) as in previous studies [2,8,43]), the recovered Green's functions for the cylindrical models differ from the previously obtained Green's functions by a factor Dr a = 1/(2πCr c ). The form of the Green's function given here thus simplifies the reciprocity identity as by absorbing location dependent coefficients (a proof is given in Appendix 5).
Our comparison between the tapering and cylindrical models indicates that the parabolic geometry tends to increase somatic voltage responses driven by distal dendritic inputs and that this increase can be attributed to the local dendritic taper alone. The earlier work of Bird and Cuntz [4] mathematically justifies that the parabolic dendritic taper provides the optimal signal transfer from the input location to the soma in a model of a single tapering cable. Our study considered a still simplified, but nevertheless a more advanced 'Y'-shaped tree model, which indicates that this optimal signal transfer property is also applicable to more complex dendritic branching structures.
As it stands this modelling framework excludes any active dynamics that can be attributed to the voltage-gated ionic channels. Mathematical analysis of the dendritic models with active properties is often impossible due to nonlinearities in the underlying equations, with a few exceptions of dendritic models with hot spots [28,29] or simplified piecewise linear models of active spines used in frameworks such as the Spike-Diffuse-Spike (SDS) type model [7,38]. It would be interesting to extend the SDS model from the cylindrical to tapering dendritic segments and investigate the geometric effects on the speed of wave propagation at the level of a single cable as well as on branching structures. Another possible extension of this work is to go beyond a single cell model to a network level. Considering a model of spatially extended neurons coupled by electrical synapses (gap-junctions) it is possible to extend the work in [39] and derive the node factors for gap-junctional boundary conditions located on tapering dendritic segments. The Green's function for the whole network can then be efficiently computed using the method of local point matching [43]. Such extensions would allow us to rigorously investigate the effects of dendritic tapers on the voltage dynamics at the level of a single cell and electrically coupled neuronal networks with the dendritic and somatic membrane models being passive, quasi-active or active. Finally, the proposed framework is compatible with stochastic cable theory [40], allowing one to expand a recent work of Gowers et al. [13] of calculating the firing rates of branching neurons with tapering dendrites. (88c) For simplicity we assume that the cross section of the cable is always a disc, and thus the membrane surface area of this section is Substituting Eqs. (88a), (88b) and (88c) into Eq. (87) and taking the limit Δ → 0, we obtain Without loss of generality we assume that an injected current I inj (t) is applied to a small section (x 0 , x 0 + Δ). In the considered limit Δ → 0 it becomes an input current at a discrete point x = x 0 , which gives At the same time the axial current I(x) is assumed to satisfy Ohm's law, that is, where is the axial resistance of the section (x, x + Δ), and is the cross-sectional area. We can thus rewrite Eq. (92) and find the axial current as which reduces to This leads to which is the generalised cable equation given in (1).
where η j (y) = 1 z j (y)φ j (Y ) . (105) To simplify the notation, we omit ω in all the functions here and we work simultaneously with the two spatial coordinates linked by the bijective mapping μ : x → X. Without loss of generality we also assume that a node is located at x = 0 (and thus X = 0) in all three types of nodes, imposing the boundary conditions listed in Sect. 3.3 on the Green's function (104). In practice any node factor should be evaluated with its local parameters, and its value is independent of the choice of local spatical coordinates.

Terminal node
Under the above assumptions the dendritic morphology is simply a semi-infinite cable. Applying the sum-over-trips framework there are only two trips, a direct trip and a reflective one. Equation (104) thus reduces to where α k , k ∈ {o, c} is the reflective node factor for killed and sealed terminals respectively. If the terminal node is killed, the boundary condition (40) is imposed, equivalently, Substituting the Green's function (106) into Eq. (107) we obtain which simply gives as all the other functions on the left hand side are positive. If the terminal node is sealed, the boundary condition (41) is imposed, which gives ∂G ∂x (0,y) = 0.
By the chain rule which implies by using the derivatives previously found in Eqs. (121a) and (121b). Thus Eq. (126) is reduced to which gives It can be easily proved that the above expressions are identical to the node factors (61) and (62).

Appendix 4: The Green's function in the 'Y'-shaped dendritic tree
To find the Green's function for the 'Y'-shaped dendritic model discussed in Sect. 4.4, we apply the sum-over-trips framework and use the method of local point matching [43]. We denote the dendritic segments by 0 (primary) and 1 (secondary), and the nodes by S (somatic), B (branching) and T (sealed terminal), and thereby identify all node factors below.
• At the somatic node, where • At the branching node, where and (135) • At the two sealed terminal nodes, All the variables should be evaluated according to the definitions in Sect. 3.5. The details are omitted here. Assuming that the output is located at the soma and applying the method of local point matching we obtain where L y = γ 0 μ 0 (y), L y-l 0 = γ 1 μ 1 (yl 0 ), and construct the Green's function using Eq. (65).

Appendix 5: The reciprocity identity for an arbitrary dendritic tree
Equation (18) of the model in the Laplace domain is the second-order linear ordinary differential equation and can thus be rewritten in the Sturm-Liouville form. The differential operator is self-adjoint, and therefore the Green's function must be symmetric in space [36], supporting the reciprocity principle G(x, y) = G(y, x) for the Green's function in an infinite single cable. Here we demonstrate that this principle is also valid on an arbitrary dendritic tree with any boundary conditions discussed in Sect. 3.3. For any trip and its reversal configuration, namely 'pirt' , we can write down the trip coefficients and find the ratio between them as A trip A pirt = α ik 1 α k 1 k 2 α k 2 k 3 · · · α k n-1 k n α k n j α jk n α k n k n-1 α k n-1 k n-2 · · · α k 2 k 1 α k 1 given that all reflective node factors have cancelled each other. Assuming k 0 = i, k n+1 = j, and segments k m and k m+1 are connected at node χ m for all m ∈ {0, 1, 2, . . . , n}, using either Eq. (58) or (62) we have In addition, any given trip and its reversal trip are identical in their lengths, regardless of any scaling (i.e.l trip =l pirt ), which gives Therefore using Eq. (47) we obtain which gives the reciprocity identity (86) either in the frequency or time domain. There is an immediate corollary from the reciprocity identity (86) stating that G ij (x, y) is continuous with respect to the input location y. To prove this we recall the continuity of membrane potentials at the nodes discussed in Sect. 3.3, which ensures that G ji (y, x) is continuous in y, and thus G ij (x, y) = G ji (y, x) is also continuous in y. Although the corollary seems to be trivial, its important consequence is the case of having the input location y at a node (e.g. the soma). In such a scenario we can place the input at y on any of the dendritic segments that are attached to that node and then consider the distance between y and the node to be zero. The corollary guarantees the uniqueness of the solution, even if we choose different attached segments in such circumstances. Therefore the framework can be applied directly to any input location.