Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Evaluation of sets of oriented and non-oriented receptive fields as local descriptors

Local descriptors are increasingly used for the task of object recognition because of their perceived robustness with respect to occlusions and to global geometrical deformations. We propose a performance criterion for a local descriptor based on the tradeoff between selectivity and invariance. In this paper, we evaluate several local descriptors with respect to selectivity and invariance. The descriptors that we evaluated are Gaussian derivatives up to the third order, gray image patches, and Laplacian-based descriptors with either three scales or one scale filters. We compare selectivity and invariance to several affine changes such as rotation, scale, brightness, and viewpoint. Comparisons have been made keeping the dimensionality of the descriptors roughly constant. The overall results indicate a good performance by the descriptor based on a set of oriented Gaussian filters. It is interesting that oriented receptive fields similar to the Gaussian derivatives as well as receptive fields similar to the Laplacian are found in primate visual cortex.

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequent in practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of the sample, thus producing a transfer of mass from some components to other ones, but preserving the total mass present in the system. This evolution is traditionally modelled as a system of ordinary di erential equations of the mass of each component. However, this kind of evolution can be decomposed into a compositional change, expressed in terms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despite of some subcompositions behaving linearly. The goal is to study the characteristics of such simplicial systems of di erential equa- tions such as linearity and stability. This is performed extracting the compositional dif ferential equations from the mass equations. Then, simplicial derivatives are expressed in coordinates of the simplex, thus reducing the problem to the standard theory of systems of di erential equations, including stability. The characterisation of stability of these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and the associated behaviour of the orbits are the main tools. For a three component system, these orbits can be plotted both in coordinates of the simplex or in a ternary diagram. A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is a radioactive decay

Addressing endogeneity in the relationship between health and employment through a simultaneous equations model

Este documento analiza la relación de doble causalidad entre salud y empleo y su comportamiento dinámico usando datos de Estados Unidos tomados del PSID (Pane Study of Income Dynamics). Este estudio usa dos variables dependientes (Estado de salud auto-reportado y Empleo), las cuales son estimadas usando un modelo probit bivariado para abordar el problema de endegeneidad presente en dicha relación. Los resultados muestran evidencia significativa de la existencia de dicha endogeneidad y del impacto positivo que tiene sobre la probabilidad de ser empleado tener un buen estado de salud y vicesersa, sin embargo, el impacto de la situación de empleo sobre el estado de salud se encuentra que no es significativa.

Polimorfismos del gen Butirilcolinesterasa responsables de reacciones adversas en pacientes consumidores de “cocaína”

La butirilcolinesterasa humana (BChE; EC 3.1.1.8) es una enzima polimórfica sintetizada en el hígado y en el tejido adiposo, ampliamente distribuida en el organismo y encargada de hidrolizar algunos ésteres de colina como la procaína, ésteres alifáticos como el ácido acetilsalicílico, fármacos como la metilprednisolona, el mivacurium y la succinilcolina y drogas de uso y/o abuso como la heroína y la cocaína. Es codificada por el gen BCHE (OMIM 177400), habiéndose identificado más de 100 variantes, algunas no estudiadas plenamente, además de la forma más frecuente, llamada usual o silvestre. Diferentes polimorfismos del gen BCHE se han relacionado con la síntesis de enzimas con niveles variados de actividad catalítica. Las bases moleculares de algunas de esas variantes genéticas han sido reportadas, entre las que se encuentra las variantes Atípica (A), fluoruro-resistente del tipo 1 y 2 (F-1 y F-2), silente (S), Kalow (K), James (J) y Hammersmith (H). En este estudio, en un grupo de pacientes se aplicó el instrumento validado Lifetime Severity Index for Cocaine Use Disorder (LSI-C) para evaluar la gravedad del consumo de “cocaína” a lo largo de la vida. Además, se determinaron Polimorfismos de Nucleótido Simple (SNPs) en el gen BCHE conocidos como responsables de reacciones adversas en pacientes consumidores de “cocaína” mediante secuenciación del gen y se predijo el efecto delos SNPs sobre la función y la estructura de la proteína, mediante el uso de herramientas bio-informáticas. El instrumento LSI-C ofreció resultados en cuatro dimensiones: consumo a lo largo de la vida, consumo reciente, dependencia psicológica e intento de abandono del consumo. Los estudios de análisis molecular permitieron observar dos SNPs codificantes (cSNPs) no sinónimos en el 27.3% de la muestra, c.293A>G (p.Asp98Gly) y c.1699G>A (p.Ala567Thr), localizados en los exones 2 y 4, que corresponden, desde el punto de vista funcional, a la variante Atípica (A) [dbSNP: rs1799807] y a la variante Kalow (K) [dbSNP: rs1803274] de la enzima BChE, respectivamente. Los estudios de predicción In silico establecieron para el SNP p.Asp98Gly un carácter patogénico, mientras que para el SNP p.Ala567Thr, mostraron un comportamiento neutro. El análisis de los resultados permite proponer la existencia de una relación entre polimorfismos o variantes genéticas responsables de una baja actividad catalítica y/o baja concentración plasmática de la enzima BChE y algunas de las reacciones adversas ocurridas en pacientes consumidores de cocaína.

