Synchronization processes in nonlinear systems and their relation to cortical oscillatory dynamics

This thesis was focused on theoretical models of synchronization to cortical dynamics as measured by magnetoencephalography (MEG). Dynamical systems theory was used in both identifying relevant variables for brain coordination and also in devising methods for their quantification. We presented a method for studying interactions of linear and chaotic neuronal sources using MEG beamforming techniques. We showed that such sources can be accurately reconstructed in terms of their location, temporal dynamics and possible interactions. Synchronization in low-dimensional nonlinear systems was studied to explore specific correlates of functional integration and segregation. In the case of interacting dissimilar systems, relevant coordination phenomena involved generalized and phase synchronization, which were often intermittent. Spatially-extended systems were then studied. For locally-coupled dissimilar systems, as in the case of cortical columns, clustering behaviour occurred. Synchronized clusters emerged at different frequencies and their boundaries were marked through oscillation death. The macroscopic mean field revealed sharp spectral peaks at the frequencies of the clusters and broader spectral drops at their boundaries. These results question existing models of Event Related Synchronization and Desynchronization. We re-examined the concept of the steady-state evoked response following an AM stimulus. We showed that very little variability in the AM following response could be accounted by system noise. We presented a methodology for detecting local and global nonlinear interactions from MEG data in order to account for residual variability. We found crosshemispheric nonlinear interactions of ongoing cortical rhythms concurrent with the stimulus and interactions of these rhythms with the following AM responses. Finally, we hypothesized that holistic spatial stimuli would be accompanied by the emergence of clusters in primary visual cortex resulting in frequency-specific MEG oscillations. Indeed, we found different frequency distributions in induced gamma oscillations for different spatial stimuli, which was suggestive of temporal coding of these spatial stimuli. Further, we addressed the bursting character of these oscillations, which was suggestive of intermittent nonlinear dynamics. However, we did not observe the characteristic-3/2 power-law scaling in the distribution of interburst intervals. Further, this distribution was only seldom significantly different to the one obtained in surrogate data, where nonlinear structure was destroyed. In conclusion, the work presented in this thesis suggests that advances in dynamical systems theory in conjunction with developments in magnetoencephalography may facilitate a mapping between levels of description int he brain. this may potentially represent a major advancement in neuroscience.

Circannual rhythms of maturation in the female rainbow trout and their entrainment by changes in photoperiod

The importance of endogenous rhythms in the photoperiodic control of the annual reproduction cycle in female rainbow trout was investigated. The effect of photoperiod regimes on the different stages of maturation was assessed by recording the timing of ovulation and from quantifying associated changes in serum oestradiol-17,testosterone and total calcium. Maintained under constant 6L:18D and constant temperature for up to four years, rainbow trout exhibited an endogenous rhythm of maturation with a periodicity of approximately one year. This rhythm of maturation appears to be driven by an autonomous circannual oscillator or clock which can be dissociated from the neuroendocrine mechanisms controlling gonadal maturation. Under conditions of constant 18L:6D or LL the periodicity of the maturation rhythm was 5.5-6 months; it is suggested that this periodicity may be caused by a splitting or uncoupling of at least two circannual clocks involved in the control of maturation. Abrupt changes in the length of the photoperiod act as a zeitgeber to entrain the endogenous rhythm of maturation. Whether the timing of maturation is advanced or delayed depends primarily on the direction of the change in photoperiod and its timing in relation to the phase of the rhythm, with the magnitude of the alteration in photoperiod having only a supplementary effect. The effect of specific changes in photoperiod on the entrainment of the maturation cycle can be described in terms of a phase-response curve. Photic information is transduced, probably by the pineal gland, into a daily rhythm of melatonin; exposure of rainbow trout to skeleton and resonance photoperiod regimes indicated that daylength measurement is effected by endogenous circadian clock(s) rather than by hour-glass mechanisms. A gating mechanism is closely associated with the circannual clock which determines the timing of onset of maturation in virgin female rainbow trout, only allowing fish that have attained a threshold stage of development to undergo gonadal maturation. Collectively the results support the hypothesis that the female rainbow trout exhibits an endogenous circannual rhythm of maturation which can be entrained by changes in photoperiod.

