title: A Makespan Lower Bound for the Scheduling of the Tiled Cholesky Factorization based on ALAP scheduling

abstract: Due to the advent of multicore architectures and massive parallelism, the tiled Cholesky factorization algorithm has recently received plenty of attention and is often referenced by practitioners as a case study. It is also implemented in mainstream dense linear algebra libraries. However, we note that theoretical study of the parallelism of this algorithm is currently lacking. In this paper, we present new theoretical results about the tiled Cholesky factorization in the context of a parallel homogeneous model without communication costs. We use standard flop-based weights for the tasks. For a $t$-by-$t$ matrix, we know that the critical path of the tiled Cholesky algorithm is $9t-10$ and that the weight of all tasks is $t^3$. In this context, we prove that no schedule with less than $0.185 t^2$ processing units can finish in a time less than the critical path. In perspective, a naive bound gives $0.11 t^2.$ We then give a schedule which needs less than $0.25 t^2+0.16t+3$ processing units to complete in the time of the critical path. In perspective, a naive schedule gives $0.50 t^2.$ In addition, given a fixed number of processing units, $p$, we give a lower bound on the execution time as follows: $$\max( \frac{t^{3}}{p}, \frac{t^{3}}{p} - 3\frac{t^2}{p} + 6\sqrt{2p} - 7 , 9t-10).$$ The interest of the latter formula lies in the middle term. Our results stem from the observation that the tiled Cholesky factorization is much better behaved when we schedule it with an ALAP (As Late As Possible) heuristic than an ASAP (As Soon As Possible) heuristic. We also provide scheduling heuristics which match closely the lower bound on execution time. We believe that our theoretical results will help practical scheduling studies. Indeed, our results enable to better characterize the quality of a practical schedule with respect to an optimal schedule.


title: The $\rho \gamma^* \to \pi (\rho)$ transition form factors in the Perturbative QCD factorization approach

abstract: In this paper, we studied the $\rho \gamma^* \to \pi$ and $\rho\gamma^*\to \rho$ transition processes and made the calculations for the $\rho\pi$ transition form factor $Q^4 F_{\rho\pi}(Q^2)$ and the $\rho$ meson electromagnetic form factors, $F_{\rm LL, LT,TT}(Q^2)$ and $F_{1,2,3}(Q^2)$, by employing the perturbative QCD (PQCD) factorization approach. For the $\rho \gamma^* \to \pi$ transition, we found that the contribution to form factor $Q^4 F_{\rho\pi}(Q)$ from the term proportional to the distribution amplitude combination $\phi^T_{\rho}(x_1)\phi^P_{\pi}(x_2)$ is absolutely dominant, and the PQCD predictions for both the size and the $Q^2$-dependence of this form factor $Q^4 F_{\rho\pi}(Q^2)$ agree well with those from the extended ADS/QCD models or the light-cone QCD sum rule. For the $\rho \gamma^* \to \rho$ transition and in the region of $Q^2\geq 3$ GeV$^2$, further more, we found that the PQCD predictions for the magnitude and their $Q^2$-dependence of the $F_1(Q^2)$ and $F_2(Q^2)$ form factors agree well with those from the QCD sum rule, while the PQCD prediction for $F_3(Q^2)$ is much larger than the one from the QCD sum rule.


title: Results of measurements of an environment neutron background at BNO INR RAS objects with the helium proportional counter

abstract: A method of measurements of the environmental neutron background at the Baksan Neutrino Observatory of the INR RAS are described. Measurements were done by using of a proportional counter filled with mixture of Ar(2 at)+$^3$He(4 at). The results obtained at the surface and the underground laboratory of the BNO INR RAS are presented. It is shown that a neutron background in the underground laboratory at the 4900 m w.e. depth is decreased by $\sim 260$ times without any special shield in a comparison with the Earth surface. A neutron flux density in the 5-1323.5~cm air height region is constant within the determination error and equal to $(7.1\pm0.1_{\rm{stat}}\pm0.3_{\rm{syst}})\times10^{-3}$ s$^{-1}\cdot$cm$^{-2}$.


title: The asymptotics of the Struve function ${\bf H}_\nu(z)$ for large complex order and argument

abstract: We re-examine the asymptotic expansion of the Struve function ${\bf H}_\nu(z)$ for large complex values of $\nu$ and $z$ satisfying $|\arg\,\nu|\leq\pi/2$ and $|\arg\,z|<\pi/2$. Watson's analysis covers only the case of $\nu$ and $z$ of the same phase with $\nu/z$ in the intervals $(0,1)$ and $(1,\infty)$. The domains in the complex $\nu/z$-plane where the expansion takes on different forms are obtained.


title: Optimal convergence for adaptive IGA boundary element methods for weakly-singular integral equations

abstract: In a recent work, we analyzed a weighted-residual error estimator for isogeometric boundary element methods in 2D and proposed an adaptive algorithm which steers the local mesh-refinement of the underlying partition as well as the multiplicity of the knots. In the present work, we give a mathematical proof that this algorithm leads to convergence even with optimal algebraic rates. Technical contributions include a novel mesh-size function which also monitors the knot multiplicity as well as inverse estimates for NURBS in fractional-order Sobolev norms.


title: Canonical quantization of electromagnetic field in the presence of nonlinear anisotropic magnetodielectric medium with spatial-temporal dispersion

abstract: Modeling a nonlinear anisotropic magnetodielectric medium with spatial-temporal dispersion by two continuum collections of three dimensional harmonic oscillators, a fully canonical quantization of the electromagnetic field is demonstrated in the presence of such a medium. Some coupling tensors of various ranks are introduced that couple the magnetodielectric medium with the electromagnetic field. The polarization and magnetization fields of the medium are defined in terms of the coupling tensors and the oscillators modeling the medium. The electric and magnetic susceptibility tensors of the medium are obtained in terms of the coupling tensors. It is shown that the electric field satisfy an integral equation in frequency domain. The integral equation is solved by an iteration method and the electric field is found up to an arbitrary accuracy.


title: On the topology of a boolean representable simplicial complex

abstract: It is proved that fundamental groups of boolean representable simplicial complexes are free and the rank is determined by the number and nature of the connected components of their graph of flats for dimension $\geq 2$. In the case of dimension 2, it is shown that boolean representable simplicial complexes have the homotopy type of a wedge of spheres of dimensions 1 and 2. Also in the case of dimension 2, necessary and sufficient conditions for shellability and being sequentially Cohen-Macaulay are determined. Complexity bounds are provided for all the algorithms involved.


title: Canonical quantization of electromagnetic field in the presence of absorbing bi-anisotropic multilayer magnetodielectric media

abstract: A bi-anisotropic magnetodielectric medium is modeled by two independent set of three dimensional harmonic oscillators .A fully canonical quantization of electromagnetic field is demonstrated in the presence of a bi-anisotropic magnetodielectric medium. The electric and magnetic polarization fields of the medium are obtained in terms of the dynamical variable modeling the medium. The Heisenberg equations of the system are solved for a multilayer bi-anisotropic magnetodielectric medium.


title: Multifractal Flexibly Detrended Fluctuation Analysis

abstract: Multifractal time series analysis is a approach that shows the possible complexity of the system. Nowadays, one of the most popular and the best methods for determining multifractal characteristics is Multifractal Detrended Fluctuation Analysis (MFDFA). However, it has some drawback. One of its core elements is detrending of the series. In the classical MFDFA a trend is estimated by fitting a polynomial of degree $m$ where $m=const$. We propose that the degree $m$ of a polynomial was not constant ($m\neq const$) and its selection was ruled by an established criterion. Taking into account the above amendment, we examine the multifractal spectra both for artificial and real-world mono- and the multifractal time series. Unlike classical MFDFA method, obtained singularity spectra almost perfectly reflects the theoretical results and for real time series we observe a significant right side shift of the spectrum.


title: Models for SIMP dark matter and dark photon

abstract: We give a review on the SIMP paradigm and discuss a consistent model for SIMP dark mesons in the context of a dark QCD with flavor symmetry. The $Z'$-portal interaction is introduced being compatible with stable dark mesons and is responsible for making the SIMP dark mesons remain in kinetic equilibrium with the SM during the freeze-out process. The SIMP parameter space of the $Z'$ gauge boson can be probed by future collider and direct detection experiments.


