Irreversible spherical model and its stationary entropy production rate

The nonequilibrium stationary state of an irreversible spherical model is investigated on hypercubic lattices. The model is defined by Langevin equations similar to the reversible case, but with asymmetric transition rates. In spite of being irreversible, we have succeeded in finding an explicit form for the stationary probability distribution, which turns out to be of the Boltzmann-Gibbs type. This enables one to evaluate the exact form of the entropy production rate at the stationary state, which is non-zero if the dynamical rules of the transition rates are asymmetric.

Exact field-driven interface dynamics in the two-dimensional stochastic Ising model with helicoidal boundary conditions

We investigate the interface dynamics of the two-dimensional stochastic Ising model in an external field under helicoidal boundary conditions. At sufficiently low temperatures and fields, the dynamics of the interface is described by an exactly solvable high-spin asymmetric quantum Hamiltonian that is the infinitesimal generator of the zero range process. Generally, the critical dynamics of the interface fluctuations is in the Kardar-Parisi-Zhang universality class of critical behavior. We remark that a whole family of RSOS interface models similar to the Ising interface model investigated here can be described by exactly solvable restricted high-spin quantum XXZ-type Hamiltonians. (C) 2012 Elsevier B.V. All rights reserved.

STRUCTURE AND DYNAMICS: THE TRANSITION FROM NONEQUILIBRIUM TO EQUILIBRIUM IN INTEGRATE-AND-FIRE DYNAMICS

The transient and equilibrium properties of dynamics unfolding in complex systems can depend critically on specific topological features of the underlying interconnections. In this work, we investigate such a relationship with respect to the integrate-and-fire dynamics emanating from a source node and an extended network model that allows control of the small-world feature as well as the length of the long-range connections. A systematic approach to investigate the local and global correlations between structural and dynamical features of the networks was adopted that involved extensive simulations (one and a half million cases) so as to obtain two-dimensional correlation maps. Smooth, but diverse surfaces of correlation values were obtained in all cases. Regarding the global cases, it has been verified that the onset avalanche time (but not its intensity) can be accurately predicted from the structural features within specific regions of the map (i.e. networks with specific structural properties). The analysis at local level revealed that the dynamical features before the avalanches can also be accurately predicted from structural features. This is not possible for the dynamical features after the avalanches take place. This is so because the overall topology of the network predominates over the local topology around the source at the stationary state.

Mean-Field critical behavior and ergodicity break in a nonequilibrium one-dimensional RSOS growth model

We investigate the nonequilibrium roughening transition of a one-dimensional restricted solid-on-solid model by directly sampling the stationary probability density of a suitable order parameter as the surface adsorption rate varies. The shapes of the probability density histograms suggest a typical Ginzburg-Landau scenario for the phase transition of the model, and estimates of the "magnetic" exponent seem to confirm its mean-field critical behavior. We also found that the flipping times between the metastable phases of the model scale exponentially with the system size, signaling the breaking of ergodicity in the thermodynamic limit. Incidentally, we discovered that a closely related model not considered before also displays a phase transition with the same critical behavior as the original model. Our results support the usefulness of off-critical histogram techniques in the investigation of nonequilibrium phase transitions. We also briefly discuss in the appendix a good and simple pseudo-random number generator used in our simulations.

Alfven eigenmode properties and dynamics in the TJ-II stellarator

