Invaded cluster dynamics for frustrated models

The invaded cluster (IC) dynamics introduced by Machta et al. [Phys. Rev. Lett. 75, 2792 (1995)] is extended to the fully frustrated Ising model on a square lattice. The properties of the dynamics that exhibits numerical evidence of self-organized criticality are studied. The fluctuations in the IC dynamics are shown to be intrinsic of the algorithm and the fluctuation-dissipation theorem is no longer valid. The relaxation time is found to be very short and does not present a critical size dependence.

Damage spreading transition in glasses: A probe for the ruggedness of the configurational landscape

We consider damage spreading transitions in the framework of mode-coupling theory. This theory describes relaxation processes in glasses in the mean-field approximation which are known to be characterized by the presence of an exponentially large number of metastable states. For systems evolving under identical but arbitrarily correlated noises, we demonstrate that there exists a critical temperature T0 which separates two different dynamical regimes depending on whether damage spreads or not in the asymptotic long-time limit. This transition exists for generic noise correlations such that the zero damage solution is stable at high temperatures, being minimal for maximal noise correlations. Although this dynamical transition depends on the type of noise correlations, we show that the asymptotic damage has the good properties of a dynamical order parameter, such as (i) independence of the initial damage; (ii) independence of the class of initial condition; and (iii) stability of the transition in the presence of asymmetric interactions which violate detailed balance. For maximally correlated noises we suggest that damage spreading occurs due to the presence of a divergent number of saddle points (as well as metastable states) in the thermodynamic limit consequence of the ruggedness of the free-energy landscape which characterizes the glassy state. These results are then compared to extensive numerical simulations of a mean-field glass model (the Bernasconi model) with Monte Carlo heat-bath dynamics. The freedom of choosing arbitrary noise correlations for Langevin dynamics makes damage spreading an interesting tool to probe the ruggedness of the configurational landscape.

Magnetic relaxation of a one-dimensional model for small particle systems with dipolar interactions: Monte Carlo simulation

Atribution as a function of the time are analyzed and this study leads to a deeper knowledge of the microscopic processes involved in the magnetic relaxation

Dynamics of driven three-dimensional thin films: from hdrophilic to superhydrophobic substrates

We study the forced displacement of a thin film of fluid in contact with vertical and inclined substrates of different wetting properties, that range from hydrophilic to hydrophobic, using the lattice-Boltzmann method. We study the stability and pattern formation of the contact line in the hydrophilic and superhydrophobic regimes, which correspond to wedge-shaped and nose-shaped fronts, respectively. We find that contact lines are considerably more stable for hydrophilic substrates and small inclination angles. The qualitative behavior of the front in the linear regime remains independent of the wetting properties of the substrate as a single dispersion relation describes the stability of both wedges and noses. Nonlinear patterns show a clear dependence on wetting properties and substrate inclination angle. The effect is quantified in terms of the pattern growth rate, which vanishes for the sawtooth pattern and is finite for the finger pattern. Sawtooth shaped patterns are observed for hydrophilic substrates and low inclination angles, while finger-shaped patterns arise for hydrophobic substrates and large inclination angles. Finger dynamics show a transient in which neighboring fingers interact, followed by a steady state where each finger grows independently.

Dynamics of Fréedericksz transition in a fluctuating magnetic field

Fréedericksz transition under twist deformation in a nematic layer is discussed when the magnetic field has a random component. A dynamical model which includes the thermal fluctuations of the system is presented. The randomness of the field produces a shift of the instability point. Beyond this instability point the time constant characteristic of the approach to the stationary stable state decreases because of the field fluctuations. The opposite happens for fields smaller than the critical one. The decay time of an unstable state, calculated as a mean first-passage time, is also decreased by the field fluctuations.

Evaporation and coarsening dynamics with open boundaries

We present a study of the evaporation dynamics of a substance undergoing a coarsening process. The system is modeled by the Cahn-Hilliard equation with absorbing boundaries. We have found that the dynamics, although of a diffusive nature, is much slower than the usual one without coarsening. Analytical and simulation results are in reasonable agreement.

State-to-state reaction probabilities within the quantum transition state framework

Rigorous quantum dynamics calculations of reaction rates and initial state-selected reaction probabilities of polyatomic reactions can be efficiently performed within the quantum transition state concept employing flux correlation functions and wave packet propagation utilizing the multi-configurational time-dependent Hartree approach. Here, analytical formulas and a numerical scheme extending this approach to the calculation of state-to-state reaction probabilities are presented. The formulas derived facilitate the use of three different dividing surfaces: two dividing surfaces located in the product and reactant asymptotic region facilitate full state resolution while a third dividing surface placed in the transition state region can be used to define an additional flux operator. The eigenstates of the corresponding thermal flux operator then correspond to vibrational states of the activated complex. Transforming these states to reactant and product coordinates and propagating them into the respective asymptotic region, the full scattering matrix can be obtained. To illustrate the new approach, test calculations study the D + H2(ν, j) → HD(ν′, j′) + H reaction for J = 0.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

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.

