Landscape structure affects dispersal in the greater white-toothed shrew: inference between genetic and simulated ecological distances

Animal dispersal in a fragmented landscape depends on the complex interaction between landscape structure and animal behavior. To better understand how individuals disperse, it is important to explicitly represent the properties of organisms and the landscape in which they move. A common approach to modelling dispersal includes representing the landscape as a grid of equal sized cells and then simulating individual movement as a correlated random walk. This approach uses a priori scale of resolution, which limits the representation of all landscape features and how different dispersal abilities are modelled. We develop a vector-based landscape model coupled with an object-oriented model for animal dispersal. In this spatially explicit dispersal model, landscape features are defined based on their geographic and thematic properties and dispersal is modelled through consideration of an organism's behavior, movement rules and searching strategies (such as visual cues). We present the model's underlying concepts, its ability to adequately represent landscape features and provide simulation of dispersal according to different dispersal abilities. We demonstrate the potential of the model by simulating two virtual species in a real Swiss landscape. This illustrates the model's ability to simulate complex dispersal processes and provides information about dispersal such as colonization probability and spatial distribution of the organism's path

Structure and Stability of 3He Droplets

We have studied the structure of 3He droplets at zero temperature using a density functional approach plus a configuration interaction calculation in an harmonic oscillator major shell. The most salient feature of open shell drops is that the valence atoms couple their spins to the maximum value compatible with Pauli's principle, building a large magnetic moment. We have determined that 29 atoms constitute the smallest self-bound droplet.

Energy and structure of dilute hard- and soft-sphere gases

The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, such as the variational correlated theory, the Bogoliubov model, and the uniform limit approximation, valid in the weak-interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. Jastrow-type correlation provides a good description of the systems, both hard- and soft-spheres, if the hypernetted chain energy functional is freely minimized and the resulting Euler equation is solved. The study of the soft-sphere potentials confirms the appearance of a dependence of the energy on the shape of the potential at gas paremeter values of x~0.001. For quantities other than the energy, such as the radial distribution functions and the momentum distributions, the dependence appears at any value of x. The occurrence of a maximum in the radial distribution function, in the momentum distribution, and in the excitation spectrum is a natural effect of the correlations when x increases. The asymptotic behaviors of the functions characterizing the structure of the systems are also investigated. The uniform limit approach is very easy to implement and provides a good description of the soft-sphere gas. Its reliability improves when the interaction weakens.

Functional interaction between the estrogen receptor and CTF1: analysis of the vitellogenin gene B1 promoter in yeast.

Eukaryotic gene expression depends on a complex interplay between the transcriptional apparatus and chromatin structure. We report here a yeast model system for investigating the functional interaction between the human estrogen receptor (hER) and CTF1, a member of the CTF/NFI transcription factor family. We show that a CTF1-fusion protein and the hER transactivate a synthetic promoter in yeast in a synergistic manner. This interaction requires the proline-rich transactivation domain of CTF1. When the natural estrogen-dependent vitellogenin B1 promoter is tested in yeast, CTF1 and CTF1-fusion proteins are unable to activate transcription, and no synergy is observed between hER, which activates the B1 promoter, and these factors. Chromatin structure analysis on this promoter reveals positioned nucleosomes at -430 to -270 (+/-20 bp) and at -270 to - 100 (+/-20 bp) relative to the start site of transcription. The positions of the nucleosomes remain unchanged upon hormone-dependent transcriptional activation of the promoter, and the more proximal nucleosome appears to mask the CTF/NFI site located at - 101 to -114. We conclude that a functional interaction of hER with the estrogen response element located upstream of a basal promoter occurs in yeast despite the nucleosomal organization of this promoter, whereas the interaction of CTF1 with its target site is apparently precluded by a nucleosome.

Relativistic mean-field interaction with density-dependent meson-nucleon vertices based on microscopical calculations

Although ab initio calculations of relativistic Brueckner theory lead to large scalar isovector fields in nuclear matter, at present, successful versions of covariant density functional theory neglect the interactions in this channel. A new high-precision density functional DD-MEδ is presented which includes four mesons, σ, ω, δ, and ρ, with density-dependent meson-nucleon couplings. It is based to a large extent on microscopic ab initiocalculations in nuclear matter. Only four of its parameters are determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, in particular the density dependence of the meson-nucleon vertices, are adjusted to nonrelativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass mp*−mn* derived from relativistic Brueckner theory is used to determine the coupling strength of the δ meson and its density dependence.

