Análisis de la actividad neuronal y cerebral mediante técnicas de redes complejas

El cerebro humano es probablemente uno de los sistemas más complejos a los que nos enfrentamos en la actualidad, si bien es también uno de los más fascinantes. Sin embargo, la compresión de cómo el cerebro organiza su actividad para llevar a cabo tareas complejas es un problema plagado de restos y obstáculos. En sus inicios la neuroimagen y la electrofisiología tenían como objetivo la identificación de regiones asociadas a activaciones relacionadas con tareas especificas, o con patrones locales que variaban en el tiempo dada cierta actividad. Sin embargo, actualmente existe un consenso acerca de que la actividad cerebral tiene un carácter temporal multiescala y espacialmente extendido, lo que lleva a considerar el cerebro como una gran red de áreas cerebrales coordinadas, cuyas conexiones funcionales son continuamente creadas y destruidas. Hasta hace poco, el énfasis de los estudios de la actividad cerebral funcional se han centrado en la identidad de los nodos particulares que forman estas redes, y en la caracterización de métricas de conectividad entre ellos: la hipótesis subyacente es que cada nodo, que es una representación mas bien aproximada de una región cerebral dada, ofrece a una única contribución al total de la red. Por tanto, la neuroimagen funcional integra los dos ingredientes básicos de la neuropsicología: la localización de la función cognitiva en módulos cerebrales especializados y el rol de las fibras de conexión en la integración de dichos módulos. Sin embargo, recientemente, la estructura y la función cerebral han empezado a ser investigadas mediante la Ciencia de la Redes, una interpretación mecánico-estadística de una antigua rama de las matemáticas: La teoría de grafos. La Ciencia de las Redes permite dotar a las redes funcionales de una gran cantidad de propiedades cuantitativas (robustez, centralidad, eficiencia, ...), y así enriquecer el conjunto de elementos que describen objetivamente la estructura y la función cerebral a disposición de los neurocientíficos. La conexión entre la Ciencia de las Redes y la Neurociencia ha aportado nuevos puntos de vista en la comprensión de la intrincada anatomía del cerebro, y de cómo las patrones de actividad cerebral se pueden sincronizar para generar las denominadas redes funcionales cerebrales, el principal objeto de estudio de esta Tesis Doctoral. Dentro de este contexto, la complejidad emerge como el puente entre las propiedades topológicas y dinámicas de los sistemas biológicos y, específicamente, en la relación entre la organización y la dinámica de las redes funcionales cerebrales. Esta Tesis Doctoral es, en términos generales, un estudio de cómo la actividad cerebral puede ser entendida como el resultado de una red de un sistema dinámico íntimamente relacionado con los procesos que ocurren en el cerebro. Con este fin, he realizado cinco estudios que tienen en cuenta ambos aspectos de dichas redes funcionales: el topológico y el dinámico. De esta manera, la Tesis está dividida en tres grandes partes: Introducción, Resultados y Discusión. En la primera parte, que comprende los Capítulos 1, 2 y 3, se hace un resumen de los conceptos más importantes de la Ciencia de las Redes relacionados al análisis de imágenes cerebrales. Concretamente, el Capitulo 1 está dedicado a introducir al lector en el mundo de la complejidad, en especial, a la complejidad topológica y dinámica de sistemas acoplados en red. El Capítulo 2 tiene como objetivo desarrollar los fundamentos biológicos, estructurales y funcionales del cerebro, cuando éste es interpretado como una red compleja. En el Capítulo 3, se resumen los objetivos esenciales y tareas que serán desarrolladas a lo largo de la segunda parte de la Tesis. La segunda parte es el núcleo de la Tesis, ya que contiene los resultados obtenidos a lo largo de los últimos cuatro años. Esta parte está dividida en cinco Capítulos, que contienen una versión detallada de las publicaciones llevadas a cabo durante esta Tesis. El Capítulo 4 está relacionado con la topología de las redes funcionales y, específicamente, con la detección y cuantificación de los nodos mas importantes: aquellos denominados “hubs” de la red. En el Capítulo 5 se muestra como las redes funcionales cerebrales pueden ser vistas no como una única red, sino más bien como una red-de-redes donde sus componentes tienen que coexistir en una situación de balance funcional. De esta forma, se investiga cómo los hemisferios cerebrales compiten para adquirir centralidad en la red-de-redes, y cómo esta interacción se mantiene (o no) cuando se introducen fallos deliberadamente en la red funcional. El Capítulo 6 va un paso mas allá al considerar las redes funcionales como sistemas vivos. En este Capítulo se muestra cómo al analizar la evolución de la topología de las redes, en vez de tratarlas como si estas fueran un sistema estático, podemos caracterizar mejor su estructura. Este hecho es especialmente relevante cuando se quiere tratar de encontrar diferencias entre grupos que desempeñan una tarea de memoria, en la que las redes funcionales tienen fuertes fluctuaciones. En el Capítulo 7 defino cómo crear redes parenclíticas a partir de bases de datos de actividad cerebral. Este nuevo tipo de redes, recientemente introducido para estudiar las anormalidades entre grupos de control y grupos anómalos, no ha sido implementado nunca en datos cerebrales y, en este Capítulo explico cómo hacerlo cuando se quiere evaluar la consistencia de la dinámica cerebral. Para concluir esta parte de la Tesis, el Capítulo 8 se centra en la relación entre las propiedades topológicas de los nodos dentro de una red y sus características dinámicas. Como mostraré más adelante, existe una relación entre ellas que revela que la posición de un nodo dentro una red está íntimamente correlacionada con sus propiedades dinámicas. Finalmente, la última parte de esta Tesis Doctoral está compuesta únicamente por el Capítulo 9, el cual contiene las conclusiones y perspectivas futuras que pueden surgir de los trabajos expuestos. En vista de todo lo anterior, espero que esta Tesis aporte una perspectiva complementaria sobre uno de los más extraordinarios sistemas complejos frente a los que nos encontramos: El cerebro humano. ABSTRACT The human brain is probably one of the most complex systems we are facing, thus being a timely and fascinating object of study. Characterizing how the brain organizes its activity to carry out complex tasks is highly non-trivial. While early neuroimaging and electrophysiological studies typically aimed at identifying patches of task-specific activations or local time-varying patterns of activity, there has now been consensus that task-related brain activity has a temporally multiscale, spatially extended character, as networks of coordinated brain areas are continuously formed and destroyed. Up until recently, though, the emphasis of functional brain activity studies has been on the identity of the particular nodes forming these networks, and on the characterization of connectivity metrics between them, the underlying covert hypothesis being that each node, constituting a coarse-grained representation of a given brain region, provides a unique contribution to the whole. Thus, functional neuroimaging initially integrated the two basic ingredients of early neuropsychology: localization of cognitive function into specialized brain modules and the role of connection fibres in the integration of various modules. Lately, brain structure and function have started being investigated using Network Science, a statistical mechanics understanding of an old branch of pure mathematics: graph theory. Network Science allows endowing networks with a great number of quantitative properties, thus vastly enriching the set of objective descriptors of brain structure and function at neuroscientists’ disposal. The link between Network Science and Neuroscience has shed light about how the entangled anatomy of the brain is, and how cortical activations may synchronize to generate the so-called functional brain networks, the principal object under study along this PhD Thesis. Within this context, complexity appears to be the bridge between the topological and dynamical properties of biological systems and, more specifically, the interplay between the organization and dynamics of functional brain networks. This PhD Thesis is, in general terms, a study of how cortical activations can be understood as the output of a network of dynamical systems that are intimately related with the processes occurring in the brain. In order to do that, I performed five studies that encompass both the topological and the dynamical aspects of such functional brain networks. In this way, the Thesis is divided into three major parts: Introduction, Results and Discussion. In the first part, comprising Chapters 1, 2 and 3, I make an overview of the main concepts of Network Science related to the analysis of brain imaging. More specifically, Chapter 1 is devoted to introducing the reader to the world of complexity, specially to the topological and dynamical complexity of networked systems. Chapter 2 aims to develop the biological, topological and functional fundamentals of the brain when it is seen as a complex network. Next, Chapter 3 summarizes the main objectives and tasks that will be developed along the forthcoming Chapters. The second part of the Thesis is, in turn, its core, since it contains the results obtained along these last four years. This part is divided into five Chapters, containing a detailed version of the publications carried out during the Thesis. Chapter 4 is related to the topology of functional networks and, more specifically, to the detection and quantification of the leading nodes of the network: the hubs. In Chapter 5 I will show that functional brain networks can be viewed not as a single network, but as a network-of-networks, where its components have to co-exist in a trade-off situation. In this way, I investigate how the brain hemispheres compete for acquiring the centrality of the network-of-networks and how this interplay is maintained (or not) when failures are introduced in the functional network. Chapter 6 goes one step beyond by considering functional networks as living systems. In this Chapter I show how analyzing the evolution of the network topology instead of treating it as a static system allows to better characterize functional networks. This fact is especially relevant when trying to find differences between groups performing certain memory tasks, where functional networks have strong fluctuations. In Chapter 7 I define how to create parenclitic networks from brain imaging datasets. This new kind of networks, recently introduced to study abnormalities between control and anomalous groups, have not been implemented with brain datasets and I explain in this Chapter how to do it when evaluating the consistency of brain dynamics. To conclude with this part of the Thesis, Chapter 8 is devoted to the interplay between the topological properties of the nodes within a network and their dynamical features. As I will show, there is an interplay between them which reveals that the position of a node in a network is intimately related with its dynamical properties. Finally, the last part of this PhD Thesis is composed only by Chapter 9, which contains the conclusions and future perspectives that may arise from the exposed results. In view of all, I hope that reading this Thesis will give a complementary perspective of one of the most extraordinary complex systems: The human brain.

