Methods in general model localization

The aim of this study was to evaluate and test methods which could improve local estimates of a general model fitted to a large area. In the first three studies, the intention was to divide the study area into sub-areas that were as homogeneous as possible according to the residuals of the general model, and in the fourth study, the localization was based on the local neighbourhood. According to spatial autocorrelation (SA), points closer together in space are more likely to be similar than those that are farther apart. Local indicators of SA (LISAs) test the similarity of data clusters. A LISA was calculated for every observation in the dataset, and together with the spatial position and residual of the global model, the data were segmented using two different methods: classification and regression trees (CART) and the multiresolution segmentation algorithm (MS) of the eCognition software. The general model was then re-fitted (localized) to the formed sub-areas. In kriging, the SA is modelled with a variogram, and the spatial correlation is a function of the distance (and direction) between the observation and the point of calculation. A general trend is corrected with the residual information of the neighbourhood, whose size is controlled by the number of the nearest neighbours. Nearness is measured as Euclidian distance. With all methods, the root mean square errors (RMSEs) were lower, but with the methods that segmented the study area, the deviance in single localized RMSEs was wide. Therefore, an element capable of controlling the division or localization should be included in the segmentation-localization process. Kriging, on the other hand, provided stable estimates when the number of neighbours was sufficient (over 30), thus offering the best potential for further studies. Even CART could be combined with kriging or non-parametric methods, such as most similar neighbours (MSN).

Correlation of Oxygen Storage Capacity and Structural Distortion in Transition-Metal-, Noble-Metal-, and Rare-Earth-Ion-Substituted CeO2 from First Principles Calculation

Oxygen storage/release (OSC) capacity is an important feature common to all three-way catalysts to combat harmful exhaust emissions. To understand the mechanism of improved OSC for doped CeO2, we undertook the structural investigation by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), H-2-TPR (temperature-programmed hydrogen reduction) and density functional theoretical (DFT) calculations of transition-metal-, noble-metal-, and rare-earth (RE)-ion-substituted ceria. In this report, we present the relationship between the OSC and structural changes induced by the dopant ion in CeO2. Transition metal and noble metal ion substitution in ceria greatly enhances the reducibility of Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu, Pd, Pt, Ru), whereas rare-earth-ion-substituted Ce(1-x)A(x)O(2-delta) (A = La, Y) have very little effect in improving the OSC. Our simulated optimized structure shows deviation in cation oxygen bond length from ideal bond length of 2.34 angstrom (for CeO2). For example, our theoretical calculation for Ce28Mn4O62 structure shows that Mn-O bonds are in 4 + 2 coordination with average bond lengths of 2.0 and 3.06 angstrom respectively. Although the four short Mn-O bond lengths spans the bond distance region of Mn2O3, the other two Mn-O bonds are moved to longer distances. The dopant transition and noble metal ions also affects Ce coordination shell and results in the formation of longer Ce-O bonds as well. Thus longer cation oxygen bonds for both dopant and host ions results in enhanced synergistic reduction of the solid solution. With Pd ion substitution in Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu) further enhancement in OSC is observed in H-2-TPR. This effect is reflected in our model calculations by the presence of still longer bonds compared to the model without Pd ion doping. The synergistic effect is therefore due to enhanced reducibility of both dopant and host ion induced due to structural distortion of fluorite lattice in presence of dopant ion. For RE ions (RE = Y, La), our calculations show very little deviation of bonds lengths from ideal fluorite structure. The absence of longer Y-O/La-O and Ce-O bonds make the structure much less susceptible to reduction.

Identification of collagenolytic MMP networks as potential biomarkers in progressive chronic periodontitis subjects and MMP-8 null allele model

