title: Latin Polytopes

abstract: Latin squares are well studied combinatorial objects. In this paper we generalize the concept and propose new objects like Latin triangles, free Latin squares, Latin tetrahedra, free Latin cubes, etc. We start with a classic definition of Latin squares followed by one based on the concept of latinized board. A Latin square appears then as a combinatorial design whose points are geometric and whose lines (the rows and columns) are invariant under the symmetries of the square. The generalization that follows proceeds by 1. broadening this geometric symmetry 2. considering more general configurations of points and 3. admitting lines that intersect more freely. The resulting concept is the Latin board. Finally, we particularize Latin boards to define Latin polytopes, Latin polygons and Latin polyhedra.


title: On linearly related sequences of difference derivatives of discrete orthogonal polynomials

abstract: Let D_v the difference operator and q-difference operators defined by D_\omega p(x) = \frac{p(x+\omega)-p(x)}{\omega} and D_q p(x) = \frac{p(qx)-p(x)}{(q-1)x}, respectively. Let U and V be two moment regular linear functionals and let (P_n)_n and Q_n)_n be their corresponding orthogonal polynomial sequences (OPS). We discuss an inverse problem in the theory of discrete orthogonal polynomials involving the above two OPS assuming that their difference derivatives $D_\nu$ of higher orders m and k (resp.) are connected by a linear algebraic structure relation such as $$ \sum_{i=0}^M a_{i,n} D_\nu^m P_{n+m-i}(x) = \sum_{i=0}^N b_{i,n} D_\nu^k Q_{n+k-i}(x), \quad n\geq 0, $$ where $M,N,m,k=0,1,2,... Under certain conditions, we prove that U and V are related by a rational factor \c{c} (in the distributional sense). Moreover, when m\neq k then both U and V are D_v-semiclassical functionals. This leads us to the concept of (M,N)-D_v-coherent pair of order (m,k) extending to the discrete case several previous works. As an application we consider the OPS with respect to a certain following Sobolev-type discrete inner product.


title: The statistics of frictional families

abstract: We develop a theoretical description for mechanically stable frictional packings in terms of the difference between the total number of contacts required for isostatic packings of frictionless disks and the number of contacts in frictional packings, $m=N_c^0-N_c$. The saddle order $m$ represents the number of unconstrained degrees of freedom that a static packing would possess if friction were removed. Using a novel numerical method that allows us to enumerate disk packings for each $m$, we show that the probability to obtain a packing with saddle order $m$ at a given static friction coefficient $\mu$, $P_m(\mu)$, can be expressed as a power-series in $\mu$. Using this form for $P_m(\mu)$, we quantitatively describe the dependence of the average contact number on friction coefficient for static disk packings obtained from direct simulations of the Cundall-Strack model for all $\mu$ and $N$.


title: Noncommutative covering projections and $K$-homology

abstract: If $X$ is a topological space then there is a natural homomorphism $\pi_1(X)\rightarrow K_1(X)$ from a fundamental group to a $K_1$-homology group. Covering projections depend of fundamental group. So $K_1$-homology groups are interrelated with covering projections. This article is concerned with a noncommutative analogue of this interrelationship.


title: Permutation symmetry in spinor quantum gases: selection rules, conservation laws, and correlations

abstract: Many-body systems of identical arbitrary-spin particles, with separable spin and spatial degrees of freedom, are considered. Their eigenstates can be classified by Young diagrams, corresponding to non-trivial permutation symmetries (beyond the conventional paradigm of symmetric--antisymmetric states). The present work obtains (a) selection rules for additional non-separable (dependent on spins and coordinates) $k$-body interactions: the Young diagrams, associated with the initial and the final states of a transition, can differ by relocation of no more than $k$ boxes between their rows; and (b) correlation rules: eigenstate-averaged local correlations of $k$ particles vanish if $k$ exceeds the number of columns (for bosons) or rows (for fermions) in the associated Young diagram. It also elucidates the physical meaning of the quantities conserved due to permutation symmetry --- in 1929, Dirac identified those with characters of the symmetric group --- relating them to experimentally observable correlations of several particles. The results provide a way to control the formation of entangled states belonging to multidimensional non-Abelian representations of the symmetric group. These states can find applications in quantum computation and metrology.

categories: cond-mat.quant-gas math-ph math.MP physics.atom-ph quant-ph


title: Interpreting $Z_c(3900)$ and $Z_c(4025)/Z_c(4020)$ as charged tetraquark states

abstract: In the framework of color flux-tube model with a four-body confinement potential, the lowest charged tetraquark states $[Qq][\bar{Q}'\bar{q}']~(Q=c,b,q=u,d,s)$ are studied by using the variational method, Gaussian expansion method. The results indicate that some compact resonance states can be formed, the states can not decay into two color singlet mesons $Q\bar{q}'$ and $\bar{Q}'q$ through the breakdown and recombination of color flux tubes but into $Q\bar{Q}'$ and $q\bar{q}'$. The four-body confinement potential is an crucial dynamical mechanism for the formation of states, The decay mechanism is similar to that of compound nucleus and therefore the states should be called "color confined, multi-quark resonance" states. The newly observed charged states $Z_c(3900)$ and $Z_c(4025)/Z_c(4020)$ can be accommodated in the color flux-tube model and can be interpreted as the $S$-wave tetraquark states $[cu][{\bar{c}\bar{d}}]$ with quantum numbers $I=1$ and $J=1$ and 2, respectively.


title: Uniqueness and stability result for Cauchy's equation of motion for a certain class of hyperelastic materials

abstract: We consider Cauchy's equation of motion for hyperelastic materials. The solution of this nonlinear initial-boundary value problem is the vector field which discribes the displacement which a particle of this material perceives when exposed to stress and external forces. This equation is of greatest relevance when investigating the behaviour of elastic, anisotropic composites and for the detection of defects in such materials from boundary measurements. This is why results on unique solvability and continuous dependence from the initial values are of large interest in materials research and structural health monitoring. In this article we present such a result, provided that reasonable smoothness assumptions for the displacement field and the boundary of the domain are satisfied, for a certain class of hyperelastic materials where the first Piola-Kirchhoff tensor is written as a conic combination of finitely many given tensors.


title: UNLocBoX: A MATLAB convex optimization toolbox for proximal-splitting methods

abstract: Convex optimization is an essential tool for machine learning, as many of its problems can be formulated as minimization problems of specific objective functions. While there is a large variety of algorithms available to solve convex problems, we can argue that it becomes more and more important to focus on efficient, scalable methods that can deal with big data. When the objective function can be written as a sum of "simple" terms, proximal splitting methods are a good choice. UNLocBoX is a MATLAB library that implements many of these methods, designed to solve convex optimization problems of the form $\min_{x \in \mathbb{R}^N} \sum_{n=1}^K f_n(x).$ It contains the most recent solvers such as FISTA, Douglas-Rachford, SDMM as well a primal dual techniques such as Chambolle-Pock and forward-backward-forward. It also includes an extensive list of common proximal operators that can be combined, allowing for a quick implementation of a large variety of convex problems.


title: Model performance indicators ERP systems

abstract: Implementation process ERP is complex and expensive process. Typically always be faced with many failures. Successfully implemented in an organization has many challenges. Organizations in the deployment and success of the system depends on several factors.One of the key factors in the successful deployment of systems methodology is the implementation process. Methodology has several indicators for successful implementation of ERP systems, we have examined. And indicators for each of the methodologies have identified. The proposed method is also an important indicator of the success of security controls and indicators to be monitored and controlled.


title: Covering spaces of character varieties

abstract: Let F be a finitely generated discrete group. Given a covering map H to G of Lie groups with G either compact or complex reductive, there is an induced covering map Hom(F, H) to Hom(F, G). We show that when the fundamental group of G is torsion-free and F is free, free Abelian, or the fundamental group of a closed Riemann surface M of genus g, this map induces a covering map between the corresponding moduli spaces of representations. We give conditions under which this map is actually the universal covering, leading to new information regarding fundamental groups of these moduli spaces. Let pi be the fundamental group of M. As an application, we show that for g>0, the stable moduli space Hom(pi, SU)/SU is homotopy equivalent to infinite complex projective space. In the Appendix by Ho and Liu, it is shown show that there is a bijection between the number of connected components of Hom(pi, G) and the fundamental group of [G,G] for all complex connected reductive Lie groups G.


title: Floquet calculation of harmonic generation from hydrogen molecular ions in strong monochromatic laser fields

abstract: We present Floquet calculations of high harmonic generation (HHG) for the lowest two electronic states of the $\mbox{H}_2^+$ ion by strong continuous-wave laser fields. We solve the non-Hermitian matrix problem to get accurate solutions to the periodic time-dependent Schr\"odinger equation (TDSE) by applying a pseudospectral representation combined with a complex absorbing potential method. This represents an alternative approach to direct TDSE solutions to obtain the harmonic spectra for the ion. We compare our HHG rates for the lower and upper states of $\mbox{H}^{+}_{2}$, which correspond to the gerade and ungerade ground states in the field-free case, with previously obtained results in the literature. We show that the enhancement of the ionization rates at the critical internuclear separation $R_{c}\approx 8\,au$ plays some role in the appearance of very strong harmonic orders $n=5-11$ at $\lambda = 1064\,nm$ and $n=5-9$ at $\lambda = 800\,nm$ and intensity $I=10^{14}\,W/cm^{2}$.


title: Appearance of an accretion disk perturbed by fractional Brownian Motion density

abstract: This paper aims to investigate/map the effects that perturbations applied to an accretion disk might produce on the registered Light Curves (LC). The case of accretion disks around supermassive active black holes (AGNs) is studied with the goal to explain some of the statistical properties of the observed IntraDay Variability (IDV). The region producing optical IDV is perturbed by allowing it to develop a mass density of a fractional Brownian Motion-like type. The light curves and spectral slopes are calculated and compared to observational data for different Hurst parameters. The spectral slopes of the simulated light curves vary in the range $(0.4,2.5)$. The agreement with observational data shows that a magnetized disk subjected to stochastic perturbations can produce some of the features observed in the light curves.