El enfoque integrado en los programas de regeneración urbana de barrios desfavorecidos: una visión relacional

Actualmente se reconoce de forma generalizada que la problemática de los barrios más desfavorecidos presenta un trasfondo estructural de magnitud multidimensional que implica a múltiples actores, que requiere respuestas de carácter multiescalar, multidimensional y multiagente y un alto nivel de cooperación institucional y de coordinación entre la diversidad de actores, que pone en cuestión los modelos de gobernanza y gestión urbana tradicionales. Frente a los retos e implicaciones que plantea la noción de espacio como construcción social y los problemas de coordinación detectados que limitan el alcance de la regeneración urbana integrada, en esta investigación se pregunta por los principales elementos que pueden facilitar la dinamización y coordinación de los procesos y se plantea que es preciso aprovechar el propio proceso participativo para reforzar el papel de estos elementos y generar mejores condiciones para la colaboración en un proceso de aprendizaje práctico de formas de análisis e intervención integrada y gestión relacional, que supere los niveles de información y consulta de la práctica participativa habitual. Se centra por tanto la atención en los factores que facilitan la dinamización y coordinación de los procesos (elementos dinamizadores y articuladores), en el nivel de participación social y cooperación institucional en los distintos momentos de la planificación de la estrategia entendida como un proceso continuo de reflexión-acción-reflexión y en aquellas condiciones del modelo de participación que pueden propiciar relaciones estratégicas y mejorar la capacidad de interacción y acción del sistema de actores y relaciones en el proceso de planificación y gestión, que lleva a cuestionar en qué medida los soportes (foros de participación y nodos y ejes de actividad y relación) y el tipo de técnicas (técnicas de dinamización y técnicas de identificación y activación de oportunidades) de los modelos de participación de experiencias concretas pueden fomentar formas de interacción que mejoren la comunicación y la movilización de recursos. Para analizar esta situación en la práctica, se desarrolla un marco conceptual, metodológico y operativo a partir de la contribución de conceptos, teorías y modelos basados en una visión relacional y multidimensional de las dinámicas de desarrollo local integrado. Esta metodología se valida en cuatro casos de referencia a partir de la información de documentos oficiales y de las visiones de la diversidad de actores implicados recogidas a través de una primera consulta a expertos (13) y entrevistas semiestructuradas a representantes de la diversidad de actores implicados (75), que se ha procesado por medio de tres herramientas conceptuales de análisis complementarias desarrolladas a lo largo de la investigación apoyadas en la perspectiva relacional de la reinterpretación actual del espacio y las escalas: el ‘Sistema de acción local’, el ‘Diagrama de momentos’ del proceso de planificación y gestión, y la ‘Matriz de caracterización’ del modelo de participación. La evaluación y la comparativa de los casos según el método de análisis evidencian la importancia del papel ejercido por los elementos articuladores (actores y organismos de enlace, espacios de encuentro y acción, y los significados y visiones compartidas) y dinamizadores (el capital social y las capacidades institucionales propiciadas en redes de colaboración) a la hora de facilitar la colaboración entre las distintas instituciones responsables, reforzar las relaciones entre actores, grupos y entidades del barrio, y mejorar los problemas de convivencia vecinal, y en definitiva, de mejorar las posibilidades de los diferentes actores urbanos de implicarse de manera activa. Además, las fórmulas alternativas que se han puesto en marcha en los distintos casos frente a las dificultades que presenta la administración a la hora de trabajar por proyectos y el resto de actores de colaborar de forma constructiva, han supuesto un importante proceso de aprendizaje y práctica en formas de participación y colaboración. ABSTRACT It is widely accepted that the current situation of multi-deprived neighborhoods has a structural and a multidimensional basis involving multiple actors that calls for multiscale, multidimensional and multiagent solutions and requires a high level of institutional cooperation and a great deal of coordination between different actors. This calls into question traditional governance and management models. In order to cope with the challenges and implications of the idea of the social construction of space and these coordination requirements that compromise the scope of integrated urban regeneration strategies, this research focuses on the main elements and factors that facilitate the mobilization and coordination of the regeneration processes, and recognizes the importance of participation processes facilitating better social, communication and collaboration skills in a practical learning process of integrated analysis and urban regeneration interventions and relational governance mechanisms going beyond mere information and consultation levels. Thus, the study focuses on the principal factors which facilitate the mobilization and coordination of integrated area regeneration, the level of social participation and institutional cooperation at different planning moments in a continuous process of reflection-action-reflection, and the conditions under which the participation process may promote strategic relationships and improve the performance of the system of actors and relationships in the planning and management process. This questions to what extent the supports (participation arenas and activity nodes and axis) and the type of techniques (participation encouragement and opportunities identification and activation techniques) considered in participation models had led to a level of interaction that has improved communication and resource mobilization capabilities. In order to analyze this situation in practice, a conceptual, methodological and operational framework is developed, which considers the contributions of the review of concepts, theories and models based on a relational and multidimensional view of local development dynamics. This methodology is validated in four representative case studies in view of official documents and the contributions of diverse stakeholders gathered through a first consultation to experts (13) and semi-structured interviews to representatives of the diversity of actors involved (75). This data has been processed by three complementary analytical tools developed during the research, based on the relational perspective of the current reinterpretation of space and scales: the 'Local Action System', the 'Diagram of planning moments', and the 'Characterization matrix of the participation model'. The comparative evaluation of the four case studies carried out using the analytical method developed during the project, shows the importance of the availability and the role of linkage elements (actors and organizations, collective spaces and shared visions) and drivers (social capital and institutional capacities generated through collaborative networks) in facilitating the collaboration among all responsible institutions, in promoting and strengthening relations between actors, groups and neighborhood organizations, and in improving the conviviality and social cohesion conditions, and ultimately in fostering the possibilities of diversity of urban actors to engage and participate actively in the transformation process, increasing accordingly the options for change. In addition, despite all the difficulties, these practices have implemented alternative formulas addressing the limitations of the administration in working by projects and not by competences and of the stakeholders in collaborating in a constructive way. The development and implementation of these participatory formulas has been an important learning process.

