Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.

Can we observe collective neuronal activity from macroscopic aggregate signals?

The fundamental problem faced by noninvasive neuroimaging techniques such as EEG/MEG1 is to elucidate functionally important aspects of the microscopic neuronal network dynamics from macroscopic aggregate measurements. Due to the mixing of the activities of large neuronal populations in the observed macroscopic aggregate, recovering the underlying network that generates the signal in the absence of any additional information represents a considerable challenge. Recent MEG studies have shown that macroscopic measurements contain sufficient information to allow the differentiation between patterns of activity, which are likely to represent different stimulus-specific collective modes in the underlying network (Hadjipapas, A., Adjamian, P., Swettenham, J.B., Holliday, I.E., Barnes, G.R., 2007. Stimuli of varying spatial scale induce gamma activity with distinct temporal characteristics in human visual cortex. NeuroImage 35, 518–530). The next question arising in this context is whether aspects of collective network activity can be recovered from a macroscopic aggregate signal. We propose that this issue is most appropriately addressed if MEG/EEG signals are to be viewed as macroscopic aggregates arising from networks of coupled systems as opposed to aggregates across a mass of largely independent neural systems. We show that collective modes arising in a network of simulated coupled systems can be indeed recovered from the macroscopic aggregate. Moreover, we show that nonlinear state space methods yield a good approximation of the number of effective degrees of freedom in the network. Importantly, information about hidden variables, which do not directly contribute to the aggregate signal, can also be recovered. Finally, this theoretical framework can be applied to experimental MEG/EEG data in the future, enabling the inference of state dependent changes in the degree of local synchrony in the underlying network.

Sensory and cortical activation of distinct glial cell subtypes in the somatosensory thalamus of young rats

The rodent ventrobasal (VB) thalamus receives sensory inputs from the whiskers and projects to the cortex, from which it receives reciprocal excitatory afferents. Much is known about the properties and functional roles of these glutamatergic inputs to thalamocortical neurons in the VB, but no data are available on how these afferents can affect thalamic glial cells. In this study, we used combined electrophysiological recordings and intracellular calcium ([Ca(2+)](i)) imaging to investigate glial cell responses to synaptic afferent stimulation. VB thalamus glial cells can be divided into two groups based on their [Ca(2+)](i) and electrophysiological responses to sensory and corticothalamic stimulation. One group consists of astrocytes, which stain positively for S100B and preferentially load with SR101, have linear current-voltage relations and low input resistance, show no voltage-dependent [Ca(2+)](i) responses, but express mGluR5-dependent [Ca(2+)](i) transients following stimulation of the sensory and/or corticothalamic excitatory afferent pathways. Cells of the other glial group, by contrast, stain positively for NG2, and are characterized by high input resistance, the presence of voltage-dependent [Ca(2+)](i) elevations and voltage-gated inward currents. There were no synaptically induced [Ca(2+)](i) elevations in these cells under control conditions. These results show that thalamic glial cell responses to synaptic input exhibit different properties to those of thalamocortical neurons. As VB astrocytes can respond to synaptic stimulation and signal to neighbouring neurons, this glial cell organization may have functional implications for the processing of somatosensory information and modulation of behavioural state-dependent thalamocortical network activities.

Study of Queuing Systems with a Generalized Departure Process

This work was supported by the Bulgarian National Science Fund under grant BY-TH-105/2005.

Limit Theorems for Regenerative Excursion Processes

This work is supported by Bulgarian NFSI, grant No. MM–704/97

Branching Processes in Autoregressive Random Environment

2000 Mathematics Subject Classification: 60J80, 60K05.

Selective Coupling between Theta Phase and Neocortical Fast Gamma Oscillations during REM-Sleep in Mice

TORT, A. B. L. ; SCHEFFER-TEIXEIRA, R ; Souza, B.C. ; DRAGUHN, A. ; BRANKACK, J. . Theta-associated high-frequency oscillations (110-160 Hz) in the hippocampus and neocortex. Progress in Neurobiology , v. 100, p. 1-14, 2013.

Bumps, breathers, and waves in a neural network with spike frequency adaptation

In this Letter we introduce a continuum model of neural tissue that include the effects of so-called spike frequency adaptation (SFA). The basic model is an integral equation for synaptic activity that depends upon the non-local network connectivity, synaptic response, and firing rate of a single neuron. A phenomenological model of SFA is examined whereby the firing rate is taken to be a simple state-dependent threshold function. As in the case without SFA classical Mexican-Hat connectivity is shown to allow for the existence of spatially localized states (bumps). Importantly an analysis of bump stability using recent Evans function techniques shows that bumps may undergo instabilities leading to the emergence of both breathers and traveling waves. Moreover, a similar analysis for traveling pulses leads to the conditions necessary to observe a stable traveling breather. Direct numerical simulations both confirm our theoretical predictions and illustrate the rich dynamic behavior of this model, including the appearance of self-replicating bumps.