title: Nonstandard functional interpretations and categorical models

abstract: Recently, the second author, Briseid and Safarik introduced nonstandard Dialectica, a functional interpretation that is capable of eliminating instances of familiar principles of nonstandard arithmetic - including overspill, underspill, and generalisations to higher types - from proofs. We show that, under few metatheoretical assumptions, the properties of this interpretation are mirrored by first order logic in a constructive sheaf model of nonstandard arithmetic due to Moerdijk, later developed by Palmgren. In doing so, we also draw some new connections between nonstandard principles, and principles that are rejected by strict constructivism. Furthermore, we introduce a variant of the Diller-Nahm interpretion with two different kinds of quantifiers (with and without computational meaning), similar to Hernest's light Dialectica interpretation, and show that one can obtain nonstandard Dialectica from this by weakening the computational content of the existential quantifiers -- a process we call herbrandisation. We also define a constructive sheaf model mirroring this new functional interpretation and show that the process of herbrandisation has a clear meaning in terms of these sheaf models.


title: Signal to Noise Ratio in Lensless Compressive Imaging

abstract: We analyze the signal to noise ratio (SNR) in a lensless compressive imaging (LCI) architecture. The architecture consists of a sensor of a single detecting element and an aperture assembly of an array of programmable elements. LCI can be used in conjunction with compressive sensing to capture images in a compressed form of compressive measurements. In this paper, we perform SNR analysis of the LCI and compare it with imaging with a pinhole or a lens. We will show that the SNR in the LCI is independent of the image resolution, while the SNR in either pinhole aperture imaging or lens aperture imaging decreases as the image resolution increases. Consequently, the SNR in the LCI is much higher if the image resolution is large enough.


title: Maslov indices, Poisson brackets, and singular differential forms

abstract: Maslov indices are integers that appear in semiclassical wave functions and quantization conditions. They are often notoriously difficult to compute. We present methods of computing the Maslov index that rely only on typically elementary Poisson brackets and simple linear algebra. We also present a singular differential form, whose integral along a curve gives the Maslov index of that curve. The form is closed but not exact, and transforms by an exact differential under canonical transformations. We illustrate the method with the $6j$-symbol, which is important in angular momentum theory and in quantum gravity.


title: micrOMEGAs : a code for the calculation of Dark Matter properties in generic models of particle interaction

abstract: These lecture notes describe the micrOMEGAs code for the calculation of Dark Matter observables in extensions of the standard model.


title: Earth and Moon orbital anomalies

abstract: A time-dependent gravitational constant or mass would correctly describe the suspected increasing of both: the Astronomical unit and the eccentricity of the Lunar orbit around the Earth.


title: Analytic Approach to the Approximate Solution of the Independent DGLAP Evolution Equations with Respect to the Hard-Pomeron Behavior

abstract: We show that it is possible to use hard-Pomeron behavior to the gluon distribution and singlet structure function at low $x$. We derive a second-order independent differential equation for the gluon distribution and the singlet structure function. In this approach, both singlet quarks and gluons have the same high-energy behavior at small $x$. These equations are derived from the next-to-leading order DGLAP evolution equations. All results can be consistently described in the framework of perturbative QCD, which shows an increase of gluon distribution and singlet structure functions as $x$ decreases.


title: Detecting Memory and Structure in Human Navigation Patterns Using Markov Chain Models of Varying Order

abstract: One of the most frequently used models for understanding human navigation on the Web is the Markov chain model, where Web pages are represented as states and hyperlinks as probabilities of navigating from one page to another. Predominantly, human navigation on the Web has been thought to satisfy the memoryless Markov property stating that the next page a user visits only depends on her current page and not on previously visited ones. This idea has found its way in numerous applications such as Google's PageRank algorithm and others. Recently, new studies suggested that human navigation may better be modeled using higher order Markov chain models, i.e., the next page depends on a longer history of past clicks. Yet, this finding is preliminary and does not account for the higher complexity of higher order Markov chain models which is why the memoryless model is still widely used. In this work we thoroughly present a diverse array of advanced inference methods for determining the appropriate Markov chain order. We highlight strengths and weaknesses of each method and apply them for investigating memory and structure of human navigation on the Web. Our experiments reveal that the complexity of higher order models grows faster than their utility, and thus we confirm that the memoryless model represents a quite practical model for human navigation on a page level. However, when we expand our analysis to a topical level, where we abstract away from specific page transitions to transitions between topics, we find that the memoryless assumption is violated and specific regularities can be observed. We report results from experiments with two types of navigational datasets (goal-oriented vs. free form) and observe interesting structural differences that make a strong argument for more contextual studies of human navigation in future work.


title: Operation of a planar-electrode ion-trap array with adjustable RF electrodes

abstract: One path to realizing systems of trapped atomic ions suitable for large-scale quantum computing and simulation is to create a two-dimensional array of ion traps. Interactions between nearest-neighbouring ions could then be turned on and off by tuning the ions' relative positions and frequencies. We demonstrate and characterize the operation of a planar-electrode ion-trap array. Driving the trap with a network of phase-locked radio-frequency (RF) resonators which provide independently variable voltage amplitudes we vary the position and motional frequency of a 40Ca+ ion in two dimensions within the trap array. With suitable miniaturization of the trap structure, this provides a viable architecture for large-scale quantum simulations.


title: Trace Formula For Two Variables

abstract: A natural generalization of Krein's theorem to a pair of commuting tuples $\left(H_1^0,H_2^0\right)$ and $\left(H_1,H_2\right)$ of bounded self-adjoint operators in a separable Hilbert space $\mathcal{H}$ with $H_j-H_j^0 = V_j\in \mathcal{B}_2(\mathcal{H})$(set of all Hilbert-Schmidt operators on $\mathcal{H}$) for $j=1,2,$ leads to a Stokes-like formula under trace. A major ingredient in the proof is the finite-dimensional approximation result for commuting self-adjoint n-tuples of operators, a generalization of Weyl-von Neumann-Berg's theorem.


title: Astrophysical weak-interaction rates for selected $A=20$ and $A=24$ nuclei

abstract: We have evaluated the electron capture rates on $^{20}$Ne, $^{20}$F, $^{24}$Mg, $^{24}$Na and the $\beta$ decay rates for $^{20}$F and $^{24}$Na at temperature and density conditions relevant for the late-evolution stages of stars with $M=8$-12 M$_\odot$. The rates are based on recent experimental data and large-scale shell model calculations. We show that the electron capture rates on $^{20}$Ne, $^{24}$Mg and the $^{20}$F, $^{24}$Na $\beta$-decay rates are based on data in this astrophysical range, except for the capture rate on $^{20}$Ne, which we predict to have a dominating contribution from the second-forbidden transition between the $^{20}$Ne and $^{20}$F ground states in the density range $\log \rho Y_e (\mathrm{g~cm}^{-3}) = 9.3$-9.6. The dominance of a few individual transitions allows us to present the various rates by analytical expressions at the relevant astrophysical conditions. We also derive the screening corrections to the rates.


title: A Game Theoretic Analysis of Collaboration in Wikipedia

abstract: Peer production projects such as Wikipedia or open-source software development allow volunteers to collectively create knowledge based products. The inclusive nature of such projects poses difficult challenges for ensuring trustworthiness and combating vandalism. Prior studies in the area deal with descriptive aspects of peer production, failing to capture the idea that while contributors collaborate, they also compete for status in the community and for imposing their views on the product. In this paper we investigate collaborative authoring in Wikipedia where contributors append and overwrite previous contributions to a page. We assume that a contributors goal is to maximize ownership of content sections such that content owned (or originated) by her survived the most recent revision of the page. We model contributors interactions to increase their content ownership as a noncooperative game where a players utility is associated with content owned and cost is a function of effort expended. Our results capture several real life aspects of contributors interactions within peer production projects. We show that at the Nash equilibrium there is an inverse relationship between the effort required to make a contribution and the survival of a contributors content. In other words majority of the content that survives is necessarily contributed by experts who expend relatively less effort than non experts. An empirical analysis of Wikipedia articles provides support for our models predictions. Implications for research and practice are discussed in the context of trustworthy collaboration as well as vandalism.


title: The Mathematics of Fels Sculptures

abstract: We give a purely mathematical interpretation and construction of sculptures rendered by one of the authors, known herein as Fels sculptures. We also show that the mathematical framework underlying Ferguson's sculpture, {\it The Ariadne Torus}, may be considered a special case of the more general constructions presented here. More general discussions are also presented about the creation of such sculptures whether they be virtual or in higher dimensional space.


title: Sequential Model-Based Ensemble Optimization

abstract: One of the most tedious tasks in the application of machine learning is model selection, i.e. hyperparameter selection. Fortunately, recent progress has been made in the automation of this process, through the use of sequential model-based optimization (SMBO) methods. This can be used to optimize a cross-validation performance of a learning algorithm over the value of its hyperparameters. However, it is well known that ensembles of learned models almost consistently outperform a single model, even if properly selected. In this paper, we thus propose an extension of SMBO methods that automatically constructs such ensembles. This method builds on a recently proposed ensemble construction paradigm known as agnostic Bayesian learning. In experiments on 22 regression and 39 classification data sets, we confirm the success of this proposed approach, which is able to outperform model selection with SMBO.


title: Extensions and traces of functions of bounded variation on metric spaces

abstract: In the setting of a metric space equipped with a doubling measure and supporting a Poincar\'e inequality, and based on results by Bj\"orn and Shanmugalingam (2007), we show that functions of bounded variation can be extended from any bounded uniform domain to the whole space. Closely related to extensions is the concept of boundary traces, which have previously been studied by Hakkarainen et al. (2014). On spaces that satisfy a suitable locality condition for sets of finite perimeter, we establish some basic results for the traces of functions of bounded variation. Our analysis of traces also produces novel results on the behavior of functions of bounded variation in their jump sets.