Coupled Nonlinear Ginzburg-Landau and Mechanics Model for Martensitic Transformations in Polycrystals

Las transformaciones martensíticas (MT) se definen como un cambio en la estructura del cristal para formar una fase coherente o estructuras de dominio multivariante, a partir de la fase inicial con la misma composición, debido a pequeños intercambios o movimientos atómicos cooperativos. En el siglo pasado se han descubierto MT en diferentes materiales partiendo desde los aceros hasta las aleaciones con memoria de forma, materiales cerámicos y materiales inteligentes. Todos muestran propiedades destacables como alta resistencia mecánica, memoria de forma, efectos de superelasticidad o funcionalidades ferroicas como la piezoelectricidad, electro y magneto-estricción etc. Varios modelos/teorías se han desarrollado en sinergia con el desarrollo de la física del estado sólido para entender por qué las MT generan microstructuras muy variadas y ricas que muestran propiedades muy interesantes. Entre las teorías mejor aceptadas se encuentra la Teoría Fenomenológica de la Cristalografía Martensítica (PTMC, por sus siglas en inglés) que predice el plano de hábito y las relaciones de orientación entre la austenita y la martensita. La reinterpretación de la teoría PTMC en un entorno de mecánica del continuo (CM-PTMC) explica la formación de los dominios de estructuras multivariantes, mientras que la teoría de Landau con dinámica de inercia desentraña los mecanismos físicos de los precursores y otros comportamientos dinámicos. La dinámica de red cristalina desvela la reducción de la dureza acústica de las ondas de tensión de red que da lugar a transformaciones débiles de primer orden en el desplazamiento. A pesar de las diferencias entre las teorías estáticas y dinámicas dado su origen en diversas ramas de la física (por ejemplo mecánica continua o dinámica de la red cristalina), estas teorías deben estar inherentemente conectadas entre sí y mostrar ciertos elementos en común en una perspectiva unificada de la física. No obstante las conexiones físicas y diferencias entre las teorías/modelos no se han tratado hasta la fecha, aun siendo de importancia crítica para la mejora de modelos de MT y para el desarrollo integrado de modelos de transformaciones acopladas de desplazamiento-difusión. Por lo tanto, esta tesis comenzó con dos objetivos claros. El primero fue encontrar las conexiones físicas y las diferencias entre los modelos de MT mediante un análisis teórico detallado y simulaciones numéricas. El segundo objetivo fue expandir el modelo de Landau para ser capaz de estudiar MT en policristales, en el caso de transformaciones acopladas de desplazamiento-difusión, y en presencia de dislocaciones. Comenzando con un resumen de los antecedente, en este trabajo se presentan las bases físicas de los modelos actuales de MT. Su capacidad para predecir MT se clarifica mediante el ansis teórico y las simulaciones de la evolución microstructural de MT de cúbicoatetragonal y cúbicoatrigonal en 3D. Este análisis revela que el modelo de Landau con representación irreducible de la deformación transformada es equivalente a la teoría CM-PTMC y al modelo de microelasticidad para predecir los rasgos estáticos durante la MT, pero proporciona una mejor interpretación de los comportamientos dinámicos. Sin embargo, las aplicaciones del modelo de Landau en materiales estructurales están limitadas por su complejidad. Por tanto, el primer resultado de esta tesis es el desarrollo del modelo de Landau nolineal con representación irreducible de deformaciones y de la dinámica de inercia para policristales. La simulación demuestra que el modelo propuesto es consistente fcamente con el CM-PTMC en la descripción estática, y también permite una predicción del diagrama de fases con la clásica forma ’en C’ de los modos de nucleación martensítica activados por la combinación de temperaturas de enfriamiento y las condiciones de tensión aplicada correlacionadas con la transformación de energía de Landau. Posteriomente, el modelo de Landau de MT es integrado con un modelo de transformación de difusión cuantitativa para elucidar la relajación atómica y la difusión de corto alcance de los elementos durante la MT en acero. El modelo de transformaciones de desplazamiento y difusión incluye los efectos de la relajación en borde de grano para la nucleación heterogenea y la evolución espacio-temporal de potenciales de difusión y movilidades químicas mediante el acoplamiento de herramientas de cálculo y bases de datos termo-cinéticos de tipo CALPHAD. El modelo se aplica para estudiar la evolución microstructural de aceros al carbono policristalinos procesados por enfriamiento y partición (Q&P) en 2D. La microstructura y la composición obtenida mediante la simulación se comparan con los datos experimentales disponibles. Los resultados muestran el importante papel jugado por las diferencias en movilidad de difusión entre la fase austenita y martensita en la distibución de carbono en las aceros. Finalmente, un modelo multi-campo es propuesto mediante la incorporación del modelo de dislocación en grano-grueso al modelo desarrollado de Landau para incluir las diferencias morfológicas entre aceros y aleaciones con memoria de forma con la misma ruptura de simetría. La nucleación de dislocaciones, la formación de la martensita ’butterfly’, y la redistribución del carbono después del revenido son bien representadas en las simulaciones 2D del estudio de la evolución de la microstructura en aceros representativos. Con dicha simulación demostramos que incluyendo las dislocaciones obtenemos para dichos aceros, una buena comparación frente a los datos experimentales de la morfología de los bordes de macla, la existencia de austenita retenida dentro de la martensita, etc. Por tanto, basado en un modelo integral y en el desarrollo de códigos durante esta tesis, se ha creado una herramienta de modelización multiescala y multi-campo. Dicha herramienta acopla la termodinámica y la mecánica del continuo en la macroescala con la cinética de difusión y los modelos de campo de fase/Landau en la mesoescala, y también incluye los principios de la cristalografía y de la dinámica de red cristalina en la microescala. ABSTRACT Martensitic transformation (MT), in a narrow sense, is defined as the change of the crystal structure to form a coherent phase, or multi-variant domain structures out from a parent phase with the same composition, by small shuffles or co-operative movements of atoms. Over the past century, MTs have been discovered in different materials from steels to shape memory alloys, ceramics, and smart materials. They lead to remarkable properties such as high strength, shape memory/superelasticity effects or ferroic functionalities including piezoelectricity, electro- and magneto-striction, etc. Various theories/models have been developed, in synergy with development of solid state physics, to understand why MT can generate these rich microstructures and give rise to intriguing properties. Among the well-established theories, the Phenomenological Theory of Martensitic Crystallography (PTMC) is able to predict the habit plane and the orientation relationship between austenite and martensite. The re-interpretation of the PTMC theory within a continuum mechanics framework (CM-PTMC) explains the formation of the multivariant domain structures, while the Landau theory with inertial dynamics unravels the physical origins of precursors and other dynamic behaviors. The crystal lattice dynamics unveils the acoustic softening of the lattice strain waves leading to the weak first-order displacive transformation, etc. Though differing in statics or dynamics due to their origins in different branches of physics (e.g. continuum mechanics or crystal lattice dynamics), these theories should be inherently connected with each other and show certain elements in common within a unified perspective of physics. However, the physical connections and distinctions among the theories/models have not been addressed yet, although they are critical to further improving the models of MTs and to develop integrated models for more complex displacivediffusive coupled transformations. Therefore, this thesis started with two objectives. The first one was to reveal the physical connections and distinctions among the models of MT by means of detailed theoretical analyses and numerical simulations. The second objective was to expand the Landau model to be able to study MTs in polycrystals, in the case of displacive-diffusive coupled transformations, and in the presence of the dislocations. Starting with a comprehensive review, the physical kernels of the current models of MTs are presented. Their ability to predict MTs is clarified by means of theoretical analyses and simulations of the microstructure evolution of cubic-to-tetragonal and cubic-to-trigonal MTs in 3D. This analysis reveals that the Landau model with irreducible representation of the transformed strain is equivalent to the CM-PTMC theory and microelasticity model to predict the static features during MTs but provides better interpretation of the dynamic behaviors. However, the applications of the Landau model in structural materials are limited due its the complexity. Thus, the first result of this thesis is the development of a nonlinear Landau model with irreducible representation of strains and the inertial dynamics for polycrystals. The simulation demonstrates that the updated model is physically consistent with the CM-PTMC in statics, and also permits a prediction of a classical ’C shaped’ phase diagram of martensitic nucleation modes activated by the combination of quenching temperature and applied stress conditions interplaying with Landau transformation energy. Next, the Landau model of MT is further integrated with a quantitative diffusional transformation model to elucidate atomic relaxation and short range diffusion of elements during the MT in steel. The model for displacive-diffusive transformations includes the effects of grain boundary relaxation for heterogeneous nucleation and the spatio-temporal evolution of diffusion potentials and chemical mobility by means of coupling with a CALPHAD-type thermo-kinetic calculation engine and database. The model is applied to study for the microstructure evolution of polycrystalline carbon steels processed by the Quenching and Partitioning (Q&P) process in 2D. The simulated mixed microstructure and composition distribution are compared with available experimental data. The results show that the important role played by the differences in diffusion mobility between austenite and martensite to the partitioning in carbon steels. Finally, a multi-field model is proposed by incorporating the coarse-grained dislocation model to the developed Landau model to account for the morphological difference between steels and shape memory alloys with same symmetry breaking. The dislocation nucleation, the formation of the ’butterfly’ martensite, and the redistribution of carbon after tempering are well represented in the 2D simulations for the microstructure evolution of the representative steels. With the simulation, we demonstrate that the dislocations account for the experimental observation of rough twin boundaries, retained austenite within martensite, etc. in steels. Thus, based on the integrated model and the in-house codes developed in thesis, a preliminary multi-field, multiscale modeling tool is built up. The new tool couples thermodynamics and continuum mechanics at the macroscale with diffusion kinetics and phase field/Landau model at the mesoscale, and also includes the essentials of crystallography and crystal lattice dynamics at microscale.

