Studying Avalanches in the ground-state of the two-dimensional random field Ising model driven by an external field

We study the exact ground state of the two-dimensional random-field Ising model as a function of both the external applied field B and the standard deviation ¿ of the Gaussian random-field distribution. The equilibrium evolution of the magnetization consists in a sequence of discrete jumps. These are very similar to the avalanche behavior found in the out-of-equilibrium version of the same model with local relaxation dynamics. We compare the statistical distributions of magnetization jumps and find that both exhibit power-law behavior for the same value of ¿. The corresponding exponents are compared.

Optical-model potential for electron and positron elastic scattering by atoms

An optical-model potential for systematic calculations of elastic scattering of electrons and positrons by atoms and positive ions is proposed. The electrostatic interaction is determined from the Dirac-Hartree-Fock self-consistent atomic electron density. In the case of electron projectiles, the exchange interaction is described by means of the local-approximation of Furness and McCarthy. The correlation-polarization potential is obtained by combining the correlation potential derived from the local density approximation with a long-range polarization interaction, which is represented by means of a Buckingham potential with an empirical energy-dependent cutoff parameter. The absorption potential is obtained from the local-density approximation, using the Born-Ochkur approximation and the Lindhard dielectric function to describe the binary collisions with a free-electron gas. The strength of the absorption potential is adjusted by means of an empirical parameter, which has been determined by fitting available absolute elastic differential cross-section data for noble gases and mercury. The Dirac partial-wave analysis with this optical-model potential provides a realistic description of elastic scattering of electrons and positrons with energies in the range from ~100 eV up to ~5 keV. At higher energies, correlation-polarization and absorption corrections are small and the usual static-exchange approximation is sufficiently accurate for most practical purposes.

KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

Lattice-gas model of particles with orientational and positional degrees of freedom: Mean-field treatment

We consider a lattice-gas model of particles with internal orientational degrees of freedom. In addition to antiferromagnetic nearest-neighbor (NN) and next-nearest-neighbor (NNN) positional interactions we also consider NN and NNN interactions arising from the internal state of the particles. The system then shows positional and orientational ordering modes with associated phase transitions at Tp and To temperatures at which long-range positional and orientational ordering are, respectively, lost. We use mean-field techniques to obtain a general approach to the study of these systems. By considering particular forms of the orientational interaction function we study coupling effects between both phase transitions arising from the interplay between orientational and positional degrees of freedom. In mean-field approximation coupling effects appear only for the phase transition taking place at lower temperatures. The strength of the coupling depends on the value of the long-range order parameter that remains finite at that temperature.

Characterization of the HeCo mutant mouse: a new model of subcortical band heterotopia associated with seizures and behavioral deficits.

In human, neuronal migration disorders are commonly associated with developmental delay, mental retardation, and epilepsy. We describe here a new mouse mutant that develops a heterotopic cortex (HeCo) lying in the dorsolateral hemispheric region, between the homotopic cortex (HoCo) and subcortical white matter. Cross-breeding demonstrated an autosomal recessive transmission. Birthdating studies and immunochemistry for layer-specific markers revealed that HeCo formation was due to a transit problem in the intermediate zone affecting both radially and tangentially migrating neurons. The scaffold of radial glial fibers, as well as the expression of doublecortin is not altered in the mutant. Neurons within the HeCo are generated at a late embryonic age (E18) and the superficial layers of the HoCo have a correspondingly lower cell density and layer thickness. Parvalbumin immunohistochemistry showed the presence of gamma-aminobutyric acidergic cells in the HeCo and the mutant mice have a lowered threshold for the induction of epileptic seizures. The mutant showed a developmental delay but, in contrast, memory function was relatively spared. Therefore, this unique mouse model resembles subcortical band heterotopia observed in human. This model represents a new and rare tool to better understand cortical development and to investigate future therapeutic strategies for refractory epilepsy.

Closure of the Monte Carlo dynamical equation in the spherical Sherrington-Kirkpatrick model

We study the analytical solution of the Monte Carlo dynamics in the spherical Sherrington-Kirkpatrick model using the technique of the generating function. Explicit solutions for one-time observables (like the energy) and two-time observables (like the correlation and response function) are obtained. We show that the crucial quantity which governs the dynamics is the acceptance rate. At zero temperature, an adiabatic approximation reveals that the relaxational behavior of the model corresponds to that of a single harmonic oscillator with an effective renormalized mass.

Magnetic relaxation of a one-dimensional model for small particle systems with dipolar interactions: Monte Carlo simulation

