Network dynamics in the healthy and epileptic developing brain

Electroencephalography (EEG) allows recording of cortical activity at high temporal resolution. EEG recordings can be summarised along different dimensions using network-level quantitative measures, e.g. channel-to-channel correlation, or band power distributions across channels. These reveal network patterns that unfold over a range of different time scales and can be tracked dynamically. Here we describe the dynamics of network-state transitions in EEG recordings of spontaneous brain activity in normally developing infants and infants with severe early infantile epileptic encephalopathies (n=8, age: 1-8 months). We describe differences in measures of EEG dynamics derived from band power, and correlation-based summaries of network-wide brain activity. We further show that EEGs from different patient groups and controls may be distinguishable based on a small set of the novel quantitative measures introduced here, which describe dynamic network state switching. Quantitative measures related to the sharpness of switching from one correlation pattern to another show the largest differences between groups. These findings reveal that the early epileptic encephalopathies are associated with characteristic dynamic features at the network level. Quantitative network-based analyses like the one presented here may in future inform the clinical use of quantitative EEG for diagnosis

Extracting the transition network of epileptic seizure onset

In presurgical monitoring, focal seizure onset is visually assessed from intracranial electroencephalogram (EEG), typically based on the selection of channels that show the strongest changes in amplitude and frequency. As epileptic seizure dynamics is increasingly considered to reflect changes in potentially distributed neural networks, it becomes important to also assess the interrelationships between channels. We propose a workflow to quantitatively extract the nodes and edges contributing to the seizure onset using an across-seizure scoring. We propose a quantification of the consistency of EEG channel contributions to seizure onset within a patient. The workflow is exemplified using recordings from patients with different degrees of seizure-onset consistency. We propose a data-driven analysis method to investigate the transitions to focal-onset seizures from invasive recordings. We employ a combination of uni- and bivariate quantification and a quality score to extract the channels and channel pairs that optimally display the transition in invasive recordings. This allows us to identify the consistency of anatomical contributions involved in seizure onset across multiple seizures of a patient. We hope our approach leads us closer to a better understanding of the complex transition dynamics in refractory epilepsy

The Evolution of Universe with th B-I Type Phantom Scalar Field

We considered the phantom cosmology with a lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1+\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$, which is original from the nonlinear Born-Infeld type scalar field with the lagrangian $\displaystyle L=\frac{1}{\eta}[1-\sqrt{1-\eta g^{\mu\nu}\phi_{, \mu}\phi_{, \nu}}]-u(\phi)$. This cosmological model can explain the accelerated expansion of the universe with the equation of state parameter $w\leq-1$. We get a sufficient condition for a arbitrary potential to admit a late time attractor solution: the value of potential $u(X_c)$ at the critical point $(X_c,0)$ should be maximum and large than zero. We study a specific potential with the form of $u(\phi)=V_0(1+\frac{\phi}{\phi_0})e^{(-\frac{\phi}{\phi_0})}$ via phase plane analysis and compute the cosmological evolution by numerical analysis in detail. The result shows that the phantom field survive till today (to account for the observed late time accelerated expansion) without interfering with the nucleosynthesis of the standard model(the density parameter $\Omega_{\phi}\simeq10^{-12}$ at the equipartition epoch), and also avoid the future collapse of the universe.Comment: 17 pages, 10 figures,typos corrected, references added,figures added and enriched, title changed, main result remaine

Generalized Sagnac Effect

