A continuous strain-space model of viral evolution within a host

Viruses rapidly evolve, and HIV in particular is known to be one of the fastest evolving human viruses. It is now commonly accepted that viral evolution is the cause of the intriguing dynamics exhibited during HIV infections and the ultimate success of the virus in its struggle with the immune system. To study viral evolution, we use a simple mathematical model of the within-host dynamics of HIV which incorporates random mutations. In this model, we assume a continuous distribution of viral strains in a one-dimensional phenotype space where random mutations are modelled by di ffusion. Numerical simulations show that random mutations combined with competition result in evolution towards higher Darwinian fitness: a stable traveling wave of evolution, moving towards higher levels of fi tness, is formed in the phenoty space.

Reaction-diffusion wave fronts: Multigeneration biological species under climate change

A generalization of reaction-diffusion models to multigeneration biological species is presented. It is based on more complex random walks than those in previous approaches. The new model is developed analytically up to infinite order. Our predictions for the speed agree to experimental data for several butterfly species better than existing models. The predicted dependence for the speed on the number of generations per year allows us to explain the change in speed observed for a specific invasion

Effect of expiratory flow limitation on respiratory mechanical impedance: a model study

Large phasic variations of respiratory mechanical impedance (Zrs) have been observed during induced expiratory flow limitation (EFL) (M. Vassiliou, R. Peslin, C. Saunier, and C. Duvivier. Eur. Respir. J. 9: 779-786, 1996). To clarify the meaning of Zrs during EFL, we have measured from 5 to 30 Hz the input impedance (Zin) of mechanical analogues of the respiratory system, including flow-limiting elements (FLE) made of easily collapsible rubber tubing. The pressures upstream (Pus) and downstream (Pds) from the FLE were controlled and systematically varied. Maximal flow (Vmax) increased linearly with Pus, was close to the value predicted from wave-speed theory, and was obtained for Pus-Pds of 4-6 hPa. The real part of Zin started increasing abruptly with flow (V) >85%Vmax and either further increased or suddenly decreased in the vicinity of V¿max. The imaginary part of Zin decreased markedly and suddenly above 95%Vmax. Similar variations of Zin during EFL were seen with an analogue that mimicked the changes of airway transmural pressure during breathing. After pressure andV measurements upstream and downstream from the FLE were combined, the latter was analyzed in terms of a serial (Zs) and a shunt (Zp) compartment. Zs was consistent with a large resistance and inertance, and Zp with a mainly elastic element having an elastance close to that of the tube walls. We conclude that Zrs data during EFL mainly reflect the properties of the FLE.

Regular wave propagation out of noise in chemical active media

A pacemaker, regularly emitting chemical waves, is created out of noise when an excitable photosensitive Belousov-Zhabotinsky medium, strictly unable to autonomously initiate autowaves, is forced with a spatiotemporal patterned random illumination. These experimental observations are also reproduced numerically by using a set of reaction-diffusion equations for an activator-inhibitor model, and further analytically interpreted in terms of genuine coupling effects arising from parametric fluctuations. Within the same framework we also address situations of noise-sustained propagation in subexcitable media.

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.

s-wave charmed baryon resonances from a coupled-channel approach with heavy quark symmetry

We study charmed baryon resonances that are generated dynamically within a unitary meson-baryon coupled-channel model that treats the heavy pseudoscalar and vector mesons on equal footing as required by heavy-quark symmetry. It is an extension of recent SU(4) models with t-channel vector-meson exchanges to an SU(8) spin-flavor scheme, but differs considerably from the SU(4) approach in how the strong breaking of the flavor symmetry is implemented. Some of our dynamically generated states can be readily assigned to recently observed baryon resonances, while others do not have a straightforward identification and require the compilation of more data as well as an extension of the model to d-wave meson-baryon interactions and p-wave coupling in the neglected s- and u-channel diagrams. Of several novelties, we find that the Delta c(2595), which emerged as a ND quasibound state within the SU(4) approaches, becomes predominantly a ND* quasibound state in the present SU(8) scheme.

Expanding scroll rings and negative tension turbulence in a model of excitable media