Atomic structure of the two intermediate phase glasses SiSe4 and GeSe4

The microscopic origin of the intermediate phase in two prototypical covalently bonded AxB1-x network glass forming systems, where A=Ge or Si, B=Se, and 0=x=1, was investigated by combining neutron diffraction with first-principles molecular-dynamics methods. Specifically, the structure of glassy GeSe4 and SiSe4 was examined, and the calculated total structure factor and total pair-correlation function for both materials are in good agreement with experiment. The structure of both glasses differs markedly from a simple model comprising undefective AB4 corner-sharing tetrahedra in which all A atoms are linked by B2 dimers. Instead, edge-sharing tetrahedra occur and the twofold coordinated Se atoms form three distinct structural motifs, namely, Se-Se2, Se-SeGe (or Se-SeSi), and Se-Ge2 (or Se-Si2). This identifies several of the conformations that are responsible for the structural variability in GexSe1-x and SixSe1-x glasses, a quantity that is linked to the finite width of the intermediate phase window.

Intracavity dynamics in high-power mode-locked fiber lasers

A theoretical model is developed which characterizes the intracavity pulse evolutions in high-power fiber lasers. It is shown that experimentally observed dynamics of the key pulse parameters can be described by a reduced model of ordinary differential equations. Critical in driving the intracavity dynamics is the amplitude and phase modulations generated by the discrete elements in the laser. The theory gives a simple geometrical description of the intracavity dynamics and possible operation modes of the laser cavity. Furthermore, it provides a simple and efficient method for optimizing the performance of complex multiparametric laser systems.

Cortical dynamics and synchronization related to multiple target consolidation under rapid-serial-visual-presentation conditions

The present report reviews behavioural, electroencephalographic, and especially magnetoencephalographic findings on the cortical mechanisms underlying attentional processes that separate targets from distractors and that ensure durable target representations for goal-directed action. A common way of investigation is to observe the system’s overt and covert behaviour when capacity limitations are reached. Here we focus on the aspect of temporally enhanced processing load, namely on performance deficits occurring under rapid-serial-visual-presentation (RSVP) conditions. The most prominent of these deficits is the so-called “attentional blink” (AB) effect. We first report MEG findings with respect to the time course of activation that shows modulations around 300 ms after target onset which reflect demands and success of target consolidation. Then, findings regarding long-range inter-area phase synchronization are reported that are hypothesized to mediate communication within the attentional network. Changes in synchronization reflect changes in the attentional demands of the task and are directly related to behavioural performance. Furthermore, enhanced vigilance of the system elicits systematically increased synchronization indices. A hypothetical framework is sketched out that aims at explaining limitations in multiple target consolidation under RSVP conditions.

Quantifying the complexity of chaos in multibasin multidimensional dynamics of molecular systems

The simulated classical dynamics of a small molecule exhibiting self-organizing behavior via a fast transition between two states is analyzed by calculation of the statistical complexity of the system. It is shown that the complexity of molecular descriptors such as atom coordinates and dihedral angles have different values before and after the transition. This provides a new tool to identify metastable states during molecular self-organization. The highly concerted collective motion of the molecule is revealed. Low-dimensional subspaces dynamics is found sensitive to the processes in the whole, high-dimensional phase space of the system. © 2004 Wiley Periodicals, Inc.

Dynamic complexity of chaotic transitions in high-dimensional classical dynamics:Leu-Enkephalin folding

