Representation errors and retrievals in linear and nonlinear data assimilation

This article shows how one can formulate the representation problem starting from Bayes’ theorem. The purpose of this article is to raise awareness of the formal solutions,so that approximations can be placed in a proper context. The representation errors appear in the likelihood, and the different possibilities for the representation of reality in model and observations are discussed, including nonlinear representation probability density functions. Specifically, the assumptions needed in the usual procedure to add a representation error covariance to the error covariance of the observations are discussed,and it is shown that, when several sub-grid observations are present, their mean still has a representation error ; socalled ‘superobbing’ does not resolve the issue. Connection is made to the off-line or on-line retrieval problem, providing a new simple proof of the equivalence of assimilating linear retrievals and original observations. Furthermore, it is shown how nonlinear retrievals can be assimilated without loss of information. Finally we discuss how errors in the observation operator model can be treated consistently in the Bayesian framework, connecting to previous work in this area.

Diagonalization and Jordan Normal Form – motivation through Maple

Following an introduction to the diagonalization of matrices, one of the more difficult topics for students to grasp in linear algebra is the concept of Jordan normal form. In this note, we show how the important notions of diagonalization and Jordan normal form can be introduced and developed through the use of the computer algebra package Maple®.

Communicability across evolving networks

Many natural and technological applications generate time ordered sequences of networks, defined over a fixed set of nodes; for example time-stamped information about ‘who phoned who’ or ‘who came into contact with who’ arise naturally in studies of communication and the spread of disease. Concepts and algorithms for static networks do not immediately carry through to this dynamic setting. For example, suppose A and B interact in the morning, and then B and C interact in the afternoon. Information, or disease, may then pass from A to C, but not vice versa. This subtlety is lost if we simply summarize using the daily aggregate network given by the chain A-B-C. However, using a natural definition of a walk on an evolving network, we show that classic centrality measures from the static setting can be extended in a computationally convenient manner. In particular, communicability indices can be computed to summarize the ability of each node to broadcast and receive information. The computations involve basic operations in linear algebra, and the asymmetry caused by time’s arrow is captured naturally through the non-mutativity of matrix-matrix multiplication. Illustrative examples are given for both synthetic and real-world communication data sets. We also discuss the use of the new centrality measures for real-time monitoring and prediction.

On rank, underdetermined systems and consistency

The concepts of rank, underdetermined systems and consistency in linear algebra are discussed in the context of a puzzle. The article begins with a specific example, moving on to a generalization of the example and then to the general n x n case. As well as providing a solution of the puzzle, the article aims to provide students with a greater understanding of these abstract ideas in linear algebra through the study of the puzzle.

In vitro fermentation of linear and alpha-1,2-branched dextrans by the human fecal microbiota

The role of structure and molecular weight in fermentation selectivity in linear α-1,6 dextrans and dextrans with α-1,2 branching was investigated. Fermentation by gut bacteria was determined in anaerobic, pH-controlled fecal batch cultures after 36 h. Inulin (1%, wt/vol), which is a known prebiotic, was used as a control. Samples were obtained at 0, 10, 24, and 36 h of fermentation for bacterial enumeration by fluorescent in situ hybridization and short-chain fatty acid analyses. The gas production of the substrate fermentation was investigated in non-pH-controlled, fecal batch culture tubes after 36 h. Linear and branched 1-kDa dextrans produced significant increases in Bifidobacterium populations. The degree of α-1,2 branching did not influence the Bifidobacterium populations; however, α-1,2 branching increased the dietary fiber content, implying a decrease in digestibility. Other measured bacteria were unaffected by the test substrates except for the Bacteroides-Prevotella group, the growth levels of which were increased on inulin and 6- and 70-kDa dextrans, and the Faecalibacterium prausnitzii group, the growth levels of which were decreased on inulin and 1-kDa dextrans. A considerable increase in short-chain fatty acid concentration was measured following the fermentation of all dextrans and inulin. Gas production rates were similar among all dextrans tested but were significantly slower than that for inulin. The linear 1-kDa dextran produced lower total gas and shorter time to attain maximal gas production compared to those of the 70-kDa dextran (branched) and inulin. These findings indicate that dextrans induce a selective effect on the gut flora, short-chain fatty acids, and gas production depending on their length.

Linear nucleation of polymers

