Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

Physicochemical characterization of nanoparticles formed between DNA and phosphorylcholine substituted chitosans

The interactions between phosphorylcholine-substituted chitosans (PC-CH) and calf-thymus DNA (ct-DNA) were investigated focusing on the effects of the charge ratio, the pH, and phosphorylcholine content on the size and stability of the complexes using the ethidium bromide fluorescence assay, gel electrophoresis, dynamic light scattering. and fluorescence microscopy. The size and colloidal stability of deacetylated chitosan (CH/DNA) and PC-CH/DNA complexes were strongly dependent on phosphorylcholine content, charge ratios, and pH. The interaction strengths were evaluated from ethidium bromide fluorescence, and at N/P ratios higher than 5.0, no DNA release was observed in any synthesized PC-CH/DNA polyplexes by gel electrophoresis. The PC-CH/DNA polyplexes exhibited a higher resistance to aggregation compared to deacetylated chitosan (CH) at neutral pH. At low pH values highly charged chitosan and its phosphorylcholine derivatives had strong binding affinity with DNA, whereas at higher pH Values CH formed large aggregates and only C-CH derivatives were able to form small nanoparticles with hydrodynamic radii varying from 100 to 150 nm. Nanoparticles synthesized at low ionic strength with PC-CH derivatives containing moderate degrees of substitution (DS = 20% and 40%) remained stable for weeks. Photomicroscopies also confirmed that rhodamine-labeled PC(40)CH derivative nanoparticles presented higher colloidal stability than those synthesized using deacetylated chitosan. Accordingly, due to their improved physicochemical properties these phosphorylcholine-modified chitosans provide new perspectives for controlling the properties of polyplexes. (C) 2009 Elsevier Inc. All rights reserved.

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

A regra de Bent contextualiza a força da ligação de hidrogênio em clusters trimoleculares

This article presents a theoretical study of the molecular properties of trimolecular clusters of CnHm∙∙∙HCN∙∙∙HX formed by the ϖ∙∙∙H and n∙∙∙H hydrogen bonds. The interaction strengths of these interactions are in line with the variations in s-character, and independently, the red-shift rise whether stronger or weaker bound systems are carried out. This behavior was justified via NBO analysis and supported by Bent´s rule, wherein the greater variations in s-character of X are in good agreement with larger red-shifts and vice-versa. To conclude, the refinement of the supermolecule approach and NBO binding energies also corroborate in this regard.

Anisotropy of magnetoviscous effect in structure-forming ferrofluids

The magnetoviscous effect, change in viscosity with change in magnetic field strength, and the anisotropy of magnetoviscous effect, change in viscosity with orientation of magnetic field, have been a focus of interest since four decades. A satisfactory understanding of the microscopic origin of anisotropy of magnetoviscous effect in magnetic fluids is still a matter of debate and a field of intense research. Here, we present an extensive simulation study to understand the relation between the anisotropy of magnetoviscous effect and the underlying change in micro-structures of ferrofluids. Our results indicate that field-induced chain-like structures respond very differently depending on their orientation relative to the direction of an externally applied shear flow, which leads to a pronounced anisotropy of viscosity. In this work, we focus on three exemplary values of dipolar interaction strengths which correspond to weak, intermediate and strong interactions between dipolar colloidal particles. We compare our simulation results with an experimental study on cobalt-based ferrofluids as well as with an existing theoretical model called the chain model. A non-monotonic behaviour in the anisotropy of magnetoviscous effect is observed with increasing dipolar interaction strength and is explained in terms of micro-structure formation.

The role of predator overlap in the robustness and extinction of a four species predator-prey network

Predators and preys often form species networks with asymmetric patterns of interaction. We study the dynamics of a four species network consisting of two weakly connected predator-prey pairs. We focus our analysis on the effects of the cross interaction between the predator of the first pair and the prey of the second pair. This is an example where the predator overlap, which is the proportion of predators that a given prey shares with other preys, is not uniform across the network due to asymmetries in patterns of interaction. We explore the behavior of the system under different interaction strengths and study the dynamics of survival and extinction. In particular, we consider situations in which the four species have initial populations lower than their long-term equilibrium, simulating catastrophic situations in which their abundances are reduced due to human action or environmental change. We show that, under these reduced initial conditions, and depending on the strength of the cross interaction, the populations tend to oscillate before re-equilibrating, disturbing the community equilibrium and sometimes reaching values that are only a small fraction of the equilibrium population, potentially leading to their extinction. We predict that, contrary to one`s intuition, the most likely scenario is the extinction of the less predated preys. (C) 2010 Elsevier B.V. All rights reserved.

Dynamical self-stabilization of the Mott insulator: Time evolution of the density and entanglement entropy of out-of-equilibrium cold fermion gases

The time evolution of the out-of-equilibrium Mott insulator is investigated numerically through calculations of space-time-resolved density and entropy profiles resulting from the release of a gas of ultracold fermionic atoms from an optical trap. For adiabatic, moderate and sudden switching-off of the trapping potential, the out-of-equilibrium dynamics of the Mott insulator is found to differ profoundly from that of the band insulator and the metallic phase, displaying a self-induced stability that is robust within a wide range of densities, system sizes and interaction strengths. The connection between the entanglement entropy and changes of phase, known for equilibrium situations, is found to extend to the out-of-equilibrium regime. Finally, the relation between the system`s long time behavior and the thermalization limit is analyzed. Copyright (C) EPLA, 2011

Thomas-Fermi ground state of dipolar fermions in a circular storage ring