Experiments were conducted to study light propagation in a light waveguide loop consisting of linearly and circularly moving segments. We found that any segment of the loop contributes to the total phase difference between two counterpropagating light beams in the loop. The contribution is proportional to a product of the moving velocity v and the projection of the segment length Deltal on the moving direction, Deltaphi=4pivDeltal/clambda. It is independent of the type of motion and the refractive index of waveguides. The finding includes the Sagnac effect of rotation as a special case and suggests a new fiber optic sensor for measuring linear motion with nanoscale sensitivity.Comment: 3 pages (including 3 figures

Density of states near the Mott-Hubbard transition in the limit of large dimensions

The zero temperature Mott-Hubbard transition as a function of the Coulomb repulsion U is investigated in the limit of large dimensions. The behavior of the density of states near the transition at U=U_c is analyzed in all orders of the skeleton expansion. It is shown that only two transition scenarios are consistent with the skeleton expansion for U<U_c: (i) The Mott-Hubbard transition is "discontinuous" in the sense that in the density of states finite spectral weight is redistributed at U_c. (ii) The transition occurs via a point at U=U_c where the system is neither a Fermi liquid nor an insulator.Comment: 4 pages, 1 figure; revised version accepted for publication in Phys. Rev. Let

A computational study of stimulus driven epileptic seizure abatement

This is the final version of the article. Available from Public Library of Science via the DOI in this record.Active brain stimulation to abate epileptic seizures has shown mixed success. In spike-wave (SW) seizures, where the seizure and background state were proposed to coexist, single-pulse stimulations have been suggested to be able to terminate the seizure prematurely. However, several factors can impact success in such a bistable setting. The factors contributing to this have not been fully investigated on a theoretical and mechanistic basis. Our aim is to elucidate mechanisms that influence the success of single-pulse stimulation in noise-induced SW seizures. In this work, we study a neural population model of SW seizures that allows the reconstruction of the basin of attraction of the background activity as a four dimensional geometric object. For the deterministic (noise-free) case, we show how the success of response to stimuli depends on the amplitude and phase of the SW cycle, in addition to the direction of the stimulus in state space. In the case of spontaneous noise-induced seizures, the basin becomes probabilistic introducing some degree of uncertainty to the stimulation outcome while maintaining qualitative features of the noise-free case. Additionally, due to the different time scales involved in SW generation, there is substantial variation between SW cycles, implying that there may not be a fixed set of optimal stimulation parameters for SW seizures. In contrast, the model suggests an adaptive approach to find optimal stimulation parameters patient-specifically, based on real-time estimation of the position in state space. We discuss how the modelling work can be exploited to rationally design a successful stimulation protocol for the abatement of SW seizures using real-time SW detection.This work was supported by the EPSRC (EP/K026992/1), the BBSRC, the DTC for Systems Biology (University of Manchester), and the Nanyang Technological University Singapore. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

A large Hilbert space QRPA and RQRPA calculation of neutrinoless double beta decay

A large Hilbert space is used for the calculation of the nuclear matrix elements governing the light neutrino mass mediated mode of neutrinoless double beta decay of Ge76, Mo100, Cd116, Te128 and Xe136 within the proton-neutron quasiparticle random phase approximation (pn-QRPA) and the renormalized QRPA with proton-neutron pairing (full-RQRPA) methods. We have found that the nuclear matrix elements obtained with the standard pn-QRPA for several nuclear transitions are extremely sensitive to the renormalization of the particle-particle component of the residual interaction of the nuclear hamiltonian. Therefore the standard pn-QRPA does not guarantee the necessary accuracy to allow us to extract a reliable limit on the effective neutrino mass. This behaviour, already known from the calculation of the two-neutrino double beta decay matrix elements, manifests itself in the neutrinoless double-beta decay but only if a large model space is used. The full-RQRPA, which takes into account proton-neutron pairing and considers the Pauli principle in an approximate way, offers a stable solution in the physically acceptable region of the particle-particle strength. In this way more accurate values on the effective neutrino mass have been deduced from the experimental lower limits of the half-lifes of neutrinoless double beta decay.Comment: 19 pages, RevTex, 1 Postscript figur

Quantum Monte Carlo calculation of the finite temperature Mott-Hubbard transition

We present clear numerical evidence for the coexistence of metallic and insulating dynamical mean field theory(DMFT) solutions in a half-filled single-band Hubbard model with bare semicircular density of states at finite temperatures. Quantum Monte Carlo(QMC) method is used to solve the DMFT equations. We discuss important technical aspects of the DMFT-QMC which need to be taken into account in order to obtain the reliable results near the coexistence region. Among them are the critical slowing down of the iterative solutions near phase boundaries, the convergence criteria for the DMFT iterations, the interpolation of the discretized Green's function and the reduction of QMC statistical and systematic errors. Comparison of our results with those of other numerical methods is presented in a phase diagram.Comment: 4 pages, 5 figure

Local-Ansatz Approach with Momentum Dependent Variational Parameters to Correlated Electron Systems

A new wavefunction which improves the Gutzwiller-type local ansatz method has been proposed to describe the correlated electron system. The ground-state energy, double occupation number, momentum distribution function, and quasiparticle weight have been calculated for the half-filled band Hubbard model in infinite dimensions. It is shown that the new wavefunction improves the local-ansatz approach (LA) proposed by Stollhoff and Fulde. Especially, calculated momentum distribution functions show a reasonable momentum dependence. The result qualitatively differs from those obtained by the LA and the Gutzwiller wavefunction. Furthermore, the present approach combined with the projection operator method CPA is shown to describe quantitatively the excitation spectra in the insulator regime as well as the critical Coulomb interactions for a gap formation in infinite dimensions.Comment: To be published in Phys. Soc. Jpn. 77 No.11 (2008

Mott-Hubbard transition in infinite dimensions

We calculate the zero-temperature gap and quasiparticle weight of the half-filled Hubbard model with a random dispersion relation. After extrapolation to the thermodynamic limit, we obtain reliable bounds on these quantities for the Hubbard model in infinite dimensions. Our data indicate that the Mott-Hubbard transition is continuous, i.e., that the quasiparticle weight becomes zero at the same critical interaction strength at which the gap opens.Comment: 4 pages, RevTeX, 5 figures included with epsfig Final version for PRL, includes L=14 dat