Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature-induced grain coarsening. Despite this large-scale trend, a considerable grain size scatter occurs on the outcrop-scale indicating local influence of second-order effects such as thermal perturbations, fluid flow and second-phase particles. Second-phase particles, whose sizes range from nano- to the micron-scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron-scale second-phase particles, or to samples containing a large number of nano-scale particles. The obtained data set suggests that the second phases induce a temperature-controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D 1⁄4 C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second-phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian-type distributions of the pure samples. In contrast, fluid-enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable d18O and d13C isotope ratios in fluid-affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid-induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1.To summarize, considerable grain size variations of up to one order of magnitude can locally result from second-order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.

Nuevo método de identificación usando la respuesta escalón para sistemas de fase no mínima de segundo orden

Este artículo presenta un nuevo método de identificación para sistemas de fase no mínima basado en la respuesta escalón. El enfoque propuesto provee un modelo aproximado de segundo orden evitando diseños experimentales complejos. El método propuesto es un algoritmo de identificación cerrado basado en puntos característicos de la respuesta escalón de sistemas de fase no mínima de segundo orden. Él es validado usando diferentes modelos lineales. Ellos tienen respuesta inversa entre 3,5% y 38% de la respuesta en régimen permanente. En simulaciones, ha sido demostrado que resultados satisfactorios pueden ser obtenidos usando el procedimiento de identificación propuesto, donde los parámetros identificados presentan errores relativos medios, menores que los obtenidos mediante el método de Balaguer.

Comprendre et maîtriser le passage de type I à type II de puits quantiques d'In(x)Ga(1-x)As(y)Sb(1-y) sur substrat de GaSb

Les antimoniures sont des semi-conducteurs III-V prometteurs pour le développement de dispositifs optoélectroniques puisqu'ils ont une grande mobilité d'électrons, une large gamme spectrale d'émission ou de détection et offrent la possibilité de former des hétérostructures confinées dont la recombinaison est de type I, II ou III. Bien qu'il existe plusieurs publications sur la fabrication de dispositifs utilisant un alliage d'In(x)Ga(1-x)As(y)Sb(1-y) qui émet ou détecte à une certaine longueur d'onde, les détails, à savoir comment sont déterminés les compositions et surtout les alignements de bande, sont rarement explicites. Très peu d'études fondamentales sur l'incorporation d'indium et d'arsenic sous forme de tétramères lors de l'épitaxie par jets moléculaires existent, et les méthodes afin de déterminer l'alignement des bandes des binaires qui composent ces alliages donnent des résultats variables. Un modèle a été construit et a permis de prédire l'alignement des bandes énergétiques des alliages d'In(x)Ga(1-x)As(y)Sb(1-y) avec celles du GaSb pour l'ensemble des compositions possibles. Ce modèle tient compte des effets thermiques, des contraintes élastiques et peut aussi inclure le confinement pour des puits quantiques. De cette manière, il est possible de prédire la transition de type de recombinaison en fonction de la composition. Il est aussi montré que l'indium ségrègue en surface lors de la croissance par épitaxie par jets moléculaires d'In(x)Ga(1-x)Sb sur GaSb, ce qui avait déjà été observé pour ce type de matériau. Il est possible d'éliminer le gradient de composition à cette interface en mouillant la surface d'indium avant la croissance de l'alliage. L'épaisseur d'indium en surface dépend de la température et peut être évaluée par un modèle simple simulant la ségrégation. Dans le cas d'un puits quantique, il y aura une seconde interface GaSb sur In(x)Ga(1-x)Sb où l'indium de surface ira s'incorporer. La croissance de quelques monocouches de GaSb à basse température immédiatement après la croissance de l'alliage permet d'incorporer rapidement ces atomes d'indium et de garder la seconde interface abrupte. Lorsque la composition d'indium ne change plus dans la couche, cette composition correspond au rapport de flux d'atomes d'indium sur celui des éléments III. L'arsenic, dont la source fournit principalement des tétramères, ne s'incorpore pas de la même manière. Les tétramères occupent deux sites en surface et doivent interagir par paire afin de créer des dimères d'arsenic. Ces derniers pourront alors être incorporés dans l'alliage. Un modèle de cinétique de surface a été élaboré afin de rendre compte de la diminution d'incorporation d'arsenic en augmentant le rapport V/III pour une composition nominale d'arsenic fixe dans l'In(x)Ga(1-x)As(y)Sb(1-y). Ce résultat s'explique par le fait que les réactions de deuxième ordre dans la décomposition des tétramères d'arsenic ralentissent considérablement la réaction d'incorporation et permettent à l'antimoine d'occuper majoritairement la surface. Cette observation montre qu'il est préférable d'utiliser une source de dimères d'arsenic, plutôt que de tétramères, afin de mieux contrôler la composition d'arsenic dans la couche. Des puits quantiques d'In(x)Ga(1-x)As(y)Sb(1-y) sur GaSb ont été fabriqués et caractérisés optiquement afin d'observer le passage de recombinaison de type I à type II. Cependant, celui-ci n'a pas pu être observé puisque les spectres étaient dominés par un niveau énergétique dans le GaSb dont la source n'a pu être identifiée. Un problème dans la source de gallium pourrait être à l'origine de ce défaut et la résolution de ce problème est essentielle à la continuité de ces travaux.