Electronic structure and properties of Cu2O

The structural and electronic properties of Cu2O have been investigated using the periodic Hartree-Fock method and a posteriori density-functional corrections. The lattice parameter, bulk modulus, and elastic constants have been calculated. The electronic structure of and bonding in Cu2O are analyzed and compared with x-ray photoelectron spectroscopy spectra, showing a good agreement for the valence-band states. To check the quality of the calculated electron density, static structure factors and Compton profiles have been calculated, showing a good agreement with the available experimental data. The effective electron and hole masses have been evaluated for Cu2O at the center of the Brillouin zone. The calculated interaction energy between the two interpenetrated frameworks in the cuprite structure is estimated to be around -6.0 kcal/mol per Cu2O formula. The bonding between the two independent frameworks has been analyzed using a bimolecular model and the results indicate an important role of d10-d10 type interactions between copper atoms.

Complex regional pain syndrome type I affects brain structure in prefrontal and motor cortex.

The complex regional pain syndrome (CRPS) is a rare but debilitating pain disorder that mostly occurs after injuries to the upper limb. A number of studies indicated altered brain function in CRPS, whereas possible influences on brain structure remain poorly investigated. We acquired structural magnetic resonance imaging data from CRPS type I patients and applied voxel-by-voxel statistics to compare white and gray matter brain segments of CRPS patients with matched controls. Patients and controls were statistically compared in two different ways: First, we applied a 2-sample ttest to compare whole brain white and gray matter structure between patients and controls. Second, we aimed to assess structural alterations specifically of the primary somatosensory (S1) and motor cortex (M1) contralateral to the CRPS affected side. To this end, MRI scans of patients with left-sided CRPS (and matched controls) were horizontally flipped before preprocessing and region-of-interest-based group comparison. The unpaired ttest of the "non-flipped" data revealed that CRPS patients presented increased gray matter density in the dorsomedial prefrontal cortex. The same test applied to the "flipped" data showed further increases in gray matter density, not in the S1, but in the M1 contralateral to the CRPS-affected limb which were inversely related to decreased white matter density of the internal capsule within the ipsilateral brain hemisphere. The gray-white matter interaction between motor cortex and internal capsule suggests compensatory mechanisms within the central motor system possibly due to motor dysfunction. Altered gray matter structure in dorsomedial prefrontal cortex may occur in response to emotional processes such as pain-related suffering or elevated analgesic top-down control.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Anisotropic fluid with SU(2) type structure in general relativity: A model of localized matter

A model of anisotropic fluid with three perfect fluid components in interaction is studied. Each fluid component obeys the stiff matter equation of state and is irrotational. The interaction is chosen to reproduce an integrable system of equations similar to the one associated to self-dual SU(2) gauge fields. An extension of the BelinskyZakharov version of the inverse scattering transform is presented and used to find soliton solutions to the coupled Einstein equations. A particular class of solutions that can be interpreted as lumps of matter propagating in empty space-time is examined.

Neutron structure of type-III antifreeze protein allows the reconstruction of AFP-ice interface.

Antifreeze proteins (AFPs) inhibit ice growth at sub-zero temperatures. The prototypical type-III AFPs have been extensively studied, notably by X-ray crystallography, solid-state and solution NMR, and mutagenesis, leading to the identification of a compound ice-binding surface (IBS) composed of two adjacent ice-binding sections, each which binds to particular lattice planes of ice crystals, poisoning their growth. This surface, including many hydrophobic and some hydrophilic residues, has been extensively used to model the interaction of AFP with ice. Experimentally observed water molecules facing the IBS have been used in an attempt to validate these models. However, these trials have been hindered by the limited capability of X-ray crystallography to reliably identify all water molecules of the hydration layer. Due to the strong diffraction signal from both the oxygen and deuterium atoms, neutron diffraction provides a more effective way to determine the water molecule positions (as D(2) O). Here we report the successful structure determination at 293 K of fully perdeuterated type-III AFP by joint X-ray and neutron diffraction providing a very detailed description of the protein and its solvent structure. X-ray data were collected to a resolution of 1.05 Å, and neutron Laue data to a resolution of 1.85 Å with a "radically small" crystal volume of 0.13 mm(3). The identification of a tetrahedral water cluster in nuclear scattering density maps has allowed the reconstruction of the IBS-bound ice crystal primary prismatic face. Analysis of the interactions between the IBS and the bound ice crystal primary prismatic face indicates the role of the hydrophobic residues, which are found to bind inside the holes of the ice surface, thus explaining the specificity of AFPs for ice versus water.

