Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution. (C) 2009 American Institute of Physics. [doi:10.1063/1.3253297]

Generalized Sznajd model for opinion propagation

In the last decade the Sznajd model has been successfully employed in modeling some properties and scale features of both proportional and majority elections. We propose a version of the Sznajd model with a generalized bounded confidence rule-a rule that limits the convincing capability of agents and that is essential to allow coexistence of opinions in the stationary state. With an appropriate choice of parameters it can be reduced to previous models. We solved this model both in a mean-field approach (for an arbitrary number of opinions) and numerically in a Barabaacutesi-Albert network (for three and four opinions), studying the transient and the possible stationary states. We built the phase portrait for the special cases of three and four opinions, defining the attractors and their basins of attraction. Through this analysis, we were able to understand and explain discrepancies between mean-field and simulation results obtained in previous works for the usual Sznajd model with bounded confidence and three opinions. Both the dynamical system approach and our generalized bounded confidence rule are quite general and we think it can be useful to the understanding of other similar models.

Role of noise in population dynamics cycles

Noise is an intrinsic feature of population dynamics and plays a crucial role in oscillations called phase-forgetting quasicycles by converting damped into sustained oscillations. This function of noise becomes evident when considering Langevin equations whose deterministic part yields only damped oscillations. We formulate here a consistent and systematic approach to population dynamics, leading to a Fokker-Planck equation and the associate Langevin equations in accordance with this conceptual framework, founded on stochastic lattice-gas models that describe spatially structured predator-prey systems. Langevin equations in the population densities and predator-prey pair density are derived in two stages. First, a birth-and-death stochastic process in the space of prey and predator numbers and predator-prey pair number is obtained by a contraction method that reduces the degrees of freedom. Second, a van Kampen expansion in the inverse of system size is then performed to get the Fokker-Planck equation. We also study the time correlation function, the asymptotic behavior of which is used to characterize the transition from the cyclic coexistence of species to the ordinary coexistence.

Power-law creep behavior of a semiflexible chain

Rheological properties of adherent cells are essential for their physiological functions, and microrheological measurements on living cells have shown that their viscoelastic responses follow a weak power law over a wide range of time scales. This power law is also influenced by mechanical prestress borne by the cytoskeleton, suggesting that cytoskeletal prestress determines the cell's viscoelasticity, but the biophysical origins of this behavior are largely unknown. We have recently developed a stochastic two-dimensional model of an elastically joined chain that links the power-law rheology to the prestress. Here we use a similar approach to study the creep response of a prestressed three-dimensional elastically jointed chain as a viscoelastic model of semiflexible polymers that comprise the prestressed cytoskeletal lattice. Using a Monte Carlo based algorithm, we show that numerical simulations of the chain's creep behavior closely correspond to the behavior observed experimentally in living cells. The power-law creep behavior results from a finite-speed propagation of free energy from the chain's end points toward the center of the chain in response to an externally applied stretching force. The property that links the power law to the prestress is the chain's stiffening with increasing prestress, which originates from entropic and enthalpic contributions. These results indicate that the essential features of cellular rheology can be explained by the viscoelastic behaviors of individual semiflexible polymers of the cytoskeleton.

Time correlation function in systems with two coexisting biological species

We study a stochastic lattice model describing the dynamics of coexistence of two interacting biological species. The model comprehends the local processes of birth, death, and diffusion of individuals of each species and is grounded on interaction of the predator-prey type. The species coexistence can be of two types: With self-sustained coupled time oscillations of population densities and without oscillations. We perform numerical simulations of the model on a square lattice and analyze the temporal behavior of each species by computing the time correlation functions as well as the spectral densities. This analysis provides an appropriate characterization of the different types of coexistence. It is also used to examine linked population cycles in nature and in experiment.

Frustration and anomalous behavior in the Bell-Lavis model of liquid water

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core,models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.

Elastic Maier-Saupe-Zwanzig model and some properties of nematic elastomers

We introduce a simple mean-field lattice model to describe the behavior of nematic elastomers. This model combines the Maier-Saupe-Zwanzig approach to liquid crystals and an extension to lattice systems of the Warner-Terentjev theory of elasticity, with the addition of quenched random fields. We use standard techniques of statistical mechanics to obtain analytic solutions for the full range of parameters. Among other results, we show the existence of a stress-strain coexistence curve below a freezing temperature, analogous to the P-V diagram of a simple fluid, with the disorder strength playing the role of temperature. Below a critical value of disorder, the tie lines in this diagram resemble the experimental stress-strain plateau and may be interpreted as signatures of the characteristic polydomain-monodomain transition. Also, in the monodomain case, we show that random fields may soften the first-order transition between nematic and isotropic phases, provided the samples are formed in the nematic state.

Importance of granular structure in the initial conditions for the elliptic flow

We show the effects of the granular structure of the initial conditions of a hydrodynamic description of high-energy nucleus-nucleus collisions on some observables, especially on the elliptic-flow parameter upsilon(2). Such a structure enhances production of isotropically distributed high-p(T) particles, making upsilon(2) smaller there. Also, it reduces upsilon(2) in the forward and backward regions where the global matter density is smaller and, therefore, where such effects become more efficacious.

Biaxial strain-induced suppression of spinodal decomposition in GaMnAs and GaCrAs

The thermodynamic properties of the magnetic semiconductors GaMnAs and GaCrAs are studied under biaxial strain. The calculations are based on the projector augmented wave method combined with the generalized quasichemical approach to treat the disorder and composition effects. Considering the influence of biaxial strain, we find a tendency to the suppression of binodal decomposition mainly for GaMnAs under compressive strain. For a substrate with a lattice constant 5% smaller than the one of GaAs, for GaMnAs, the solubility limit increases up to 40%. Thus, the strain can be a useful tool for tailoring magnetic semiconductors to the formation or not of embedded nanoclusters. (C) 2010 American Institute of Physics. [doi:10.1063/1.3448025]

