Response curves of deterministic and probabilistic cellular automata in one and two dimensions

One of the most important problems in the theory of cellular automata (CA) is determining the proportion of cells in a specific state after a given number of time iterations. We approach this problem using patterns in preimage sets - that is, the set of blocks which iterate to the desired output. This allows us to construct a response curve - a relationship between the proportion of cells in state 1 after niterations as a function of the initial proportion. We derive response curve formulae for many two-dimensional deterministic CA rules with L-neighbourhood. For all remaining rules, we find experimental response curves. We also use preimage sets to classify surjective rules. In the last part of the thesis, we consider a special class of one-dimensional probabilistic CA rules. We find response surface formula for these rules and experimental response surfaces for all remaining rules.

Electrospinning and characterization of self-assembled inclusion complexies

L’électrofilage est une technique permettant de fabriquer des fibres polymériques dont le diamètre varie entre quelques nanomètres et quelques microns. Ces fibres ont donc un rapport surface/volume très élevé. Les fibres électrofilées pourraient trouver des applications dans le relargage de médicaments et le génie tissulaire, comme membranes et capteurs chimiques, ou dans les nanocomposites et dispositifs électroniques. L’électrofilage était initialement utilisé pour préparer des toiles de fibres désordonnées, mais il est maintenant possible d’aligner les fibres par l’usage de collecteurs spéciaux. Cependant, il est important de contrôler non seulement l’alignement macroscopique des fibres mais aussi leur orientation au niveau moléculaire puisque l’orientation influence les propriétés mécaniques, optiques et électriques des polymères. Les complexes moléculaires apparaissent comme une cible de choix pour produire des nanofibres fortement orientées. Dans les complexes d’inclusion d’urée, les chaînes polymères sont empilées dans des canaux unidimensionnels construits à partir d’un réseau tridimensionnel de molécules d’urée liées par des ponts hydrogène. Ainsi, les chaînes polymère sonts très allongées à l’échelle moléculaire. Des nanofibres du complexe PEO-urée ont été préparées pour la première fois par électrofilage de suspensions et de solutions. Tel qu’attendu, une orientation moléculaire inhabituellement élevée a été observée dans ces fibres. De tels complexes orientés pourraient être utilisés à la fois dans des études fondamentales et dans la préparation de matériaux hiérarchiquement structurés. La méthode d’électrofilage peut parfois aussi être utilisée pour préparer des matériaux polymériques métastables qui ne peuvent pas être préparés par des méthodes conventionnelles. Ici, l’électrofilage a été utilisé pour préparer des fibres des complexes stables (α) et "métastables" (β) entre le PEO et l’urée. La caractérisation du complexe β, qui était mal connu, révèle un rapport PEO:urée de 12:8 appartenant au système orthorhombique avec a = 1.907 nm, b = 0.862 nm et c = 0.773 nm. Les chaînes de PEO sont orientées selon l’axe de la fibre. Leur conformation est significativement affectée par les ponts hydrogène. Une structure en couches a été suggérée pour la forme β, plutôt que la structure conventionnelle en canaux adoptée par la forme α. Nos résultats indiquent que le complexe β est thermodynamiquement stable avant sa fonte et peut se transformer en forme α et en PEO liquide par un processus de fonte et recristallisation à 89 ºC. Ceci va dans le sens contraire aux observations faites avec le complexe β obtenu par trempe du complexe α fondu. En effet, le complexe β ainsi obtenu est métastable et contient des cristaux d’urée. Il peut subir une transition de phases cinétique solide-solide pour produire du complexe α dans une vaste gamme de températures. Cette transition est induite par un changement de conformation du PEO et par la formation de ponts hydrogène intermoléculaires entre l’urée et le PEO. Le diagramme de phases du système PEO-urée a été tracé sur toute la gamme de compositions, ce qui a permis d’interpréter la formation de plusieurs mélanges qui ne sont pas à l’équilibre mais qui sont été observés expérimentalement. La structure et le diagramme de phases du complexe PEO-thiourée, qui est aussi un complexe très mal connu, ont été étudiés en détail. Un rapport molaire PEO :thiourée de 3:2 a été déduit pour le complexe, et une cellule monoclinique avec a = 0.915 nm, b = 1.888 nm, c = 0.825 nm et β = 92.35º a été déterminée. Comme pour le complexe PEO-urée de forme β, une structure en couches a été suggérée pour le complexe PEO-thiourée, dans laquelle les molécules de thiourée seraient disposées en rubans intercalés entre deux couches de PEO. Cette structure en couches pourrait expliquer la température de fusion beaucoup plus faible des complexes PEO-thiourée (110 ºC) et PEO-urée de forme β (89 ºC) en comparaison aux structures en canaux du complexe PEO-urée de forme α (143 ºC).

Photoacoustic thermal characterization of Al2O3–Ag ceramic nanocomposites

