Analysis of inverse stochastic resonance and the long-term firing of Hodgkin–Huxley neurons with Gaussian white noise

In order to explain the occurrence of a minimum in firing rate which occurs for certain mean input levels μ$\mu$ as noise level σ$\sigma$ increases (inverse stochastic resonance, ISR) in Hodgkin–Huxley (HH) systems, we analyze the underlying transitions from a stable equilibrium point to limit cycle and vice-versa. For a value of μ$\mu$ at which ISR is pronounced, properties of the corresponding stable equilibrium point are found. A linearized approximation around this point has oscillatory solutions from whose maxima spikes tend to occur. A one dimensional diffusion is also constructed for small noise. Properties of the basin of attraction of the limit cycle (spike) are investigated heuristically. Long term trials of duration 500000 ms are carried out for values of σ$\sigma$ from 0 to 2.0. The graph of mean spike count versus σ$\sigma$ is divided into 4 regions _R1,…,_R4$R_1,\dots ,R_4$, where _R3$R_3$ contains the minimum associated with ISR. In _R1$R_1$ transitions to the basin of attraction of the rest point are not observed until a small critical value of σ=_σc1$\sigma ={\sigma }_{c_1}$ is reached, at the beginning of _R2$R_2$. The sudden decline in firing rate when σ$\sigma$ is just greater than _σc1${\sigma }_{c_1}$ implies that there is only a small range of noise levels 0<σ<_σc1$0<\sigma <{\sigma }_{c_1}$ where repetitive spiking is safe from annihilation by noise. The firing rate remains small throughout _R3$R_3$. At a larger critical value σ=_σc2$\sigma ={\sigma }_{c_2}$ which signals the beginning of _R4$R_4$, the probability of transitions from the basin of attraction of the equilibrium point to that of the limit cycle apparently becomes greater than zero and the spike rate thereafter increases with increasing σ$\sigma$. The quantitative scheme underlying the ISR curve is outlined in terms of the properties of exit time random variables. In the final subsection, several statistical properties of the main random variables associated with long term spiking activity are given, including distributions of exit times from the two relevant basins of attraction and the interspike interval.     

Non-Abelian fusion rules from an Abelian system

We demonstrate the emergence of non-Abelian fusion rules for excitations of a two dimensional lattice model built out of Abelian degrees of freedom. It can be considered as an extension of the usual toric code model on a two dimensional lattice augmented with matter fields. It consists of the usual C(Z_p)$\mathbb{C}\left({\mathbb{Z}}_p\right)$ gauge degrees of freedom living on the links together with matter degrees of freedom living on the vertices. The matter part is described by a n$n$ dimensional vector space which we call H_n$H_n$. The Z_p${\mathbb{Z}}_p$ gauge particles act on the vertex particles and thus H_n$H_n$ can be thought of as a C(Z_p)$\mathbb{C}\left({\mathbb{Z}}_p\right)$ module. An exactly solvable model is built with operators acting in this Hilbert space. The vertex excitations for this model are studied and shown to obey non-Abelian fusion rules. We will show this for specific values of n$n$ and p$p$, though we believe this feature holds for all n>p$n>p$. We will see that non-Abelian anyons of the quantum double of C(S₃)$\mathbb{C}\left(S_3\right)$ are obtained as part of the vertex excitations of the model with n=6$n=6$ and p=3$p=3$. Ising anyons are obtained in the model with n=4$n=4$ and p=2$p=2$. The n=3$n=3$ and p=2$p=2$ case is also worked out as this is the simplest model exhibiting non-Abelian fusion rules. Another common feature shared by these models is that the ground states have a higher symmetry than Z_p${\mathbb{Z}}_p$. This makes them possible candidates for realizing quantum computation.     

Three PT-symmetric Hamiltonians with completely different spectra

