On the relation between K-means and PLSA

Non-negative matrix factorization [5](NMF) is a well known tool for unsupervised machine learning. It can be viewed as a generalization of the K-means clustering, Expectation Maximization based clustering and aspect modeling by Probabilistic Latent Semantic Analysis (PLSA). Specifically PLSA is related to NMF with KL-divergence objective function. Further it is shown that K-means clustering is a special case of NMF with matrix L2 norm based error function. In this paper our objective is to analyze the relation between K-means clustering and PLSA by examining the KL-divergence function and matrix L2 norm based error function.

The evolution of complexity in social organization-A model using dominance-subordinate behavior in two social wasp species

Dominance and subordinate behaviors are important ingredients in the social organizations of group living animals. Behavioral observations on the two eusocial species Ropalidia marginata and Ropalidia cyathiformis suggest varying complexities in their social systems. The queen of R. cyathiformis is an aggressive individual who usually holds the top position in the dominance hierarchy although she does not necessarily show the maximum number of acts of dominance, while the R. marginata queen rarely shows aggression and usually does not hold the top position in the dominance hierarchy of her colony. In R. marginata, more workers are involved in dominance-subordinate interactions as compared to R. cyathiformis. These differences are reflected in the distribution of dominance-subordinate interactions among the hierarchically ranked individuals in both the species. The percentage of dominance interactions decreases gradually with hierarchical ranks in R. marginata while in R. cyathiformis it first increases and then decreases. We use an agent-based model to investigate the underlying mechanism that could give rise to the observed patterns for both the species. The model assumes, besides some non-interacting individuals, the interaction probabilities of the agents depend on their pre-differentiated winning abilities. Our simulations show that if the queen takes up a strategy of being involved in a moderate number of dominance interactions, one could get the pattern similar to R. cyathiformis, while taking up the strategy of very low interactions by the queen could lead to the pattern of R. marginata. We infer that both the species follow a common interaction pattern, while the differences in their social organization are due to the slight changes in queen as well as worker strategies. These changes in strategies are expected to accompany the evolution of more complex societies from simpler ones.

Quick recovery in transactional flash

Maintaining metadata consistency is a critical issue in designing a filesystem. Although satisfactory solutions are available for filesystems residing on magnetic disks, these solutions may not give adequate performance for filesystems residing on flash devices. Prabhakaran et al. have designed a metadata consistency mechanism specifically for flash chips, called Transactional Flash1]. It uses cyclic commit mechanism to provide transactional abstractions. Although significant improvement over usual journaling techniques, this mechanism has certain drawbacks such as complex protocol and necessity to read whole flash during recovery, which slows down recovery process. In this paper we propose addition of thin journaling layer on top of Transactional Flash to simplify the protocol and speed up the recovery process. The simplified protocol named Quick Recovery Cyclic Commit (QRCC) uses journal stored on NOR flash for recovery. Our evaluations on actual raw flash card show that journal writes add negligible penalty compared to original Transactional Flash's write performance, while quick recovery is facilitated by journal in case of failures.

Warped alternatives to Froggatt-Nielsen models

We consider the Randall-Sundrum (RS) setup to be a theory of flavor, as an alternative to Froggatt-Nielsen models instead of as a solution to the hierarchy problem. The RS framework is modified by taking the low-energy brane to be at the grand unified theory (GUT) scale. This also alleviates constraints from flavor physics. Fermion masses and mixing angles are fit at the GUT scale. The ranges of the bulk mass parameters are determined using a chi(2) fit taking into consideration the variation in O(1) parameters. In the hadronic sector, the heavy top quark requires large bulk mass parameters localizing the right-handed top quark close to the IR brane. Two cases of neutrino masses are considered: (a) Planck scale lepton number violation and (b) Dirac neutrino masses. Contrary to the case of weak scale RS models, both these cases give reasonable fits to the data, with the Planck scale lepton number violation fitting slightly better compared to the Dirac case. In the supersymmetric version, the fits are not significantly different except for the variation in tan beta. If the Higgs superfields and the supersymmetry breaking spurion are localized on the same brane, then the structure of the sfermion masses are determined by the profiles of the zero modes of the hypermultiplets in the bulk. Trilinear terms have the same structure as the Yukawa matrices. The resultant squark spectrum is around similar to 2-3 TeV required by the light Higgs mass to be around 125 GeV and to satisfy the flavor violating constraints.

