title: Hot Outflows in Galaxy Clusters

abstract: The gas-phase metallicity distribution has been analyzed for the hot atmospheres of 29 galaxy clusters using {\it Chandra X-ray Observatory} observations. All host brightest cluster galaxies (BCGs) with X-ray cavity systems produced by radio AGN. We find high elemental abundances projected preferentially along the cavities of 16 clusters. The metal-rich plasma was apparently lifted out of the BCGs with the rising X-ray cavities (bubbles) to altitudes between twenty and several hundred kiloparsecs. A relationship between the maximum projected altitude of the uplifted gas (the "iron radius") and jet power is found with the form $R_{\rm Fe} \propto P_{\rm jet}^{0.45}$. The estimated outflow rates are typically tens of solar masses per year but exceed $100 ~\rm M_\odot ~yr^{-1}$ in the most powerful AGN. The outflow rates are 10% to 20% of the cooling rates, and thus alone are unable to offset a cooling inflow. Nevertheless, hot outflows effectively redistribute the cooling gas and may play a significant role at regulating star formation and AGN activity in BCGs and presumably in giant elliptical galaxies. The metallicity distribution overall can be complex, perhaps due to metal rich gas returning in circulation flows or being blown around in the hot atmospheres. Roughly 15% of the work done by the cavities is expended lifting the metal-enriched gas, implying their nuclear black holes have increased in mass by at least $\sim 10^7$ M$_\odot$ to $10^9$ M$_\odot$. Finally, we show that hot outflows can account for the broad, gas-phase metallicity distribution compared to the stellar light profiles of BCGs, and we consider a possible connection between hot outflows and cold molecular gas flows discovered in recent ALMA observations.


title: A note on (anti-)self dual quasi Yamabe gradient soliton

abstract: In this note we prove that a (anti-)self dual quasi Yamabe soliton with positive sectional curvature is rotationally symmetric. This generalizes a recent result of G. Huang and H. Li in dimension four. Whence, (anti-) self dual gradient Yamabe solitons with positive sectional curvature is rotationally symmetric. We also prove that half conformally flat gradient Yamabe soliton has a special warped product structure provided that the potential function has no critical point.


title: Lattice Boltzmann simulations of a two-dimensional Fermi gas at unitarity

abstract: We present fully nonlinear dissipative fluid dynamics simulations of a trapped two-dimensional Fermi gas at unitarity using a Lattice Boltzmann algorithm. We are able to simulate non-harmonic trapping potentials, temperature-dependent viscosities as well as a discretized version of the ballistic (non-interacting) behavior. Our approach lends itself to direct comparison with experimental data, opening up the possibility of a precision determination of transport coefficients in the unitary Fermi gas. Furthermore, we predict the presence of a non-hydrodynamic component in the quadrupole mode, which should be observable experimentally.


title: Fluctuations and the rapidity dependence of charged particles spectra in fixed centrality bins in p$A$ collisions

abstract: We argue that large fluctuations in the saturation momentum are necessary to explain the ATLAS and ALICE data on pA collisions measured at the LHC. Using a form for the distribution of fluctuations motivated by theoretical studies of the non-linear evolution equations for the Color Glass Condensate, we find a remarkably good agreement between theory and the measured distributions. If the saturation momentum fluctuates, we argue that the cross section for a proton probe should also fluctuate, consistent with previous observations. Finally we discuss these results and their possible consistency with what is known about multiplicity fluctuations in such collisions, and attempts to describe the ridge effect from properties of the initial state.


title: Cross-correlation cosmic shear with the SDSS and VLA FIRST surveys

abstract: We measure the cosmic shear power spectrum on large angular scales by cross-correlating the shapes of ~9 million galaxies measured in the optical SDSS survey with the shapes of ~2.7x10^5 radio galaxies measured by the overlapping VLA-FIRST survey. Our measurements span the multipole range 10 < l < 130, corresponding to angular scales 2deg < {\theta} < 20deg. On these scales, the shear maps from both surveys suffer from significant systematic effects that prohibit a measurement of the shear power spectrum from either survey alone. Conversely we demonstrate that a power spectrum measured by cross-correlating the two surveys is unbiased. We measure an E-mode power spectrum from the data that is inconsistent with zero signal at the 99\% confidence (~2.7{\sigma}) level. The odd-parity B-mode signal and the EB cross- correlation are both found to be consistent with zero (within 1{\sigma}). These constraints are obtained after a careful error analysis that accounts for uncertainties due to cosmic variance, random galaxy shape noise and shape measurement errors, as well as additional errors associated with the observed large-scale systematic effects in the two surveys. Our constraints are consistent with the expected signal in the concordance cosmological model assuming recent estimates of the cosmological parameters from the Planck satellite, and literature values for the median redshifts of the SDSS and FIRST galaxy populations. The cross-power spectrum approach described in this paper represents a powerful technique for mitigating shear systematics and will be ideal for extracting robust results, with the exquisite control of systematics required, from future cosmic shear surveys with the SKA, LSST, Euclid and WFIRST-AFTA.


title: The Phase Diagram of Scalar Field Theory on the Fuzzy Disc

abstract: Using a recently developed bootstrapping method, we compute the phase diagram of scalar field theory on the fuzzy disc with quartic even potential. We find three distinct phases with second and third order phase transitions between them. In particular, we find that the second order phase transition happens approximately at a fixed ratio of the two coupling constants defining the potential. We compute this ratio analytically in the limit of large coupling constants. Our results qualitatively agree with previously obtained numerical results.


title: Yang-Baxter operators need quantum entanglement to distinguish knots

abstract: Any solution to the Yang-Baxter equation yields a family of representations of braid groups. Under certain conditions, identified by Turaev, the appropriately normalized trace of these representations yields a link invariant. Any Yang-Baxter solution can be interpreted as a two-qudit quantum gate. Here we show that if this gate is non-entangling, then the resulting invariant of knots is trivial. We thus obtain a general connection between topological entanglement and quantum entanglement, as suggested by Kauffman et al.


title: Towards single face shortest vertex-disjoint paths in undirected planar graphs

abstract: Given $k$ pairs of terminals $\{(s_{1}, t_{1}), \ldots, (s_{k}, t_{k})\}$ in a graph $G$, the min-sum $k$ vertex-disjoint paths problem is to find a collection $\{Q_{1}, Q_{2}, \ldots, Q_{k}\}$ of vertex-disjoint paths with minimum total length, where $Q_{i}$ is an $s_i$-to-$t_i$ path between $s_i$ and $t_i$. We consider the problem in planar graphs, where little is known about computational tractability, even in restricted cases. Kobayashi and Sommer propose a polynomial-time algorithm for $k \le 3$ in undirected planar graphs assuming all terminals are adjacent to at most two faces. Colin de Verdiere and Schrijver give a polynomial-time algorithm when all the sources are on the boundary of one face and all the sinks are on the boundary of another face and ask about the existence of a polynomial-time algorithm provided all terminals are on a common face. We make progress toward Colin de Verdiere and Schrijver's open question by giving an $O(kn^5)$ time algorithm for undirected planar graphs when $\{(s_{1}, t_{1}), \ldots, (s_{k}, t_{k})\}$ are in counter-clockwise order on a common face.


title: High degrees of random recursive trees

abstract: For $n\ge 1$, let $T_n$ be a random recursive tree on the vertex set $[n]=\{1,\ldots,n\}$. Let $\mathrm{deg}_{T_n}(v)$ be the degree of vertex $v$ in $T_n$, that is, the number of children of $v$ in $T_n$. Devroye and Lu showed that the maximum degree $\Delta_n$ of $T_n$ satisfies $\Delta_n/\lfloor \log_2 n\rfloor \to 1$ almost surely; Goh and Schmutz showed distributional convergence of $\Delta_n - \lfloor \log_2 n \rfloor$ along suitable subsequences. In this work we show how a version of Kingman's coalescent can be used to access much finer properties of the degree distribution in $T_n$. For any $i\in \mathbb{Z}$, let $X_i^{(n)}=|\{v\in [n]: \mathrm{deg}_{T_n}(v)= \lfloor \log n\rfloor +i\}|$. Also, let $\mathcal{P}$ be a Poisson point process on $\mathbb{R}$ with rate function $\lambda(x)=2^{-x}\cdot \ln 2$. We show that, up to lattice effects, the vectors $(X_i^{(n)},\, i\in \mathbb{Z})$ converge weakly in distribution to $(\mathcal{P}[i,i+1),\, i\in \mathbb{Z})$. We also prove asymptotic normality of $X_i^{(n)}$ when $i=i(n) \to -\infty$ slowly, and obtain precise asymptotics for $\mathbb{P}(\Delta_n - \log_2 n > i)$ when $ i(n) \to \infty$ and $i(n)/\log n$ is not too large. Our results recover and extends the previous results on maximal and near-maximal degrees in random recursive trees.


title: On locating statistics in the world of finding out

abstract: This paper attempts to situate statistics in relation to qualitative research methods and other means of "finding out". It compares and contrasts aspects of qualitative research methods and statistical inquiry and attempts to answer the question of whether and how elements of qualitative research methods should be included in statistics teaching.


title: The Starburst-AGN connection: quenching the fire and feeding the beast