Laser-induced nondestructive photoacoustic (PA) technique has been employed to determine the thermal diffusivity of nanometal (Ag) dispersed ceramic alumina matrix sintered at different temperatures. The thermal diffusivity values are evaluated by knowing the transition frequency from the amplitude spectrum of PA signal using the one-dimensional heat flow model of Rosencwaig and Gersho. Analysis of the data shows that heat transport and hence the thermal diffusivity value is greatly affected by the influence of incorporation of foreign atom. It is also seen that sintering temperature affects the thermal diffusivity value in a substantial manner. The results are interpreted in terms of variation in porosity and carrier-assisted heat transport mechanism in nanometal dispersed ceramics.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Ab Initio study of magnetic interactions in the KCuF3 and K2CuF4 low-dimensional Systems

The ab initio cluster model approach has been used to study the electronic structure and magnetic coupling of KCuF3 and K2CuF4 in their various ordered polytype crystal forms. Due to a cooperative Jahn-Teller distortion these systems exhibit strong anisotropies. In particular, the magnetic properties strongly differ from those of isomorphic compounds. Hence, KCuF3 is a quasi-one-dimensional (1D) nearest neighbor Heisenberg antiferromagnet whereas K2CuF4 is the only ferromagnet among the K2MF4 series of compounds (M=Mn, Fe, Co, Ni, and Cu) behaving all as quasi-2D nearest neighbor Heisenberg systems. Different ab initio techniques are used to explore the magnetic coupling in these systems. All methods, including unrestricted Hartree-Fock, are able to explain the magnetic ordering. However, quantitative agreement with experiment is reached only when using a state-of-the-art configuration interaction approach. Finally, an analysis of the dependence of the magnetic coupling constant with respect to distortion parameters is presented.

From a discrete to a continuum model of cell dynamics in one dimension

Multiscale modeling is emerging as one of the key challenges in mathematical biology. However, the recent rapid increase in the number of modeling methodologies being used to describe cell populations has raised a number of interesting questions. For example, at the cellular scale, how can the appropriate discrete cell-level model be identified in a given context? Additionally, how can the many phenomenological assumptions used in the derivation of models at the continuum scale be related to individual cell behavior? In order to begin to address such questions, we consider a discrete one-dimensional cell-based model in which cells are assumed to interact via linear springs. From the discrete equations of motion, the continuous Rouse [P. E. Rouse, J. Chem. Phys. 21, 1272 (1953)] model is obtained. This formalism readily allows the definition of a cell number density for which a nonlinear "fast" diffusion equation is derived. Excellent agreement is demonstrated between the continuum and discrete models. Subsequently, via the incorporation of cell division, we demonstrate that the derived nonlinear diffusion model is robust to the inclusion of more realistic biological detail. In the limit of stiff springs, where cells can be considered to be incompressible, we show that cell velocity can be directly related to cell production. This assumption is frequently made in the literature but our derivation places limits on its validity. Finally, the model is compared with a model of a similar form recently derived for a different discrete cell-based model and it is shown how the different diffusion coefficients can be understood in terms of the underlying assumptions about cell behavior in the respective discrete models.

Stabilization of a helical water chain in a metal-organic host of a trinuclear Schiff base complex

Three heterometallic trinuclear Schiff base complexes, [{GuL(1)(H2O)}(2)Ni(CN)(4)]center dot 4H(2)O (1), [{CuL2(H2O)}(2)Ni(CN)(4)] (2), and [{CuL3(H2O)}(2)Ni(CN)(4)] (3) (HL1 = 7-amino-4-methyl-5-azahept-3-en-2-one, HL2 = 7-methylamino-4-methyl-5-azahept-3-en-2-one, and HL3 = 7-dimethylamino-4-methyl-5-azahept-3-en-2-one), were synthesized. All three complexes were characterized by elemental analysis, IR and UV spectroscopies, and thermal analysis. Two of them (1 and 3) were also characterized by single crystal X-ray crystallography. Complex 1 forms a hydrogen-bonded one-dimensional metal-organic framework that stabilizes a helical water chain into its cavity, but when any of the amine hydrogen atoms of the Schiff base are replaced by methyl groups, as in L 2 and L 3, the water chain, vanishes, showing explicitly the importance of the host-guest H-bonding interactions for the stabilization of a water cluster.

Facile synthesis of Cu(II) complexes of monocondensed N,N,N donor Schiff base ligands: crystal structure, spectroscopic and magnetic properties