We discuss three Hamiltonians, each with a central-field part H₀$H_0$ and a PT-symmetric perturbation igz$igz$. When H₀$H_0$ is the isotropic Harmonic oscillator the spectrum is real for all g$g$ because H$H$ is isospectral to H₀+g²/2$H_0+g^2/2$. When H₀$H_0$ is the Hydrogen atom then infinitely many eigenvalues are complex for all g$g$. If the potential in H₀$H_0$ is linear in the radial variable r$r$ then the spectrum of H$H$ exhibits real eigenvalues for 0<g<g_c$0<g<g_c$ and a PT phase transition at g_c$g_c$.     

Transition between SU(4) and SU(2) Kondo effect

Motivated by experiments in nanoscopic systems, we study a generalized Anderson, which consist of two spin degenerate doublets hybridized to a singlet by the promotion of an electron to two conduction bands, as a function of the energy separation δ$\delta$ between both doublets. For δ=0$\delta =0$ or very large, the model is equivalent to a one-level SU(N$N$) Anderson model, with N=4$N=4$ and 2 respectively. We study the evolution of the spectral density for both doublets (ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$ and ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$) and their width in the Kondo limit as δ$\delta$ is varied, using the non-crossing approximation (NCA). As δ$\delta$ increases, the peak at the Fermi energy in the spectral density (Kondo peak) splits and the density of the doublet of higher energy ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$ shifts above the Ferrmi energy. The Kondo temperature T_K (determined by the half-width at half maximum of the Kondo peak in density of the doublet of lower energy ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$) decreases dramatically. The variation of T_K with δ$\delta$ is reproduced by a simple variational calculation.     

Thermodynamics properties of an anisotropic Ising antiferromagnet in both external longitudinal and transverse fields

We have studied the anisotropic two-dimensional nearest-neighbor Ising model with competitive interactions in both uniform longitudinal field H$H$ and transverse magnetic field Ω$\Omega$. Using the effective-field theory (EFT) with correlation in cluster with N=1$N=1$ spin we calculate the thermodynamic properties as a function of temperature with values H$H$ and Ω$\Omega$ fixed. The model consists of ferromagnetic interaction _Jx$J_x$ in the x$x$ direction and antiferromagnetic interaction _Jy$J_y$ in the y$y$ direction, and it is found that for H/_Jy∈[0,2]$H/J_y\in [0,2]$ the system exhibits a second-order phase transition. The thermodynamic properties are obtained for the particular case of λ=_Jx/_Jy=1$\lambda =J_x/J_y=1$ (isotropic square lattice).     

The H2+ molecular ion: Low-lying states

Matching for a wavefunction the WKB expansion at large distances and Taylor expansion at small distances leads to a compact, few-parametric uniform approximation found in Turbiner and Olivares-Pilon (2011). The ten low-lying eigenstates of H₂⁺ of the quantum numbers (n,m,Λ,±)$(n,m,\Lambda ,\pm )$  with n=m=0$n=m=0$ at Λ=0,1,2$\Lambda =0,1,2$, with n=1$n=1$, m=0$m=0$ and n=0$n=0$, m=1$m=1$ at Λ=0$\Lambda =0$ of both parities are explored for all interproton distances R$R$. For all these states this approximation provides the relative accuracy ?10⁻⁵$?10^{-5}$ (not less than 5 s.d.) locally, for any real coordinate x$x$ in eigenfunctions, when for total energy E(R)$E\left(R\right)$ it gives 10-11 s.d. for R∈[0,50]$R\in [0,50]$  a.u. Corrections to the approximation are evaluated in the specially-designed, convergent perturbation theory. Separation constants are found with not less than 8 s.d. The oscillator strength for the electric dipole transitions E1$E1$ is calculated with not less than 6 s.d. A dramatic dip in the E1$E1$ oscillator strength f_{1sσg−3pσu}$f_{1s{\sigma }_g-3p{\sigma }_u}$ at R∼R_eq$R\sim R_{eq}$ is observed. The magnetic dipole and electric quadrupole transitions are calculated for the first time with not less than 6 s.d. in oscillator strength. For two lowest states (0,0,0,±)$(0,0,0,\pm )$ (or, equivalently, 1sσ_g$1s{\sigma }_g$ and 2pσ_u$2p{\sigma }_u$ states) the potential curves are checked and confirmed in the Lagrange mesh method within 12 s.d. Based on them the Energy Gap between 1sσ_g$1s{\sigma }_g$ and 2pσ_u$2p{\sigma }_u$ potential curves is approximated with modified Pade Re^−R[Pade(8/7)](R)$R{e}^{-R}\left[Pade(8/7)\right]\left(R\right)$ with not less than 4-5 figures at R∈[0,40]$R\in [0,40]$ a.u. Sum of potential curves E_1sσg+E_2pσu$E_{1s{\sigma }_g}+E_{2p{\sigma }_u}$ is approximated by Pade 1/R[Pade(5/8)](R)$1/R\left[Pade(5/8)\right]\left(R\right)$ in R∈[0,40]$R\in [0,40]$ a.u. with not less than 3-4 figures.     

