Anisotropic seismic tomography with the reversible jump Markov chain Monte Carlo

The established isotropic tomographic models show the features of subduction zones in terms of seismic velocity anomalies, but they are generally subjected to the generation of artifacts due to the lack of anisotropy in forward modelling. There is evidence for the significant influence of seismic anisotropy in the mid-upper mantle, especially for boundary layers like subducting slabs. As consequence, in isotropic models artifacts may be misinterpreted as compositional or thermal heterogeneities. In this thesis project the application of a trans-dimensional Metropolis-Hastings method is investigated in the context of anisotropic seismic tomography. This choice arises as a response to the important limitations introduced by traditional inversion methods which use iterative procedures of optimization of a function object of the inversion. On the basis of a first implementation of the Bayesian sampling algorithm, the code is tested with some cartesian two-dimensional models, and then extended to polar coordinates and dimensions typical of subduction zones, the main focus proposed for this method. Synthetic experiments with increasing complexity are realized to test the performance of the method and the precautions for multiple contexts, taking into account also the possibility to apply seismic ray-tracing iteratively. The code developed is tested mainly for 2D inversions, future extensions will allow the anisotropic inversion of seismological data to provide more realistic imaging of real subduction zones, less subjected to generation of artifacts.

Relativistic Hydrodynamics in Isotropic and Anisotropic Cosmological Backgrounds

A seguito di recenti risultati nel campo dell’Astrofisica, con questo elaborato ci si propone di approfondire il ruolo della viscosità in una classe di modelli cosmologici. La discussione ha lo scopo di generalizzare delle tecniche applicate con successo nello studio dei fluidi ideali a sistemi dove anche la viscosità é un parametro che necessita di essere largamente preso in considerazione. Nello specifico, ci si serve di strumenti computazionali e geometrici e di teoria delle equazioni differenziali alle derivate parziali per comprendere, da un punto di vista matematico, come associare le Equazioni di Einstein ai fluidi in specifici background cosmologici. Questo elaborato parte dallo studio dei fluidi ideali in un background di tipo FLRW al fine di generalizzarlo a fluidi viscosi che scorrono in un background di tipo Bianchi I. Infine si indicano alcuni problemi ancora aperti relativi al caso dissipativo e le possibili strategie risolutive di tali questioni.

Polyphenol Extraction Optimisation From Ceylon Gooseberry (dovyalis Hebecarpa) Pulp.

Originally from Asia, Dovyalis hebecarpa is a dark purple/red exotic berry now also produced in Brazil. However, no reports were found in the literature about phenolic extraction or characterisation of this berry. In this study we evaluate the extraction optimisation of anthocyanins and total phenolics in D. hebecarpa berries aiming at the development of a simple and mild analytical technique. Multivariate analysis was used to optimise the extraction variables (ethanol:water:acetone solvent proportions, times, and acid concentrations) at different levels. Acetone/water (20/80 v/v) gave the highest anthocyanin extraction yield, but pure water and different proportions of acetone/water or acetone/ethanol/water (with >50% of water) were also effective. Neither acid concentration nor time had a significant effect on extraction efficiency allowing to fix the recommended parameters at the lowest values tested (0.35% formic acid v/v, and 17.6 min). Under optimised conditions, extraction efficiencies were increased by 31.5% and 11% for anthocyanin and total phenolics, respectively as compared to traditional methods that use more solvent and time. Thus, the optimised methodology increased yields being less hazardous and time consuming than traditional methods. Finally, freeze-dried D. hebecarpa showed high content of target phytochemicals (319 mg/100g and 1,421 mg/100g of total anthocyanin and total phenolic content, respectively).

Influence Of Substrate Steps On The Catalytic Properties Of Pt Layers: The Ethanol Electrooxidation Reaction.

The ethanol oxidation reaction (EOR) is investigated on Pt/Au(hkl) electrodes. The Au(hkl) single crystals used belong to the [n(111)x(110)] family of planes. Pt is deposited following the galvanic exchange of a previously deposited Cu monolayer using a Pt(2+) solution. Deposition is not epitaxial and the defects on the underlying Au(hkl) substrates are partially transferred to the Pt films. Moreover, an additional (100)-step-like defect is formed, probably as a result of the strain resulting from the Pt and Au lattice mismatch. Regarding the EOR, both vicinal Pt/Au(hkl) surfaces exhibit a behavior that differs from that expected for stepped Pt; for instance, the smaller the step density on the underlying Au substrate, the greater the ability to break the CC bond in the ethanol molecule, as determined by in situ Fourier transform infrared spectroscopy measurements. Also, we found that the acetic acid production is favored as the terrace width decreases, thus reflecting the inefficiency of the surface array to cleave the ethanol molecule.