abstract: The merger of two spiral galaxies is believed to be one of the main channels for the production of elliptical and early-type galaxies. In the process, the system becomes an (ultra) luminous infrared galaxy, or (U)LIRG, that morphs to a quasar, to a K+A galaxy, and finally to an early-type galaxy. The time scales for this metamorphosis are only loosely constrained by observations. In particular, the K+A phase should follow immediately after the QSO phase during which the dust and gas remaining from the (U)LIRG phase are expelled by the AGN. An intermediate class of QSOs with K+A spectral signatures, the post-starburst QSOs or PSQ, may represent the transitional phase between QSOs and K+As. We have compiled a sample of 72 {bona fide} $z<0.5$ PSQ from the SDSS DR7 QSO catalogue. We find the intermediate age populations in this sample to be on average significantly weaker and metal poorer than their putative descendants, the K+A galaxies. The typical spectral energy distribution of PSQ is well fitted by three components: starlight; an obscured power-law; and a hot dust component required to reproduce the mid-IR fluxes. From the slope and bolometric luminosity of the power-law component we estimate typical masses and accretion rates of the AGN, but we find little evidence of powerful radio-loud or strong X-ray emitters in our sample. This may indicate that the power-law component originates in a nuclear starburst rather than in an AGN, as expected if the bulk of their young stars are still being formed, or that the AGN is still heavily enshrouded in dust and gas. We find that both alternatives are problematic and that more and better optical, X-ray, and mm-wave observations are needed to elucidate the evolutionary history of PSQ.


title: On graded irreducible representations of Leavitt path algebras

abstract: Using the E-algebraic systems, various graded irreducible representations of a Leavitt path algebra L of a graph E over a field K are constructed. The concept of a Laurent vertex is introduced and it is shown that the minimal graded left ideals of L are generated by Laurent vertices and/or line points leading to a detailed description of the graded socle of L. Following this, a complete characterization is obtained of the Leavitt path algebras over which every graded irreducible representation is finitely presented. A useful result is that the irreducible representation V_[p] induced by infinite path tail-equivalent to an infinite path p is graded if and only if p is an irrational infinite path. We also show that every one-sided ideal of L is graded if and only if the graph E contains no cycles. Supplementing a theorem of one of the co-authors that every Leavitt path algebra is graded von Neumann regular, we describe the graded self-injective Leavitt path algebras. These turn out to be direct sums of matrix rings of arbitrary size over K and k[x,x^-1].


title: XMM-Newton observations of SGR 1806-20 over seven years following the 2004 Giant Flare

abstract: We report on the study of 14 XMM-Newton observations of the magnetar SGR 1806-20 spread over a period of 8 years, starting in 2003 and extending to 2011. We find that in mid 2005, a year and a half after a giant flare (GF), the torques on the star increased to the largest value yet seen, with a long term average rate between 2005 and 2011 of $\lvert\dot{\nu}\rvert\approx1.35\times10^{-11}$ Hz s$^{-1}$, an order of magnitude larger than its historical level measured in 1995. The pulse morphology of the source is complex in the observations following the GF, while its pulsed-fraction remained constant at about $7\%$ in all observations. Spectrally, the combination of a black-body (BB) and power-law (PL) components is an excellent fit to all observations. The BB and PL fluxes increased by a factor of 2.5 and 4, respectively, while the spectra hardened, in concordance with the 2004 major outburst that preceded the GF. The fluxes decayed exponentially back to quiescence with a characteristic time-scale of $\tau\sim1.5$ yrs, although they did not reach a constant value until at least 3.5 years later (2009). The long-term timing and spectral behavior of the source point to a decoupling between the mechanisms responsible for their respective behavior. We argue that low level seismic activity causing small twists in the open field lines can explain the long lasting large torques on the star, while the spectral behavior is due to a twist imparted onto closed field lines after the 2004 large outburst.


title: Practical Run-time Checking via Unobtrusive Property Caching

abstract: The use of annotations, referred to as assertions or contracts, to describe program properties for which run-time tests are to be generated, has become frequent in dynamic programing languages. However, the frameworks proposed to support such run-time testing generally incur high time and/or space overheads over standard program execution. We present an approach for reducing this overhead that is based on the use of memoization to cache intermediate results of check evaluation, avoiding repeated checking of previously verified properties. Compared to approaches that reduce checking frequency, our proposal has the advantage of being exhaustive (i.e., all tests are checked at all points) while still being much more efficient than standard run-time checking. Compared to the limited previous work on memoization, it performs the task without requiring modifications to data structure representation or checking code. While the approach is general and system-independent, we present it for concreteness in the context of the Ciao run-time checking framework, which allows us to provide an operational semantics with checks and caching. We also report on a prototype implementation and provide some experimental results that support that using a relatively small cache leads to significant decreases in run-time checking overhead.


title: Initial-state cold nuclear matter energy loss effects on inclusive jet production in p+A collisions at RHIC and LHC

abstract: Recent measurements of the centrality and rapidity dependence of single inclusive jet production in p+Pb collisions at the LHC have revealed large and non-trivial nuclear modification of the production cross section for this process. In this paper, we explore to what extent such nuclear modification can be understood by the framework of standard cold nuclear matter effects, in particular initial-state cold nuclear matter energy loss. We demonstrate quantitatively that theoretical calculations which include medium-induced radiative corrections can describe rather reasonably the attenuation of the jet production yields in the large transverse momentum region in d+Au collisions at RHIC and p+Pb collisions at the LHC for central to semi-central collisions. We further show that the observed scaling behavior of the nuclear modification factor as a function of the total jet energy $p_T\, {\cosh} (y)$ for various rapidity intervals has a natural explanation in the picture of cold nuclear matter energy loss. On the other hand, the observed enhancement in peripheral collisions is not described in this picture and could have a different origin.


title: Parallel simulation for the ultra-short laser pulses' propagation in air

abstract: A parallel 2D+1 split-step Fourier method with Crank-Nicholson scheme running on multi-core shared memory architectures is developed to study the propagation of ultra-short high-intensity laser pulses in air. The parallel method achieves a near linear speed-up with results for the efficiency of more than 95% on a 24-core machine. This method is of great potential application in studying the long-distance propagation of the ultra-short high intensity laser pulses.


title: Tweeting Over The Border: An Empirical Study of Transnational Migration in San Diego and Tijuana

abstract: Sociological studies on transnational migration are often based on surveys or interviews, an expensive and time consuming approach. On the other hand, the pervasiveness of mobile phones and location aware social networks has introduced new ways to understand human mobility patterns at a national or global scale. In this work, we leverage geo located information obtained from Twitter as to understand transnational migration patterns between two border cities (San Diego, USA and Tijuana, Mexico). We obtained 10.9 million geo located tweets from December 2013 to January 2015. Our method infers human mobility by inspecting tweet submissions and user's home locations. Our results depict a trans national community structure that exhibits the formation of a functional metropolitan area that physically transcends international borders. These results show the potential for re analysing sociology phenomena from a technology based empirical perspective.


title: Estimation and uncertainty of reversible Markov models

abstract: Reversibility is a key concept in Markov models and Master-equation models of molecular kinetics. The analysis and interpretation of the transition matrix encoding the kinetic properties of the model relies heavily on the reversibility property. The estimation of a reversible transition matrix from simulation data is therefore crucial to the successful application of the previously developed theory. In this work we discuss methods for the maximum likelihood estimation of transition matrices from finite simulation data and present a new algorithm for the estimation if reversibility with respect to a given stationary vector is desired. We also develop new methods for the Bayesian posterior inference of reversible transition matrices with and without given stationary vector taking into account the need for a suitable prior distribution preserving the meta- stable features of the observed process during posterior inference. All algorithms here are implemented in the PyEMMA software - http://pyemma.org - as of version 2.0.

categories: physics.chem-ph math-ph math.MP physics.comp-ph physics.data-an stat.CO


title: Design Contracts For Networked Automation Systems Co-design

abstract: Networked automation systems (NAS) are characterized by confluence of control, computation, communication and Information (C3I) technologies. Design decisions of one domain are affected by the constraints posed by others. Reliable NAS design should address the requirements of the system, and simultaneously meet the constraints posed by other domains and this is called co-design in literature. Co-design requires clear definition of interfaces among these domains. Control design in NAS is affected by the timing imperfections posed by other domains. In this investigation, we first study the different sources of timing imperfections in NAS, and classify them based on their occurrence. The concept of jitter is used to define the timing imperfections induced by various system components. Using this analysis, we classify the jitter based on their behavior and domain of occurrence. Our analysis shows that the jitter induced in NAS can be classified based on domain as- hardware, software and communication. Next, we use this analysis to model the jitter from the components of NAS. Modeling timing imperfections helps in capturing the interfaces among the domains, and we use the concept of design contracts to capture the interfaces. Design contracts describe the semantic mapping among the domains and are specified using the jitter margins. Implementing design contracts requires knowledge of the jitter margin and, the results from control theory are used to this extent.


title: Finite Dimensionality of the Global Attractor for the Solutions to 3D Primitive Equations with Viscosity

abstract: A new method is presented to prove finiteness of the fractal and Hausdorff dimensions of the global attractor for the strong solutions to the 3D Primitive Equations with viscosity, which is applicable to even more general situations than the recent result of [7] in the sense that it removes all extra technical conditions imposed by previous analyses. More specifically, for finiteness of the dimensions of the global attractor, we only need the heat source $Q\in L^2$ which is exactly the condition for the existence of global strong solutions and the existence of the global attractor of these solutions; while the best existing result, which was obtained very recently in [7], still needs the extra condition that $\partial_z Q\in L^2 $ for finiteness of the dimensions of the global attractor. Moreover, the new method can be applied to cases with more complicated boundary conditions which present essential difficulties for previous methods.


title: Local observability of state variables and parameters in nonlinear modeling quantified by delay reconstruction

abstract: Features of the Jacobian matrix of the delay coordinates map are exploited for quantifying the robustness and reliability of state and parameter estimations for a given dynamical model using an observed time series. Relevant concepts of this approach are introduced and illustrated for discrete and continuous time systems employing a filtered H\'enon map and a R\"ossler system.


title: Massive MIMO in Real Propagation Environments: Do All Antennas Contribute Equally?