Aplicación integral en construcción de materiales eco-eficientes con incorporación de cenizas de biomasa

El empleo de biomasa como combustible para la generación de bio-energía va en aumento en la actualidad, debido a su impacto medioambiental nulo en cuanto a las emisiones de CO2. Por lo tanto la generación de cenizas de biomasa, residuo de la producción de esta energía, constituye un problema medioambiental con un claro impacto social y económico. Este tipo de ceniza tiene contenidos en óxidos que la hacen atractiva para su empleo como sustituto parcial del cemento Portland, lo cual proporciona una salida eco-eficiente a este residuo, reduciendo al mismo tiempo la emisión de gases de efecto invernadero asociada a la fabricación del cemento. Esta investigación se centra en el desarrollo de nuevos e innovadores materiales base-cemento eco-eficientes que incorporan ceniza de biomasa para su aplicación integral en construcción. Para ello, se emplea una ceniza de biomasa (CB) procedente de un combustor de lecho fluidizado, cuya biomasa de combustión es principalmente restos de corteza de eucalipto, suministrada por el grupo ENCE-Navia (Asturias). El trabajo desarrollado en la presente tesis doctoral, tiene como primera fase la caracterización de esta ceniza y el análisis de viabilidad de su valorización en materiales base-cemento. Dentro de este análisis, se propone la activación de la ceniza CB mediante tratamiento hidrotermal (TH) en diferentes condiciones de medio activante, temperatura y tiempo de proceso, con el objetivo de favorecer la formación de fases hidratadas que potencien la valorización de la ceniza en el campo de los materiales de construcción. Como fase hidratada de interés se obtiene la fase tobermorita (Ca2.25(Si3O7.5(OH)1.5)(H2O)), precursora del gel C-S-H, responsable del desarrollo de resistencias mecánicas en los materiales base-cemento. El proceso de TH se optimiza para la síntesis más eficiente de esta fase. El estudio posterior de las propiedades mecánicas y micro-estructurales de pastas de cemento eco-eficientes que incorporan la ceniza CB y la ceniza tratada hidrotermalmente, CB-TH, confirma una mayor viabilidad de incorporación de la ceniza CB como sustituto parcial del cemento Portland. Como siguiente paso en el desarrollo de estos innovadores materiales base-cemento eco-eficientes se amplía el estudio multi-escalar de los materiales que incorporan CB mediante diferentes ensayos físico-mecánicos y de durabilidad. Los resultados indican que la presencia de la ceniza de biomasa no tiene efectos negativos sobre las propiedades físicas de los morteros eco-eficientes estudiados. Sin embargo, la adición de CB proporciona una mejor durabilidad del material al producir modificaciones de la microestructura que dificultan el transporte de agentes agresivos. Por otro lado, los morteros con un 10 y 20% de sustitución parcial de cemento por la ceniza de biomasa CB (CB-10 y CB-20) presentan una resistencia a compresión de 53.3 y 50.5 MPa a 28 días de curado, respectivamente. Estos morteros son comparables con un cemento Portland tradicional tipo CEM I de clase de resistencia 42.5 R. Por último, y con el fin de proporcionar la apertura de estos nuevos cementos eco-eficientes al mercado en el campo de los materiales de construcción, se estudian propiedades concretas relacionadas con diferentes tipos de aplicaciones. Concretamente se estudian en detalle las propiedades relativas a la aplicación en baldosas de mortero y los resultados indican unas prestaciones del material eco-eficiente con incorporación de CB similares o mejoradas con respecto al cemento Portland. Se analiza también la viabilidad de aplicación estructural de los cementos eco-eficientes desarrollados mediante el estudio de la adherencia al acero, que resulta similar a la del material de referencia. En cuanto a los resultados de extracción y caracterización de la fase acuosa de los poros, en todas las matrices eco-eficientes se obtiene un pH que garantiza la pasivación de la armadura. Sin embargo, el alto contenido en cloruros de dicha fase acuosa sugiere la conveniencia de realizar un análisis más detallado para la aplicación de los nuevos materiales eco-eficientes en hormigón armado. Se comprueba que todas las matrices que incorporan CB en porcentajes entre un 10 y un 90%, se pueden considerar adecuadas como nuevos materiales de construcción más eco-eficientes en aplicaciones con distintos niveles de exigencias mecánicas y sin problemas ambientales asociados con procesos de lixiviación. Con el presente trabajo de investigación se completan los objetivos iniciales de la tesis, con la obtención de nuevos e innovadores materiales base-cemento eco-eficientes que incorporan cenizas de biomasa (CB) con aplicación integral en el campo de la construcción. ABSTRACT The use of biomass as a fuel for the generation of bio-energy is increasing nowadays, due to its zero environmental impact in terms of CO2 emissions. Therefore the generation of biomass ash, a by-product of this energy, is an environmental problem with a clear social and economic impact. This type of ash contains oxides that make it attractive to be used as a partial replacement of Portland cement, providing an eco-efficient solution to this residue, while reducing the emission of greenhouse gases associated with the production of cement. The present research is focused on the development of new and innovative eco-efficient cement-based materials that incorporate biomass ash for their comprehensive application in construction. For this purpose a biomass ash (CB) is used from a fluidized bed forest combustor mainly fed with the bark of eucalyptus trees, provided by the ENCE-Navia (Asturias) group. The work includes in the first stage the characterization of the raw materials and the analysis of viability of their valorization in cement-based materials. Within this analysis, the activation of the ash is proposed by hydrothermal treatment (HT) in different conditions of activation medium, temperature and process duration, aiming an enhanced formation of hydrated phases to improve the ash valorization in the construction materials field. As an interesting hydrated phase, the tobermorite (Ca2.25(Si3O7.5(OH)1.5)(H2O)) is obtained from the process. This phase is considered as a precursor of the gel C-S-H, responsible for the development of mechanical strength in cement-based materials. HT process is optimized for the most efficient synthesis of tobermorite. The analysis of mechanical and microstructural properties of eco-efficient cement pastes incorporating CB ash and hydrothermally treated ash, CB-TH, confirms an improved viability of incorporation of CB ash as a partial replacement for Portland cement in the case. As a next step in the development of these innovative eco-efficient cement-based materials, a multiscale study of the materials that incorporate CB by different physical-mechanical and durability tests is carried out. The results indicate that the presence of biomass ash does not give rise to negative effects on the physical properties of the eco-efficient mortars analyzed. Nevertheless, the addition of CB produces a better durability performance due to microstructural modifications that hinder the transport of aggressive agents through the material. Moreover, mortars with a 10% and 20% of partial substitution of cement by the CB biomass ash (CB-10 and CB-20) show a compressive resistance of 53.3 and 50.5 MPa at 28 days of curing, respectively. These mortars are comparable to an ordinary Portland cement type CEM I with a resistance class of 42.5R. Finally, and in order to provide the opening of these new eco-efficient cement to the market in the field of construction materials, certain properties specifically related to different types of applications are studied. Among these, the properties concerning the application in mortar tiles are analyzed and the results indicate a similar, or even better performance of the eco-efficient mortar that incorporates CB, with respect to Portland cement. The viability of structural application of the developed eco-efficient cement is also performed considering the study of the adhesion to steel, with results similar to those of the reference material. Regarding the results of extraction and analysis of the aqueous phase of the pores, a pH value guaranteeing reinforcement passivation is obtained for all the eco-efficient matrices. However, high chloride content is obtained suggesting the suitability of a more detailed study to evaluate the application of these new eco-efficient materials in reinforced concrete. It is established that all the matrices incorporating CB in percentages between 10 and 90% may be considered adequate as new more eco-efficient construction materials in applications with different levels of mechanical demand and without environmental problems associated to leaching processes. In this research the initial objectives of the thesis are fulfilled by obtaining new and innovative eco-efficient cement-based materials that incorporate biomass ashes (CB) with comprehensive application in the construction field.

