The condensation and ordering of models of empty liquids

We consider a simple model consisting of particles with four bonding sites ("patches"), two of type A and two of type B, on the square lattice, and investigate its global phase behavior by simulations and theory. We set the interaction between B patches to zero and calculate the phase diagram as the ratio between the AB and the AA interactions, epsilon(AB)*, varies. In line with previous work, on three-dimensional off-lattice models, we show that the liquid-vapor phase diagram exhibits a re-entrant or "pinched" shape for the same range of epsilon(AB)*, suggesting that the ratio of the energy scales - and the corresponding empty fluid regime - is independent of the dimensionality of the system and of the lattice structure. In addition, the model exhibits an order-disorder transition that is ferromagnetic in the re-entrant regime. The use of low-dimensional lattice models allows the simulation of sufficiently large systems to establish the nature of the liquid-vapor critical points and to describe the structure of the liquid phase in the empty fluid regime, where the size of the "voids" increases as the temperature decreases. We have found that the liquid-vapor critical point is in the 2D Ising universality class, with a scaling region that decreases rapidly as the temperature decreases. The results of simulations and theoretical analysis suggest that the line of order-disorder transitions intersects the condensation line at a multi-critical point at zero temperature and density, for patchy particle models with a re-entrant, empty fluid, regime. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3657406]

Biaxial Nematic Order in the Hard-boomerang Fluid

We consider a fluid of hard boomerangs, each composed of two hard spherocylinders joined at their ends at an angle Psi. The resulting particle is nonconvex and biaxial. The occurence of nematic order in such a system has been investigated using Straley's theory, which is a simplificaton of Onsager's second-virial treatment of long hard rods, and by bifurcation analysis. The excluded volume of two hard boomerangs has been approximated by the sum of excluded volumes of pairs of constituent spherocylinders, and the angle-dependent second-virial coefficient has been replaced by a low-order interpolating function. At the so-called Landau point, Psi(Landau)approximate to 107.4 degrees, the fluid undergoes a continuous transition from the isotropic to a biaxial nematic (B) phase. For Psi not equal Psi(Landau) ordering is via a first-order transition into a rod-like uniaxial nematic phase (N(+)) if Psi > Psi(Landau), or a plate-like uniaxial nematic (N(-)) phase if Psi < Psi(Landau). The B phase is separated from the N(+) and N(-) phases by two lines of continuous transitions meeting at the Landau point. This topology of the phase diagram is in agreement with previous studies of spheroplatelets and biaxial ellipsoids. We have checked the accuracy of our theory by performing numerical calculations of the angle-dependent second virial coefficient, which yields Psi(Landau)approximate to 110 degrees for very long rods, and Psi(Landau)approximate to 90 degrees for short rods. In the latter case, the I-N transitions occur at unphysically high packing fractions, reflecting the inappropriateness of the second-virial approximation in this limit.

Reply to "Comment on 'Effect of polydispersity on the ordering transition of adsorbed self- assembled rigid rods'"

We comment on the nature of the ordering transition of a model of equilibrium polydisperse rigid rods on the square lattice, which is reported by Lopez et al. to exhibit random percolation criticality in the canonical ensemble, in sharp contrast to (i) our results of Ising criticality for the same model in the grand canonical ensemble [Phys. Rev. E 82, 061117 (2010)] and (ii) the absence of exponent(s) renormalization for constrained systems with logarithmic specific-heat anomalies predicted on very general grounds by Fisher [Phys. Rev. 176, 257 (1968)].

Chaos and crises in a model for cooperative hunting: A symbolic dynamics approach

In this work we investigate the population dynamics of cooperative hunting extending the McCann and Yodzis model for a three-species food chain system with a predator, a prey, and a resource species. The new model considers that a given fraction sigma of predators cooperates in prey's hunting, while the rest of the population 1-sigma hunts without cooperation. We use the theory of symbolic dynamics to study the topological entropy and the parameter space ordering of the kneading sequences associated with one-dimensional maps that reproduce significant aspects of the dynamics of the species under several degrees of cooperative hunting. Our model also allows us to investigate the so-called deterministic extinction via chaotic crisis and transient chaos in the framework of cooperative hunting. The symbolic sequences allow us to identify a critical boundary in the parameter spaces (K, C-0) and (K, sigma) which separates two scenarios: (i) all-species coexistence and (ii) predator's extinction via chaotic crisis. We show that the crisis value of the carrying capacity K-c decreases at increasing sigma, indicating that predator's populations with high degree of cooperative hunting are more sensitive to the chaotic crises. We also show that the control method of Dhamala and Lai [Phys. Rev. E 59, 1646 (1999)] can sustain the chaotic behavior after the crisis for systems with cooperative hunting. We finally analyze and quantify the inner structure of the target regions obtained with this control method for wider parameter values beyond the crisis, showing a power law dependence of the extinction transients on such critical parameters.