Leu-Enkephalin in explicit water is simulated using classical molecular dynamics. A ß-turn transition is investigated by calculating the topological complexity (in the "computational mechanics" framework [J. P. Crutchfield and K. Young, Phys. Rev. Lett., 63, 105 (1989)]) of the dynamics of both the peptide and the neighbouring water molecules. The complexity of the atomic trajectories of the (relatively short) simulations used in this study reflect the degree of phase space mixing in the system. It is demonstrated that the dynamic complexity of the hydrogen atoms of the peptide and almost all of the hydrogens of the neighbouring waters exhibit a minimum precisely at the moment of the ß-turn transition. This indicates the appearance of simplified periodic patterns in the atomic motion, which could correspond to high-dimensional tori in the phase space. It is hypothesized that this behaviour is the manifestation of the effect described in the approach to molecular transitions by Komatsuzaki and Berry [T. Komatsuzaki and R.S. Berry, Adv. Chem. Phys., 123, 79 (2002)], where a "quasi-regular" dynamics at the transition is suggested. Therefore, for the first time, the less chaotic character of the folding transition in a realistic molecular system is demonstrated. © Springer-Verlag Berlin Heidelberg 2006.

Bridging large and small scales of water models using hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework

This thesis presents a two-dimensional water model investigation and development of a multiscale method for the modelling of large systems, such as virus in water or peptide immersed in the solvent. We have implemented a two-dimensional ‘Mercedes Benz’ (MB) or BN2D water model using Molecular Dynamics. We have studied its dynamical and structural properties dependence on the model’s parameters. For the first time we derived formulas to calculate thermodynamic properties of the MB model in the microcanonical (NVE) ensemble. We also derived equations of motion in the isothermal–isobaric (NPT) ensemble. We have analysed the rotational degree of freedom of the model in both ensembles. We have developed and implemented a self-consistent multiscale method, which is able to communicate micro- and macro- scales. This multiscale method assumes, that matter consists of the two phases. One phase is related to micro- and the other to macroscale. We simulate the macro scale using Landau Lifshitz-Fluctuating Hydrodynamics, while we describe the microscale using Molecular Dynamics. We have demonstrated that the communication between the disparate scales is possible without introduction of fictitious interface or approximations which reduce the accuracy of the information exchange between the scales. We have investigated control parameters, which were introduced to control the contribution of each phases to the matter behaviour. We have shown, that microscales inherit dynamical properties of the macroscales and vice versa, depending on the concentration of each phase. We have shown, that Radial Distribution Function is not altered and velocity autocorrelation functions are gradually transformed, from Molecular Dynamics to Fluctuating Hydrodynamics description, when phase balance is changed. In this work we test our multiscale method for the liquid argon, BN2D and SPC/E water models. For the SPC/E water model we investigate microscale fluctuations which are computed using advanced mapping technique of the small scales to the large scales, which was developed by Voulgarakisand et. al.

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

The structure and dynamics of saxicolous lichen communities