Active control of waves in a cochlear model with subpartitions.

Multiscale asymptotic methods developed previously to study macromechanical wave propagation in cochlear models are generalized here to include active control of a cochlear partition having three subpartitions, the basilar membrane, the reticular lamina, and the tectorial membrane. Activation of outer hair cells by stereocilia displacement and/or by lateral wall stretching result in a frequency-dependent force acting between the reticular lamina and basilar membrane. Wavelength-dependent fluid loads are estimated by using the unsteady Stokes' equations, except in the narrow gap between the tectorial membrane and reticular lamina, where lubrication theory is appropriate. The local wavenumber and subpartition amplitude ratios are determined from the zeroth order equations of motion. A solvability relation for the first order equations of motion determines the subpartition amplitudes. The main findings are as follows: The reticular lamina and tectorial membrane move in unison with essentially no squeezing of the gap; an active force level consistent with measurements on isolated outer hair cells can provide a 35-dB amplification and sharpening of subpartition waveforms by delaying dissipation and allowing a greater structural resonance to occur before the wave is cut off; however, previously postulated activity mechanisms for single partition models cannot achieve sharp enough tuning in subpartitioned models.

Impact of artifact correction methods on R-R interbeat signals to quantifying heart rate variability (HRV) according to linear and nonlinear methods.