Effects of colored noise and noise delay on a calcium oscillation system

As a calcium oscillations system is in steady state, the effects of colored noise and noise delay on the system is investigated using stochastic simulation methods. The results indicate that: (1) the colored noise can induce coherence bi-resonance phenomenon. (2) there exist three peaks in the R–_τ0$R\text{–}{\tau }_0$ (R$R$ is the reciprocal coefficient of variance, and _τ0${\tau }_0$ is the self-correlation time of the colored noise) curves. For the same noise intensity Q=1$Q=1$, the Gaussian colored noise can induce calcium spikes but the white noise cannot do this. (3) the delay time can improve noise induced spikes regularity as _τ0${\tau }_0$ is small, and R$R$ has a significant minimum with increasing τ$\tau$ as _τ0${\tau }_0$ is large. (4) large values of ζ$\zeta$ reduce noise induced spikes regularity.     

Is space-time symmetry a suitable generalization of parity-time symmetry?

We discuss space-time symmetric Hamiltonian operators of the form H=H₀+igH^′$H=H_0+ig{H}^{\prime }$, where H₀$H_0$ is Hermitian and g$g$ real. H₀$H_0$ is invariant under the unitary operations of a point group G$G$ while H^′$H^{\prime }$ is invariant under transformation by elements of a subgroup G^′$G^{\prime }$ of G$G$. If G$G$ exhibits irreducible representations of dimension greater than unity, then it is possible that H$H$ has complex eigenvalues for sufficiently small nonzero values of g$g$. In the particular case that H$H$ is parity-time symmetric then it appears to exhibit real eigenvalues for all 0<g<g_c$0<g<g_c$, where g_c$g_c$ is the exceptional point closest to the origin. Point-group symmetry and perturbation theory enable one to predict whether H$H$ may exhibit real or complex eigenvalues for g>0$g>0$. We illustrate the main theoretical results and conclusions of this paper by means of two- and three-dimensional Hamiltonians exhibiting a variety of different point-group symmetries.     

Quantum spherical model with competing interactions

We analyse the phase diagram of a quantum mean spherical model in terms of the temperature T$T$, a quantum parameter g$g$, and the ratio p=−_J2/_J1$p=-J_2/J_1$, where _J1>0$J_1>0$ refers to ferromagnetic interactions between first-neighbour sites along the d$d$ directions of a hypercubic lattice, and _J2<0$J_2<0$ is associated with competing antiferromagnetic interactions between second neighbours along m≤d$m\le d$ directions. We regain a number of known results for the classical version of this model, including the topology of the critical line in the g=0$g=0$ space, with a Lifshitz point at p=1/4$p=1/4$, for d>2$d>2$, and closed-form expressions for the decay of the pair correlations in one dimension. In the T=0$T=0$ phase diagram, there is a critical border, _gc=_gc(p)$g_c=g_c\left(p\right)$ for d≥2$d\ge 2$, with a singularity at the Lifshitz point if d<(m+4)/2$d<(m+4)/2$. We also establish upper and lower critical dimensions, and analyse the quantum critical behavior in the neighborhood of p=1/4$p=1/4$.     