title: Isgur-Wise functions, Bjorken-Uraltsev Sum Rules and their Lorentz group interpretation, In Memoriam Nikolai Uraltsev

abstract: In the heavy quark limit of QCD, using the Operator Product Expansion and the non-forward amplitude, as proposed by Nikolai Uraltsev, we formulate sum rules that generalize Bjorken and Uraltsev sum rules. We recover the Uraltsev lower bound for the slope of the Isgur-Wise (IW) function, that we generalize to higher derivatives. We show that these results have a clear interpretation in terms of the Lorentz group, since the IW function is given by an overlap between the initial and final light clouds, related by Lorentz transformations. Both the Lorentz group and the Sum Rules approach are equivalent. Moreover, we formulate an integral representation of the IW function with a positive measure. Inverting this integral formula, we obtain the measure in terms of the IW function, allowing to formulate criteria to decide if a given ansatz for the IW function is compatible or not with the sum rule constraints. We compare these theoretical constraints to some forms proposed in the literature.


title: Transversals as generating sets in finitely generated groups

abstract: We explore transversals of finite index subgroups of finitely generated groups. We show that when $H$ is a subgroup of a rank $n$ group $G$ and $H$ has index at least $n$ in $G$ then we can construct a left transversal for $H$ which contains a generating set of size $n$ for $G$, and that the construction is algorithmic when $G$ is finitely presented. We also show that, in the case where $G$ has rank $n \leq3$, there is a simultaneous left-right transversal for $H$ which contains a generating set of size $n$ for $G$. We finish by showing that if $H$ is a subgroup of a rank $n$ group $G$ with index less than $3 \cdot 2^{n-1}$, and $H$ contains no primitive elements of $G$, then $H$ is normal in $G$ and $G/H \cong C_{2}^{n}$.


title: More on Asymptotic behaviors of 9j Coefficients

abstract: Previouslywe noted in numerical calculations that a certain unitary ninej coefficient U9j =((jj)^{2j} (jj)^{2j} | (jj)^{2j} (jj) ^{(2j-2)})^{I} decreases with increasing j and for small I the decrease is of the form C j^{m} e^{-\alpha j} . The exponential decay dominates. It is here shown that \alpha = 4 ln(2) and m=3/2.

categories: nucl-th physics.atm-clus physics.atom-ph quant-ph


title: Surgery along star shaped plumbings and exotic smooth structures on 4-manifolds

abstract: We define a new 4-dimensional symplectic cut and paste operation which is analogous to Fintushel and Stern's rational blow-down. We use this operation to produce multiple constructions of symplectic smoothly exotic complex projective space blown-up eight times, seven times, and six times. We also show how this operation can be used in conjunction with knot surgery to construct an infinite family of minimal exotic smooth structures on the complex projective space blown-up seven times.


title: Equations of Motion for Variational Electrodynamics

abstract: We extend the variational problem of Wheeler-Feynman electrodynamics by putting the electromagnetic functional in a local space of absolutely continuous trajectories possessing a derivative (velocities) of bounded variation. Generalizing the calculus of variations for extrema with a finite number of velocity discontinuities (breaking points), we prove that the critical-point-conditions for the two-body problem in the extended local space are Euler-Lagrange equations holding Lebesgue-almost-everywhere plus the generalized Weierstrass-Erdmann conditions that (i) the partial momenta must be absolutely continuous functions and (ii) the Legendre transforms of the partial Lagrangians (i.e, the partial energies) must be absolutely continuous functions.


title: Mode density and frequency extraction in the delta Scuti star HD 50844

abstract: We consider the high mode density reported in the delta Scuti star HD 50844 observed by CoRoT. Using simulations, we find that extracting frequencies down to a given false alarm probability by means of successive prewhitening leads to a gross over-estimate of the number of frequencies in a star. This is due to blending of the peaks in the periodogram due to the finite duration of the time series. Prewhitening is equivalent to adding a frequency to the data which is carefully chosen to interfere destructively with a given frequency in the data. Since the frequency extracted from a blended peak is not quite correct, the interference is not destructive with the result that many additional fictitious frequencies are added to the data. In data with very high signal-to-noise, such as the CoRoT data, these spurious frequencies are highly significant. Continuous prewhitening thus causes a cascade of spurious frequencies which leads to a much larger estimate of the mode density than is actually the case. The results reported for HD 50844 are consistent with this effect. Direct comparison of the power in the raw periodogram in this star with that in delta Scuti stars observed by Kepler shows that HD 50844 has a typical mode density.


title: A Measurement-based Analysis of the Energy Consumption of Data Center Servers

abstract: Energy consumption is a growing issue in data centers, impacting their economic viability and their public image. In this work we empirically characterize the power and energy consumed by different types of servers. In particular, in order to understand the behavior of their energy and power consumption, we perform measurements in different servers. In each of them, we exhaustively measure the power consumed by the CPU, the disk, and the network interface under different configurations, identifying the optimal operational levels. One interesting conclusion of our study is that the curve that defines the minimal CPU power as a function of the load is neither linear nor purely convex as has been previously assumed. Moreover, we find that the efficiency of the various server components can be maximized by tuning the CPU frequency and the number of active cores as a function of the system and network load, while the block size of I/O operations should be always maximized by applications. We also show how to estimate the energy consumed by an application as a function of some simple parameters, like the CPU load, and the disk and network activity. We validate the proposed approach by accurately estimating the energy of a map-reduce computation in a Hadoop platform.


title: On the vanishing of Hochster's theta invariant

abstract: Hochster's theta invariant is defined for a pair of finitely generated modules on a hypersurface ring having only an isolated singularity. Up to a sign, it agrees with the Euler invariant of a pair of matrix factorizations. Working over the complex numbers, Buchweitz and van Straten have established an interesting connection between Hochster's theta invariant and the classical linking form on the link of the singularity. In particular, they establish the vanishing of the theta invariant if the hypersurface is even-dimensional by exploiting the fact that the (reduced) cohomology of the Milnor fiber is concentrated in odd degrees in this situation. In this paper, we give purely algebraic versions of some of these results. In particular, we establish the vanishing of the theta invariant for isolated hypersurface singularities of odd dimension in characteristic p > 0, under some mild extra assumptions. This confirms, in a large number of cases, a conjecture of Hailong Dao.


title: Partition Function in One, Two and Three Spatial Dimensions from Effective Lagrangian Field Theory

abstract: The systematic effective Lagrangian method was first formulated in the context of the strong interaction: chiral perturbation theory (CHPT) is the effective theory of Quantum Chromodynamics (QCD). It was then pointed out that the method can be transferred to the nonrelativistic domain -- in particular, to describe the low-energy properties of ferromagnets. Interestingly, whereas for Lorentz-invariant systems the effective Lagrangian method fails in one spatial dimension ($d_s$=1), it perfectly works for nonrelativistic systems in $d_s$=1. In the present brief review, we give an outline of the method and then focus on the partition function for ferromagnetic spin chains, ferromagnetic films and ferromagnetic crystals up to three loops in the perturbative expansion -- an accuracy never achieved by conventional condensed matter methods. We then compare ferromagnets in $d_s$=1,2,3 with the behavior of QCD at low temperatures by considering the pressure and the order parameter. The two apparently very different systems (ferromagnets and QCD) are related from a universal point of view based on the spontaneously broken symmetry. In either case, the low-energy dynamics is described by an effective theory containing Goldstone bosons as basic degrees of freedom.


title: Efficient real frequency solver for dynamical mean field theory

abstract: We here present how a self-consistent solution of the dynamical mean field theory equations can be obtained using exact diagonalization of an Anderson impurity model with accuracies comparable to those found using renormalization group or quantum Monte Carlo methods. We show how one can solve a correlated quantum impurity coupled to several hundred uncorrelated bath sites, using a restricted active basis set. The number of bath sites determines the resolution of the obtained spectral function, which consists of peaks roughly spaced by the band width divided by the number of bath sites. The self-consistency cycles are fully performed on the real frequency axis and expressed as numerical stable matrix operations. The same impurity solver has been used on Ligand Field and finite size cluster calculations and is capable of treating involved Hamiltonians including the full rotational invariant Coulomb interaction, spin-orbit coupling and low symmetry crystal-fields. The proposed method allows for the calculation of a variety of correlation functions at little extra cost.


title: Associative Memories Based on Multiple-Valued Sparse Clustered Networks

abstract: Associative memories are structures that store data patterns and retrieve them given partial inputs. Sparse Clustered Networks (SCNs) are recently-introduced binary-weighted associative memories that significantly improve the storage and retrieval capabilities over the prior state-of-the art. However, deleting or updating the data patterns result in a significant increase in the data retrieval error probability. In this paper, we propose an algorithm to address this problem by incorporating multiple-valued weights for the interconnections used in the network. The proposed algorithm lowers the error rate by an order of magnitude for our sample network with 60% deleted contents. We then investigate the advantages of the proposed algorithm for hardware implementations.


title: Large Deviations for stationary probabilities of a family of continuous time Markov chains via Aubry-Mather theory