Option pricing in market models driven by telegraph processes with jumps

Esta tesis está dividida en dos partes: en la primera parte se presentan y estudian los procesos telegráficos, los procesos de Poisson con compensador telegráfico y los procesos telegráficos con saltos. El estudio presentado en esta primera parte incluye el cálculo de las distribuciones de cada proceso, las medias y varianzas, así como las funciones generadoras de momentos entre otras propiedades. Utilizando estas propiedades en la segunda parte se estudian los modelos de valoración de opciones basados en procesos telegráficos con saltos. En esta parte se da una descripción de cómo calcular las medidas neutrales al riesgo, se encuentra la condición de no arbitraje en este tipo de modelos y por último se calcula el precio de las opciones Europeas de compra y venta.

Lógica en la toma de decisiones de inversión en el sector de extracción de petróleo y gas natural

Este documento se centra en la presentación de información y análisis de la misma a la hora de establecer la manera en que empresas del sector de extracción de gas natural y generación de energía a base de dicho recurso, toman decisiones en cuanto a inversión, centrándose en la lógica que usan a la hora de emprender este proceso. Esto debido a la constante necesidad de establecer procesos que permitan tomar decisiones más acertadas, incluyendo todas las herramientas posibles para tal fin. La lógica es una de estas herramientas, pues permite encadenar factores con el fin de obtener resultados positivos. Por tal razón, se hace importante conocer el uso de esta herramienta, teniendo en cuentas de qué manera y en que contextos es usada. Con el fin de tener una mayor orientación, este estudio estará centrado en un sector específico, el cual es el de la extracción de petróleo y gas natural. Lo anterior entendiendo la necesidad existente de fundamentación teórica que permita establecer de manera clara la forma apropiada de tomar decisiones en un sector tan diverso y complejo como lo es el mencionado. El contexto empresarial actual exige una visión global, no basada en la lógica lineal causal que hoy se tiene como referencia. El sector de extracción de petróleo y gas natural es un ejemplo particular en cuanto a la manera en cuanto se toman decisiones en inversión, puesto que en su mayoría son empresas de capital intensivo, las cuales mantienen un flujo elevado de recursos monetarios.

Why labor income shares seem to be constant?

The common assumptions that labor income share does not change over time or across countries and that factor income shares are equal to the elasticity of output with respect to factors have had important implications for economic theory. However, there are various theoretical reasons why the elasticity of output with respect to reproducible factors should be correlated with the stage of development. In particular, the behavior of international trade and capital flows and the existence of factor saving innovations imply such a correlation. If this correlation exists and if factor income shares are equal to the elasticity of output with respect to factors then the labor income share must be negatively correlated with the stage of development. We propose an explanation for why labor income share has no correlation with income per capita: the existence of a labor intensive sector which produces non tradable goods.

Energy saving innovations, non-exhaustible sources of energy and long-run: what would happen if we run out of oil?

We formulate and solve a model of factor saving technological improvement considering three factors of production: labor, capital and energy. The productive activities have three main characteristics: first, in order to use capital goods firms need energy; second, there are two sources of energy: non-exhaustible and exhaustible; third, capital goods can be of different qualities and the quality of these goods can be changed along two dimensions -reducing the need of energy or changing the source of energy used in the production process. The economy goes through three stages of development after industrialization. In the first, firms make use of exhaustible energy and the e¢ ciency in the use of energy is constant. In the second stage, as the price of energy grows the e¢ ciency in its use is increased. In the third stage, the price of exhaustible sources is so high that firms have incentives to use non-exhaustible sources of energy. During this stage the price of energy is constant. In this set up, the end of the oil age has level effects on consumption and output but it does not cause the collapse of the economic system.