abstract: Massive MIMO can greatly increase both spectral and transmit-energy efficiency. This is achieved by allowing the number of antennas and RF chains to grow very large. However, the challenges include high system complexity and hardware energy consumption. Here we investigate the possibilities to reduce the required number of RF chains, by performing antenna selection. While this approach is not a very effective strategy for theoretical independent Rayleigh fading channels, a substantial reduction in the number of RF chains can be achieved for real massive MIMO channels, without significant performance loss. We evaluate antenna selection performance on measured channels at 2.6 GHz, using a linear and a cylindrical array, both having 128 elements. Sum-rate maximization is used as the criterion for antenna selection. A selection scheme based on convex optimization is nearly optimal and used as a benchmark. The achieved sum-rate is compared with that of a very simple scheme that selects the antennas with the highest received power. The power-based scheme gives performance close to the convex optimization scheme, for the measured channels. This observation indicates a potential for significant reductions of massive MIMO implementation complexity, by reducing the number of RF chains and performing antenna selection using simple algorithms.


title: A Warm Restart Strategy for Solving Sudoku by Sparse Optimization Methods

abstract: This paper is concerned with the popular Sudoku problem. We proposed a warm restart strategy for solving Sudoku puzzles, based on the sparse optimization technique. Furthermore, we defined a new difficulty level for Sudoku puzzles. The efficiency of the proposed method is tested using a dataset of Sudoku puzzles, and the numerical results show that the accurate recovery rate can be enhanced from 84%+ to 99%+ using the L1 sparse optimization method.


title: Origin of the large anisotropic g-factor of holes in bismuth

abstract: The ratio of the Zeeman splitting to the cyclotron energy ($M=\Delta E_{\rm Z}/\hbar \omega_{\rm c}$) for hole-like carriers in bismuth has been quantified with a great precision by many experiments performed during the past five decades. It exceeds 2 when the magnetic field is along the trigonal axis and vanishes in the perpendicular configuration. Theoretically, however, $M$ is expected to be isotropic and equal to unity in a two-band Dirac model. We argue that a solution to this half-a-century-old puzzle can be found by extending the $k\cdot p$ theory to multiple bands. Our model not only gives a quantitative account of magnitude and anisotropy of $M$ for hole-like carriers in bismuth, but also explains its contrasting evolution with antimony doping pressure, both probed by new experiments reported here. The present results have important implications for the magnitude and anisotropy of $M$ in other systems with strong spin-orbit coupling.


title: Interior Nucleation of Martensite in a Cubic-to-Tetragonal Phase Transformation

abstract: Using a variational model based on non-linear elasticity we investigate whether in a cubic-to-tetragonal phase transformation it is energetically preferable to nucleate martensite within austenite. Under minimal growth assumptions on the free energy density $W$ away from the wells, we derive explicit upper bounds on $W^{qc}({\mathbb{I}})$, i.e. on the macroscopic free energy density of a region that has been macroscopically deformed by the identity map. The bounds only depend on material parameters, the temperature difference and the growth of $W$ away from the wells. By comparing $W^{qc}({\mathbb{I}})$ and $W({\mathbb{I}})$ we conclude that nucleation is always energetically preferable and are able to provide quantitative upper bounds for the gain in negative energy due to nucleation.


title: Fast Bidirectional Probability Estimation in Markov Models

abstract: We develop a new bidirectional algorithm for estimating Markov chain multi-step transition probabilities: given a Markov chain, we want to estimate the probability of hitting a given target state in $\ell$ steps after starting from a given source distribution. Given the target state $t$, we use a (reverse) local power iteration to construct an `expanded target distribution', which has the same mean as the quantity we want to estimate, but a smaller variance -- this can then be sampled efficiently by a Monte Carlo algorithm. Our method extends to any Markov chain on a discrete (finite or countable) state-space, and can be extended to compute functions of multi-step transition probabilities such as PageRank, graph diffusions, hitting/return times, etc. Our main result is that in `sparse' Markov Chains -- wherein the number of transitions between states is comparable to the number of states -- the running time of our algorithm for a uniform-random target node is order-wise smaller than Monte Carlo and power iteration based algorithms; in particular, our method can estimate a probability $p$ using only $O(1/\sqrt{p})$ running time.


title: Personalized PageRank Estimation and Search: A Bidirectional Approach

abstract: We present new algorithms for Personalized PageRank estimation and Personalized PageRank search. First, for the problem of estimating Personalized PageRank (PPR) from a source distribution to a target node, we present a new bidirectional estimator with simple yet strong guarantees on correctness and performance, and 3x to 8x speedup over existing estimators in experiments on a diverse set of networks. Moreover, it has a clean algebraic structure which enables it to be used as a primitive for the Personalized PageRank Search problem: Given a network like Facebook, a query like "people named John", and a searching user, return the top nodes in the network ranked by PPR from the perspective of the searching user. Previous solutions either score all nodes or score candidate nodes one at a time, which is prohibitively slow for large candidate sets. We develop a new algorithm based on our bidirectional PPR estimator which identifies the most relevant results by sampling candidates based on their PPR; this is the first solution to PPR search that can find the best results without iterating through the set of all candidate results. Finally, by combining PPR sampling with sequential PPR estimation and Monte Carlo, we develop practical algorithms for PPR search, and we show via experiments that our algorithms are efficient on networks with billions of edges.


title: On-line Survival Analysis of Power Electronic Converters Using Step Noise-Cox Processes

abstract: This paper is focused on survival analysis of electrical components. The main goal of this paper is to develop a method for on-line estimation of the Mean Time To Failure (MTTF) of electrical components under dynamic stress levels. The proposed method models the variations of the stress levels as a stochastic process. Hence, a stochastic failure rate function can be developed for each electrical component. Later, this function is used as the underlying rate of a doubly stochastic Poisson process (known as Cox processes). Furthermore, this Cox process is used for on-line estimation of the Mean Residual Life (MRL) using the observed stress levels. The proposed method provides a good estimate of the age and life expectancy of each component. An experimental case study is provided to demonstrate the proposed method.


title: Angular distribution in s-channel formation of the pentaquark $\Theta^+$ baryon

abstract: Using the DIANA data on the charge-exchange reaction $K^+n \rightarrow pK^0$ on a bound neutron, in which the s-channel formation of the pentaquark baryon $\Theta^+(1538)$ has been observed, we analyze the dependence of the background-subtracted $\Theta^+ \rightarrow pK^0$ signal on the $K^0$ emission angle in the $pK^0$ rest frame. In order to describe the observed $\cos\Theta_K^\mathrm{cms}$ distribution, invoking the interference between the nonresonant s-wave and the $\Theta^+$-mediated higher-wave contributions to the amplitude of the charge-exchange reaction is required at a 2.8$\sigma$ level. The spin--parity assignment of 1/2$^-$ for the $\Theta^+$ baryon is ruled out at a statistical level of 2.9 standard deviations. A physically-meaningful selection in $\cos\Theta_K^\mathrm{cms}$ based on the observed angular dependence of the $\Theta^+ \rightarrow pK^0$ signal allows to boost the statistical significance of the signal up to 7.1 standard deviations. This is far in excess of previously reported signals and renders the $\Theta^+$ existence more credible.


title: Epidemic Spreading in Random Rectangular Networks

abstract: The use of network theory to model disease propagation on populations introduces important elements of reality to the classical epidemiological models. The use of random geometric graphs (RGG) is one of such network models that allows for the consideration of spatial properties on disease propagation. In certain real-world scenarios -like in the analysis of a disease propagating through plants- the shape of the plots and fields where the host of the disease is located may play a fundamental role on the propagation dynamics. Here we consider a generalization of the RGG to account for the variation of the shape of the plots/fields where the hosts of a disease are allocated. We consider a disease propagation taking place on the nodes of a random rectangular graph (RRG) and we consider a lower bound for the epidemic threshold of a Susceptible-Infected-Susceptible (SIS) or Susceptible-Infected-Recovered (SIR) model on these networks. Using extensive numerical simulations and based on our analytical results we conclude that (ceteris paribus) the elongation of the plot/field in which the nodes are distributed makes the network more resilient to the propagation of a disease due to the fact that the epidemic threshold increases with the elongation of the rectangle. These results agree with accumulated empirical evidence and simulation results about the propagation of diseases on plants in plots/fields of the same area and different shapes.


title: Fast phase randomisation via two-folds

abstract: A two-fold is a singular point on the discontinuity surface of a piecewise-smooth vector field, at which the vector field is tangent to the discontinuity surface on both sides. If an orbit passes through an invisible two-fold (also known as a Teixeira singularity) before settling to regular periodic motion, then the phase of that motion cannot be determined from initial conditions, and in the presence of small noise the asymptotic phase of a large number of sample solutions is highly random. In this paper we show how the probability distribution of the asymptotic phase depends on the global nonlinear dynamics. We also show how the phase of a smooth oscillator can be randomised by applying a simple discontinuous control law that generates an invisible two-fold. We propose that such a control law can be used to desynchronise a collection of oscillators, and that this manner of phase randomisation is fast compared to existing methods (which use fixed points as phase singularities) because there is no slowing of the dynamics near a two-fold.


title: Sequence Heterogeneity Accelerates Protein Search for Targets on DNA

abstract: The process of protein search for specific binding sites on DNA is fundamentally important since it marks the beginning of all major biological processes. We present a theoretical investigation that probes the role of DNA sequence symmetry, heterogeneity and chemical composition in the protein search dynamics. Using a discrete-state stochastic approach with a first-passage events analysis, which takes into account the most relevant physical-chemical processes, a full analytical description of the search dynamics is obtained. It is found that, contrary to existing views, the protein search is generally faster on DNA with more heterogeneous sequences. In addition, the search dynamics might be affected by the chemical composition near the target site. The physical origins of these phenomena are discussed. Our results suggest that biological processes might be effectively regulated by modifying chemical composition, symmetry and heterogeneity of a genome.


title: Two-component integrable generalizations of Burgers equations with nondiagonal linearity

abstract: Two-component second and third-order Burgers type systems with nondiagonal constant matrix of leading order terms are classified for higher symmetries. New symmetry integrable systems with their master symmetries are obtained. Some third order systems are observed to possess conservation laws. Bi-Poisson structures of systems possessing conservation laws are given.