abstract: We consider a family of continuous time symmetric random walks indexed by $k\in \mathbb{N}$, $\{X_k(t),\,t\geq 0\}$. For each $k\in \mathbb{N}$ the matching random walk take values in the finite set of states $\Gamma_k=\frac{1}{k}(\mathbb{Z}/k\mathbb{Z})$ which is a subset of the unitary circle. The stationary probability for such process converges to the uniform distribution on the circle, when $k\to \infty$. We disturb the system considering a fixed $C^2$ potential $V: \mathbb{S}^1 \to \mathbb{R}$ and we will denote by $V_k$ the restriction of $V$ to $\Gamma_k$. Then, we define a non-stochastic semigroup generated by the matrix $k\,\, L_k + k\,\, V_k$, where $k\,\, L_k $ is the infinifesimal generator of $\{X_k(t),\,t\geq 0\}$. From the continuous time Perron's Theorem one can normalized such semigroup, and, then we get another stochastic semigroup which generates a continuous time Markov Chain taking values on $\Gamma_k$. The stationary probability vector for such Markov Chain is denoted by $\pi_{k,V}$. We assume that the maximum of $V$ is attained in a unique point $x_0$ of $\mathbb{S}^1$, and from this will follow that $\pi_{k,V}\to \delta_{x_0}$. Our main goal is to analyze the large deviation principle for the family $\pi_{k,V}$, when $k \to\infty$. The deviation function $I^V$, which is defined on $ \mathbb{S}^1$, will be obtained from a procedure based on fixed points of the Lax-Oleinik operator and Aubry-Mather theory.


title: Measures of disturbance and incompatibility for quantum measurements

abstract: We propose a class of incompatibility measures for quantum observables based on quantifying the effect of a measurement of one observable on the statistics of the outcomes of another. Specifically, for a pair of observables $A$ and $B$ with purely discrete spectra, we compare the following two probability distributions: one resulting from a measurement of $A$ followed by a measurement of $B$ on a given state, and the other obtained from a measurement of $B$ alone on the same state. We show that maximizing the distance between these two distributions over all states yields a valid measure of the incompatibility of observables $A$ and $B$, which is zero if and only if they commute and is strictly greater than zero (and less than or equal to one) otherwise. For finite dimensional systems, we obtain a tight upper bound on the incompatibility of any pair of observables and show that the bound is attained when the observables are totally non-degenerate and associated with mutually unbiased bases. In the process, we also establish an important connection between the incompatibility of a pair of observables and the maximal disturbances due to their measurements. Finally, we indicate how these measures of incompatibility and disturbance can be extended to the more general class of non-projective measurements. In particular, we obtain a non-trivial upper bound on the incompatibility of one L\"uders instrument with another.


title: New equidistribution estimates of Zhang type

abstract: We prove distribution estimates for primes in arithmetic progressions to large smooth squarefree moduli, with respect to congruence classes obeying Chinese Remainder Theorem conditions, obtaining an exponent of distribution $\frac{1}{2} + \frac{7}{300}$.


title: A Study on the Optimal Implementation of Statistical Multiplexing in DVB Distribution Systems

abstract: The paper presents an overview of the main methods used to improve the efficiency of DVB systems, based on multiplexing, through a study on the impact of the multiplexing methods used in DVB, having as a final goal a better usage of the data capacity and the possibility to insert new services into the original DVB Transport Stream. This study revealed that not all DVB providers are using statistical multiplexing. Based on this study, we were able to propose a method to improve the original DVB stream, originated from DVB-S or DVB-T providers. This method is proposing the detection of null packets, removal and reinserting a new service, with a VBR content. The method developed in this research can be implemented even in optimized statistical multiplexing systems, due to a residual use of null packets for data rate adjustment. There is no need to have access in the original stream multiplexer, since the method allows the implementation on the fly, near to the end user. The proposed method is proposed to be applied in DVB-S to DVB-C translation, using the computing power of a PC or in a FPGA implementation.


title: Observation of quantum interference in the plasmonic Hong-Ou-Mandel effect

abstract: We report direct evidence of the bosonic nature of surface plasmon polaritons (SPPs) in a scattering-based beamsplitter. A parametric down-conversion source is used to produce two indistinguishable photons, each of which is converted into a SPP on a metal-stripe waveguide and then made to interact through a semi-transparent Bragg mirror. In this plasmonic analog of the Hong-Ou-Mandel experiment, we measure a coincidence dip with a visibility of 72%, a key signature that SPPs are bosons and that quantum interference is clearly involved.


title: Accuracy tests of radiation schemes used in hot Jupiter global circulation models

abstract: The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.


title: Embeddability in the 3-sphere is decidable

abstract: We show that the following algorithmic problem is decidable: given a $2$-dimensional simplicial complex, can it be embedded (topologically, or equivalently, piecewise linearly) in $\mathbf{R}^3$? By a known reduction, it suffices to decide the embeddability of a given triangulated 3-manifold $X$ into the 3-sphere $S^3$. The main step, which allows us to simplify $X$ and recurse, is in proving that if $X$ can be embedded in $S^3$, then there is also an embedding in which $X$ has a short meridian, i.e., an essential curve in the boundary of $X$ bounding a disk in $S^3\setminus X$ with length bounded by a computable function of the number of tetrahedra of $X$.


title: Nonlinear photoacoustic wavefront shaping (PAWS) for single speckle-grain optical focusing in scattering media

abstract: Non-invasively focusing light into strongly scattering media, such as biological tissue, is highly desirable but challenging. Recently, wavefront shaping technologies guided by ultrasonic encoding or photoacoustic sensing have been developed to address this limitation. So far, these methods provide only acoustic diffraction-limited optical focusing. Here, we introduce nonlinear photoacoustic wavefront shaping (PAWS), which achieves optical diffraction-limited (i.e. single-speckle-grain) focusing in scattering media. We develop an efficient dual-pulse excitation approach to generate strong nonlinear photoacoustic (PA) signals based on the Grueneisen memory effect. These nonlinear PA signals are used as feedback to guide iterative wavefront optimization. By maximizing the amplitude of the nonlinear PA signal, light is effectively focused to a single optical speckle grain. Experimental results demonstrate a clear optical focus on the scale of 5-7 micrometers, which is ~10 times smaller than the acoustic focus in linear dimension, with an enhancement factor of ~6000 in peak fluence. This technology has the potential to provide highly confined strong optical focus deep in tissue for microsurgery of Parkinson's disease and epilepsy or single-neuron imaging and optogenetic activation.


title: A real space auxiliary field approach to the BCS-BEC crossover

abstract: The BCS to BEC crossover in attractive Fermi systems is a prototype of weak to strong coupling evolution in many body physics. While extensive numerical results are available, and several approximate methods have been developed, most of these schemes are unsuccessful in the presence of spatial inhomogeneity. Such situations call for a real space approach that can handle large spatial scales and retain the crucial thermal fluctuations. With this in mind, we present comprehensive results of a real space auxiliary field approach to the BCS to BEC crossover in the attractive Hubbard model in two dimensions. The scheme reproduces the Hartree-Fock-Bogoliubov ground state, and leads to a $T_c$ scale that agrees with quantum Monte Carlo estimates to within a few percent. We provide results on the $T_c$, amplitude and phase fluctuations, density of states, and the momentum resolved spectral function over the entire interaction and temperature window. We suggest how the method generalises successfully to the presence of disorder, trapping, and population imbalance.


title: Almost square Banach spaces

abstract: We single out and study a natural class of Banach spaces -- almost square Banach spaces. In an almost square space we can find, given a finite set $x_1,x_2,\ldots,x_N$ in the unit sphere, a unit vector $y$ such that $\|x_i-y\|$ is almost one. These spaces have duals that are octahedral and finite convex combinations of slices of the unit ball of an almost square space have diameter 2. We provide several examples and characterizations of almost square spaces. We prove that non-reflexive spaces which are M-ideals in their biduals are almost square. We show that every separable space containing a copy of $c_0$ can be renormed to be almost square. A local and a weak version of almost square spaces are also studied.


title: Observation of the Meissner effect with ultracold atoms in bosonic ladders

abstract: We report on the observation of the Meissner effect in bosonic flux ladders of ultracold atoms. Using artificial gauge fields induced by laser-assisted tunneling, we realize arrays of decoupled ladder systems that are exposed to a uniform magnetic field. By suddenly decoupling the ladders and projecting into isolated double wells, we are able to measure the currents on each side of the ladder. For large coupling strengths along the rungs of the ladder, we find a saturated maximum chiral current corresponding to a full screening of the artificial magnetic field. For lower coupling strengths, the chiral current decreases in good agreement with expectations of a vortex lattice phase. Our work marks the first realization of a low-dimensional Meissner effect and, furthermore, it opens the path to exploring interacting particles in low dimensions exposed to a uniform magnetic field.


title: Lifshitz Quartic Quasitopological Black Holes

abstract: In this paper we shall elucidate some of the effects of the quartic quasitopological term for Lifshitz-symmetric black holes. The field equations of this theory are difficult to solve exactly; here we will use numerical solutions both to verify previous exact solutions for quartic quasitopological AdS black holes as well as to examine new quasitopological Lifshitz-symmetric black hole solutions, in order to determine the effect of the quartic coupling parameter on the black hole's thermodynamic behavior. We shall find that the quartic parameter controls solutions very similarly to the cubic parameter, allowing for the construction of a theory with another free parameter which may find meaning in the phase transition behaviour of a gauge/gravity context.


title: Atomic form factor for twisted vortex photons interacting with atoms

abstract: The relatively new atomic form factor for twisted (vortex) beams, which carry orbital angular momentum (OAM), is considered and compared to the conventional atomic form factor for plane wave beams that carry only spin angular momentum (SAM). Since the vortex symmetry of a twisted photon is more complex that that of a plane-wave, evaluation of the atomic form factor is also more complex for twisted photons. On the other hand, the twisted photon has additional parameters, including the OAM quantum number, $\ell$, the nodal radial number, $p$, and the Rayleigh range, $z_R$ that determines the cone angle of the vortex. This Rayleigh range may be used as a variable parameter to control, in new ways, the interaction of twisted photons with matter. Here we address: i) normalization of the vortex atomic form factor, ii) displacement of target atoms away from the center of the beam vortex, and iii) formulation of transition probabilities for a variety of photon-atom processes. We attend to features related to new experiments that can test the range of validity and accuracy of calculations of these variations of the atomic form factor. Using the absolute square of the form factor for vortex beams, we introduce a vortex factor that can be directly measured.


title: Markov bridges: SDE representation