Four new copper(II) complexes, [((CuLN3)-N-1)(2)](ClO4)(2) (1), [(CuL2 N-3)(2)](ClO4)(2) (2), [CuL3(N-3)ClO4)](n) (3) and [CuL4(mu-1,1-N-3)(mu-1,3-N-3)(ClO4)](n) (4) where L-1 = N-1-pyridin-2-yl-methylene-propane-1,3-diamine, L-2 = N-1-(1-pyridin-2-yl-ethylidene)propane-1,3-diamine, L-3 =N-1-(1-pyridin-2-yl-ethylidene)ethane-1,2-diamine and L-4=N-1-(1-pyridin-2-yl-ethylidene)propane-1,2-diamine are four tridentate N,N,N donor Schiff base ligands, have been derived and structurally characterized by X-ray crystallography. Compounds 1 and 2 consist of double basal-apical end-on (EO) azide bridged dinuclear Cu-II complexes with square-pyramidal geometry. In complex 3 the square planar mononuclear [CuL3 (N-3)] units are linked by weakly coordinated perchlorate ions in the axial positions of Cu-II to form a one-dimensional chain. Two such chains are connected by hydrogen bonds involving perchlorate ions and azide groups. Compound 4 consists of 1-D chains in which the Cu-II ions with a square-pyramidal geometry are alternately bridged by single EO and end-to-end (EE) azido ligands, both adopting a basal-apical disposition. Variable temperature (300-2 K) magnetic susceptibility measurements and magnetization measurements at 2 K have been performed. The results reveal that complexes 1 and 2 are antiferromagnetically coupled through azido bridges (J= -12.18 +/- 0.09 and -4.43 +/- 0.1 cm(-1) for 1 and 2, respectively). Complex 3 shows two different magnetic interactions through the two kinds of hydrogen bonds; one is antiferromagnetic (J(1) = - 9.69 +/- 0.03 cm(-1)) and the other is ferromagnetic (J(2) = 1.00 +/- 0.01 cm(-1)). From a magnetic point of view complex 4 is a ferromagnetic dinuclear complex (J= 1.91 +/- 0.01 cm(-1)) coupled through the EO bridge only. The coupling through the EE bridge is practically nil as the N(azido)-Cu-II (axial) distance (2.643 angstrom) is too long. (C) 2006 Elsevier Ltd. All rights reserved.

Zinc(II) mediated synthesis of an N-substituted 2-(2-pyridyl)imidazole from a 1,2-diketone and 2-(aminomethyl)pyridine and its ligational behaviour

Reaction of anhydrous ZnCl2 with the 1:2 condensate (L) of benzil and 2-(aminomethyl)pyridine in methanol gives monomeric ZnL'Cl-2 (1) where L' is 2-[(4,5-diphenyl-2-pyridin-2-yl-1H-imidazol-1-yl)-methyl]pyridine. In the X-ray crystal structure, 1 is found to contain tetrahedral zinc with an N2Cl2 coordination sphere and the N-substituent methylpyridine fragment hanging as a free arm. A tentative mechanism is proposed for the zinc mediated conversion of L-->L'. Demetallation of 1 by the action of aqueous NaOH yields L' in the free state. When L' is reacted with Zn(ClO4)(2).6H(2)O in a 1:2 molar proportion, [Zn(L')(2)](n)(ClO4)(2n).(H2O)(n/2).(CH2Cl2)(n/2) (2) is obtained. The zinc atom in 2, as revealed by X-ray crystallography, has a trigonal bipyramidal N-5 coordination sphere. There are two independent ligands in the asymmetric unit of 2. One of them bonds only to one zinc atom in a bidentate mode with the N-substituent methylpyridine hanging free while the other ligand binds to two different zinc atoms in a tridentate fashion, employing the N-substituent methylpyridine nitrogen atom to form the polymeric one-dimensional chain cation.

Novel charge transfer supramolecular assemblies with Keggin anions and 2-amino-5-nitropyridine

Several novel compounds with the non-linear optical chromophore 2-amino-5-nitropyridine (2A5NP) and Keggin polyoxoanions (alpha-isomers), having the general formula (2A5NP)(m)H-n[XM12O40]center dot xH(2)O, M = Mo, W, were synthesised. Compounds were obtained with X = P, n = 3, m = 3 and 4 and X = Si, n = m = 4 ( x = 2 - 6). Thus, for each of the anions [PMo12O40](3-) and [PW12O40](3-) two different compounds were obtained, with the same anion and organic counterpart but with a different stoichiometric ratio. These presented different charge transfer properties and thermal stability. All compounds were characterised by spectroscopic and analytical techniques. The single crystal X-ray diffraction structure of (2A5NP)(4)H-3[PMo12O40]center dot 2.5H(2)O center dot 0.5C(2)H(5)OH showed that the water solvent molecules and the organic chromophores are assembled via infinite one-dimensional chains of hydrogen bonds with formation of open channels, which accommodate [ PMo12O40] 3- and ethanol solvent molecules.

Structural variations in Ni(II) complexes of salen type di-Schiff base ligands

