title: Joint Latency and Cost Optimization for Erasure-coded Data Center Storage

abstract: Modern distributed storage systems offer large capacity to satisfy the exponentially increasing need of storage space. They often use erasure codes to protect against disk and node failures to increase reliability, while trying to meet the latency requirements of the applications and clients. This paper provides an insightful upper bound on the average service delay of such erasure-coded storage with arbitrary service time distribution and consisting of multiple heterogeneous files. Not only does the result supersede known delay bounds that only work for a single file or homogeneous files, it also enables a novel problem of joint latency and storage cost minimization over three dimensions: selecting the erasure code, placement of encoded chunks, and optimizing scheduling policy. The problem is efficiently solved via the computation of a sequence of convex approximations with provable convergence. We further prototype our solution in an open-source, cloud storage deployment over three geographically distributed data centers. Experimental results validate our theoretical delay analysis and show significant latency reduction, providing valuable insights into the proposed latency-cost tradeoff in erasure-coded storage.


title: Observational implications of mattergenesis during inflation

abstract: The observed baryon asymmetry, as well as potentially an asymmetry in the dark matter sector, can be produced through dissipative particle production during inflation. A distinctive feature of this mechanism is the generation of matter isocurvature perturbations that are fully (anti-)correlated with the dominant adiabatic curvature perturbations. We show that chaotic warm inflation models yield anti-correlated isocurvature modes that may partially or even completely screen the contribution of primordial gravity waves to the CMB temperature power spectrum. The tensor-to-scalar ratio inferred from the latter may thus be parametrically smaller than the one deduced from B-mode polarization maps, which is particularly relevant in the light of the recently announced results of the BICEP2 experiment.


title: A weighty interpretation of the Galactic Centre excess

abstract: Previous attempts at explaining the gamma-ray excess near the Galactic Centre have focussed on dark matter annihilation directly into Standard Model particles. This results in a preferred dark matter mass of 30-40 GeV (if the annihilation is into b quarks) or 10 GeV (if it is into leptons). Here we show that the gamma-ray excess is also consistent with heavier dark matter particles; in models of secluded dark matter, dark matter with mass up to 76 GeV provides a good fit to the data. This occurs if the dark matter first annihilates to an on-shell particle that subsequently decays to Standard Model particles through a portal interaction. This is a generic process that works in models with annihilation, semi-annihilation or both. We explicitly demonstrate this in a model of hidden vector dark matter with an SU(2) gauge group in the hidden sector.


title: Logic and operator algebras

abstract: The most recent wave of applications of logic to operator algebras is a young and rapidly developing field. This is a snapshot of the current state of the art.


title: Wireless Josephson Amplifier

abstract: Josephson junction parametric amplifiers are playing a crucial role in the readout chain in superconducting quantum information experiments. However, their integration with current 3D cavity implementations poses the problem of transitioning between waveguide, coax cables and planar circuits. Moreover, Josephson amplifiers require auxiliary microwave components, like directional couplers and/or hybrids, that are sources of spurious losses and impedance mismatches that limit measurement efficiency and amplifier tunability. We have developed a new wireless architecture for these parametric amplifiers that eliminates superfluous microwave components and interconnects. This greatly simplifies their assembly and integration into experiments. We present an experimental realization of such a device operating in the $9-11$ GHz band with about $100$ MHz of amplitude gain-bandwidth product, on par with devices mounted in conventional sample holders. The simpler impedance environment presented to the amplifier also results in increased amplifier tunability.


title: Monge-Kantorovich norms on spaces of vector measures

abstract: One considers Hilbert space valued measures on the Borel sets of a compact metric space. A natural numerical valued integral of vector valued continuous functions with respect to vector valued functions is defined. Using this integral, different norms (we called them Monge-Kantorovich norm, modified Monge-Kantorovich norm and Hanin norm) on the space of measures are introduced, generalizing the theory of (weak) convergence for probability measures on metric spaces. These norms introduce new (equivalent) metrics on the initial compact metric space.


title: The Panchromatic Hubble Andromeda Treasury. VI. The reliability of far-ultraviolet flux as a star formation tracer on sub-kpc scales

abstract: We have used optical observations of resolved stars from the Panchromatic Hubble Andromeda Treasury (PHAT) to measure the recent (< 500 Myr) star formation histories (SFHs) of 33 FUV-bright regions in M31. The region areas ranged from ~$10^4$ to $10^6$ pc$^2$, which allowed us to test the reliability of FUV flux as a tracer of recent star formation on sub-kpc scales. The star formation rates (SFRs) derived from the extinction-corrected observed FUV fluxes were, on average, consistent with the 100-Myr mean SFRs of the SFHs to within the 1$\sigma$ scatter. Overall, the scatter was larger than the uncertainties in the SFRs and particularly evident among the smallest regions. The scatter was consistent with an even combination of discrete sampling of the initial mass function and high variability in the SFHs. This result demonstrates the importance of satisfying both the full-IMF and the constant-SFR assumptions for obtaining precise SFR estimates from FUV flux. Assuming a robust FUV extinction correction, we estimate that a factor of 2.5 uncertainty can be expected in FUV-based SFRs for regions smaller than $10^5$ pc$^2$, or a few hundred pc. We also examined ages and masses derived from UV flux under the common assumption that the regions are simple stellar populations (SSPs). The SFHs showed that most of the regions are not SSPs, and the age and mass estimates were correspondingly discrepant from the SFHs. For those regions with SSP-like SFHs, we found mean discrepancies of 10 Myr in age and a factor of 3 to 4 in mass. It was not possible to distinguish the SSP-like regions from the others based on integrated FUV flux.


title: Dynamic and Multi-functional Labeling Schemes

abstract: We investigate labeling schemes supporting adjacency, ancestry, sibling, and connectivity queries in forests. In the course of more than 20 years, the existence of $\log n + O(\log \log)$ labeling schemes supporting each of these functions was proven, with the most recent being ancestry [Fraigniaud and Korman, STOC '10]. Several multi-functional labeling schemes also enjoy lower or upper bounds of $\log n + \Omega(\log \log n)$ or $\log n + O(\log \log n)$ respectively. Notably an upper bound of $\log n + 5\log \log n$ for adjacency+siblings and a lower bound of $\log n + \log \log n$ for each of the functions siblings, ancestry, and connectivity [Alstrup et al., SODA '03]. We improve the constants hidden in the $O$-notation. In particular we show a $\log n + 2\log \log n$ lower bound for connectivity+ancestry and connectivity+siblings, as well as an upper bound of $\log n + 3\log \log n + O(\log \log \log n)$ for connectivity+adjacency+siblings by altering existing methods. In the context of dynamic labeling schemes it is known that ancestry requires $\Omega(n)$ bits [Cohen, et al. PODS '02]. In contrast, we show upper and lower bounds on the label size for adjacency, siblings, and connectivity of $2\log n$ bits, and $3 \log n$ to support all three functions. There exist efficient adjacency labeling schemes for planar, bounded treewidth, bounded arboricity and interval graphs. In a dynamic setting, we show a lower bound of $\Omega(n)$ for each of those families.


title: Shiva++: An Enhanced Graph based Ontology Matcher

abstract: With the web getting bigger and assimilating knowledge about different concepts and domains, it is becoming very difficult for simple database driven applications to capture the data for a domain. Thus developers have come out with ontology based systems which can store large amount of information and can apply reasoning and produce timely information. Thus facilitating effective knowledge management. Though this approach has made our lives easier, but at the same time has given rise to another problem. Two different ontologies assimilating same knowledge tend to use different terms for the same concepts. This creates confusion among knowledge engineers and workers, as they do not know which is a better term then the other. Thus we need to merge ontologies working on same domain so that the engineers can develop a better application over it. This paper shows the development of one such matcher which merges the concepts available in two ontologies at two levels; 1) at string level and 2) at semantic level; thus producing better merged ontologies. We have used a graph matching technique which works at the core of the system. We have also evaluated the system and have tested its performance with its predecessor which works only on string matching. Thus current approach produces better results.


title: Information Theoretic Analysis of Concurrent Information Transfer and Power Gain

abstract: In this paper, we analyze the fundamental trade-off between information transfer and power gain by means of an information-theoretic framework in communications circuits. This analysis is of interest as many of today's applications require that maximum information and maximum signal power are extracted (or transferred) through the circuit at the same time for further processing so that a compromise concerning the signal spectral shape as well as the matching network has to be found. To this end, the optimization framework is applied to a two-port circuit, which is used as an abstraction for a broadband amplifier. Thereby, we characterize the involved Pareto bound by considering different optimization problems. The first one aims at optimizing the input power spectral density (PSD) as well as the source and load admittances, whereas the second approach assumes the PSD to be fixed and uniformly distributed within a fixed bandwidth and optimizes the source and load admittances only. Moreover, we will show that additional matching networks may help to improve the trade-off.


title: Learning in Repeated Games: Human Versus Machine

abstract: While Artificial Intelligence has successfully outperformed humans in complex combinatorial games (such as chess and checkers), humans have retained their supremacy in social interactions that require intuition and adaptation, such as cooperation and coordination games. Despite significant advances in learning algorithms, most algorithms adapt at times scales which are not relevant for interactions with humans, and therefore the advances in AI on this front have remained of a more theoretical nature. This has also hindered the experimental evaluation of how these algorithms perform against humans, as the length of experiments needed to evaluate them is beyond what humans are reasonably expected to endure (max 100 repetitions). This scenario is rapidly changing, as recent algorithms are able to converge to their functional regimes in shorter time-scales. Additionally, this shift opens up possibilities for experimental investigation: where do humans stand compared with these new algorithms? We evaluate humans experimentally against a representative element of these fast-converging algorithms. Our results indicate that the performance of at least one of these algorithms is comparable to, and even exceeds, the performance of people.


title: Multiparticle SYM equations of motion and pure spinor BRST blocks