abstract: Let $X$ be a Markov process taking values in $\mathbf{E}$ with continuous paths and transition function $(P_{s,t})$. Given a measure $\mu$ on $(\mathbf{E}, \mathscr{E})$, a Markov bridge starting at $(s,\varepsilon_x)$ and ending at $(T^*,\mu)$ for $T^* <\infty$ has the law of the original process starting at $x$ at time $s$ and conditioned to have law $\mu$ at time $T^*$. We will consider two types of conditioning: a) {\em weak conditioning} when $\mu$ is absolutely continuous with respect to $P_{s,t}(x,\cdot)$ and b) {\em strong conditioning} when $\mu=\varepsilon_z$ for some $z \in \mathbf{E}$. The main result of this paper is the representation of a Markov bridge as a solution to a stochastic differential equation (SDE) driven by a Brownian motion in a diffusion setting. Under mild conditions on the transition density of the underlying diffusion process we establish the existence and uniqueness of weak and strong solutions of this SDE.


title: An interpretive conjecture for physics beyond the standard models: generalized complementarity

abstract: Our interpretive conjecture is inspired by the epistemology due to Ferdinand Gonseth (1890-1975) who interpreted complementarity as the relationship between profound and apparent reality horizons. It consists, on the one hand, on enlarging the scope of quantum theory to the most profound reality horizon, namely a triply quantum theory of gravitation that would be able to take into account simultaneously as elementary quanta the Planck's constant, the Planck's space-time area and the Boltzmann constant, and, on the other hand, on interpreting in terms of generalized complementarity three doubly quantum schemata taking into account by pairs, these three elementary quanta and form the apparent reality horizon.


title: New Stellar Models - Boon or Bane?

abstract: The impact of new stellar evolution models with rotation on the predictions of population synthesis models is discussed. Massive rotating stars have larger convective cores than their non-rotating counterparts, and their outer layers are chemically enriched due to increased mixing. Together, these two effects lead to hotter and more luminous stars, in particular during later evolutionary phases. As a result, stellar populations containing massive stars are predicted to become more luminous for a given mass and to emit more ionizing photons. Depending on the assumed rotation velocity, rotation causes profound changes in the properties of young stellar populations. These changes are most noticeable at later evolutionary phases and at shorter wavelengths of the spectral energy distribution. Most strikingly, the Lyman continuum luminosity increases by up to a factor of five in O- and Wolf-Rayet stars. Care is required when comparing these models to observations, and some fine-tuning of the models is still required before recalibrations of star-formation indicators should be attempted.


title: Generating function of the tilings of Aztec rectangle with holes

abstract: We consider a generating function of the domino tilings of an Aztec rectangle with several boundary unit squares removed. Our generating function involves two statistics: the rank of the tiling and half number of vertical dominoes as in the Aztec diamond theorem by Elkies, Kuperberg, Larsen and Propp. In addition, our work deduces a combinatorial explanation for an interesting connection between the number of lozenge tilings of a semihexagon and the number of domino tilings of an Aztec rectangle.


title: A Characterisation of Strong Integer Additive Set-Indexers of Graphs

abstract: An integer additive set-indexer is defined as an injective function $f:V(G)\rightarrow 2^{\mathbb{N}_0}$ such that the induced function $g_f:E(G) \rightarrow 2^{\mathbb{N}_0}$ defined by $g_f (uv) = f(u)+ f(v)$ is also injective, where $f(u)+f(v)$ is the sumset of $f(u)$ and $f(v)$. If $g_f(uv)=k~\forall~uv\in E(G)$, then $f$ is said to be a $k$-uniform integer additive set-indexers. An integer additive set-indexer $f$ is said to be a strong integer additive set-indexer if $|g_f(uv)|=|f(u)|.|f(v)|~\forall ~ uv\in E(G)$. We already have some characteristics of the graphs which admit strong integer additive set-indexers. In this paper, we study the characteristics of certain graph classes, graph operations and graph products that admit strong integer additive set-indexers.


title: Exploring the role of the charm quark in the $\Delta I=1/2$ rule

abstract: We study the dependence on the charm quark mass of the leading-order low-energy constants of the $\Delta S=1$ effective Hamiltonian, with the aim of elucidating the role of the charm mass scale in the $\Delta I=1/2$ rule for $K\to\pi\pi$ decay. To that purpose, finite-volume Chiral Perturbation Theory predictions are matched to QCD simulations, performed in the quenched approximation with overlap fermions and $m_u=m_d=m_s$. Light quark masses range between a few MeV up to around one third of the physical strange mass, while charm masses range between $m_u$ and a few hundred MeV. Novel variance reduction techniques are used to obtain a signal for penguin contractions in correlation functions involving four-fermion operators. The important role played by the subtractions required to construct renormalised amplitudes for $m_c \neq m_u$ is discussed in detail. We find evidence that the moderate enhancement of the $\Delta I=1/2$ amplitude previously found in the GIM limit $m_c=m_u$ increases only slightly as $m_c$ abandons the light quark regime. Hints of a stronger enhancement for even higher values of $m_c$ are also found, but their confirmation requires a better understanding of the subtraction terms.


title: Supersymmetry Breaking and the Cosmological Constant

abstract: I review three attempts to explain the small value of the cosmological constant, and their connection to SUSY breaking. They are The String Landscape, Supersymmetric Large Extra Dimensions (SLED), and the Holographic Space-time Formalism invented by Fischler and myself.


title: Sharp Estimates for Norms of Functions from Conjugate Classes in Metrics C and L

abstract: Sharp estimates for C - and L - norms of functions that are conjugate with functions from the classes W^rH_\omega of periodic functions having prescribed concave majorant of moduli of continuity, as well as sharp estimates for the differences of such functions are obtained.


title: A new perspective on least squares under convex constraint

abstract: Consider the problem of estimating the mean of a Gaussian random vector when the mean vector is assumed to be in a given convex set. The most natural solution is to take the Euclidean projection of the data vector on to this convex set; in other words, performing "least squares under a convex constraint." Many problems in modern statistics and statistical signal processing theory are special cases of this general situation. Examples include the lasso and other high-dimensional regression techniques, function estimation problems, matrix estimation and completion, shape-restricted regression, constrained denoising, linear inverse problems, etc. This paper presents three general results about this problem, namely, (a) an exact computation of the main term in the estimation error by relating it to expected maxima of Gaussian processes (existing results only give upper bounds), (b) a theorem showing that the least squares estimator is always admissible up to a universal constant in any problem of the above kind and (c) a counterexample showing that least squares estimator may not always be minimax rate-optimal. The result from part (a) is then used to compute the error of the least squares estimator in two examples of contemporary interest.


title: Variance reduction with practical all-to-all lattice propagators

abstract: We discuss all-to-all quark propagator techniques in two (related) contexts within Lattice QCD: the computation of closed quark propagators, and applications to the so-called "eye diagrams" appearing in the computation of non-leptonic kaon decay amplitudes. Combinations of low-mode averaging and diluted stochastic volume sources that yield optimal signal-to-noise ratios for the latter problem are developed. We also apply a recently proposed probing algorithm to compute directly the diagonal of the inverse Dirac operator, and compare its performance with that of stochastic methods. At fixed computational cost the two procedures yield comparable signal-to-noise ratios, but probing has practical advantages which make it a promising tool for a wide range of applications in Lattice QCD.


title: Large-Scale Asymmetries in the Transitional Disks of SAO 206462 and SR 21

abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations in the dust continuum (690 GHz, 0.45 mm) and 12CO J=6-5 spectral line emission, of the transitional disks surrounding the stars SAO 206462 and SR 21. These ALMA observations resolve the dust-depleted disk cavities and extended gaseous disks, revealing large-scale asymmetries in the dust emission of both disks. We modeled these disks structures with a ring and an azimuthal gaussian, where the azimuthal gaussian is motivated by the steady-state vortex solution from Lyra & Lin (2013). Compared to recent observations of HD 142527, Oph IRS 48, and LkHa 330, these are low-contrast (< 2) asymmetries. Nevertheless, a ring alone is not a good fit, and the addition of a vortex prescription describes these data much better. The asymmetric component encompasses 15% and 28% of the total disk emission in SAO 206462 and SR 21 respectively, which corresponds to a lower limit of 2 MJup of material within the asymmetry for both disks. Although the contrast in the dust asymmetry is low, we find that the turbulent velocity inside it must be large (~20% of the sound speed) in order to drive these azimuthally wide and radially narrow vortex-like structures. We obtain residuals from the ring and vortex fitting that are still significant, tracing non-axisymmetric emission in both disks. We compared these submillimeter observations with recently published H-band scattered light observations. For SR 21 the scattered light emission is distributed quite differently from the submillimeter continuum emission, while for SAO 206462 the submillimeter residuals are suggestive of spiral-like structure similar to the near-IR emission.


title: Evidence of resonant mode coupling and the relationship between low and high frequencies in a rapidly rotating A star

abstract: In the theory of resonant mode coupling, the parent and child modes are directly related in frequency and phase. The oscillations present in the fast rotating Delta Scuti star KIC 8054146 allow us to test the most general and generic aspects of such a theory. The only direct way to separate the parent and coupled (child) modes is to examine the correlations in amplitude variability between the different frequencies. For the dominant family of related frequencies, only a single mode and a triplet are the origins of nine dominant frequency peaks ranging from 2.93 to 66.30 cycles per day (as well as dozens of small-amplitude combination modes and a predicted and detected third high-frequency triplet). The mode-coupling model correctly predicts the large amplitude variations of the coupled modes as a product of the amplitudes of the parent modes, while the phase changes are also correctly modeled. This differs from the behavior of 'normal' combination frequencies in that the amplitudes are three orders of magnitude larger and may exceed even the amplitudes of the parent modes. We show that two dominant low frequencies at 5.86 and 2.93 cycles per day in the gravity-mode region are not harmonics of each other, and their properties follow those of the almost equidistant high-frequency triplet. We note that the previously puzzling situation of finding two strong peaks in the low-frequency region related by nearly a factor of two in frequency has been seen in other Delta Scuti stars as well.