Measuring and Controlling the Chaotic Motion of Profits

The study of economic systems has generated deep interest in exploring the complexity of chaotic motions in economy. Due to important developments in nonlinear dynamics, the last two decades have witnessed strong revival of interest in nonlinear endogenous business chaotic models. The inability to predict the behavior of dynamical systems in the presence of chaos suggests the application of chaos control methods, when we are more interested in obtaining regular behavior. In the present article, we study a specific economic model from the literature. More precisely, a system of three ordinary differential equations gather the variables of profits, reinvestments and financial flow of borrowings in the structure of a firm. Firstly, using results of symbolic dynamics, we characterize the topological entropy and the parameter space ordering of kneading sequences, associated with one-dimensional maps that reproduce significant aspects of the model dynamics. The analysis of the variation of this numerical invariant, in some realistic system parameter region, allows us to quantify and to distinguish different chaotic regimes. Finally, we show that complicated behavior arising from the chaotic firm model can be controlled without changing its original properties and the dynamics can be turned into the desired attracting time periodic motion (a stable steady state or into a regular cycle). The orbit stabilization is illustrated by the application of a feedback control technique initially developed by Romeiras et al. [1992]. This work provides another illustration of how our understanding of economic models can be enhanced by the theoretical and numerical investigation of nonlinear dynamical systems modeled by ordinary differential equations.

Torre solar: estudo do rendimento e viabilidade de centrais térmicas solares para produção de energia eléctrica de pequena e média dimensão

O aproveitamento da radiação solar representa um recurso energético extremamente benéfico, quer no âmbito energético, quer no ambiental, contribuindo para a redução das emissões de gases nocivos para a atmosfera. Portugal apresenta uma radiação solar total média anual bastante elevada, colocando-se entre um dos países que apresentam melhores potencialidades para o aproveitamento da energia solar. A torre solar consiste numa estufa com uma determinada extensão com uma chaminé no seu centro e o seu funcionamento baseia-se no aquecimento do ar que circula por baixo da estufa, sendo expelido pela chaminé. Nesta tese é abordado o tema da torre solar e os princípios físicos inerentes ao seu funcionamento. Foi estudado e descrito o método de cálculo de diversos parâmetros e resultados associados ao funcionamento da torre solar. Elaborou-se uma folha de cálculo para obtenção dos valores de simulações de torres com diversas dimensões, tecendo-se conclusões quanto aos resultados e às variações dos mesmos, consoante as alterações de dimensão dos elementos que a constituem. Foram descritos os vários elementos que constituem uma torre solar, bem como as suas características e tipologias. Efectuou-se um estudo com aplicação a um caso real, para se tecer algumas conclusões e comentários relativamente à viabilidade de uma torre solar para a situação em causa. Por fim, com base em todo o trabalho desenvolvido e abordado, foi possível tecer-se algumas conclusões quanto à viabilidade das torres solares.

Structural characterization and immunogenicity in wildtype and immune tolerant mice of degraded recombinant human interferon alpha2b