Chronic periodontitis results from a complex aetiology, including the formation of a subgingival biofilm and the elicitation of the host s immune and inflammatory response. The hallmark of chronic periodontitis is alveolar bone loss and soft periodontal tissue destruction. Evidence supports that periodontitis progresses in dynamic states of exacerbation and remission or quiescence. The major clinical approach to identify disease progression is the tolerance method, based on sequential probing. Collagen degradation is one of the key events in periodontal destructive lesions. Matrix metalloproteinase (MMP)-8 and MMP-13 are the primary collagenolytic MMPs that are associated with the severity of periodontal inflammation and disease, either by a direct breakdown of the collagenised matrix or by the processing of non-matrix bioactive substrates. Despite the numerous host mediators that have been proposed as potential biomarkers for chronic periodontitis, they reflect inflammation rather than the loss of periodontal attachment. The aim of the present study was to determine the key molecular MMP-8 and -13 interactions in gingival crevicular fluid (GCF) and gingival tissue from progressive periodontitis lesions and MMP-8 null allele mouse model. In study (I), GCF and gingival biopsies from active and inactive sites of chronic periodontitis patients, which were determined clinically by the tolerance method, and healthy GCF were analysed for MMP-13 and tissue inhibitor of matrix metalloproteinases (TIMP)-1. Chronic periodontitis was characterised by increased MMP-13 levels and the active sites showed a tendency of decreased TIMP-1 levels associated with increments of MMP-13 and total protein concentration compared to inactive sites. In study (II), we investigated whether MMP-13 activity was associated with TIMP-1, bone collagen breakdown through ICTP levels, as well as the activation rate of MMP-9 in destructive lesions. The active sites demonstrated increased GCF ICTP levels as well as lowered TIMP-1 detection along with elevated MMP-13 activity. MMP-9 activation rate was enhanced by MMP-13 in diseased gingival tissue. In study (III), we analysed the potential association between the levels, molecular forms, isoenzyme distribution and degree of activation of MMP-8, MMP-14, MPO and the inhibitor TIMP-1 in GCF from periodontitis progressive patients at baseline and after periodontal therapy. A positive correlation was found for MPO/MMP-8 and their levels associated with progression episodes and treatment response. Because MMP-8 is activated by hypochlorous acid in vitro, our results suggested an interaction between the MPO oxidative pathway and MMP-8 activation in GCF. Finally, in study (IV), on the basis of the previous finding that MMP-8-deficient mice showed impaired neutrophil responses and severe alveolar bone loss, we aimed to characterise the detection patterns of LIX/CXCL5, SDF-1/CXCL12 and RANKL in P. gingivalis-induced experimental periodontitis and in the MMP-8-/- murine model. The detection of neutrophil-chemoattractant LIX/CXCL5 was restricted to the oral-periodontal interface and its levels were reduced in infected MMP-8 null mice vs. wild type mice, whereas the detection of SDF-1/CXCL12 and RANKL in periodontal tissues increased in experimentally-induced periodontitis, irrespectively from the genotype. Accordingly, MMP-8 might regulate LIX/CXCL5 levels by undetermined mechanisms, and SDF-1/CXCL12 and RANKL might promote the development and/or progression of periodontitis.

Exact eigenstates of the Pariser-Parr-Pople model for anthracene

The pi-electronic structure of anthracene is discussed by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical PM3 calculations. Symmetry adaptation of the 2.8 million singlet valence-bond (VB) diagrams is explicitly demonstrated for D2h and electron-hole symmetry. Standard PPP parameters provide a comprehensive fit to one- and two-photon anthracene spectra and intensities up to the strong 1 B-1(3u)-absorption at 5.24 eV, the 10th excited state in the dense correlated spectrum, and indicate a reassignment of two-photon absorptions. The singlet-triplet gap and fine-structure constants also agree with experiment. Fully-relaxed PM3 geometries are obtained for the anthracene ground state and for singlet, triplet, and charged bipolarons. The PM3 bond lengths correlate well with PPP bond orders for the idealized structure. Single-determinantal PM3 excitation and relaxation energies for bipolarons are consistent with exact PPP results and contrast all-valence electron with pi-electron calculations. Several correlation effects are noted in the rich pi-spectra of anthracene in connection with improved PPP modeling of conjugated molecules and polymers.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Correlation of seasonal changes in sperm output with endocrinological changes in the adult male bonnet monkey, Macaca radiata