title: Proof of a conjecture on `plateaux' phenomenon of graph Laplacian eigenvalues

abstract: Let $G$ be a simple graph. A pendant path of $G$ is a path such that one of its end vertices has degree $1$, the other end has degree $\ge3$, and all the internal vertices have degree $2$. Let $p_k(G)$ be the number of pendant paths of length $k$ of $G$, and $q_k(G)$ be the number of vertices with degree $\ge3$ which are an end vertex of some pendant paths of length $k$. Motivated by the problem of characterizing dendritic trees, N. Saito and E. Woei conjectured that any graph $G$ has some Laplacian eigenvalue with multiplicity at least $p_k(G)-q_k(G)$. We prove a more general result for both Laplacian and signless Laplacian eigenvalues from which the conjecture follows.


title: Networks, Dynamic Factors, and the Volatility Analysis of High-Dimensional Financial Series

abstract: We consider weighted directed networks for analysing, over the period 2000-2013, the interdependencies between volatilities of a large panel of stocks belonging to the S\&P100 index. In particular, we focus on the so-called {\it Long-Run Variance Decomposition Network} (LVDN), where the nodes are stocks, and the weight associated with edge $(i,j)$ represents the proportion of $h$-step-ahead forecast error variance of variable $i$ accounted for by variable $j$'s innovations. To overcome the curse of dimensionality, we decompose the panel into a component driven by few global, market-wide, factors, and an idiosyncratic one modelled by means of a sparse vector autoregression (VAR) model. Inversion of the VAR together with suitable identification restrictions, produces the estimated network, by means of which we can assess how {\it systemic} each firm is.~Our analysis demonstrates the prominent role of financial firms as sources of contagion, especially during the~2007-2008 crisis.


title: On the uniqueness of the Gauss-Bonnet-Chern formula (after Gilkey-Park-Sekigawa)

abstract: On an oriented Riemannian manifold, the Gauss-Bonnet-Chern formula asserts that the Pfaffian of the metric represents, in de Rham cohomology, the Euler class of the tangent bundle. Hence, if the underlying manifold is compact, the integral of the Pfaffian is a topological invariant; namely, the Euler characteristic of the manifold. In this paper we refine a result originally due to Gilkey that characterizes this formula as the only (non-trivial) integral of a differential invariant that is independent of the underlying metric.


title: Bibliometric-Enhanced Information Retrieval: 3rd International BIR Workshop

abstract: The BIR workshop brings together experts in Bibliometrics and Information Retrieval. While sometimes perceived as rather loosely related, these research areas share various interests and face similar challenges. Our motivation as organizers of the BIR workshop stemmed from a twofold observation. First, both communities only partly overlap, albeit sharing various interests. Second, it will be profitable for both sides to tackle some of the emerging problems that scholars face today when they have to identify relevant and high quality literature in the fast growing number of electronic publications available worldwide. Bibliometric techniques are not yet used widely to enhance retrieval processes in digital libraries, although they offer value-added effects for users. Information professionals working in libraries and archives, however, are increasingly confronted with applying bibliometric techniques in their services. The first BIR workshop in 2014 set the research agenda by introducing each group to the other, illustrating state-of-the-art methods, reporting on current research problems, and brainstorming about common interests. The second workshop in 2015 further elaborated these themes. This third BIR workshop aims to foster a common ground for the incorporation of bibliometric-enhanced services into scholarly search engine interfaces. In particular we will address specific communities, as well as studies on large, cross-domain collections like Mendeley and ResearchGate. This third BIR workshop addresses explicitly both scholarly and industrial researchers.


title: Anderson localization in optical lattices with correlated disorder

abstract: We study the Anderson localization of atomic gases exposed to simple-cubic optical lattices with a superimposed disordered speckle pattern. The two mobility edges in the first band and the corresponding critical filling factors are determined as a function of the disorder strength, ranging from vanishing disorder up to the critical disorder intensity where the two mobility edges merge and the whole band becomes localized. Our theoretical analysis is based both on continuous-space models which take into account the details of the spatial correlation of the speckle pattern, and also on a simplified tight-binding model with an uncorrelated distribution of the on-site energies. The mobility edges are computed via the analysis of the energy-level statistics, and we determine the universal value of the ratio between consecutive level spacings at the mobility edge. We analyze the role of the spatial correlation of the disorder, and we also discuss a qualitative comparison with available experimental data for interacting atomic Fermi gases measured in the moderate interaction regime.


title: Optical and Near-Infrared Polarimetry for a Highly Dormant Comet 209P/LINEAR

abstract: We conducted an optical and near-infrared polarimetric observation of the highly dormant Jupiter-Family Comet, 209P/LINEAR. Because of its low activity, we were able to determine the linear polarization degrees of the coma dust particles and nucleus independently, that is $P_n$=30.3$^{+1.3}_{-0.9}$% at $\alpha$=92.2$^\circ$ and $P_n$=31.0$^{+1.0}_{-0.7}$% at $\alpha$=99.5$^\circ$ for the nucleus, and $P_c$=28.8$^{+0.4}_{-0.4}$% at $\alpha$=92.2$^\circ$ and 29.6$^{+0.3}_{-0.3}$% at $\alpha$=99.5$^\circ$ for the coma. We detected no significant variation in $P$ at the phase angle coverage of 92.2$^\circ$-99.5$^\circ$, which may imply that the obtained polarization degrees are nearly at maximum in the phase-polarization curves. By fitting with an empirical function, we obtained the maximum values of linear polarization degrees $P_\mathrm{max}$=30.8% for the nucleus and $P_\mathrm{max}$=29.6% for the dust coma. The $P_\mathrm{max}$ of the dust coma is consistent with those of dust-rich comets. The low geometric albedo of $P_v$=0.05 was derived from the slope-albedo relationship and was associated with high $P_\mathrm{max}$. We examined $P_\mathrm{max}$-albedo relations between asteroids and 209P, and found that the so-called Umov law seems to be applicable on this cometary surface.


title: Optimal growth trajectories with finite carrying capacity

abstract: We consider the problem of finding optimal strategies that maximize the average growth-rate of multiplicative stochastic processes. For a geometric Brownian motion the problem is solved through the so-called Kelly criterion, according to which the optimal growth rate is achieved by investing a constant given fraction of resources at any step of the dynamics. We generalize these finding to the case of dynamical equations with finite carrying capacity, which can find applications in biology, mathematical ecology, and finance. We formulate the problem in terms of a stochastic process with multiplicative noise and a non-linear drift term that is determined by the specific functional form of carrying capacity. We solve the stochastic equation for two classes of carrying capacity functions (power laws and logarithmic), and in both cases compute optimal trajectories of the control parameter. We further test the validity of our analytical results using numerical simulations.


title: Monic monomial representations I Gorenstein-projective modules

abstract: For a $k$-algebra $A$, a quiver $Q$, and an ideal $I$ of $kQ$ generated by monomial relations, let $\Lambda: = A\otimes_k kQ/I$. We introduce the monic representations of $(Q, I)$ over $A$. We give properties of the structural maps of monic representations, and prove that the category ${\rm mon}(Q, I, A)$ of the monic representations of $(Q, I)$ over $A$ is a resolving subcategory of ${\rm rep}(Q, I, A)$. We introduce the condition ${\rm(G)}$. The main result claims that a $\m$-module is Gorenstein-projective if and only if it is a monic module satisfying ${\rm(G)}$. As consequences, the monic $\m$-modules are exactly the projective $\m$-modules if and only if $A$ is semisimple; and they are exactly the Gorenstein-projective $\m$-modules if and only if $A$ is selfinjective, and if and only if ${\rm mon}(Q, I, A)$ is Frobenius.


title: Chiral phase transition and Schwinger mechanism in a pure electric field

abstract: We systematically study the chiral symmetry breaking and restoration in the presence of a pure electric field in the Nambu--Jona-Lasinio (NJL) model at finite temperature and baryon chemical potential. In addition, we also study the effect of the chiral phase transition on the charged pair production due to the Schwinger mechanism. For these purposes, a general formalism for parallel electric and magnetic fields is developed at finite temperature and chemical potential for the first time. In the pure electric field limit $B\rightarrow0$, we compute the order parameter, the transverse-to-longitudinal ratio of the Goldstone mode velocities, and the Schwinger pair production rate as functions of the electric field. The inverse catalysis effect of the electric field to chiral symmetry breaking is recovered. And the Goldstone mode is find to disperse anisotropically such that the transverse velocity is always smaller than the longitudinal one, especially at nonzero temperature and baryon chemical potential. As expected, the quark-pair production rate is greatly enhanced by the chiral symmetry restoration.


title: Particle-based simulations of steady-state mass transport at high P\'eclet numbers

abstract: Conventional approaches for simulating steady-state distributions of particles under diffusive and advective transport at high P\'eclet numbers involve solving the diffusion and advection equations in at least two dimensions. Here, we present an alternative computational strategy by combining a particle-based rather than a field-based approach with the initialisation of particles in proportion to their flux. This method allows accurate prediction of the steady state and is applicable even at high P\'eclet numbers where traditional particle-based Monte-Carlo methods starting from randomly initialised particle distributions fail. We demonstrate that generating a flux of particles according to a predetermined density and velocity distribution at a single fixed time and initial location allows for accurate simulation of mass transport under flow. Specifically, upon initialisation in proportion to their flux, these particles are propagated individually and detected by summing up their Monte-Carlo trajectories in predefined detection regions. We demonstrate quantitative agreement of the predicted concentration profiles with the results of experiments performed with fluorescent particles in microfluidic channels under continuous flow. This approach is computationally advantageous and readily allows non-trivial initial distributions to be considered. In particular, this method is highly suitable for simulating advective and diffusive transport in microfluidic devices.


title: Perfect independent sets with respect to infinitely many relations

abstract: We prove a result on perfect cliques with respect to countably many G-delta relations on a complete metric space. As an application, we show that a Polish group contains a free subgroup generated by a perfect set as long as it contains any uncountable free subgroup. This answers a recent question of G{\l}\c{a}b and Strobin.


title: Comprehensive indicator comparisons intelligible to non-experts: The case of two SNIP versions

abstract: A framework is proposed for comparing different types of bibliometric indicators, introducing the notion of an Indicator Comparison Report. It provides a comprehensive overview of the main differences and similarities of indicators. The comparison shows both the strong points and the limitations of each of the indicators at stake, rather than over-promoting one indicator and ignoring the benefits of alternative constructs. It focuses on base notions, assumptions, and application contexts, which makes it more intelligible to non-experts. As an illustration, a comparison report is presented for the original and the modified SNIP (Source Normalized Impact per Paper) indicator of journal citation impact.


title: On full text download and citation distributions in scientific-scholarly journals

abstract: A statistical analysis of full text downloads of articles in Elseviers ScienceDirect covering all disciplines reveals large differences in download frequencies, their skewness, and their correlation with Scopus-based citation counts, between disciplines, journals, and document types. Download counts tend to be two orders of magnitude higher and less skewedly distributed than citations. A mathematical model based on the sum of two exponentials does not adequately capture monthly download counts. The degree of correlation at the article level within a journal is similar to that at the journal level in the discipline covered by that journal, suggesting that the differences between journals are to a large extent discipline specific. Despite the fact that in all study journals download and citation counts per article positively correlate, little overlap may exist between the set of articles appearing in the top of the citation distribution and that with the most frequently downloaded ones. Usage and citation leaks, bulk downloading, differences between reader and author populations in a subject field, the type of document or its content, differences in obsolescence patterns between downloads and citations, different functions of reading and citing in the research process, all provide possible explanations of differences between download and citation distributions.


title: Characterization of diagonally dominant H-matrices

abstract: We first show that sufficient conditions for a diagonally dominant matrix to be a nonsingular one (and also an H-matrix), obtained independently by Shivakumar and Chew in 1974, and Farid in 1995, are equivalent. Then we simplify the characterization of diagonally dominant H-matrices obtained by Huang in 1995, and using it prove that the Shivakumar-Chew-Farid sufficient condition for a diagonally dominant matrix to be an H-matrix, is also necessary.


title: Altmetrics as traces of the computerization of the research process

abstract: I propose a broad, multi-dimensional conception of altmetrics, namely as traces of the computerization of the research process. Computerization should be conceived in its broadest sense, including all recent developments in ICT and software, taking place in society as a whole. I distinguish four aspects of the research process: the collection of research data and development of research methods; scientific information processing; communication and organization; and, last but not least, research assessment. I will argue that in each aspect, computerization plays a key role, and metrics are being developed to describe this process. I propose to label the total collection of such metrics as Altmetrics. I seek to provide a theoretical foundation of altmetrics, based on notions developed by Michael Nielsen in his monograph Reinventing Discovery: The New Era of Networked Science. Altmetrics can be conceived as tools for the practical realization of the ethos of science and scholarship in a computerized or digital age.


title: Efficient tensor tomography in fan-beam coordinates

abstract: We propose a thorough analysis of the tensor tomography problem on the Euclidean unit disk parameterized in fan-beam coordinates. This includes, for the inversion of the Radon transform over functions, using another range characterization first appearing in [Pestov-Uhlmann, IMRN 2004] to enforce in a fast way classical moment conditions at all orders. When considering direction-dependent integrands (e.g., tensors), a problem where injectivity no longer holds, we propose a suitable representative (other than the traditionally sought-after solenoidal candidate) to be reconstructed, as well as an efficient procedure to do so. Numerical examples illustrating the method are provided at the end.


title: Periodic transit and variability search with simultaneous systematics filtering: Is it worth it?

abstract: By using subsets of the HATNet and K2 (Kepler two-wheel) Campaign 1 databases, we examine the effectiveness of filtering out systematics from photometric time series while simultaneously searching for periodic signals. We carry out tests to recover simulated sinusoidal and transit signals added to time series with both real and artificial noise. We find that the simple (and more traditional) method that performs correction for systematics first and signal search thereafter, produces higher signal recovery rates on the average, while also being substantially faster than the simultaneous method. Independently of the method of search, once the signal is found, a far less time consuming full-fledged model, incorporating both the signal and systematics, must be employed to recover the correct signal shape. As a by-product of the tests on the K2 data, we find that for longer period sinusoidal signals the detection rate decreases (after an optimum value is reached) as the number of light curves used for systematics filtering increases. The decline of the detection rate is observable in both methods of filtering, albeit the simultaneous method performs better in the regime of relative high template number. We suspect that the observed phenomenon is linked to the increased role of low amplitude intrinsic stellar variability in the space-based data. This assumption is also supported by the substantially higher stability of the detection rates for transit signals against the increase of the template number.


title: Set families with a forbidden pattern

abstract: A balanced pattern of order $2d$ is an element $P \in \{+,-\}^{2d}$, where both signs appear $d$ times. Two sets $A,B \subset [n]$ form $P$-pattern, which we denote by $\operatorname{pat}(A,B) = P$, if $A\triangle B = \{j_1,\ldots ,j_{2d}\}$ with $1\leq j_1<\cdots < j_{2d}\leq n$ and $\{i\in [2d]: P_i = + \} = \{i\in [2d]: j_i \in A \setminus B\}$. We say ${\cal A} \subset {\cal P}[n]$ is $P$-free if $\operatorname{pat}(A,B)\neq P$ for all $A,B \in {\cal A}$. We consider the following extremal question: how large can a family ${\cal A} \subset {\cal P}[n]$ be if ${\cal A} $ is $P$-free? We prove a number of results on the sizes of such families. In particular, we show that for some fixed $c>0$, if $P$ is a $d$-balanced pattern with $d < c \log \log n$ then $|{\cal A}| = o(2^n)$. We then give stronger bounds in the cases when (i) $P$ consists of $d$ $+$ signs, followed by $d$ $-$ signs and (ii) $P$ consists of alternating signs. In both cases, if $d = o(\sqrt n)$ then $|{\cal A} | = o(2^n)$. In the case of (i), this is tight. .


title: Water deuteration and ortho-to-para nuclear spin ratio of H2 in molecular clouds formed via the accumulation of HI gas

abstract: We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.


title: Gyro-viscosity and linear dispersion relations in pair-ion magnetized plasmas

abstract: A fluid theory has been developed by taking account of gyro-viscosity to study wave propagation characteristics in a homogeneous pair-ion magnetized plasma with a cylindrical symmetry. The exact dispersion relations derived by the Hankel-Fourier transformation are shown comparable with those observed in the experiment by Oohara and co-workers. The gyro-viscosity is responsible for the change in propagation characteristics of the ion cyclotron wave from forward to backward by suppressing the effect of the thermal pressure which normally causes the forward nature of dispersion. Although the experiment has been already explained by a kinetic theory by the present authors, the kinetic derivations are so involved because of exact particle orbits in phase space, finite Larmor radius effects, and higher order ion cyclotron resonances. The present fluid theory provides a simple and transparent structure to the dispersion relations since the gyro-viscosity is renormalized into the ion cyclotron frequency which itself indicates the backward nature of dispersion. The usual disadvantage of a fluid theory, which treats only fundamental modes of eigen-waves excited in a system and is not able to describe higher harmonics that a kinetic theory does, is compensated by simple derivations and clear picture based on the renormalization of the gyro-viscosity.


title: Integrality Gaps and Approximation Algorithms for Dispersers and Bipartite Expanders

abstract: We study the problem of approximating the quality of a disperser. A bipartite graph $G$ on $([N],[M])$ is a $(\rho N,(1-\delta)M)$-disperser if for any subset $S\subseteq [N]$ of size $\rho N$, the neighbor set $\Gamma(S)$ contains at least $(1-\delta)M$ distinct vertices. Our main results are strong integrality gaps in the Lasserre hierarchy and an approximation algorithm for dispersers. \begin{enumerate} \item For any $\alpha>0$, $\delta>0$, and a random bipartite graph $G$ with left degree $D=O(\log N)$, we prove that the Lasserre hierarchy cannot distinguish whether $G$ is an $(N^{\alpha},(1-\delta)M)$-disperser or not an $(N^{1-\alpha},\delta M)$-disperser. \item For any $\rho>0$, we prove that there exist infinitely many constants $d$ such that the Lasserre hierarchy cannot distinguish whether a random bipartite graph $G$ with right degree $d$ is a $(\rho N, (1-(1-\rho)^d)M)$-disperser or not a $(\rho N, (1-\Omega(\frac{1-\rho}{\rho d + 1-\rho}))M)$-disperser. We also provide an efficient algorithm to find a subset of size exact $\rho N$ that has an approximation ratio matching the integrality gap within an extra loss of $\frac{\min\{\frac{\rho}{1-\rho},\frac{1-\rho}{\rho}\}}{\log d}$. \end{enumerate} Our method gives an integrality gap in the Lasserre hierarchy for bipartite expanders with left degree~$D$. $G$ on $([N],[M])$ is a $(\rho N,a)$-expander if for any subset $S\subseteq [N]$ of size $\rho N$, the neighbor set $\Gamma(S)$ contains at least $a \cdot \rho N$ distinct vertices. We prove that for any constant $\epsilon>0$, there exist constants $\epsilon'<\epsilon,\rho,$ and $D$ such that the Lasserre hierarchy cannot distinguish whether a bipartite graph on $([N],[M])$ with left degree $D$ is a $(\rho N, (1-\epsilon')D)$-expander or not a $(\rho N, (1-\epsilon)D)$-expander.


title: Integral Images: Efficient Algorithms for Their Computation and Storage in Resource-Constrained Embedded Vision Systems

abstract: The integral image, an intermediate image representation, has found extensive use in multi-scale local feature detection algorithms, such as Speeded-Up Robust Features (SURF), allowing fast computation of rectangular features at constant speed, independent of filter size. For resource-constrained real-time embedded vision systems, computation and storage of integral image presents several design challenges due to strict timing and hardware limitations. Although calculation of the integral image only consists of simple addition operations, the total number of operations is large owing to the generally large size of image data. Recursive equations allow substantial decrease in the number of operations but require calculation in a serial fashion. This paper presents two new hardware algorithms that are based on the decomposition of these recursive equations, allowing calculation of up to four integral image values in a row-parallel way without significantly increasing the number of operations. An efficient design strategy is also proposed for a parallel integral image computation unit to reduce the size of the required internal memory (nearly 35% for common HD video). Addressing the storage problem of integral image in embedded vision systems, the paper presents two algorithms which allow substantial decrease (at least 44.44%) in the memory requirements. Finally, the paper provides a case study that highlights the utility of the proposed architectures in embedded vision systems.


title: Dehn twists on Kauffman bracket skein algebras

abstract: We give an explicit formula for the action of the Dehn twist along a simple closed curve in a compact connected oriented surface on the completion of the filtered skein modules. To do this, we introduce filtrations of the Kauffman bracket skein algebra and the Kauffman bracket skein modules on the surface.


title: Proactive Eavesdropping via Jamming for Rate Maximization over Rayleigh Fading Channels

abstract: Instead of against eavesdropping, this letter proposes a new paradigm in wireless security by studying how a legitimate monitor (e.g., government agencies) efficiently eavesdrops a suspicious wireless communication link. The suspicious transmitter controls its communication rate over Rayleigh fading channels to maintain a target outage probability at the receiver, and the legitimate monitor can successfully eavesdrop only when its achievable rate is no smaller than the suspicious communication rate. We propose a proactive eavesdropping via jamming approach to maximize the average eavesdropping rate, where the legitimate monitor sends jamming signals with optimized power control to moderate the suspicious communication rate.


title: Molecular Cloud Fragmentation and Core Collapse

abstract: I review some steps in the conversion of molecular cloud gas into stars and planets, with an emphasis in this presentation on the early stage molecular cloud fragmentation that leads to elongated filaments/ribbons. Magnetic fields can play a crucial role in all stages and need to be invoked particularly for early stage fragmentation as well as in late core collapse where it may control disk formation. I also review some elements of hydrodynamic modeling of disk evolution.


title: Memory-Efficient Design Strategy for a Parallel Embedded Integral Image Computation Engine

abstract: In embedded vision systems, parallel computation of the integral image presents several design challenges in terms of hardware resources, speed and power consumption. Although recursive equations significantly reduce the number of operations for computing the integral image, the required internal memory becomes prohibitively large for an embedded integral image computation engine for increasing image sizes. With the objective of achieving high-throughput with minimum hardware resources, this paper proposes a memory-efficient design strategy for a parallel embedded integral image computation engine. Results show that the design achieves nearly 35% reduction in memory for common HD video.


title: Alfven wave coupled with flow-driven fluid instability in interpenetrating plasmas

abstract: The Alfven wave is analyzed in case of one quasineutral plasma propagating with some constant speed $v_0$ through another static quasineutral plasma. A dispersion equation is derived describing the Alfven wave coupled with the flow driven mode $\omega= k v_0$ and solutions are discussed analytically and numerically. The usual solutions for two oppositely propagating Alfv\'en waves are substantially modified due to the flowing plasma. More profound is modification of the solution propagating in the negative direction with respect to the magnetic field and the plasma flow. For a large enough flow speed (exceeding the Alfven speed in the static plasma), this negative solution may become non-propagating, with frequency equal to zero. In this case it represents a spatial variation of the electromagnetic field. For greater flow speed it becomes a forward mode, and it may merge with the positive one. This merging of the two modes represents the starting point for a flow-driven instability, with two complex-conjugate solutions. The Alfven wave in interpenetrating plasmas is thus modified and coupled with the flow-driven mode and this coupled mode is shown to be growing when the flow speed is large enough. The energy for the instability is macroscopic kinetic energy of the flowing plasma. The dynamics of plasma particles caused by such a coupled wave still remains similar to the ordinary Alfven wave. This means that well-known stochastic heating by the Alfv\'en wave may work, and this should additionally support the potential role of the Alfven wave in the coronal heating.


title: Overlapping group logistic regression with applications to genetic pathway selection

abstract: Discovering important genes that account for the phenotype of interest has long been challenging in genomewide expression analysis. Analyses such as Gene Set Enrichment Analysis (GSEA) that incorporate pathway information have become widespread in hypothesis testing, but pathway-based approaches have been largely absent from regression methods due to the challenges of dealing with overlapping pathways and the resulting lack of available software. The R package grpreg is widely used to fit group lasso and other group-penalized regression models; in this study, we develop an extension, grpregOverlap, to allow for overlapping group structure using the latent variable approach proposed by Jacob et al. (2009). We compare this approach to the ordinary lasso and to GSEA using both simulated and real data. We find that incorporation of prior pathway information substantially improves the accuracy of gene expression classifiers, and we shed light on several ways in which hypothesis-testing approaches such as GSEA differ from regression approaches with respect to the analysis of pathway data.


title: Rapid Online Analysis of Local Feature Detectors and Their Complementarity

abstract: A vision system that can assess its own performance and take appropriate actions online to maximize its effectiveness would be a step towards achieving the long-cherished goal of imitating humans. This paper proposes a method for performing an online performance analysis of local feature detectors, the primary stage of many practical vision systems. It advocates the spatial distribution of local image features as a good performance indicator and presents a metric that can be calculated rapidly, concurs with human visual assessments and is complementary to existing offline measures such as repeatability. The metric is shown to provide a measure of complementarity for combinations of detectors, correctly reflecting the underlying principles of individual detectors. Qualitative results on well-established datasets for several state-of-the-art detectors are presented based on the proposed measure. Using a hypothesis testing approach and a newly-acquired, larger image database, statistically-significant performance differences are identified. Different detector pairs and triplets are examined quantitatively and the results provide a useful guideline for combining detectors in applications that require a reasonable spatial distribution of image features. A principled framework for combining feature detectors in these applications is also presented. Timing results reveal the potential of the metric for online applications.


title: Positivity of Intersection Multiplicity Over a Two-Dimensional Base

abstract: We show that Serre's Intersection Multiplicity Conjecture holds for a formal power series ring A over a complete, two-dimensional regular local ring R. From this, we deduce the corresponding result for the local rings of any scheme X which is a smooth extension of a regular, two-dimensional base Y. We also investigate the connection between intersection multiplicity and transversality in the unramified setting via a local analysis on the blowup.


title: The Curie-Weiss Model of SOC in Higher Dimension

abstract: We build and study a multidimensional version of the Curie-Weiss model of self-organized criticality we have designed in arXiv:1301.6911. For symmetric distributions satisfying some integrability condition, we prove that the sum $S_n$ of the randoms vectors in the model has a typical critical behaviour. The fluctuations are of order $n^{3/4}$ and the limiting law has a density proportional to the exponential of a fourth-degree polynomial.


title: $a_{0}(980)$ as a companion pole of $a_{0}(1450)$

abstract: Light scalar hadrons are often understood as dynamically generated resonances. These arise as `companion poles' in the propagators of $q\bar{q}$ seed states when accounting for meson-loop contributions to the self-energies of the latter. Following this idea, we demonstrate that for the scalar--isovector state $a_{0}(1450)$ the full one-loop propagator has two poles: a pole of the seed state $a_{0}(1450)$ and a companion pole corresponding to $a_{0}(980)$. The positions of these poles are studied by varying the relative coupling strength between the non-derivative and derivative parts of the interactions.