Physical Layer Network Coding for the K-User Multiple Access Relay Channel

We propose a Physical layer Network Coding (PNC) scheme for the K-user wireless Multiple Access Relay Channel, in which K source nodes want to transmit messages to a destination node D with the help of a relay node R. The proposed scheme involves (i) Phase 1 during which the source nodes alone transmit and (ii) Phase 2 during which the source nodes and the relay node transmit. At the end of Phase 1, the relay node decodes the messages of the source nodes and during Phase 2 transmits a many-to-one function of the decoded messages. To counter the error propagation from the relay node, we propose a novel decoder which takes into account the possibility of error events at R. It is shown that if certain parameters are chosen properly and if the network coding map used at R forms a Latin Hypercube, the proposed decoder offers the maximum diversity order of two. Also, it is shown that for a proper choice of the parameters, the proposed decoder admits fast decoding, with the same decoding complexity order as that of the reference scheme based on Complex Field Network Coding (CFNC). Simulation results indicate that the proposed PNC scheme offers a large gain over the CFNC scheme.

Top polarization in stop production at the LHC

We survey the expected polarization of the top produced in the decay of a scalar top quark, t -> t((chi) over tildei)(0), i = 1 - 2. The phenomenology is quite interesting, since the expected polarization depends both on the mixing in the stop and neutralino sectors and on the mass differences between the stop and the neutralino. We find that a mixed stop behaves almost like a right-handed stop due to the larger hypercharge that enters the stop/top/gaugino coupling and that these polarisation effects disappear, when m((t) over tilde1) approximate to m(t) +m((chi) over tildei)(0). After a discussion on the expected top polarization from the decay of a scalar top quark, we focus on the interplay of polarization and kinematics at the LHC. We discuss different probes of the top polarization in terms of lab-frame observables. We find that these observables faithfully reflect the polarization of the parent top-quark, but also have a non-trivial dependence on the kinematics of the stop production and decay process. In addition, we illustrate the effect of top polarization on the energy and transverse momentum of the decay lepton in the laboratory frame. Our results show that both spectra are softened substantially in case of a negatively polarized top, particularly for a large mass difference between the stop and the neutralino. Thus, the search strategies, and the conclusions that can be drawn from them, depends not just on the mass difference m((t) over tilde) - m((chi) over tildei)(0) due to the usual kinematic effects but also on the effects of top polarization on the decay kinematics the extent of which depends in turn on the said mass difference.

Inference for the Component and System Lifetime Distribution of a k-unit Parallel System Based on System Data

In this paper, we consider the inference for the component and system lifetime distribution of a k-unit parallel system with independent components based on system data. The components are assumed to have identical Weibull distribution. We obtain the maximum likelihood estimates of the unknown parameters based on system data. The Fisher information matrix has been derived. We propose -expectation tolerance interval and -content -level tolerance interval for the life distribution of the system. Performance of the estimators and tolerance intervals is investigated via simulation study. A simulated dataset is analyzed for illustration.

On k-stellated and k-stacked spheres