Análisis teórico del contacto plasma-superficie y sus aplicaciones industriales

Actualmente, la física de plasmas constituye una parte importante de la investigación en física que está siendo desarrollada. Su campo de aplicación varía desde el estudio de plasmas interestelares y cósmicos, como las estrellas, las nebulosas, el medio intergaláctico, etc.; hasta aplicaciones más terrenales como la producción de microchips o los dispositivos de iluminación. Resulta particularmente interesante el estudio del contacto de una superficie metálica con un plasma. Siendo la razón que, la dinámica de la interfase formada entre un plasma imperturbado y una superficie metálica, resulta de gran importancia cuando se trata de estudiar problemas como: la implantación iónica en una oblea de silicio, el grabado por medio de plasmas, la carga de una aeronave cuando atraviesa la ionosfera y la diagnosis de plasmas mediante sondas de Langmuir. El uso de las sondas de Langmuir está extendido a través de multitud de aplicaciones tecnológicas e industriales como método de diagnosis de plasmas. Algunas de estas aplicaciones han sido mencionadas justo en el párrafo anterior. Es más, su uso también es muy popular en la investigación en física de plasmas, por ser una de las pocas técnicas de diagnosis que proporciona información local sobre el plasma. El equipamiento donde es habitualmente implementado varía desde plasmas de laboratorio de baja temperatura hasta plasmas de fusión en dispositivos como tokamaks o stellerators. La geometría más popular de este tipo de sondas es cilíndrica, y la principal magnitud que se usa para diagnosticar el plasma es la corriente recogida por la sonda cuando se encuentra polarizada a un cierto potencial. Existe un interes especial en diagnosticar por medio de la medida de la corriente iónica recogida por la sonda, puesto que produce una perturbación muy pequeña del plasma en comparación con el uso de la corriente electrónica. Dada esta popularidad, no es de extrañar que grandes esfuerzos se hayan realizado en la consecución de un modelo teórico que explique el comportamiento de una sonda de Langmuir inmersa en un plasma. Hay que remontarse a la primera mitad del siglo XX para encontrar las primeras teorías que permiten diagnosticar parámetros del plasma mediante la medida de la corriente iónica recogida por la sonda de Langmuir. Desde entonces, las mejoras en estos modelos y el desarrollo de otros nuevos ha sido una constante en la investigación en física de plasmas. No obstante, todavía no está claro como los iones se aproximan a la superficie de la sonda. Las dos principales, a la par que opuestas, aproximaciones al problema que están ampliamente aceptadas son: la radial y la orbital; siendo el problema que ambas predicen diferentes valores para la corriente iónica. Los experimentos han arrojado resultados de acuerdo con ambas teorías, la radial y la orbital; y lo que es más importante, una transición entre ambos ha sido recientemente observada. La mayoría de los logros conseguidos a la hora de comprender como los iones caen desde el plasma hacia la superficie de la sonda, han sido llevados a cabo en el campo de la dinámica de fluidos o la teoría cinética. Por otra parte, este problema puede ser abordado mediante el uso de simulaciones de partículas. La principal ventaja de las simulaciones de partículas sobre los modelos de fluidos o cinéticos es que proporcionan mucha más información sobre los detalles microscópicos del movimiento de las partículas, además es relativamente fácil introducir interacciones complejas entre las partículas. No obstante, estas ventajas no se obtienen gratuitamente, ya que las simulaciones de partículas requieren grandísimos recursos. Por esta razón, es prácticamente obligatorio el uso de técnicas de procesamiento paralelo en este tipo de simulaciones. El vacío en el conocimiento de las sondas de Langmuir, es el que motiva nuestro trabajo. Nuestra aproximación, y el principal objetivo de este trabajo, ha sido desarrollar una simulación de partículas que nos permita estudiar el problema de una sonda de Langmuir inmersa en un plasma y que está negativamente polarizada con respecto a éste. Dicha simulación nos permitiría estudiar el comportamiento de los iones en los alrededores de una sonda cilíndrica de Langmuir, así como arrojar luz sobre la transición entre las teorías radiales y orbitales que ha sido observada experimentalmente. Justo después de esta sección introductoria, el resto de la tesis está dividido en tres partes tal y como sigue: La primera parte está dedicada a establecer los fundamentos teóricos de las sondas de Langmuir. En primer lugar, se realiza una introducción general al problema y al uso de sondas de Langmuir como método de diagnosis de plasmas. A continuación, se incluye una extensiva revisión bibliográfica sobre las diferentes teorías que proporcionan la corriente iónica recogida por una sonda. La segunda parte está dedicada a explicar los detalles de las simulaciones de partículas que han sido desarrolladas a lo largo de nuestra investigación, así como los resultados obtenidos con las mismas. Esta parte incluye una introducción sobre la teoría que subyace el tipo de simulaciones de partículas y las técnicas de paralelización que han sido usadas en nuestros códigos. El resto de esta parte está dividido en dos capítulos, cada uno de los cuales se ocupa de una de las geometrías consideradas en nuestras simulaciones (plana y cilíndrica). En esta parte discutimos también los descubrimientos realizados relativos a la transición entre el comportamiento radial y orbital de los iones en los alrededores de una sonda cilíndrica de Langmuir. Finalmente, en la tercera parte de la tesis se presenta un resumen del trabajo realizado. En este resumen, se enumeran brevemente los resultados de nuestra investigación y se han incluido algunas conclusiones. Después de esto, se enumeran una serie de perspectivas futuras y extensiones para los códigos desarrollados.