Complex scale-free networks with tunable power-law exponent and clustering

We introduce a network evolution process motivated by the network of citations in the scientific literature. In each iteration of the process a node is born and directed links are created from the new node to a set of target nodes already in the network. This set includes m$m$ “ambassador” nodes and l$l$ of each ambassador’s descendants where m$m$ and l$l$ are random variables selected from any choice of distributions _pl$p_l$ and _qm$q_m$. The process mimics the tendency of authors to cite varying numbers of papers included in the bibliographies of the other papers they cite. We show that the degree distributions of the networks generated after a large number of iterations are scale-free and derive an expression for the power-law exponent. In a particular case of the model where the number of ambassadors is always the constant m$m$ and the number of selected descendants from each ambassador is the constant l$l$, the power-law exponent is (2l+1)/l$(2l+1)/l$. For this example we derive expressions for the degree distribution and clustering coefficient in terms of l$l$ and m$m$. We conclude that the proposed model can be tuned to have the same power law exponent and clustering coefficient of a broad range of the scale-free distributions that have been studied empirically.     

The Bethe lattice treatment of sound attenuation for a spin- 3/2 Ising model

The sound attenuation phenomena is investigated for a spin- 3/2 Ising model on the Bethe lattice in terms of the recursion relations by using the Onsager theory of irreversible thermodynamics. The dependencies of sound attenuation on the temperature (T$T$), frequency (w$w$), Onsager coefficient (γ$\gamma$) and external magnetic field (H$H$) near the second-order (_Tc)$\left(T_c\right)$ and first-order (_Tt)$\left(T_t\right)$ phase transition temperatures are examined for given coordination numbers q$q$ on the Bethe lattice. It is assumed that the sound wave couples to the order-parameter fluctuations which decay mainly via the order-parameter relaxation process, thus two relaxation times are obtained and which are used to obtain an expression for the sound attenuation coefficient (α)$\left(\alpha \right)$. Our investigations revealed that only one peak is obtained near _Tt$T_t$ and three peaks are found near _Tc$T_c$ when the Onsager coefficient is varied at a given constant frequency for q=3$q=3$. Fixing the Onsager coefficient and varying the frequency always leads to two peaks for q=3,4$q=3,4$ and 6 near _Tc$T_c$. The sound attenuation peaks are observed near _Tt$T_t$ at lower values of external magnetic field, but as it increases the sound attenuation peaks decrease and eventually disappear.     

A new entropy based on a group-theoretical structure

A multi-parametric version of the nonadditive entropy S_q$S_q$ is introduced. This new entropic form, denoted by S_a,b,r$S_{a,b,r}$, possesses many interesting statistical properties, and it reduces to the entropy S_q$S_q$ for b=0$b=0$, a=r:=1−q$a=r:=1-q$ (hence Boltzmann–Gibbs entropy S_BG$S_{BG}$ for b=0$b=0$, a=r→0$a=r\to 0$). The construction of the entropy S_a,b,r$S_{a,b,r}$ is based on a general group-theoretical approach recently proposed by one of us, Tempesta (2016). Indeed, essentially all the properties of this new entropy are obtained as a consequence of the existence of a rational group law, which expresses the structure of S_a,b,r$S_{a,b,r}$ with respect to the composition of statistically independent subsystems. Depending on the choice of the parameters, the entropy S_a,b,r$S_{a,b,r}$ can be used to cover a wide range of physical situations, in which the measure of the accessible phase space increases say exponentially with the number of particles N$N$ of the system, or even stabilizes, by increasing N$N$, to a limiting value.     

Fluxmetric and magnetometric demagnetizing factors for cylinders