Electronic structure and properties of Cu2O

The structural and electronic properties of Cu2O have been investigated using the periodic Hartree-Fock method and a posteriori density-functional corrections. The lattice parameter, bulk modulus, and elastic constants have been calculated. The electronic structure of and bonding in Cu2O are analyzed and compared with x-ray photoelectron spectroscopy spectra, showing a good agreement for the valence-band states. To check the quality of the calculated electron density, static structure factors and Compton profiles have been calculated, showing a good agreement with the available experimental data. The effective electron and hole masses have been evaluated for Cu2O at the center of the Brillouin zone. The calculated interaction energy between the two interpenetrated frameworks in the cuprite structure is estimated to be around -6.0 kcal/mol per Cu2O formula. The bonding between the two independent frameworks has been analyzed using a bimolecular model and the results indicate an important role of d10-d10 type interactions between copper atoms.

Electronic structure and magnetic interactions of the spin-chain compounds Ca2CuO3 and Sr2CuO3

The results are presented of a combined periodic and cluster model approach to the electronic structure and magnetic interactions in the spin-chain compounds Ca2CuO3 and Sr2CuO3. An extended t-J model is presented that includes in-chain and interchain hopping and magnetic interaction processes with parameters extracted from ab initio calculations. For both compounds, the in-chain magnetic interaction is found to be around -240 meV, larger than in any of the other cuprates reported in the literature. The interchain magnetic coupling is found to be weakly antiferromagnetic, -1 meV. The effective in-chain hopping parameters are estimated to be ~650 meV for both compounds, whereas the value of the interchain hopping parameter is 30 meV for Sr2CuO3 and 40 meV for Ca2CuO3, in line with the larger interchain distance in the former compound. These effective parameters are shown to be consistent with expressions recently suggested for the Néel temperature and the magnetic moments, and with relations that emerge from the t-J model Hamiltonian. Next, we investigate the physical nature of the band gap. Periodic calculations indicate that an interpretation in terms of a charge-transfer insulator is the most appropriate one, in contrast to the suggestion of a covalent correlated insulator recently reported in the literature.

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.

Unravelling the interaction between tectonic and sedimentary processes during lithospheric thinning in the Alpine Tethys margins

The discovery of exhumed continental mantle and hyper-extended crust in present-day magma-poor rifted margins is at the origin of a paradigm shift within the research field of deep-water rifted margins. It opened new questions about the strain history of rifted margins and the nature and composition of sedimentary, crustal and mantle rocks in rifted margins. Thanks to the benefit of more than one century of work in the Alps and access to world-class outcrops preserving the primary relationships between sediments and crustal and mantle rocks from the fossil Alpine Tethys margins, it is possible to link the subsidence history and syn-rift sedimentary evolution with the strain distribution observed in the crust and mantle rocks exposed in the distal rifted margins. In this paper, we will focus on the transition from early to late rifting that is associated with considerable crustal thinning and a reorganization of the rift system. Crustal thinning is at the origin of a major change in the style of deformation from high-angle to low-angle normal faulting which controls basin-architecture, sedimentary sources and processes and the nature of basement rocks exhumed along the detachment faults in the distal margin. Stratigraphic and isotopic ages indicate that this major change occurred in late Sinemurian time, involving a shift of the syn-rift sedimentation toward the distal domain associated with a major reorganization of the crustal structure with exhumation of lower and middle crust. These changes may be triggered by mantle processes, as indicated by the infiltration of MOR-type magmas in the lithospheric mantle, and the uplift of the Brianconnais domain. Thinning and exhumation of the crust and lithosphere also resulted in the creation of new paleogeographic domains, the Proto Valais and Liguria-Piemonte domains. These basins show a complex, 3D temporal and spatial evolution that might have evolved, at least in the case of the Liguria-Piemonte basin, in the formation of an embryonic oceanic crust. The re-interpretation of the rift evolution and the architecture of the distal rifted margins in the Alps have important implications for the understanding of rifted margins worldwide, but also for the paleogeographic reconstruction of the Alpine domain and its subsequent Alpine compressional overprint.