In the analysis of heart rate variability (HRV) are used temporal series that contains the distances between successive heartbeats in order to assess autonomic regulation of the cardiovascular system. These series are obtained from the electrocardiogram (ECG) signal analysis, which can be affected by different types of artifacts leading to incorrect interpretations in the analysis of the HRV signals. Classic approach to deal with these artifacts implies the use of correction methods, some of them based on interpolation, substitution or statistical techniques. However, there are few studies that shows the accuracy and performance of these correction methods on real HRV signals. This study aims to determine the performance of some linear and non-linear correction methods on HRV signals with induced artefacts by quantification of its linear and nonlinear HRV parameters. As part of the methodology, ECG signals of rats measured using the technique of telemetry were used to generate real heart rate variability signals without any error. In these series were simulated missing points (beats) in different quantities in order to emulate a real experimental situation as accurately as possible. In order to compare recovering efficiency, deletion (DEL), linear interpolation (LI), cubic spline interpolation (CI), moving average window (MAW) and nonlinear predictive interpolation (NPI) were used as correction methods for the series with induced artifacts. The accuracy of each correction method was known through the results obtained after the measurement of the mean value of the series (AVNN), standard deviation (SDNN), root mean square error of the differences between successive heartbeats (RMSSD), Lomb\'s periodogram (LSP), Detrended Fluctuation Analysis (DFA), multiscale entropy (MSE) and symbolic dynamics (SD) on each HRV signal with and without artifacts. The results show that, at low levels of missing points the performance of all correction techniques are very similar with very close values for each HRV parameter. However, at higher levels of losses only the NPI method allows to obtain HRV parameters with low error values and low quantity of significant differences in comparison to the values calculated for the same signals without the presence of missing points.

Um modelo multiescala concorrente para representar o processo de fissuração do concreto.