Selective Coupling between Theta Phase and Neocortical Fast Gamma Oscillations during REM-Sleep in Mice

TORT, A. B. L. ; SCHEFFER-TEIXEIRA, R ; Souza, B.C. ; DRAGUHN, A. ; BRANKACK, J. . Theta-associated high-frequency oscillations (110-160 Hz) in the hippocampus and neocortex. Progress in Neurobiology , v. 100, p. 1-14, 2013.

Exotic dynamics in a firing rate model of neural tissue with threshold accommodation

Many of the equations describing the dynamics of neural systems are written in terms of firing rate functions, which themselves are often taken to be threshold functions of synaptic activity. Dating back to work by Hill in 1936 it has been recognized that more realistic models of neural tissue can be obtained with the introduction of state-dependent dynamic thresholds. In this paper we treat a specific phenomenological model of threshold accommodation that mimics many of the properties originally described by Hill. Importantly we explore the consequences of this dynamic threshold at the tissue level, by modifying a standard neural field model of Wilson-Cowan type. As in the case without threshold accommodation classical Mexican-Hat connectivity is shown to allow for the existence of spatially localized states (bumps) in both one and two dimensions. Importantly an analysis of bump stability in one dimension, using recent Evans function techniques, shows that bumps may undergo instabilities leading to the emergence of both breathers and traveling waves. Moreover, a similar analysis for traveling pulses leads to the conditions necessary to observe a stable traveling breather. In the regime where a bump solution does not exist direct numerical simulations show the possibility of self-replicating bumps via a form of bump splitting. Simulations in two space dimensions show analogous localized and traveling solutions to those seen in one dimension. Indeed dynamical behavior in this neural model appears reminiscent of that seen in other dissipative systems that support localized structures, and in particular those of coupled cubic complex Ginzburg-Landau equations. Further numerical explorations illustrate that the traveling pulses in this model exhibit particle like properties, similar to those of dispersive solitons observed in some three component reaction-diffusion systems. A preliminary account of this work first appeared in S Coombes and M R Owen, Bumps, breathers, and waves in a neural network with spike frequency adaptation, Physical Review Letters 94 (2005), 148102(1-4).

As Santas Casas da Misericórdia na República Brasileira 1922-1945

O principal objectivo desta dissertação é o de conhecer um pouco melhor o processo de implementação das Santas Casas de Misericórdia no Brasil, dando especial ênfase à sua expansão durante o período de consolidação da República brasileira, mais concretamente entre 1922 a 1945. A necessária contextualização levou a pesquisa sobre as Misericórdias até à fase colonial e imperial do Brasil, acabando por demonstrar que as mesmas se fortaleceram no segmento de assistência médica, durante o período em análise, tomando o Estado brasileiro dependente das suas actividades. Este trabalho discute ainda o imaginário social da caridade e filantropia e a forma como tais preceitos configuraram a assistência médico-social no país. ABSTRACT; This dissertation aims to better know the implementation process of the Santas Casas de Misericórdia in Brazil, highlighting their expansion during the Republic, mainly between 1922 and 1945. For a better historical contextualization the study explores the Brazil's colonial and imperial phases, demonstrating that the Misericórdias progressively strength their power in the medical assistance segment, becoming the State dependent of their activities. The dissertation also discusses the philanthropy and charity's collective social imagery, as well as the way in which such concepts shaped the medico-social assistance in the country.

Engineering a Quantum Many-body Hamiltonian with Trapped Ions

While fault-tolerant quantum computation might still be years away, analog quantum simulators offer a way to leverage current quantum technologies to study classically intractable quantum systems. Cutting edge quantum simulators such as those utilizing ultracold atoms are beginning to study physics which surpass what is classically tractable. As the system sizes of these quantum simulators increase, there are also concurrent gains in the complexity and types of Hamiltonians which can be simulated. In this work, I describe advances toward the realization of an adaptable, tunable quantum simulator capable of surpassing classical computation. We simulate long-ranged Ising and XY spin models which can have global arbitrary transverse and longitudinal fields in addition to individual transverse fields using a linear chain of up to 24 Yb+ 171 ions confined in a linear rf Paul trap. Each qubit is encoded in the ground state hyperfine levels of an ion. Spin-spin interactions are engineered by the application of spin-dependent forces from laser fields, coupling spin to motion. Each spin can be read independently using state-dependent fluorescence. The results here add yet more tools to an ever growing quantum simulation toolbox. One of many challenges has been the coherent manipulation of individual qubits. By using a surprisingly large fourth-order Stark shifts in a clock-state qubit, we demonstrate an ability to individually manipulate spins and apply independent Hamiltonian terms, greatly increasing the range of quantum simulations which can be implemented. As quantum systems grow beyond the capability of classical numerics, a constant question is how to verify a quantum simulation. Here, I present measurements which may provide useful metrics for large system sizes and demonstrate them in a system of up to 24 ions during a classically intractable simulation. The observed values are consistent with extremely large entangled states, as much as ~95% of the system entangled. Finally, we use many of these techniques in order to generate a spin Hamiltonian which fails to thermalize during experimental time scales due to a meta-stable state which is often called prethermal. The observed prethermal state is a new form of prethermalization which arises due to long-range interactions and open boundary conditions, even in the thermodynamic limit. This prethermalization is observed in a system of up to 22 spins. We expect that system sizes can be extended up to 30 spins with only minor upgrades to the current apparatus. These results emphasize that as the technology improves, the techniques and tools developed here can potentially be used to perform simulations which will surpass the capability of even the most sophisticated classical techniques, enabling the study of a whole new regime of quantum many-body physics.

