Notas del primer seminario de integrabilidad de la Universitat Politècnica de Catalunya

Estas notas corresponden a las exposiciones presentadas en el \emph{Primer Seminario de Integrabilidad}, dentro de lo que se denomina \emph{Aula de Sistemas Din\'amicos}. Durante este evento se realizaron seis conferencias, todas presentadas por miembros del grupo de Sistemas Din\'amicos de la UPC. El programa desarrollado fue el siguiente:\\\begin{center}AULA DE SISTEMAS DIN\'AMICOS\end{center}\begin{center}\texttt{http://www.ma1.upc.es/recerca/seminaris/aulasd-cat.html}\end{center}\begin{center}SEMINARIO DE INTEGRABILIDAD\end{center}\begin{center}Martes 29 y Mi\'ercoles 30 de marzo de 2005\\Facultad de Matem\'aticas y Estad\'{\i}stica, UPC\\Aula: Seminario 1\end{center}\bigskip\begin{center}PROGRAMA Y RES\'UMENES\end{center}{\bf Martes 29 de marzo}\begin{itemize}\item15:30. Juan J. Morales-Ruiz. \emph{El problema de laintegrabilidad en Sistemas Din\'amicos}\medskip {\bf Resumen.} En esta presentaci\'on se pretende dar unaidea de conjunto, pero sin entrar en detalles, sobre las diversasnociones de integrabilidad, asociadas a nombres de matem\'aticostan ilustres como Liouville, Galois-Picard-Vessiot, Lie, Darboux,Kowalevskaya, Painlev\'e, Poincar\'e, Kolchin, Lax, etc. Adem\'astambi\'en mencionaremos la revoluci\'on que supuso en los a\~nossesenta del siglo pasado el descubrimiento de Gardner, Green,Kruskal y Miura sobre un nuevo m\'etodo para resolver en algunoscasos determinadas ecuaciones en derivadas parciales. \medskip\item16:00. David G\'omez-Ullate. \emph{Superintegrabilidad, pares deLax y modelos de $N-$cuerpos en el plano}\medskip{\bf Resumen.} Introduciremos algunas t\'ecnicas cl\'asicas paraconstruir modelos de N-cuerpos integrables, como los pares de Laxo la din\'amica de los ceros de un polinomio. Revisaremos lanoci\'on de integrabilidad Liouville y superintegrabilidad, ydiscutiremos un nuevo m\'etodo debido a F. Calogero para contruirmodelos de N-cuerpos en el plano con muchas \'orbitasperi\'odicas. La exposici\'on se acompa\~nar\'a de animaciones delmovimiento de los cuerpos, y se plantear\'an algunos problemasabiertos.\medskip\item17:00. Pausa\medskip\item17:30. Yuri Fedorov. \emph{An\'alisis de Kovalevskaya--Painlev\'ey Sistemas Algebraicamente Integrables}\medskip{\bf Resumen.} Muchos sistemas integrables poseen una propiedadremarcable: todas sus soluciones son funciones meromorfas deltiempo como una variable compleja. Tal comportamiento, que serefiere como propiedad de Kovalevskaya-Painleve (KP) y que se usafrecuentemente como una ensayo de integrabilidad, no es accidentaly tiene unas ra\'{\i}ces geom\'etricas profundas. En esta charladescribiremos una clase de tales sistemas (conocidos como lossistemas algebraicamente integrables) y subrayaremos suspropiedades geom\'etricas principales que permiten predecir laestructura de las soluciones complejas y adem\'as encontrarlasexpl\'{\i}citamente. Eso lo ilustraremos con algunos sistemas dela mec\'anica cl\'asica. Tambi\'en mencionaremos unasgeneralizaciones \'utiles de la noci\'on de integrabilidadalgebraica y de la propiedad KP.\end{itemize}\medskip{\bf Mi\'ercoles 30 de marzo}\begin{itemize}\item 15:30. Rafael Ram\'{\i}rez-Ros. \emph{El m\'etodo de Poincar\'e}\medskip{\bf Resumen.} Dado un sistema Hamiltoniano aut\'onomo cercano acompletamente integrable Poincar\'e prob\'o que, en general, noexiste ninguna integral primera adicional uniforme en elpar\'ametro de perturbaci\'on salvo el propio Hamiltoniano.Esbozaremos las ideas principales del m\'etodo de prueba ycomentaremos algunas extensiones y generalizaciones.\newpage\item16:30. Chara Pantazi. \emph{El M\'etodo de Darboux}\medskip{\bf Resumen.} Darboux, en 1878, present\'o su m\'etodo paraconstruir integrales primeras de campos vectoriales polinomialesutilizando sus curvas invariantes algebraicas. En estaexposici\'on presentaremos algunas extensiones del m\'etodocl\'asico de Darboux y tambi\'en algunas aplicaciones.\medskip\item17:30. Pausa\medskip\item18:00. Juan J. Morales-Ruiz. \emph{M\'etodos recientes paradetectar la no integrabilidad}\medskip{\bf Resumen.} En 1982 Ziglin utiliza la estructura de laecuaci\'on en variaciones de Poincar\'e (sobre una curva integralparticular) como una herramienta fundamental para detectar la nointegrabilidad de un sistema Hamiltoniano. En esta charla sepretende dar una idea de esta aproximaci\'on a la nointegrabilidad, junto con t\'ecnicas m\'as recientes queinvolucran la teor\'{\i}a de Galois de ecuaciones diferencialeslineales, haciendo \'enfasis en los ejemplos m\'as que en lateor\'{\i}a general. Ilustraremos estos m\'etodos con resultadossobre la no integrabilidad de algunos problemas de $N$ cuerpos enMec\'anica Celeste.\end{itemize}

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

Alteration of neuromuscular function in squash

The alteration in neuromuscular function of knee extensor muscles was characterised after a squash match in 10 trained players. Maximal voluntary contraction (MVC) and surface EMG activity of vastus lateralis (VL) and vastus medialis (VM) muscles were measured before and immediately after a 1-h squash match. M-wave and twitch contractile properties were analysed following single stimuli. MVC declined (280.5+/-46.8 vs. 233.6+/-35.4 Nm, -16%; P<0.001) after the exercise and this was accompanied by an impairment of central activation, as attested by decline in voluntary activation (76.7+/-10.4 vs. 71.3+/-9.6%, -7%; P<0.05) and raw EMG activity of the two vastii (-17%; P<0.05), whereas RMS/M decrease was lesser (VL: -5%; NS and VM: -12%; P=0.10). In the fatigued state, no significant changes in M-wave amplitude (VL: -9%; VM: -5%) or duration were observed. Following exercise, the single twitch was characterised by lower peak torque (-20%; P<0.001) as well as shorter half-relaxation time (-13%; P<0.001) and reduced maximal rate of twitch tension development (-23%; P<0.001) and relaxation (-17%; P<0.05). A 1-h squash match play caused peripheral fatigue by impairing excitation-contraction coupling, whereas sarcolemmal excitability seems well preserved. Our results also emphasise the role of central activation failure as a possible mechanism contributing to the torque loss observed in knee extensors. Physical conditioners should consider these effects when defining their training programs for squash players.

Free-breathing T2 mapping at 3T for the monitoring of cardiac allograft rejection: initial results

rejection can lead to loss of function. Histological reading of endomyocardial biopsy remains the "gold standard" for guiding immunosuppression, despite its methodological limitations (sampling error and interobserver variability). The measurement of the T2 relaxation time has been suggested for detection of allograft rejection, on the pathophysiological basis that the T2 relaxation time prolongs with local edema resulting from acute allograft rejection. Using breath-held cardiac magnetic resonance T2 mapping at 1.5 T, Usman et al. (CircCardiovascImaging2012) detected moderate allograft rejection (grade 2R, ISHLT 2004). With modern immunosuppression grade 2R rejection has become a rare event, but the need remains for a technique that permits the discrimination of absent (grade 0R) and mild rejection (grade 1R). We therefore investigated whether an increase of magnetic field strength to 3T and the use of real-time navigator-gated respiration compensation allow for an increase in the sensitivity of T2 relaxation time detection that is necessary to achieve this discrimination. Methods: Eighteen patients received EMB (Tan et al., ArchPatholLabMed2007) and cardiac T2 mapping on the same day. Reading of T2 maps was blinded to the histological results. For final analysis, 3 cases with known 2R rejection at the time of T2 mapping were added, yielding 21 T2 mapping sessions. A respiration-navigator-gated radial gradient-recalled-echo pulse sequence (resolution 1.17 mm2, matrix 2562, trigger time 3 heartbeats, T2 preparation duration TET2 Prep = 60/30/0 ms) was applied to obtain 3 short-axis T2 maps (van Heeswijk et al., JACCCardiovascImaging2012), which were segmented according to AHA guidelines (Cerqueira et al, Circulation2001). The highest segmental T2 values were grouped according to histological rejection grade and differences were analyzed by Student's t-test, except for the non-blinded cases with 2R rejection. The degree of discrimination was determined using the Spearman's ranked correlation test. Results: The high-quality T2 maps allowed for visual differentiation of the rejection degrees (Figure 1), and the correlation of T2 mapping with the histological grade of acute cellular rejection was significant (Spearman's r = 0.56, p = 0.007). The 0R (n = 15) and 1R (n = 3) degrees demonstrated significantly different T2 values (46.9 ± 5.0 and 54.3 ± 3.0 ms, p = 0.02, Figure 2). Cases with 2R rejection showed clear T2 elevation (T2 = 60.3 ± 16.2 ms). Conclusions: This pilot study demonstrates that non-invasive free-breathing cardiac T2 mapping at 3T discriminates between no and mild cardiac allograft rejection. Confirmation of these encouraging results in a larger cohort should consider a study able to show equivalency or superiority of T2 mapping.

Colloid cysts of the third ventricle: correlation of MR and CT findings with histology and chemical analysis.

Eight patients with colloid cysts of the third ventricle were examined with CT and MR. In six, surgical resection was performed and the material was subjected to histologic evaluation; the concentrations of trace elements were determined by particle-induced X-ray emission. Stereotaxic aspiration was performed in two. The investigation showed that colloid cysts are often iso- or hypodense relative to brain on CT (5/8), but sometimes have a center of increased density. Increased density did not correlate with increased concentration of calcium or other metals but did not correlate with high cholesterol content. Colloid cysts appear more heterogeneous on MR (6/8) than on CT (3/8), despite a homogeneous appearance at histology. High signal on short TR/TE sequences is correlated with a high cholesterol content. A marked shortening of the T2 relaxation time is often noticed in the central part of the cyst. Analysis of trace elements showed that this phenomenon is not related to the presence of metals with paramagnetic effects. Our analysis of the contents of colloid cysts does not support the theory that differing metallic concentrations are responsible for differences in MR signal intensity or CT density. We did find that increased CT density and high MR signal correlated with high cholesterol content.

An active strain electromechanical model for cardiac tissue

We propose a finite element approximation of a system of partial differential equations describing the coupling between the propagation of electrical potential and large deformations of the cardiac tissue. The underlying mathematical model is based on the active strain assumption, in which it is assumed that a multiplicative decomposition of the deformation tensor into a passive and active part holds, the latter carrying the information of the electrical potential propagation and anisotropy of the cardiac tissue into the equations of either incompressible or compressible nonlinear elasticity, governing the mechanical response of the biological material. In addition, by changing from an Eulerian to a Lagrangian configuration, the bidomain or monodomain equations modeling the evolution of the electrical propagation exhibit a nonlinear diffusion term. Piecewise quadratic finite elements are employed to approximate the displacements field, whereas for pressure, electrical potentials and ionic variables are approximated by piecewise linear elements. Various numerical tests performed with a parallel finite element code illustrate that the proposed model can capture some important features of the electromechanical coupling, and show that our numerical scheme is efficient and accurate.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Investigation of interaction between native and impurity defects in ZnSe

Zinc selenide is a prospective material for optoelectronics. The fabrication of ZnSe­based light-emitting diodes is hindered by complexity of p-type doping of the component materials. The interaction between native and impurity defects, the tendency of doping impurity to form associative centres with native defects and the tendency to self-compensation are the main factors impeding effective control of the value and type of conductivity. The thesis is devoted to the study of the processes of interaction between native and impurity defects in zinc selenide. It is established that the Au impurity has the most prominent amphoteric properties in ZnSe among Cu, Ag and Au impurities, as it forms a great number of both Au; donors and Auz„ acceptors. Electrical measurements show that Ag and Au ions introduced into vacant sites of the Zn sublattice form simple single-charged Agz„+ and Auzn+ states with d1° electron configuration, while Cu ions can form both single-charged Cuz„ (d1) and double-charged Cuzr`+ (d`o) centres. Amphoteric properties of Ag and Au transition metals stimulated by time are found for the first time from both electrical and luminescent measurements. A model that explains the changes in electrical and luminescent parameters by displacement of Ag ions into interstitial sites due to lattice deformation forces is proposed. Formation of an Ag;-donor impurity band in ZnSe samples doped with Ag and stored at room temperature is also studied. Thus, the properties of the doped samples are modified due to large lattice relaxation during aging. This fact should be taken into account in optoelectronic applications of doped ZnSe and related compounds.

Interplay between connexin40 and nitric oxide signaling during hypertension.

Connexins (Cxs) and endothelial nitric oxide synthase (eNOS) contribute to the adaptation of endothelial and smooth muscle cells to hemodynamic changes. To decipher the in vivo interplay between these proteins, we studied Cx40-null mice, a model of renin-dependent hypertension which displays an altered endothelium-dependent relaxation of the aorta because of reduced eNOS levels. These mice, which were either untreated or subjected to the 1-kidney, 1-clip (1K1C) procedure, a model of volume-dependent hypertension, were compared with control mice submitted to either the 1K1C or the 2-kidney, 1-clip (2K1C) procedure, a model of renin-dependent hypertension. All operated mice became hypertensive and featured hypertrophy and altered Cx expression of the aorta. The combination of volume- and renin-dependent hypertension in Cx40-/- 1K1C mice raised blood pressure and cardiac weight index. Under these conditions, all aortas showed increased levels of Cx40 in endothelial cells and of both Cx37 and Cx45 in smooth muscle cells. In the wild-type 1K1C mice, the interactions between Cx40 and Cx37 with eNOS were enhanced, resulting in increased NO release. The Cx40-eNOS interaction could not be observed in mice lacking Cx40, which also featured decreased levels of eNOS. In these animals, the volume overload caused by the 1K1C procedure resulted in increased phosphorylation of eNOS and in a higher NO release. The findings provide evidence that Cx40 and Cx37 play an in vivo role in the regulation of eNOS.

Hienopaperimassojen lujuusominaisuudet nopeassa vetokuormituksessa

Työssä verrattiin koivu-, akaasia- ja eucalyptussellujen käyttökelpoisuutta hienopapereiden kuituraaka-aineena. Kirjallisuusosassa todettiin radan hallinnan paperikoneiden avoimissa vienneissä riippuvan käytetyn geometrian lisäksi lujuus- ja relaksaatio-ominaisuuksista. Kuituverkoston käyttäytymiseen jännityksen alaisena vaikuttavat kuitudimensiot, kosteus ja lämpötila sekä kuituverkostoon kohdistetun jännityksen nopeus ja määrä. Relaksaationopeus ja vetolujuus kasvavat kuivilla papereilla vetonopeuden lisääntyessä. Kosteuspitoisuuden kasvattaminen alentaa puukuiduissa olevien polymeerien lasisiirtymälämpötilaa, jonka seurauksena vetolujuus ja relaksaatiokireys laskevat voimakkaasti. Kosteuspitoisuuden kasvaessa murtovenymä kasvaa lähes lineaarisesti ja repäisylujuus sekä murtotyö saavuttavat maksiminsa tietyssä kosteuspitoisuudessa. Kokeellisessa osassa keskityttiin hienopaperimassojen lujuus- ja relaksaatiokäyttäytymisen selvittämiseen nopeassa vetokuormituksessa. Lisäksi määritettiin laatu-, massa- ja rakenneominaisuuksia valituille koepisteille. Muuttujina kokeissa olivat massojen kuiva-ainepitoisuudet ja jauhatusolosuhteet sekä havusellun osuus hienopaperimassoissa.CSF-tasoon 350 ml jauhetuista näytteistä parhaat lujuus ja relaksaatio-ominaisuudet olivat koivulla ja heikoimmat akaasialla. Erot koepisteiden välillä korostuivat pienellä havusellun määrällä, mutta kaventuivat huomattavasti havusellun määrää lisättäessä. Samaan vetolujuuteen jauhettaessa massojen erot poistuvat kokonaan. Kuivilla näytteillä löydettiin erinomainen korrelaatio myötölujuuden ja relaksaatiokireyden välille. Puristinkuivien näytteiden relaksaatiokireyksiä voidaan kokeiden valossa ennustaa parhaiten kuivien näytteiden vetolujuuksista. Myös elastisten venymien osuuksille ja kuituseinämien paksuuksille löydettiin selvä yhteys.Eucalyptus- ja akaasiamassojen erinomaisuus hienopapereiden raaka-aineena korostui niiden optisissa ominaisuuksissa, erityisesti korkeana valonsirontana. Verrattaessa samassa vetolujuudessa ja relaksaatiokireydessä valonsirontakertoimien arvoja havaittiin akaasian olevan paras koivun jäädessä heikoimmaksi. Lisäksi akaasian ja eucalyptuksen kapeat kuitujakaumat ovat edullisia painokoneessa värin tasaisen imeytymisen kannalta. Akaasian pienet ja taipuisat kuidut antavat paperille tasaisen pinnan ja siten painatuksessa tasaisen painoalustan.

Sherali-Adams Relaxations and Indistinguishability in Counting Logics

Two graphs with adjacency matrices $\mathbf{A}$ and $\mathbf{B}$ are isomorphic if there exists a permutation matrix $\mathbf{P}$ for which the identity $\mathbf{P}^{\mathrm{T}} \mathbf{A} \mathbf{P} = \mathbf{B}$ holds. Multiplying through by $\mathbf{P}$ and relaxing the permutation matrix to a doubly stochastic matrix leads to the linear programming relaxation known as fractional isomorphism. We show that the levels of the Sherali--Adams (SA) hierarchy of linear programming relaxations applied to fractional isomorphism interleave in power with the levels of a well-known color-refinement heuristic for graph isomorphism called the Weisfeiler--Lehman algorithm, or, equivalently, with the levels of indistinguishability in a logic with counting quantifiers and a bounded number of variables. This tight connection has quite striking consequences. For example, it follows immediately from a deep result of Grohe in the context of logics with counting quantifiers that a fixed number of levels of SA suffice to determine isomorphism of planar and minor-free graphs. We also offer applications in both finite model theory and polyhedral combinatorics. First, we show that certain properties of graphs, such as that of having a flow circulation of a prescribed value, are definable in the infinitary logic with counting with a bounded number of variables. Second, we exploit a lower bound construction due to Cai, Fürer, and Immerman in the context of counting logics to give simple explicit instances that show that the SA relaxations of the vertex-cover and cut polytopes do not reach their integer hulls for up to $\Omega(n)$ levels, where $n$ is the number of vertices in the graph.

Spatiotemporal variability in radial growth of trees in the Central Pyrenees: climatic and site influences.

To understand how tree growth has responded to recent climate warming, an understanding of the tree-climate-site complex is necessary. To achieve this, radial growth variability among 204 trees established before 1850 was studied in relation to both climatic and site factors. Seventeen forest stands were sampled in the Spanish Central Pyrenees. Three species were studied: Pinus uncinata, Abies alba, and Pinus sylvestris. For each tree, a ring-width residual chronology was built. All trees cross-dated well, indicating a common influence of the regional climate. For the 1952-1993 period, the radial growth of all species, especially P. uncinata, was positively correlated with warm Novembers during the year before ring formation and warm Mays of the year the annual ring formed. Differences in species-stand elevation modulated the growth-climate associations. Radial growth in P. uncinata at high elevation sites was reduced when May temperatures were colder and May precipitation more abundant. In the 20th century, two contrasting periods in radial growth were observed: one (1900-1949) with low frequency of narrow and wide rings, low mean annual sensitivity, and low common growth variation; and another (1950-1994) with the reverse characteristics. The increased variability in radial growth since the 1950s was observed for all species and sites, which suggests a climatic cause. The low shared variance among tree chronologies during the first half of the 20th century may result from a"relaxation" of the elevation gradient, allowing local site conditions to dominate macroclimatic influence. These temporal trends may be related to the recently reported increase of climatic variability and warmer conditions. This study emphasizes the need to carefully assess the relationships between radial growth and site conditions along ecological gradients to improve dendroclimatic reconstructions.

From caging to rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory

By means of computer simulations and solution of the equations of the mode coupling theory (MCT),we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question,provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This dis-agreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.

1C.06: Ambulatory pulse pressure is negatively associated with excretions of urinary caffeine and its metabolites

OBJECTIVE: Systolic blood pressure (BP) has been associated with urinary caffeine and its metabolites such as paraxanthine and theophylline. Caffeine and caffeine metabolites could influence arterial pulse pressure (PP) via sympathomimetic effects, smooth muscle relaxation, and phosphodiesterase inhibition. The purpose of this analysis was to explore the association of ambulatory PP with urinary caffeine and its related metabolites in a large population-based sample. DESIGN AND METHOD: Families were randomly selected from the general population of three Swiss cities (2009-2013). Ambulatory BP monitoring was conducted using validated Diasys Integra devices. PP was defined as the difference between the systolic and diastolic ambulatory BP. Urinary caffeine, paraxanthine, theophylline, and theobromine excretions were measured in 24 h urine using ultra-high performance liquid chromatography tandem mass spectrometry. Urinary excretions were log-transformed to satisfy regression assumptions. We used linear mixed models to explore the associations of urinary caffeine and caffeine metabolite excretions with 24-hour, day- and night-time PP while adjusting for major confounders. RESULTS: The 836 participants (48.9% men) included in this analysis had mean (±SD) age of 47.8 (±17.5), and mean 24-hour systolic and diastolic BP of 120.1 mmHg (±13.9) and 78.0 (±8.6). Except theobromine, log transformed urinary caffeine and caffeine metabolite excretions were associated negatively with 24-hour, daytime and night-time ambulatory PP. 24-hour, daytime, and night-time ambulatory PP decreased by -0.804 mmHg (SE, 0.209), -0.749 (0.215), and -0.968 (0.243) (all P values <0.005), for each doubling excretion of caffeine. Strong negative associations with night-time ambulatory PP were observed for paraxanthine and theophylline.(Figure is included in full-text article.) CONCLUSIONS: : The negative associations of PP with caffeine, paraxanthine, and theophylline excretions suggest that caffeine and its metabolites do lower BP, possibly by modifying arterial stiffness.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.