Alfven eigenmodes (AEs) were studied in neutral beam injection (NBI) heated plasmas in the TJ-II stellarator using a heavy ion beam probe (HIBP) in the core, and by Langmuir (LP) and Mirnov probes (MP) at the edge. AEs were detected over the whole plasma radius by the HIBP with a spatial resolution of about 1 cm. AE-induced oscillations were detected in the plasma density n(e), electric potential phi and poloidal magnetic field B-pol with frequencies 50 kHz < f(AE) < 300 kHz. The LP, MP and HIBP data showed a high level of coherency for specific branches of AEs. Poloidal mode wave-vectors k(theta), mode numbers m (m < 8) and propagation velocities V-theta similar to 30 km s(-1) were detected for various branches of AEs, having different radial locations. When the density rose due to NBI fuelling, the AE frequency decreased as predicted by the Alfven law f(AE) similar to n(e)(-1/2). During the AE frequency decay the following new AE features were observed: (i) the poloidal wave-vector k(theta) and mode number m remained constant, (ii) the cross-phases between the oscillations in B-pol, n(e) and electric potential remained constant, having an individual value for each AE branch, (iii) V-theta decreased proportional to the AE frequency. The interaction of the AEs with the bulk (thermal) plasma resulted in clearly pronounced quasi-coherent peaks in the electrostatic turbulent particle flux spectra. Various AE branches exhibited different contributions to the particle flux: outward, inward and also zero, depending on the phase relations between the oscillations in E-pol and n(e), which are specific for each branch. A comparison with MHD mode modelling indicated that some of the more prominent frequency branches can be identified as radially extended helical AEs.

Analysis of energy dissipation in stirred suspension polymerisation reactors using computational fluid dynamics

This work evaluates the spatial distribution of normalised rates of droplet breakage and droplet coalescence in liquidliquid dispersions maintained in agitated tanks at operation conditions normally used to perform suspension polymerisation reactions. Particularly, simulations are performed with multiphase computational fluid dynamics (CFD) models to represent the flow field in liquidliquid styrene suspension polymerisation reactors for the first time. CFD tools are used first to compute the spatial distribution of the turbulent energy dissipation rates (e) inside the reaction vessel; afterwards, normalised rates of droplet breakage and particle coalescence are computed as functions of e. Surprisingly, multiphase simulations showed that the rates of energy dissipation can be very high near the free vortex surfaces, which has been completely neglected in previous works. The obtained results indicate the existence of extremely large energy dissipation gradients inside the vessel, so that particle breakage occurs primarily in very small regions that surround the impeller and the free vortex surface, while particle coalescence takes place in the liquid bulk. As a consequence, particle breakage should be regarded as an independent source term or a boundary phenomenon. Based on the obtained results, it can be very difficult to justify the use of isotropic assumptions to formulate particle population balances in similar systems, even when multiple compartment models are used to describe the fluid dynamic behaviour of the agitated vessel. (C) 2011 Canadian Society for Chemical Engineering

Quantum-Critical Spin Dynamics in Quasi-One-Dimensional Antiferromagnets

By means of nuclear spin-lattice relaxation rate T-1(-1), we follow the spin dynamics as a function of the applied magnetic field in two gapped quasi-one-dimensional quantum antiferromagnets: the anisotropic spin-chain system NiCl2-4SC(NH2)(2) and the spin-ladder system (C5H12N)(2)CuBr4. In both systems, spin excitations are confirmed to evolve from magnons in the gapped state to spinons in the gapless Tomonaga-Luttinger-liquid state. In between, T-1(-1) exhibits a pronounced, continuous variation, which is shown to scale in accordance with quantum criticality. We extract the critical exponent for T-1(-1), compare it to the theory, and show that this behavior is identical in both studied systems, thus demonstrating the universality of quantum-critical behavior.

Oocyte meiotic-stage-specific differences in spindle depolymerization in response to temperature changes monitored with polarized field microscopy and immunocytochemistry

