On cluster approximations of cellular automata

In this thesis I examine one commonly used class of methods for the analytic approximation of cellular automata, the so-called local cluster approximations. This class subsumes the well known mean-field and pair approximations, as well as higher order generalizations of these. While a straightforward method known as Bayesian extension exists for constructing cluster approximations of arbitrary order on one-dimensional lattices (and certain other cases), for higher-dimensional systems the construction of approximations beyond the pair level becomes more complicated due to the presence of loops. In this thesis I describe the one-dimensional construction as well as a number of approximations suggested for higher-dimensional lattices, comparing them against a number of consistency criteria that such approximations could be expected to satisfy. I also outline a general variational principle for constructing consistent cluster approximations of arbitrary order with minimal bias, and show that the one-dimensional construction indeed satisfies this principle. Finally, I apply this variational principle to derive a novel consistent expression for symmetric three cell cluster frequencies as estimated from pair frequencies, and use this expression to construct a quantitatively improved pair approximation of the well-known lattice contact process on a hexagonal lattice.

Soluautomaattien approksimaatioista

Tässä tutkielmassa tarkastelen soluautomaatteja ja niiden yleisesti käytettyj ä keskivertokenttä- ("mean field") ja pariapproksimaatioita. Määrittelen soluautomaatin käsitteen yleisellä tasolla ja kuvailen kolme eri tyypillisesti käytettyä soluautomaattityyppiä, joista kahdesta esimerkkeinä mainitsen John Conwayn tunnetun "Game of Life" -soluautomaatin sekä epidemiologian alalla yleisesti käytetyn hilakontaktiprosessin. Tutkielman toisessa osassa esittelen eri tapoja, joilla soluautomaattien käyttäytymistä voidaan analyyttisesti approksimoida ja karakterisoida, mukaanlukien niin kutsuttu soluautomaatin aktiviteetti (λ) sekä yllämainitut keskivertokenttä- ja pariapproksimaatiot. Kahta jälkimmäistä sovellan hilakontaktiprosessiin, ja vertailen näin saatuja tuloksia.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

Phonon-mediated two-magnon light scattering from ferromagnetic insulators

An indirect mechanism of light scattering from spin-waves in ferromagnetic insulators via two-magnon one-phonon process is proposed. Following linear response theory, an expression has been derived for the differential scattering cross-section in the mean-field-approximation.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Novel phase-transition behavior in an aqueous electrolyte solution

We have investigated the near-critical behavior of the susceptibility of a ternary liquid mixture of 3-methylpyridine. water, and sodium bromide as a function of the salt concentration. The susceptibility was determined from light-scattering measurements performed at a scattering angle of 90 degrees in the one-phase region near the locus of lower consolute points. A sharp crossover from asymptotic Ising behavior to mean-field behavior has been observed at concentrations ranging from 8 to 16.5 mass% NaBr. The range of asymptotic Ising behavior shrinks with increasing salt concentration and vanishes at a NaBr concentration of about 17 mass%. where complete mean-field-like behavior of the susceptibility is observed. A simultaneous pronounced increase in the background scattering at concentrations above 15 mass%, as well as a dip in the critical locus at 17 mass % NaBr, suggests that this phenomenon can be interpreted as mean-field tricritical behavior associated with the formation of a microheterogeneous phase due to clustering of the molecules and ions. An analogy with tri critical behavior observed in polymer solutions as well as the possibility of a charge-density-wave phase is also discussed. In addition, we, have observed a third soap-like phase an the liquid-liquid interface in several binary and ternary liquid mixtures.

Modeling server-unreliability in closed queuing-networks

A method is presented to model server unreliability in closed queuing networks. Breakdowns and repairs of servers, assumed to be time-dependent, are modeled using virtual customers and virtual servers in the system. The problem is thus converted into a closed queue with all reliable servers and preemptive resume priority centers. Several recent preemptive priority approximations and an approximation of the one proposed are used in the analysis. This method has approximately the same computational requirements as that of mean-value analysis for a network of identical dimensions and is therefore very efficient

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some of the current issues of the solar simulations are put forward.

Search on a hypercubic lattice using a quantum random walk. I. d > 2

Random walks describe diffusion processes, where movement at every time step is restricted to only the neighboring locations. We construct a quantum random walk algorithm, based on discretization of the Dirac evolution operator inspired by staggered lattice fermions. We use it to investigate the spatial search problem, that is, to find a marked vertex on a d-dimensional hypercubic lattice. The restriction on movement hardly matters for d > 2, and scaling behavior close to Grover's optimal algorithm (which has no restriction on movement) can be achieved. Using numerical simulations, we optimize the proportionality constants of the scaling behavior, and demonstrate the approach to that for Grover's algorithm (equivalent to the mean-field theory or the d -> infinity limit). In particular, the scaling behavior for d = 3 is only about 25% higher than the optimal d -> infinity value.

Metal-Insulator Transition in the Hubbard Model on a Triangular Lattice