The bioconcentration factor of pesticides in potatoes.

We presented a model that estimates the bioconcentration factor (BCF) of pesticides in potatoes supposing that the pesticide in the soil solution is absorbed by the potato by passive diffusion, following Fick?s second law. The pesticides in the model are nonionic organic substances, traditionally used in potato crops that degrade in the soil according to a first order kinetic equation. This presents an expression that relates BCF with the pesticide elimination rate by the potato, with the pesticide accumulation rate within the potato, with the rate of growth of the potato and with the pesticide degradation rate in the soil. BCF was estimated supposing steady state equilibrium of the quotient between the pesticide concentration in the potato and the pesticide concentration in the soil solution. It is suggested that a negative correlation exists between the pesticide BCF and the soil sorption partition coefficient. The model was built based on the work of Trapp et al. (2007), [Diffusion of PAH in Potato and Carrot Slices and Applications for a Potato Model] in which an expression to calculate the diffusivity of persistent organic substances in potatoes is presented. The model consists in adding to the expression of Trapp et al. (2007) the hypothesis that the pesticide degrades in the soil. The value of BCF suggests which pesticides should be monitored in potatoes.

Phase transitions in gene networks evolved under different selection rules

Mathematical models of gene regulation are a powerful tool for understanding the complex features of genetic control. While various modeling efforts have been successful at explaining gene expression dynamics, much less is known about how evolution shapes the structure of these networks. An important feature of gene regulatory networks is their stability in response to environmental perturbations. Regulatory systems are thought to have evolved to exist near the transition between stability and instability, in order to have the required stability to environmental fluctuations while also being able to achieve a wide variety of functions (corresponding to different dynamical patterns). We study a simplified model of gene network evolution in which links are added via different selection rules. These growth models are inspired by recent work on `explosive' percolation which shows that when network links are added through competitive rather than random processes, the connectivity phase transition can be significantly delayed, and when it is reached, it appears to be first order (discontinuous, e.g., going from no failure at all to large expected failure) instead of second order (continuous, e.g., going from no failure at all to very small expected failure). We find that by modifying the traditional framework for networks grown via competitive link addition to capture how gene networks evolve to avoid damage propagation, we also see significant delays in the transition that depend on the selection rules, but the transitions always appear continuous rather than `explosive'.