Objective: To compare the polymerization status of mouse oocyte spindles exposed to various temperatures at various stages of meiosis. Design: Experimental animal study. Setting: University animal laboratory. Animal(s): CF1 mice. Intervention(s): Immature oocytes matured to metaphase I (MI), telophase I (TI), and metaphase II (MII) were incubated at 37 degrees C (control), room temperature (RT), or 4 degrees C for 0, 10, 30, and 60 minutes. Spindle analysis subsequently was performed using polarized field microscopy and immunocytochemistry. Spindles of TI and MII oocytes that underwent vitrification and warming were analyzed also by immunocytochemistry. Main Outcome Measure(s): Detection of polymerized meiotic spindles. Result(s): At RT, and after 60 minutes at 4 degrees C, a significant time-dependent decrease in the percentage of polymerized meiotic spindles was observed in MI and MII oocytes, but not in TI oocytes. The polymerization of TI spindles at 4 degrees C was similar to that of TI spindles at 4 degrees C that underwent vitrification and warming. Conclusion(s): Significant differences in the microtubule dynamics of MI, TI, and MII oocytes incubated at different temperatures were observed. In particular, meiotic spindles in TI oocytes exhibited less depolymerization than did metaphase spindles. (Fertil Steril (R) 2012; 97: 714-9. (C) 2012 by American Society for Reproductive Medicine.)

Dissipation effects in random transverse-field Ising chains

We study the effects of Ohmic, super-Ohmic, and sub-Ohmic dissipation on the zero-temperature quantum phase transition in the random transverse-field Ising chain by means of an (asymptotically exact) analytical strong-disorder renormalization-group approach. We find that Ohmic damping destabilizes the infinite-randomness critical point and the associated quantum Griffiths singularities of the dissipationless system. The quantum dynamics of large magnetic clusters freezes completely, which destroys the sharp phase transition by smearing. The effects of sub-Ohmic dissipation are similar and also lead to a smeared transition. In contrast, super-Ohmic damping is an irrelevant perturbation; the critical behavior is thus identical to that of the dissipationless system. We discuss the resulting phase diagrams, the behavior of various observables, and the implications to higher dimensions and experiments.

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Holographic field theory models of dark energy in interaction with dark matter

We discuss two Lagrangian interacting dark energy models in the context of the holographic principle. The potentials of the interacting fields are constructed. The models are compared with CMB distance information, baryonic acoustic oscillations, lookback time and the Constitution supernovae sample. For both models, the results are consistent with a nonvanishing interaction in the dark sector of the Universe and the sign of coupling is consistent with dark energy decaying into dark matter, alleviating the coincidence problem-with more than 3 standard deviations of confidence for one of them. However, this is because the noninteracting holographic dark energy model is a bad fit to the combination of data sets used in this work as compared to the cosmological constant with cold dark matter model, so that one needs to introduce the interaction in order to improve this model.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Statistical physics approach to dendritic computation: The excitable-wave mean-field approximation

We analytically study the input-output properties of a neuron whose active dendritic tree, modeled as a Cayley tree of excitable elements, is subjected to Poisson stimulus. Both single-site and two-site mean-field approximations incorrectly predict a nonequilibrium phase transition which is not allowed in the model. We propose an excitable-wave mean-field approximation which shows good agreement with previously published simulation results [Gollo et al., PLoS Comput. Biol. 5, e1000402 (2009)] and accounts for finite-size effects. We also discuss the relevance of our results to experiments in neuroscience, emphasizing the role of active dendrites in the enhancement of dynamic range and in gain control modulation.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

Effect of external fields in Axelrod's model of social dynamics

The study of the effects of spatially uniform fields on the steady-state properties of Axelrod's model has yielded plenty of counterintuitive results. Here, we reexamine the impact of this type of field for a selection of parameters such that the field-free steady state of the model is heterogeneous or multicultural. Analyses of both one- and two-dimensional versions of Axelrod's model indicate that the steady state remains heterogeneous regardless of the value of the field strength. Turning on the field leads to a discontinuous decrease on the number of cultural domains, which we argue is due to the instability of zero-field heterogeneous absorbing configurations. We find, however, that spatially nonuniform fields that implement a consensus rule among the neighborhood of the agents enforce homogenization. Although the overall effects of the fields are essentially the same irrespective of the dimensionality of the model, we argue that the dimensionality has a significant impact on the stability of the field-free homogeneous steady state.