New path integral simulation algorithms and their application to creep in the quantum sine-Gordon chain

A path integral simulation algorithm which includes a higher-order Trotter approximation (HOA)is analyzed and compared to an approach which includes the correct quantum mechanical pair interaction (effective Propagator (EPr)). It is found that the HOA algorithmconverges to the quantum limit with increasing Trotter number P as P^{-4}, while the EPr algorithm converges as P^{-2}.The convergence rate of the HOA algorithm is analyzed for various physical systemssuch as a harmonic chain,a particle in a double-well potential, gaseous argon, gaseous helium and crystalline argon. A new expression for the estimator for the pair correlation function in the HOA algorithm is derived. A new path integral algorithm, the hybrid algorithm, is developed.It combines an exact treatment of the quadratic part of the Hamiltonian and thehigher-order Trotter expansion techniques.For the discrete quantum sine-Gordon chain (DQSGC), it is shown that this algorithm works more efficiently than all other improved path integral algorithms discussed in this work. The new simulation techniques developed in this work allow the analysis of theDQSGC and disordered model systems in the highly quantum mechanical regime using path integral molecular dynamics (PIMD)and adiabatic centroid path integral molecular dynamics (ACPIMD).The ground state phonon dispersion relation is calculated for the DQSGC by the ACPIMD method.It is found that the excitation gap at zero wave vector is reduced by quantum fluctuations. Two different phases exist: One phase with a finite excitation gap at zero wave vector, and a gapless phase where the excitation gap vanishes.The reaction of the DQSGC to an external driving force is analyzed at T=0.In the gapless phase the system creeps if a small force is applied, and in the phase with a gap the system is pinned. At a critical force, the systems undergo a depinning transition in both phases and flow is induced. The analysis of the DQSGC is extended to models with disordered substrate potentials. Three different cases are analyzed: Disordered substrate potentials with roughness exponent H=0, H=1/2,and a model with disordered bond length. For all models, the ground state phonon dispersion relation is calculated.

Structural and strain analysis in the Late Paleozoic coastal accretionary wedge of central Chile

Abstract In this study structural and finite strain data are used to explore the tectonic evolution and the exhumation history of the Chilean accretionary wedge. The Chilean accretionary wedge is part of a Late Paleozoic subduction complex that developed during subduction of the Pacific plate underneath South America. The wedge is commonly subdivided into a structurally lower Western Series and an upper Eastern Series. This study shows the progressive development of structures and finite strain from the least deformed rocks in the eastern part of the Eastern Series of the accretionary wedge to higher grade schist of the Western Series at the Pacific coast. Furthermore, this study reports finite-strain data to quantify the contribution of vertical ductile shortening to exhumation. Vertical ductile shortening is, together with erosion and normal faulting, a process that can aid the exhumation of high-pressure rocks. In the east, structures are characterized by upright chevron folds of sedimentary layering which are associated with a penetrative axial-plane foliation, S1. As the F1 folds became slightly overturned to the west, S1 was folded about recumbent open F2 folds and an S2 axial-plane foliation developed. Near the contact between the Western and Eastern Series S2 represents a prominent subhorizontal transposition foliation. Towards the structural deepest units in the west the transposition foliation became progressively flat lying. Finite-strain data as obtained by Rf/Phi and PDS analysis in metagreywacke and X-ray texture goniometry in phyllosilicate-rich rocks show a smooth and gradual increase in strain magnitude from east to west. There are no evidences for normal faulting or significant structural breaks across the contact of Eastern and Western Series. The progressive structural and strain evolution between both series can be interpreted to reflect a continuous change in the mode of accretion in the subduction wedge. Before ~320-290 Ma the rocks of the Eastern Series were frontally accreted to the Andean margin. Frontal accretion caused horizontal shortening and upright folds and axial-plane foliations developed. At ~320-290 Ma the mode of accretion changed and the rocks of the Western Series were underplated below the Andean margin. This basal accretion caused a major change in the flow field within the wedge and gave rise to vertical shortening and the development of the penetrative subhorizontal transposition foliation. To estimate the amount that vertical ductile shortening contributed to the exhumation of both units finite strain is measured. The tensor average of absolute finite strain yield Sx=1.24, Sy=0.82 and Sz=0.57 implying an average vertical shortening of ca. 43%, which was compensated by volume loss. The finite strain data of the PDS measurements allow to calculate an average volume loss of 41%. A mass balance approximates that most of the solved material stays in the wedge and is precipitated in quartz veins. The average of relative finite strain is Sx=1.65, Sy=0.89 and Sz=0.59 indicating greater vertical shortening in the structurally deeper units. A simple model which integrates velocity gradients along a vertical flow path with a steady-state wedge is used to estimate the contribution of deformation to ductile thinning of the overburden during exhumation. The results show that vertical ductile shortening contributed 15-20% to exhumation. As no large-scale normal faults have been mapped the remaining 80-85% of exhumation must be due to erosion.

Neural network activity in the neonatal acute slice, slice culture and cell culture

Information processing and storage in the brain may be presented by the oscillations and cell assemblies. Here we address the question of how individual neurons associate together to assemble neural networks and present spontaneous electrical activity. Therefore, we dissected the neonatal brain at three different levels: acute 1-mm thick brain slice, cultured organotypic 350-µm thick brain slice and dissociated neuronal cultures. The spatio-temporal properties of neural activity were investigated by using a 60-channel Micro-electrode arrays (MEA), and the cell assemblies were studied by using a template-matching method. We find local on-propagating as well as large- scale propagating spontaneous oscillatory activity in acute slices, spontaneous network activity characterized by synchronized burst discharges in organotypic cultured slices, and autonomous bursting behaviour in dissociated neuronal cultures. Furthermore, repetitive spike patterns emerge after one week of dissociated neuronal culture and dramatically increase their numbers as well as their complexity and occurrence in the second week. Our data indicate that neurons can self-organize themselves, assembly to a neural network, present spontaneous oscillations, and emerge spatio-temporal activation patterns. The spontaneous oscillations and repetitive spike patterns may serve fundamental functions for information processing and storage in the brain.