Methods for Modeling and Analyzing Concurrent Software

Concurrent software executes multiple threads or processes to achieve high performance. However, concurrency results in a huge number of different system behaviors that are difficult to test and verify. The aim of this dissertation is to develop new methods and tools for modeling and analyzing concurrent software systems at design and code levels. This dissertation consists of several related results. First, a formal model of Mondex, an electronic purse system, is built using Petri nets from user requirements, which is formally verified using model checking. Second, Petri nets models are automatically mined from the event traces generated from scientific workflows. Third, partial order models are automatically extracted from some instrumented concurrent program execution, and potential atomicity violation bugs are automatically verified based on the partial order models using model checking. Our formal specification and verification of Mondex have contributed to the world wide effort in developing a verified software repository. Our method to mine Petri net models automatically from provenance offers a new approach to build scientific workflows. Our dynamic prediction tool, named McPatom, can predict several known bugs in real world systems including one that evades several other existing tools. McPatom is efficient and scalable as it takes advantage of the nature of atomicity violations and considers only a pair of threads and accesses to a single shared variable at one time. However, predictive tools need to consider the tradeoffs between precision and coverage. Based on McPatom, this dissertation presents two methods for improving the coverage and precision of atomicity violation predictions: 1) a post-prediction analysis method to increase coverage while ensuring precision; 2) a follow-up replaying method to further increase coverage. Both methods are implemented in a completely automatic tool.

Credit spreads and monetary policy in a small open economy: an analysis of the colombian economy

This paper estimates Bejarano and Charry (2014)’s small open economy with financial frictions model for the Colombian economy using Bayesian estimation techniques. Additionally, I compute the welfare gains of implementing an optimal response to credit spreads into an augmented Taylor rule. The main result is that a reaction to credit spreads does not imply significant welfare gains unless the economic disturbances increases its volatility, like the disruption implied by a financial crisis. Otherwise its impact over the macroeconomic variables is null.

Customer Loyalty toward grocery retailing service: a multidimensional approach to costumer perceived value

Structured abstract Purpose: To deepen, in grocery retail context, the roles of consumer perceived value and consumer satisfaction, as antecedents’ dimensions of customer loyalty intentions. Design/Methodology/approach: Also employing a short version (12-items) of the original 19-item PERVAL scale of Sweeney & Soutar (2001), a structural equation modeling approach was applied to investigate statistical properties of the indirect influence on loyalty of a reflective second order customer perceived value model. The performance of three alternative estimation methods was compared through bootstrapping techniques. Findings: Results provided i) support for the use of the short form of the PERVAL scale in measuring consumer perceived value; ii) the influence of the four highly correlated independent latent predictors on satisfaction was well summarized by a higher-order reflective specification of consumer perceived value; iii) emotional and functional dimensions were determinants for the relationship with the retailer; iv) parameter’s bias with the three methods of estimation was only significant for bootstrap small sample sizes. Research limitations:/implications: Future research is needed to explore the use of the short form of the PERVAL scale in more homogeneous groups of consumers. Originality/value: Firstly, to indirectly explain customer loyalty mediated by customer satisfaction it was adopted a recent short form of PERVAL scale and a second order reflective conceptualization of value. Secondly, three alternative estimation methods were used and compared through bootstrapping and simulation procedures.

Constructive development of the solutions of linear equations in introductory ordinary differential equations

The solution of linear ordinary differential equations (ODEs) is commonly taught in first year undergraduate mathematics classrooms, but the understanding of the concept of a solution is not always grasped by students until much later. Recognising what it is to be a solution of a linear ODE and how to postulate such solutions, without resorting to tables of solutions, is an important skill for students to carry with them to advanced studies in mathematics. In this study we describe a teaching and learning strategy that replaces the traditional algorithmic, transmission presentation style for solving ODEs with a constructive, discovery based approach where students employ their existing skills as a framework for constructing the solutions of first and second order linear ODEs. We elaborate on how the strategy was implemented and discuss the resulting impact on a first year undergraduate class. Finally we propose further improvements to the strategy as well as suggesting other topics which could be taught in a similar manner.