Atribution as a function of the time are analyzed and this study leads to a deeper knowledge of the microscopic processes involved in the magnetic relaxation

MODELING OF SOIL LOAD-BEARING CAPACITY AS A FUNCTION OF SOIL MECHANICAL RESISTANCE TO PENETRATION

Estimation of soil load-bearing capacity from mathematical models that relate preconsolidation pressure (σp) to mechanical resistance to penetration (PR) and gravimetric soil water content (U) is important for defining strategies to prevent compaction of agricultural soils. Our objective was therefore to model the σp and compression index (CI) according to the PR (with an impact penetrometer in the field and a static penetrometer inserted at a constant rate in the laboratory) and U in a Rhodic Eutrudox. The experiment consisted of six treatments: no-tillage system (NT); NT with chiseling; and NT with additional compaction by combine traffic (passing 4, 8, 10, and 20 times). Soil bulk density, total porosity, PR (in field and laboratory measurements), U, σp, and CI values were determined in the 5.5-10.5 cm and 13.5-18.5 cm layers. Preconsolidation pressure (σp) and CI were modeled according to PR in different U. The σp increased and the CI decreased linearly with increases in the PR values. The correlations between σp and PR and PR and CI are influenced by U. From these correlations, the soil load-bearing capacity and compaction susceptibility can be estimated by PR readings evaluated in different U.

An agency-oriented model to explain vine-growing specialization in the province of Barcelona (Catalonia, Spain) in the mid-nineteenth century

[cat] En aquest treball presentem un model per explicar el procés d’especialització vitícola assolit als municipis de la província de Barcelona, a mitjans del s. XIX,que cerca entendre com va sorgir històricament un avantatge comparatiu fruit d’un procés que esdevindria un dels punts de partida del procés d’industrialització a Catalunya. Els resultats confirmen els papers jugats pel impuls “Boserupià” de la població en un context d’intensificació de l’ús de la terra, i d’un impuls del mercat “Smithià” en un context d’expansió de la demanda per part de les economies atlàntiques. També es posa de manifest la importància de les dotacions agro-ecològiques i les condicions socioinstitucionals relacionades amb la desigualtat d’ingrés. La difusió de la vinya donà com a resultat unes comunitats rurals menys desiguals fins al 1820, tot i que aquesta desigualtat augmentà de nou a partir d'aleshores.

What uses are mating types? The "developmental switch" model.

Why mating types exist at all is subject to much debate. Among hypotheses, mating types evolved to control organelle transmission during sexual reproduction, or to prevent inbreeding or same-clone mating. Here I review data from a diversity of taxa (including ciliates, algae, slime molds, ascomycetes, and basidiomycetes) to show that the structure and function of mating types run counter the above hypotheses. I argue instead for a key role in triggering developmental switches. Genomes must fulfill a diversity of alternative programs along the sexual cycle. As a haploid gametophyte, an individual may grow vegetatively (through haploid mitoses), or initiate gametogenesis and mating. As a diploid sporophyte, similarly, it may grow vegetatively (through diploid mitoses) or initiate meiosis and sporulation. Only diploid sporophytes (and not haploid gametophytes) should switch on the meiotic program. Similarly, only haploid gametophytes (not sporophytes) should switch on gametogenesis and mating. And they should only do so when other gametophytes are ready to do the same in the neighborhood. As argued here, mating types have evolved primarily to switch on the right program at the right moment.

Cholesterol transport from late endosomes to the Golgi regulates t-SNARE trafficking, assembly, and function

Cholesterol regulates plasma membrane (PM) association and functioning of syntaxin-4 and soluble N-ethylmaleimide-sensitive fusion protein 23 (SNAP23) in the secretory pathway. However, the molecular mechanism and cellular cholesterol pools that determine the localization and assembly of these target membrane SNAP receptors (t-SNAREs) are largely unknown. We recently demonstrated that high levels of annexin A6 (AnxA6) induce accumulation of cholesterol in late endosomes, thereby reducing cholesterol in the Golgi and PM. This leads to an impaired supply of cholesterol needed for cytosolic phospholipase A2 (cPLA2) to drive Golgi vesiculation and caveolin transport to the cell surface. Using AnxA6-overexpressing cells as a model for cellular cholesterol imbalance, we identify impaired cholesterol egress from late endosomes and diminution of Golgi cholesterol as correlating with the sequestration of SNAP23/syntaxin-4 in Golgi membranes. Pharmacological accumulation of late endosomal cholesterol and cPLA2 inhibition induces a similar phenotype in control cells with low AnxA6 levels. Ectopic expression of Niemann-Pick C1 (NPC1) or exogenous cholesterol restores the location of SNAP23 and syntaxin-4 within the PM. Importantly, AnxA6-mediated mislocalization of these t-SNAREs correlates with reduced secretion of cargo via the SNAP23/syntaxin-4¿dependent constitutive exocytic pathway. We thus conclude that inhibition of late endosomal export and Golgi cholesterol depletion modulate t-SNARE localization and functioning along the exocytic pathway.

