Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

A reverse crystal Engineering approach

The enormous impact of crystal engineering in modern solid state chemistry takes advantage from the connection between a typical basic science field and the word engineering. Regrettably, the engineering aspect of organic or metal organic crystalline materials are limited, so far, to descriptive structural features, sometime entangled with topological aspects, but only rarely with true material design. This should include not only the fabrication and structural description at micro- and nano-scopic level of the solids, but also a proper reverse engineering, a fundamental discipline for engineers. Translated into scientific language, the reverse crystal engineering refers to a dedicated and accurate analysis of how the building blocks contribute to generate a given material property. This would enable a more appropriate design of new crystalline material. We propose here the application of reverse crystal engineering to optical properties of organic and metal organic framework structures, applying the distributed atomic polarizability approach that we have extensively investigated in the past few years[1,2].

Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

HOMO Stabilisation in π-Extended Dibenzotetrathiafulvalene Derivatives for Their Application in Organic Field-Effect Transistors

Three new organic semiconductors, in which either two methoxy units are directly linked to a dibenzotetrathiafulvalene (DB-TTF) central core and a 2,1,3-chalcogendiazole is fused on the one side, or four methoxy groups are linked to the DB-TTF, have been synthesised as active materials for organic field-effect transistors (OFETs). Their electrochemical behaviour, electronic absorption and fluorescence emission as well as photoinduced intramolecular charge transfer were studied. The electron-withdrawing 2,1,3-chalcogendiazole unit significantly affects the electronic properties of these semiconductors, lowering both the HOMO and LUMO energy levels and hence increasing the stability of the semiconducting material. The solution-processed single-crystal transistors exhibit high performance with a hole mobility up to 0.04 cm2 V−1 s−1 as well as good ambient stability.

Charge density analysis for crystal engineering

This review reports on the application of charge density analysis in the field of crystal engineering, which is one of the most growing and productive areas of the entire field of crystallography. While methods to calculate or measure electron density are not discussed in detail, the derived quantities and tools, useful for crystal engineering analyses, are presented and their applications in the recent literature are illustrated. Potential developments and future perspectives are also highlighted and critically discussed.

Dissolution rate spectra data of calcite single crystal and micrite material

The large discrepancy between field and laboratory measurements of mineral reaction rates is a long-standing problem in earth sciences, often attributed to factors extrinsic to the mineral itself. Nevertheless, differences in reaction rate are also observed within laboratory measurements, raising the possibility of intrinsic variations as well. Critical insight is available from analysis of the relationship between the reaction rate and its distribution over the mineral surface. This analysis recognizes the fundamental variance of the rate. The resulting anisotropic rate distributions are completely obscured by the common practice of surface area normalization. In a simple experiment using a single crystal and its polycrystalline counterpart, we demonstrate the sensitivity of dissolution rate to grain size, results that undermine the use of "classical" rate constants. Comparison of selected published crystal surface step retreat velocities (Jordan and Rammensee, 1998) as well as large single crystal dissolution data (Busenberg and Plummer, 1986) provide further evidence of this fundamental variability. Our key finding highlights the unsubstantiated use of a single-valued "mean" rate or rate constant as a function of environmental conditions. Reactivity predictions and long-term reservoir stability calculations based on laboratory measurements are thus not directly applicable to natural settings without a probabilistic approach. Such a probabilistic approach must incorporate both the variation of surface energy as a general range (intrinsic variation) as well as constraints to this variation owing to the heterogeneity of complex material (e.g., density of domain borders). We suggest the introduction of surface energy spectra (or the resulting rate spectra) containing information about the probability of existing rate ranges and the critical modes of surface energy.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.

Elastic (stress-strain) halo associated to ion-induced nano-tracks in lithium niobate: role of crystal anisotropy

The elastic strain/stress fields (halo) around a compressed amorphous nano-track (core) caused by a single high-energy ion impact on LiNbO3 are calculated. A method is developed to approximately account for the effects of crystal anisotropy of LiNbO3 (symmetry 3m) on the stress fields for tracks oriented along the crystal axes (X, Y or Z). It only considers the zero-order (axial) harmonic contribution to the displacement field in the perpendicular plane and uses effective Poisson moduli for each particular orientation. The anisotropy is relatively small; however, it accounts for some differential features obtained for irradiations along the crystallographic axes X, Y and Z. In particular, the irradiation-induced disorder (including halo) and the associated surface swelling appear to be higher for irradiations along the X- or Y-axis in comparison with those along the Z-axis. Other irradiation effects can be explained by the model, e.g. fracture patterns or the morphology of pores after chemical etching of tracks. Moreover, it offers interesting predictions on the effect of irradiation on lattice parameters

Influence of Crystal Quality on the Excitation and Propagation of Surface and Bulk Acoustic Waves in Polycrystalline AlN Films

We investigate the excitation and propagation of acoustic waves in polycrystalline aluminum nitride films along the directions parallel and normal to the c-axis. Longitudinal and transverse propagations are assessed through the frequency response of surface acoustic wave and bulk acoustic wave devices fabricated on films of different crystal qualities. The crystalline properties significantly affect the electromechanical coupling factors and acoustic properties of the piezoelectric layers. The presence of misoriented grains produces an overall decrease of the piezoelectric activity, degrading more severely the excitation and propagation of waves traveling transversally to the c-axis. It is suggested that the presence of such crystalline defects in c-axis-oriented films reduces the mechanical coherence between grains and hinders the transverse deformation of the film when the electric field is applied parallel to the surface.

Electrooptic characterization of tunable cylindrical liquid crystal lenses

In this work, one-dimensional arrays of cylindrical adaptive liquid crystal lenses were manufactured and characterized; and test devices were filled with nematic liquid crystal. Comb interdigitated electrodes were designed as a mask pattern for the control electrode on the top glass substrates. A radial graded refractive index along each microsized lens was achieved by fabricating a layer of high resistance sheet deposited as a control electrode. These tunable lenses were switched by applying amplitude and frequency optimized waveforms on the control electrode. Phase profiles generated by the radial electric field distribution on each lens were measured by a convectional interferometric technique.

Optical Nonlinearities at the Interface between Glass and Liquid Crystal

In this paper, the optical behavior of a nonlinear interface is studied. The nonlinear medium has been a nematic liquid crystal, namely MBBA, and the nonlinear one, glasses of different types (F-10 and F-2) depending on the experimental needs. The anchoring forces at the boundary have been found to inhibit the action of the evanescent field in the case of total internal reflection. Most of observed nonlinearities are due to thermal effects. As a consequence, liquid crystals do not seem to be good candidates for total internal reflection optical bistability.

Optical nonlinearities at the interface between liquid crystal and glass: behavior of the transmitted beam.

During the past years a great interest has been devoted to the study of possible applications of non-linear interfaces, mainly in the field of Optical Bistability. Several papers have been published in this field, and some of them dealing with liquid crystals as non-linear material.

Instabilities in hybrid liquid crystal optical bistable devices

In this paper we report some of the experimental results that can be obtained in the field of hybrid optical bistable devices when liquid crystals are employed as non linear materials. The advantages with respect to other materials are the very low voltages and power needed, compatibles with I.C.'s levels.

Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.