Adsorció d'or i zinc amb resines impregnades XAD-2. Superfície d'Equilibri, un nou concepte per a l'adsorció

Dins dels processos de recuperació de metalls de dissolucions diluïdes s'ha realitzat un estudi del procés d'extracció d'or i de zinc mitjançant resines amberlite XAD-2 impregnades amb sulfur de triisobutil fosfina (TIBPS) i àcid di-(2-etilhexil) fosfòric (DEHPA) respectivament. S'ha realitzat un estudi de l'equilibri de l'adsorció d'espècies metàl·liques d'aquests metalls amb les resines indicades anteriorment. Amb la metodologia emprada per a la determinació dels punts d'equilibri dels experiments en batch i en columna, s'ha vist que una única isoterma no podia descriure el fenomen global d'equilibri i que en funció de la metodologia emprada s'obtenien isotermes diferents. Es va introduir una nova variable per poder explicar el fenomen observat, i per tant, amb aquesta nova variable l'equació de la isoterma es converteix amb l'equació d'una supèrfície que s'ha definit com a Superfície d'Equilibri. S'han determinat les equacions de les Suprfícies d'Equilibri dels sistemes d'adsorció estudiats (Au(III) TIBPS/XAD-2 i Zn(II) DEHPA/XAD-2) observan una bona coincidència de tots els punts d'equilibri obtinguts sobre la superfície, així com, un bon ajust de totes les isotermes obtingudes en funció de les diferents metodologies emprades sobre les respectives superfícies d'equilibri. Aquest nou concepte generalitza el concepte d'isoterma d'un procés d'adsorció. Fimalment, s'ha plantejat un model matemàtic d'adsorció per a determinar el coeficient efectiu de difusió (De) i el coeficient de transferància de matèria (kf) per ambdós sistemes d'adsorció estudiats mitjançant l'aplicació del model de difusió de sòlid homogeni (HSDM), utilitzant com a condició de contorn en el model la isoterma de Langmuir obtinguda mitjançant els experiments en columna de llit fix i emprant també l'equació obtinguda mitjançant el nou concepte de Superfície d'Equilibri. Els resultats obtinguts són molt satisfactoris, per tant, es pot concloure que la Superfície d'Equilibri és una bona eina per a descriure l'equilibri en els processos d'adsorció d'or i zinc amb les resines amberlite XAD-2 impregnades amb TIBPS i DEHPA respectivament.

Sexual presence and intentional dynamics: deviant and non-deviant sexual self-regulation from the first person stance

Child molesters (n=13) and sexually non-deviant subjects (n=29) were immersed with virtual characters depicting relevant sexual features while their sexual arousal and gaze behaviour were assessed to characterize their sexual preferences and intentional dynamics. Sexual arousal was measured using circumferential penile plethysmography (PPG). Gaze behaviour dynamics was derived from average gaze radial angular deviation (GRAD) and GRAD coefficient of variation (GRADCV). Results show distinct sexual arousal profiles according to sexual preferences and point towards the existence of specific gaze behaviour dynamics guided by sexual intentions. Theoretical interpretations are based on the ecological psychology of J.J. Gibson and the integrated theory of sexual offending (Ward, 2009; Ward & Beech, 2006). Theoretical underpinnings coming from these approaches are advocated as being especially well suited to explain how virtual reality can help probing into child molesters’ phenomenology as lived from the first-person stance.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Understanding mixing efficiency in the oceans. Do the nonlinearities of the equation of state matter?