Hyperhoneycomb Iridate β-li_{2}iro_{3} As A Platform For Kitaev Magnetism.

A complex iridium oxide β-Li_{2}IrO_{3} crystallizes in a hyperhoneycomb structure, a three-dimensional analogue of honeycomb lattice, and is found to be a spin-orbital Mott insulator with J_{eff}=1/2 moment. Ir ions are connected to the three neighboring Ir ions via Ir-O_{2}-Ir bonding planes, which very likely gives rise to bond-dependent ferromagnetic interactions between the J_{eff}=1/2 moments, an essential ingredient of Kitaev model with a spin liquid ground state. Dominant ferromagnetic interaction between J_{eff}=1/2 moments is indeed confirmed by the temperature dependence of magnetic susceptibility χ(T) which shows a positive Curie-Weiss temperature θ_{CW}∼+40  K. A magnetic ordering with a very small entropy change, likely associated with a noncollinear arrangement of J_{eff}=1/2 moments, is observed at T_{c}=38  K. With the application of magnetic field to the ordered state, a large moment of more than 0.35  μ_{B}/Ir is induced above 3 T, a substantially polarized J_{eff}=1/2 state. We argue that the close proximity to ferromagnetism and the presence of large fluctuations evidence that the ground state of hyperhoneycomb β-Li_{2}IrO_{3} is located in close proximity of a Kitaev spin liquid.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

New species of thorny catfish, genus Leptodoras Boulenger (Siluriformes: Doradidae), from Tapajós and Xingu basins, Brazil

A new species of thorny catfish (family Doradidae) is described as a member of genus Leptodoras Boulenger by having a modified oral hood and first gill arch with enlarged accessory lamellae extending well onto medial face of gill filaments. The new species is distinguished by three characteristics unique within Leptodoras: gas bladder moderately sized (not reduced) with simple walls (diverticula absent), and paired bony capsules on anteriormost vertebrae reduced to paired cup-like laminar ossifications separated by triangular septum. The new species was discovered from material recently collected in the Tapajós basin at Serra do Cachimbo and in the headwaters of rio Xingu. The species description is supplemented by a detailed osteological description, discussion on its placement among congeners, and a re-evaluation of the diagnosis of Leptodoras and its relationships with sister taxon Anduzedoras oxyrhynchus. In addition, the distributions of several species of Leptodoras are expanded based on recent collecting efforts and identification of museum specimens not previously studied.

A positional FEM Formulation for geometrical non-linear analysis of shells

This work presents a fully non-linear finite element formulation for shell analysis comprising linear strain variation along the thickness of the shell and geometrically exact description for curved triangular elements. The developed formulation assumes positions and generalized unconstrained vectors as the variables of the problem, not displacements and finite rotations. The full 3D Saint-Venant-Kirchhoff constitutive relation is adopted and, to avoid locking, the rate of thickness variation enhancement is introduced. As a consequence, the second Piola-Kirchhoff stress tensor and the Green strain measure are employed to derive the specific strain energy potential. Curved triangular elements with cubic approximation are adopted using simple notation. Selected numerical simulations illustrate and confirm the objectivity, accuracy, path independence and applicability of the proposed technique.

Evolution of magnetic properties and crystallographic texture in electrical steel with large plastic deformation

Deformation leads to a hardening of steel due to an increase in the density of dislocations and a reduction in their mobility, giving rise to a state of elevated residual stresses in the crystal lattice. In the microstructure, one observes an increase in the contribution of crystalline orientations which are unfavorable to the magnetization, as seen, for example, by a decrease in B(50), the magnetic flux density at a field of 50 A/cm. The present study was carried out with longitudinal strips of fully processed non-oriented (NO) electrical steel, with deformations up to 70% resulting from cold rolling in the longitudinal direction. With increasing plastic deformation, the value of B(50) gradually decreases until it reaches a minimum value, where it remains even for larger deformations. On the other hand, the coercive field H(c) continually increases. Magnetometry results and electron backscatter diffraction results are compared and discussed. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3560895]

Trapezoidal cross-sectional influence on FinFET threshold voltage and corner effects

Fin field effect transistors (FinFETS) are silicon-on-insulator (SOI) transistors with three-dimensional structures. As a result of some fabrication-process limitations (as nonideal anisotropic overetch) some FinFETs have inclined surfaces, which results in trapezoidal cross sections instead of rectangular sections, as expected. This geometric alteration results in some device issues, like carrier profile, threshold voltage, and corner effects. This work analyzes these consequences based on three-dimensional numeric simulation of several dual-gate and triple-gate FinFETs. The simulation results show that the threshold voltage depends on the sidewall inclination angle and that this dependence varies according to the body doping level. The corner effects also depend on the inclination angle and doping level. (C) 2008 The Electrochemical Society.