title: Essential enhancements revisited

abstract: In 1991 Aizenman and Grimmett claimed that any `essential enhancement' of site or bond percolation on a lattice lowers the critical probability, an important result with many implications, such as strict inequalities between critical probabilities on suitable pairs of lattices. Their proof has two parts, one probabilistic and one combinatorial. In this paper we point out that a key combinatorial lemma, for which they provide only a figure as proof, is false. We prove an alternative form of the lemma, and thus the enhancement result, in the special cases of site percolation on the square, triangular and cubic lattices, and for bond percolation on $Z^d$, $d\ge 2$. The general case remains open, even for site percolation on $Z^d$, $d\ge 4$.


title: Results on identified particle production in pp, p-Pb and Pb-Pb collisions measured with ALICE at the LHC

abstract: Using the unique capabilities of the ALICE detectors for particle identification, different measurements have been performed to study the properties of the hot and dense matter created in the Pb-Pb collisions at sqrt(sNN)= 2.76 TeV. The analysis of the p-Pb data at sqrt(sNN) = 5.02 TeV reveals that the suppression of high pt hadrons observed in heavy nuclei collisions can not be explained as due to initial state effects. The systems created in the p-Pb collisions do not show evidence of jet quenching but, surprisingly, exhibit characteristics of flow. In this paper a review of the main results on identified particle production measured in different colliding systems is presented, data are also compared to models.


title: Cognitive Aging as Interplay between Hebbian Learning and Criticality

abstract: Cognitive ageing seems to be a story of global degradation. As one ages there are a number of physical, chemical and biological changes that take place. Therefore it is logical to assume that the brain is no exception to this phenomenon. The principle purpose of this project is to use models of neural dynamics and learning based on the underlying principle of self-organised criticality, to account for the age related cognitive effects. In this regard learning in neural networks can serve as a model for the acquisition of skills and knowledge in early development stages i.e. the ageing process and criticality in the network serves as the optimum state of cognitive abilities. Possible candidate mechanisms for ageing in a neural network are loss of connectivity and neurons, increase in the level of noise, reduction in white matter or more interestingly longer learning history and the competition among several optimization objectives. In this paper we are primarily interested in the affect of the longer learning history on memory and thus the optimality in the brain. Hence it is hypothesized that prolonged learning in the form of associative memory patterns can destroy the state of criticality in the network. We base our model on Tsodyks and Markrams [49] model of dynamic synapses, in the process to explore the effect of combining standard Hebbian learning with the phenomenon of Self-organised criticality. The project mainly consists of evaluations and simulations of networks of integrate and fire-neurons that have been subjected to various combinations of neural-level ageing effects, with the aim of establishing the primary hypothesis and understanding the decline of cognitive abilities due to ageing, using one of its important characteristics, a longer learning history.


title: On the existence of the $\it NN$-decoupled dibaryon $d^\ast_1$(1956)

abstract: We present strong experimental evidence for the existence of a nonstrange $\it NN$-decoupled dibaryon with a mass of about 1956 MeV stable against strong decay, which is a very likely candidate for long-sought multiquark hadrons with the baryon number B=2. We start by presenting the first evidence for the existence of this dibaryon called $d^{\ast }_{1}$(1956), that was found in the energy spectrum of coincident photons emitted at $\pm 90^0$ from the reaction $pp \to pp \gamma\gamma$ at an energy of 216 MeV measured at JINR. We then show its signatures found in several experimental energy spectra of single photons and invariant mass spectra of photon pairs from photon production processes induced by nucleon-nucleon and nucleon-nucleus collisions at intermediate energies.


title: The nature of the Roberge-Weiss transition in $N_f=2$ QCD with Wilson fermions

abstract: At imaginary values of the quark chemical potential $\mu$, Quantum Chromodynamics shows an interesting phase structure due to an exact center, or Roberge-Weiss (RW), symmetry. This can be used to constrain QCD at real $\mu$, where the sign problem prevents Monte Carlo simulations of the lattice theory. In previous studies of this region with staggered fermions it was found that the RW endpoint, where the center transition changes from first-order to a crossover, depends non-trivially on the quark mass: for high and low masses, it is a triple point connecting to the deconfinement and chiral transitions, respectively, changing to a second-order endpoint for intermediate mass values. These parameter regions are separated by tricritical points. Here we present a confirmation of these findings using Wilson fermions on $N_\tau=4$ lattices. In addition, our results provide a successful quantitative check for a heavy quark effective lattice theory at finite density.


title: Varieties of Picard rank one as components of ample divisors

abstract: Let $\mathcal{V}$ be an integral normal complex projective variety of dimension $n\geq 3$ and denote by $\mathcal{L}$ an ample line bundle on $\mathcal{V}$. By imposing that the linear system $|\mathcal{L}|$ contains an element $A=A_{1}+...+A_{r}, r\geq 1$, where all the $A_{i}$'s are distinct effective Cartier divisors with Pic$(A_i)=\mathbb{Z}$, we show that such a $\mathcal{V}$ is as special as the components $A_i$ of $A\in |\mathcal{L}|$. After making a list of some consequences about the positivity of the components $A_i$, we characterize pairs $(\mathcal{V}, \mathcal{L})$ as above when either $A_1\cong\mathbb{P}^{n-1}$ and Pic$(A_j)=\mathbb{Z}$ for $j=2,...,r,$ or $\mathcal{V}$ is smooth and each $A_i$ is a variety of small degree with respect to $[H_i]_{A_i}$, where $[H_i]_{A_i}$ is the restriction to $A_i$ of a suitable line bundle $H_i$ on $\mathcal{V}$.


title: Exchange-induced crystallization of soft core bosons

abstract: We study the phase diagram of a two-dimensional assembly of bosons interacting via a soft core repulsive pair potential of varying strength, and compare it to that of the equivalent system in which particles are regarded as distinguishable. We show that quantum-mechanical exchanges stabilize a "cluster crystal" phase in a wider region of parameter space than predicted by calculations in which exchanges are neglected. This physical effect is diametrically opposite to that which takes place in hard core Bose systems such as $^4$He, wherein exchanges strengthen the fluid phase. It is underlain in the cluster crystal phase of soft core bosons by the free energy gain associated to the formation of local Bose-Einstein condensates.


title: The Three-body problem and the shape sphere

abstract: [This is an expository article. I have submitted it to the American Mathematical Monthly.] The three-body problem defines a dynamics on the space of triangles in the plane. The shape sphere is the moduli space of oriented similarity classes of planar triangles and lies inside shape space, a Euclidean 3-space parametrizing oriented congruence classes of triangles. We derive and investigate the geometry and dynamics induced on these spaces by the three-body problem. We present two theorems concerning the three-body problem whose discovery was made through the shape space perspective


title: Quantum Control of Molecular Gas Hydrodynamics

abstract: We demonstrate that strong impulsive gas heating or heating suppression at standard temperature and pressure can occur from coherent rotational excitation or de-excitation of molecular gases using a sequence of non-ionizing laser pulses. For the case of excitation, subsequent collisional decoherence of the ensemble leads to gas heating significantly exceeding that from plasma absorption under the same laser focusing conditions. In both cases, the macroscopic hydrodynamics of the gas can be finely controlled with ~40 fs temporal sensitivity.


title: Monotone quantities involving a weighted $\sigma_k$ integral along inverse curvature flows

abstract: We give a family of monotone quantities along smooth solutions to the inverse curvature flows in Euclidean spaces. We also derive a related geometric inequality for closed hypersurfaces with positive k-th mean curvature.


title: On the theoretic and practical merits of the banding estimator for large covariance matrices

abstract: This paper considers the banding estimator proposed in Bickel and Levina (2008) for estimation of large covariance matrices. We prove that the banding estimator achieves rate-optimality under the operator norm, for a class of approximately banded covariance matrices, improving the existing results in Bickel and Levina (2008). In addition, we propose a Stein's unbiased risk estimate (Sure)-type approach for selecting the bandwidth for the banding estimator. Simulations indicate that the Sure-tuned banding estimator outperforms competing estimators.


title: The sign of the logistic regression coefficient

abstract: Let Y be a binary random variable and X a scalar. Let $\hat\beta$ be the maximum likelihood estimate of the slope in a logistic regression of Y on X with intercept. Further let $\bar x_0$ and $\bar x_1$ be the average of sample x values for cases with y=0 and y=1, respectively. Then under a condition that rules out separable predictors, we show that sign($\hat\beta$) = sign($\bar x_1-\bar x_0$). More generally, if $x_i$ are vector valued then we show that $\hat\beta=0$ if and only if $\bar x_1=\bar x_0$. This holds for logistic regression and also for more general binary regressions with inverse link functions satisfying a log-concavity condition. Finally, when $\bar x_1\ne \bar x_0$ then the angle between $\hat\beta$ and $\bar x_1-\bar x_0$ is less than ninety degrees in binary regressions satisfying the log-concavity condition and the separation condition, when the design matrix has full rank.


title: Near-IR Direct Detection of Water Vapor in Tau Boo b

abstract: We use high dynamic range, high-resolution L-band spectroscopy to measure the radial velocity variations of the hot Jupiter in the tau Bootis planetary system. The detection of an exoplanet by the shift in the stellar spectrum alone provides a measure of the planet's minimum mass, with the true mass degenerate with the unknown orbital inclination. Treating the tau Boo system as a high flux ratio double-lined spectroscopic binary permits the direct measurement of the planet's true mass as well as its atmospheric properties. After removing telluric absorption and cross-correlating with a model planetary spectrum dominated by water opacity, we measure a 6-sigma detection of the planet at K_p = 111 +- 5 km/s, with a 1-sigma upper limit on the spectroscopic flux ratio of 10^-4. This radial velocity leads to a planetary orbital inclination of i = 45+3-4degrees and a mass of M_P = 5.90+0.35-0.20 M_ Jup. We report the first detection of water vapor in the atmosphere of a non-transiting hot Jupiter, tau Boo b.


title: Generalized Moment Method for Gap Estimation and Quantum Monte Carlo Level Spectroscopy

abstract: We formulate a convergent sequence for the energy gap estimation in the worldline quantum Monte Carlo method. The ambiguity left in the conventional gap calculation for quantum systems is eliminated. Our estimation will be unbiased in the low-temperature limit and also the error bar is reliably estimated. The level spectroscopy from quantum Monte Carlo data is developed as an application of the unbiased gap estimation. From the spectral analysis, we precisely determine the Kosterlitz-Thouless quantum phase-transition point of the spin-Peierls model. It is stablished that the quantum phonon with a finite frequency is essential to the critical theory governed by the antiadiabatic limit, i.e., the $k=1$ SU(2) Wess-Zumino-Witten model.


title: A simple mechanism for the anti-glitch observed in magnetar 1E 2259+586

abstract: In this short communication we present a simple internal mechanism that accounts for the recently observed anti-glitch in magnetar 1E 2259+586. We propose that the decay of an internal toroidal magnetic field component would de-stabilize an originally stable prolate star configuration. Then, the subsequent rearrangement of the stellar structure would give rise to a "more spherical" configuration, resulting in a sudden spin-down of the whole star. We present here some order of magnitude calculation to give confidence to this scenario, using the simplest analytical stellar model and let more detailed calculations for a more technical future paper.


title: The rise of SN2014J in the nearby galaxy M82

abstract: We report on the discovery of SN2014J in the nearby galaxy M82. Given its proximity, it offers the best opportunity to date to study a thermonuclear supernova over a wide range of the electromagnetic spectrum. The first set of optical, near-IR and mid-IR observations of SN2014J, orchestrated by the intermediate Palomar Transient Factory (iPTF), show that SN2014J is a spectroscopically normal Type Ia supernova, albeit exhibiting high-velocity features in its spectrum and heavily reddened by dust in the host galaxy. Our earliest detections start just hours after the fitted time of explosion. We use high-resolution optical spectroscopy to analyze the dense intervening material and do not detect any evolution in the resolved absorption features during the lightcurve rise. Similarly to other highly reddened Type Ia supernovae, a low value of total-to-selective extinction, Rv < 2, provides the best match to our observations. We also study pre-explosion optical and near-IR images from HST with special emphasis on the sources nearest to the SN location.


title: RADIA: RNA and DNA Integrated Analysis for Somatic Mutation Detection

abstract: The detection of somatic single nucleotide variants is a crucial component to the characterization of the cancer genome. Mutation calling algorithms thus far have focused on comparing the normal and tumor genomes from the same individual. In recent years, it has become routine for projects like The Cancer Genome Atlas (TCGA) to also sequence the tumor RNA. Here we present RADIA (RNA and DNA Integrated Analysis), a method that combines the patient-matched normal and tumor DNA with the tumor RNA to detect somatic mutations. The inclusion of the RNA increases the power to detect somatic mutations, especially at low DNA allelic frequencies. By integrating the DNA and RNA, we are able to rescue back calls that would be missed by traditional mutation calling algorithms that only examine the DNA. RADIA was developed for the identification of somatic mutations using both DNA and RNA from the same individual. We demonstrate high sensitivity (84%) and very high specificity (98% and 99%) in real data from endometrial carcinoma and lung adenocarcinoma from TCGA. Mutations with both high DNA and RNA read support have the highest validation rate of over 99%. We also introduce a simulation package that spikes in artificial mutations to real data, rather than simulating sequencing data from a reference genome. We evaluate sensitivity on the simulation data and demonstrate our ability to rescue back calls at low DNA allelic frequencies by including the RNA. Finally, we highlight mutations in important cancer genes that were rescued back due to the incorporation of the RNA. Software available at https://github.com/aradenbaugh/radia/


title: Interval scheduling and colorful independent sets

abstract: Numerous applications in scheduling, such as resource allocation or steel manufacturing, can be modeled using the NP-hard Independent Set problem (given an undirected graph and an integer k, find a set of at least k pairwise non-adjacent vertices). Here, one encounters special graph classes like 2-union graphs (edge-wise unions of two interval graphs) and strip graphs (edge-wise unions of an interval graph and a cluster graph), on which Independent Set remains NP-hard but admits constant-ratio approximations in polynomial time. We study the parameterized complexity of Independent Set on 2-union graphs and on subclasses like strip graphs. Our investigations significantly benefit from a new structural "compactness" parameter of interval graphs and novel problem formulations using vertex-colored interval graphs. Our main contributions are: 1. We show a complexity dichotomy: restricted to graph classes closed under induced subgraphs and disjoint unions, Independent Set is polynomial-time solvable if both input interval graphs are cluster graphs, and is NP-hard otherwise. 2. We chart the possibilities and limits of effective polynomial-time preprocessing (also known as kernelization). 3. We extend Halld\'orsson and Karlsson (2006)'s fixed-parameter algorithm for Independent Set on strip graphs parameterized by the structural parameter "maximum number of live jobs" to show that the problem (also known as Job Interval Selection) is fixed-parameter tractable with respect to the parameter k and generalize their algorithm from strip graphs to 2-union graphs. Preliminary experiments with random data indicate that Job Interval Selection with up to fifteen jobs and 5*10^5 intervals can be solved optimally in less than five minutes.


title: Monodisperse cluster crystals: classical and quantum dynamics

abstract: We study the phases and dynamics of a gas of monodisperse particles interacting via soft-core potentials in two spatial dimensions, which is of interest for soft-matter colloidal systems and quantum atomic gases. Using exact theoretical methods, we demonstrate that the equilibrium low-temperature classical phase simultaneously breaks continuous translational symmetry and dynamic space-time homogeneity, whose absence is usually associated with out-of-equilibrium glassy phenomena. This results in an exotic self-assembled cluster crystal with coexisting liquid-like long-time dynamical properties, which corresponds to a classical analog of supersolid behavior. We demonstrate that the effects of quantum fluctuations and bosonic statistics on cluster-glassy crystals are separate and competing: zero-point motion tends to destabilize crystalline order, which can be restored by bosonic statistics.


title: First Dynamic Computation of Synchrotron Emission from the Cygnus A Radio Cavity; Evidence for Electron Pair Plasma in Cavity

abstract: Cosmic rays, thermal gas and magnetic fields in FRII radio cavities are assumed to come entirely from winds flowing from just behind the jet shocks. Combining analytic and computational methods, it is shown that the computed radio-electron energy distribution and synchrotron emissivity spectra everywhere in the Cygnus A radio cavity agrees with radio observations of the Cygnus A lobes. The magnetic field energy density is small everywhere and evolves passively in the post-shock wind. Most synchrotron emission arises in recent post-shock material as it flows back along the radio cavity wall. Because it experienced less adiabatic expansion, the magnetic field in this young backflow is larger than elsewhere in the radio lobe, explaining the observed radio synchrotron limb-brightening. The boundary backflow decelerates due to small cavity pressure gradients, causing large-scale fields perpendicular to the backflow (and synchrotron emission) to grow exponentially unlike observations. However, if the field is random on subgrid (sub-kpc) scales, the computed field reproduces both the magnitude and slowly decreasing radio synchrotron emissivity observed along the backflow. The radio synchrotron spectrum and image computed with a small-scale random field agree with VLA observations. The total relativistic energy density in the post-jet shock region required in computations to inflate the radio cavity matches the energy density of relativistic electrons observed in the post-shock region of Cygnus A. This indicates that the component in the jet and cavity that dominates the dynamical evolution is a relativistic pair plasma.


title: Advanced spectral analysis of ionospheric waves observed with sparse arrays

abstract: This paper presents a case study from a single, six-hour observing period to illustrate the application of techniques developed for interferometric radio telescopes to the spectral analysis of observations of ionospheric fluctuations with sparse arrays. We have adapted the deconvolution methods used for making high dynamic range images of cosmic sources with radio arrays to making comparably high dynamic range maps of spectral power of wavelike ionospheric phenomena. In the example presented here, we have used observations of the total electron content (TEC) gradient derived from Very Large Array (VLA) observations of synchrotron emission from two galaxy clusters at 330 MHz as well as GPS-based TEC measurements from a sparse array of 33 receivers located within New Mexico near the VLA. We show that these techniques provide a significant improvement in signal to noise (S/N) of detected wavelike structures by correcting for both measurement inaccuracies and wavefront distortions. This is especially true for the GPS data when combining all available satellite/receiver pairs, which probe a larger physical area and likely have a wider variety of measurement errors than in the single-satellite case. In this instance, we found the peak S/N of the detected waves was improved by more than an order of magnitude. The data products generated by the deconvolution procedure also allow for a reconstruction of the fluctuations as a two-dimensional waveform/phase screen that can be used to correct for their effects.


title: AMS-02 data confronts acceleration of cosmic ray secondaries in nearby sources

abstract: We revisit the model proposed earlier to account for the observed increase in the positron fraction in cosmic rays with increasing energy, in the light of new data from the Alpha Magnetic Spectrometer (AMS-02) experiment. The model accounts for the production and acceleration of secondary electrons and positrons in nearby supernova remnants which results in an additional, harder component that becomes dominant at high energies. By fitting this to AMS-02 data we can calculate the expected concomitant rise of the boron-to-carbon ratio, as well as of the fraction of antiprotons. If these predictions are confirmed by the forthcoming AMS-02 data it would conclusively rule out all other proposed explanations, in particular dark matter annihilations or decays.


title: Network Traffic Anomaly Detection

abstract: This paper presents a tutorial for network anomaly detection, focusing on non-signature-based approaches. Network traffic anomalies are unusual and significant changes in the traffic of a network. Networks play an important role in today's social and economic infrastructures. The security of the network becomes crucial, and network traffic anomaly detection constitutes an important part of network security. In this paper, we present three major approaches to non-signature-based network detection: PCA-based, sketch-based, and signal-analysis-based. In addition, we introduce a framework that subsumes the three approaches and a scheme for network anomaly extraction. We believe network anomaly detection will become more important in the future because of the increasing importance of network security.


title: Fully Constrained Majorana Neutrino Mass Matrices using $\Sigma(72\times3)$

abstract: In 2002, two neutrino mixing ansatze having trimaximally mixed middle ($\nu_2$) columns, namely tri-chi-maximal mixing ($\text{T}\chi\text{M}$) and tri-phi-maximal mixing ($\text{T}\phi\text{M}$), were proposed. It was recently shown that $\text{T}\chi\text{M}$ with $\chi=\pm \frac{\pi}{16}$ as well as $\text{T}\phi\text{M}$ with $\phi = \pm \frac{\pi}{16}$ leads to the solution, $\sin^2 \theta_{13} = \frac{2}{3} \sin^2 \frac{\pi}{16}$, consistent with the latest measurements of the reactor mixing angle, $\theta_{13}$. To obtain $\text{T}\chi\text{M}_{(\chi=\pm \frac{\pi}{16})}$ and $\text{T}\phi\text{M}_{(\phi=\pm \frac{\pi}{16})}$, we utilised the type I see-saw framework with fully constrained Majorana neutrino mass matrices. These mass matrices also resulted in a relation among the neutrino masses, $m_1:m_2:m_3=\frac{\left(2+\sqrt{2}\right)}{1+\sqrt{2(2+\sqrt{2})}}:1:\frac{\left(2+\sqrt{2}\right)}{-1+\sqrt{2(2+\sqrt{2})}}$. In this paper we construct a flavour model based on the discrete group $\Sigma(72\times3)$ and obtain the aforementioned results. A Majorana neutrino mass matrix (a symmetric $3\times3$ matrix with 6 complex degrees of freedom) is conveniently mapped into a flavon field transforming as the complex 6 dimensional representation of $\Sigma(72\times3)$. Specific vacuum alignments of the flavons are used to arrive at the desired mass matrices.


title: Robust Satisfiability of Systems of Equations

abstract: We study the problem of \emph{robust satisfiability} of systems of nonlinear equations, namely, whether for a given continuous function $f:\,K\to\mathbb{R}^n$ on a~finite simplicial complex $K$ and $\alpha>0$, it holds that each function $g:\,K\to\mathbb{R}^n$ such that $\|g-f\|_\infty \leq \alpha$, has a root in $K$. Via a reduction to the extension problem of maps into a sphere, we particularly show that this problem is decidable in polynomial time for every fixed $n$, assuming $\dim K \le 2n-3$. This is a substantial extension of previous computational applications of \emph{topological degree} and related concepts in numerical and interval analysis. Via a reverse reduction we prove that the problem is undecidable when $\dim K\ge 2n-2$, where the threshold comes from the \emph{stable range} in homotopy theory. For the lucidity of our exposition, we focus on the setting when $f$ is piecewise linear. Such functions can approximate general continuous functions, and thus we get approximation schemes and undecidability of the robust satisfiability in other possible settings.


title: The Informed Sampler: A Discriminative Approach to Bayesian Inference in Generative Computer Vision Models

abstract: Computer vision is hard because of a large variability in lighting, shape, and texture; in addition the image signal is non-additive due to occlusion. Generative models promised to account for this variability by accurately modelling the image formation process as a function of latent variables with prior beliefs. Bayesian posterior inference could then, in principle, explain the observation. While intuitively appealing, generative models for computer vision have largely failed to deliver on that promise due to the difficulty of posterior inference. As a result the community has favoured efficient discriminative approaches. We still believe in the usefulness of generative models in computer vision, but argue that we need to leverage existing discriminative or even heuristic computer vision methods. We implement this idea in a principled way with an "informed sampler" and in careful experiments demonstrate it on challenging generative models which contain renderer programs as their components. We concentrate on the problem of inverting an existing graphics rendering engine, an approach that can be understood as "Inverse Graphics". The informed sampler, using simple discriminative proposals based on existing computer vision technology, achieves significant improvements of inference.


title: Decomposition of random graphs into complete bipartite graphs

abstract: We consider the problem of partitioning the edge set of a graph $G$ into the minimum number $\tau(G)$ of edge-disjoint complete bipartite subgraphs. We show that for a random graph $G$ in $G(n,p)$, for $p$ is a constant no greater than $1/2$, almost surely $\tau(G)$ is between $n- c(\ln_{1/p} n)^{3+\epsilon}$ and $n - 2\ln_{1/(1-p)} n$ for any positive constants $c$ and $\epsilon$.


title: Voting for Committees in Agreeable Societies

abstract: We examine the following voting situation. A committee of $k$ people is to be formed from a pool of n candidates. The voters selecting the committee will submit a list of $j$ candidates that they would prefer to be on the committee. We assume that $j \leq k < n$. For a chosen committee, a given voter is said to be satisfied by that committee if her submitted list of $j$ candidates is a subset of that committee. We examine how popular is the most popular committee. In particular, we show there is always a committee that satisfies a certain fraction of the voters and examine what characteristics of the voter data will increase that fraction.


title: Fair Division and Redistricting

abstract: Recently, Landau, Reid and Yershov provided a novel solution to the problem of redistricting. Instead of trying to ensure fairness by restricting the shape of the possible maps or by assigning the power to draw the map to nonbiased entities, the solution ensures fairness by balancing competing interests against each other. This kind of solution is an example of what are known as "fair division" solutions--- such solutions involve the preferences of all parties and are accompanied by rigorous guarantees of a specified well-defined notion of fairness. In this expository article, we give an introduction to the ideas of fair division in the context of this redistricting solution. Through examples and discussion we clarify how fair division methods can play an important role in a realistic redistricting solution by introducing an interactive step that incorporates a certain kind of fairness that can be used in concert with, and not a substitute for, other necessary or desired criteria for a good redistricting solution.


title: A comparison of Euclidean metrics and their application in statistical inferences in the spike train space

abstract: Statistical analysis and inferences on spike trains are one of the central topics in neural coding. It is of great interest to understand the underlying distribution and geometric structure of given spike train data. However, a fundamental obstacle is that the space of all spike trains is not an Euclidean space, and non-Euclidean metrics have been commonly used in the literature to characterize the variability and pattern in neural observations. Over the past few years, two Euclidean-like metrics were independently developed to measure distance in the spike train space. An important benefit of these metrics is that the spike train space will be suitable for embedding in Euclidean spaces due to their Euclidean properties. In this paper, we systematically compare these two metrics on theory, properties, and applications. Because of its Euclidean properties, one of these metrics has been further used in defining summary statistics (i.e. mean and variance) and conducting statistical inferences in the spike train space. Here we provide equivalent definitions using the other metric and show that consistent statistical inferences can be conducted. We then apply both inference frameworks in a neural coding problem for a recording in geniculate ganglion stimulated by different tastes. It is found that both frameworks achieve desirable results and provide useful new tools in statistical inferences in neural spike train space.


title: Shadowing correction to the gluon distribution behavior at small x

abstract: We determined the saturation exponent of the gluon distribution using the solution of the QCD nonlinear Dokshitzer-Gribov-Lipatov-Altarelli-parisi (NLDGLAP) evolution equation at small $x$. The very small $x$ behavior of the gluon distribution is obtained by solving the Gribov, Levin, Ryskin, Mueller and Qiu (GLR-MQ) evolution equation with the nonlinear shadowing term incorporated. The form of initial condition for the equation is determined. We find, with decreasing $x$, the emergence of a singular behavior and the eventual taming (at $R=5\hspace{0.1cm}GeV^{-1}$) and the essential taming (at $R=2\hspace{0.1cm}GeV^{-1}$) of this singular behavior by shadowing term. The nonlinear gluon density functions are calculated and compared with the results for the integrated gluon density from the Balitsky- Kovchegov(BK) equation for the different values of $Q^{2}$. It is shown, that the results for the gluon density function are comparable with the results obtained from BK equation solution. Also we show that for each $x$, the $Q^{2}$ dependence of the data is well described by gluon shadowing corrections to GLR-MQ equation. The resulting analytic expression allow us to predict the logarithmic derivative $\frac{{\partial}F^{s}_{2}(x,Q^{2})}{{\partial}lnQ^{2}}$ and to compare the results with H1 data and a QCD analysis fit.


title: On Statistics of the Riemann Zeros Differences

abstract: Numerical study of the distribution of the Riemann zeros differences following the work [1] shows the significance of the function for which the prime sum expression is proposed. Computational results related to this definition explored with various prime cut-offs.


title: Bulge Growth and Quenching since z = 2.5 in CANDELS/3D-HST

abstract: Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.


title: The Big Problems in Star Formation: the Star Formation Rate, Stellar Clustering, and the Initial Mass Function

abstract: Star formation lies at the center of a web of processes that drive cosmic evolution: generation of radiant energy, synthesis of elements, formation of planets, and development of life. Decades of observations have yielded a variety of empirical rules about how it operates, but at present we have no comprehensive, quantitative theory. In this review I discuss the current state of the field of star formation, focusing on three central questions: what controls the rate at which gas in a galaxy converts to stars? What determines how those stars are clustered, and what fraction of the stellar population ends up in gravitationally-bound structures? What determines the stellar initial mass function, and does it vary with star-forming environment? I use these three question as a lens to introduce the basics of star formation, beginning with a review of the observational phenomenology and the basic physical processes. I then review the status of current theories that attempt to solve each of the three problems, pointing out links between them and opportunities for theoretical and numerical work that crosses the scale between them. I conclude with a discussion of prospects for theoretical progress in the coming years.