title: Robust Non-linear Wiener-Granger Causality For Large High-dimensional Data

abstract: Wiener-Granger causality is a widely used framework of causal analysis for temporally resolved events. We introduce a new measure of Wiener-Granger causality based on kernelization of partial canonical correlation analysis with specific advantages in the context of large high-dimensional data. The introduced measure is able to detect non-linear and non-monotonous signals, is designed to be immune to noise, and offers tunability in terms of computational complexity in its estimations. Furthermore, we show that, under specified conditions, the introduced measure can be regarded as an estimate of conditional mutual information (transfer entropy). The functionality of this measure is assessed using comparative simulations where it outperforms other existing methods. The paper is concluded with an application to climatological data.


title: Dynamical measurement of the stellar surface density of face-on galaxies

abstract: The DiskMass survey recently provided measurements of the vertical velocity dispersions of disk stars in a sample of nearly face-on galaxies. By setting the disk scale-heights to be equal to those of edge-on galaxies with similar scale-lengths, it was found that these disks must be sub-maximal, with surprisingly low K-band mass-to-light ratios of the order of $M_\star/L_K \simeq 0.3 M_\odot/L_\odot$. This study made use of a simple relation between the disk surface density and the measured velocity dispersion and scale height of the disk, neglecting the shape of the rotation curve and the dark matter contribution to the vertical force, which can be especially important in the case of sub-maximal disks. Here, we point out that these simplifying assumptions led to an overestimation of the stellar mass-to-light ratios. Relaxing these assumptions, we compute even lower values than previously reported for the mass-to-light ratios, with a median $M_\star/L_K \simeq 0.18 M_\odot/L_\odot$, where 14 galaxies have $M_\star/L_K < 0.11$. Invoking prolate dark matter halos made only a small difference to the derived $M_\star/L_K$, although extreme prolate halos ($q>1.5$ for the axis ratios of the potential) might help. The cross-terms in the Jeans equation are also generally negligible. These deduced K-band stellar mass-to-light ratios are even more difficult to reconcile with stellar population synthesis models than the previously reported ones.