We have examined the monthly variations in sperm output and attempted to correlate the profiles of endocrine hormones secreted with the sperm counts throughout the ,year in the adult male bonnet monkey. As previously reported, there was a distinct spurt in sperm output beginning September through December months. A concomitant increase in serum testosterone and prolactin concentrations were also noted during September through November (mid and post-monsoon season). Although there was a marked increase in gonadotropin releasing hormone stimulated testosterone secretion, the peak testosterone concentrations post gonadotropin releasing hormone injection did not vary significantly (P>0.05) throughout the year. Basal serum follicle stimulating hormone concentrations did not vary significantly (P>0.05) during April to June months compared to September-November months. Serum inhibin concentration remained unaltered throughout the year, except in the month of March. The results of this study provide evidence for annual rhythms in prolactin and testosterone secretion and a distinct seasonality in the sperm output of the adult male bonnet monkey, but the pituitary responsiveness to exogenous gonadotropin releasing hormone remains unaltered throughout the year. Because of the existence of seasonality as noted in the present study, future studies which utilize the adult male bonnet monkey as an experimental model need to take into consideration the seasonal effects on reproductive function in this species.

Molecular Dynamics Investigation of Structure and Transport in the K2O-2SiO2 System Using a Partial Charge Based Model Potential

Potassium disilicate glass and melt have been investigated by using a new partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygenrich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.

Molecular dynamics Investigation of structure and transport in the K2O-2SiO2 system using a partial charge based model potential

Potassium disilicate glass and melt have been investigated by using anew partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygen-rich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.

Perturbative growth of cosmological clustering .2. The two-point correlation

We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest order nonlinear correction to the two-point correlation and the pair velocity for Gaussian initial conditions in a critical density matter-dominated cosmological model. We compare our results with the results obtained using the hydrodynamic equations that neglect pressure and find that the two match, indicating that there are no effects of multistreaming at this order of perturbation. We analytically study the effect of small scales on the large scales by calculating the nonlinear correction for a Dirac delta function initial two-point correlation. We find that the induced two-point correlation has a x(-6) behavior at large separations. We have considered a class of initial conditions where the initial power spectrum at small k has the form k(n) with 0 < n less than or equal to 3 and have numerically calculated the nonlinear correction to the two-point correlation, its average over a sphere and the pair velocity over a large dynamical range. We find that at small separations the effect of the nonlinear term is to enhance the clustering, whereas at intermediate scales it can act to either increase or decrease the clustering. At large scales we find a simple formula that gives a very good fit for the nonlinear correction in terms of the initial function. This formula explicitly exhibits the influence of small scales on large scales and because of this coupling the perturbative treatment breaks down at large scales much before one would expect it to if the nonlinearity were local in real space. We physically interpret this formula in terms of a simple diffusion process. We have also investigated the case n = 0, and we find that it differs from the other cases in certain respects. We investigate a recently proposed scaling property of gravitational clustering, and we find that the lowest order nonlinear terms cause deviations from the scaling relations that are strictly valid in the linear regime. The approximate validity of these relations in the nonlinear regime in l(T)-body simulations cannot be understood at this order of evolution.

Sponge phase transitions from a lattice model

We use Monte Carlo simulations to obtain thermodynamic functions and correlation functions in a lattice model we propose for sponge phases. We demonstrate that the surface-density correlation function dominates the scattering only along the symmetric-sponge (SS) to asymmetric-sponge (AS) phase boundary but not the boundary between the sponge-with-free-edges (SFE) and symmetric-sponge phases. At this second thermodynamic transition the scattering is dominated instead by an edge-density (or seam-density) correlation function. This prediction provides an unambiguous diagnostic for experiments in search of the SS-SFE transition.

Simulation of ENSO-Related Surface Winds in the Tropical Pacific by an Atmospheric General Circulation Model Forced by Observed Sea Surface Temperatures