abstract: In this paper a multiparticle generalization of linearized ten-dimensional super Yang--Mills superfields is proposed. Their equations of motions are shown to take the same form as in the single-particle case, supplemented by contact terms. A recursive construction of these superfields is inspired by the iterated OPEs among massless vertex operators in the pure spinor formalism. An enlarged set of BRST-covariant pure spinor blocks is then defined in a streamlined fashion and combined to multiparticle vertex operators. The latter can be used to universally describe tree-level subdiagrams in the perturbative open and closed superstring, regardless of the loop order. The inherent symmetries of the multiparticle superfields are reproduced by structure constants of the gauge group, hinting a natural appearance of the BCJ-duality between color and kinematics in the fundamentals of super Yang--Mills theory. We present one-loop applications where known scalar cohomology objects are systematically recomputed and a novel vector cohomology particularly relevant to the closed string is constructed for arbitrary multiplicity.


title: Between 2- and 3-colorability

abstract: We consider the question of the existence of homomorphisms between $G_{n,p}$ and odd cycles when $p=c/n,\,1<c\leq 4$. We show that for any positive integer $\ell$, there exists $\epsilon=\epsilon(\ell)$ such that if $c=1+\epsilon$ then w.h.p. $G_{n,p}$ has a homomorphism from $G_{n,p}$ to $C_{2\ell+1}$ so long as its odd-girth is at least $2\ell+1$. On the other hand, we show that if $c=4$ then w.h.p. there is no homomorphism from $G_{n,p}$ to $C_5$. Note that in our range of interest, $\chi(G_{n,p})=3$ w.h.p., implying that there is a homomorphism from $G_{n,p}$ to $C_3$.


title: Neighborhoods on the Grasmannian of marginals with bounded isotropic constant

abstract: We show that for any isotropic log-concave probability measure $\mu$ on $\mathbb R^n$, for every $\varepsilon > 0$, every $1 \leq k \leq \sqrt{n}$ and any $E \in G_{n,k}$ there exists $F \in G_{n,k}$ with $d(E,F) < \varepsilon$ and $L_{\pi_F\mu} < C/\varepsilon$.


title: Dependent Lindeberg CLT - Finite Dimensional for Empirical Processes of Cluster Functionals

abstract: Drees and Rootz\'en [2010] have proven central limit theorems (CLT) for empirical processes of extreme values cluster functionals built from $\beta$-mixing processes. The problem with this family of $\beta$-mixing processes is that it is quite restrictive, as has been shown by Andrews [1984]. We expand this result to a more general dependent processes family, known as weakly dependent processes in the sense of Doukhan and Louhichi [1999], but in finite-dimensional convergence (fidis). We show an example where the application of the CLT-fidis is sufficient in several cases, including a small simulation of the extremogram introduced by Davis and Mikosch [2009] to confirm the efficacy of our result.


title: N-body Simulations of Terrestrial Planet Formation under the Influence of a Hot Jupiter

abstract: We investigate the formation of multiple-planet systems in the presence of a hot Jupiter using extended N-body simulations that are performed simultaneously with semi-analytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of a hot Jupiter. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal-mass; the largest planet constitutes more than 50 percent of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper (Ogihara et al. 2013), we have found a new physical mechanism of induced migration of the hot Jupiter, which is called a crowding-out. If the hot Jupiter opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the hot Jupiter remains. We also discuss angular momentum transfer between the planets and disk.


title: Spectral gap estimates for some block matrices

abstract: We estimate the size of the spectral gap at zero for some Hermitian block matrices. Included are quasi-definite matrices, quasi-semidefinite matrices (the closure of the set of the quasi-definite matrices) and some related block matrices which need not belong to either of these classes. Matrices of such structure arise in quantum models of possibly disordered systems with supersymmetry or graphene like symmetry. Some of the results immediately extend to infinite dimension.


title: Strong Dependence of the Inner Edge of the Habitable Zone on Planetary Rotation Rate

abstract: Planetary rotation rate is a key parameter in determining atmospheric circulation and hence the spatial pattern of clouds. Since clouds can exert a dominant control on planetary radiation balance, rotation rate could be critical for determining mean planetary climate. Here we investigate this idea using a three-dimensional general circulation model with a sophisticated cloud scheme. We find that slowly rotating planets (like Venus) can maintain an Earth-like climate at nearly twice the stellar flux as rapidly rotating planets (like Earth). This suggests that many exoplanets previously believed to be too hot may actually be habitable, depending on their rotation rate. The explanation for this behavior is that slowly rotating planets have a weak Coriolis force and long daytime illumination, which promotes strong convergence and convection in the substellar region. This produces a large area of optically thick clouds, which greatly increases the planetary albedo. In contrast, on rapidly rotating planets a much narrower belt of clouds form in the deep tropics, leading to a relatively low albedo. A particularly striking example of the importance of rotation rate suggested by our simulations is that a planet with modern Earth's atmosphere, in Venus' orbit, and with modern Venus' (slow) rotation rate would be habitable. This would imply that if Venus went through a runaway greenhouse, it had a higher rotation rate at that time.


title: Plasmons in a Superlattice of Fullerenes or Metallic Shells

abstract: A theory for the collective plasma excitations in a linear periodic array of spherical two-dimensional electron gases (S2DEGs) is presented. This is a simple model for an ultra thin and narrow microribbon of fullerenes or metallic shells. Coulomb coupling between electrons located on the same sphere and on different spheres is included in the random-phase approximation (RPA). Electron hopping between spheres is neglected in these calculations. The resulting plasmon-dispersion equation is solved numerically. Results are presented for a superlattice of single-wall S2DEGs as a function of the wave vector. The plasmon dispersions are obtained for different spherical separations. We show that the one-dimensional translational symmetry of the lattice is maintained in the plasmon spectrum. Additionally, we compare the plasmon dispersion when the superlatice direction is parallel or perpendicular to the axis of quantization. However, because of anisotropy in the Coulomb matrix elements, there is anticrossing in the plasmon dispersion only in the direction perpendicular to the quantization axis. The S2DEG may serve as a simple model for fullerenes, when their energy bands are far apart.


title: Kinetic theory and numerical simulations of two-species coagulation

abstract: In this work we study the stochastic process of two-species coagulation. This process consists in the aggregation dynamics taking place in a ring. Particles and clusters of particles are set in this ring and they can move either clockwise or counterclockwise. They have a probability to aggregate forming larger clusters when they collide with another particle or cluster. We study the stochastic process both analytically and numerically. Analytically, we derive a kinetic theory which approximately describes the process dynamics. One of our strongest assumptions in this respect is the so called well-stirred limit, that allows neglecting the appearance of spatial coordinates in the theory, so this becomes effectively reduced to a zeroth dimensional model. We determine the long time behavior of such a model, making emphasis in one special case in which it displays self-similar solutions. In particular these calculations answer the question of how the system gets ordered, with all particles and clusters moving in the same direction, in the long time. We compare our analytical results with direct numerical simulations of the stochastic process and both corroborate its predictions and check its limitations. In particular, we numerically confirm the ordering dynamics predicted by the kinetic theory and explore properties of the realizations of the stochastic process which are not accessible to our theoretical approach.

categories: math.AP cond-mat.stat-mech math-ph math.MP nlin.AO


title: A Generalized Cut-Set Bound for Deterministic Multi-Flow Networks and its Applications

abstract: We present a new outer bound for the sum capacity of general multi-unicast deterministic networks. Intuitively, this bound can be understood as applying the cut-set bound to concatenated copies of the original network with a special restriction on the allowed transmit signal distributions. We first study applications to finite-field networks, where we obtain a general outer-bound expression in terms of ranks of the transfer matrices. We then show that, even though our outer bound is for deterministic networks, a recent result relating the capacity of AWGN KxKxK networks and the capacity of a deterministic counterpart allows us to establish an outer bound to the DoF of KxKxK wireless networks with general connectivity. This bound is tight in the case of the "adjacent-cell interference" topology, and yields graph-theoretic necessary and sufficient conditions for K DoF to be achievable in general topologies.


title: The solar abundance problem: the effect of the turbulent kinetic flux on the solar envelope model

abstract: Recent 3D-simulations have shown that the turbulent kinetic flux (TKF) is significant. We discuss the effects of TKF on the size of convection zone and find that the TKF may help to solve the solar abundance problem. The solar abundance problem is that, with new abundances, the solar convection zone depth, sound speed in the radiative interior, the helium abundance and density in the convective envelope are not in agreement with helioseismic inversions. We have done Monte Carlo simulations on solar convective envelope models with different profile of TKF to test the effects. The solar abundance problem is revealed in the standard solar convective envelope model with AGSS09 composition, which shows significant differences ($\rm{\sim 10 %}$) on density from the helioseicmic inversions, but the differences in the model with old composition GN93 is small ($\rm{\sim 0.5 %}$). In the testing models with different imposed TKF, it is found that the density profile is sensitive to the value of TKF at the base of convective envelope and insensitive to the structure of TKF in the convection zone. Required value of turbulent kinetic luminosity at the base is about $\rm{-13%\sim-19%L_{\odot}}$. Comparing with the 3D-simulations, this value is plausible. This study is for the solar convective envelope only. The evolutionary solar models with TKF are required for investigating its effects on the solar interior structure below the convection zone and the whole solar abundance problem, but the profile of TKF in the overshoot region is needed.


title: Tight bounds for learning a mixture of two gaussians

abstract: We consider the problem of identifying the parameters of an unknown mixture of two arbitrary $d$-dimensional gaussians from a sequence of independent random samples. Our main results are upper and lower bounds giving a computationally efficient moment-based estimator with an optimal convergence rate, thus resolving a problem introduced by Pearson (1894). Denoting by $\sigma^2$ the variance of the unknown mixture, we prove that $\Theta(\sigma^{12})$ samples are necessary and sufficient to estimate each parameter up to constant additive error when $d=1.$ Our upper bound extends to arbitrary dimension $d>1$ up to a (provably necessary) logarithmic loss in $d$ using a novel---yet simple---dimensionality reduction technique. We further identify several interesting special cases where the sample complexity is notably smaller than our optimal worst-case bound. For instance, if the means of the two components are separated by $\Omega(\sigma)$ the sample complexity reduces to $O(\sigma^2)$ and this is again optimal. Our results also apply to learning each component of the mixture up to small error in total variation distance, where our algorithm gives strong improvements in sample complexity over previous work. We also extend our lower bound to mixtures of $k$ Gaussians, showing that $\Omega(\sigma^{6k-2})$ samples are necessary to estimate each parameter up to constant additive error.