Fluxmetric and magnetometric demagnetizing factors, _Nf$N_{\mathrm{f}}$ and _Nm$N_{\mathrm{m}}$, for cylinders along the axial direction are numerically calculated as functions of material susceptibility χ$\chi$ and the ratio γ$\gamma$ of length to diameter. The results have an accuracy better than 0.1% with respect to min(_Nf,m,1-_Nf,m)$\min (N_{\mathrm{f},\mathrm{m}},1-N_{\mathrm{f},\mathrm{m}})$ and are tabulated in the range of 0.01?γ?500$0.01?\gamma ?500$ and -1?χ<∞$-1?\chi <\infty$. _Nm$N_{\mathrm{m}}$ along the radial direction is evaluated with a lower accuracy from _Nm$N_{\mathrm{m}}$ along the axis and tabulated in the range of 0.01?γ?1$0.01?\gamma ?1$ and -1?χ<∞$-1?\chi <\infty$. Some previous results are discussed and several applications are explained based on the new results.     

Rigidity of noncompact complete Bach-flat manifolds

Let (M,g)$(M,g)$ be a noncompact complete Bach-flat manifold with positive Yamabe constant. We prove that (M,g)$(M,g)$ is flat if (M,g)$(M,g)$ has zero scalar curvature and sufficiently small _L2$L_2$ bound of curvature tensor. When (M,g)$(M,g)$ has nonconstant scalar curvature, we prove that (M,g)$(M,g)$ is conformal to the flat space if (M,g)$(M,g)$ has sufficiently small _L2$L_2$ bound of curvature tensor and _L4/3$L_{4/3}$ bound of scalar curvature.     

Singular unitarity in “quantization commutes with reduction”

Let M$M$ be a connected compact quantizable Kähler manifold equipped with a Hamiltonian action of a connected compact Lie group G$G$. Let M//G=?⁻¹(0)/G=_M0$M//G={?}^{-1}\left(0\right)/G=M_0$ be the symplectic quotient at value 0 of the moment map ?$?$. The space _M0$M_0$ may in general not be smooth. It is known that, as vector spaces, there is a natural isomorphism between the quantum Hilbert space over _M0$M_0$ and the G$G$-invariant subspace of the quantum Hilbert space over M$M$. In this paper, without any regularity assumption on the quotient _M0$M_0$, we discuss the relation between the inner products of these two quantum Hilbert spaces under the above natural isomorphism; we establish asymptotic unitarity to leading order in Planck’s constant of a modified map of the above isomorphism under a “metaplectic correction” of the two quantum Hilbert spaces.     

Reduction of generalized complex structures

We study reduction of generalized complex structures. More precisely, we investigate the following question. Let J$J$ be a generalized complex structure on a manifold M$M$, which admits an action of a Lie group G$G$ preserving J$J$. Assume that _M0$M_0$ is a G$G$-invariant smooth submanifold and the G$G$-action on _M0$M_0$ is proper and free so that _MG?_M0/G$M_G?M_0/G$ is a smooth manifold. Under what condition does J$J$ descend to a generalized complex structure on _MG$M_G$? We describe a sufficient condition for the reduction to hold, which includes the Marsden–Weinstein reduction of symplectic manifolds and the reduction of the complex structures in Kähler manifolds as special cases. As an application, we study reduction of generalized Kähler manifolds.     

Comparative study between a two-dimensional anisotropic Heisenberg antiferromagnet with easy-axis single-ion anisotropy and one with easy-axis exchange anisotropy

Generally, in literature, easy-axis single ion anisotropy and easy-axis exchange anisotropy was treated in indistinct way. In this work we propose to perform a comparative study of the effects of these two easy-axis anisotropies on the behavior of the magnetization and the critical temperature (_Tc)$(T_{\mathrm{c}})$ in the 2D classical Heisenberg antiferromagnetic model. In order to study the low-temperature thermodynamics of this model, we should consider the contribution of anisotropic spin waves, using a self-consistent harmonic approximation (SCHA) theory. We compare the predictions of SCHA with numerical simulations on L×L$L\times L$ square lattices using Monte Carlo (MC) simulations, which include effects due to all thermodynamically allowed excitations. Our SCHA results are in good agreement with our MC simulations results and have shown that the strong K$K$ limit gives two different Ising-like behavior. In the exchange anisotropic case, the dependence of _Tc$T_{\mathrm{c}}$ on anisotropic parameter K$K$ becomes linear and in the single-ion anisotropic case, _Tc$T_{\mathrm{c}}$ becomes independent of K$K$. Also, using MC simulations and finite size scaling, we show that the critical exponents in the two anisotropic case are compatible with the 2D Ising values α=0.125$\alpha =0.125$ and γ=1.75$\gamma =1.75$.     