The authors present the simulation of the tropical Pacific surface wind variability by a low-resolution (R15 horizontal resolution and 18 vertical levels) version of the Center for Ocean-Land-Atmosphere Interactions, Maryland, general circulation model (GCM) when forced by observed global sea surface temperature. The authors have examined the monthly mean surface winds acid precipitation simulated by the model that was integrated from January 1979 to March 1992. Analyses of the climatological annual cycle and interannual variability over the Pacific are presented. The annual means of the simulated zonal and meridional winds agree well with observations. The only appreciable difference is in the region of strong trade winds where the simulated zonal winds are about 15%-20% weaker than observed, The amplitude of the annual harmonics are weaker than observed over the intertropical convergence zone and the South Pacific convergence zone regions. The amplitudes of the interannual variation of the simulated zonal and meridional winds are close to those of the observed variation. The first few dominant empirical orthogonal functions (EOF) of the simulated, as well as the observed, monthly mean winds are found to contain a targe amount of high-frequency intraseasonal variations, While the statistical properties of the high-frequency modes, such as their amplitude and geographical locations, agree with observations, their detailed time evolution does not. When the data are subjected to a 5-month running-mean filter, the first two dominant EOFs of the simulated winds representing the low-frequency EI Nino-Southern Oscillation fluctuations compare quite well with observations. However, the location of the center of the westerly anomalies associated with the warm episodes is simulated about 15 degrees west of the observed locations. The model simulates well the progress of the westerly anomalies toward the eastern Pacific during the evolution of a warm event. The simulated equatorial wind anomalies are comparable in magnitude to the observed anomalies. An intercomparison of the simulation of the interannual variability by a few other GCMs with comparable resolution is also presented. The success in simulation of the large-scale low-frequency part of the tropical surface winds by the atmospheric GCM seems to be related to the model's ability to simulate the large-scale low-frequency part of the precipitation. Good correspondence between the simulated precipitation and the highly reflective cloud anomalies is seen in the first two EOFs of the 5-month running means. Moreover, the strong correlation found between the simulated precipitation and the simulated winds in the first two principal components indicates the primary role of model precipitation in driving the surface winds. The surface winds simulated by a linear model forced by the GCM-simulated precipitation show good resemblance to the GCM-simulated winds in the equatorial region. This result supports the recent findings that the large-scale part of the tropical surface winds is primarily linear.

Density-matrix renormalization-group study of low-lying excitations of polyacene within a Pariser-Parr-Pople model

We have carried out symmetrized density-matrix renormalization-group calculations to study the nature of excited states of long polyacene oligomers within a Pariser-Parr-Pople Hamiltonian. We have used the C-2 symmetry, the electron-hole symmetry, and the spin parity of the system in our calculations. We find that there is a crossover in the lowest dipole forbidden two-photon state and the lowest dipole allowed excited state with size of the oligomer. In the long system limit, the two-photon state lies below the lowest dipole allowed excited state. The triplet state lies well below the two-photon state and energetically does not correspond to its description as being made up of two triplets. These results are in agreement with the general trends in linear conjugated polymers. However, unlike in linear polyenes wherein the two-photon state is a localized excitation, we find that in polyacenes, the two-photon excitation is spread out over the system. We have doped the systems with a hole and an electron and have calculated the charge excitation gap. Using the charge gap and the optical gap, we estimate the binding energy of the 1(1)B(-) exciton to be 2.09 eV. We have also studied doubly doped polyacenes and find that the bipolaron in these systems, to be composed of two separated polarons, as indicated by the calculated charge-density profile and charge-charge correlation function. We have studied bond orders in various states in order to get an idea of the excited state geometry of the system. We find that the ground state, the triplet state, the dipole allowed state, and the polaron excitations correspond to lengthening of the rung bonds in the interior of the oligomer while the two-photon excitation corresponds to the rung bond lengths having two maxima in the system.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 2. The fluctuating-force model