title: Nonharmonic Phonons in $\alpha$-Iron at High Temperatures

abstract: Phonon densities of states (DOS) of bcc $\alpha$-$^{57}$Fe were measured from room temperature through the 1044K Curie transition and the 1185K fcc $\gamma$-Fe phase transition using nuclear resonant inelastic x-ray scattering. At higher temperatures all phonons shift to lower energies (soften) with thermal expansion, but the low transverse modes soften especially rapidly above 700K, showing strongly nonharmonic behavior that persists through the magnetic transition. Interatomic force constants for the bcc phase were obtained by iteratively fitting a Born-von K\'{a}rm\'{a}n model to the experimental phonon spectra using a genetic algorithm optimization. The second-nearest-neighbor fitted axial force constants weakened significantly at elevated temperatures. An unusually large nonharmonic behavior is reported, which increases the vibrational entropy and accounts for a contribution of 35 meV/atom in the free energy at high temperatures. The nonharmonic contribution to the vibrational entropy follows the thermal trend of the magnetic entropy, and may be coupled to magnetic excitations. A small change in vibrational entropy across the $\alpha$-$\gamma$ structural phase transformation is also reported.


title: Real multiplication revisited

abstract: It is proved that the Hilbert class field of a real quadratic field ${\Bbb Q}(\sqrt{D})$ modulo a power $m$ of the conductor $f$ is generated by the Fourier coefficients of the Hecke eigenform for a congruence subgroup of level $fD$.


title: Clinching Auctions Beyond Hard Budget Constraints

abstract: Constraints on agent's ability to pay play a major role in auction design for any setting where the magnitude of financial transactions is sufficiently large. Those constraints have been traditionally modeled in mechanism design as \emph{hard budget}, i.e., mechanism is not allowed to charge agents more than a certain amount. Yet, real auction systems (such as Google AdWords) allow more sophisticated constraints on agents' ability to pay, such as \emph{average budgets}. In this work, we investigate the design of Pareto optimal and incentive compatible auctions for agents with \emph{constrained quasi-linear utilities}, which captures more realistic models of liquidity constraints that the agents may have. Our result applies to a very general class of allocation constraints known as polymatroidal environments, encompassing many settings of interest such as multi-unit auctions, matching markets, video-on-demand and advertisement systems. Our design is based Ausubel's \emph{clinching framework}. Incentive compatibility and feasibility with respect to ability-to-pay constraints are direct consequences of the clinching framework. Pareto-optimality, on the other hand, is considerably more challenging, since the no-trade condition that characterizes it depends not only on whether agents have their budgets exhausted or not, but also on prices {at} which the goods are allocated. In order to get a handle on those prices, we introduce novel concepts of dropping prices and saturation. These concepts lead to our main structural result which is a characterization of the tight sets in the clinching auction outcome and its relation to dropping prices.


title: The variety of three-dimensional real Jordan algebras

abstract: In this paper, we study the variety $Jor_{3}$ of three-dimensional Jordan algebras over the field of real numbers. We establish the list of $26$ non-isomorphic Jordan algebras and describe the irreducible components of $Jor_{3}$ proving that it is the union of Zariski closure of the orbits of $8$ rigid algebras.


title: A Geometric Distance Oracle for Large Real-World Graphs

abstract: Many graph processing algorithms require determination of shortest-path distances between arbitrary numbers of node pairs. Since computation of exact distances between all node-pairs of a large graph, e.g., 10M nodes and up, is prohibitively expensive both in computational time and storage space, distance approximation is often used in place of exact computation. In this paper, we present a novel and scalable distance oracle that leverages the hyperbolic core of real-world large graphs for fast and scalable distance approximation. We show empirically that the proposed oracle significantly outperforms prior oracles on a random set of test cases drawn from public domain graph libraries. There are two sets of prior work against which we benchmark our approach. The first set, which often outperforms other oracles, employs embedding of the graph into low dimensional Euclidean spaces with carefully constructed hyperbolic distances, but provides no guarantees on the distance estimation error. The second set leverages Gromov-type tree contraction of the graph with the additive error guaranteed not to exceed $2\delta\log{n}$, where $\delta$ is the hyperbolic constant of the graph. We show that our proposed oracle 1) is significantly faster than those oracles that use hyperbolic embedding (first set) with similar approximation error and, perhaps surprisingly, 2) exhibits substantially lower average estimation error compared to Gromov-like tree contractions (second set). We substantiate our claims through numerical computations on a collection of a dozen real world networks and synthetic test cases from multiple domains, ranging in size from 10s of thousand to 10s of millions of nodes.


title: A generalization of Kuo condensation

abstract: Kuo introduced his 4-point condensation in 2003 for bipartite planar graphs. In 2006 Kuo generalized this 4-point condensation to planar graphs that are not necessarily bipartite. His formula expressed the product between the number of perfect matching of the original graph $G$ and that of the subgraph obtained from $G$ by removing the four distinguished vertices as a Pfaffian of order 4, whose entries are numbers of perfect matchings of subgraphs of $G$ obtained by removing various pairs of vertices chosen from among the four distinguished ones. The compelling elegance of this formula is inviting of generalization. Kuo generalized it to $2k$ points under the special assumption that the subgraph obtained by removing some subset of the $2k$ vertices has precisely one perfect matching. In this paper we prove that the formula holds in the general case. We also present a couple of applications.


title: Temporal Evolution of the Magnetic Topology of the NOAA Active Region 11158

abstract: We studied the temporal evolution of the magnetic topology of the active region (AR) 11158 based on the reconstructed three-dimensional magnetic fields in the corona. The \nlfff\ extrapolation method was applied to the 12 minutes cadence data obtained with the \hmi\ (HMI) onboard the \sdo\ (SDO) during five days. By calculating the squashing degree factor Q in the volume, the derived quasi-separatrix layers (QSLs) show that this AR has an overall topology, resulting from a magnetic quadrupole, including an hyperbolic flux tube (HFT) configuration which is relatively stable at the time scale of the flare ($\sim 1-2$ hours). A strong QSL, which corresponds to some highly sheared arcades that might be related to the formation of a flux rope, is prominent just before the M6.6 and X2.2 flares, respectively. These facts indicate the close relationship between the strong QSL and the high flare productivity of AR 11158. In addition, with a close inspection of the topology, we found a small-scale HFT which has an inverse tear-drop structure above the aforementioned QSL before the X2.2 flare. It indicates the existence of magnetic flux rope at this place. Even though a global configuration (HFT) is recognized in this AR, it turns out that the large-scale HFT only plays a secondary role during the eruption. In final, we dismiss a trigger based on the breakout model and highlight the central role of the flux rope in the related eruption.


title: Valley Order and Loop Currents in Graphene on Hexagonal Boron Nitride

abstract: In this letter, we examine the role of Coulomb interactions in the emergence of macroscopically ordered states in graphene supported on hexagonal boron nitride substrates. Due to incommensuration effects with the substrate, graphene can develop gapped low energy modes that spatially conform into a triangular superlattice of quantum rings. In the presence of these modes, we show that Coulomb interactions lead to spontaneous formation of chiral loop currents in bulk and to macroscopic spin-valley order at zero temperature. We show that this exotic state breaks time reversal symmetry and can be detected with interferometry and polar Kerr measurements.


title: On an Enneper-Weierstrass-type representation of constant Gaussian curvature surfaces in $3$-dimensional hyperbolic space