title: Strong hypercontractivity and logarithmic Sobolev inequalities on stratified complex Lie groups

abstract: We show that for a hypoelliptic Dirichlet form operator A on a stratified complex Lie group, if the logarithmic Sobolev inequality holds, then a holomorphic projection of A is strongly hypercontractive in the sense of Janson. This extends previous results of Gross to a setting in which the operator A is not holomorphic.


title: Modeling and Simulation for Fluid-Rotating Structure Interaction

abstract: In this paper, we study a dynamic fluid-structure interaction (FSI) model for an elastic structure that is immersed and spinning in the fluid. We develop a linear constitutive model to describe the motion of a rotational elastic structure which is suitable for the application of arbitrary Lagrangian-Eulerian (ALE) method in FSI simulation. Additionally, a novel ALE mapping method is designed to generate the moving fluid mesh while the deformable structure spins in a non-axisymmetric fluid channel. The structure velocity is adopted as the principle unknown to form a monolithic saddle-point system together with fluid velocity and pressure. We discretize the nonlinear saddle-point system with mixed finite element method and Newton's linearization, and prove that the derived saddle-point problem is well-posed. The developed methodology is applied to a self-defined elastic structure and a realistic hydro-turbine under a prescribed angular velocity. Both illustrate the satisfactory numerical results of an elastic structure that is deforming and rotating while interacting with the fluid. The numerical validation is also conducted to demonstrate the modeling consistency.


title: On the Principal Permanent Rank Characteristic Sequences of Graphs and Digraphs

abstract: The principal permanent rank characteristic sequence is a binary sequence $r_0 r_1 \ldots r_n$ where $r_k = 1$ if there exists a principal square submatrix of size $k$ with nonzero permanent and $r_k = 0$ otherwise, and $r_0 = 1$ if there is a zero diagonal entry. A characterization is provided for all principal permanent rank sequences obtainable by the family of nonnegative matrices as well as the family of nonnegative symmetric matrices. Constructions for all realizable sequences are provided. Results for skew-symmetric matrices are also included.