title: Contractible edges in 3-connected graphs that preserve a minor

abstract: Let $G$ be a $3$-connected graph with a $3$-connected (or sufficiently small) simple minor $H$. We establish that $G$ has a forest $F$ with at least $\left\lceil(|G|-|H|+1)/2\right\rceil$ edges such that $G/e$ is $3$-connected with an $H$-minor for each $e\in E(F)$. Moreover, we may pick $F$ with $|G|-|H|$ edges provided $G$ is triangle-free. These results are sharp. Our result generalizes a previous one by Ando et. al., which establishes that a $3$-connected graph $G$ has at least $\left\lceil|G|/2\right\rceil$ contractible edges. As another consequence, each triangle-free $3$-connected graph has an spanning tree of contractible edges. Our results follow from a more general theorem on graph minors, a splitter theorem, which is also established here.


title: Higgs relaxation and the matter-antimatter asymmetry of the universe

abstract: The recent measurement of the Higgs boson mass implies a relatively slow rise of the Standard Model Higgs potential at large scales, and a possible second minimum at even larger scales. Consequently, the Higgs field may develop a large vacuum expectation value during inflation. The relaxation of the Higgs field from its large postinflationary value to the minimum of the effective potential represents an important stage in the evolution of the universe. During this epoch, the time-dependent Higgs condensate can create an effective chemical potential for the lepton number, leading to a generation of the lepton asymmetry in the presence of some large right-handed Majorana neutrino masses. Electroweak sphalerons redistribute this asymmetry between leptons and baryons. Higgs relaxation leptogenesis can explain the observed matter-antimatter asymmetry of the universe even if the Standard Model is valid up to the scale of inflation, and any new physics is suppressed by that high scale.


title: Parabolic Anderson model in a dynamic random environment: random conductances

abstract: The parabolic Anderson model is defined as the partial differential equation \partial u(x,t)/\partial t = \kappa\Delta u(x,t) + \xi(x,t)u(x,t), x\in\Z^d, t\geq 0, where \kappa \in [0,\infty) is the diffusion constant, \Delta is the discrete Laplacian, and \xi is a dynamic random environment that drives the equation. The initial condition u(x,0)=u_0(x), x\in\Z^d, is taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate 2d\kappa, split into two at rate \xi \vee 0, and die at rate (-\xi) \vee 0. Our focus is on the Lyapunov exponents \lambda_p(\kappa) = \lim_{t\to\infty} \frac{1}{t} \log \E([u(0,t)]^p)^{1/p}, p \in \N, and \lambda_0(\kappa) = \lim_{t\to\infty} \frac{1}{t}\log u(0,t). We investigate what happens when \kappa\Delta is replaced by \Delta^\cK, where \cK = \{\mathcal{K}(x,y)\colon\,x,y\in\Z^d,\,x \sim y\} is a collection of random conductances between neighbouring sites replacing the constant conductances \kappa in the homogeneous model. We show that the associated annealed Lyapunov exponents are given by the formula \lambda_p(\cK) = \sup\{\lambda_p(\kappa) \colon\,\kappa\in\Supp(\cK)\}, where \Supp(\cK) is the set of values taken by the \cK-field. We also show that for the associated quenched Lyapunov exponent this formula only provides a lower bound. Our proof is valid for three classes of reversible \xi, and for all \cK satisfying a certain clustering property, namely, there are arbitrarily large balls where \cK is almost constant and close to any value in \Supp(\cK). What our result says is that the Lyapunov exponents are controlled by those pockets of \cK where the conductances are close to the value that maximises the growth in the homogeneous setting.


title: Valuations and Frobenius

abstract: The behavior of the Frobenius map is investigated for valuation rings of prime characteristic. We show that valuation rings are always F-pure. We introduce a generalization of the notion of strong F-regularity, which we call F-pure regularity, and show that a valuation ring is F-pure regular if and only if it is Noetherian. For valuations on function fields, we show that the Frobenius map is finite if and only if the valuation is divisorial; in this case the valuation ring is Frobenius split. For Noetherian valuation rings in function fields, we show that the valuation ring is Frobenius split if and only if Frobenius is finite, or equivalently, if and only if the valuation ring is excellent.


title: Structural Resolution for Logic Programming

abstract: We introduce a Three Tier Tree Calculus (3TC) that defines in a systematic way three tiers of tree structures underlying proof search in logic programming. We use 3TC to define a new -- structural -- version of resolution for logic programming.


title: Using coarse GPS data to quantify city-scale transportation system resilience to extreme events

abstract: This article proposes a method to quantitatively measure the resilience of transportation systems using GPS data from taxis. The granularity of the GPS data necessary for this analysis is relatively coarse; it only requires coordinates for the beginning and end of trips, the metered distance, and the total travel time. The method works by computing the historical distribution of pace (normalized travel times) between various regions of a city and measuring the pace deviations during an unusual event. This method is applied to a dataset of nearly 700 million taxi trips in New York City, which is used to analyze the transportation infrastructure resilience to Hurricane Sandy. The analysis indicates that Hurricane Sandy impacted traffic conditions for more than five days, and caused a peak delay of two minutes per mile. Practically, it identifies that the evacuation caused only minor disruptions, but significant delays were encountered during the post-disaster reentry process. Since the implementation of this method is very efficient, it could potentially be used as an online monitoring tool, representing a first step toward quantifying city scale resilience with coarse GPS data.


title: Implicit Large Eddy Simulation of a wingtip vortex at $Re_c = 1.2\cdot 10^6$

abstract: In this article we present recent developments in numerical methods for performing a Large Eddy Simulation (LES) of the formation and evolution of a wingtip vortex. The development of these vortices in the near wake, in combination with the large Reynolds numbers present in these cases, make these types of test cases particularly challenging to investigate numerically. We first give an overview of the Spectral Vanishing Viscosity--implicit LES (SVV-iLES) solver that is used to perform the simulations, and highlight techniques that have been adopted to solve various numerical issues that arise when studying such cases. To demonstrate the method's viability, we present results from numerical simulations of flow over a NACA 0012 profile wingtip at $Re_c = 1.2\cdot 10^6$ and compare them against experimental data, which is to date the highest Reynolds number achieved for a LES that has been correlated with experiments for this test case. Our model correlates favorably with experiment, both for the characteristic jetting in the primary vortex and pressure distribution on the wing surface. The proposed method is of general interest for the modeling of transitioning vortex dominated flows over complex geometries.


title: Large Complex Correlated Wishart Matrices: The Pearcey Kernel and Expansion at the Hard Edge

abstract: We study the eigenvalue behaviour of large complex correlated Wishart matrices near an interior point of the limiting spectrum where the density vanishes (cusp point), and refine the existing results at the hard edge as well. More precisely, under mild assumptions for the population covariance matrix, we show that the limiting density vanishes at generic cusp points like a cube root, and that the local eigenvalue behaviour is described by means of the Pearcey kernel if an extra decay assumption is satisfied. As for the hard edge, we show that the density blows up like an inverse square root at the origin. Moreover, we provide an explicit formula for the $1/N$ correction term for the fluctuation of the smallest random eigenvalue.


title: Spectral and temporal properties of the ultra-luminous X-ray pulsar in M82 from 15 years of Chandra observations and analysis of the pulsed emission using NuSTAR

abstract: The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit.


title: Risk Quantification in Stochastic Simulation under Input Uncertainty

abstract: When simulating a complex stochastic system, the behavior of output response depends on input parameters estimated from finite real-world data, and the finiteness of data brings input uncertainty into the system. The quantification of the impact of input uncertainty on output response has been extensively studied. Most of the existing literature focuses on providing inferences on the mean response at the true but unknown input parameter, including point estimation and confidence interval construction. Risk quantification of mean response under input uncertainty often plays an important role in system evaluation and control, because it provides inferences on extreme scenarios of mean response in all possible input models. To the best of our knowledge, it has rarely been systematically studied in the literature. In this paper, first we introduce risk measures of mean response under input uncertainty, and propose a nested Monte Carlo simulation approach to estimate them. Then we develop asymptotical properties such as consistency and asymptotic normality for the proposed nested risk estimators. We further study the associated budget allocation problem for efficient nested risk simulation, and finally use a sharing economy example to illustrate the importance of accessing and controlling risk due to input uncertainty.


title: Three-Dimensional Interface Roughness in Layered Semiconductor Structures and Its Effects on Intersubband Transitions

abstract: A general model for treating the effects of three dimensional interface roughness (IFR) in layered semiconductor structures has been derived and experimentally verified. Configurational averaging of the IFR potential produces an effective grading potential in the out-of-plane direction, greatly altering the energy spectrum of the structures. IFR scattering self-energy is also derived for the general case; when IFR is strong, its scattering effect is shown to dominate over phonon interaction and impurity scattering. When applied to intersubband transitions, the theoretical predictions explain the experimental observation of the anomalous energy shift and unusual broadening of the ISB transitions in III-Nitride thin-layered superlattices.


title: Ward Identities for Scale and Special Conformal Transformations in Inflation

abstract: We derive the general Ward identities for scale and special conformal transformations in theories of single field inflation. Our analysis is model independent and based on symmetry considerations alone. The identities we obtain are valid to all orders in the slow roll expansion. For special conformal transformations, the Ward identities include a term which is non-linear in the fields that arises due to a compensating spatial reparametrization. Some observational consequences are also discussed.


title: Low Scale Composite Higgs Model and 1.8 $\sim$ 2 TeV Diboson Excess