Transient patterns in nematic liquid crystals: Domain-wall dynamics

We study the dynamics of the late stages of the Fréedericksz transition in which a periodic transient pattern decays to a final homogeneous state. A stability analysis of an unstable stationary pattern is presented, and equations for the evolution of the domain walls are obtained. Using results of previous theories, we analyze the effect that the specific dynamics of the problem, incorporating hydrodynamic couplings, has on the expected logarithmic domain growth law.

Antagonism of neuromuscular blockade but not muscle relaxation affects depth of anaesthesia

RÉSUMÉ Introduction L'effet des agents myorelaxants ainsi que des anticholinestérases sur la profondeur d'anesthésie a été étudié avec des résultats contradictoires. C'est pourquoi nous avons évalué l'effet de l'atracurium et de la néostigmine sur le BIS (bispectral index) ainsi que sur les potentiels auditives évoqués (middle-latency auditory evoked potentials, A-Line® autoregressive index [AAI]). Méthodes Après avoir obtenu l'accord du comité d'éthique local, nous avons étudié 40 patients ayant donné leur consentement écrit, ASA I-II, âgé de 18-69 ans. L'anesthésie générale a consisté en anesthésie intra-veineuse à objectif de concentration avec du propofol et du remifentanil. La fonction de la jonction neuromusculaire était monitorée en continu au moyen d'un électromyographe. Le BIS et l'AAI ont été enregistrés en continu. Après avoir atteint des valeurs stables au niveau du BIS, les patients ont été attribués à deux groupes par randomisation. Les patients du groupe 1 ont reçu 0.4 mg kg-1 d'atracurium et 5 minutes plus tard le même volume de NaCI 0.9%, dans le groupe 2 la séquence d'injection était inversée, le NaCI 0.9% en premier et l'atracurium en deuxième. Au moment où le premier « twitch » d'un train de quatre atteignait 10% de l'intensité avant la relaxation, les patients ont été randomisés une deuxième fois. Les patients du groupe N ont reçu 0.04 mg kg-1 de néostigmine et 0.01 rn9 kg-1 de glycopyrrolate alors que le groupe contrôle (G) ne recevait que 0.01 mg kg-] de glycopyrrolate. Résultats : L 'injection d'atracurium ou de NaCI 0.9% n'a pas eu d'effet sur le BIS ou l'AAI. Après l'injection de néostigmine avec glycopyrrolate, le BIS et I `AAI a augmenté de manière significative (changement maximal moyen du BIS 7.1 ± 7.5, P< 0.001, de l'AAI 9.7 ± 10.5, P< 0.001). Suite à l'injection de glycopyrrolate seule, le BIS et l'AAI a augmenté également (changement maximal moyen du BIS 2.2 ± 3.4, P< 0.008, de l'AAI 3.5 ± 5.7, P< 0.012), mais cette augmentation était significativement moins importante que dans le groupe N (P< 0.012 pour le BIS, P< 0.027 pour l'AAI). Conclusion Ces résultats laissent supposer que la néostigmine peut altérer la profondeur de l'anesthésie. La diminution de la profondeur d'anesthésie enregistrée par le BIS et l'AAI correspond probablement à une réapparition brusque d'une stimulation centrale liée à la proprioception. Au contraire, lors de la curarisation, le tonus musculaire diminue de manière beaucoup plus progressive, pouvant ainsi expliquer l'absence d'effet sur la profondeur d'anesthésie. ABSTRACT Background. Conflicting effects of neuromuscular blocking drugs and anticholinesterases on depth of anaesthesia have been reported. Therefore we evaluated the effect of atracurium and neostigmine on bispectral index (BIS) and middle-latency auditory evoked potentials (AAI). Methods. We studied 40 patients (ASA I-II) aged 18-69 yr. General anaesthesia consisted of propofol and remifentanil by target-controlled infusion and neuromuscular function was monitored by electromyography. When BIS reached stable values, patients were randomly assigned to one of two groups. Group I received atracurium 0.4 mg kg-1 and, 5 min later, the same volume of NaCl 0.9%; group 2 received saline first and then atracurium. When the first twitch of a train of four reached 10% of control intensity, patients were again randomized: one group (N) received neostigmine 0.04 mg kg-1 and glycopyrrolate 0.01 mg kg-1, and the control group (G) received only glycopyrrolate. Results. Injection of atracurium or NaCl 0.9% had no effect on BIS or AAI. After neostigmine¬glycopyrrolate, BIS and AAI increased significantly (mean maximal change of BIS 7.1 [SD 7.5], P<0.001; mean maximal change of AAI 9.7 [10.5], P<0.001). When glycopyrrolate was injected alone BIS and AAI also increased (mean maximal change of BIS 2.2 [3.4], P=0.008; mean maximal change of AAI 3.5 [5.7], P=0.012), but this increase was significantly less than in group N (P=0.012 for BIS; P=0.027 for AAI). Conclusion. These data suggest that neostigmine alters the state of propofol-remifentanil anaesthesia and may enhance recovery.