There exist two central measures of turbulent mixing in turbulent stratified fluids that are both caused by molecular diffusion: 1) the dissipation rate D(APE) of available potential energy APE; 2) the turbulent rate of change Wr, turbulent of background gravitational potential energy GPEr. So far, these two quantities have often been regarded as the same energy conversion, namely the irreversible conversion of APE into GPEr, owing to the well known exact equality D(APE)=Wr, turbulent for a Boussinesq fluid with a linear equation of state. Recently, however, Tailleux (2009) pointed out that the above equality no longer holds for a thermally-stratified compressible, with the ratio ξ=Wr, turbulent/D(APE) being generally lower than unity and sometimes even negative for water or seawater, and argued that D(APE) and Wr, turbulent actually represent two distinct types of energy conversion, respectively the dissipation of APE into one particular subcomponent of internal energy called the "dead" internal energy IE0, and the conversion between GPEr and a different subcomponent of internal energy called "exergy" IEexergy. In this paper, the behaviour of the ratio ξ is examined for different stratifications having all the same buoyancy frequency N vertical profile, but different vertical profiles of the parameter Υ=α P/(ρCp), where α is the thermal expansion coefficient, P the hydrostatic pressure, ρ the density, and Cp the specific heat capacity at constant pressure, the equation of state being that for seawater for different particular constant values of salinity. It is found that ξ and Wr, turbulent depend critically on the sign and magnitude of dΥ/dz, in contrast with D(APE), which appears largely unaffected by the latter. These results have important consequences for how the mixing efficiency should be defined and measured in practice, which are discussed.

Diffusion of the synthetic pyrethroid permethrin into bed-sediments

Bed-sediments are a sink for many micro-organic contaminants in aquatic environments. The impact of toxic contaminants on benthic fauna often depends on their spatial distribution, and the fate of the parent compounds and their metabolites. The distribution of a synthetic pyrethroid, permethrin, a compound known to be toxic to aquatic invertebrates, was studied using river bed-sediments in lotic flume channels. trans/cis-Permethrin diagnostic ratios were used to quantify the photoisomerization of the trans isomer in water. Rates were affected by the presence of sediment particles and colloids when compared to distilled water alone. Two experiments in dark/light conditions with replicate channels were undertaken using natural sediment, previously contaminated with permethrin, to examine the effect of the growth of an algal biofilm at the sediment-water interface on diffusive fluxes of permethrin into the sediment. After 42 days, the bulk water was removed, allowing a fine sectioning of the sediment bed (i.e., every mm down to 5 mm and then 5-10 mm, then every 10 mm down to 50 mm). Permethrin was detected in all cases down to a depth of 5-10 mm, in agreement with estimates by the Millington and Quirk model, and measurements of concentrations in pore water produced a distribution coefficient (K-d) for each section, High K-d's were observed for the top layers, mainly as a result of high organic matter and specific surface area. Concentrations in the algal biofilm measured at the end of the experiment under light conditions, and increases in concentration in the top 1 mm of the sediment, demonstrated that algal/bacterial biofilm material was responsible for high K-d's at the sediment surface, and for the retardation of permethrin diffusion. This specific partition of permethrin to fine sediment particles and algae may enhance its threat to benthic invertebrates. In addition,the analysis of trans/cis-permethrin isomer ratios in sediment showed greater losses of trans-permethrin in the experiment under light conditions, which may have also resulted from enhanced biological activity at the sediment surface.

The effect of aqueous diffusion on the fractionation of chlorine and bromine stable isotopes

Diffusive isotopic fractionation factors are important in order to understand natural processes and have practical application in radioactive waste storage and carbon dioxide sequestration. We determined the isotope fractionation factors and the effective diffusion coefficients of chloride and bromide ions during aqueous diffusion in polyacrylamide gel. Diffusion was determined as functions of temperature, time and concentration. The effect of temperature is relatively large on the diffusion coefficient (D) but only small on isotope fractionation. For chlorine, the ratio, D-35cl/D-37cl varied from 1.00128 +/- 0.00017 (1 sigma) at 2 degrees C to 1.00192 +/- 0.00015 at 80 degrees C. For bromine, D-79Br/D-81Br varied from 1.00098 +/- 0.00009 at 2 degrees C to 1.0064 +/- 0.00013 at 21 degrees C and 1.00078 +/- 0.00018 (1 sigma) at 80 degrees C. There were no significant effects on the isotope fractionation due to concentration. The lack of sensitivity of the diffusive isotope fractionation to anything at the most common temperatures (0 to 30 C) makes it particularly valuable for application to understanding processes in geological environments and an important natural tracer in order to understand fluid transport processes. (C) 2009 Elsevier Ltd. All rights reserved.