abstract: We consider a simple solution to explain the recent diboson excess observed by ATALS and CMS Collaborations in models with custodial symmetry $SU(2)_L \times SU(2)_R \rightarrow SU(2)_c$. The $SU(2)_L$ triplet vector boson $\rho$ with mass range of $1.8 \sim 2$ TeV would be produced through the Drell-Yan process with sizable diboson decay branching to account for the excess. The other $SU(2)_L \times SU(2)_R$ bidoublet axial vector boson $a$ would cancel all deviations of electroweak obervables induced by $\rho$ even if the SM fermions mix with some heavy vector like (composite) fermions which couple to $\rho$ ("non-universally partially composite"), therefore allows arbitrary couplings between each SM fermion and $\rho$. We present our model in the "General Composite Higgs" framework with $SO(5) \times U(1)_X \rightarrow SO(4) \times U(1)_X$ breaking at scale $f$ and demand the first Weinberg sum rule and positive gauge boson form factors as the theoretical constraints. We find that our model can fit the diboson excess very well if the left-handed SM light quarks, charged leptons and tops have zero, zero/moderately small and moderate/large composite components for reasonable values of $g_\rho$ and $f$. The correlation between tree level $S$ parameter and the $h \rightarrow Z \gamma$ suggest a large $a$ contribution to $h \rightarrow Z \gamma$ and it is indeed a $\mathcal{O}(1)$ effect in our parameter space which provides a strong hint for our scenario if this diboson excess is confirmed by the $13 \sim 14$ TeV LHC Run II.


title: The Total Acquisition Number of the Randomly Weighted Path

abstract: There exists a significant body of work on determining the acquisition number $a_t(G)$ of various graphs when the vertices of those graphs are each initially assigned a unit weight. We determine properties of the acquisition number of the path, star, complete, complete bipartite, cycle, and wheel graphs for variations on this initial weighting scheme, with the majority of our work focusing on the expected acquisition number of randomly weighted graphs. In particular, we bound the expected acquisition number $E(a_t(P_n))$ of the $n$-path when $n$ distinguishable "units" of integral weight, or chips, are randomly distributed across its vertices between $0.242n$ and $0.375n$. With computer support, we improve it by showing that $E(a_t(P_n))$ lies between $0.29523n$ and $0.29576n$. We then use subadditivity to show that the limiting ratio $\lim E(a_t(P_n))/n$ exists, and simulations reveal more exactly what the limiting value equals. The Hoeffding-Azuma inequality is used to prove that the acquisition number is tightly concentrated around its expected value. Additionally, in a different context, we offer a non-optimal acquisition protocol algorithm for the randomly weighted path and exactly compute the expected size of the resultant residual set.


title: Practical Selection of SVM Supervised Parameters with Different Feature Representations for Vowel Recognition

abstract: It is known that the classification performance of Support Vector Machine (SVM) can be conveniently affected by the different parameters of the kernel tricks and the regularization parameter, C. Thus, in this article, we propose a study in order to find the suitable kernel with which SVM may achieve good generalization performance as well as the parameters to use. We need to analyze the behavior of the SVM classifier when these parameters take very small or very large values. The study is conducted for a multi-class vowel recognition using the TIMIT corpus. Furthermore, for the experiments, we used different feature representations such as MFCC and PLP. Finally, a comparative study was done to point out the impact of the choice of the parameters, kernel trick and feature representations on the performance of the SVM classifier


title: An Empirical Comparison of SVM and Some Supervised Learning Algorithms for Vowel recognition

abstract: In this article, we conduct a study on the performance of some supervised learning algorithms for vowel recognition. This study aims to compare the accuracy of each algorithm. Thus, we present an empirical comparison between five supervised learning classifiers and two combined classifiers: SVM, KNN, Naive Bayes, Quadratic Bayes Normal (QDC) and Nearst Mean. Those algorithms were tested for vowel recognition using TIMIT Corpus and Mel-frequency cepstral coefficients (MFCCs).


title: Divergence-free Approximate Riemann Solver for the Quasi-neutral Two-fluid Plasma Model

abstract: A numerical method for the quasi-neutral two-fluid (QNTF) plasma model is described. The basic equations are ion and electron fluid equations and the Maxwell equations without displacement current. The neglect of displacement current is consistent with the assumption of charge neutrality. It thus reduces to the ideal magnetohydrodynamic (MHD) equations in the long wavelength limit, but the two-fluid effect appearing at ion and electron inertial scales is fully taken into account. It is shown that the basic equations may be rewritten in a form that has formally the same structure as the MHD equations. The total mass, momentum, and energy are all written in the conservative form. A new three-dimensional numerical simulation code has been developed for the QNTF equations. The HLL (Harten-Lax-van Leer) approximate Riemann solver combined with the upwind constrained transport (UCT) scheme is applied. The method was originally developed for MHD (Londrillo & Del Zanna, 2004), but works quite well for the present model as well. The simulation code is able to capture sharp multidimensional discontinuities as well as dispersive waves arising from the two-fluid effect at small scales without producing $\nabla \cdot \mathbf{B}$ errors. It is well known that conventional Hall-MHD codes often suffer a numerical stability issue associated with short wavelength whistler waves. On the other hand, since finite electron inertia introduces an upper bound to the phase speed of whistler waves in the present model, our code is free from the issue even without explicit dissipation terms or implicit time integration. Numerical experiments have confirmed that there is no need to resolve characteristic time scales such as plasma frequency or cyclotron frequency for numerical stability. Consequently, the QNTF model offers a better alternative to the Hall-MHD or fully electromagnetic two-fluid models.


title: Robust speech recognition using consensus function based on multi-layer networks

abstract: The clustering ensembles mingle numerous partitions of a specified data into a single clustering solution. Clustering ensemble has emerged as a potent approach for ameliorating both the forcefulness and the stability of unsupervised classification results. One of the major problems in clustering ensembles is to find the best consensus function. Finding final partition from different clustering results requires skillfulness and robustness of the classification algorithm. In addition, the major problem with the consensus function is its sensitivity to the used data sets quality. This limitation is due to the existence of noisy, silence or redundant data. This paper proposes a novel consensus function of cluster ensembles based on Multilayer networks technique and a maintenance database method. This maintenance database approach is used in order to handle any given noisy speech and, thus, to guarantee the quality of databases. This can generates good results and efficient data partitions. To show its effectiveness, we support our strategy with empirical evaluation using distorted speech from Aurora speech databases.


title: The ground state construction of bilayer graphene

abstract: We consider a model of half-filled bilayer graphene, in which the three dominant Slonczewski-Weiss-McClure hopping parameters are retained, in the presence of short range interactions. Under a smallness assumption on the interaction strength $U$ as well as on the inter-layer hopping $\epsilon$, we construct the ground state in the thermodynamic limit, and prove its analyticity in $U$, uniformly in $\epsilon$. The interacting Fermi surface is degenerate, and consists of eight Fermi points, two of which are protected by symmetries, while the locations of the other six are renormalized by the interaction, and the effective dispersion relation at the Fermi points is conical. The construction reveals the presence of different energy regimes, where the effective behavior of correlation functions changes qualitatively. The analysis of the crossover between regimes plays an important role in the proof of analyticity and in the uniform control of the radius of convergence. The proof is based on a rigorous implementation of fermionic renormalization group methods, including determinant estimates for the renormalized expansion.

categories: math-ph cond-mat.stat-mech cond-mat.str-el math.MP


title: Incorporating Belief Function in SVM for Phoneme Recognition

abstract: The Support Vector Machine (SVM) method has been widely used in numerous classification tasks. The main idea of this algorithm is based on the principle of the margin maximization to find an hyperplane which separates the data into two different classes.In this paper, SVM is applied to phoneme recognition task. However, in many real-world problems, each phoneme in the data set for recognition problems may differ in the degree of significance due to noise, inaccuracies, or abnormal characteristics; All those problems can lead to the inaccuracies in the prediction phase. Unfortunately, the standard formulation of SVM does not take into account all those problems and, in particular, the variation in the speech input. This paper presents a new formulation of SVM (B-SVM) that attributes to each phoneme a confidence degree computed based on its geometric position in the space. Then, this degree is used in order to strengthen the class membership of the tested phoneme. Hence, we introduce a reformulation of the standard SVM that incorporates the degree of belief. Experimental performance on TIMIT database shows the effectiveness of the proposed method B-SVM on a phoneme recognition problem.


title: The positive piezoconductive effect in graphene

abstract: As the thinnest conductive and elastic material, graphene is expected to play a crucial role in post-Moore era. Besides applications on electronic devices, graphene has shown great potential for nano-electromechanical systems. While interlayer interactions play a key role in modifying the electronic structures of layered materials, no attention has been given to their impact on electromechanical properties. Here we report the positive piezoconductive effect observed in suspended bi- and multi-layer graphene. The effect is highly layer number dependent and shows the most pronounced response for tri-layer graphene. The effect, and its dependence on the layer number, can be understood as resulting from the strain-induced competition between interlayer coupling and intralayer transport, as confirmed by the numerical calculations based on the non-equilibrium Green's function method. Our results enrich the understanding of graphene and point to layer number as a powerful tool for tuning the electromechanical properties of graphene for future applications.


title: Centrifugal deterministic lateral displacement separation system

abstract: This work investigates the migration of spherical particles of different sizes in a centrifuge-driven deterministic lateral displacement (c-DLD) device. Specifically, we use a scaled-up model to study the motion of suspended particles through a square array of cylindrical posts under the action of centrifugation. Experiments show that separation of particles by size is possible depending on the orientation of the driving acceleration with respect to the array of posts (forcing angle). We focus on the fractionation of binary suspensions and measure the separation resolution at the outlet of the device for different forcing angles. We found excellent resolution at intermediate forcing angles, when large particles are locked to move at small migration angles but smaller particles follow the forcing angle more closely. Finally, we show that reducing the initial concentration (number) of particles, approaching the dilute limit of single particles, leads to increased resolution in the separation.


title: The challenges of SVM optimization using Adaboost on a phoneme recognition problem

abstract: The use of digital technology is growing at a very fast pace which led to the emergence of systems based on the cognitive infocommunications. The expansion of this sector impose the use of combining methods in order to ensure the robustness in cognitive systems.


title: Towards Input Device Satisfaction Through Hand Anthropometry