abstract: For all $k\in]0,1[$, we construct a canonical bijection between the space of ramified coverings of the sphere and the space of complete immersed surfaces in $3$-dimensional hyperbolic space of finite area and of constant extrinsic curvature equal to $k$. We show, furthermore, that this bijection restricts to a homeomorphism over each stratum of the space of ramified coverings of the sphere.


title: Secure Degrees of Freedom of the MIMO Multiple Access Channel with Multiple unknown Eavesdroppers

abstract: We investigate the secure degrees of freedom (SDoF) of a two-transmitter Gaussian multiple access channel with multiple antennas at the transmitters, the legitimate receiver with the existence of an unknown number of eavesdroppers each with a number of antennas less than or equal to a known value $N_E$. The channel matrices between the legitimate transmitters and the receiver are available everywhere, while the legitimate pair does not know the eavesdroppers' channels matrices. We provide the exact sum SDoF for the considered system. A new comprehensive upperbound is deduced and a new achievable scheme based on utilizing jamming is exploited. We prove that Cooperative Jamming is SDoF optimal even without the instantaneous eavesdropper CSI available at the transmitters.


title: Dynamical behavior connection of the gluon distribution and the proton structure function at small $x$

abstract: We make a critical study of the relationship between the singlet structure function $F_{2}^{S}$ and the gluon distribution $G(x,Q^{2})$ proposed in the past two decades, which is frequently used to extract the gluon distribution from the proton structure function. We show that a simple relation is not generally valid in the simplest state. We complete this relation by using a Laplace transform method and hard-pomeron behavior at LO and NLO at small $x$. Our study shows that this relation is dependent on the splitting functions and initial conditions at $Q^{2} = Q^{2}_{0}$ and running coupling constant at NLO. The resulting analytic expression allows us to predict the proton structure function with respect to the gluon distributions and to compare the results with H1 data and a QCD analysis fit. Comparisons with other results are made and predictions for the proposed best approach are also provided.


title: Efficient Semidefinite Branch-and-Cut for MAP-MRF Inference

abstract: We propose a Branch-and-Cut (B&C) method for solving general MAP-MRF inference problems. The core of our method is a very efficient bounding procedure, which combines scalable semidefinite programming (SDP) and a cutting-plane method for seeking violated constraints. In order to further speed up the computation, several strategies have been exploited, including model reduction, warm start and removal of inactive constraints. We analyze the performance of the proposed method under different settings, and demonstrate that our method either outperforms or performs on par with state-of-the-art approaches. Especially when the connectivities are dense or when the relative magnitudes of the unary costs are low, we achieve the best reported results. Experiments show that the proposed algorithm achieves better approximation than the state-of-the-art methods within a variety of time budgets on challenging non-submodular MAP-MRF inference problems.


title: Mapping to a Reference Genome Structure

abstract: To support comparative genomics, population genetics, and medical genetics, we propose that a reference genome should come with a scheme for mapping each base in any DNA string to a position in that reference genome. We refer to a collection of one or more reference genomes and a scheme for mapping to their positions as a reference structure. Here we describe the desirable properties of reference structures and give examples. To account for natural genetic variation, we consider the more general case in which a reference genome is represented by a graph rather than a set of phased chromosomes; the latter is treated as a special case.


title: Categorical Bockstein sequences

abstract: We construct the reduction of an exact category with a twist functor with respect to an element of its graded center in presence of an exact-conservative forgetful functor annihilating this central element. The procedure allows, e.g., to recover the abelian/exact category of modular representations of a finite group from the exact category of its l-adic representations. The construction uses matrix factorizations in a nontraditional way. We obtain the Bockstein long exact sequences for the Ext groups in the exact categories produced by reduction. Our motivation comes from the theory of Artin-Tate motives and motivic sheaves with finite coefficients, and our key techniques generalize those of Section 4 in arXiv:1006.4343


title: On the MacWilliams Identity for Classical and Quantum Convolutional Codes

abstract: The weight generating functions associated with convolutional codes (CCs) are based on state space realizations or the weight adjacency matrices (WAMs). The MacWilliams identity for CCs on the WAMs was first established by Gluesing- Luerssen and Schneider in the case of minimal encoders, and generalized by Forney. We consider this problem in the viewpoint of constraint codes and obtain a simple and direct proof of this MacWilliams identity in the case of minimal encoders. For our purpose, we choose a different representation for the exact weight generating function (EWGF) of a block code, by defining it as a linear combination of orthonormal vectors in Dirac bra-ket notation. This representation provides great flexibility so that general split weight generating functions and their MacWilliams identities can be easily obtained from the MacWilliams identity for EWGFs. As a result, we also obtain the MacWilliams identity for the input-parity weight adjacency matrices of a systematic convolutional code and its dual. Finally, paralleling the development of the classical case, we establish the MacWilliams identity for quantum convolutional codes.


title: Notes on a conjecture of Manoussakis concerning Hamilton cycles in digraphs

abstract: In 1992, Manoussakis conjectured that a strongly 2-connected digraph $D$ on $n$ vertices is hamiltonian if for every two distinct pairs of independent vertices $x,y$ and $w,z$ we have $d(x)+d(y)+d(w)+d(z)\geq 4n-3$. In this note we show that $D$ has a Hamilton path, which gives an affirmative evidence supporting this conjecture.


title: Resonant bands, Aomoto complex, and real 4-nets

abstract: The resonant band is a useful notion for the computation of the nontrivial monodromy eigenspaces of the Milnor fiber of a real line arrangement. In this article, we develop the resonant band description for the cohomology of the Aomoto complex. As an application, we prove that real 4-nets do not exist.


title: Linear Turan numbers of r-uniform linear cycles and related Ramsey numbers

abstract: An $r$-uniform hypergraph is called an $r$-graph. A hypergraph is linear if every two edges intersect in at most one vertex. Given a linear $r$-graph $H$ and a positive integer $n$, the linear Tur\'an number $ex_L(n,H)$ is the maximum number of edges in a linear $r$-graph $G$ that does not contain $H$ as a subgraph. For each $\ell\geq 3$, let $C^r_\ell$ denote the $r$-uniform linear cycle of length $\ell$, which is an $r$-graph with edges $e_1,\ldots, e_\ell$ such that $\forall i\in [\ell-1]$, $|e_i\cap e_{i+1}|=1$, $|e_\ell\cap e_1|=1$ and $e_i\cap e_j=\emptyset$ for all other pairs $\{i,j\}, i\neq j$. For all $r\geq 3$ and $\ell\geq 3$, we show that there exist positive constants $c_{m,r}$ and $c'_{m,r}$, depending only $m$ and $r$, such that $ex_L(n,C^r_{2m})\leq c_{m,r} n^{1+\frac{1}{m}}$ and $ex_L(n,C^r_{2m+1})\leq c'_{m,r} n^{1+\frac{1}{m}}$. This answers a question of Kostochka, Mubayi, and Verstra\"ete. For even cycles, our result extends the result of Bondy and Simonovits on the Tur\'an numbers of even cycles to linear hypergraphs. Using our results on linear Tur\'an numbers we also obtain bounds on the cycle-complete hypergraph Ramsey numbers. We show that there are positive constants $a_{m,r}$ and $b_{m,r}$, depending only on $m$ and $r$, such that $R(C^r_{2m}, K^r_t)\leq a_{m,r} (\frac{t}{\ln t})^\frac{m}{m-1}$ and $R(C^r_{2m+1}, K^r_t)\leq b_{m,r} t^\frac{m}{m-1}$.


title: Spherical harmonics with maximal Lp (2<p<=6) norm growth

abstract: In this paper, we show that there exists a positive density subsequence of orthonormal spherical harmonics which achieves the maximal Lp norm growth for 2<p<=6, therefore giving an example of a Riemannian surface supporting such subsequence of eigenfunctions. This answers the question proposed by Sogge and Zelditch (arXiv:1011.0215). Furthermore, we provide an explicit lower bound on the density in this example.


title: A decade of dark matter searches with ground-based Cherenkov telescopes

abstract: In the general scenario of Weakly Interacting Massive Particles (WIMP), dark matter (DM) can be observed via astrophysical gamma rays because photons are produced in various DM annihilation or decay processes, either as broad-band or line emission, or because of the secondary processes of charged particles in the final stages of the annihilations or the decays. The energy range of the former processes is accessible by current ground-based Imaging Atmospheric Cherenkov telescopes (IACTs, like H.E.S.S., MAGIC and VERITAS). The strengths of this technique are: a) the expected DM gamma-ray spectra show peculiar features like bumps, spikes and cutoff that make them clearly distinguishable from the smoother astrophysical spectra, b) the expected DM spectrum is universal and therefore by observing two or more DM targets with the same spectrum, a clear identification (besides detection) of DM would be enabled. The role of IACTs may gain more importance in the future as the results from the LHC may hint to a DM particle with mass at the TeV or above, where the IACTs sensitivity is unsurpassed by other experiments. In this contribution, a review of the search for DM with the current generation of IACT will be presented.


title: The Formation of Uranus & Neptune: Challenges and Implications For Intermediate-Mass Exoplanets

abstract: In this paper we investigate the formation of Uranus and Neptune, according to the core-nucleated accretion model, considering formation locations ranging from 12 to 30 AU from the Sun, and with various disk solid-surface densities and core accretion rates. It is shown that in order to form Uranus-like and Neptune-like planets in terms of final mass and solid-to-gas ratio, very specific conditions are required. We also show that when recently proposed high solid accretion rates are assumed, along with solid surface densities about 10 times those in the minimum-mass solar nebula, the challenge in forming Uranus and Neptune at large radial distances is no longer the formation timescale, but is rather finding agreement with the final mass and composition of these planets. In fact, these conditions are more likely to lead to gas-giant planets. Scattering of planetesimals by the forming planetary core is found to be an important effect at the larger distances. Our study emphasizes how (even slightly) different conditions in the protoplanetary disk and the birth environment of the planetary embryos can lead to the formation of very different planets in terms of final masses and compositions (solid-to-gas ratios), which naturally explains the large diversity of intermediate-mass exoplanets.


title: Compressive Joint Angular-Frequency Power Spectrum Estimation

abstract: We introduce a new compressive power spectrum estimation approach in both frequency and direction of arrival (DOA). Wide-sense stationary signals produced by multiple uncorrelated sources are compressed in both the time and spatial domain where the latter compression is implemented by activating only some of the antennas in the underlying uniform linear array (ULA).We sample the received signal at every active antenna at sub-Nyquist rate, compute both the temporal and spatial correlation functions between the sub-Nyquist rate samples, and apply least squares to reconstruct the full-blown two-dimensional power spectrum matrix where the rows and columns correspond to the frequencies and the angles, respectively. This is possible under the full column rank condition of the system matrices and without applying any sparsity constraint on the signal statistics. Further, we can estimate the DOAs of the sources by locating the peaks of the angular power spectrum. We can theoretically estimate the frequency bands and the DOAs of more uncorrelated sources than active sensors using sub-Nyquist sampling.


title: Training-Free Non-Intrusive Load Monitoring of Electric Vehicle Charging with Low Sampling Rate

abstract: Non-intrusive load monitoring (NILM) is an important topic in smart-grid and smart-home. Many energy disaggregation algorithms have been proposed to detect various individual appliances from one aggregated signal observation. However, few works studied the energy disaggregation of plug-in electric vehicle (EV) charging in the residential environment since EVs charging at home has emerged only recently. Recent studies showed that EV charging has a large impact on smart-grid especially in summer. Therefore, EV charging monitoring has become a more important and urgent missing piece in energy disaggregation. In this paper, we present a novel method to disaggregate EV charging signals from aggregated real power signals. The proposed method can effectively mitigate interference coming from air-conditioner (AC), enabling accurate EV charging detection and energy estimation under the presence of AC power signals. Besides, the proposed algorithm requires no training, demands a light computational load, delivers high estimation accuracy, and works well for data recorded at the low sampling rate 1/60 Hz. When the algorithm is tested on real-world data recorded from 11 houses over about a whole year (total 125 months worth of data), the averaged error in estimating energy consumption of EV charging is 15.7 kwh/month (while the true averaged energy consumption of EV charging is 208.5 kwh/month), and the averaged normalized mean square error in disaggregating EV charging load signals is 0.19.