We introduce the class Sigma(k)(d) of k-stellated (combinatorial) spheres of dimension d (0 <= k <= d + 1) and compare and contrast it with the class S-k(d) (0 <= k <= d) of k-stacked homology d-spheres. We have E-1(d) = S-1(d), and Sigma(k)(d) subset of S-k(d) ford >= 2k-1. However, for each k >= 2 there are k-stacked spheres which are not k-stellated. For d <= 2k - 2, the existence of k-stellated spheres which are not k-stacked remains an open question. We also consider the class W-k(d) (and K-k(d)) of simplicial complexes all whose vertex-links belong to Sigma(k)(d - 1) (respectively, S-k(d - 1)). Thus, W-k(d) subset of K-k(d) for d >= 2k, while W-1(d) = K-1(d). Let (K) over bar (k)(d) denote the class of d-dimensional complexes all whose vertex-links are k-stacked balls. We show that for d >= 2k + 2, there is a natural bijection M -> (M) over bar from K-k(d) onto (K) over bar (k)(d + 1) which is the inverse to the boundary map partial derivative: (K) over bar (k)(d + 1) -> (K) over bar (k)(d). Finally, we complement the tightness results of our recent paper, Bagchi and Datta (2013) 5], by showing that, for any field F, an F-orientable (k + 1)-neighbourly member of W-k(2k + 1) is F-tight if and only if it is k-stacked.

Formation energies and the stability of the oxides of K

The most common valencies associated with K and O atoms are 1+ and 2-. As a result, one expects K2O to be the oxide of potassium which is the most stable with respect to its constituents. Calculating the formation energy within electronic structure calculations using hybrid functionals, one finds that K2O2 has the largest formation energy, implying the largest stability of this oxide of potassium with respect to its constituents. This is traced to the presence of oxygen dimers in the K2O2 structure which interact strongly resulting in a larger formation energy compared to the more ionic K2O.

Bipartite Powers of k-chordal Graphs

Let k be an integer and k >= 3. A graph G is k-chordal if G does not have an induced cycle of length greater than k. From the definition it is clear that 3-chordal graphs are precisely the class of chordal graphs. Duchet proved that, for every positive integer m, if G m is chordal then so is G(m+2). Brandst `` adt et al. in Andreas Brandsadt, Van Bang Le, and Thomas Szymczak. Duchet- type theorems for powers of HHD- free graphs. Discrete Mathematics, 177(1- 3): 9- 16, 1997.] showed that if G m is k - chordal, then so is G(m+2). Powering a bipartite graph does not preserve its bipartitedness. In order to preserve the bipartitedness of a bipartite graph while powering Chandran et al. introduced the notion of bipartite powering. This notion was introduced to aid their study of boxicity of chordal bipartite graphs. The m - th bipartite power G(m]) of a bipartite graph G is the bipartite graph obtained from G by adding edges (u; v) where d G (u; v) is odd and less than or equal to m. Note that G(m]) = G(m+1]) for each odd m. In this paper we show that, given a bipartite graph G, if G is k-chordal then so is G m], where k, m are positive integers with k >= 4

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.

Temperature dependent magnetic ordering and electrical transport behavior of nano zinc ferrite from 20 to 800 K

The nano ZnFe2O4 compound was prepared by eco-friendly hydrothermal method. The characterization of the sample for its structure, morphology and composition were done by powder X-ray diffraction (PXRD), scanning electron microscopy (SEM), dynamic light scattering, Fourier transform infrared spectroscopy, zeta surface profiler and UV-Visible spectroscopy studies. The PXRD measurement reveals that the compound shows spinel cubic phase belong Fd (3) over barm (227) space group. Morphology of the compound from SEM and surface profile shows nearly spherical agglomerated particles with well defined grains and grain boundaries. The material shows the semiconducting behavior with E-g of 2.3 eV at room temperature (RT). The variation in the magnetic ordering was observed for wide range of temperature. The compound behaves like a soft magnetic material with ferrimagnetic at various temperatures except at RT. Both magnetic and EPR studies supports the superparamagnetic behavior of the the sample. The DC conductivity, dielectric and AC conductivity behavior of the 1000 degrees C pellets sintered for 2 h shows good frequency dependent transport properties. The present study facilitate in selecting the suitable materials for the nanoelectronics and spintronic applications. (C) 2013 Elsevier B.V. All rights reserved.

Studies on rheocasting using cooling slope