Anelastic spectroscopy study of the metal-insulator transition of Nd(1-x)Eu(x)NiO(3)

Measurements are presented of the complex dynamic Young's modulus of NdNiO(3) and Nd(0.65)Eu(0.35)NiO(3) through the metal-insulator transition (MIT). Upon cooling, the modulus presents a narrow dip at the MIT followed by an abrupt stiffening of similar to 6%. The anomaly is reproducible between cooling and heating in Nd(0.65)Eu(0.35)NiO(3) but appears only as a slow stiffening during cooling in undoped NdNiO(3), in conformance with the fact that the MIT in RNiO(3) changes from strongly first order to second order when the mean R size is decreased. The elastic anomaly seems not to be associated with the antiferromagnetic transition, which is distinct from the MIT in Nd(0.65)Eu(0.35)NiO(3). It is concluded that the steplike stiffening is due to the disappearance or freezing of dynamic Jahn-Teller (JT) distortions through the MIT, where the JT active Ni(3+) is disproportionated into alternating Ni(3+delta) and Ni(3-delta). The fluctuating octahedral JT distortion necessary to justify the observed jump in the elastic modulus is estimated as similar to 3% but does not have a role in determining the MIT, since the otherwise-expected precursor softening is not observed.

Optical second harmonic generation in the centrosymmetric magnetic semiconductors EuTe and EuSe

The magnetic europium chalcogenide semiconductors EuTe and EuSe are investigated by the spectroscopy of second harmonic generation (SHG) in the vicinity of the optical band gap formed by transitions involving the 4f and 5d electronic orbitals of the magnetic Eu(2+) ions. In these materials with centrosymmetric crystal lattice the electric-dipole SHG process is symmetry forbidden so that no signal is observed in zero magnetic field. Signal appears, however, in applied magnetic field with the SHG intensity being proportional to the square of magnetization. The magnetic field and temperature dependencies of the induced SHG allow us to introduce a type of nonlinear optical susceptibility determined by the magnetic-dipole contribution in combination with a spontaneous or induced magnetization. The experimental results can be described qualitatively by a phenomenological model based on a symmetry analysis and are in good quantitative agreement with microscopic model calculations accounting for details of the electronic energy and spin structure.

Resistively detected NMR of the nu=1 quantum Hall state: A tilted magnetic field study

Previous resistively detected NMR (RDNMR) studies on the nu approximate to 1 quantum Hall state have reported a ""dispersionlike"" line shape and extremely short nuclear-spin-lattice relaxation times, observations which have been attributed to the formation of a skyrme lattice. Here we examine the evolution of the RDNMR line shape and nuclear-spin relaxation for Zeeman: Coulomb energy ratios ranging from 0.012 to 0.036. According to theory, suppression of the skyrme crystal, along with the associated Goldstone mode nuclear-spin-relaxation mechanism, is expected at the upper end of this range. However, we find that the anomalous line shape persists at high Zeeman energy, and only a modest decrease in the RDNMR-detected nuclear-spin-relaxation rate is observed.

Zero-phonon emission and magnetic polaron parameters in EuTe

A phonon structure in the photoluminescence of EuTe was discovered, with a well-defined zero-phonon emission line (ZPL). The ZPL redshifts linearly with the intensity of applied magnetic field, indicating spin relaxation of the photoexcited electron, and saturates at a lower magnetic field than the optical absorption bandgap, which is attributed to formation of magnetic polarons. From the difference in these saturation fields, the zero-field polaron binding energy and radius are estimated to be 43 meV and 3.2 (in units of the EuTe lattice parameter), respectively. (C) 2011 American Institute of Physics. [doi:10.1063/1.3634030]

EPR study of local symmetry sites of Ce(3+) in Pb(1-x)Ce(x)A (A=S, Se, and Te)

The local site symmetry of Ce(3+) ions in the diluted magnetic semiconductors Pb(1-x)Ce(x)A (A=S, Se, and Te) has been investigated by electron-paramagnetic resonance (EPR). The experiments were carried out on single crystals with cerium concentration x ranging from 0.001 to 0.035. The isotropic line due to Ce(3+) ions located at the substitutional Pb cation site with octahedral symmetry was observed for all the studied samples. We determined the effective Lande factors to be g=1.333, 1.364, and 1.402 for A=S, Se, and Te, respectively. The small difference with the predicted Lande factor g of 10/7 for the Gamma(7) (J=5/2) ground state was attributed to crystal-field admixture. In addition, EPR lines from Ce(3+) ions located at sites with small distortion from the original octahedral symmetry were also observed. Two distinct sites with axial distortion along the < 001 > crystallographic direction were identified and a third signal in the spectrum was attributed to sites with the cubic symmetry distorted along the < 110 > direction. The distortion at these distinct Ce sites is attributed to Pb lattice vacancies near the cerium ions that compensate for its donor activity.

Higher Moments of Net Proton Multiplicity Distributions at RHIC

200 GeV corresponding to baryon chemical potentials (mu(B)) between 200 and 20 MeV. Our measurements of the products kappa sigma(2) and S sigma, which can be related to theoretical calculations sensitive to baryon number susceptibilities and long-range correlations, are constant as functions of collision centrality. We compare these products with results from lattice QCD and various models without a critical point and study the root s(NN) dependence of kappa sigma(2). From the measurements at the three beam energies, we find no evidence for a critical point in the QCD phase diagram for mu(B) below 200 MeV.