abstract: We collected the hand anthropometric data of 91 respondents to come up with a Filipino-based measurement to determine the suitability of an input device for a digital equipment, the standard PC keyboard. For correlation purposes, we also collected other relevant information like age, height, province of origin, and gender, among others. We computed the percentiles for each finger to classify various finger dimensions and identify length-specific anthropometric cut-points. We compared the percentiles of each finger dimension against the actual length of the longest key combinations when correct finger placement is used for typing, to determine whether the standard PC keyboard is fit for use by our sampled population. Our analysis shows that the members of the population with hand dimensions at extended position below 75th percentile and at 99th percentile are the ones who would most likely not reach the longest key combination for the left and the right hands, respectively. Using machine vision and image processing techniques, we automated the anthropometric process and compared the accuracy of its measurements to that of manual process'. We compared the measurement generated by our automated anthropometric process with the measurements using the manual one and we found out that they have a very minimal absolute difference. The data collected from this study could be used in other studies such as determining a good design for mobile and other handheld devices, or input devices other than keyboard. The automated method that we developed could be used to easily measure hand dimensions given a digital image of the hand and could be extended for measuring the entire human body for various other applications.


title: Yang-Baxter relation planar algebras

abstract: We give a new type of Schur-Weyl duality for the representations of a family of quantum subgroups and their centralizer algebra. We define and classify singly-generated, Yang-Baxter relation planar algebras. We present the skein theoretic construction of a new parameterized planar algebra. We construct infinitely many new subfactors and unitary fusion categories, and compute their trace formula as a closed-form expression, in terms of Young diagrams.


title: Inversion of the elliptical Radon transform arising in migration imaging using the regular Radon transform

abstract: In recent years, many types of elliptical Radon transforms that integrate functions over various sets of ellipses/ellipsoids have been considered, relating to studies in bistatic synthetic aperture radar, ultrasound reflection tomography, radio tomography, and migration imaging. In this article, we consider the transform that integrates a given function in $\mathbf R^n$ over a set of ellipses (when $n=2$) or ellipsoids of rotation (when $n\geq 3$) with foci restricted to a hyperplane. We show a relation between this elliptical Radon transform and the regular Radon transform, and provide the inversion formula for the elliptical Radon transform using this relation. Numerical simulations are performed to demonstrate the suggested algorithms in two dimensions, and these simulations are also provided in this article.


title: Elastic Net Procedure for Partially Linear Models

abstract: Variable selection plays an important role in the high-dimensional data analysis. However the high-dimensional data often induces the strongly correlated variables problem. In this paper, we propose Elastic Net procedure for partially linear models and prove the group effect of its estimate. By a simulation study, we show that the strongly correlated variables problem can be better handled by the Elastic Net procedure than Lasso, ALasso and Ridge. Based on an empirical analysis, we can get that the Elastic Net procedure is particularly useful when the number of predictors $p$ is much bigger than the sample size $n$.


title: Influence of clustering and hadron potentials on the rapidity distribution of protons from the UrQMD model

abstract: The Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model supplemented by potentials for both pre-formed hadrons and confined baryons (called UrQMD/M) are used to describe rapidity distributions of both the E895 proton data at AGS energies and the NA49 net proton data at SPS energies. With the help of a coalescence afterburner using only one parameter set of ($R_0$, $P_0$)=(3.8 fm, 0.3 GeV$/$c), both sets of experimental data can be described fairly well except for a small discrepancy seen for the net protons at mid-rapidity from heavy ion collisions (HICs) at high SPS energies. Furthermore, in contrast to the logarithmic dependence with beam energy at SIS energies there are still about 10$\%$ of protons in clusters from central HICs at the beam energy of $80$ GeV$/$nucleon.


title: Lithium Intercalation in Graphene/MoS2 Composites: First-Principles Insights

abstract: As a storage material for Li-ion batteries, graphene/molybdenum disulfide (Gr/MoS2) composites have been intensively studied in experiments. But the relevant theoretical works from first-principles are lacking. In the current work, van-der-Waals-corrected density functional theory calculations are performed to investigate the interaction of Li in Gr/MoS2 composites. Three interesting features are revealed for the intercalated Gr/Li(n)/MoS2 composites (n = 1 to 9). One is the reason for large Li storage capacity of Gr/MoS2: due to the binding energies per Li atom increase with the increasing number of intercalated Li atoms. Secondly, the band gap opening of Gr is found, and the band gap is enlarged with the increasing number of intercalated Li atoms, up to 160 meV with nine Li; hence these results suggest an efficient way to tune the band gap of graphene. Thirdly, the Dirac cone of Gr always preserve for different number of ionic bonded Li atoms.


title: The Promise of Quantum Simulation

abstract: Quantum simulation promises to be one of the primary application of quantum computers, should one be constructed. This article briefly summarizes the history quantum simulation in light of the recent result of Wang and coworkers demonstrating calculation of the ground and excited states for a HeH+ molecule, and concludes with a discussion of why this and other recent progress in the field suggests that quantum simulation of quantum chemistry has a bright future.


title: Artificially constructed plasmaron and plasmon-exciton molecule in 2D-metals

abstract: Resonant optical excitation was used to create a macroscopic non-equilibrium ensemble of dark excitons with unprecedentedly long lifetime in a two-dimensional (2D) electron system placed in a quantizing magnetic field. Exotic three-particle and four-particle states, plasmarons and plasmon-exciton molecules, coupled with the surrounding electrons through the collective plasma oscillations are engineered. Plasmarons and plasmon-exciton molecules are manifested as new features in the recombination spectra of non-equilibrium systems. Magnetoplasmarons and magnetotrions in a two-dimensional Hall insulator.


title: The abundance of 28Si32S, 29Si32S, 28Si34S, and 30Si32S in the inner layers of the envelope of IRC+10216

abstract: We present high spectral resolution mid-IR observations of SiS towards the C-rich AGB star IRC+10216 carried out with the Texas Echelon-cross-Echelle Spectrograph mounted on the NASA Infrared Telescope Facility. We have identified 204 ro-vibrational lines of 28Si32S, 26 of 29Si32S, 20 of 28Si34S, and 15 of 30Si32S in the frequency range 720-790 cm-1. These lines belong to bands v=1-0, 2-1, 3-2, 4-3, and 5-4, and involve rotational levels with Jlow<90. About 30 per cent of these lines are unblended or weakly blended and can be partially or entirely fitted with a code developed to model the mid-IR emission of a spherically symmetric circumstellar envelope composed of expanding gas and dust. The observed lines trace the envelope at distances to the star <35R* (~0.7 arcsec). The fits are compatible with an expansion velocity of 1+2.5(r/R*-1) km/s between 1 and 5R*, 11 km/s between 5 and 20R*, and 14.5 km/s outwards. The derived abundance profile of 28Si32S with respect to H2 is 4.9e-6 between the stellar photosphere and 5R*, decreasing linearly to 1.6e-6 at 20R* and to 1.3e-6 at 50R*. 28Si32S seems to be rotationally under LTE in the region of the envelope probed with our observations and vibrationally out of LTE in most of it. There is a red-shifted emission excess in the 28Si32S lines of band v=1-0 that cannot be found in the lines of bands v=2-1, 3-2, 4-3, and 5-4. This excess could be explained by an enhancement of the vibrational temperature around 20R* behind the star. The derived isotopic ratios 28Si/29Si, and 32S/34S are 17 and 14, compatible with previous estimates.


title: Quantum data hiding in the presence of noise

abstract: When classical or quantum information is broadcast to separate receivers, there exist codes that encrypt the encoded data such that the receivers cannot recover it when performing local operations and classical communication, but they can decode reliably if they bring their systems together and perform a collective measurement. This phenomenon is known as quantum data hiding and hitherto has been studied under the assumption that noise does not affect the encoded systems. With the aim of applying the quantum data hiding effect in practical scenarios, here we define the data-hiding capacity for hiding classical information using a quantum channel. Using this notion, we establish a regularized upper bound on the data hiding capacity of any quantum broadcast channel, and we prove that coherent-state encodings have a strong limitation on their data hiding rates. We then prove a lower bound on the data hiding capacity of channels that map the maximally mixed state to the maximally mixed state (we call these channels "mictodiactic"---they can be seen as a generalization of unital channels when the input and output spaces are not necessarily isomorphic) and argue how to extend this bound to generic channels and to more than two receivers.


title: Multi-band Eilenberger theory of superconductivity: Systematic low-energy projection

abstract: We propose the general multi-band quasiclassical Eilenberger theory of superconductivity to describe quasiparticle excitations in inhomogeneous systems. With the use of low-energy projection matrix, the $M$-band quasiclassical Eilenberger equations are systematically obtained from $N$-band Gor'kov equations. Here $M$ is the internal degrees of freedom in the bands crossing the Fermi energy and $N$ is the degree of freedom in a model. Our framework naturally includes inter-band off-diagonal elements of Green's functions, which have usually been neglected in previous multi-band quasiclassical frameworks. The resultant multi-band Eilenberger and Andreev equations are similar to the single-band ones, except for multi-band effects. The multi-band effects can exhibit the non-locality and the anisotropy in the mapped systems. Our framework can be applied to an arbitrary Hamiltonian (e.g. a tight-binding Hamiltonian derived by the first-principle calculation). As examples, we use our framework in various kinds of systems, such as noncentrosymmetric superconductor CePt$_{3}$Si, three-orbital model for Sr$_{2}$RuO$_{4}$, heavy fermion CeCoIn$_{5}$/YbCoIn$_{5}$ superlattice, a topological superconductor with the strong spin-orbit coupling Cu$_{x}$Bi$_{2}$Se$_{3}$, and a surface system on a topological insulator.


title: On a singular minimizing problem