The saxicolous lichen vegetation on Ordovician slate rock at the mouth of the River Dovey, South Merionethshire, Wales was described in relation to several environmental variables which include aspect, slope angle, light intensity, rock porosity, rock microtopography and rock stability. Each of the measured environmental variables was shown to influence the lichen vegetation. A number of groups of species which were characteristic of certain environments were described. The data from the saxicolous lichen communities were analysed using multivariate analysis. Qualitative and quantitative data were ordinated, the qualitative data being easier to interpret ecologically, and site number (which reflects distance from the sea and altitude), rock porosity and light intensity were shown to be important environmental variables. A classification of the data was also carried out. The results of the ordination and classification were combined together and a model constructed which describes saxicolous lichen vegetation. A method which uses the model as an aid to the design and interpretation of field experiments is described. The model is applied to an experiment which investigates the effect on growth of transplanting four saxicolous lichens to different aspects. Growth was inhibited in Physcia orbicularis and Parmelia conspersa on rock surfaces of northwest aspect compared with growth on rock surfaces of southeast aspect. Growth was inhibited in Parmelia glabratula ssp. fuliginosa on rock surfaces of southeast aspect compared with rock surfaces of northwesr aspect. The growth of Parmelia saxatilis was similar at both southeast and northwesr aspects. Growth inhibition or stimulation in thalli of Physcia orbicularis, Parmelia conspersa and Parmelia glabratula ssp. fuliginosa after transplantation was consistent with the predictions of the model while the results for Parmelia saxatilis were not as expected. There was evidence that the frequency of Parmelia conspersa and Parmelia glabratula at a site is related to an effect of the environment on the growth of the thalli. There was also evidence that the frequency of Physcia orbicularis at a site is related to an effect of the environment on the establishment phase of the thalli and for the competitive exclusion of Parmelia saxatilis thalli from southeast facing rock surfaces. The distribution of lichens in relation to height on nine rock surfaces was investigated. It was suggested that the distribution of the lichens was influenced by microclimatic factors which are related to height on the rock, environmental variables which are associated with the rock substratum (e.g. rock porosity and rock microtopography) and by historical factors. The pattern of one crustose and one foliose lichen on four rock surfaces of different aspect and slope was investigated. On the vertically inclined surface the density of small thalli of Buellia aethalea and Parmelia glabratula ssp fuliginosa was correlated with the microtopography of the surface in transects horizontally across the rock surface but not in transects vertically down the rock surface. there were consitent differences in the scale and intensity of pattern horizontally and vertically and also a decrease in the intensity of pattern vertically as the slope of the rock surface decreased. These results were consistent with the suggestion of a gradient of microclimatic factors up the rock. The differences in the scale and intensity of pattern in different size classes in the population were consistent with the changes in pattern with time which have been shown to occur during succession in sand dune and salt marsh vegetation. The relationship between thallus size and height on a rock surface and between the radial growth rate and location of a thallus on a rock surface were investigated. Thalli of Parmelia glabratula ssp. fuliginosa were larger at the top of the rock surface than at the bottom and the data were consistent with the suggestion that the colonisation of the rock surface began at the top and, in time, spread downwards. The radial growth rate of the thalli could not be related to variation in slope, porosity, microtopography or directly to height on the rock but could be related to the horizontal location of the thalli on the rock. These results were consistent with the suggestion that here is a gradient of microclimatic factors across the rock surface which is also modified by height on the rock surface. The succession of lichen communities was described by relating the vegetation to rock porosity, rock microtopography, species diversity and rock stability. An initial stage dominated by crustose lichens leads to communities dominated by crustose, foliose and fruticose species. In the late stages of the succession on some rock surfaces crustose species again become dominant. The occurrence of the climax state and cyclic vegetation change in lichen communities are discussed. A mthod of estimating the age structure of a lichen population by relating thallus size to growth rate is described. The sources of error in the method are discussed in some detail and several refinements suggested to increase the accuracy of the method. The population dynamics of Parmelia glabratula ssp. fuliginosa was investigated by applying life tables to the age structures of eight different populations. The data were consistent with a period of relatively constant recruitment of thalli into the populations. Mortality in lichen populations was divided into deaths which occur after fragmentation of the thallus and deaths which occur after catastrophic environmental events. THe data suggest that the rate of fragmenting death is dependent on the age of the thallus while the rate of catastrophic death is dependent on the number of thalli established in an age class. A comparison of the numbers of thalli in each age class in the eight populations suggested that population density is controlled firstly, by climate and secondly, by variables related to the local rock surface environment. The rate of fragmenting death is related to the diversity of the community and the influence of diversity together with environmental variables in fluctuating or cyclic changes in population number.

Dynamics of cholesteric liquid crystals in the presence of a random magnetic field:stochastic dynamics of cholesteric liquid crystal

Based on dynamic renormalization group techniques, this letter analyzes the effects of external stochastic perturbations on the dynamical properties of cholesteric liquid crystals, studied in presence of a random magnetic field. Our analysis quantifies the nature of the temperature dependence of the dynamics; the results also highlight a hitherto unexplored regime in cholesteric liquid crystal dynamics. We show that stochastic fluctuations drive the system to a second-ordered Kosterlitz-Thouless phase transition point, eventually leading to a Kardar-Parisi-Zhang (KPZ) universality class. The results go beyond quasi-first order mean-field theories, and provides the first theoretical understanding of a KPZ phase in distorted nematic liquid crystal dynamics.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

Real-time high-resolution heterodyne-based measurements of spectral dynamics in fibre lasers