title: Local Rank Modulation for Flash Memories II

abstract: Local rank modulation scheme was suggested recently for representing information in flash memories in order to overcome drawbacks of rank modulation. For $0 < s\leq t\leq n$ with $s$ divides $n$, an $(s,t,n)$-LRM scheme is a local rank modulation scheme where the $n$ cells are locally viewed cyclically through a sliding window of size $t$ resulting in a sequence of small permutations which requires less comparisons and less distinct values. The gap between two such windows equals to $s$. In this work, encoding, decoding, and asymptotic enumeration of the $(1,3,n)$-LRM scheme is studied. The techniques which are suggested have some generalizations for $(1,t,n)$-LRM, $t > 3$, but the proofs will become more complicated. The enumeration problem is presented also as a purely combinatorial problem. Finally, we prove the conjecture that the size of a constant weight $(1,2,n)$-LRM Gray code with weight two is at most $2n$.


title: An Elementary Proof That Rationally Isometric Quadratic Forms Are Isometric

abstract: Let $R$ be a valuation ring with fraction field $K$ and $2\in R^\times$. We give an elementary proof of the following known result: Two unimodular quadratic forms over $R$ are isometric over $K$ if and only if they are isometric over $R$. Our proof does not use Witt's Cancelation Theorem and yields an explicit algorithm to construct an isometry over $R$ from a given isometry over $K$. The statement actually holds for hermitian forms over valuated involutary division rings, provided mild assumptions. A python implementation of the algorithm derived from the proof can be found on the author's home page.


title: The Betti numbers for a family of solvable Lie algebras

abstract: We give a characterization of symplectic quadratic Lie algebras that their Lie algebra of inner derivations has an invertible derivation. A family of symplectic quadratic Lie algebras is introduced to illustrate this situation. Finally, we calculate explicitly the Betti numbers of a family of solvable Lie algebras in two ways: using the cohomology of quadratic Lie algebras and applying a Pouseele's result on extensions of the one-dimensional Lie algebra by Heisenberg Lie algebras


title: Consequences of Deformation of the Heisenberg Algebra

abstract: In this paper we will demonstrate that like the existence of a minimum measurable length, the existence of a maximum measurable momentum, also influence all quantum mechanical systems. Beyond the simple one dimensional case, the existence of a maximum momentum will induce non-local corrections to the first quantized Hamiltonian. However, these non-local corrections can be effectively treated as local corrections by using the theory of harmonic extensions of functions. We will also analyses the second quantization of this deformed first quantized theory. Finally, we will analyses the gauge symmetry corresponding to this deformed theory.


title: Abelian and non-abelian cohomology

abstract: We place the representation variety in the broader context of abelian and nonabelian cohomology. We outline the equivalent constructions of the moduli spaces of flat bundles, of smooth integrable connections, and of holomorphic integrable connections over a compact Kaehler manifold. In addition, we describe the moduli space of Higgs bundles and how it relates to the representation variety. We attempt to avoid abstraction, but strive to present and clarify the unifying ideas underlying the theory.


title: The second Hilbert coefficients and the sectional genera of ideals

abstract: This paper gives an upper bound of the second Hilbert coefficients ${\mathrm{e}}_2(I)$ of $\mathfrak{m}$-primary ideals $I$ in a Cohen-Macaulay local ring $(A,{\mathfrak{m}})$ in terms of the sectional genera of $I$.


title: Anisotropic plasma at finite $U(1)$ chemical potential

abstract: We present a type IIB supergravity solution dual to a spatially anisotropic $\mathcal{N}=4$ super Yang-Mills plasma at finite $U(1)$ chemical potential and finite temperature. The effective five-dimensional gravitational theory is a consistent Einstein-Maxwell-dilaton-Axion truncation of the gauged supergravity. We obtain the solutions both numerically and analytically. We study the phase structure and thermodynamic instabilities of the solution, and find new instabilities independent of the renormalization scheme.


title: Clustering via Mode Seeking by Direct Estimation of the Gradient of a Log-Density

abstract: Mean shift clustering finds the modes of the data probability density by identifying the zero points of the density gradient. Since it does not require to fix the number of clusters in advance, the mean shift has been a popular clustering algorithm in various application fields. A typical implementation of the mean shift is to first estimate the density by kernel density estimation and then compute its gradient. However, since good density estimation does not necessarily imply accurate estimation of the density gradient, such an indirect two-step approach is not reliable. In this paper, we propose a method to directly estimate the gradient of the log-density without going through density estimation. The proposed method gives the global solution analytically and thus is computationally efficient. We then develop a mean-shift-like fixed-point algorithm to find the modes of the density for clustering. As in the mean shift, one does not need to set the number of clusters in advance. We empirically show that the proposed clustering method works much better than the mean shift especially for high-dimensional data. Experimental results further indicate that the proposed method outperforms existing clustering methods.


title: Using Covert Topological Information for Defense Against Malicious Attacks on DC State Estimation

abstract: Accurate state estimation is of paramount importance to maintain the power system operating in a secure and efficient state. The recently identified coordinated data injection attacks to meter measurements can bypass the current security system and introduce errors to the state estimates. The conventional wisdom to mitigate such attacks is by securing meter measurements to evade malicious injections. In this paper, we provide a novel alternative to defend against false-data injection attacks using covert power network topological information. By keeping the exact reactance of a set of transmission lines from attackers, no false data injection attack can be launched to compromise any set of state variables. We first investigate from the attackers' perspective the necessary condition to perform injection attack. Based on the arguments, we characterize the optimal protection problem, which protects the state variables with minimum cost, as a well-studied Steiner tree problem in a graph. Besides, we also propose a mixed defending strategy that jointly considers the use of covert topological information and secure meter measurements when either method alone is costly or unable to achieve the protection objective. A mixed integer linear programming (MILP) formulation is introduced to obtain the optimal mixed defending strategy. To tackle the NP-hardness of the problem, a tree pruning-based heuristic is further presented to produce an approximate solution in polynomial time. The advantageous performance of the proposed defending mechanisms is verified in IEEE standard power system testcases.


title: Robust antiferromagnetism preventing superconductivity in pressurized Ba0.61K0.39Mn2Bi2

abstract: Here, we report that K-doped BaMn2Bi2 shows no experimental evidence of superconductivity down to 1.5 K under pressures up to 35.6 GPa, however, a tetragonal to an orthorhombic phase transition is observed at pressure of 20 GPa. Theoretical calculations for the tetragonal and orthorhombic phases, on basis of our high-pressure XRD data, find that the AFM order is robust in both of the phases in pressurized Ba0.61K0.39Mn2Bi2. Our experimental and theoretical results suggest that the K-doped BaMn2Bi2 belongs to a strong Hunds AFM metal with a hybridization of localized spin electrons and itinerant electrons, and that its robust AFM order essentially prevents the emergence of superconductivity.


title: On Multicorns and Unicorns II: Bifurcations in Spaces of Antiholomorphic Polynomials

abstract: The multicorns are the connectedness loci of unicritical antiholomorphic polynomials $\bar{z}^d + c$. We investigate the structure of boundaries of hyperbolic components: we prove that the structure of bifurcations from hyperbolic components of even period is as one would expect for maps that depend holomorphically on a complex parameter (for instance, as for the Mandelbrot set; in this setting, this is a non-obvious fact), while the bifurcation structure at hyperbolic components of odd period is very different. In particular, the boundaries of odd period hyperbolic components consist only of parabolic parameters, and there are bifurcations between hyperbolic components along entire arcs, but only of bifurcation ratio $2$. We also count the number of hyperbolic components of any period of the multicorns. Since antiholomorphic polynomials depend only real-analytically on the parameters, most of the techniques used in this paper are quite different from the ones used to prove the corresponding results in a holomorphic setting.


title: Tchebychev Polynomial Approximations for $m^{th}$ Order Boundary Value Problems

abstract: Higher order boundary value problems (BVPs) play an important role modeling various scientific and engineering problems. In this article we develop an efficient numerical scheme for linear $m^{th}$ order BVPs. First we convert the higher order BVP to a first order BVP. Then we use Tchebychev orthogonal polynomials to approximate the solution of the BVP as a weighted sum of polynomials. We collocate at Tchebychev clustered grid points to generate a system of equations to approximate the weights for the polynomials. The excellency of the numerical scheme is illustrated through some examples.


title: Optimal receiver for binary coherent signals

abstract: One of the most fascinating aspects of quantum mechanics is the principle impossibility of deterministic errorless discrimination of nonorthogonal signals, such as coherent states. On the one hand, it prevents perfect cloning of quantum states and enables secure communication. On the other hand, it makes a grand challenge to reach the ultimate measurement precision. Although the minimum possible error rate (the Helstrom bound) has been known for almost five decades, there is no practical way to achieve it. Developing the realistic optimal measurement strategies to attain the Helstrom bound is of utmost importance for high-precision applications, long-distance free-space and optical fiber communication, gravitational wave detection, optical sensing in biology and medicine, to name a few. In this work, we show an optimal receiver for coherent states which admits a relatively simple technological implementation. The receiver is based on multichannel splitting of the signal, followed by feed-forward signal displacement and photon detection.


title: Sustaining IT PMOs during Cycles of Global Recession

abstract: Growth in the number of PMOs established by the industry over last decade and ever growing body of literature on PMO related research in academia is a clear indication that there is very clear interest of researchers, practitioners and industries across the globe to understand and explore value propositions of PMO. However, there is still a lack of consensus on many critical aspects of PMOs. While there are many PMOs being established, but there are also many being closed and disbanded, which is definitely a matter of concern. In industry environment, a narrow majority of PMOs are well-regarded by their organizations and are seen as contributing business value, many of the others are still struggling to show value for money and some are failing, causing a high mortality rate among PMOs. This paper is the result of a study undertaken to get a deeper understanding of factors that may be causing mortality and failure of PMOs. Post Implementation Reviews of 4-failed & 3-challenged PMOs in IT-Industry were carried out with concerned Project Managers & PMO-staff, using grounded theory research method, with support from the concerned enterprise from IT-Industry.