We discuss the results of an extensive mean-field investigation of the half-filled Hubbard model on a triangular lattice at zero temperature. At intermediate U we find a first-order metal-insulator transition from an incommensurate spiral magnetic metal to a semiconducting state with a commensurate linear spin density wave ordering stabilized by the competition between the kinetic energy and the frustrated nature of the magnetic interaction. At large U the ground state is that of a classical triangular antiferromagnet within our approximation. In the incommensurate spiral metallic phase the Fermi surface has parts in which the wave function renormalization Z is extremely small. The evolution of the Fermi surface and the broadening of the quasi-particle band along with the variation of the plasma frequency and a charge stiffness constant with U/t are discussed.

Spiral states in the square-lattice Hubbard model

We present a variety of physical implications of a mean-field theory for spiral spin-density-wave states in the square-lattice Hubbard model for small deviations from half filling. The phase diagram with the paramagnetic metal, two spiral (semimetallic) states, and ferromagnet is calculated. The momentum distribution function and the (quasiparticle) density of states are discussed. There is a significant broadening of the quasiparticle bands when the antiferromagnetic insulator is doped. The evolution of the Fermi surface and the variation of the plasma frequency and a charge-stiffness constant with U/t and δ are calculated. The connection to results based on the Schwinger-boson-slave-fermion formalism is made.

The shear dynamo problem for small magnetic Reynolds numbers

We study large-scale kinematic dynamo action due to turbulence in the presence of a linear shear flow in the low-conductivity limit. Our treatment is non-perturbative in the shear strength and makes systematic use of both the shearing coordinate transformation and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed to have low magnetic Reynolds number (Re-m), but could have arbitrary fluid Reynolds number. The equation for the magnetic fluctuations is expanded perturbatively in the small quantity, Re-m. Our principal results are as follows: (i) the magnetic fluctuations are determined to the lowest order in Rem by explicit calculation of the resistive Green's function for the linear shear flow; (ii) the mean electromotive force is then calculated and an integro-differential equation is derived for the time evolution of the mean magnetic field. In this equation, velocity fluctuations contribute to two different kinds of terms, the 'C' and 'D' terms, respectively, in which first and second spatial derivatives of the mean magnetic field, respectively, appear inside the space-time integrals; (iii) the contribution of the D term is such that its contribution to the time evolution of the cross-shear components of the mean field does not depend on any other components except itself. Therefore, to the lowest order in Re-m, but to all orders in the shear strength, the D term cannot give rise to a shear-current-assisted dynamo effect; (iv) casting the integro-differential equation in Fourier space, we show that the normal modes of the theory are a set of shearing waves, labelled by their sheared wavevectors; (v) the integral kernels are expressed in terms of the velocity-spectrum tensor, which is the fundamental dynamical quantity that needs to be specified to complete the integro-differential equation description of the time evolution of the mean magnetic field; (vi) the C term couples different components of the mean magnetic field, so they can, in principle, give rise to a shear-current-type effect. We discuss the application to a slowly varying magnetic field, where it can be shown that forced non-helical velocity dynamics at low fluid Reynolds number does not result in a shear-current-assisted dynamo effect.

Theory of the non-Fermi-liquid transition point in the two-impurity Kondo model

We present an explicit solution of the problem of two coupled spin-1/2 impurities, interacting with a band of conduction electrons. We obtain an exact effective bosonized Hamiltonian, which is then treated by two different methods (low-energy theory and mean-field approach). Scale invariance is explicitly shown at the quantum critical point. The staggered susceptibility behaves like ln(T(K)/T) at low T, whereas the magnetic susceptibility and [S1.S2] are well behaved at the transition. The divergence of C(T)/T when approaching the transition point is also studied. The non-Fermi-liquid (actually marginal-Fermi-liquid) critical point is shown to arise because of the existence of anomalous correlations, which lead to degeneracies between bosonic and fermionic states of the system. The methods developed in this paper are of interest for studying more physically relevant models, for instance, for high-T(c) cuprates.

A Cache coherence protocol for MIN-based multiprocessors

In this paper we present a cache coherence protocol for multistage interconnection network (MIN)-based multiprocessors with two distinct private caches: private-blocks caches (PCache) containing blocks private to a process and shared-blocks caches (SCache) containing data accessible by all processes. The architecture is extended by a coherence control bus connecting all shared-block cache controllers. Timing problems due to variable transit delays through the MIN are dealt with by introducing Transient states in the proposed cache coherence protocol. The impact of the coherence protocol on system performance is evaluated through a performance study of three phases. Assuming homogeneity of all nodes, a single-node queuing model (phase 3) is developed to analyze system performance. This model is solved for processor and coherence bus utilizations using the mean value analysis (MVA) technique with shared-blocks steady state probabilities (phase 1) and communication delays (phase 2) as input parameters. The performance of our system is compared to that of a system with an equivalent-sized unified cache and with a multiprocessor implementing a directory-based coherence protocol. System performance measures are verified through simulation.