Conventional tools for measurement of laser spectra (e.g. optical spectrum analysers) capture data averaged over a considerable time period. However, the generation spectrum of many laser types may involve spectral dynamics whose relatively fast time scale is determined by their cavity round trip period, calling for instrumentation featuring both high temporal and spectral resolution. Such real-time spectral characterisation becomes particularly challenging if the laser pulses are long, or they have continuous or quasi-continuous wave radiation components. Here we combine optical heterodyning with a technique of spatiooral intensity measurements that allows the characterisation of such complex sources. Fast, round-trip-resolved spectral dynamics of cavity-based systems in real-time are obtained, with temporal resolution of one cavity round trip and frequency resolution defined by its inverse (85 ns and 24 MHz respectively are demonstrated). We also show how under certain conditions for quasi-continuous wave sources, the spectral resolution could be further increased by a factor of 100 by direct extraction of phase information from the heterodyned dynamics or by using double time scales within the spectrogram approach.

Control of the in-cavity dynamics in mode-locked fibre lasers

At the level of fundamental research, fibre lasers provide convenient and reproducible experimental settings for the study of a variety of nonlinear dynamical processes, while at the applied research level, pulses with different and optimised features – e.g., in terms of pulse duration, temporal and/or spectral intensity profile, energy, repetition rate and emission bandwidth – are sought with the general constraint of developing efficient cavity architectures. In this talk, we review our recent progress on the realisation of different regimes of pulse generation in passively mode-locked fibre lasers through control of the in-cavity propagation dynamics. We report on the possibility to achieve both parabolic self-similar and triangular pulse shaping in a mode-locked fibre laser via adjustment of the net normal dispersion and integrated gain of the cavity [1]. We also show that careful control of the gain/loss parameters of a net-normal dispersion laser cavity provides the means of achieving switching among Gaussian pulse, dissipative soliton and similariton pulse solutions in the cavity [2,3]. Furthermore, we report on our recent theoretical and experimental studies of pulse shaping by inclusion of an amplitude and phase spectral filter into the cavity of a laser. We numerically demonstrate that a mode-locked fibre laser can operate in dif- ferent pulse-generation regimes, including parabolic, flattop and triangular waveform generations, depending on the amplitude profile of the in-cavity spectral filter [4]. An application of technique using a flat-top spectral filter is demonstrated to achieve the direct generation of sinc-shaped optical Nyquist pulses of high quality and of a widely tuneable bandwidth from the laser [5]. We also report on a recently-developed versa- tile erbium-doped fibre laser, in which conventional soliton, dispersion-managed soli- ton (stretched-pulse) and dissipative soliton mode-locking regimes can be selectively and reliably targeted by programming different group-velocity dispersion profiles and bandwidths on an in-cavity programmable filter [6]. References: 1. S. Boscolo and S. K. Turitsyn, Phys. Rev. A 85, 043811 (2012). 2. J. Peng et al., Phys. Rev. A 86, 033808 (2012). 3. J. Peng, Opt. Express 24, 3046-3054 (2016). 4. S. Boscolo, C. Finot, H. Karakuzu, and P. Petropoulos, Opt. Lett. 39, 438-441 (2014). 5. S. Boscolo, C. Finot, and S. K. Turitsyn, IEEE Photon. J. 7, 7802008 (2015). 6. J. Peng and S. Boscolo, Sci. Rep. 6, 25995 (2016).

Mode-locking intracavity dynamics:modeling nonlinear power and energy fluctuations

Modern high-power, pulsed lasers are driven by strong intracavity fluctuations. Critical in driving the intracavity dynamics is the nontrivial phase profiles generated and their periodic modification from either nonlinear mode-coupling, spectral filtering or dispersion management. Understanding the theoretical origins of the intracavity fluctuations helps guide the design, optimization and construction of efficient, high-power and high-energy pulsed laser cavities. Three specific mode-locking component are presented for enhancing laser energy: waveguide arrays, spectral filtering and dispersion management. Each component drives a strong intracavity dynamics that is captured through various modeling and analytic techniques.