The use of high-melting fibres as linear nuclei for quiescent polymeric melts is instrumental in providing the superior mechanical properties of polymeric self-composites. It also has inherent advantages in the elucidation of fundamental aspects of polymeric crystallization and self-organization, not least in allowing systematic microscopic studies of polymeric crystallization from nucleation through to the growth interface. This has demonstrated explicitly that lamellae develop in two distinct ways, for slower and faster growth, depending on whether fold packing has or has not time to order before the next molecular layer is added with only the former leading to banded growth in linear polyethylene. Other gains in understanding concern cellulation and morphological instability, internuclear interference, isothermal lamellar thickening and banded growth being a consequence of the partial relief of initial surface stress. (c) 2006 Elsevier Ltd. All rights reserved.

Coriolis interactions about x-y axes in symmetric tops

Selection rules and matrix elements are derived for Coriolis interactions between vibrational levels due to rotation about (x, y) axes in symmetric top molecules. The theory is developed in detail for the case of interaction between an A1 and an E species vibrational level in a C3v molecule; perturbations to both the positions and the intensities of the rovibration transitions in the spectrum are considered. A computer program has been written which calculates exactly the perturbed spectrum of two interacting rovibration bands according to this model, the results being presented directly by a graph plotter connected to the computer. This has been used to interpret perturbations observed in two pairs of interacting fundamentals in the spectrum of CH3F (ν2 - ν5 and ν3 - ν6) and one pair in CD3Cl (ν2 - ν5). The resulting analysis of the observed spectrum leads to new values for some vibration-rotation interaction constants and also leads to a unique determination of the sign relationship between the dipole moment derivatives in each pair of interacting normal vibrations. These sign relations are summarized in Figs. 8, 12, and 15.

Potential energy surface intersections in triatomic molecules

The different types of surface intersection which may occur in linear configurations of triatomic molecules are reviewed, particularly with regard to the way in which the degeneracy is split as the molecule bends. The Renner-Teller effect in states of symmetry Π, Δ, Φ, etc., and intersections between Σ and Π, Σ and Δ, and Π and Δ states are discussed. A general method of modelling such intersecting potential surfaces is proposed, as a development of the model previously used by Murrell and Carter and co-workers for single-valued surfaces. Some of the lower energy surfaces of H2O, NH2, O3, C3, and HNO are discussed as examples.

Facile strategies for the synthesis and crystallization of linear trinuclear nickel(II)-Schiff base, complexes with carboxylate bridges: Tuning of coordination geometry and magnetic properties

Three new linear trinuclear nickel(II) complexes, [Ni-3(salpen)(2)(OAc)(2)(H2O)(2)]center dot 4H(2)O (1) (OAc = acetate, CH3COO-), [Ni-3(salpen)(2)(OBz)(2)] (2) (OBz=benzoate, PhCOO-) and [Ni-3(salpen)(2)(OCn)(2)(CH3CN)(2)] (4) (OCn = cinnamate, PhCH=CHCOO-), H(2)salpen = tetradentate ligand, N,N'-bis(salicylidene)-1,3-pentanediamine have been synthesized and characterized structurally and magnetically. The choice of solvent for growing single crystal was made by inspecting the morphology of the initially obtained solids with the help of SEM study. The magnetic properties of a closely related complex, [Ni-3(salpen)(2)(OPh)(2)(EtOH)] (3) (OPh = phenyl acetate, PhCH2COO-) whose structure and solution properties have been reported recently, has also been studied here. The structural analyses reveal that both phenoxo and carboxylate bridging are present in all the complexes and the three Ni(II) atoms remain in linear disposition. Although the Schiff base ligand and the syn-syn bridging bidentate mode of the carboxylate group remain the same in complexes 1-4, the change of alkyl/aryl group of the carboxylates brings about systematic variations between six- and five-coordination in the geometry of the terminal Ni(II) centres of the trinuclear units. The steric demand as well as hydrophobic nature of the alkyl/aryl group of the carboxylate is found to play a crucial role in the tuning of the geometry. Variable-temperature (2-300 K) magnetic susceptibility measurements show that complexes 1-4 are antiferromagnetically coupled (J = -3.2(1), -4.6(1). -3.2(1) and -2.8(1) cm(-1) in 1-4, respectively). Calculations of the zero-field splitting parameter indicate that the values of D for complexes 1-4 are in the high range (D = +9.1(2), +14.2(2), +9.8(2) and +8.6(1) cm(-1) for 1-4, respectively). The highest D value of +14.2(2) and +9.8(2) cm(-1) for complexes 2 and 3, respectively, are consistent with the pentacoordinated geometry of the two terminal nickel(II) ions in 2 and one terminal nickel(II) ion in 3. (C) 2009 Elsevier Ltd. All rights reserved.