Comment on “Dynamics of thermal growth of silicon oxide films on Si”

The paper commented on here R. M. C. de Almeida, S. Gonçalves, I. J. R. Baumvol and F. C. Stedile Phys. Rev. B 61 12992 (2000) claims that the Deal and Grove model of oxidation is unable to describe the kinetics in the thin oxide regime due to two main simplifications: (a) the steady-state assumption and (b) the abrupt Si∕SiO2 interface assumption. Although reasonably good fits are obtained without these simplifications, it will be shown that the values of the kinetic parameters are not reliable and that the solutions given for different partial pressures are erroneous. Finally, it will be shown that the correct solution of their model is unable to predict the oxidation rate enhancement observed in the thin oxide regime and that the predicted width of the interface compatible with the Deal and Grove rate constants is too large

Development of Asphalt Dynamic Modulus Master Curve Using Falling Weight Deflectometer Measurements, TR-659

The asphalt concrete (AC) dynamic modulus (|E*|) is a key design parameter in mechanistic-based pavement design methodologies such as the American Association of State Highway and Transportation Officials (AASHTO) MEPDG/Pavement-ME Design. The objective of this feasibility study was to develop frameworks for predicting the AC |E*| master curve from falling weight deflectometer (FWD) deflection-time history data collected by the Iowa Department of Transportation (Iowa DOT). A neural networks (NN) methodology was developed based on a synthetically generated viscoelastic forward solutions database to predict AC relaxation modulus (E(t)) master curve coefficients from FWD deflection-time history data. According to the theory of viscoelasticity, if AC relaxation modulus, E(t), is known, |E*| can be calculated (and vice versa) through numerical inter-conversion procedures. Several case studies focusing on full-depth AC pavements were conducted to isolate potential backcalculation issues that are only related to the modulus master curve of the AC layer. For the proof-of-concept demonstration, a comprehensive full-depth AC analysis was carried out through 10,000 batch simulations using a viscoelastic forward analysis program. Anomalies were detected in the comprehensive raw synthetic database and were eliminated through imposition of certain constraints involving the sigmoid master curve coefficients. The surrogate forward modeling results showed that NNs are able to predict deflection-time histories from E(t) master curve coefficients and other layer properties very well. The NN inverse modeling results demonstrated the potential of NNs to backcalculate the E(t) master curve coefficients from single-drop FWD deflection-time history data, although the current prediction accuracies are not sufficient to recommend these models for practical implementation. Considering the complex nature of the problem investigated with many uncertainties involved, including the possible presence of dynamics during FWD testing (related to the presence and depth of stiff layer, inertial and wave propagation effects, etc.), the limitations of current FWD technology (integration errors, truncation issues, etc.), and the need for a rapid and simplified approach for routine implementation, future research recommendations have been provided making a strong case for an expanded research study.

Energy relaxation in nonlinear one-dimensional lattices

We study energy relaxation in thermalized one-dimensional nonlinear arrays of the Fermi-Pasta-Ulam type. The ends of the thermalized systems are placed in contact with a zero-temperature reservoir via damping forces. Harmonic arrays relax by sequential phonon decay into the cold reservoir, the lower-frequency modes relaxing first. The relaxation pathway for purely anharmonic arrays involves the degradation of higher-energy nonlinear modes into lower-energy ones. The lowest-energy modes are absorbed by the cold reservoir, but a small amount of energy is persistently left behind in the array in the form of almost stationary low-frequency localized modes. Arrays with interactions that contain both a harmonic and an anharmonic contribution exhibit behavior that involves the interplay of phonon modes and breather modes. At long times relaxation is extremely slow due to the spontaneous appearance and persistence of energetic high-frequency stationary breathers. Breather behavior is further ascertained by explicitly injecting a localized excitation into the thermalized arrays and observing the relaxation behavior.

From caging to rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory

By means of computer simulations and solution of the equations of the mode coupling theory (MCT),we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question,provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This dis-agreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.