Three di-Schiff-base ligands, N,N'-bis(salicylidene)-1,3-propanediamine (H(2)Salpn), N,N'-bis(salicylidene)-1,3-pentanedianiine (H(2)Salpen) and N,N'-bis(salicylidine)-ethylenediamine (H(2)Salen) react with Ni(SCN)(2). 4H(2)O in 2:3 molar ratios to form the complexes; mononuclear [Ni(HSalpn)(NCS)(H2O)]center dot H2O (1a), trinuclear [{Ni(Salpen)}(2)Ni(NCS)(2)] (2b) and trinuclear [{Ni(Salen)}(2)Ni(NCS)(2)] (3) respectively. All the complexes have been characterized by elemental analyses, IR and UV-VIS spectra, and room temperature magnetic susceptibility measurements. The structures of la and 2b have been confirmed by X-ray single crystal analysis. In complex la, the Ni(II) atom is coordinated equatorially by the tetradentate, mononegative Schiff-base, HSalpn. Axial coordination of isothiocyanate group and a water molecule completes its octahedral geometry. The hydrogen atom attached to one of the oxygen atoms of the Schiff base is involved in a very strong hydrogen bond with a neighboring unit to form a centrosymmetric dimer. In 2b, two square planar [Ni(Salpen)] units act as bide mate oxygen donor ligands to a central Ni(II) which is also coordinated by two mutually cis N-bonded thiocyanate ligands to complete its distorted octahedral geometry. Complex 3 possesses a similar structure to that of 2b. A dehydrated form of la and a hydrated form of 2b have been obtained and characterized. The importance of electronic and steric factors in the variation of the structures is discussed. (c) 2007 Elsevier Ltd. All rights reserved.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

The power of NMR. Part 1: the beginnings

Nuclear mnagnetic resonance (NMR) spectroscopy involves the excitation of nuclei by electromagnetic radiation in the radio-frequency range of the electromagnetic spectrum. For a nucleus to absorb energy from radiowaves in this way, it must hve the quantum mechanical property of spin. A spinning nucleus, such as that of the hydrogen atom, will dopt one f only two possible states when placed in a magnetic field. (In NMR, the hydrogen nucleus is often referred to as a proton, and is given the abbreviation 1H.) Az the strength of the magnetic field is increased, there is a proportional increase in the energy 'gap' between these two states. We can predic the resonant frequency at which any spinning nucleus will absorb energy from radio-frequency radiation as it jumps from the lower energy state to the upper state.

Use of a reduced Schiff-Base ligand to prepare novel chloro-bridged chains of rare Cu(II) trinuclear complexes with mixed azido/oxo and chloro/oxo bridges

Two mixed bridged one-dimensional (1D) polynuclear complexes, [Cu3L2(mu(1,1)-N-3)(2)(mu-Cl)Cl](n) (1) and {[Cu3L2(mu-Cl)(3)Cl]center dot 0.46CH(3)OH}(n), (2), have been synthesized using the tridentate reduced Schiff-base ligand HL (2-[(2-dimethylamino-ethylamino)-methyl]-phenol). The complexes have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. In both complexes the basic trinuclear angular units are joined together by weak chloro bridges to form a 1D chain. The trinuclear structure of 1 is composed of two terminal square planar [Cu(L)(mu(1,1)-N-3)] units connected by a central Cu(II) atom through bridging nitrogen atoms of end-on azido ligands and the phenoxo oxygen atom of the tridentate ligand. These four coordinating atoms along with a chloride ion form a distorted trigonal bipyramidal geometry around the central Cu(II). The structure of 2 is similar; the only difference being a Cl bridge replacing the mu(1,1)-N-3 bridge in the trinuclear unit. The magnetic properties of both trinuclear complexes can be very well reproduced with a simple linear symmetrical trimer model (H = JS(i)S(i+1)) with only one intracluster exchange coupling (J) including a weak intertrimer interaction (.j) reproduced with the molecular field approximation. This model provides very satisfactory fits for both complexes in the whole temperature range with the following parameters: g = 2.136(3), J = 93.9(3) cm(-1) and zj= -0.90(3) cm(-1) (z = 2) for 1 and g = 2.073(7), J = -44.9(4) cm(-1) and zJ = -1.26(6) cm(-1) (z = 2) for 2.

Four-dimensional variational data assimilation for high resolution nested models

Four-dimensional variational data assimilation (4D-Var) is used in environmental prediction to estimate the state of a system from measurements. When 4D-Var is applied in the context of high resolution nested models, problems may arise in the representation of spatial scales longer than the domain of the model. In this paper we study how well 4D-Var is able to estimate the whole range of spatial scales present in one-way nested models. Using a model of the one-dimensional advection–diffusion equation we show that small spatial scales that are observed can be captured by a 4D-Var assimilation, but that information in the larger scales may be degraded. We propose a modification to 4D-Var which allows a better representation of these larger scales.