Evidence for oxygen-island formation on Al(111): Cluster-model theory and x-ray photoelectron spectroscopy

The O 1s x-ray photoelectron spectroscopy spectrum for Al(111)/O at 300 K shows two components whose behavior as a function of time and variation of detection angle are consistent with either (a) a surface species represented by the higher binding-energy (BE) component and a subsurface species represented by the lower BE component, or (b) small close-packed oxygen islands with the interior atoms represented by the lower BE component and the perimeter atoms by the higher BE component. We have modeled both situations using ab initio Hartree-Fock wave functions for clusters of Al and O atoms. For an O atom in a threefold site, it was found that a below-surface position gave a higher O 1s BE than an above-surface position, incompatible with interpretation (a). This change in the O 1s BE could arise because the bond for O to Al may have a more covalent character when the O is below the surface than when it is above the surface. We present evidence consistent with this view. An O adatom island with all the O atoms in threefold sites gives calculated O 1s BE's which are significantly higher for the perimeter O atoms. Further, the results for an isolated O island without the Al substrate present also give higher BE¿s for the perimeter atoms. Both these results are consistent with interpretation (b). Published scanning-tunneling-microscopy data supports the suggestion that the chemisorbed state consists of small, close-packed islands, whereas the presence of two vibrational modes in high-resolution electron-energy-loss spectroscopy data has been interpreted as representing surface and subsurface oxygen atoms. In light of the present results, we suggest that a vibrational interpretation in terms of interior and perimeter adatoms should be considered.

Ab initio valence-bond cluster model for ionic solids: Alkaline-earth oxides

A linear M-O-M (M=metal, O=oxygen) cluster embedded in a Madelung field, and also including the quantum effects of the neighboring ions, is used to represent the alkaline-earth oxides. For this model an ab initio wave function is constructed as a linear combination of Slater determinants written in an atomic orbital basis set, i.e., a valence-bond wave function. Each valence-bond determinant (or group of determinants) corresponds to a resonating valence-bond structure. We have obtained ab initio valence-bond cluster-model wave functions for the electronic ground state and the excited states involved in the optical-gap transitions. Numerical results are reasonably close to the experimental values. Moreover, the model contains the ionic model as a limiting case and can be readily extended and improved.

Effective t-J model Hamiltonian parameters of monolayered cuprate superconductors from ab initio electronic structure calculations

The magnetic coupling constant of selected cuprate superconductor parent compounds has been determined by means of embedded cluster model and periodic calculations carried out at the same level of theory. The agreement between both approaches validates the cluster model. This model is subsequently employed in state-of-the-art configuration interaction calculations aimed to obtain accurate values of the magnetic coupling constant and hopping integral for a series of superconducting cuprates. Likewise, a systematic study of the performance of different ab initio explicitly correlated wave function methods and of several density functional approaches is presented. The accurate determination of the parameters of the t-J Hamiltonian has several consequences. First, it suggests that the appearance of high-Tc superconductivity in existing monolayered cuprates occurs with J/t in the 0.20¿0.35 regime. Second, J/t=0.20 is predicted to be the threshold for the existence of superconductivity and, third, a simple and accurate relationship between the critical temperatures at optimum doping and these parameters is found. However, this quantitative electronic structure versus Tc relationship is only found when both J and t are obtained at the most accurate level of theory.

Obesity and the incidence of chronic diseases: a seemingly unrelated probit approach

Western societies can reduce avoidable mortality and morbidity by better understanding the relationship between obesity and chronic disease. This paper examines the joint determinants of obesity and of heart disease, diabetes, hypertension, and elevated cholesterol. It analyzes a broadly representative Spanish dataset, the 1999 Survey on Disabilities, Impairments and Health Status, using a health production theoretical framework together with a seemingly unrelated probit model approach that controls for unobserved heterogeneity and endogeneity. Its findings provide suggestive evidence of a positive and significant, although specification-dependent, association between obesity and the prevalence of chronic illness