Este trabalho propõe uma técnica de modelagem multiescala concorrente do concreto considerando duas escalas distintas: a mesoescala, onde o concreto é modelado como um material heterogêneo, e a macroescala, na qual o concreto é tratado como um material homogêneo. A heterogeneidade da estrutura mesoscópica do concreto é idealizada considerando três fases distintas, compostas pelos agregados graúdos e argamassa (matriz), estes considerados materiais homogêneos, e zona de transição interfacial (ZTI), tratada como a parte mais fraca entre as três fases. O agregado graúdo é gerado a partir de uma curva granulométrica e posicionado na matriz de forma aleatória. Seu comportamento mecânico é descrito por um modelo constitutivo elástico-linear, devido a sua maior resistência quando comparado com as outras duas fases do concreto. Elementos finitos contínuos com alta relação de aspecto em conjunto com um modelo constitutivo de dano são usados para representar o comportamento não linear do concreto, decorrente da iniciação de fissuras na ZTI e posterior propagação para a matriz, dando lugar à formação de macrofissuras. Os elementos finitos de interface com alta relação de aspecto são inseridos entre todos os elementos regulares da matriz e entre os da matriz e agregados, representando a ZTI, tornando-se potenciais caminhos de propagação de fissuras. No estado limite, quando a espessura do elemento de interface tende a zero (h ?0) e, consequentemente, a relação de aspecto tende a infinito, estes elementos apresentam a mesma cinemática da aproximação contínua de descontinuidades fortes (ACDF), sendo apropriados para representar a formação de descontinuidades associados a fissuras, similar aos modelos coesivos. Um modelo de dano à tração é proposto para representar o comportamento mecânico não linear das interfaces, associado à formação de fissuras, ou até mesmo ao eventual fechamento destas. A fim de contornar os problemas causados pela malha de elementos finitos de transição entre as malhas da macro e da mesoescala, que, em geral, apresentam diferenças expressivas 5 de refinamento, utiliza-se uma técnica recente de acoplamento de malhas não conformes. Esta técnica é baseada na definição de elementos finitos de acoplamento (EFAs), os quais são capazes de estabelecer a continuidade de deslocamento entre malhas geradas de forma completamente independentes, sem aumentar a quantidade total de graus de liberdade do problema, podendo ser utilizados tanto para acoplar malhas não sobrepostas quanto sobrepostas. Para tornar possível a análise em multiescala em casos nos quais a região de localização de deformações não pode ser definida a priori, propõe-se uma técnica multiescala adaptativa. Nesta abordagem, usa-se a distribuição de tensões da escala macroscópica como um indicador para alterar a modelagem das regiões críticas, substituindo-se a macroescala pela mesoescala durante a análise. Consequentemente, a malha macroscópica é automaticamente substituída por uma malha mesoscópica, onde o comportamento não linear está na iminência de ocorrer. Testes numéricos são desenvolvidos para mostrar a capacidade do modelo proposto de representar o processo de iniciação e propagação de fissuras na região tracionada do concreto. Os resultados numéricos são comparados com os resultados experimentais ou com aqueles obtidos através da simulação direta em mesoescala (SDM).

Caracterização macro e micromorfológica do solo para compreensão de processos erosivos lineares, topossequência Manacá, São Carlos - SP

A presente pesquisa tem como objetivo compreender a dinâmica de comportamento do solo sob escala macro e micromorfológica visualizados em topossequência, no que concerne aos agentes morfológicos que condicionam e contribuem para deflagração de processos erosivos. A área de estudo está inserida na sub-bacia hidrográfica do Laranja Azeda localizada na região centro-leste do estado de São Paulo, no município de São Carlos/SP, e têm fundamental importância por pertencer à bacia hidrográfica do Ribeirão Feijão, importante manancial urbano para a cidade. O planejamento de uso e ocupação adequados aos fatores físicos que compõe a dinâmica desta paisagem são essenciais visando a conservação e preservação dos recursos hídricos ali existentes, onde a expressiva ocorrência de processos erosivos são objetos de preocupação, já que estes podem causar assoreamento de rios e reservatórios. Utilizando uma metodologia multiescalar para seleção da área de pesquisa em detalhe e compreensão da organização e dinâmica da cobertura pedológica, foram utilizados os procedimentos propostos pela Análise Estrutural da Cobertura Pedológica e conceitos e técnicas da micromorfologia de solos. Verifica-se que a distribuição dos solos na Topossequência Manacá está estritamente correlacionada à transformação vertical do materialde origem em solo, em cuja vertente existe uma diferenciação litológica que condiciona a morfologia diferenciada, tanto em escala macromorfológica quanto micromorfológica. O terço superior e médio da vertente está associado à depósitos colúvio-eluvionaresda Formação Itaqueri, onde desenvolve-se um Latossolo Vermelho Amarelo. Já o terço inferior da vertente corresponde a um solo formado a partir dos arenitos da Formação Botucatu, sendo enquadrado enquanto Neossolo Quartzarênico. Com o auxílio técnicas de análise bidimensional de imagens retiradas das lâminas delgadas de solo, foi possível visualizar e quantificar a macroposidade ao longo da vertente, importante atributo morfológico que controla os fluxos de água e são agentes condicionantes para o desenvolvimento de processos erosivos. Conclui-se que a ocorrência de voçorocas no terço médio inferior da vertente é a materialização em forma de processos erosivos deste comportamento diferencial da massa do solo, onde portanto, na Topossequência Manacá a busca de equilíbrio dinâmico na vertente é induzida pela dinâmica genética evolutiva das formações geológicas que sustentam a paisagem, desencadeada em processos erosivos que tendem a progredir em desequilíbrio, a depender do manejo estabelecido para o local.

Dois ensaios em finanças