title: $n$-widths and Approximation theory on Compact Riemannian Manifolds

abstract: We determine upper asymptotic estimates of Kolmogorov and linear $n$-widths of unit balls in Sobolev and Besov norms in $L_{p}$-spaces on compact Riemannian manifolds. The proofs rely on estimates for the near-diagonal localization of the kernels of elliptic operators. We also summarize some of our previous results about approximations by eigenfunctions of elliptic operators on compact homogeneous manifolds.


title: Second derivative estimates for uniformly elliptic operators on Riemannian manifolds

abstract: In this paper, we obtain a uniform $W^{2,\varepsilon}$-estimate of solutions to the fully nonlinear uniformly elliptic equations on Riemannian manifolds with a lower bound of sectional curvature using the ABP method.


title: Multiresolution analysis on compact Riemannian manifolds

abstract: In the chapter "Multiresolution Analysis on Compact Riemannian Manifolds" Isaac Pesenson describes multiscale analysis, sampling, interpolation and approximation of functions defined on manifolds. His main achievements are: construction on manifolds of bandlimited and space-localized frames which have Parseval property and construction of variational splines on manifolds. Such frames and splines enable multiscale analysis on arbitrary compact manifolds, and they already found a number of important applications (statistics, CMB, crystallography) related to such manifolds as two-dimensional sphere and group of its rotations.


title: Signal of doubly charged Higgs at $e^+e^-$ colliders

abstract: Masses and signals of the production of Doubly charged Higgses(DCH) in the framework of the supersymmetric reduced minimal 3-3-1 model(SUSYRM331) are investigated. In the DCH sector, weprove that there always exists a region of the parameter space where the mass of the lightest DCH is in order of $\mathcal{O}(100)$ GeV even when all other new particles are very heavy.The lightest DCH mainly decays to two same-sign leptons while the dominated decay channels of the heavy DCHs are those decay to heavy particles. We analyze each production cross section for $e^+e^- \rightarrow H^{++} H^{--}$ as a function of a few kinematic variables, which are useful to discuss the creation of the DCHs in the $e^+e^-$ colliders as a signal of new physics beyond the Standard Model. The numerical study shows that the cross sections for creating the lightest DCH can reach values of few pbs. The two other DCHs are too heavy, beyond the observable range of experiments. The lightest DCH may be detected by the International Linear Collider (ILC) or the Compact LInear Collider (CLIC) by searching its decay to a same-sign charged lepton pair.


title: Stochastic nonlinear Schr\"odinger equations

abstract: This paper is devoted to the well-posedness of stochastic nonlinear Schr\"odinger equations in the energy space H1(Rd), which is a natural continuation of our recent work [1]. We consider both focusing and defocusing nonlinearities and prove global well-posedness in H1(Rd), including also the pathwise continuous dependence on initial conditions, with exponents exactly the same as in the deterministic case. In particular, this work improves earlier results in [4]. Moreover, the local existence, uniqueness and blowup alternative are also established for the energy-critical case. The approach presented here is mainly based on the rescaling approach already used in [1] to study the L2 case and also on the Strichartz estimates established in [12] for large perturbations of the Laplacian.


title: Existence of minimizers for Kohn-Sham within the Local Spin Density Approximation