Neural network activity in the neonatal rat barrel cortex in vivo

Coordinated patterns of electrical activity are important for the early development of sensory systems. The spatiotemporal dynamics of these early activity patterns and the role of the peripheral sensory input for their generation are essentially unknown. There are two projects in this thesis. In project1, we performed extracellular multielectrode recordings in the somatosensory cortex of postnatal day 0 to 7 rats in vivo and observed three distinct patterns of synchronized oscillatory activity. (1) Spontaneous and periphery-driven spindle bursts of 1–2 s in duration and ~10 Hz in frequency occurred approximately every 10 s. (2) Spontaneous and sensory-driven gamma oscillations of 150–300 ms duration and 30–40 Hz in frequency occurred every 10–30 s. (3) Long oscillations appeared only every ~20 min and revealed the largest amplitude (250–750 µV) and longest duration (>40 s). These three distinct patterns of early oscillatory activity differently synchronized the neonatal cortical network. Whereas spindle bursts and gamma oscillations did not propagate and synchronized a local neuronal network of 200–400 µm in diameter, long oscillations propagated with 25–30 µm/s and synchronized 600-800 µm large ensembles. All three activity patterns were triggered by sensory activation. Single electrical stimulation of the whisker pad or tactile whisker activation elicited neocortical spindle bursts and gamma activity. Long oscillations could be only evoked by repetitive sensory stimulation. The neonatal oscillatory patterns in vivo depended on NMDAreceptor-mediated synaptic transmission and gap junctional coupling. Whereas spindle bursts and gamma oscillations may represent an early functional columnar-like pattern, long oscillations may serve as a propagating activation signal consolidating these immature neuronal networks. In project2, Using voltage-sensitive dye imaging and simultaneous multi-channel extracellular recordings in the barrel cortex and somatosensory thalamus of newborn rats in vivo, we found that spontaneous and whisker stimulation induced activity patterns were restricted to functional cortical columns already at the day of birth. Spontaneous and stimulus evoked cortical activity consisted of gamma oscillations followed by spindle bursts. Spontaneous events were mainly generated in the thalamus or by spontaneous whisker movements. Our findings indicate that during early developmental stages cortical networks self-organize in ontogenetic columns via spontaneous gamma oscillations triggered by the thalamus or sensory periphery.

Granulite facies metamorphism and partial melting processes in the Ivrea Zone, Northern Italy

The Ivrea Zone in northern Italy has been the focus of numerous petrological, geochemical and structural studies. It is commonly inferred to represent an almost complete section through the mid to lower continental crust, in which metamorphism and partial melting of the abundant metapelites was the result of magmatic underplating by a large volume of mantle-derived magma. This study concerns amphibolite and granulite facies metamorphism in the Ivrea Zone with focus on metapelites and metapsammites/metagreywackes from Val Strona di Omegna and metapelites from Val Sesia and Val Strona di Postua, with the aim to better constrain their metamorphic evolution as well as their pressure and temperature conditions via phase equilibria modelling.rnrnIn Val Strona di Omegna, the metapelites show a structural and mineralogical change from mica-schists with the common assemblage bi-mu-sill-pl-q-ilm ± liq at the lowest grades, through metatexitic migmatites (g-sill-bi-ksp-pl-q-ilm-liq) at intermediate grades, to complex diatexitic migmatites (g-sill-ru-bi-ksp-pl-q-ilm-liq) at the highest grades. Within this section several mappable isograds occur, including the first appearance of K-feldspar in the metapelites, the first appearance of orthopyroxene in the metabasites and the disappearance of prograde biotite from the metapelites. The inferred onset of partial melting in the metapelites occurs around Massiola. The prograde suprasolidus evolution of the metapelites is consistent with melting via the breakdown of first muscovite then biotite. Maximum modelled melt fractions of 30–40 % are predicted at the highest grade. The regional metamorphic field gradient in Val Strona di Omegna is constrained to range from conditions of 3.5–6.5 kbar at T = 650–730 °C to P > 9 kbar at T > 900 °C. The peak P–T estimates, particularly for granulite facies conditions, are significantly higher (around 100 °C) than those of most previous studies. In Val Sesia and Val Strona di Postua to the south the exposure is more restricted. P–T estimates for the metapelites are 750–850 °C and 5–6.5 kbar in Val Sesia and approximately 800–900 °C and 5.5–7 kbar in Val Strona di Postua. These results show similar temperatures but lower pressure than metapelites in Val Strona di Omegna. Metapelites in Val Sesia in contact with the Mafic Complex exhibit a metatexitic structure, while in Val Strona di Postua diatexitic structures occur. Further, metapelites at the contact with the Mafic Complex contain cordierite (± spinel) that overprint the regional metamorphic assemblages and are interpreted to have formed during contact metamorphism related to intrusion of the Mafic Complex. The lower pressures in the high-grade rocks in Val Sesia and Val Strona di Postua are consistent with some decompression from the regional metamorphic peak prior to the intrusion of the Mafic Complex, suggesting the rocks followed a clockwise P–T path. In contrast, the metapelites in Val Strona di Omegna, especially in the granulite facies, do not contain any cordierite or any evidence for a contact metamorphic overprint. The extrapolated granulite facies mineral isograds are cut by the rocks of the Mafic Complex to the south. Therefore, the Mafic Complex cannot have caused the regional metamorphism and it is unlikely that the Mafic Complex occurs in Val Strona di Omegna.