Purpose - To study the influence of protein structure on the immunogenicity in wildtype and immune tolerant mice of well-characterized degradation products of recombinant human interferon alpha2b (rhIFNα2b). Methods - RhIFNα2b was degraded by metal catalyzed oxidation (M), crosslinking with glutaraldehyde (G), oxidation with hydrogen peroxide (H) and incubation in a boiling water bath (B). The products were characterized with UV absorption, circular dichroism and fluorescence spectroscopy, gel permeation chromatography, reversed-phase HPLC, SDS-PAGE, Western blotting and mass spectrometry. The immunogenicity of the products was evaluated in wildtype mice and in transgenic mice immune tolerant for hIFNα2. Serum antibodies were detected by ELISA or surface plasmon resonance. Results - M-rhIFNα2b contained covalently aggregated rhIFNα2b with three methionines partly oxidized to methionine sulfoxides. G-rhIFNα2b contained covalent aggregates and did not show changes in secondary structure. H-rhIFNα2b was only chemically changed with four partly oxidized methionines. B-rhIFNα2b was largely unfolded and heavily aggregated. Native (N) rhIFNα2b was immunogenic in the wildtype mice but not in the transgenic mice, showing that the latter were immune tolerant for rhIFNα2b. The antirhIFNα2b antibody levels in the wildtype mice depended on the degradation product: M-rhIFNα2b > H-rhIFNα2b ~ N-rhIFNα2b >> B-rhIFNα2b; G-rhIFNα2b did not induce anti-rhIFNα2b antibodies. In the transgenic mice, only M-rhIFNα2b could break the immune tolerance. Conclusions - RhIFNα2b immunogenicity is related to its structural integrity. Moreover, the immunogenicity of aggregated rhIFNα2b depends on the structure and orientation of the constituent protein molecules and/or on the aggregate size.

Structural characterization and immunogenicity in wild-type and immune tolerant mice of degraded recombinant human interferon Alpha2b

Purpose: This study was conducted to study the influence of protein structure on the immunogenicity in wild-type and immune tolerant mice of well-characterized degradation products of recombinant human interferon alpha2b (rhIFNα2b). Methods: RhIFNα2b was degraded by metal-catalyzed oxidation (M), cross-linking with glutaraldehyde (G), oxidation with hydrogen peroxide (H), and incubation in a boiling water bath (B). The products were characterized with UV absorption, circular dichroism and fluorescence spectroscopy, gel permeation chromatography, reverse-phase high-pressure liquid chromatography, sodium dodecyl sulfate polyacrylamide gel electrophoresis, Western blotting, and mass spectrometry. The immunogenicity of the products was evaluated in wild-type mice and in transgenic mice immune tolerant for hIFNα2. Serum antibodies were detected by enzyme-linked immunosorbent assay or surface plasmon resonance. Results: M-rhIFNα2b contained covalently aggregated rhIFNα2b with three methionines partly oxidized to methionine sulfoxides. G-rhIFNα2b contained covalent aggregates and did not show changes in secondary structure. H-rhIFNα2b was only chemically changed with four partly oxidized methionines. B-rhIFNα2b was largely unfolded and heavily aggregated. Nontreated (N) rhIFNα2b was immunogenic in the wild-type mice but not in the transgenic mice, showing that the latter were immune tolerant for rhIFNα2b. The anti-rhIFNα2b antibody levels in the wild-type mice depended on the degradation product: M-rhIFNα2b > H-rhIFNα2b ∼ N-rhIFNα2b ≫ B-rhIFNα2b; G-rhIFNα2b did not induce anti-rhIFNα2b antibodies. In the transgenic mice, only M-rhIFNα2b could break the immune tolerance. Conclusions: RhIFNα2b immunogenicity is related to its structural integrity. Moreover, the immunogenicity of aggregated rhIFNα2b depends on the structure and orientation of the constituent protein molecules and/or on the aggregate size.

Analysis of functionally graded sandwich plate structures with piezoelectric skins, using B-spline finite strip method

Functionally graded materials are composite materials wherein the composition of the constituent phases can vary in a smooth continuous way with a gradation which is function of its spatial coordinates. This characteristic proves to be an important issue as it can minimize abrupt variations of the material properties which are usually responsible for localized high values of stresses, and simultaneously providing an effective thermal barrier in specific applications. In the present work, it is studied the static and free vibration behaviour of functionally graded sandwich plate type structures, using B-spline finite strip element models based on different shear deformation theories. The effective properties of functionally graded materials are estimated according to Mori-Tanaka homogenization scheme. These sandwich structures can also consider the existence of outer skins of piezoelectric materials, thus achieving them adaptive characteristics. The performance of the models, are illustrated through a set of test cases. (C) 2012 Elsevier Ltd. All rights reserved.

Topological Complexity and Predictability in the Dynamics of a Tumor Growth Model with Shilnikov's Chaos

