A Statistical Model of Current Loops and Magnetic Monopoles

We formulate a natural model of loops and isolated vertices for arbitrary planar graphs, which we call the monopole-dimer model. We show that the partition function of this model can be expressed as a determinant. We then extend the method of Kasteleyn and Temperley-Fisher to calculate the partition function exactly in the case of rectangular grids. This partition function turns out to be a square of a polynomial with positive integer coefficients when the grid lengths are even. Finally, we analyse this formula in the infinite volume limit and show that the local monopole density, free energy and entropy can be expressed in terms of well-known elliptic functions. Our technique is a novel determinantal formula for the partition function of a model of isolated vertices and loops for arbitrary graphs.

A consistent model of isolated force-free magnetic arches

The magnetic flux tube concentrating strong magnetic field is the basic configuration of magneticfield in the solar atmosphere. In the present paper, the equilibrium of isolated magnetic flux tube inthe solar atmosphere is discussed. In the viewpoint of mathematics, the boundary condition is nonlinearand the position of boundary needs to be determined by the physical condition although the equation ofmagnetic potential is linear for the linear force-free field. Analytical solutions to the arches of bothuniform circular cross-section and non-uniform cross section have been obtained. The results show thatthe nonlinear problem may have or not have any solution according to different azimuthal components of the magnetic field; the number of solutions to the nonlinear problem is four at most, and two in some cases. In the present paper, the analytical solutions to the approximations of both fat and slender arches are given in detail, and the general features of magnetic arch structure are shown.

Kinetic model of continuous-wave flow chemical lasers

In this paper an analysis of the kinetic theory of the continuous-wave flow chemical lasers(CWFCL) is presented with emphasis being laid on the effects of inhomogeneous broadeningon CWFCL's performance. The results obtained are applicable to the case where laser fre-quency is either coincident or incoincident with that of the eenter of the line shape. This rela-tion has been,compared with that of the rate model in common use. These two models are almostidentical as the broadening parameter η is larger than 1. The smaller the value of η, thegreater the difference between the results of these two models will be. For fixed η, the dif-ferences between fhe results of the two models increase with the increase of the frequencyshift parameter ξ. When η is about less than 0.2. the kinetic model can predict exactly the in-homogeneous broadening effects,while the rate model cannot.

A Model of Ischemia-Induced Neuroblast Activation in the Adult Subventricular Zone

12 p.

Early astrocytic atrophy in the entorhinal cortex of a triple transgenic animal model of Alzheimer's disease