In the present work, a cooling channel is employed to produce semi-solid A356 alloy slurry. To understand the transport process involved, a 3D non-isothermal, multiphase volume averaging model has been developed for simulation of the semi-solid slurry generation process in the cooling channel. For simulation purpose, the three phases considered are the parent melt, the nearly spherical grains and air as separated but highly coupled interpenetrating continua. The conservation equations of mass, momentum, energy and species have been solved for each phase and the thermal and mechanical interactions (drag force) among the phases have been considered using appropriate model. The superheated liquid alloy is poured at the top of the cooling slope/channel, where specified velocity inlet boundary condition is used in the model, and allowed to flow along gravity through the channel. The melt loses its superheat and becomes semisolid up to the end of cooling channel due to the evolving -Al grains with decreasing temperature. The air phase forms a definable air/liquid melt interface, i.e. free surface, due its low density. The results obtained from the present model includes volume fractions of three different phases considered, grain evolution, grain growth rate, size and distribution of solid grains. The effect of key process variables such as pouring temperature, slope angle of the cooling channel and cooling channel wall temperature on temperature distribution, velocity distribution, grain formation and volume fraction of different phases are also studied. The results obtained from the simulations are validated by microstructure study using SEM and quantitative image analysis of the semi-solid slurry microstructure obtained from the experimental set-up.

Electronic structure of Nd1-xYxMnO3 from Mn K edge absorption spectroscopy and DFT methods

The electronic structure of Nd1-xYxMnO3 (x-0-0.5) is studied using x-ray absorption near-edge structure (XANES) spectroscopy at the Mn K-edge along with the DFT-based LSDA+U and real space cluster calculations. The main edge of the spectra does not show any variation with doping. The pre-edge shows two distinct features which appear well-separated with doping. The intensity of the pre-edge decreases with doping. The theoretical XANES were calculated using real space multiple scattering methods which reproduces the entire experimental spectra at the main edge as well as the pre-edge. Density functional theory calculations are used to obtain the Mn 4p, Mn 3d and O 2p density of states. For x=0, the site-projected density of states at 1.7 eV above Fermi energy shows a singular peak of unoccupied e(g) (spin-up) states which is hybridized Mn 4p and O 2p states. For x=0.5, this feature develops at a higher energy and is highly delocalized and overlaps with the 3d spin-down states which changes the pre-edge intensity. The Mn 4p DOS for both compositions, show considerable difference between the individual p(x), p(y) and p(z)), states. For x=0.5, there is a considerable change in the 4p orbital polarization suggesting changes in the Jahn-Teller effect with doping. (C) 2013 Elsevier Ltd. All rights reserved.

Probing the flavor violating scalar top quark signal at the LHC

The Large Hadron Collider (LHC) has completed its run at 8 TeV with the experiments ATLAS and CMS having collected about 25 fb(-1) of data each. Discovery of a light Higgs boson coupled with lack of evidence for supersymmetry at the LHC so far, has motivated studies of supersymmetry in the context of naturalness with the principal focus being the third generation squarks. In this work, we analyze the prospects of the flavor violating decay mode (t) over tilde (1) -> c chi(0)(1) at 8 and 13 TeV center-of-mass energy at the LHC. This channel is also relevant in the dark matter context for the stop-coannihilation scenario, where the relic density depends on the mass difference between the lighter stop quark ((t) over tilde (1)) and the lightest neutralino (chi(0)(1)) states. This channel is extremely challenging to probe, especially for situations when the mass difference between the lighter stop quark and the lightest neutralino is small. Using certain kinematical properties of signal events we find that the level of backgrounds can be reduced substantially. We find that the prospect for this channel is limited due to the low production cross section for top squarks and limited luminosity at 8 TeV, but at the 13 TeV LHC with 100 fb(-1) luminosity, it is possible to probe top squarks with masses up to similar to 450 GeV. We also discuss how the sensitivity could be significantly improved by tagging charm jets.