title: A Historical Analysis of the Field of OR/MS using Topic Models

abstract: This study investigates the content of the published scientific literature in the fields of operations research and management science (OR/MS) since the early 1950s. Our study is based on 80,757 published journal abstracts from 37 of the leading OR/MS journals. We have developed a topic model, using Latent Dirichlet Allocation (LDA), and extend this analysis to reveal the temporal dynamics of the field, journals, and topics. Our analysis shows the generality or specificity of each of the journals, and we identify groups of journals with similar content, which are both consistent and inconsistent with intuition. We also show how journals have become more or less unique in their scope. A more detailed analysis of each journals' topics over time shows significant temporal dynamics, especially for journals with niche content. This study presents an observational, yet objective, view of the published literature from OR/MS that would be of interest to authors, editors, journals, and publishers. Furthermore, this work can be used by new entrants to the fields of OR/MS to understand the content landscape, as a starting point for discussions and inquiry of the field at large, or as a model for other fields to perform similar analyses.


title: Limit lamination theorem for H-disks

abstract: In this paper we prove a theorem concerning lamination limits of sequences of compact disks $M_n$ embedded in $\mathbb{R}^3$ with constant mean curvature $H_n$, when the boundaries of these disks tend to infinity. This theorem generalizes to the non-zero constant mean curvature case Theorem 0.1 by Colding and Minicozzi in [8]. We apply this theorem to prove the existence of a chord arc result for compact disks embedded in $\mathbb{R}^3$ with constant mean curvature; this chord arc result generalizes Theorem 0.5 by Colding and Minicozzi in [9] for minimal disks.


title: Assessing The Performance Bounds Of Local Feature Detectors: Taking Inspiration From Electronics Design Practices

abstract: Since local feature detection has been one of the most active research areas in computer vision, a large number of detectors have been proposed. This has rendered the task of characterizing the performance of various feature detection methods an important issue in vision research. Inspired by the good practices of electronic system design, a generic framework based on the improved repeatability measure is presented in this paper that allows assessment of the upper and lower bounds of detector performance in an effort to design more reliable and effective vision systems. This framework is then employed to establish operating and guarantee regions for several state-of-the art detectors for JPEG compression and uniform light changes. The results are obtained using a newly acquired, large image database (15092 images) with 539 different scenes. These results provide new insights into the behavior of detectors and are also useful from the vision systems design perspective.


title: Performance Characterization of Image Feature Detectors in Relation to the Scene Content Utilizing a Large Image Database

abstract: Selecting the most suitable local invariant feature detector for a particular application has rendered the task of evaluating feature detectors a critical issue in vision research. No state-of-the-art image feature detector works satisfactorily under all types of image transformations. Although the literature offers a variety of comparison works focusing on performance evaluation of image feature detectors under several types of image transformation, the influence of the scene content on the performance of local feature detectors has received little attention so far. This paper aims to bridge this gap with a new framework for determining the type of scenes, which maximize and minimize the performance of detectors in terms of repeatability rate. Several state-of-the-art feature detectors have been assessed utilizing a large database of 12936 images generated by applying uniform light and blur changes to 539 scenes captured from the real world. The results obtained provide new insights into the behaviour of feature detectors.


title: Distributions of Long-Lived Radioactive Nuclei Provided by Star Forming Environments

abstract: Radioactive nuclei play an important role in planetary evolution by providing an internal heat source, which affects planetary structure and helps facilitate plate tectonics. A minimum level of nuclear activity is thought to be necessary --- but not sufficient --- for planets to be habitable. Extending previous work that focused on short-lived nuclei, this paper considers the delivery of long-lived radioactive nuclei to circumstellar disks in star forming regions. Although the long-lived nuclear species are always present, their abundances can be enhanced through multiple mechanisms. Most stars form in embedded cluster environments, so that disks can be enriched directly by intercepting ejecta from supernovae within the birth clusters. In addition, molecular clouds often provide multiple episodes of star formation, so that nuclear abundances can accumulate within the cloud; subsequent generations of stars can thus receive elevated levels of radioactive nuclei through this distributed enrichment scenario. This paper calculates the distribution of additional enrichment for $^{40}$K, the most abundant of the long-lived radioactive nuclei. We find that distributed enrichment is more effective than direct enrichment. For the latter mechanism, ideal conditions lead to about 1 in 200 solar systems being directly enriched in $^{40}$K at the level inferred for the early solar nebula (thereby doubling the abundance). For distributed enrichment from adjacent clusters, about 1 in 80 solar systems are enriched at the same level. Distributed enrichment over the entire molecular cloud is more uncertain, but can be even more effective.


title: BSM Lessons from the SM Higgs

abstract: In this talk I will review the implications on Standard Model (SM) and Beyond the SM (BSM) theory of the experimental exploration of the scalar sector. Given that the Higgs discovery has been the most important achievement, I will start with a general overview of its dramatic impact on the physics of fundamental interaction. Next, I will describe how scalar sector results, including the measurements of the Higgs couplings and mass and the negative searches for extra scalars, constrain concrete BSM scenarios related with Naturalness. In particular, I will discuss composite Higgs and low energy supersymmetry, both in its minimal and next-to-minimal implementations.


title: Periods of magnetic field variations in the Ap star $\gamma$ Equulei (HD 201601)

abstract: This paper presents a series of 95 new measurements of the longitudinal (effective) magnetic field $B_e$ of the Ap star $\gamma$ Equ (HD 201601). Observations were obtained at the coud\'e focus of the 1-m reflector at the Special Astrophysical Observatory (SAO RAS) in Russia over a time period of 4190 days (more than 11 years). We compiled a long record of $B_e$ points, adding our measurements to all published data. The time series of magnetic data consists of 395 $B_e$ points extending for 24488 days, or over 67 years. Various methods of period determination were examined for the case in which the length of the observed time series is rather short and amounts only to ~69 percent of the period. We argue that the fitting of a sine wave to the observed $B_e$ points by least squares yields the most reliable period in the case of $\gamma$ Equ. Therefore, the best period for long-term magnetic variations of $\gamma$ Equ, and hence the rotational period, is $P_{\rm rot}=35462.5 \pm 1149$ days $= 97.16 \pm 3.15$ years.


title: Constraints on the neutrino parameters by future cosmological 21cm line and precise CMB polarization observations (PhD thesis, The Graduate University for Advanced Studies (SOKENDAI))

abstract: Observations of the 21 cm line radiation coming from the epoch of reionization have a great capacity to study the cosmological growth of the Universe. Also, CMB polarization produced by gravitational lensing has a large amount of information about the growth of matter fluctuations at late time. In this thesis, we investigate their sensitivities to the impact of neutrino property on the growth of density fluctuations, such as the total neutrino mass, the neutrino mass hierarchy, the effective number of neutrino species (extra radiation), and the lepton asymmetry of our Universe. We will show that by combining the precise CMB polarization observations with Square Kilometer Array (SKA) we can measure the impact of non-zero neutrino mass on the growth of density fluctuation, and determine the neutrino mass hierarchy at 2 sigma level if the total neutrino mass is smaller than 0.1 eV. Additionally, we will show that by using these combinations we can constrain the lepton asymmetry better than big-bang nucleosynthesis (BBN). Besides we discuss constraints on that in the presence of some extra radiation, and show that the 21 cm line observations can substantially improve the constraints obtained by CMB alone, and allow us to distinguish the effects of the lepton asymmetry from those of extra radiation.


title: Satisfiability of Short Circuit Logic