Dynamical systems modeling tumor growth have been investigated to determine the dynamics between tumor and healthy cells. Recent theoretical investigations indicate that these interactions may lead to different dynamical outcomes, in particular to homoclinic chaos. In the present study, we analyze both topological and dynamical properties of a recently characterized chaotic attractor governing the dynamics of tumor cells interacting with healthy tissue cells and effector cells of the immune system. By using the theory of symbolic dynamics, we first characterize the topological entropy and the parameter space ordering of kneading sequences from one-dimensional iterated maps identified in the dynamics, focusing on the effects of inactivation interactions between both effector and tumor cells. The previous analyses are complemented with the computation of the spectrum of Lyapunov exponents, the fractal dimension and the predictability of the chaotic attractors. Our results show that the inactivation rate of effector cells by the tumor cells has an important effect on the dynamics of the system. The increase of effector cells inactivation involves an inverse Feigenbaum (i.e. period-halving bifurcation) scenario, which results in the stabilization of the dynamics and in an increase of dynamics predictability. Our analyses also reveal that, at low inactivation rates of effector cells, tumor cells undergo strong, chaotic fluctuations, with the dynamics being highly unpredictable. Our findings are discussed in the context of tumor cells potential viability.

Water-Based Cellulose Liquid Crystal System Investigated by Rheo-NMR

Water-based cellulose cholesteric liquid crystalline phases at rest can undergo structural changes induced by shear flow. This reflects on the deuterium spectra recorded when the system is investigated by rheo-nuclear magnetic resonance (rheo-NMR) techniques. In this work, the model system hydroxypropylcellulose (HPC)+water is revisited using rheo-NMR to clarify unsettled points regarding its behavior under shear and in relaxation. The NMR spectra allow the identification of five different stable ordering states, within shear and relaxation, which are well integrated in a mesoscopic picture of the system's structural evolution under shear and relaxation. This picture emerging from the large body of studies available for this system by other experimental techniques, accounts well for the NMR data and is in good agreement with the three distinct regions of steady shear flow recognized for some lyotropic LC polymers. Shear rates in between 0.1 and 1.0 s(-1) where investigated using a Taylor-Couette flow and deuterated water was used as solvent for the deuterium NMR (DNMR) analysis.

Ordenamento do território na cidade de Lisboa: o caso da Alta de Lisboa

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Edificações

Does independent component analysis play a role in unmixing hyperspectral data?

Independent component analysis (ICA) has recently been proposed as a tool to unmix hyperspectral data. ICA is founded on two assumptions: 1) the observed spectrum vector is a linear mixture of the constituent spectra (endmember spectra) weighted by the correspondent abundance fractions (sources); 2)sources are statistically independent. Independent factor analysis (IFA) extends ICA to linear mixtures of independent sources immersed in noise. Concerning hyperspectral data, the first assumption is valid whenever the multiple scattering among the distinct constituent substances (endmembers) is negligible, and the surface is partitioned according to the fractional abundances. The second assumption, however, is violated, since the sum of abundance fractions associated to each pixel is constant due to physical constraints in the data acquisition process. Thus, sources cannot be statistically independent, this compromising the performance of ICA/IFA algorithms in hyperspectral unmixing. This paper studies the impact of hyperspectral source statistical dependence on ICA and IFA performances. We conclude that the accuracy of these methods tends to improve with the increase of the signature variability, of the number of endmembers, and of the signal-to-noise ratio. In any case, there are always endmembers incorrectly unmixed. We arrive to this conclusion by minimizing the mutual information of simulated and real hyperspectral mixtures. The computation of mutual information is based on fitting mixtures of Gaussians to the observed data. A method to sort ICA and IFA estimates in terms of the likelihood of being correctly unmixed is proposed.

Does Independent Component Analysis Play a Role in Unmixing Hyperspectral Data?

Chapter in Book Proceedings with Peer Review First Iberian Conference, IbPRIA 2003, Puerto de Andratx, Mallorca, Spain, JUne 4-6, 2003. Proceedings

Parallel implementation of vertex component analysis for hyperspectral endmember extraction

International Conference with Peer Review 2012 IEEE International Conference in Geoscience and Remote Sensing Symposium (IGARSS), 22-27 July 2012, Munich, Germany