Properties of reachable sets in the sub-Lorentzian geometry

The aim of this paper is to develop local theory of future timelike, nonspacelike and null reachable sets from a given point _q0$q_0$ in the sub-Lorentzian geometry. In particular, we prove that if U$U$ is a normal neighbourhood of _q0$q_0$ then the three reachable sets, computed relative to U$U$, have identical interiors and boundaries with respect to U$U$. Further, among other things, we show that for Lorentzian metrics on contact distributions on R²ⁿ⁺¹${\mathbb{R}}^{2n+1}$, n≥1$n\ge 1$, the boundary of reachable sets from _q0$q_0$ is, in a neighbourhood of _q0$q_0$, made up of null future directed curves starting from _q0$q_0$. Every such curve has only a finite number of non-smooth points; smooth pieces of every such curve are Hamiltonian geodesics. For general sub-Lorentzian structures, contrary to the Lorentzian case, timelike curves may appear on the boundary. It turns out that such curves are always Goh curves. We also generalize a classical result on null Lorentzian geodesics: every null future directed Hamiltonian sub-Lorentzian geodesic initiating at _q0$q_0$ is contained, at least to a certain moment of time, in the boundary of the reachable set from _q0$q_0$.     

On the parallel transport of the Ricci curvatures

Geometrical characterizations are given for the tensor R⋅S$R\cdot S$, where S$S$ is the Ricci tensor   of a (semi-)Riemannian manifold (M,g)$(M,g)$ and R$R$ denotes the curvature operator   acting on S$S$ as a derivation, and of the Ricci Tachibana tensor  _∧g⋅S${\wedge }_g\cdot S$, where the natural metrical operator  _∧g${\wedge }_g$ also acts as a derivation on S$S$. As a combination, the Ricci curvatures   associated with directions on M$M$, of which the isotropy determines that M$M$ is Einstein, are extended to the Ricci curvatures of Deszcz   associated with directions and planes on M$M$, and of which the isotropy determines that M$M$ is Ricci pseudo-symmetric in the sense of Deszcz.     

Post-breakthrough scaling in reservoir field simulation

We study the oil displacement and production behavior in an isothermal thin layered reservoir model subjected to water flooding. We use the CMG’s (Computer Modelling Group  ) numerical simulators to solve mass balance equations. The influences of the viscosity ratio (m≡_μoil/_μwater$m\equiv {\mu }_{oil}/{\mu }_{water}$) and the inter-well (injector-producer) distance r$r$ on the oil production rate C(t)$C\left(t\right)$ and the breakthrough time _tbr$t_{br}$ are investigated. Two types of reservoir configuration are used, namely one with random porosities and another with a percolation cluster structure. We observe that the breakthrough time follows a power-law of m$m$ and r$r$, _tbr∝r^αm^β$t_{br}\propto r^{\alpha }{m}^{\beta }$, with α=1.8$\alpha =1.8$ and β=−0.25$\beta =-0.25$ for the random porosity type, and α=1.0$\alpha =1.0$ and β=−0.2$\beta =-0.2$ for the percolation cluster type. Moreover, our results indicate that the oil production rate is a power law of time. In the percolation cluster type of reservoir, we observe that P(t)∝t^γ$P\left(t\right)\propto t^{\gamma }$, with γ=−1.81$\gamma =-1.81$, where P(t)$P\left(t\right)$ is the time derivative of C(t)$C\left(t\right)$. The curves related to different values of m$m$ and r$r$ may be collapsed suggesting a universal behavior for the oil production rate.