Scroll waves in excitable media, described by the Barkley model, are studied. In the parameter region of weak excitability, negative tension of wave filaments is found. It leads to expansion of scroll rings and instability of wave filaments. A circular filament tends to stretch, bend, loop, and produce an expanding tangle that fills up the volume. The filament does not undergo fragmentation before it touches the boundaries. Statistical properties of such Winfree turbulence of scroll waves are numerically investigated.

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.

Measures of ionicity of alkaline-earth oxides from the analysis of ab initio cluster wave functions

We present an analysis of the M-O chemical bonding in the binary oxides MgO, CaO, SrO, BaO, and Al2O3 based on ab initio wave functions. The model used to represent the local environment of a metal cation in the bulk oxide is an MO6 cluster which also includes the effect of the lattice Madelung potential. The analysis of the wave functions for these clusters leads to the conclusion that all the alkaline-earth oxides must be regarded as highly ionic oxides; however, the ionic character of the oxides decreases as one goes from MgO, almost perfectly ionic, to BaO. In Al2O3 the ionic character is further reduced; however, even in this case, the departure from the ideal, fully ionic, model of Al3+ is not exceptionally large. These conclusions are based on three measures, a decomposition of the Mq+-Oq- interaction energy, the number of electrons associated to the oxygen ions as obtained from a projection operator technique, and the analysis of the cation core-level binding energies. The increasing covalent character along the series MgO, CaO, SrO, and BaO is discussed in view of the existing theoretical models and experimental data.

Ground and excited states of KNiF3: An ab initio cluster-model approach

Finite cluster models and a variety of ab initio wave functions have been used to study the electronic structure of bulk KNiF3. Several electronic states, including the ground state and some charge-transfer excited states, have been considered. The study of the cluster-model wave functions has permitted an understanding of the nature of the chemical bond in the electronic ground state. This is found to be highly ionic and the different ionic and covalent contributions to the bonding have been identified and quantified. Finally, we have studied the charge-transfer excited states leading to the optical gap and have found that calculated and experimental values are in good agreement. The wave functions corresponding to these excited states have also been analyzed and show that although KNiF3 may be described as a ligand-to-metal charge-transfer insulator there is a strong configuration mixing with the metal-to-metal charge-transfer states.

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.

Regular wave propagation out of noise in chemical active media

A pacemaker, regularly emitting chemical waves, is created out of noise when an excitable photosensitive Belousov-Zhabotinsky medium, strictly unable to autonomously initiate autowaves, is forced with a spatiotemporal patterned random illumination. These experimental observations are also reproduced numerically by using a set of reaction-diffusion equations for an activator-inhibitor model, and further analytically interpreted in terms of genuine coupling effects arising from parametric fluctuations. Within the same framework we also address situations of noise-sustained propagation in subexcitable media.

Wave-height hazard analysis in Eastern Coast of Spain : Bayesian approach using generalized Pareto distribution

Standard practice of wave-height hazard analysis often pays little attention to the uncertainty of assessed return periods and occurrence probabilities. This fact favors the opinion that, when large events happen, the hazard assessment should change accordingly. However, uncertainty of the hazard estimates is normally able to hide the effect of those large events. This is illustrated using data from the Mediterranean coast of Spain, where the last years have been extremely disastrous. Thus, it is possible to compare the hazard assessment based on data previous to those years with the analysis including them. With our approach, no significant change is detected when the statistical uncertainty is taken into account. The hazard analysis is carried out with a standard model. Time-occurrence of events is assumed Poisson distributed. The wave-height of each event is modelled as a random variable which upper tail follows a Generalized Pareto Distribution (GPD). Moreover, wave-heights are assumed independent from event to event and also independent of their occurrence in time. A threshold for excesses is assessed empirically. The other three parameters (Poisson rate, shape and scale parameters of GPD) are jointly estimated using Bayes' theorem. Prior distribution accounts for physical features of ocean waves in the Mediterranean sea and experience with these phenomena. Posterior distribution of the parameters allows to obtain posterior distributions of other derived parameters like occurrence probabilities and return periods. Predictives are also available. Computations are carried out using the program BGPE v2.0