The Effect of Spatial Scale on Predicting Time Series: A Study on Epidemiological System Identification

A susceptible-infective-recovered (SIR) epidemiological model based on probabilistic cellular automaton (PCA) is employed for simulating the temporal evolution of the registered cases of chickenpox in Arizona, USA, between 1994 and 2004. At each time step, every individual is in one of the states S, I, or R. The parameters of this model are the probabilities of each individual (each cell forming the PCA lattice ) passing from a state to another state. Here, the values of these probabilities are identified by using a genetic algorithm. If nonrealistic values are allowed to the parameters, the predictions present better agreement with the historical series than if they are forced to present realistic values. A discussion about how the size of the PCA lattice affects the quality of the model predictions is presented. Copyright (C) 2009 L. H. A. Monteiro et al.

System Identification and Prediction of Dengue Fever Incidence in Rio de Janeiro

Identification, prediction, and control of a system are engineering subjects, regardless of the nature of the system. Here, the temporal evolution of the number of individuals with dengue fever weekly recorded in the city of Rio de Janeiro, Brazil, during 2007, is used to identify SIS (susceptible-infective-susceptible) and SIR (susceptible-infective-removed) models formulated in terms of cellular automaton (CA). In the identification process, a genetic algorithm (GA) is utilized to find the probabilities of the state transition S -> I able of reproducing in the CA lattice the historical series of 2007. These probabilities depend on the number of infective neighbors. Time-varying and non-time-varying probabilities, three different sizes of lattices, and two kinds of coupling topology among the cells are taken into consideration. Then, these epidemiological models built by combining CA and GA are employed for predicting the cases of sick persons in 2008. Such models can be useful for forecasting and controlling the spreading of this infectious disease.

Estimation and prediction of parameters and breeding values in soybean using REML/BLUP and Least Squares

The aim of this study was to compare REML/BLUP and Least Square procedures in the prediction and estimation of genetic parameters and breeding values in soybean progenies. F(2:3) and F(4:5) progenies were evaluated in the 2005/06 growing season and the F(2:4) and F(4:6) generations derived thereof were evaluated in 2006/07. These progenies were originated from two semi-early, experimental lines that differ in grain yield. The experiments were conducted in a lattice design and plots consisted of a 2 m row, spaced 0.5 m apart. The trait grain yield per plot was evaluated. It was observed that early selection is more efficient for the discrimination of the best lines from the F(4) generation onwards. No practical differences were observed between the least square and REML/BLUP procedures in the case of the models and simplifications for REML/BLUP used here.

System-Size Independence of Directed Flow Measured at the BNL Relativistic Heavy-Ion Collider

We measure directed flow (v(1)) for charged particles in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62.4 GeV, as a function of pseudorapidity (eta), transverse momentum (p(t)), and collision centrality, based on data from the STAR experiment. We find that the directed flow depends on the incident energy but, contrary to all available model implementations, not on the size of the colliding system at a given centrality. We extend the validity of the limiting fragmentation concept to v(1) in different collision systems, and investigate possible explanations for the observed sign change in v(1)(p(t)).

Centrality dependence of charged hadron and strange hadron elliptic flow from root s(NN)=200 GeVAu+Au collisions

We present STAR results on the elliptic flow upsilon(2) Of charged hadrons, strange and multistrange particles from,root s(NN) = 200 GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC). The detailed study of the centrality dependence of upsilon(2) over a broad transverse momentum range is presented. Comparisons of different analysis methods are made in order to estimate systematic uncertainties. To discuss the nonflow effect, we have performed the first analysis Of upsilon(2) with the Lee-Yang zero method for K(S)(0) and A. In the relatively low PT region, P(T) <= 2 GeV/c, a scaling with m(T) - m is observed for identified hadrons in each centrality bin studied. However, we do not observe nu 2(p(T))) scaled by the participant eccentricity to be independent of centrality. At higher PT, 2 1 <= PT <= 6 GeV/c, V2 scales with quark number for all hadrons studied. For the multistrange hadron Omega, which does not suffer appreciable hadronic interactions, the values of upsilon(2) are consistent with both m(T) - m scaling at low p(T) and number-of-quark scaling at intermediate p(T). As a function ofcollision centrality, an increase of p(T)-integrated upsilon(2) scaled by the participant eccentricity has been observed, indicating a stronger collective flow in more central Au+Au collisions.