abstract: The logical connectives typically found in programming languages are similar to their mathematical counterparts, yet different due to their short-circuit behaviour -- when evaluating them, the second argument is only evaluated if the first argument is not sufficient to determine the result. Combined with the possibility of side-effects, this creates a different type of logic called Short Circuit Logic. A greater theoretical understanding of this logic can lead to more efficient programming and faster program execution. In this thesis, formula satisfiability in the context of Short Circuit Logic is discussed. A formal definition of evaluation based on valuation algebras is presented, alongside an alternative definition based on valuation paths. The accompanying satisfiability and `path-satisfiability' are then proven to be equivalent, and an implementation of path-satisfiability is given. Although five types of valuation algebras can be discerned, there are only three corresponding types of valuation paths. From this, conclusions are drawn about satisfiability and side-effects; the manner in which side-effects alter truth values is relevant when analysing satisfiability, but the side-effects themselves are not.


title: Statistical theory of light nucleus reaction and application to $^9$Be(p, xn) reaction

abstract: A statistical theory of light nucleus reaction (STLN) is proposed to describe both neutron and light charged particle induced nuclear reactions with 1p-shell light nuclei involved. The dynamic of STLN is described by the unified Hauser-Feshbach and exciton model, of which the angular momentum and parity conservations are considered in equilibrium and pre-equilibrium processes. The Coulomb barriers of the incident and outgoing charged particles, which seriously influence the open reaction channels, could be reasonably considered in the incident channel and the different outgoing channels. In kinematics, the recoiling effects in various emission processes are taken strictly into account. Taking $^9$Be(p, xn) reaction as an example, we calculate the double-differential cross sections of outgoing neutrons and charged particles using PUNF code in the frame of STLN. The calculated results agree very well with the existing experimental neutron double-differential cross sections at $E_p=18$ MeV, and indicate that PUNF code is a powerful tool to set up file-6 in the reaction data library for the light charged particle induced nuclear reactions with 1p-shell light nuclei involved.


title: The Role of the Pauli-Lubanski Vector for the Dirac, Weyl, Proca, Maxwell, and Fierz-Pauli Equations

abstract: We analyze basic relativistic wave equations for the classical fields, such as Dirac's equation, Weyl's two-component equation for massless neutrinos, and the Proca, Maxwell, and Fierz-Pauli equations, from the viewpoint of the Pauli-Lubanski vector and the Casimir operators of the Poincare group. In general, in this group-theoretical approach, the above wave equations arise in certain overdetermined forms, which can be reduced to the conventional ones by a Gaussian elimination. A connection between the spin of a particle/field and consistency of the corresponding overdetermined system is emphasized in the massless case.


title: Mg line formation in late-type stellar atmospheres: II. Calculations in a grid of 1D models

abstract: Mg is the alpha element of choice for Galactic population and chemical evolution studies as it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from LTE need to be accounted for. Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in late-type stars. These can be used e.g., to correct LTE spectra and abundances. Using the model atom built and tested in the preceding paper in this series, we performed non-LTE radiative transfer calculations in a grid of 3945 stellar 1D atmospheric models. We used a sub-grid of 86 models to explore the propagation of errors in the recent atomic collision calculations to the radiative transfer results. We obtained departure coefficients for all the levels and equivalent widths (in LTE and non-LTE) for all the radiative transitions included in the "final" model atom of Osorio et al.. We present and describe our results and show some examples of applications of the data. The errors due to uncertainties in the collisional data are investigated and tabulated. The results for equivalent widths and departure coefficients are made freely available. Giants tend to have negative abundance corrections while dwarfs have positive, though small, corrections. Error analysis results show that uncertainties related to the atomic collision data are typically of order 0.01 dex or less, although for few stellar models in specific lines uncertainties can be as large as 0.03 dex. As these errors are less than or of the same order as typical corrections, we expect that we can use these results to extract Mg abundances from high quality spectra more reliably than from classical LTE analysis.


title: Self-Consistent Sources for Integrable Equations via Deformations of Binary Darboux Transformations

abstract: We reveal the origin and structure of self-consistent source extensions of integrable equations from the perspective of binary Darboux transformations. They arise via a deformation of the potential that is central in this method. As examples, we obtain in particular matrix versions of self-consistent source extensions of the sine-Gordon, nonlinear Schrodinger, KdV, Boussinesq, KP, Davey-Stewartson, two-dimensional Toda lattice and discrete KP systems. We also recover a (2+1)-dimensional version of the Yajima-Oikawa system from a deformation of the pKP hierarchy. By construction, these systems are accompanied by a hetero binary Darboux transformation, which generates solutions of such a system from a solution of the source-free system and additionally solutions of an associated linear system and its adjoint. The essence of all this is encoded in universal equations in the framework of bidifferential calculus.


title: Viral epidemiology of the adult Apis Mellifera infested by the Varroa destructor mite

abstract: The ectoparasitic mite Varroa destructor has become one of the major worldwide threats for apiculture. Varroa destructor attacks the honey bee Apis mellifera weakening its host by sucking hemolymph. However, the damage to bee colonies is not strictly related to the parasitic action of the mite but it derives, above all, from its action as vector increasing the trasmission of many viral diseases such as acute paralysis (ABPV) and deformed wing viruses (DWV), that are considered among the main causes of CCD (Colony Collapse Disorder). In this work we discuss an SI model that describes how the presence of the mite affects the epidemiology of these viruses on adult bees. We characterize the system behavior, establishing that ultimately either only healthy bees survive, or the disease becomes endemic and mites are wiped out. Another dangerous alternative is the Varroa invasion scenario with the extinction of healthy bees. The final possible configuration is the coexistence equilibrium in which honey bees share their infected hive with mites. The analysis is in line with some observed facts in natural honey bee colonies. Namely, these diseases are endemic. Further, if the mite population is present, necessarily the viral infection occurs. The findings of this study indicate that a low horizontal transmission rate of the virus among honey bees in beehives will help in protecting bee colonies from Varroa infestation and viral epidemics.


title: Digitally reconfigurable complex 2D Dual lattice structure by optical phase engineering

abstract: We present a method to combine two periodic lattice wave-fields to generate a complex dual lattice wave-field which could be employed for microfabrication of corresponding 2D dual lattice structures. Since the addition of two periodic lattice wave-fields is coherent in nature the resultant dual lattice structure is highly dependent on the relative phase difference between constituent wave-fields. We show that it is possible to have control over the dual lattice pattern by precisely controlling this relative phase difference. This control is enabled by making use of digitally addressable phase only spatial light modulator (SLM). We provide the computational method for calculation of corresponding phase mask to be displayed on the SLM and also verify the results experimentally by employing a simple 4f Fourier filter based geometry. The method is completely scalable and reconfigurable in terms of the choice of periodic lattice wave-fields and has the potential to form gradient phase masks which could be useful for fabrication of graded index optical components.


title: Distributed saddle-point subgradient algorithms with Laplacian averaging

abstract: We present distributed subgradient methods for min-max problems with agreement constraints on a subset of the arguments of both the convex and concave parts. Applications include constrained minimization problems where each constraint is a sum of convex functions in the local variables of the agents. In the latter case, the proposed algorithm reduces to primal-dual updates using local subgradients and Laplacian averaging on local copies of the multipliers associated to the global constraints. For the case of general convex-concave saddle-point problems, our analysis establishes the convergence of the running time-averages of the local estimates to a saddle point under periodic connectivity of the communication digraphs. Specifically, choosing the gradient step-sizes in a suitable way, we show that the evaluation error is proportional to $1/\sqrt{t}$, where $t$ is the iteration step. We illustrate our results in simulation for an optimization scenario with nonlinear constraints coupling the decisions of agents that cannot communicate directly.


title: Loss-Gain Equalized Reconfigurable Phaser for Dynamic Radio Analog Signal Processing (R-ASP)

abstract: We present a loss-gain equalized reconfigurable phaser for dynamic radio analog signal processing (R-ASP). Such a phaser provides real-time tunable group delay response with all-pass transmission. We propose a lumped loss-gain implementation, where tuning and equalization are mostly easily achieved. A theoretical study derives the transfer function and the fundamental characteristics of the device. The phaser is finally experimentally demonstrated, first using a single loss-gain pair and finally a three cascaded loss-gain pair structure with full reconfigurability , where up-chirp and down-chirp group delays are shown for illustration. It is expected that this phaser will find wide applications in radio analog signal processing (R-ASP) systems requiring dynamic adaptability.


title: Interplay of curvature-induced micro- and nanodomain structures in multicomponent lipid bilayers

abstract: We discuss different mechanisms for curvature-induced domain formation in multicomponent lipid membranes and present a theoretical model that allows us to study the interplay between the domains. The model represents the membrane by two coupled monolayers, which each carry an additional order parameter field describing the local lipid composition. The spontaneous curvature of each monolayer is coupled to the local composition, moreover, the lipid compositions on opposing monolayers are coupled to each other. Using this model, we calculate the phase behavior of the bilayer in mean-field approximation. The resulting phase diagrams are surprisingly complex and reveal a variety of phases and phase transitions, including a decorated microdomain phase where nanodomains are aligned along the microdomain boundaries. Our results suggest that external membrane tension can be used to control the lateral organization of nanodomains (which might be associated with lipid "rafts") in a multicomponent lipid bilayer.


title: Second-order perturbation theory: problems on large scales

abstract: In general-relativistic perturbation theory, a point mass accelerates away from geodesic motion due to its gravitational self-force. Because the self-force is small, one can often approximate the motion as geodesic. However, it is well known that self-force effects accumulate over time, making the geodesic approximation fail on long timescales. It is less well known that this failure at large times translates to a failure at large distances as well. At second perturbative order, two large-distance pathologies arise: spurious secular growth and infrared-divergent retarded integrals. Both stand in the way of practical computations of second-order self-force effects. Utilizing a simple flat-space scalar toy model, I develop methods to overcome these obstacles. The secular growth is tamed with a multiscale expansion that captures the system's slow evolution. The divergent integrals are eliminated by matching to the correct retarded solution at large distances. I also show how to extract conservative self-force effects by taking local-in-time "snapshots" of the global solution. These methods are readily adaptable to the physically relevant case of a point mass orbiting a black hole.


title: Thermal transport in a 2D stressed nanostructure with mass gradient

abstract: Inspired by some recent molecular dynamics (MD) simulations and experiments on suspended graphene nanoribbons, we study a simplified model where the atoms are disposed in a rectangular lattice coupled by nearest neighbor interactions which are quadratic in the interatomic distance. The system has a mechanical strain, and the border atoms are coupled to Langevin thermal baths. Atom masses vary linearly in the longitudinal direction, modeling an isotope or doping distribution. This asymmetry and tension modify thermal properties. Although the atomic interaction is quadratic, the potential is anharmonic in the coordinates. By direct MD simulations and solving Fokker--Planck equations at low temperatures, we can better understand the role of anharmonicities in thermal rectification. We observe an increasing thermal current with an increasing applied mechanical tension. The temperatures and thermal currents vary along the transverse direction. This effect can be useful to establish which parts of the system are more sensitive to thermal damage. We also study thermal rectification as a function of strain and system size.


title: Coupled electricity and magnetism in solids: multiferroics and beyond

abstract: The interplay of electricity and magnetism, one of the cornerstones of modern physics, takes a special form in solids in such phenomena as magnetoelectricity and the possibility of multiferroic behaviour. In this paper I give a short survey of the main notions of this field, paying special attention to microscopic aspects. Some related phenomena, such as electric activity of magnetic domain walls, etc., are also shortly discussed.


title: Top-k Query Processing on Encrypted Databases with Strong Security Guarantees

abstract: Privacy concerns in outsourced cloud databases have become more and more important recently and many efficient and scalable query processing methods over encrypted data have been proposed. However, there is very limited work on how to securely process top-k ranking queries over encrypted databases in the cloud. In this paper, we focus exactly on this problem: secure and efficient processing of top-k queries over outsourced databases. In particular, we propose the first efficient and provable secure top-k query processing construction that achieves adaptively CQA security. We develop an encrypted data structure called EHL and describe several secure sub-protocols under our security model to answer top-k queries. Furthermore, we optimize our query algorithms for both space and time efficiency. Finally, in the experiments, we empirically analyze our protocol using real world datasets and demonstrate that our construction is efficient and practical.


title: Network Flows that Solve Linear Equations

abstract: We study distributed network flows as solvers in continuous time for the linear algebraic equation $\mathbf{z}=\mathbf{H}\mathbf{y}$. Each node $i$ has access to a row $\mathbf{h}_i^{\rm T}$ of the matrix $\mathbf{H}$ and the corresponding entry $z_i$ in the vector $\mathbf{z}$. The first "consensus + projection" flow under investigation consists of two terms, one from standard consensus dynamics and the other contributing to projection onto each affine subspace specified by the $\mathbf{h}_i$ and $z_i$. The second "projection consensus" flow on the other hand simply replaces the relative state feedback in consensus dynamics with projected relative state feedback. Without dwell-time assumption on switching graphs as well as without positively lower bounded assumption on arc weights, we prove that all node states converge to a common solution of the linear algebraic equation, if there is any. The convergence is global for the "consensus + projection" flow while local for the "projection consensus" flow in the sense that the initial values must lie on the affine subspaces. If the linear equation has no exact solutions, we show that the node states can converge to a ball around the least squares solution whose radius can be made arbitrarily small through selecting a sufficiently large gain for the "consensus + projection" flow under fixed bidirectional graphs. Semi-global convergence to approximate least squares solutions is demonstrated for general switching directed graphs under suitable conditions. It is also shown that the "projection consensus" flow drives the average of the node states to the least squares solution with complete graph. Numerical examples are provided as illustrations of the established results.


title: Remark on the N-barrier method for a class of autonomous elliptic systems

abstract: In this note, we aim to extend the previous work on an N-barrier maximum principle (\cite{hung2015n,hung2015maximum}) to a more general class of systems of two equations. Moreover, an N-barrier maximum principle for systems of three equations is established.


title: Optimizing the Convergence Rate of the Quantum Consensus: A Discrete Time Model

abstract: Motivated by the recent advances in the field of quantum computing, quantum systems are modelled and analyzed as networks of decentralized quantum nodes which employ distributed quantum consensus algorithms for coordination. In the literature, both continuous and discrete time models have been proposed for analyzing these algorithms. This paper aims at optimizing the convergence rate of the discrete time quantum consensus algorithm over a quantum network with $N$ qudits. The induced graphs are categorized in terms of the partitions of integer $N$ by arranging them as the Schreier graphs. It is shown that the original optimization problem reduces to optimizing the Second Largest Eigenvalue Modulus (SLEM) of the weight matrix. Exploiting the Specht module representation of partitions of $N$, the Aldous' conjecture is generalized to all partitions (except ($N$)) in the Hasse diagram of integer $N$. Based on this result, it is shown that the spectral gap of Laplacian of all induced graphs corresponding to partitions (other than ($N$)) of $N$ are the same, while the spectral radius of the Laplacian is obtained from the feasible least dominant partition in the Hasse diagram of integer $N$. The semidefinite programming formulation of the problem is addressed analytically for $N \leq d^2 + 1$ and a wide range of topologies where closed-form expressions for the optimal results are provided. For a quantum network with complete graph topology, solution of the optimization problem based on group association schemes is provided for all values of $N$.


title: Algebraic Conditions for Generating Accurate Adjacency Arrays

abstract: Data processing systems impose multiple views on data as it is processed by the system. These views include spreadsheets, databases, matrices, and graphs. Associative arrays unify and simplify these different approaches into a common two-dimensional view of data. Graph construction, a fundamental operation in the data processing pipeline, is typically done by multiplying the incidence array representations of a graph, $\mathbf{E}_\mathrm{in}$ and $\mathbf{E}_\mathrm{out}$, to produce an adjacency matrix of the graph that can be processed with a variety of machine learning clustering techniques. This work focuses on establishing the mathematical criteria to ensure that the matrix product $\mathbf{E}_\mathrm{out}^\intercal\mathbf{E}_\mathrm{in}$ is the adjacency array of the graph. It will then be shown that these criteria are also necessary and sufficient for the remaining nonzero product of incidence arrays, $\mathbf{E}_\mathrm{in}^\intercal\mathbf{E}_\mathrm{out}$ to be the adjacency matrices of the reversed graph. Algebraic structures that comply with the criteria will be identified and discussed.


title: Stochastic feedback control of quantum transport to realize a dynamical ensemble of two nonorthogonal pure states

abstract: A Markovian open quantum system which relaxes to a unique steady state $\rho_{ss}$ of finite rank can be decomposed into a finite physically realizable ensemble (PRE) of pure states. That is, as shown by Karasik and Wiseman [Phys. Rev. Lett. 106, 020406 (2011)], in principle there is a way to monitor the environment so that in the long time limit the conditional state jumps between a finite number of possible pure states. In this paper we show how to apply this idea to the dynamics of a double quantum dot arising from the feedback control of quantum transport, as previously considered by one of us and co-workers [Phys. Rev. B 84, 085302 (2011)]. Specifically, we consider the limit where the system can be described as a qubit, and show that while the control scheme can always realize a two-state PRE, in the incoherent tunneling regime there are infinitely many PREs compatible with the dynamics that cannot be so realized. For the two-state PREs that are realized, we calculate the counting statistics and see a clear distinction between the coherent and incoherent regimes.


title: Universal Meshes for the Simulation of Brittle Fracture and Moving Boundary Problems

abstract: Universal meshes have recently appeared in the literature as a compu- tationally efficient and robust paradigm for the generation of conforming simpli- cial meshes for domains with evolving boundaries. The main idea behind a univer- sal mesh is to immerse the moving boundary in a background mesh (the universal mesh), and to produce a mesh that conforms to the moving boundary at any given time by adjusting a few of elements of the background mesh. In this manuscript we present the application of universal meshes to the simulation of brittle fracturing. To this extent, we provide a high level description of a crack propagation algorithm and showcase its capabilities. Alongside universal meshes for the simulation of brit- tle fracture, we provide other examples for which universal meshes prove to be a powerful tool, namely fluid flow past moving obstacles. Lastly, we conclude the manuscript with some remarks on the current state of universal meshes and future directions.


title: Sustainability-Aware Cloud Computing Using Virtual Carbon Tax

abstract: In this paper, a solution for sustainable cloud system is proposed and then implemented on a real testbed. The solution composes of optimization of a profit model and introduction of virtual carbon tax to limit environmental footprint of the cloud. The proposed multi-criteria optimizer of the cloud system suggests new optimum CPU frequencies for CPU-cores when the local grid energy mix or the cloud workload changes. The cloud system is implemented on a blade system, and proper middlewares are developed to interact with the blades. The experimental results show that it is possible to significantly decrease the targeted environmental footprint of the system and keep it profitable.


title: Piezoelectric energy harvesting from colored fat-tailed fluctuations: An electronic analogy

abstract: Aiming to optimize piezoelectric energy harvesting from strongly colored fat-tailed fluctuations, we have recently studied the performance of a monostable inertial device under a noise whose statistics depends on a parameter $q$ (bounded for $q<1$, Gaussian for $q=1$, fat-tailed for $q>1$). We have studied the interplay between the potential shape (interpolating between square-well and harmonic-like behaviors) and the noise's statistics and spectrum, and showed that its output power grows as $q$ increases above 1. We now report a real experiment on an electronic analog of the proposed system, which sheds light on its operating principle.


title: Three-Dimensional Network Model for Coupling~of~Fracture and Mass Transport in Quasi-Brittle Geomaterials

abstract: Dual three-dimensional networks of structural and transport elements were combined to model the effect of fracture on mass transport in quasi-brittle geomaterials. Element connectivity of the structural network, representing elasticity and fracture, was defined by the Delaunay tessellation of a random set of points. The connectivity of transport elements within the transport network was defined by the Voronoi tessellation of the same set of points. A new discretisation strategy for domain boundaries was developed to apply boundary conditions for the coupled analyses. The properties of transport elements were chosen to evolve with the crack opening values of neighbouring structural elements. Through benchmark comparisons involving non-stationary transport and fracture, the proposed dual network approach was shown to be objective with respect to element size and orientation.


title: A Local Perspective on Community Structure in Multilayer Networks

abstract: The analysis of multilayer networks is among the most active areas of network science, and there are now several methods to detect dense "communities" of nodes in multilayer networks. One way to define a community is as a set of nodes that trap a diffusion-like dynamical process (usually a random walk) for a long time. In this view, communities are sets of nodes that create bottlenecks to the spreading of a dynamical process on a network. We analyze the local behavior of different random walks on multiplex networks (which are multilayer networks in which different layers correspond to different types of edges) and show that they have very different bottlenecks that hence correspond to rather different notions of what it means for a set of nodes to be a good community. This has direct implications for the behavior of community-detection methods that are based on these random walks.

categories: cs.SI math.PR nlin.AO physics.data-an physics.soc-ph


title: Reheating-era leptogenesis

abstract: We propose a novel leptogenesis scenario at the reheating era. Our setup is minimal in the sense that, in addition to the standard model Lagrangian, we only consider an inflaton and higher dimensional operators. The lepton number asymmetry is produced not by the decay of a heavy particle, but by the scattering between the standard model particles. After the decay of an inflaton, the model is described within the standard model with higher dimensional operators. The Sakharov's three conditions are satisfied by the following way. The violation of the lepton number is realized by the dimension-5 operator. The complex phase comes from the dimension-6 four lepton operator. The universe is out of equilibrium before the reheating is completed. It is found that the successful baryogenesis is realized for the wide range of parameters, the inflaton mass and reheating temperature, depending on the cutoff scale. Since we only rely on the effective Lagrangian, our scenario can be applicable to all mechanisms to generate neutrino Majorana masses.


title: On the behavior of diffusion processes with traps

abstract: We consider processes that coincide with a given diffusion process outside a finite collection of domains. In each of the domains, there is, additionally, a large drift directed towards the interior of the domain. We describe the limiting behavior of the processes as the magnitude of the drift tends to infinity, and thus the domains become trapping, with the time to exit the domains being exponentially large. In particular, in exponential time scales, metastable distributions between the trapping regions are considered.


title: Fra\"iss\'e limits in functional analysis

abstract: We provide a unified approach to Fra\"{\i}ss\'{e} limits in functional analysis, including the Gurarij space, the Poulsen simplex, and their noncommutative analogs. We obtain in this general framework many known and new results about the Gurarij space and the Poulsen simplex, and at the same time establish their noncommutative analogs. Particularly, we construct noncommutative analogs of universal operators in the sense of Rota.


title: Causal Falling Rule Lists

abstract: A causal falling rule list (CFRL) is a sequence of if-then rules that specifies heterogeneous treatment effects, where (i) the order of rules determines the treatment effect subgroup a subject belongs to, and (ii) the treatment effect decreases monotonically down the list. A given CFRL parameterizes a hierarchical bayesian regression model in which the treatment effects are incorporated as parameters, and assumed constant within model-specific subgroups. We formulate the search for the CFRL best supported by the data as a Bayesian model selection problem, where we perform a search over the space of CFRL models, and approximate the evidence for a given CFRL model using standard variational techniques. We apply CFRL to a census wage dataset to identify subgroups of differing wage inequalities between men and women.


title: Monochromatic tree covers and Ramsey numbers for set-coloured graphs

abstract: We consider a generalisation of the classical Ramsey theory setting to a setting where each of the edges of the underlying host graph is coloured with a {\em set} of colours (instead of just one colour). We give bounds for monochromatic tree covers in this setting, both for an underlying complete graph, and an underlying complete bipartite graph. We also discuss a generalisation of Ramsey numbers to our setting and propose some other new directions. Our results for tree covers in complete graphs imply that a stronger version of Ryser's conjecture holds for $k$-intersecting $r$-partite $r$-uniform hypergraphs: they have a transversal of size at most $r-k$. (Similar results have been obtained by Kir\'aly et al., see below.) However, we also show that the bound $r-k$ is not best possible in general.


title: Kloosterman sums and Maass cusp forms of half integral weight for the modular group

abstract: We estimate the sums \[ \sum_{c\leq x} \frac{S(m,n,c,\chi)}{c}, \] where the $S(m,n,c,\chi)$ are Kloosterman sums of half-integral weight on the modular group. Our estimates are uniform in $m$, $n$, and $x$ in analogy with Sarnak and Tsimerman's improvement of Kuznetsov's bound for the ordinary Kloosterman sums. Among other things this requires us to develop mean value estimates for coefficients of Maass cusp forms of weight $1/2$ and uniform estimates for $K$-Bessel integral transforms. As an application, we obtain an improved estimate for the classical problem of estimating the size of the error term in Rademacher's formula for the partition function $p(n)$.


title: Three Hours a Day: Understanding Current Teen Practices of Smartphone Application Use

abstract: Teens are using mobile devices for an increasing number of activities. Smartphones and a variety of mobile apps for communication, entertainment, and productivity have become an integral part of their lives. This mobile phone use has evolved rapidly as technology has changed and thus studies from even 2 or 3 years ago may not reflect new patterns and practices as smartphones have become more sophisticated. In order to understand current teen's practices around smartphone use, we conducted a two week, mixed-methods study with 14 diverse teens. Through voicemail diaries, interviews, and real world usage data from a logging application installed on their smartphones, we developed an understanding of the types of apps used by teens, when they use these apps, and their reasons for using specific apps in particular situations. We found that the teens in our study used their smartphones for an average of almost 3 hours per day and that two-thirds of all app use involved interacting with an average of almost 10 distinct communications applications. From our study data, we highlight key implications for the design of future mobile apps or services, specifically new social and communications-related applications that allow teens to maintain desired levels of privacy and permanence on the content that they share.


title: Source detection algorithms for dynamic contaminants based on the analysis of a hydrodynamic limit

abstract: In this work we propose and numerically analyze an algorithm for detection of a contaminant source using a dynamic sensor network. The algorithm is motivated using a global probabilistic optimization problem and is based on the analysis of the hydrodynamic limit of a discrete time evolution equation on the lattice under a suitable scaling of time and space. Numerical results illustrating the effectiveness of the algorithm are presented.


title: Bijections for planar maps with boundaries

abstract: We present bijections for planar maps with boundaries. In particular, we obtain bijections for triangulations and quadrangulations of the sphere with boundaries of prescribed lengths. For triangulations we recover the beautiful factorized formula obtained by Krikun using a (technically involved) generating function approach. The analogous formula for quadrangulations is new. We also obtain a far-reaching generalization for other face-degrees. In fact, all the known enumerative formulas for maps with boundaries are proved bijectively in the present article (and several new formulas are obtained). Our method is to show that maps with boundaries can be endowed with certain "canonical" orientations, making them amenable to the master bijection approach we developed in previous articles. As an application of our enumerative formulas, we note that they provide an exact solution of the dimer model on rooted triangulations and quadrangulations.


title: Homotopy groups of highly connected manifolds

abstract: In this paper we give a formula for the homotopy groups of $(n-1)$-connected $2n$-manifolds as a direct sum of homotopy groups of spheres in the case the $n^{th}$ Betti number is larger than $1$. We demonstrate that when the $n^{th}$ Betti number is $1$ the homotopy groups might not have such a decomposition. The techniques used in this computation also yield formulae for homotopy groups of connected sums of sphere products and CW complexes of a similar type. In all the families of spaces considered here, we establish a conjecture of J. C. Moore.


title: Percolation and coarse conformal uniformization

abstract: We formulate conjectures regarding percolation on planar triangulations suggested by assuming (quasi) invariance under coarse conformal uniformization.


title: Imaging Cyclotron Orbits of Electrons in Graphene

abstract: Electrons in graphene can travel for several microns without scattering at low temperatures, and their motion becomes ballistic, following classical trajectories. When a magnetic field B is applied perpendicular to the plane, electrons follow cyclotron orbits. Magnetic focusing occurs when electrons injected from one narrow contact focus onto a second contact located an integer number of cyclotron diameters away. By tuning the magnetic field B and electron density n in the graphene layer, we observe magnetic focusing peaks. We use a cooled scanning gate microscope to image cyclotron trajectories in graphene at 4.2 K. The tip creates a local change in density that casts a shadow by deflecting electrons flowing nearby; an image of flow can be obtained by measuring the transmission between contacts as the tip is raster scanned across the sample. On the first magnetic focusing peak, we image a cyclotron orbit that extends from one point contact to the other. In addition, we study the geometry of orbits deflected into the second point contact by the tip.


title: Learning multi-faceted representations of individuals from heterogeneous evidence using neural networks

abstract: Inferring latent attributes of people online is an important social computing task, but requires integrating the many heterogeneous sources of information available on the web. We propose learning individual representations of people using neural nets to integrate rich linguistic and network evidence gathered from social media. The algorithm is able to combine diverse cues, such as the text a person writes, their attributes (e.g. gender, employer, education, location) and social relations to other people. We show that by integrating both textual and network evidence, these representations offer improved performance at four important tasks in social media inference on Twitter: predicting (1) gender, (2) occupation, (3) location, and (4) friendships for users. Our approach scales to large datasets and the learned representations can be used as general features in and have the potential to benefit a large number of downstream tasks including link prediction, community detection, or probabilistic reasoning over social networks.


title: On the square function associated with generalized Bochner-Riesz means

abstract: We consider generalized Bochner-Riesz multipliers of the form $(1-\rho(\xi))_+^{\lambda}$ where $\rho:\mathbb{R}^2\to\mathbb{R}$ belongs to a class of rough distance functions homogeneous with respect to a nonisotropic dilation group. We prove a critical $L^4$ estimate for the associated square function, which we use to derive multiplier theorems for multipliers of the form $m\circ\rho$ where $m:\mathbb{R}\to\mathbb{C}$.


title: Magnetization dynamics in Fe$_x$Co$_{1-x}$ in presence of chemical disorder

abstract: In this paper, we present a theoretical formulation of magnetization dynamics in disordered binary alloys based on Kubo linear response theory interfaced with the combination of seamlessly three approaches; density functional based tight-binding linear muffin-tin orbitals, generalized recursion and Augmented space formalism. We apply this method to study the magnetization dynamics in chemically disordered Fe$_x$Co$_{1-x}$ ($x$ = 0.2, 0.5, 0.8) alloys. We reported that the magnon energies decrease with an increase in Co concentration. Significant magnon softening has been observed in Fe$_{20}$Co$_{80}$ at the Brillouin zone boundary. The magnon-electron scattering increases with increasing Co content which in turn modifies the hybridization between the Fe and Co atoms. This reduces the exchange energy between the atoms and soften down the magnon energy. The lowest magnon lifetime in found in Fe$_{50}$Co$_{50}$, where disorder is maximum. This clearly indicates that the damping of magnon energies in Fe$_x$Co$_{1-x}$ is governed by the hybridization between Fe and Co whereas the magnon lifetime is controlled by disorder configuration. Our atomistic spin dynamics simulations show a reasonable agreement with our theoretical approach in magnon dispersion for different alloy compositions.


title: On Termination of Polynomial Programs with Equality Conditions

abstract: We investigate the termination problem of a family of multi-path polynomial programs (MPPs), in which all assignments to program variables are polynomials, and test conditions of loops and conditional statements are polynomial equalities. We show that the set of non-terminating inputs (NTI) of such a program is algorithmically computable, thus leading to the decidability of its termination. To the best of our knowledge, the considered family of MPPs is hitherto the largest one for which termination is decidable. We present an explicit recursive function which is essentially Ackermannian, to compute the maximal length of ascending chains of polynomial ideals under a control function, and thereby obtain a complete answer to the questions raised by Seidenberg. This maximal length facilitates a precise complexity analysis of our algorithms for computing the NTI and deciding termination of MPPs. We extend our method to programs with polynomial guarded commands and show how an incomplete procedure for MPPs with inequality guards can be obtained. An application of our techniques to invariant generation of polynomial programs is further presented.