abstract: We study a minimizing problem associated with the singular problem \[ \left\{ \begin{array} [c]{ll} -\operatorname{div}\left( \left\vert \nabla u\right\vert ^{p-2}\nabla u\right) =\lambda u^{-1} & \mathrm{in\ }\Omega\\ u>0 & \mathrm{in\ }\Omega\\ u=0 & \mathrm{on\ }\partial\Omega, \end{array} \right. \] where $p>1$, $\lambda>0$ and $\Omega$ is a bounded and smooth domain of $\mathbb{R}^{N}$, $N\geq2.$ A new log-Sobolev type inequality is proved and the corresponding best constant is identifyied.


title: A stable partitioned FSI algorithm for incompressible flow and deforming beams

abstract: An added-mass partitioned (AMP) algorithm is described for solving fluid-structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier-Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler-Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet-Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.


title: Representation schemes and rigid maximal Cohen-Macaulay modules

abstract: Let k be an algebraically closed field and A be a finitely generated, centrally finite, non- negatively graded (not necessarily commutative) k-algebra. In this note we construct a representation scheme for graded maximal Cohen-Macaulay A modules. Our main application asserts that when A is commutative with an isolated singularity, for a fixed multiplicity, there are only finitely many indecomposable rigid (i.e, with no nontrivial self-extensions) MCM modules up to shifting and isomorphism. We appeal to a result by Keller, Murfet, and Van den Bergh to prove a similar result for rings that are completion of graded rings. Finally, we discuss how finiteness results for rigid MCM modules are related to recent work by Iyama and Wemyss on maximal modifying modules over compound Du Val singularities.


title: Evidence for universal relationship between the measured 1/f permittivity noise and loss tangent created by tunneling atoms

abstract: Noise from atomic tunneling systems (TSs) limit the performance of various resonant devices, ranging in application from astronomy detectors to quantum computing. Using devices containing these TSs, we study the $1/f$ permittivity noise and loss tangent in two film types containing different TS densities. The noise reveals an intrinsic value as well as the crossover to power-saturated noise. The intrinsic $1/f$ noise fits to the temperature dependence $1/T^{1+\mu}$ where $0.2\le\mu\le0.7$, which is related to previous studies and strongly interacting TS. An analysis of the noise normalized by the loss tangent and temperature is quantitatively identical for two film types, despite a factor of 5 difference in their loss tangent. Following from the broad applicability of the TS model, the data supports a universal relationship for amorphous-solid produced permittivity noise. The quantity of the observed noise particularly supports a recent model in which noise is created by weak TS-TS interactions.


title: Dynamical properties of the honeycomb-lattice Iridates ${\rm Na_2IrO_3}$

abstract: We investigate the dynamical properties of ${\rm Na_2IrO_3}$. For five effective models proposed for ${\rm Na_2IrO_3}$, we numerically calculate dynamical structure factors (DSFs) with an exact diagonalization method. An effective model obtained from $ab$ $initio$ calculations explains inelastic neutron scattering experiments adequately. We further calculate excitation modes based on linearized spin-wave theory. The spin-wave excitation of the effective models obtained by $ab$ $initio$ calculations disagrees with the low-lying excitation of DSFs. We attribute this discrepancy to the location of ${\rm Na_2IrO_3}$ in a parameter space close to the phase boundary with the Kitaev spin-liquid phase.


title: Adapting Stochastic Search For Real-time Dynamic Weighted Constraint Satisfaction

abstract: This work presents two new algorithms for performing constraint satisfaction. The first algorithm presented, DMaxWalkSat, is a constraint solver specialized for solving dynamic, weighted constraint satisfaction problems. The second algorithm, RDMaxWalkSat, is a derivative of DMaxWalkSat that has been modified into an anytime algorithm, and hence support realtime constraint satisfaction. DMaxWalkSat is shown to offer performance advantages in terms of solution quality and runtime over its parent constraint solver, MaxWalkSat. RDMaxWalkSat is shown to support anytime operation. The introduction of these algorithms brings another tool to the areas of computer science that naturally represent problems as constraint satisfaction problems, an example of which is the robust coherence algorithm.


title: Uniform Regularity Estimates in Homogenization Theory of Elliptic Systems with Lower Order Terms on the Neumann Boundary Problem

abstract: In this paper, we mainly employed the idea of the previous paper to study the sharp uniform $W^{1,p}$ estimates with $1<p\leq \infty$ for more general elliptic systems with the Neumann boundary condition on a bounded $C^{1,\eta}$ domain, arising in homogenization theory. Based on the skills developed by Z. Shen and by T. Suslina for different purposes, we also established the $L^2$ convergence rates on a bounded $C^{1,1}$ domain and a Lipschitz domain, respectively. Here we found a "rough" version of the first order correctors (see Theorem 1.3), It allows us to skip the corresponding convergence results on $\mathbb{R}^d$ that are the preconditions in T. Suslina's papers.


title: The Effects of Stoichiometry and Substitution in quasi-one-dimensional iron-chalcogenide BaFe$_{2}$S$_{3}$

abstract: The effects of off-stoichiometry and elemental substitution on electronic properties of iron-based ladder compound BaFe$_{2}$S$_{3}$ are investigated. Resistivity and magnetization are revealed to be quite sensitive to the stoichiometry of Fe atoms, and 10 \% deficiency at Fe sites reduces the antiferromagnetic transition temperature by 40 K. The antiferromagnetic transition temperature decreases even faster and collapse to zero with hole doping through 10 \% K substitution at Ba site, while the antiferromagnetic ordering phase remains with electron doping through 20 \% Co substitution at Fe site. Such electron-hole asymmetry is opposite to two-dimensional iron-based superconductors, and can be explained on the basis of both itinerant and localized electronic pictures.


title: Strongly enhanced effects of Lorentz symmetry violation in entangled Yb+ ions

abstract: Lorentz symmetry is one of the cornerstones of modern physics. However, a number of theories aiming at unifying gravity with the other fundamental interactions including string field theory suggest violation of Lorentz symmetry [1-4]. While the energy scale of such strongly Lorentz symmetry-violating physics is much higher than that currently attainable by particle accelerators, Lorentz violation may nevertheless be detectable via precision measurements at low energies [2]. Here, we carry out a systematic theoretical investigation of the sensitivity of a wide range of atomic systems to violation of local Lorentz invariance (LLI). Aim of these studies is to identify which atom shows the biggest promise to detect violation of Lorentz symmetry. We identify the Yb+ ion as an ideal system with high sensitivity as well as excellent experimental controllability. By applying quantum information inspired technology to Yb+, we expect tests of LLI violating physics in the electron-photon sector to reach levels of $10^{-23}$, five orders of magnitude more sensitive than the current best bounds [5-7]. Most importantly, the projected sensitivity of $10^{-23}$ for the Yb+ ion tests will allow for the first time to probe whether Lorentz violation is minimally suppressed at low energies for photons and electrons.


title: Native point defects in CaS: A focus on doping limit for persistent luminescence

abstract: We studied native point defects in CaS by DFT+ Hubbard U method. The effect of the localization of the d orbitals of Ca pseudopotential has been included. The Hubbard U corrected d-orbital for Ca sites are playing a role assisting the charge transfer from p orbitals of S site to the Ca site, giving both localized electron and hole states within the band gap, with an energy interval of 1.2~1.4 eV. This corresponds to the localized excitonic levels below the conduction band edge for the optical absorption. The electronic properties and formation energies of native point defects have been discussed. We found the neutral S vacancy has the lowest energy of 0.62 eV under Ca-rich limit. The Schottky defect pair defect is another dominant defect with cost of 1.51 eV per defect site from S-rich to Ca-rich chemical potential limits. The defect formation energy further summarized that the doubly positive S vacancy and doubly negative Ca vacancy are both the most stable donor-type and acceptor-type defects respectively coexisting in CaS. We also summarized a narrow doping limit energy which has been determined as 1.33 eV constantly in CaS independent to different chemical potential limits. But such doping allowed range entirely shifts from valence band vicinity toward the conduction band edge from S-rich to Ca-rich limits. Under S-rich limit, the dopable range for EF shifting corresponds to Eu2+ doping experiments, the luminescence wavelength limit is predicted to be 659 nm (1.88 eV) remarkably close to the reported results 650 nm. This gives a solid theoretical reference for the lanthanide ions doping experiments in CaS. We conclude that the defect levels combined with formation energy is an accessible way to suggest the doping energy range and an assistant to explain the activation stage of photostimulated luminescence (PSL) mechanism in lanthanide ions doped crystal materials.


title: Uniform Regularity Estimates in Homogenization Theory of Elliptic System with Lower Order terms

abstract: In this paper, we extend the uniform regularity estimates obtained by M. Avellanda and F. Lin in the paper of Compactness methods in the theory of homogenization (Comm. Pure Appl. Math. 40(1987), no.6, 803-847) to the more general second order elliptic systems in divergences form, with rapidly oscillating periodic coefficients. We establish not only sharp $W^{1,p}$ estimates, Holder estimates, Lipschitz estimates and non-tangential maximal function estimates for the Dirichelt problem on a $C^{1,\eta}$ domain, but also a sharp $O(\varepsilon)$ convergence rate in $H_0^1(\Omega)$ by virtue of the Dirichlet correctors. The well known compactness methods are not employed here, instead we construct the transformations to make full use of the corresponding results developed by M. Avellanda and F. Lin.


title: The Weihai Observatory search for close-in planets orbiting giant stars

abstract: Planets are known to orbit giant stars, yet there is a shortage of planets orbiting within ~0.5 AU (P<100 days). First-ascent giants have not expanded enough to engulf such planets, but tidal forces can bring planets to the surface of the star far beyond the stellar radius. So the question remains: are tidal forces strong enough in these stars to engulf all the missing planets? We describe a high-cadence observational program to obtain precise radial velocities of bright giants from Weihai Observatory of Shandong University. We present data on the planet host Beta Gem (HD 62509), confirming our ability to derive accurate and precise velocities; our data achieve an rms of 7.3 m/s about the Keplerian orbit fit. This planet-search programme currently receives ~100 nights per year, allowing us to aggressively pursue short-period planets to determine whether they are truly absent.


title: A Time-Multiplexed Heralded Single-Photon Source

abstract: Photons have proven to be excellent carriers of quantum information, and play essential roles in numerous quantum information processing (QIP) applications. In particular, heralded single- photon sources via spontaneous parametric-down conversion (SPDC) have been a key technology for demonstrating small-scale QIP, yet their low generation efficiency is a critical limitation for further scaling up optical QIP technology. In order to efficiently overcome the probabilistic nature of SPDC, here we demonstrate time multiplexing for up to 30 time slots of a periodically pumped heralded single-photon source, using a switchable low-loss optical storage cavity. We observe a maximum single-photon probability of 38.6% in periodic output time windows, corresponding to 6 times enhancement over a non-multiplexed source, but with no increase in the contribution of unwanted multi-photon events. Combining this time-multiplexing technique with a heralded source producing pure single-photon states should enable larger scale optical QIP systems than ever realized.


title: On Kernel Mengerian Orientations of Line Multigraphs

abstract: We present a polyhedral description of kernels in orientations of line multigraphs. Given a digraph $D$, let $FK(D)$ denote the fractional kernel polytope defined on $D$, and let ${\sigma}(D)$ denote the linear system defining $FK(D)$. A digraph $D$ is called kernel perfect if every induced subdigraph $D^\prime$ has a kernel, called kernel ideal if $FK(D^\prime)$ is integral for each induced subdigraph $D^\prime$, and called kernel Mengerian if ${\sigma} (D^\prime)$ is TDI for each induced subdigraph $D^\prime$. We show that an orientation of a line multigraph is kernel perfect iff it is kernel ideal iff it is kernel Mengerian. Our result strengthens the theorem of Borodin et al. [3] on kernel perfect digraphs and generalizes the theorem of Kiraly and Pap [7] on stable matching problem.


title: Theory of spin Hall magnetoresistance (SMR) and related phenomena

abstract: We review the recently discovered spin Hall magnetoresistance (SMR) and related effects from a theoretical point of view. The SMR is observed in bilayers of a magnetic insulator and a metal, in which spin currents aregenerated in the normal metal due to the spin Hall effect. The associated angular momentum transfer to the ferromagnetic layer and thereby the electrical resistance is modulated by the angle between the applied current and the magnetization direction. The SMR provides a convenient tool to non-invasively measure the magnetization direction and spin-transfer torque to an insulator. We introduce the minimal theoretical instruments to calculate the SMR, i.e. spin diffusion theory and quantum mechanical boundary conditions. This leads to a small set of parameters that can be fitted to experiments. We discuss the limitations of the theory as well as alternative mechanisms such as the ferromagnetic proximity effect and Rashba spin-orbit torques, and point out new developments related to the SMR.


title: FastGP: An R Package for Gaussian Processes

abstract: Despite their promise and ubiquity, Gaussian processes (GPs) can be difficult to use in practice due to the computational impediments of fitting and sampling from them. Here we discuss a short R package for efficient multivariate normal functions which uses the Rcpp and RcppEigen packages at its core. GPs have properties that allow standard functions to be sped up; as an example we include functionality for Toeplitz matrices whose inverse can be computed in O(n^2) time with methods due to Trench and Durbin (Golub & Van Loan 1996), which is particularly apt when time points (or spatial locations) of a Gaussian process are evenly spaced, since the associated covariance matrix is Toeplitz in this case. Additionally, we include functionality to sample from a latent variable Gaussian process model with elliptical slice sampling (Murray, Adams, & MacKay 2010).


title: A National Effort for Motivating Indian Students and Teachers towards Algorithmic Research

abstract: During 2008-2015, twenty-two introductory workshops on graph and geometric algorithms were organized for teachers and students (undergraduate, post-graduate and doctoral) of engineering colleges and universities at different states and union territories of India. The lectures were meant to provide exposure to the field of graph and geometric algorithms and to motivate the participants towards research. Fifty-eight professors from TIFR, IITs, IISc, IMSc, CMI, ISI Kolkata, and other institutes and universities delivered invited lectures on different topics in the design and analysis of algorithms, discrete applied mathematics, computer graphics, computer vision, and robotics. The first four workshops were funded by TIFR, BRNS and IIT Kharagpur, and the remaining workshops were funded by the NBHM. In this paper, we present the salient features of these workshops, and state our observations on the national impact of these workshops.


title: Spin-triplet superconductivity in a weak-coupling Hubbard model for the quasi-one-dimensional compound Li$_{0.9}$Mo$_6$O$_{17}$

abstract: The purple bronze Li$_{0.9}$Mo$_6$O$_{17}$ is of interest due to its quasi-one-dimensional electronic structure and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, we analyze a minimal Hubbard model for this material involving four molybdenum orbitals per unit cell near quarter filling, using asymptotically exact perturbative renormalization group methods. We find that spin triplet odd-parity superconductivity is the dominant instability. Approximate nesting properties of the two quasi-one-dimensional Fermi surfaces enhance certain second-order processes, which play crucial roles in determining the structure of the pairing gap. Notably, we find that the gap has accidental nodes, i.e. it has more sign changes than required by the point-group symmetry.


title: Theoretical Study of Coulomb Correlations and Spin-Orbit Coupling in SrIrO3

abstract: Given that energy scales associated with crystal field splitting, spin orbit coupling and coulomb correlations in iridates are comparable, hence leading to exotic properties, we investigate the physical properties of orthorhombic SrIrO3 using density functional theory. Our calculations, however, show that SrIrO 3 is a bad metal with no long range magnetic ordering, unlike its sister compounds Sr2IrO4 and Sr3Ir2O7. Moreover, despite having large band width, it appears conclusive that the larger resistivity in SrIrO 3 is due to spin orbit interactions. Besides, the effects of electron-electron correlations on its electronic structure and magnetic properties are also discussed.


title: X-ray and radio emission from Type IIn supernova SN 2010jl

abstract: We present all X-ray and radio observations of the Type IIn supernova SN 2010jl. The X-ray observations cover a period up to day 1500 with Chandra, XMM-Newton, NuSTAR and Swift-XRT. The Chandra observations after 2012 June, the XMM-Newton observation in 2013 November, and most of the Swift-XRT observations until 2014 December are presented for the first time. All the spectra can be fitted by an absorbed hot thermal model except for \chandra spectra on 2011 October and 2012 June when an additional component is needed. Although the origin of this component is uncertain, it is spatially coincident with the supernova and occurs when there are changes to the supernova spectrum in the energy range close to that of the extra component, indicating that the emission is related to the supernova. The X-ray light curve shows an initial plateau followed by a steep drop starting at day $\sim 300$. We attribute the drop to a decrease in the circumstellar density. The column density to the X-ray emission drops rapidly with time, showing that the absorption is in the vicinity of the supernova. We also present Very Large Array radio observations of SN 2010jl. Radio emission was detected from \sn from day 570 onwards. The radio light curves and spectra suggest that the radio luminosity was close to its maximum at the first detection. The velocity of the shocked ejecta derived assuming synchrotron self absorption is much less than that estimated from the optical and X-ray observations, suggesting that free-free absorption dominates.


title: Pulsating fronts for bistable on average reaction-diffusion equations in a time periodic environment

abstract: This paper is devoted to reaction-diffusion equations with bistable nonlinearities depending periodically on time. These equations admit two linearly stable states. However, the reaction terms may not be bistable at every time. These may well be a periodic combination of standard bistable and monostable nonlinearities. We are interested in a particular class of solutions, namely pulsating fronts. We prove the existence of such solutions in the case of small time periods of the nonlinearity and in the case of small perturbations of a nonlinearity for which we know there exist pulsating fronts. We also study uniqueness, monotonicity and stability of pulsating fronts.


title: Synchronization hypothesis in the Winfree model

abstract: We consider $N$ oscillators coupled by a mean field as in the Winfree model. The model is governed by two parameters: the coupling strength $\kappa$ and the spectrum width $\gamma$ of the frequencies of each oscillator. In the uncoupled regime, $\kappa=0$, each oscillator possesses its own natural frequency, and the difference between the phases of any two oscillators grows linearly in time. We say that $N$ oscillators are synchronized if the difference between any two phases is uniformly bounded in time. We identify a new hypothesis for the existence of synchronization. The domain in $(\gamma,\kappa)$ of synchronization contains coupling values that are both weak and strong. Moreover the domain is independent of the number of oscillators and the distribution of the frequencies. We give a numerical counter-example which shows that this hypothesis is necessary for the existence of synchronization.


title: Mott localization in a pure stripe antiferromagnet Rb$_{1-\delta}$Fe$_{1.5-\sigma}$S$_2$

abstract: A combination of neutron diffraction and angle-resolved photoemission spectroscopy measurements on a pure antiferromagnetic stripe Rb$_{1-\delta}$Fe$_{1.5-\sigma}$S$_2$ is reported. A neutron diffraction experiment on a powder sample shows that a 98$\%$ volume fraction of the sample is in the antiferromagnetic stripe phase with rhombic iron vacancy order and a refined composition of Rb$_{0.66}$Fe$_{1.36}$S$_{2}$, and that only 2$\%$ of the sample is in the block antiferromagnetic phase with $\sqrt{5}\times \sqrt{5}$ iron vacancy order. Furthermore, a neutron diffraction experiment on a single crystal shows that there is only a single phase with the stripe antiferromagnetic order with the refined composition of Rb$_{0.78}$Fe$_{1.35}$S$_2$, while the phase with block antiferromagnetic order is absent. Angle-resolved photoemission spectroscopy measurements on the same crystal with the pure stripe phase reveal that the electronic structure is gapped at the Fermi level with a gap larger than 0.325 eV. The data collectively demonstrates that the extra 10$\%$ iron vacancies in addition to the rhombic iron vacancy order effectively impede the formation of the block antiferromagnetic phase; the data also suggest that the stripe antiferromagnetic phase with rhombic iron vacancy order is a Mott insulator.