Developing an integrative model for understanding innovation adoption

We develop and test a theoretically-based integrative model of organizational innovation adoption. Confirmatory factor analyses using responses from 134 organizations showed that the hypothesized second-order model was a better fit to the data than the traditional model of independent factors. Furthermore, although not all elements were significant, the hypothesized model fit adoption better than the traditional model.

Decentred optical axes and aberrations along principal visual field meridians

Suggestions that peripheral imagery may affect the development of refractive error have led to interest in the variation in refraction and aberration across the visual field. It is shown that, if the optical system of the eye is rotationally symmetric about an optical axis which does not coincide with the visual axis, measurements of refraction and aberration made along the horizontal and vertical meridians of the visual field will show asymmetry about the visual axis. The departures from symmetry are modelled for second-order aberrations, refractive components and third-order coma. These theoretical results are compared with practical measurements from the literature. The experimental data support the concept that departures from symmetry about the visual axis in the measurements of crossed-cylinder astigmatism J45 and J180 are largely explicable in terms of a decentred optical axis. Measurements of the mean sphere M suggest, however, that the retinal curvature must differ in the horizontal and vertical meridians.

Numerical methods for fractional partial differential equations with Riesz space fractional derivatives

In this paper, we consider the numerical solution of a fractional partial differential equation with Riesz space fractional derivatives (FPDE-RSFD) on a finite domain. Two types of FPDE-RSFD are considered: the Riesz fractional diffusion equation (RFDE) and the Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second-order space derivative with the Riesz fractional derivative of order αset membership, variant(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first-order and second-order space derivatives with the Riesz fractional derivatives of order βset membership, variant(0,1) and of order αset membership, variant(1,2], respectively. Firstly, analytic solutions of both the RFDE and RFADE are derived. Secondly, three numerical methods are provided to deal with the Riesz space fractional derivatives, namely, the L1/L2-approximation method, the standard/shifted Grünwald method, and the matrix transform method (MTM). Thirdly, the RFDE and RFADE are transformed into a system of ordinary differential equations, which is then solved by the method of lines. Finally, numerical results are given, which demonstrate the effectiveness and convergence of the three numerical methods.

Fundamental solution and discrete random walk model for time-space fractional diffusion equation

In this paper, we consider a time-space fractional diffusion equation of distributed order (TSFDEDO). The TSFDEDO is obtained from the standard advection-dispersion equation by replacing the first-order time derivative by the Caputo fractional derivative of order α∈(0,1], the first-order and second-order space derivatives by the Riesz fractional derivatives of orders β 1∈(0,1) and β 2∈(1,2], respectively. We derive the fundamental solution for the TSFDEDO with an initial condition (TSFDEDO-IC). The fundamental solution can be interpreted as a spatial probability density function evolving in time. We also investigate a discrete random walk model based on an explicit finite difference approximation for the TSFDEDO-IC.

Describing ocular aberrations with wavefront vergence maps

A common optometric problem is to specify the eye’s ocular aberrations in terms of Zernike coefficients and to reduce that specification to a prescription for the optimum sphero-cylindrical correcting lens. The typical approach is first to reconstruct wavefront phase errors from measurements of wavefront slopes obtained by a wavefront aberrometer. This paper applies a new method to this clinical problem that does not require wavefront reconstruction. Instead, we base our analysis of axial wavefront vergence as inferred directly from wavefront slopes. The result is a wavefront vergence map that is similar to the axial power maps in corneal topography and hence has a potential to be favoured by clinicians. We use our new set of orthogonal Zernike slope polynomials to systematically analyse details of the vergence map analogous to Zernike analysis of wavefront maps. The result is a vector of slope coefficients that describe fundamental aberration components. Three different methods for reducing slope coefficients to a spherocylindrical prescription in power vector forms are compared and contrasted. When the original wavefront contains only second order aberrations, the vergence map is a function of meridian only and the power vectors from all three methods are identical. The differences in the methods begin to appear as we include higher order aberrations, in which case the wavefront vergence map is more complicated. Finally, we discuss the advantages and limitations of vergence map representation of ocular aberrations.