abstract: The purpose of this article is to extend the work by Anantharaman and Canc\`es [1], and prove the existence of minimizers for the spin-polarized Kohn-Sham model in the presence of a magnetic field within the local spin density approximation. We show that for any magnetic field that vanishes at infinity, the existence of minimizers is ensured for neutral or positively charged systems. The proof relies on classical concentration-compactness techniques.


title: The Exponentially Faster Stick-Slip Dynamics of the Peeling of an Adhesive Tape

abstract: The stick-slip dynamics is considered from the nonlinear differential-algebraic equation (DAE) point of view and the peeling dynamics is shown to be a switching differential index DAE model. In the stick-slip regime with bifurcations, the differential index can be arbitrarily high. The time scale of the peeling velocity, the algebraic variable, in this regime is shown to be exponentially faster compared to the angular velocity of the spool and/or the stretch rate of the tape. A homogenization scheme for the peeling velocity which is characterized by the bifurcations is discussed and is illustrated with numerical examples.


title: HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions

abstract: We present new analytic theory and radiative transfer computations for the atomic to molecular (HI-to-H2) transitions, and the build-up of atomic-hydrogen (HI) gas columns, in optically thick interstellar clouds, irradiated by far-ultraviolet photodissociating radiation fields. We derive analytic expressions for the total HI column densities for (1D) planar slabs, for beamed or isotropic radiation fields, from the weak- to strong-field limits, for gradual or sharp atomic to molecular transitions, and for arbitrary metallicity. Our expressions may be used to evaluate the HI column densities as functions of the radiation field intensity and the H2-dust-limited dissociation flux, the hydrogen gas density, and the metallicity-dependent H2 formation rate-coefficient and far-UV dust-grain absorption cross-section. We make the distinction between "HI-dust" and "H2-dust" opacity, and we present computations for the "universal H2-dust-limited effective dissociation bandwidth". We validate our analytic formulae with Meudon PDR code computations for the HI-to-H2 density profiles, and total HI column densities. We show that our general 1D formulae predict HI columns and H2 mass fractions that are essentially identical to those found in more complicated (and approximate) spherical (shell/core) models. We apply our theory to compute H2 mass fractions and star-formation thresholds for individual clouds in self-regulated galaxy disks, for a wide range of metallicities. Our formulae for the HI columns and H2 mass fractions may be incorporated into hydrodynamics simulations for galaxy evolution.


title: The predictable degree property, column reducedness, and minimality in multidimensional convolutional coding

abstract: Higher-dimensional analogs of the predictable degree property and column reducedness are defined, and it is proved that the two properties are equivalent. It is shown that every multidimensional convolutional code has, what is called, a minimal reduced polynomial resolution. It is uniquely determined (up to isomorphism) and leads to a number of important integer invariants of the code generalizing classical Forney's indices.


title: Anomalous Power Law Decay in Solvation Dynamics of DNA: A Mode Coupling Theory Analysis of Ion Contribution

abstract: Several time domain fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Lastly, we propose several experiments to check the prediction of the present theoretical work.


title: Ample Group Action on AS-regular Algebras and Noncommutative Graded Isolated Singularities

abstract: In this paper, we introduce a notion of ampleness of a group action $G$ on a right noetherian graded algebra $A$, and show that it is strongly related to the notion of $A^G$ to be a graded isolated singularity introduced by the second author of this paper. Moreover, if $S$ is a noetherian AS-regular algebra and $G$ is a finite ample group acting on $S$, then we will show that ${\mathcal D}^b(\operatorname{tails} S^G)\cong {\cal D}^b(\operatorname{mod} \nabla S*G)$ where $\nabla S$ is the Beilinson algebra of $S$. We will also explicitly calculate a quiver $Q_{S, G}$ such that ${\mathcal D}^b(\operatorname{tails} S^G)\cong {\mathcal D}^b(\operatorname{mod} kQ_{S, G})$ when $S$ is of dimension 2.


title: A Comparative Study of Protein Unfolding in Aqueous Urea and DMSO Solutions: Surface Polarity, Solvent Specificity and Sequence of Secondary Structure Melting

abstract: Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here we alter the solvent environment in a controlled manner by choosing two different co-solvents of water, urea and dimethyl sulphoxide (DMSO), and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha-helices and beta-sheets in these two solvents. For example, at 8M urea solution, beta-sheet parts of a protein is found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8M DMSO solution unfolds alpha helices first, followed by the separation of beta-sheets for majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solution demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and non-polar residues of a protein and the mode of molecular action of a solvent on that protein.


title: Electromagnetic response of quark-gluon plasma in heavy-ion collisions

abstract: We study the electromagnetic response of the quark-gluon plasma in AA-collisions at RHIC and LHC energies for a realistic space-time evolution of the plasma fireball. We demonstrate that for a realistic electric conductivity the electromagnetic response of the plasma is in a quantum regime when the induced electric current does not generate a classical electromagnetic field, and can only lead to a rare emission of single photons.


title: Congruence identities of regularized multiple zeta values involving a pair of index sets

abstract: Riemann zeta values are generalized to multiple zeta values (MZVs) by use of nested sums, and MZVs are generalized to regularized multiple zeta values (RMZVs) by regularization of divergent infinite series. In the present paper, we prove congruence identities of RMZVs of depth $n$ involving a pair of index sets; the congruence relation is given by the vector space spanned by MZVs of depth $n-1$ and products of MZVs. We also obtain a proof of the parity result, and a congruence sum formula for MZVs.


title: Sobre a viagem de Enrico Fermi ao Brasil em 1934

abstract: Enrico Fermi was one of the greater physicists of the XX century. In 1934, he gave several lectures in Brazil. Invited by Theodoro Ramos to work in S\~ao Paulo, he preferred to stay in Rome and went to the USA in 1938. However, Fermi recommended Gleb Wataghin to come in his place. Wataghin made history in Brazil, becoming one of the first Professors of the future S\~ao Paulo University. Besides his relevance to the History of Science, Fermi eventually leaved an indelible mark on the creation and institutionalization of national scientific research due to the indication of Wataghin. Despite this fact, very little is known about Fermi's trip to Brazil. This work tries to reconstruct the fullest possible steps of the famous Italian physicist in our lands.


title: A Spectral Model of Turnover Reduction

abstract: We give a simple explicit formula for turnover reduction when a large number of alphas are traded on the same execution platform and trades are crossed internally. We model turnover reduction via alpha correlations. Then, for a large number of alphas, turnover reduction is related to the largest eigenvalue and the corresponding eigenvector of the alpha correlation matrix.


title: Convex geodesic bicombings and hyperbolicity

abstract: A geodesic bicombing on a metric space selects for every pair of points a geodesic connecting them. We prove existence and uniqueness results for geodesic bicombings satisfying different convexity conditions. In combination with recent work by the second author on injective hulls, this shows that every word hyperbolic group acts geometrically on a proper, finite dimensional space X with a unique (hence equivariant) convex geodesic bicombing of the strongest type. Furthermore, the Gromov boundary of X is a Z-set in the closure of X, and the latter is a metrizable absolute retract, in analogy with the Bestvina--Mess theorem on the Rips complex.


title: Mass Generation in QCD -- Oscillating Quarks and Gluons

abstract: The present lecture is devoted to embedding the approximate genuine harmonic oscillator structure of valence q qbar mesons and in more detail the q q q configurations for up,down,strange - flavored baryons in Q C D for three light flavors of quark.


title: The Controlled Thermodynamic Integral for Bayesian Model Comparison

abstract: Bayesian model comparison relies upon the model evidence, yet for many models of interest the model evidence is unavailable in closed form and must be approximated. Many of the estimators for evidence that have been proposed in the Monte Carlo literature suffer from high variability. This paper considers the reduction of variance that can be achieved by exploiting control variates in this setting. Our methodology is based on thermodynamic integration and applies whenever the gradient of both the log-likelihood and the log-prior with respect to the parameters can be efficiently evaluated. Results obtained on regression models and popular benchmark datasets demonstrate a significant and sometimes dramatic reduction in estimator variance and provide insight into the wider applicability of control variates to Bayesian model comparison.


title: Two-spinor geometry and gauge freedom

abstract: Gauge freedom in quantum particle physics is shown to arise in a natural way from the geometry of two-spinors (Weyl spinors). Various related mathematical notions are reviewed, and a special ansatz of the kind "the system defines the geometry" is discussed in connection with the stated results.


title: Correlated Jamming in a Joint Source Channel Communication System

abstract: We study correlated jamming in joint source-channel communication systems. An i.i.d. source is to be communicated over a memoryless channel in the presence of a correlated jammer with non-causal knowledge of user transmission. This user-jammer interaction is modeled as a zero sum game. A set of conditions on the source and the channel is provided for the existence of a Nash equilibrium for this game, where the user strategy is uncoded transmission and the jammer strategy is i.i.d jamming. This generalizes a well-known example of uncoded communication of a Gaussian sources over Gaussian channels with additive jamming. Another example, of a Binary Symmetric source over a Binary Symmetric channel with jamming, is provided as a validation of this result.


title: Pinned modes in two-dimensional lossy lattices with local gain and nonlinearity

abstract: We introduce a system with one or two amplified nonlinear sites ("hot spots", HSs) embedded into a two-dimensional linear lossy lattice. The system describes an array of evanescently coupled optical or plasmonic waveguides, with gain applied at selected HS cores. The subject of the analysis is discrete solitons pinned to the HSs. The shape of the localized modes is found in quasi-analytical and numerical forms, using a truncated lattice for the analytical consideration. Stability eigenvalues are computed numerically, and the results are supplemented by direct numerical simulations. In the case of self-focusing nonlinearity, the modes pinned to a single HS are stable or unstable when the nonlinearity includes the cubic loss or gain, respectively. If the nonlinearity is self-defocusing, the unsaturated cubic gain acting at the HS supports stable modes in a small parametric area, while weak cubic loss gives rise to a bistability of the discrete solitons. Symmetric and antisymmetric modes pinned to a symmetric set of two HSs are considered too.


title: Topological Dynamics of Closed Subgroups of $S_\infty$

abstract: For $G$ a closed subgroup of $S_{\infty}$, we provide an explicit characterization of the greatest $G$-ambit. Using this, we provide a precise characterization of when $G$ has metrizable universal minimal flow. In particular, each such instance fits into the framework of metrizable flows developed by Kechris-Pestov-Todorcevic and Nguyen Van Th\'e; as a consequence, each $G$ with metrizable universal minimal flow has the generic point property.


title: Electrodynamic phenomena induced by a dark fluid: Analogs of pyromagnetic, piezoelectric, and striction effects

abstract: We establish a new model of coupling between a cosmic dark fluid and electrodynamic systems, based on an analogy with effects of electric and magnetic striction, piezo-electricity and piezo-magnetism, pyro-electricity and pyro-magnetism, which appear in classical electrodynamics of continuous media. Extended master equations for electromagnetic and gravitational fields are derived using Lagrange formalism. A cosmological application of the model is considered, and it is shown that a striction-type interaction between the dark energy (the main constituent of the dark fluid) and electrodynamic system provides the universe history to include the so-called unlighted epochs, during which electromagnetic waves can not propagate and thus can not scan the universe interior.


title: Investigating CXOU J163802.6-471358: a new pulsar wind nebula in the Norma region?

abstract: We present the first analysis of the extended source CXOU J163802.6--471358, which was discovered serendipitously during the {\em Chandra} X-ray survey of the Norma region of the Galactic spiral arms. The X-ray source exhibits a cometary appearance with a point source and an extended tail region. The complete source spectrum is fitted well with an absorbed power law model and jointly fitting the {\em Chandra} spectrum of the full source with one obtained from an archived {\em XMM-Newton} observation results in best fit parameters $N_{\rm H}$ $=1.5^{+0.7}_{-0.5}\times10^{23} \text{cm}{^{-2}}$ and $\Gamma=1.1^{+0.7}_{-0.6}$ (90$%$ confidence uncertainties). The unabsorbed luminosity of the full source is then $L_X\sim 4.8\times10^{33}d_{10}^2$ergs s$^{-1}$ with $d_{10}=d/10$kpc, where a distance of 10 kpc is a lower bound inferred from the large column density. The radio counterpart found for the source using data from the Molonglo Galactic Plane Survey epoch-2 (MGPS-2) shows an elongated tail offset from the X-ray emission. No infrared counterpart was found. The results are consistent with the source being a previously unknown pulsar driving a bow shock through the ambient medium.


title: Writing on a Dirty Paper in the presence of Jamming

abstract: In this paper, the problem of writing on a dirty paper in the presence of jamming is examined. We consider an AWGN channel with an additive white Gaussian state and an additive adversarial jammer. The state is assumed to be known non-causally to the encoder and the jammer but not to the decoder. The capacity of the channel in the presence of a jammer is determined. A surprising result that this capacity is equal to the capacity of a relaxed version of the problem, where the state is also known non-causally to the decoder, is proved.


title: Closing in on the Tip of the CMSSM Stau Coannihilation Strip

abstract: Near the tip of the stau coannihilation strip in the CMSSM with a neutralino LSP, the astrophysical cold dark matter density constraint forces the stau-neutralino mass difference to be small. If this mass difference is smaller than the tau mass, the stau may decay either in the outer part of an LHC detector - the `disappearing track' signature - or be sufficiently long-lived to leave the detector before decaying - the long-lived massive charged-particle signature. We combine searches for these signatures with conventional missing transverse energy searches during LHC Run 1, identifying the small remaining parts of the CMSSM stau coannihilation strip region that have not yet been excluded, and discussing how they may be explored during Run 2 of the LHC.


title: Rapid prototyping for sling design optimization

abstract: This paper deals with combination of two modern engineering methods in order to optimise the shape of a representative casting product. The product being analysed is a sling, which is used to attach pulling rope in timber transportation. The first step was 3D modelling and static stress/strain analysis using CAD/CAE software NX4. The slinger shape optimization was performed using Traction method, by means of software Optishape-TS. To define constraints for shape optimization, FEA software FEMAP was used. The mould pattern with optimized 3D shape was then prepared using Fused Deposition Modelling (FDM) Rapid prototyping method. The sling mass decreased by 20%, while signifficantly better stress distribution was achieved, with maximum stress 3.5 times less than initial value. The future researches should use 3D scanning technology in order to provide more accurate 3D model of initial part. Results of this research can be used by toolmakers in order to engage FEA/RP technology to design and manufacture lighter products with acceptable stress distribution.


title: Hints of a second explosion (a quark nova) in Cassiopeia A Supernova

abstract: We show that the explosive transition of the neutron star (NS) to a quark star (QS) (a Quark Nova) in Cassiopeia A (Cas A) a few days following the SN proper can account for several of the puzzling kinematic and nucleosynthetic features observed. The observed decoupling between Fe and 44Ti and the lack of Fe emission within 44Ti regions is expected in the QN model owing to the spallation of the inner SN ejecta by the relativistic QN neutrons. Our model predicts the 44Ti to be more prominent to the NW of the central compact object (CCO) than in the SE and little of it along the NE-SW jets, in agreement with NuStar observations. Other intriguing features of Cas A such as the lack of a pulsar wind nebula (PWN) and the reported a few percent drop of the CCO temperature over a period of 10 years are also addressed.


title: Nonhomogeneous Place-Dependent Markov Chains, Unsynchronised AIMD, and Network Utility Maximization

abstract: We present a solution of a class of network utility maximization (NUM) problems using minimal communication. The constraints of the problem are inspired less by TCP-like congestion control but by problems in the area of internet of things and related areas in which the need arises to bring the behavior of a large group of agents to a social optimum. The approach uses only intermittent feedback, no inter-agent communication, and no common clock. The proposed algorithm is a combination of the classical AIMD algorithm in conjunction with a simple probabilistic rule for the agents to respond to a capacity signal. This leads to a nonhomogeneous Markov chain and we show almost sure convergence of this chain to the social optimum.


title: Multi-Target Regression via Random Linear Target Combinations

abstract: Multi-target regression is concerned with the simultaneous prediction of multiple continuous target variables based on the same set of input variables. It arises in several interesting industrial and environmental application domains, such as ecological modelling and energy forecasting. This paper presents an ensemble method for multi-target regression that constructs new target variables via random linear combinations of existing targets. We discuss the connection of our approach with multi-label classification algorithms, in particular RA$k$EL, which originally inspired this work, and a family of recent multi-label classification algorithms that involve output coding. Experimental results on 12 multi-target datasets show that it performs significantly better than a strong baseline that learns a single model for each target using gradient boosting and compares favourably to multi-objective random forest approach, which is a state-of-the-art approach. The experiments further show that our approach improves more when stronger unconditional dependencies exist among the targets.