The EC (entorhinal cortex) is fundamental for cognitive and mnesic functions. Thus damage to this area appears as a key element in the progression of AD (Alzheimer's disease), resulting in memory deficits arising from neuronal and synaptic alterations as well as glial malfunction. In this paper, we have performed an in-depth analysis of astroglial morphology in the EC by measuring the surface and volume of the GFAP (glial fibrillary acidic protein) profiles in a triple transgenic mouse model of AD [3xTg-AD (triple transgenic mice of AD)]. We found significant reduction in both the surface and volume of GFAP-labelled profiles in 3xTg-AD animals from very early ages (1 month) when compared with non-Tg (non-transgenic) controls (48 and 54%, reduction respectively), which was sustained for up to 12 months (33 and 45% reduction respectively). The appearance of Lambda beta (amyloid beta-peptide) depositions at 12 months of age did not trigger astroglial hypertrophy; nor did it result in the close association of astrocytes with senile plaques. Our results suggest that the AD progressive cognitive deterioration can be associated with an early reduction of astrocytic arborization and shrinkage of the astroglial domain, which may affect synaptic connectivity within the EC and between the EC and other brain regions. In addition, the EC seems to be particularly vulnerable to AD pathology because of the absence of evident astrogliosis in response to A beta accumulation. Thus we can consider that targeting astroglial atrophy may represent a therapeutic strategy which might slow down the progression of AD.

TUNPOP, a computer simulation model of the yellowfin tuna population and the surface tuna fishery of the Eastern Pacific Ocean

ENGLISH: Mathematical documentation of TUNP0P, an age-structured computer simulation model of the yellowfin tuna population and surface tuna fishery of the eastern Pacific Ocean, is described. Example runs of the model are presented and discussed, and the sensitivity of the model output to changes in various parameters is examined. SPANISH: Se describe la documentación matemática de TUNP0P, un modelo computador de simulación basado en la edad de la población del atún aleta amarilla y de la pesca atunera epipelágíca del Océano Pacífico oriental. Se presentan y se discuten ejemplos de las pasadas del modelo, y se examina la sensibilidad de los resultados de salida con relación a los cambios de varios parámetros. (PDF contains 47 pages.)

Absolute measurement of roughness and lateral-correlation length of random surfaces by use of the simplified model of image-speckle contrast

We present a method of image-speckle contrast for the nonprecalibration measurement of the root-mean-square roughness and the lateral-correlation length of random surfaces with Gaussian correlation. We use the simplified model of the speckle fields produced by the weak scattering object in the theoretical analysis. The explicit mathematical relation shows that the saturation value of the image-speckle contrast at a large aperture radius determines the roughness, while the variation of the contrast with the aperture radius determines the lateral-correlation length. In the experimental performance, we specially fabricate the random surface samples with Gaussian correlation. The square of the image-speckle contrast is measured versus the radius of the aperture in the 4f system, and the roughness and the lateral-correlation length are extracted by fitting the theoretical result to the experimental data. Comparison of the measurement with that by an atomic force microscope shows our method has a satisfying accuracy. (C) 2002 Optical Society of America.

Mathematical study of complex networks : brain, internet, and power grid

The dissertation is concerned with the mathematical study of various network problems. First, three real-world networks are considered: (i) the human brain network (ii) communication networks, (iii) electric power networks. Although these networks perform very different tasks, they share similar mathematical foundations. The high-level goal is to analyze and/or synthesis each of these systems from a “control and optimization” point of view. After studying these three real-world networks, two abstract network problems are also explored, which are motivated by power systems. The first one is “flow optimization over a flow network” and the second one is “nonlinear optimization over a generalized weighted graph”. The results derived in this dissertation are summarized below.

Brain Networks: Neuroimaging data reveals the coordinated activity of spatially distinct brain regions, which may be represented mathematically as a network of nodes (brain regions) and links (interdependencies). To obtain the brain connectivity network, the graphs associated with the correlation matrix and the inverse covariance matrix—describing marginal and conditional dependencies between brain regions—have been proposed in the literature. A question arises as to whether any of these graphs provides useful information about the brain connectivity. Due to the electrical properties of the brain, this problem will be investigated in the context of electrical circuits. First, we consider an electric circuit model and show that the inverse covariance matrix of the node voltages reveals the topology of the circuit. Second, we study the problem of finding the topology of the circuit based on only measurement. In this case, by assuming that the circuit is hidden inside a black box and only the nodal signals are available for measurement, the aim is to find the topology of the circuit when a limited number of samples are available. For this purpose, we deploy the graphical lasso technique to estimate a sparse inverse covariance matrix. It is shown that the graphical lasso may find most of the circuit topology if the exact covariance matrix is well-conditioned. However, it may fail to work well when this matrix is ill-conditioned. To deal with ill-conditioned matrices, we propose a small modification to the graphical lasso algorithm and demonstrate its performance. Finally, the technique developed in this work will be applied to the resting-state fMRI data of a number of healthy subjects.

Communication Networks: Congestion control techniques aim to adjust the transmission rates of competing users in the Internet in such a way that the network resources are shared efficiently. Despite the progress in the analysis and synthesis of the Internet congestion control, almost all existing fluid models of congestion control assume that every link in the path of a flow observes the original source rate. To address this issue, a more accurate model is derived in this work for the behavior of the network under an arbitrary congestion controller, which takes into account of the effect of buffering (queueing) on data flows. Using this model, it is proved that the well-known Internet congestion control algorithms may no longer be stable for the common pricing schemes, unless a sufficient condition is satisfied. It is also shown that these algorithms are guaranteed to be stable if a new pricing mechanism is used.

Electrical Power Networks: Optimal power flow (OPF) has been one of the most studied problems for power systems since its introduction by Carpentier in 1962. This problem is concerned with finding an optimal operating point of a power network minimizing the total power generation cost subject to network and physical constraints. It is well known that OPF is computationally hard to solve due to the nonlinear interrelation among the optimization variables. The objective is to identify a large class of networks over which every OPF problem can be solved in polynomial time. To this end, a convex relaxation is proposed, which solves the OPF problem exactly for every radial network and every meshed network with a sufficient number of phase shifters, provided power over-delivery is allowed. The concept of “power over-delivery” is equivalent to relaxing the power balance equations to inequality constraints.

Flow Networks: In this part of the dissertation, the minimum-cost flow problem over an arbitrary flow network is considered. In this problem, each node is associated with some possibly unknown injection, each line has two unknown flows at its ends related to each other via a nonlinear function, and all injections and flows need to satisfy certain box constraints. This problem, named generalized network flow (GNF), is highly non-convex due to its nonlinear equality constraints. Under the assumption of monotonicity and convexity of the flow and cost functions, a convex relaxation is proposed, which always finds the optimal injections. A primary application of this work is in the OPF problem. The results of this work on GNF prove that the relaxation on power balance equations (i.e., load over-delivery) is not needed in practice under a very mild angle assumption.

Generalized Weighted Graphs: Motivated by power optimizations, this part aims to find a global optimization technique for a nonlinear optimization defined over a generalized weighted graph. Every edge of this type of graph is associated with a weight set corresponding to the known parameters of the optimization (e.g., the coefficients). The motivation behind this problem is to investigate how the (hidden) structure of a given real/complex valued optimization makes the problem easy to solve, and indeed the generalized weighted graph is introduced to capture the structure of an optimization. Various sufficient conditions are derived, which relate the polynomial-time solvability of different classes of optimization problems to weak properties of the generalized weighted graph such as its topology and the sign definiteness of its weight sets. As an application, it is proved that a broad class of real and complex optimizations over power networks are polynomial-time solvable due to the passivity of transmission lines and transformers.

A model of the olfactory bulb and beyond

The olfactory bulb of mammals aids in the discrimination of odors. A mathematical model based on the bulbar anatomy and electrophysiology is described. Simulations of the highly non-linear model produce a 35-60 Hz modulated activity, which is coherent across the bulb. The decision states (for the odor information) in this system can be thought of as stable cycles, rather than as point stable states typical of simpler neuro-computing models. Analysis shows that a group of coupled non-linear oscillators are responsible for the oscillatory activities. The output oscillation pattern of the bulb is determined by the odor input. The model provides a framework in which to understand the transformation between odor input and bulbar output to the olfactory cortex. This model can also be extended to other brain areas such as the hippocampus, thalamus, and neocortex, which show oscillatory neural activities. There is significant correspondence between the model behavior and observed electrophysiology.

It has also been suggested that the olfactory bulb, the first processing center after the sensory cells in the olfactory pathway, plays a role in olfactory adaptation, odor sensitivity enhancement by motivation, and other olfactory psychophysical phenomena. The input from the higher olfactory centers to the inhibitory cells in the bulb are shown to be able to modulate the response, and thus the sensitivity, of the bulb to odor input. It follows that the bulb can decrease its sensitivity to a pre-existing and detected odor (adaptation) while remaining sensitive to new odors, or can increase its sensitivity to discover interesting new odors. Other olfactory psychophysical phenomena such as cross-adaptation are also discussed.

Exploring the Kinetics of Domain Switching in Ferroelectrics for Structural Applications

The complex domain structure in ferroelectrics gives rise to electromechanical coupling, and its evolution (via domain switching) results in a time-dependent (i.e. viscoelastic) response. Although ferroelectrics are used in many technological applications, most do not attempt to exploit the viscoelastic response of ferroelectrics, mainly due to a lack of understanding and accurate models for their description and prediction. Thus, the aim of this thesis research is to gain better understanding of the influence of domain evolution in ferroelectrics on their dynamic mechanical response. There have been few studies on the viscoelastic properties of ferroelectrics, mainly due to a lack of experimental methods. Therefore, an apparatus and method called Broadband Electromechanical Spectroscopy (BES) was designed and built. BES allows for the simultaneous application of dynamic mechanical and electrical loading in a vacuum environment. Using BES, the dynamic stiffness and loss tangent in bending and torsion of a particular ferroelectric, viz. lead zirconate titanate (PZT), was characterized for different combinations of electrical and mechanical loading frequencies throughout the entire electric displacement hysteresis. Experimental results showed significant increases in loss tangent (by nearly an order of magnitude) and compliance during domain switching, which shows promise as a new approach to structural damping. A continuum model of the viscoelasticity of ferroelectrics was developed, which incorporates microstructural evolution via internal variables and associated kinetic relations. For the first time, through a new linearization process, the incremental dynamic stiffness and loss tangent of materials were computed throughout the entire electric displacement hysteresis for different combinations of mechanical and electrical loading frequencies. The model accurately captured experimental results. Using the understanding gained from the characterization and modeling of PZT, two applications of domain switching kinetics were explored by using Micro Fiber Composites (MFCs). Proofs of concept of set-and-hold actuation and structural damping using MFCs were demonstrated.

Tsunamis - a model of their generation and propagation

A general solution is presented for water waves generated by an arbitrary movement of the bed (in space and time) in a two-dimensional fluid domain with a uniform depth. The integral solution which is developed is based on a linearized approximation to the complete (nonlinear) set of governing equations. The general solution is evaluated for the specific case of a uniform upthrust or downthrow of a block section of the bed; two time-displacement histories of the bed movement are considered.

An integral solution (based on a linear theory) is also developed for a three-dimensional fluid domain of uniform depth for a class of bed movements which are axially symmetric. The integral solution is evaluated for the specific case of a block upthrust or downthrow of a section of the bed, circular in planform, with a time-displacement history identical to one of the motions used in the two-dimensional model.

Since the linear solutions are developed from a linearized approximation of the complete nonlinear description of wave behavior, the applicability of these solutions is investigated. Two types of non-linear effects are found which limit the applicability of the linear theory: (1) large nonlinear effects which occur in the region of generation during the bed movement, and (2) the gradual growth of nonlinear effects during wave propagation.

A model of wave behavior, which includes, in an approximate manner, both linear and nonlinear effects is presented for computing wave profiles after the linear theory has become invalid due to the growth of nonlinearities during wave propagation.

An experimental program has been conducted to confirm both the linear model for the two-dimensional fluid domain and the strategy suggested for determining wave profiles during propagation after the linear theory becomes invalid. The effect of a more general time-displacement history of the moving bed than those employed in the theoretical models is also investigated experimentally.

The linear theory is found to accurately approximate the wave behavior in the region of generation whenever the total displacement of the bed is much less than the water depth. Curves are developed and confirmed by the experiments which predict gross features of the lead wave propagating from the region of generation once the values of certain nondimensional parameters (which characterize the generation process) are known. For example, the maximum amplitude of the lead wave propagating from the region of generation has been found to never exceed approximately one-half of the total bed displacement. The gross features of the tsunami resulting from the Alaskan earthquake of 27 March 1964 can be estimated from the results of this study.

Cytokine pathway disruption in a mouse model of schizophrenia induced by Munc18-1a overexpression in the brain

Background: An accumulating body of evidence points to the significance of neuroinflammation and immunogenetics in schizophrenia, and an imbalance of cytokines in the central nervous system (CNS) has been suggested to be associated with the disorder. Munc18-overexpressing mice (Munc18-OE) have provided a model for the study of the alterations that may underlie the symptoms of subjects with schizophrenia. The aim of the present study was to elucidate the involvement of neuroinflammation and cytokine imbalance in this model. Methods: Cytokines were evaluated in the cortex and the striatum of Munc18-OE and wild-type (WT) mice by enzyme-linked immunosorbent assay (ELISA). Protein levels of specific microglia and macrophage, astrocytic and neuroinflammation markers were quantified by western blot in the cortex and the striatum of Munc18-OE and WT mice. Results: Each cytokine evaluated (Interferon-gamma (IFN-gamma), Tumor Necrosis Factor-alpha (TNF-alpha), Interleukin-2 (IL-2) and CCL2 chemokine) was present at higher levels in the striatum of Munc18-OE mice than WT. Cortical TNF-alpha and IL-2 levels were significantly lower in Munc18-OE mice than WT mice. The microglia and macrophage marker CD11b was lower in the cortexes of Munc18-OE mice than WT, but no differences were observed in the striatum. Glial Fibrillary Acidic Protein (GFAP) and Nuclear Factor-kappaB (NF-kappa B)p65 levels were not different between the groups. Interleukin-1beta (IL-1 beta) and IL-6 levels were beneath detection limits. Conclusions: The disrupted levels of cytokines detected in the brain of Munc18-OE mice was found to be similar to clinical reports and endorses study of this type for analysis of this aspect of the disorder. The lower CD11b expression in the cortex but not in the striatum of the Munc18-OE mice may reflect differences in physiological activity. The cytokine expression pattern observed in Munc18-OE mice is similar to a previously published model of schizophrenia caused by maternal immune activation. Together, these data suggest a possible role for an immune imbalance in this disorder.

Brca1 is expressed in human microglia and is dysregulated in human and animal model of ALS

Background: There is growing evidence that microglia are key players in the pathological process of amyotrophic lateral sclerosis (ALS). It is suggested that microglia have a dual role in motoneurone degeneration through the release of both neuroprotective and neurotoxic factors. Results: To identify candidate genes that may be involved in ALS pathology we have analysed at early symptomatic age (P90), the molecular signature of microglia from the lumbar region of the spinal cord of hSOD1(G93A) mice, the most widely used animal model of ALS. We first identified unique hSOD1(G93A) microglia transcriptomic profile that, in addition to more classical processes such as chemotaxis and immune response, pointed toward the potential involvement of the tumour suppressor gene breast cancer susceptibility gene 1 (Brca1). Secondly, comparison with our previous data on hSOD1(G93A) motoneurone gene profile substantiated the putative contribution of Brca1 in ALS. Finally, we established that Brca1 protein is specifically expressed in human spinal microglia and is up-regulated in ALS patients. Conclusions: Overall, our data provide new insights into the pathogenic concept of a non-cell-autonomous disease and the involvement of microglia in ALS. Importantly, the identification of Brca1 as a novel microglial marker and as possible contributor in both human and animal model of ALS may represent a valid therapeutic target. Moreover, our data points toward novel research strategies such as investigating the role of oncogenic proteins in neurodegenerative diseases.

Model-based control of a Rijke tube combustion instability

A Rijke tube is used to demonstrate model-based control of a combustion instability, where controller design is based on measurement of the unstable system. The Rijke tube used was of length 0.75m and had a grid-stabilised laminar flame in its lower half. A microphone was used as a sensor and a loudspeaker as an actuator for active control. The open loop transfer function (OLTF) required for controller design was that from the actuator to the sensor. This was measured experimentally by sending a signal with two components to the actuator. The first was a control component from an empirically designed controller, which was used to stabilise the system, thus eliminating the non-linear limit cycle. The second was a high bandwidth signal for identification of the OLTF. This approach to measuring the OLTF is generic and can be applied to large-scale combustors. The measured OLTF showed that only the fundamental mode of the tube was unstable; this was consistent with the OLTF predicted by a mathematical model of the tube, involving 1-D linear acoustic waves and a time delay heat release model. Based on the measured OLTF, a controller to stabilise the instability was designed using Nyquist techniques. This was implemented and was seen to result in an 80dB reduction in the microphone pressure spectrum. A robustness study was performed by adding an additional length to the top of the Rijke tobe. The controller was found to achieve control up to an increase in tube length of 19%. This compared favourably with the empirical controller, which lost control for an increase in tube length of less than 3%.

Angular dependence of domain wall resistivity in artificial magnetic domain structures.

We exploit the ability to precisely control the magnetic domain structure of perpendicularly magnetized Pt/Co/Pt trilayers to fabricate artificial domain wall arrays and study their transport properties. The scaling behavior of this model system confirms the intrinsic domain wall origin of the magnetoresistance, and systematic studies using domains patterned at various angles to the current flow are excellently described by an angular-dependent resistivity tensor containing perpendicular and parallel domain wall resistivities. We find that the latter are fully consistent with Levy-Zhang theory, which allows us to estimate the ratio of minority to majority spin carrier resistivities, rho downward arrow/rho upward arrow approximately 5.5, in good agreement with thin film band structure calculations.