Recent developments in the field of ultracold gases has led to the production of degenerate samples of polar molecules. These have large static electric-dipole moments, which in turn causes the molecules to interact strongly. We investigate the interaction of polar particles in waveguide geometries subject to an applied polarizing field. For circular waveguides, tilting the direction of the polarizing field creates a periodic inhomogeneity of the interparticle interaction. We explore the consequences of geometry and interaction for stability of the ground state within the Thomas-Fermi model. Certain combinations of tilt angles and interaction strengths are found to preclude the existence of a stable Thomas-Fermi ground state. The system is shown to exhibit different behavior for quasi-one-dimensional and three-dimensional trapping geometries.

Physicochemical characterization of nanoparticles formed between DNA and phosphorylcholine substituted chitosans

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mass-imbalanced three-body systems in two dimensions

We consider three-body systems in two dimensions with zero-range interactions for general masses and interaction strengths. The momentum-space Schrödinger equation is solved numerically and in the Born-Oppenheimer (BO) approximation. The BO expression is derived using separable potentials and yields a concise adiabatic potential between the two heavy particles. The BO potential is Coulomb-like and exponentially decreasing at small and large distances, respectively. While we find similar qualitative features to previous studies, we find important quantitative differences. Our results demonstrate that mass-imbalanced systems that are accessible in the field of ultracold atomic gases can have a rich three-body bound state spectrum in two-dimensional geometries. Small light-heavy mass ratios increase the number of bound states. For 87Rb-87Rb-6Li and 133Cs- 133Cs-6Li we find respectively three and four bound states. © 2013 IOP Publishing Ltd.

Síntese, caracterização e estudo físico-químico de nanopartículas formadas pela interação de derivados hidrofílicos de quitosana com DNA

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Mass-Imbalanced Three-Body Systems in 2D: Bound States and the Analytical Approach to the Adiabatic Potential

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Multimodality and flexibility of stochastic gene expression

We consider a general class of mathematical models for stochastic gene expression where the transcription rate is allowed to depend on a promoter state variable that can take an arbitrary (finite) number of values. We provide the solution of the master equations in the stationary limit, based on a factorization of the stochastic transition matrix that separates timescales and relative interaction strengths, and we express its entries in terms of parameters that have a natural physical and/or biological interpretation. The solution illustrates the capacity of multiple states promoters to generate multimodal distributions of gene products, without the need for feedback. Furthermore, using the example of a three states promoter operating at low, high, and intermediate expression levels, we show that using multiple states operons will typically lead to a significant reduction of noise in the system. The underlying mechanism is that a three-states promoter can change its level of expression from low to high by passing through an intermediate state with a much smaller increase of fluctuations than by means of a direct transition.

The swelling behaviour of polyelectrolyte networks

Being basic ingredients of numerous daily-life products with significant industrial importance as well as basic building blocks for biomaterials, charged hydrogels continue to pose a series of unanswered challenges for scientists even after decades of practical applications and intensive research efforts. Despite a rather simple internal structure it is mainly the unique combination of short- and long-range forces which render scientific investigations of their characteristic properties to be quite difficult. Hence early on computer simulations were used to link analytical theory and empirical experiments, bridging the gap between the simplifying assumptions of the models and the complexity of real world measurements. Due to the immense numerical effort, even for high performance supercomputers, system sizes and time scales were rather restricted until recently, whereas it only now has become possible to also simulate a network of charged macromolecules. This is the topic of the presented thesis which investigates one of the fundamental and at the same time highly fascinating phenomenon of polymer research: The swelling behaviour of polyelectrolyte networks. For this an extensible simulation package for the research on soft matter systems, ESPResSo for short, was created which puts a particular emphasis on mesoscopic bead-spring-models of complex systems. Highly efficient algorithms and a consistent parallelization reduced the necessary computation time for solving equations of motion even in case of long-ranged electrostatics and large number of particles, allowing to tackle even expensive calculations and applications. Nevertheless, the program has a modular and simple structure, enabling a continuous process of adding new potentials, interactions, degrees of freedom, ensembles, and integrators, while staying easily accessible for newcomers due to a Tcl-script steering level controlling the C-implemented simulation core. Numerous analysis routines provide means to investigate system properties and observables on-the-fly. Even though analytical theories agreed on the modeling of networks in the past years, our numerical MD-simulations show that even in case of simple model systems fundamental theoretical assumptions no longer apply except for a small parameter regime, prohibiting correct predictions of observables. Applying a "microscopic" analysis of the isolated contributions of individual system components, one of the particular strengths of computer simulations, it was then possible to describe the behaviour of charged polymer networks at swelling equilibrium in good solvent and close to the Theta-point by introducing appropriate model modifications. This became possible by enhancing known simple scaling arguments with components deemed crucial in our detailed study, through which a generalized model could be constructed. Herewith an agreement of the final system volume of swollen polyelectrolyte gels with results of computer simulations could be shown successfully over the entire investigated range of parameters, for different network sizes, charge fractions, and interaction strengths. In addition, the "cell under tension" was presented as a self-regulating approach for predicting the amount of swelling based on the used system parameters only. Without the need for measured observables as input, minimizing the free energy alone already allows to determine the the equilibrium behaviour. In poor solvent the shape of the network chains changes considerably, as now their hydrophobicity counteracts the repulsion of like-wise charged monomers and pursues collapsing the polyelectrolytes. Depending on the chosen parameters a fragile balance emerges, giving rise to fascinating geometrical structures such as the so-called pear-necklaces. This behaviour, known from single chain polyelectrolytes under similar environmental conditions and also theoretically predicted, could be detected for the first time for networks as well. An analysis of the total structure factors confirmed first evidences for the existence of such structures found in experimental results.