title: Hypersurface foliation approach to renormalization of ADM formulation of gravity

abstract: We carry out ADM splitting in the Lagrangian formulation and establish a procedure in which (almost) all of the unphysical components of the metric are removed by using the 4D diffeomorphism and the measure-zero 3D symmetry. The procedure introduces a constraint that corresponds to the Hamiltonian constraint of the Hamiltonian formulation, and its solution implies that the 4D dynamics admits an effective description through 3D hypersurface physics. As far as we can see, our procedure implies potential renormalizability of {the ADM formulation of} 4D Einstein gravity for which a complete gauge-fixing in the ADM formulation and hypersurface foliation of geometry are the key elements. If true, this implies that the alleged unrenormalizability of 4D Einstein gravity may be due to the presence of the unphysical fields. The procedure can straightforwardly be applied to quantization around a flat background; the Schwarzschild case seems more subtle. We discuss a potential limitation of the procedure when applying it to explicit time-dependent backgrounds.


title: Cosmological Perturbations Through a Non-Singular Ghost-Condensate/Galileon Bounce

abstract: We study the propagation of super-horizon cosmological perturbations in a non-singular bounce spacetime. The model we consider combines a ghost condensate with a Galileon term in order to induce a ghost-free bounce. Our calculation is performed in harmonic gauge, which ensures that the linearized equations of motion remain well-defined and non-singular throughout. We find that, despite the fact that near the bounce the speed of sound becomes imaginary, super-horizon curvature perturbations remain essentially constant across the bounce. In fact, we show that there is a time close to the bounce where curvature perturbations of all wavelengths are required to be momentarily exactly constant. We relate our calculations to those performed in other gauges, and comment on the relation to previous results in the literature.


title: Directional Cell Discovery in Millimeter Wave Cellular Networks

abstract: The acute disparity between increasing bandwidth demand and available spectrum, has brought millimeter wave (mmW) bands to the forefront of candidate solutions for the next-generation cellular networks. Highly directional transmissions are essential for cellular communication in these frequencies to compensate for high isotropic path loss. This reliance on directional beamforming, however, complicates initial cell search since the mobile and base station must jointly search over a potentially large angular directional space to locate a suitable path to initiate communication. To address this problem, this paper proposes a directional cell discovery procedure where base stations periodically transmit synchronization signals, potentially in time-varying random directions, to scan the angular space. Detectors for these signals are derived based on a Generalized Likelihood Ratio Test (GLRT) under various signal and receiver assumptions. The detectors are then simulated under realistic design parameters and channels based on actual experimental measurements at 28~GHz in New York City. The study reveals two key findings: (i) digital beamforming can significantly outperform analog beamforming even when the digital beamforming uses very low quantization to compensate for the additional power requirements; and (ii) omni-directional transmissions of the synchronization signals from the base station generally outperforms random directional scanning.


title: Computing periods of rational integrals

abstract: A period of a rational integral is the result of integrating, with respect to one or several variables, a rational function over a closed path. This work focuses particularly on periods depending on a parameter: in this case the period under consideration satisfies a linear differential equation, the Picard-Fuchs equation. I give a reduction algorithm that extends the Griffiths-Dwork reduction and apply it to the computation of Picard-Fuchs equations. The resulting algorithm is elementary and has been successfully applied to problems that were previously out of reach.


title: Congruences involving product of intervals and sets with small multiplicative doubling modulo a prime and applications

abstract: In the present paper we obtain new upper bound estimates for the number of solutions of the congruence $$ x\equiv y r\pmod p;\quad x,y\in \mathbb{N},\quad x,y\le H,\quad r\in\cU, $$ for certain ranges of $H$ and $|\cU|$, where $\cU$ is a subset of the field of residue classes modulo $p$ having small multiplicative doubling. We then use this estimate to show that the number of solutions of the congruence $$ x^n\equiv \lambda\pmod p; \quad x\in \N, \quad L<x<L+p/n, $$ is at most $p^{\frac{1}{3}-c}$ uniformly over positive integers $n, \lambda$ and $L$, for some absolute constant $c>0$. This implies, in particular, that if $f(x)\in \Z[x]$ is a fixed polynomial without multiple roots in $\C$, then the congruence $ x^{f(x)}\equiv 1\pmod p, \,x\in \mathbb{N}, \,x\le p,$ has at most $p^{\frac{1}{3}-c}$ solutions as $p\to\infty$, improving some recent results of Kurlberg, Luca and Shparlinski and of Balog, Broughan and Shparlinski. We use our results to show that almost all the residue classes modulo $p$ can be represented in the form $xg^y \pmod p$ with positive integers $x<p^{5/8+\varepsilon}$ and $y<p^{3/8}$. Here $g$ denotes a primitive root modulo $p$. We also prove that almost all the residue classes modulo $p$ can be represented in the form $xyzg^t \pmod p$ with positive integers $x,y,z,t<p^{1/4+\varepsilon}$.


title: Sparsity-Exploiting Moment-Based Relaxations of the Optimal Power Flow Problem

abstract: Convex relaxations of non-convex optimal power flow (OPF) problems have recently attracted significant interest. While existing relaxations globally solve many OPF problems, there are practical problems for which existing relaxations fail to yield physically meaningful solutions. This paper applies moment relaxations to solve many of these OPF problems. The moment relaxations are developed from the Lasserre hierarchy for solving generalized moment problems. Increasing the relaxation order in this hierarchy results in "tighter" relaxations at the computational cost of larger semidefinite programs. Low-order moment relaxations are capable of globally solving many small OPF problems for which existing relaxations fail. By exploiting sparsity and only applying the higher-order relaxation to specific buses, global solutions to larger problems are computationally tractable through the use of an iterative algorithm informed by a heuristic for choosing where to apply the higher-order constraints. With standard semidefinite programming solvers, the algorithm globally solves many test systems with up to 300 buses for which the existing semidefinite relaxation fails to yield globally optimal solutions.


title: Metabolic fluxes and value production

abstract: Metabolic fluxes in cells are governed by physical, physiological, and economic principles. Here I assume an optimal allocation of enzyme resources and postulate a general principle for metabolism: each enzyme must convert less valuable into more valuable metabolites to justify its own cost. The "values", called economic potentials, describe the individual contributions of metabolites to cell fitness. Local value production implies that the cost of an enzyme must be balanced by a benefit, given by the economic potential difference the catalysed reaction multiplied by the flux. Flux profiles that satisfy this principle - i.e. for which consistent potentials can be found - are called economical. Economical fluxes must lead from lower to higher economic potentials, so certain flux cycles are incompatible with any choice of economic potentials and can be excluded. To obtain economical flux profiles, non-beneficial local patterns, called futile motifs, can be systematically removed from a given flux distribution. The principle of local value production resembles thermodynamic principles and complements them in models. Here I describe a modelling framework called Value Balance Analysis (VBA) that uses the two principles and yields the same solution as enzyme cost minimisation (in kinetic models) and flux cost minimisation (in FBA). Given an economical flux distribution, kinetic models in enzyme-optimal states and with these fluxes can be constructed systematically. VBA justifies the principle of minimal fluxes and the exclusion of futile cycles, predicts enzymes that could be plausible targets for regulation, provides criteria for the usage of enzymes and pathways, and explains the choice between high-yield and low-yield flux modes.


title: Scaling Invariance of Density Functionals

abstract: Based on the homogeneity ($F[n_{\lambda m}]=\lambda^{p(m)}F[n]$) and invariance ($F[n_{\lambda m_0}]=F[n]$) properties of a functional of the electron density under uniform scaling of the coordinates in the density $n_{\lambda m}(\mathbf{r})=\lambda^{m} n(\lambda\mathbf{r}),\,(\lambda\in\mathbb{R}^+,\, m\in\mathbb{R})$, it is proven that homogeneity implies invariace and therefore all homogeneous scaling functionals have the representation $F[n]=\frac{m-m_0}{p(m)} \int_V\,\frac{\delta F[n]}{\delta n(\mathbf{r})}\,n(\mathbf{r})\,d^3r$. Also, the homogeneity ($p(m)$) and invariant ($m_0$) degrees of density functionals related to the Kohn-Sham theory are calculated. Besides, it is shown that the functional density and the electron density itself satisfy the general equation representing the local scaling invariance of a functional $\lambda \frac{d}{d\lambda} f([n_{\lambda m_0}],\mathbf{r},\mathbf{r'}) = \sum_{i=1}^3 \frac{d}{d x_i} [ x_i f([n_{\lambda m_0}],\mathbf{r},\mathbf{r'}) ] + \sum_{j=1}^3 \frac{d}{d x_j'} [ x_j' f([n_{\lambda m_0}],\mathbf{r},\mathbf{r'}) ] $. The equation simplifies for cases where the functional density depends only on the density and/or its gradient, and general forms of the solutions are provided, in particular for the non-interacting kinetic energy density is shown to take the form $t_s(n,\nabla n)= n(\mathbf{r})^{3} g[ \frac{\partial_{x_1} n(\mathbf{r})}{n(\mathbf{r})^2}, \frac{\partial_{x_2} n(\mathbf{r})}{n(\mathbf{r})^2}, \frac{\partial_{x_3} n(\mathbf{r})}{n(\mathbf{r})^2}]$ .


title: Relativistic Dynamical Collapse Model for a Scalar Field

abstract: A natural generalization of the CSL (Continuous Spontaneous Localization) theory of dynamical collapse is applied to a relativistic quantum scalar field $\phi({\bf x},t)$. It is shown that the modified Schr\"odinger equation is relativistically invariant, that the probabilities associated to all possible values of the classical scalar random field $w({\bf x},t)$ (which determines the eventual state of collapse) add up to 1, that there is no energy production out of the vacuum and, in the limit of large time, the collapse is toward eigenstates of $\phi({\bf x},0)$.