Optimal Linear Glauber Model

Contrary to the actual nonlinear Glauber model, the linear Glauber model (LGM) is exactly solvable, although the detailed balance condition is not generally satisfied. This motivates us to address the issue of writing the transition rate () in a best possible linear form such that the mean squared error in satisfying the detailed balance condition is least. The advantage of this work is that, by studying the LGM analytically, we will be able to anticipate how the kinetic properties of an arbitrary Ising system depend on the temperature and the coupling constants. The analytical expressions for the optimal values of the parameters involved in the linear are obtained using a simple Moore-Penrose pseudoinverse matrix. This approach is quite general, in principle applicable to any system and can reproduce the exact results for one dimensional Ising system. In the continuum limit, we get a linear time-dependent Ginzburg-Landau equation from the Glauber's microscopic model of non-conservative dynamics. We analyze the critical and dynamic properties of the model, and show that most of the important results obtained in different studies can be reproduced by our new mathematical approach. We will also show in this paper that the effect of magnetic field can easily be studied within our approach; in particular, we show that the inverse of relaxation time changes quadratically with (weak) magnetic field and that the fluctuation-dissipation theorem is valid for our model.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

Enhanced electroluminescence intensity of InGaN/GaN multi-quantum-wells based on Mg-doped GaN annealed in O-2

InGaN/GaN multi-quantum-well blue (461 +/- 4 nm) light emitting diodes with higher electroluminescence intensity are obtained by postgrowth thermal annealing at 720 C in O-2-ambient. Based on our first-principle total-energy calculations, we conclude that besides dissociating the Mg-H complex by forming H2O, annealing in O-2 has another positive effect on the activation of acceptor Mg in GaN. Mg can be further activated by the formation of an impurity band above the valence band maximum of host GaN from the passivated Mg-Ga-O-N complex. Our calculated ionization energy for acceptor Mg in the passivated system is about 30 meV shallower than that in pure GaN, in good agreement with previous experimental measurement. Our model can explain that the enhanced electroluminescence intensity of InGaN/GaN MQWs based on Mg-doped p-type GaN is due to a decrease in the ionization energy of Mg acceptor with the presence of oxygen. (C) 2008 American Institute of Physics.

STRUCTURAL EVOLUTION AND COMPOSITION CHANGE IN THE SURFACE REGION OF POLYPROPYLENE/CLAY NANOCOMPOSITES ANNEALED AT HIGH TEMPERATURES

A model experiment was done to clear the formation mechanism of protective layers during combustion of polypropylene (PP)/organically modified montmorillonite (OMMT) nanocomposites. The investigation was focused on the effects of annealing temperature on the structural changes and protective layer formation. The decomposition of OMMT and degradation of PP/OMMT nanocomposites were characterized by means of thermogravimetric analysis (TGA). The structural evolution and composition change in the surface region of PP/OMMT nanocomposites during heating were monitored by means of X-ray photoelectron spectroscopy (XPS), ATR-FTIR and field emission scanning electron microscopy (FESEM).

Multipartite nonlocality in a thermalized Ising spin chain

We study multipartite correlations and nonlocality in an isotropic Ising ring under transverse magnetic field at both zero and finite temperature. We highlight parity-induced differences between the multipartite Bell-like functions used in order to quantify the degree of nonlocality within a ring state and reveal a mechanism for the passive protection of multipartite quantum correlations against thermal spoiling effects that is clearly related to the macroscopic properties of the ring model.

Increasing entanglement through engineered disorder in the random Ising chain

The ground-state entanglement entropy between block of sites in the random Ising chain is studied by means of the Von Neumann entropy. We show that in presence of strong correlations between the disordered couplings and local magnetic fields the entanglement increases and becomes larger than in the ordered case. The different behavior with respect to the uncorrelated disordered model is due to the drastic change of the ground state properties. The same result holds also for the random three-state quantum Potts model.

Effect of thermal annealing on Fe40Ni38B18Mo4 thin films: modified Herzer model for magnetic evolution

Magnetic properties of nano-crystalline soft magnetic alloys have usually been correlated to structural evolution with heat treatment. However, literature reports pertaining to the study of nano-crystalline thin films are less abundant. Thin films of Fe40Ni38B18Mo4 were deposited on glass substrates under a high vacuum of ≈ 10−6 Torr by employing resistive heating. They were annealed at various temperatures ranging from 373 to 773K based on differential scanning calorimetric studies carried out on the ribbons. The magnetic characteristics were investigated using vibrating sample magnetometry. Morphological characterizations were carried out using atomic force microscopy (AFM), and magnetic force microscopy (MFM) imaging is used to study the domain characteristics. The variation of magnetic properties with thermal annealing is also investigated. From AFM and MFM images it can be inferred that the crystallization temperature of the as-prepared films are lower than their bulk counterparts. Also there is a progressive evolution of coercivity up to 573 K, which is an indication of the lowering of nano-crystallization temperature in thin films. The variation of coercivity with the structural evolution of the thin films with annealing is discussed and a plausible explanation is provided using the modified random anisotropy model

Compressible Sherrington-Kirkpatrick spin-glass model

We introduce a Sherrington-Kirkpatrick spin-glass model with the addition of elastic degrees of freedom. The problem is formulated in terms of an effective four-spin Hamiltonian in the pressure ensemble, which can be treated by the replica method. In the replica-symmetric approximation, we analyze the pressure-temperature phase diagram, and obtain expressions for the critical boundaries between the disordered and the ordered (spin-glass and ferromagnetic) phases. The second-order para-ferromagnetic border ends at a tricritical point, beyond which the transition becomes discontinuous. We use these results to make contact with the temperature-concentration phase diagrams of mixtures of hydrogen-bonded crystals.