A fluctuating-force model is developed for representing the effect of the turbulent fluid velocity fluctuations on the particle phase in a turbulent gas–solid suspension in the limit of high Stokes number, where the particle relaxation time is large compared with the correlation time for the fluid velocity fluctuations. In the model, a fluctuating force is incorporated in the equation of motion for the particles, and the force distribution is assumed to be an anisotropic Gaussian white noise. It is shown that this is equivalent to incorporating a diffusion term in the Boltzmann equation for the particle velocity distribution functions. The variance of the force distribution, or equivalently the diffusion coefficient in the Boltzmann equation, is related to the time correlation functions for the fluid velocity fluctuations. The fluctuating-force model is applied to the specific case of a Couette flow of a turbulent particle–gas suspension, for which both the fluid and particle velocity distributions were evaluated using direct numerical simulations by Goswami & Kumaran (2010). It is found that the fluctuating-force simulation is able to quantitatively predict the concentration, mean velocity profiles and the mean square velocities, both at relatively low volume fractions, where the viscous relaxation time is small compared with the time between collisions, and at higher volume fractions, where the time between collisions is small compared with the viscous relaxation time. The simulations are also able to predict the velocity distributions in the centre of the Couette, even in cases in which the velocity distribution is very different from a Gaussian distribution.

NETGEM: Network Embedded Temporal GEnerative Model for gene expression data

Background: Temporal analysis of gene expression data has been limited to identifying genes whose expression varies with time and/or correlation between genes that have similar temporal profiles. Often, the methods do not consider the underlying network constraints that connect the genes. It is becoming increasingly evident that interactions change substantially with time. Thus far, there is no systematic method to relate the temporal changes in gene expression to the dynamics of interactions between them. Information on interaction dynamics would open up possibilities for discovering new mechanisms of regulation by providing valuable insight into identifying time-sensitive interactions as well as permit studies on the effect of a genetic perturbation. Results: We present NETGEM, a tractable model rooted in Markov dynamics, for analyzing the dynamics of the interactions between proteins based on the dynamics of the expression changes of the genes that encode them. The model treats the interaction strengths as random variables which are modulated by suitable priors. This approach is necessitated by the extremely small sample size of the datasets, relative to the number of interactions. The model is amenable to a linear time algorithm for efficient inference. Using temporal gene expression data, NETGEM was successful in identifying (i) temporal interactions and determining their strength, (ii) functional categories of the actively interacting partners and (iii) dynamics of interactions in perturbed networks. Conclusions: NETGEM represents an optimal trade-off between model complexity and data requirement. It was able to deduce actively interacting genes and functional categories from temporal gene expression data. It permits inference by incorporating the information available in perturbed networks. Given that the inputs to NETGEM are only the network and the temporal variation of the nodes, this algorithm promises to have widespread applications, beyond biological systems. The source code for NETGEM is available from https://github.com/vjethava/NETGEM

Confinement-deconfinement transition and spin correlations in a generalized Kitaev model

We present a spin model, namely, the Kitaev model augmented by a loop term and perturbed by an Ising Hamiltonian, and show that it exhibits both confinement-deconfinement transitions from spin liquid to antiferromagnetic/spin-chain/ferromagnetic phases and topological quantum phase transitions between gapped and gapless spin-liquid phases. We develop a fermionic resonating-valence-bonds (RVB) mean-field theory to chart out the phase diagram of the model and estimate the stability of its spin-liquid phases, which might be relevant for attempts to realize the model in optical lattices and other spin systems. We present an analytical mean-field theory to study the confinement-deconfinement transition for large coefficient of the loop term and show that this transition is first order within such mean-field analysis in this limit. We also conjecture that in some other regimes, the confinement-deconfinement transitions in the model, predicted to be first order within the mean-field theory, may become second order via a defect condensation mechanism. Finally, we present a general classification of the perturbations to the Kitaev model on the basis of their effect on it's spin correlation functions and derive a necessary and sufficient condition, within the regime of validity of perturbation theory, for the spin correlators to exhibit a long-ranged power-law behavior in the presence of such perturbations. Our results reproduce those of Tikhonov et al. [Phys. Rev. Lett. 106, 067203 (2011)] as a special case.