In the first chapter, we test some stochastic volatility models using options on the S&P 500 index. First, we demonstrate the presence of a short time-scale, on the order of days, and a long time-scale, on the order of months, in the S&P 500 volatility process using the empirical structure function, or variogram. This result is consistent with findings of previous studies. The main contribution of our paper is to estimate the two time-scales in the volatility process simultaneously by using nonlinear weighted least-squares technique. To test the statistical significance of the rates of mean-reversion, we bootstrap pairs of residuals using the circular block bootstrap of Politis and Romano (1992). We choose the block-length according to the automatic procedure of Politis and White (2004). After that, we calculate a first-order correction to the Black-Scholes prices using three different first-order corrections: (i) a fast time scale correction; (ii) a slow time scale correction; and (iii) a multiscale (fast and slow) correction. To test the ability of our model to price options, we simulate options prices using five different specifications for the rates or mean-reversion. We did not find any evidence that these asymptotic models perform better, in terms of RMSE, than the Black-Scholes model. In the second chapter, we use Brazilian data to compute monthly idiosyncratic moments (expected skewness, realized skewness, and realized volatility) for equity returns and assess whether they are informative for the cross-section of future stock returns. Since there is evidence that lagged skewness alone does not adequately forecast skewness, we estimate a cross-sectional model of expected skewness that uses additional predictive variables. Then, we sort stocks each month according to their idiosyncratic moments, forming quintile portfolios. We find a negative relationship between higher idiosyncratic moments and next-month stock returns. The trading strategy that sells stocks in the top quintile of expected skewness and buys stocks in the bottom quintile generates a significant monthly return of about 120 basis points. Our results are robust across sample periods, portfolio weightings, and to Fama and French (1993)’s risk adjustment factors. Finally, we identify a return reversal of stocks with high idiosyncratic skewness. Specifically, stocks with high idiosyncratic skewness have high contemporaneous returns. That tends to reverse, resulting in negative abnormal returns in the following month.

Neutron-mapping polymer flow: Scattering, flow visualization, and molecular theory

Flows of complex fluids need to be understood at both macroscopic and molecular scales, because it is the macroscopic response that controls the fluid behavior, but the molecular scale that ultimately gives rise to rheological and solid-state properties. Here the flow field of an entangled polymer melt through an extended contraction, typical of many polymer processes, is imaged optically and by small-angle neutron scattering. The dual-probe technique samples both the macroscopic stress field in the flow and the microscopic configuration of the polymer molecules at selected points. The results are compared with a recent tube model molecular theory of entangled melt flow that is able to calculate both the stress and the single-chain structure factor from first principles. The combined action of the three fundamental entangled processes of reptation, contour length fluctuation, and convective constraint release is essential to account quantitatively for the rich rheological behavior. The multiscale approach unearths a new feature: Orientation at the length scale of the entire chain decays considerably more slowly than at the smaller entanglement length.

Integrating JERS-1 imaging radar and elevation models for mapping tropical rainforest communities in far North Queensland, Australia

The Wet Tropics World Heritage Area in Far North Queens- land, Australia consists predominantly of tropical rainforest and wet sclerophyll forest in areas of variable relief. Previous maps of vegetation communities in the area were produced by a labor-intensive combination of field survey and air-photo interpretation. Thus,. the aim of this work was to develop a new vegetation mapping method based on imaging radar that incorporates topographical corrections, which could be repeated frequently, and which would reduce the need for detailed field assessments and associated costs. The method employed G topographic correction and mapping procedure that was developed to enable vegetation structural classes to be mapped from satellite imaging radar. Eight JERS-1 scenes covering the Wet Tropics area for 1996 were acquired from NASDA under the auspices of the Global Rainforest Mapping Project. JERS scenes were geometrically corrected for topographic distortion using an 80 m DEM and a combination of polynomial warping and radar viewing geometry modeling. An image mosaic was created to cover the Wet Tropics region, and a new technique for image smoothing was applied to the JERS texture bonds and DEM before a Maximum Likelihood classification was applied to identify major land-cover and vegetation communities. Despite these efforts, dominant vegetation community classes could only be classified to low levels of accuracy (57.5 percent) which were partly explained by the significantly larger pixel size of the DEM in comparison to the JERS image (12.5 m). In addition, the spatial and floristic detail contained in the classes of the original validation maps were much finer than the JERS classification product was able to distinguish. In comparison to field and aerial photo-based approaches for mapping the vegetation of the Wet Tropics, appropriately corrected SAR data provides a more regional scale, all-weather mapping technique for broader vegetation classes. Further work is required to establish an appropriate combination of imaging radar with elevation data and other environmental surrogates to accurately map vegetation communities across the entire Wet Tropics.

Summarizing level-two topological relations in large spatial datasets

Summarizing topological relations is fundamental to many spatial applications including spatial query optimization. In this article, we present several novel techniques to effectively construct cell density based spatial histograms for range (window) summarizations restricted to the four most important level-two topological relations: contains, contained, overlap, and disjoint. We first present a novel framework to construct a multiscale Euler histogram in 2D space with the guarantee of the exact summarization results for aligned windows in constant time. To minimize the storage space in such a multiscale Euler histogram, an approximate algorithm with the approximate ratio 19/12 is presented, while the problem is shown NP-hard generally. To conform to a limited storage space where a multiscale histogram may be allowed to have only k Euler histograms, an effective algorithm is presented to construct multiscale histograms to achieve high accuracy in approximately summarizing aligned windows. Then, we present a new approximate algorithm to query an Euler histogram that cannot guarantee the exact answers; it runs in constant time. We also investigate the problem of nonaligned windows and the problem of effectively partitioning the data space to support nonaligned window queries. Finally, we extend our techniques to 3D space. Our extensive experiments against both synthetic and real world datasets demonstrate that the approximate multiscale histogram techniques may improve the accuracy of the existing techniques by several orders of magnitude while retaining the cost efficiency, and the exact multiscale histogram technique requires only a storage space linearly proportional to the number of cells for many popular real datasets.

The genetic structure of Australasian green turtles (Chelonia mydas): exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Phi = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.

Training Bayesian networks for image segmentation

We are concerned with the problem of image segmentation in which each pixel is assigned to one of a predefined finite number of classes. In Bayesian image analysis, this requires fusing together local predictions for the class labels with a prior model of segmentations. Markov Random Fields (MRFs) have been used to incorporate some of this prior knowledge, but this not entirely satisfactory as inference in MRFs is NP-hard. The multiscale quadtree model of Bouman and Shapiro (1994) is an attractive alternative, as this is a tree-structured belief network in which inference can be carried out in linear time (Pearl 1988). It is an hierarchical model where the bottom-level nodes are pixels, and higher levels correspond to downsampled versions of the image. The conditional-probability tables (CPTs) in the belief network encode the knowledge of how the levels interact. In this paper we discuss two methods of learning the CPTs given training data, using (a) maximum likelihood and the EM algorithm and (b) emphconditional maximum likelihood (CML). Segmentations obtained using networks trained by CML show a statistically-significant improvement in performance on synthetic images. We also demonstrate the methods on a real-world outdoor-scene segmentation task.