Time-Dependent Density Functional Molecular Orbital and Excited State Calculations on Bis(porphyrinyl)butadiynes in the Monocationic, Neutral, Monoanionic, and Dianionic Oxidation States

We report a theoretical study of the multiple oxidation states (1+, 0, 1−, and 2−) of a meso,meso-linked diporphyrin, namely bis[10,15,20-triphenylporphyrinatozinc(II)-5-yl]butadiyne (4), using Time-Dependent Density Functional Theory (TDDFT). The origin of electronic transitions of singlet excited states is discussed in comparison to experimental spectra for the corresponding oxidation states of the close analogue bis{10,15,20-tris[3‘,5‘-di-tert-butylphenyl]porphyrinatozinc(II)-5-yl}butadiyne (3). The latter were measured in previous work under in situ spectroelectrochemical conditions. Excitation energies and orbital compositions of the excited states were obtained for these large delocalized aromatic radicals, which are unique examples of organic mixed-valence systems. The radical cations and anions of butadiyne-bridged diporphyrins such as 3 display characteristic electronic absorption bands in the near-IR region, which have been successfully predicted with use of these computational methods. The radicals are clearly of the “fully delocalized” or Class III type. The key spectral features of the neutral and dianionic states were also reproduced, although due to the large size of these molecules, quantitative agreement of energies with observations is not as good in the blue end of the visible region. The TDDFT calculations are largely in accord with a previous empirical model for the spectra, which was based simplistically on one-electron transitions among the eight key frontier orbitals of the C4 (1,4-butadiyne) linked diporphyrins.

Diaminotriazine substituted diphenyl ether: reversible structural transformation and solvent dependent solid state fluorescence

The effect of molecular shape and position of hydrogen bonding functionality in the solid state structural self-assembly was investigated using diaminotriazine substituted diphenyl ether based positional isomers (1-5). The molecular shape was modulated by changing diaminotriazine position that produced channel supramolecular structures in 1, 3 and 5. There exists a direct correlation between the molecular shape and three dimensional structures; more linear molecules resulted in close-packing whereas molecules with a labyrinthine topology formed a channel structure. Supramolecular aspects pertaining to the influence of solvent of crystallization in structure formation and reversible structural transformation in solid state were also explored. 1-5 exhibited tunable solid state fluorescence (lambda(max) = 437-496 nm) depending on the diaminotriazine substitutional position and 3 showed solvent-dependent solid state fluorescence. The present study describes the generation of a supramolecular channel structure with functional properties such as tunable fluorescence by varying the position of hydrogen bond functionality and solvent of crystallization.

RNA-dependent RNA polymerase of Japanese encephalitis virus binds the initiator nucleotide GTP to form a mechanistically important pre-initiation state

Flaviviral RNA-dependent RNA polymerases (RdRps) initiate replication of the single-stranded RNA genome in the absence of a primer. The template sequence 5'-CU-3' at the 3'-end of the flaviviral genome is highly conserved. Surprisingly, flaviviral RdRps require high concentrations of the second incoming nucleotide GTP to catalyze de novo template-dependent RNA synthesis. We show that GTP stimulates de novo RNA synthesis by RdRp from Japanese encephalitis virus (jRdRp) also. Crystal structures of jRdRp complexed with GTP and ATP provide a basis for specific recognition of GTP. Comparison of the jRdRp(GTP) structure with other viral RdRp-GTP structures shows that GTP binds jRdRp in a novel conformation. Apo-jRdRp structure suggests that the conserved motif F of jRdRp occupies multiple conformations in absence of GTP. Motif F becomes ordered on GTP binding and occludes the nucleotide triphosphate entry tunnel. Mutational analysis of key residues that interact with GTP evinces that the jRdRp(GTP) structure represents a novel pre-initiation state. Also, binding studies show that GTP binding reduces affinity of RdRp for RNA, but the presence of the catalytic Mn2+ ion abolishes this inhibition. Collectively, these observations suggest that the observed pre-initiation state may serve as a check-point to prevent erroneous template-independent RNA synthesis by jRdRp during initiation.