Antiferromagnetic spherical spin-glass model

We study the thermodynamic properties and the phase diagrams of a multi-spin antiferromagnetic spherical spin-glass model using the replica method. It is a two-sublattice version of the ferromagnetic spherical p-spin glass model. We consider both the replica-symmetric and the one-step replica-symmetry-breaking solutions, the latter being the most general solution for this model. We find paramagnetic, spin-glass, antiferromagnetic and mixed or glassy antiferromagnetic phases. The phase transitions are always of second order in the thermodynamic sense, but the spin-glass order parameter may undergo a discontinuous change.

The energy landscape for solvent dynamics in electron transfer reactions: A minimalist model

Energy fluctuations of a solute molecule embedded in a polar solvent are investigated to depict the energy landscape for solvation dynamics. The system is modeled by a charged molecule surrounded by two layers of solvent dipolar molecules with simple rotational dynamics. Individual solvent molecules are treated as simple dipoles that can point toward or away from the central charge (Ising spins). Single-spin-flip Monte Carlo kinetics simulations are carried out in a two-dimensional lattice for different central charges, radii of outer shell, and temperatures. By analyzing the density of states as a function of energy and temperatures, we have determined the existence of multiple freezing transitions. Each of them can be associated with the freezing of a different layer of the solvent. (C) 2002 American Institute of Physics.

Agent-based model to rural-urban migration analysis

In this paper, we analyze the rural-urban migration phenomenon as it is usually observed in economies which are in the early stages of industrialization. The analysis is conducted by means of a statistical mechanics approach which builds a computational agent-based model. Agents are placed on a lattice and the connections among them are described via an Ising-like model. Simulations on this computational model show some emergent properties that are common in developing economies, such as a transitional dynamics characterized by continuous growth of urban population, followed by the equalization of expected wages between rural and urban sectors (Harris-Todaro equilibrium condition), urban concentration and increasing of per capita income. (c) 2005 Elsevier B.V. All rights reserved.

A Monte Carlo study of the anisotropic N=3 Ashkin-Teller model

The phase diagram of an asymmetric N = 3 Ashkin-Teller model is obtained by a numerical analysis which combines Monte Carlo renormalization group and reweighting techniques. Present results reveal several differences with those obtained by mean-field calculations and a Hamiltonian approach. In particular, we found Ising critical exponents along a line where Goldschmidt has located the Kosterlitz-Thouless multicritical point. On the other hand, we did find nonuniversal exponents along another transition line. Symmetry breaking in this case is very similar to the N = 2 case, since the symmetries associated to only two of the Ising variables are broken. However, for large values of the coupling constant ratio XW = W/K, when the only broken symmetry is of a hidden variable, we detected first-order phase transitions giving evidence supporting the existence of a multicritical point, as suggested by Goldschmidt, but in a different region of the phase diagram. © 2002 Elsevier Science B.V. All rights reserved.

Stability of Order in Solvent-Annealed Block Copolymer Thin Films

ABSTRACT: One way to produce high order in a block copolymer thin film is by solution casting a thin film and slowly evaporating the solvent in a sealed vessel. Such a solvent-annealing process is a versatile method to produce a highly ordered thin film of a block copolymer. However, the ordered structure of the film degrades over time when stored under ambient conditions. Remarkably, this aging process occurs in mesoscale thin films of polystyrene-polyisoprene triblock copolymer where the monolayer of vitrified 15 nm diameter polystyrene cylinders sink in a 20 nm thick film at 22 °C. The transformation is studied by atomic force microscopy (AFM). We describe the phenomena, characterize the aging process, and propose a semiquantitative model to explain the observations. The residual solvent effects are important but not the primary driving force for the aging process. The study may lead to effective avenue to improve order and make the morphology robust and possibly the solvent-annealing process more effective.

A macroscopic kinetic model for DNA polymerase elongation and high-fidelity nucleotide election

The enzymatically catalyzed template-directed extension of ssDNA/primer complex is an impor-tant reaction of extraordinary complexity. The DNA polymerase does not merely facilitate the insertion of dNMP, but it also performs rapid screening of substrates to ensure a high degree of fidelity. Several kinetic studies have determined rate constants and equilibrium constants for the elementary steps that make up the overall pathway. The information is used to develop a macro-scopic kinetic model, using an approach described by Ninio [Ninio J., 1987. Alternative to the steady-state method: derivation of reaction rates from first-passage times and pathway probabili-ties. Proc. Natl. Acad. Sci. U.S.A. 84, 663–667]. The principle idea of the Ninio approach is to track a single template/primer complex over time and to identify the expected behavior. The average time to insert a single nucleotide is a weighted sum of several terms, in-cluding the actual time to insert a nucleotide plus delays due to polymerase detachment from ei-ther the ternary (template-primer-polymerase) or quaternary (+nucleotide) complexes and time delays associated with the identification and ultimate rejection of an incorrect nucleotide from the binding site. The passage times of all events and their probability of occurrence are ex-pressed in terms of the rate constants of the elementary steps of the reaction pathway. The model accounts for variations in the average insertion time with different nucleotides as well as the in-fluence of G+C content of the sequence in the vicinity of the insertion site. Furthermore the model provides estimates of error frequencies. If nucleotide extension is recognized as a compe-tition between successful insertions and time delaying events, it can be described as a binomial process with a probability distribution. The distribution gives the probability to extend a primer/template complex with a certain number of base pairs and in general it maps annealed complexes into extension products.