title: Re-emergence of rapidity scale uncertainty in SCET

abstract: The artificial separation of a full-theory mode into distinct collinear and soft modes in SCET leads to divergent integrals over rapidity, which are not present in the full theory. Rapidity divergence introduces an additional scale into the problem, giving rise to its own renormalization group with respect to this new scale. Two contradicting claims exist in the literature concerning rapidity scale uncertainty. One camp has shown that the results of perturbative calculations depend on the precise choice of rapidity scale. The other has derived an all-order factorization formula with no dependence on rapidity scale, by using a form of analytic regulator to regulate rapidity divergences. We deliver a simple resolution to this controversy by deriving an alternative form of the all-order factorization formula with an analytic regulator that, despite being formally rapidity scale independent, reveals how rapidity scale dependence arises when it is truncated at a finite order in perturbation theory. With our results, one can continue to take advantage of the technical ease and simplicity of the analytic regulator approach while correctly taking into account rapidity scale dependence. As an application, we update our earlier study of WW production with jet-veto by including rapidity scale uncertainty. While the central values of the predictions are unchanged, the scale uncertainties are increased and consistency between the NLL and NNLL calculations are improved.


title: Bounds on Higgs And Top Quark Masses From The Other Degenerate Vacua Near The Planck Scale With Gravitational Contributions

abstract: Based on the weak coupling expansion of gravity, we calculate the gravitational contributions to yukawa coupling, scalar quartic coupling as well as gauge couplings with general Landau-DeWitt gauge-fixing choice and a gauge preserving (of SM gauge group) cut off regularization scheme. We find that the results depend on the Landau-DeWitt gauge-fixing parameter. Based on the two loop RGE of SM couplings with one loop full gravitational contributions in harmonic gauge, we study the constraints on the higgs and top quark mass from the requirement of existing the other degenerate vacua at the Planck-dominated region. Our numerical calculations show that nature will not develop the other degenerate vacua at the Planck-dominated region with current higgs and top quark masses. On the other hand, requiring the existence of the other degenerate vacua at the Planck-dominated region will constrain the higgs and top mass to lie at approximately 130 and 174 GeV, respectively.


title: A ~50,000 solar mass black hole in the nucleus of RGG 118

abstract: Scaling relations between black hole (BH) masses and their host galaxy properties have been studied extensively over the last two decades, and point towards co-evolution of central massive BHs and their hosts. However, these relations remain poorly constrained for BH masses below $\sim10^{6}$ M_sun. Here we present optical and X-ray observations of the dwarf galaxy RGG 118 taken with the Magellan Echellette Spectrograph on the 6.5m Clay Telescope and Chandra X-ray Observatory. Based on Sloan Digital Sky Survey spectroscopy, RGG 118 was identified as possessing narrow emission line ratios indicative of photoionization partly due to an active galactic nucleus. Our higher resolution spectroscopy clearly reveals broad H$\alpha$ emission in the spectrum of RGG 118. Using virial BH mass estimate techniques, we calculate a BH mass of $\sim50,000$ \msun. We detect a nuclear X-ray point source in RGG 118, suggesting a total accretion powered luminosity of $L=4\times10^{40}~{\rm erg~s^{-1}}$, and an Eddington fraction of $\sim1$ per cent. The BH in RGG 118 is the smallest ever reported in a galaxy nucleus and we find that it lies on the extrapolation of the $M_{\rm BH}-\sigma_{\ast}$ relation to the lowest masses yet.


title: Freeze-In Dark Matter with Displaced Signatures at Colliders

abstract: Dark matter, $X$, may be generated by new physics at the TeV scale during an early matter-dominated (MD) era that ends at temperature $T_R \ll {\rm TeV}$. Compared to the conventional radiation-dominated (RD) results, yields from both Freeze-Out and Freeze-In processes are greatly suppressed by dilution from entropy production, making Freeze-Out less plausible while allowing successful Freeze-In with a much larger coupling strength. Freeze-In is typically dominated by the decay of a particle $B$ of the thermal bath, $B \rightarrow X$. For a large fraction of the relevant cosmological parameter space, the decay rate required to produce the observed dark matter abundance leads to displaced signals at LHC and future colliders, for any $m_X$ in the range ${\rm keV} < m_X < m_B$ and for values of $m_B$ accessible to these colliders. This result applies whether the early MD era arises after conventional inflation, when $T_R$ is the usual reheat temperature, or is a generic MD era with an alternative origin. In the former case, if $m_X$ is sufficiently large to be measured from kinematics, the reheat temperature $T_R$ can be extracted. Our result is independent of the particular particle physics implementation of $B \rightarrow X$, and can occur via any operator of dimension less than 8 (4) for a post-inflation (general MD) cosmology. An interesting example is provided by DFS axion theories with TeV-scale supersymmetry and axino dark matter of mass GeV to TeV, which is typically overproduced in a conventional RD cosmology. If $B$ is the higgsino, $\tilde h$, Higgs, W and Z particles appear at the displaced decays, $\tilde h \rightarrow h \tilde a, Z \tilde a$ and $\tilde h^\pm \rightarrow W^\pm \tilde a$. The scale of axion physics, $f$, is predicted to be in the range $(3\times10^8 - 10^{12})$ GeV and, over much of this range, can be extracted from the decay length.


title: In-Situ Low-Angle Cross Sectioning: Bevel Slope Flattening due to Self-Alignment Effects

abstract: Low-angle cross sections are produced inside an Auger microprobe using the equipped depth profile ion sputter gun. Simply the sample is partly covered by a mask. Utilizing the edge of this mask the sample is sputtered with ions. Due to the shading of the mask a cross section is produced in the sample. The slope of this cross section is considerably shallower than given by the geometrical setup. This is attributed to self-alignment effects, which are due to missing sputter cascades in the transition area between sputtered and shaded sample regions and a chamfering of the mask edge. These self-alignment effects are studied here using a 104.6 nm thick SiO2 layer thermally grown on a Si substrate. In this study on one hand for a fixed ion impact angle of 15.8{\deg} as function of the sputter time several in-situ low-angle cross sections were produced. This way slope angles between an ultimate low slope angle of 0.014{\deg} and 0.085{\deg} were achieved. On the other hand for a fixed sputter time the ion impact angle was varied between 14.8{\deg}and 70.8{\deg}. For these samples cross section slope angles between 0.031{\deg}and 0.32{\deg}are observed. These results demonstrate the distinct slope flatting of in-situ cross sectioning.


title: Tides or dark matter sub-halos: Which ones are more attractive?

abstract: Young tidal dwarf galaxies (TDGs) are observed in the tidal debris of gas-rich interacting galaxies. In contrast to what is generally assumed to be the case for isolated dwarf galaxies, TDGs are not embedded in their own dark matter (DM) sub-halo. Hence, they are more sensitive to stellar feedback and could be disrupted on a short time-scale. Detailed numerical and observational studies demonstrate that isolated DM-dominated dwarf galaxies can have lifetimes of more than 10 Gyr. For TDGs that evolve in a tidal field with compressing accelerations equal to the gravitational acceleration within a DM sub-halo typical of an isolated dwarf galaxy, a similar survival time is expected. The tidal acceleration profile depends on the virial mass of the host galaxy and the distance between the TDG and its host. We analytically compare the tidal compression to the gravitational acceleration due to either cuspy or cored DM sub-halos of various virial masses. For example, the tidal field at a distance of 100 kpc to a host halo of 10^13 Msol can be as stabilizing as a 10^9 Msol DM sub-halo. By linking the tidal field to the equivalent gravitational field of a DM sub-halo, we can use existing models of isolated dwarfs to estimate the survivability of TDGs. We show that part of the unexpectedly high dynamical masses inferred from observations of some TDGs can be explained by tidal compression and hence TDGs require to contain less unobservable matter to understand their rotation curves.


title: The Sloan Digital Sky Survey Reverberation Mapping Project: Post-Starburst Signatures in Quasar Host Galaxies at z < 1

abstract: Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.


title: Temporal Evolution of the Size and Temperature of Betelgeuse's Extended Atmosphere

abstract: We use the Very Large Array (VLA) in the A configuration with the Pie Town (PT) Very Long Baseline Array (VLBA) antenna to spatially resolve the extended atmosphere of Betelgeuse over multiple epochs at 0.7, 1.3, 2.0, 3.5, and 6.1 cm. The extended atmosphere deviates from circular symmetry at all wavelengths while at some epochs we find possible evidence for small pockets of gas significantly cooler than the mean global temperature. We find no evidence for the recently reported e-MERLIN radio hotspots in any of our multi-epoch VLA/PT data, despite having sufficient spatial resolution and sensitivity at short wavelengths, and conclude that these radio hotspots are most likely interferometric artefacts. The mean gas temperature of the extended atmosphere has a typical value of 3000 K at 2 $R_{\star}$ and decreases to 1800 K at 6 $R_{\star}$, in broad agreement with the findings of the single epoch study from Lim et al. (1998). The overall temperature profile of the extended atmosphere between $2 R_{\star} \lesssim r \lesssim 6 R_{\star}$ can be described by a power law of the form $T_{\mathrm{gas}}(r) \propto r^{-0.6}$, with temporal variability of a few 100 K evident at some epochs. Finally, we present over 12 years of V band photometry, part of which overlaps our multi-epoch radio data. We find a correlation between the fractional flux density variability at V band with most radio wavelengths. This correlation is likely due to shock waves induced by stellar pulsations, which heat the inner atmosphere and ionize the more extended atmosphere through radiative means. Stellar pulsations may play an important role in exciting Betelgeuse's extended atmosphere.


title: Highly Multiplexible Thermal Kinetic Inductance Detectors for X-Ray Imaging Spectroscopy

abstract: For X-ray imaging spectroscopy, high spatial resolution over a large field of view is often as important as high energy resolution, but current X-ray detectors do not provide both in the same device. Thermal Kinetic Inductance Detectors (TKIDs) are being developed as they offer a feasible way to combine the energy resolution of transition edge sensors with pixel counts approaching CCDs and thus promise significant improvements for many X-ray spectroscopy applications. TKIDs are a variation of Microwave Kinetic Inductance Detectors (MKIDs) and share their multiplexibility: working MKID arrays with 2024 pixels have recently been demonstrated and much bigger arrays are under development. In this work, we present our first working TKID prototypes which are able to achieve an energy resolution of 75 eV at 5.9 keV, even though their general design still has to be optimized. We further describe TKID fabrication, characterization, multiplexing and working principle and demonstrate the necessity of a data fitting algorithm in order to extract photon energies. With further design optimizations we expect to be able to improve our TKID energy resolution to less than 10 eV at 5.9 keV.


title: Beurling-Type Invariant Subspaces of the Poletsky-Stessin Hardy Spaces in the Bidisc

abstract: The invariant subspaces of the Hardy space on $H^2(\mathbb{D})$ of the unit disc are very well known however in several variables the structure of the invariant subspaces of the classical Hardy spaces is not yet fully understood. In this study we examine the invariant subspace problem for Poletsky-Stessin Hardy spaces which is a natural generalization of the classical Hardy spaces to hyperconvex domains in $\mathbb{C}^n$. We showed that not all invariant subspaces of $H^{2}_{\tilde{u}}(\mathbb{D}^2)$ are of Beurling-type. To characterize the Beurling-type invariant subspaces of this space we first generalized the Lax-Halmos theorem of vector valued Hardy spaces to the vector valued Poletsky-Stessin Hardy spaces and then we give a necessary and sufficient condition for the invariant subspaces of $H^{2}_{\tilde{u}}(\mathbb{D}^2)$ to be of Beurling-type.


title: Harnack inequalities and Gaussian estimates for random walks on metric measure spaces

abstract: We characterize Gaussian estimates for transition probability of a discrete time Markov chain in terms of geometric properties of the underlying state space. In particular, we show that the following are equivalent: (1) Two sided Gaussian bounds on heat kernel (2) A scale invariant Parabolic Harnack inequality (3) Volume doubling property and a scale invariant Poincar\'{e} inequality. The underlying state space is a metric measure space, a setting that includes both manifolds and graphs as special cases. An important feature of our work is that our techniques are robust to small perturbations of the underlying space and the Markov kernel. In particular, we show the stability of the above properties under quasi-isometries. We discuss various applications and examples.


title: Global Optimality in Tensor Factorization, Deep Learning, and Beyond

abstract: Techniques involving factorization are found in a wide range of applications and have enjoyed significant empirical success in many fields. However, common to a vast majority of these problems is the significant disadvantage that the associated optimization problems are typically non-convex due to a multilinear form or other convexity destroying transformation. Here we build on ideas from convex relaxations of matrix factorizations and present a very general framework which allows for the analysis of a wide range of non-convex factorization problems - including matrix factorization, tensor factorization, and deep neural network training formulations. We derive sufficient conditions to guarantee that a local minimum of the non-convex optimization problem is a global minimum and show that if the size of the factorized variables is large enough then from any initialization it is possible to find a global minimizer using a purely local descent algorithm. Our framework also provides a partial theoretical justification for the increasingly common use of Rectified Linear Units (ReLUs) in deep neural networks and offers guidance on deep network architectures and regularization strategies to facilitate efficient optimization.


title: A giant comet-like cloud of hydrogen escaping the warm Neptune-mass exoplanet GJ 436b

abstract: Exoplanets orbiting close to their parent stars could lose some fraction of their atmospheres because of the extreme irradiation. Atmospheric mass loss primarily affects low-mass exoplanets, leading to suggest that hot rocky planets might have begun as Neptune-like, but subsequently lost all of their atmospheres; however, no confident measurements have hitherto been available. The signature of this loss could be observed in the ultraviolet spectrum, when the planet and its escaping atmosphere transit the star, giving rise to deeper and longer transit signatures than in the optical spectrum. Here we report that in the ultraviolet the Neptune-mass exoplanet GJ 436b (also known as Gliese 436b) has transit depths of 56.3 +/- 3.5% (1 sigma), far beyond the 0.69% optical transit depth. The ultraviolet transits repeatedly start ~2 h before, and end >3 h after the ~1 h optical transit, which is substantially different from one previous claim (based on an inaccurate ephemeris). We infer from this that the planet is surrounded and trailed by a large exospheric cloud composed mainly of hydrogen atoms. We estimate a mass-loss rate in the range of ~10^8-10^9 g/s, which today is far too small to deplete the atmosphere of a Neptune-like planet in the lifetime of the parent star, but would have been much greater in the past.


title: Plasmon-soliton waves in planar slot waveguides: I. Modeling

abstract: We present two complementary models to study stationary nonlinear solutions in one-dimensional plasmonic slot waveguides made of a finite-thickness nonlinear dielectric core surrounded by metal regions. The considered nonlinearity is of focusing Kerr type. In the first model, it is assumed that the nonlinear term depends only on the transverse component of the electric field and that the nonlinear refractive index change is small compared to the linear part of the refractive index. This first model allows us to describe analytically the field profiles in the whole waveguide using Jacobi elliptic special functions. It also provides a closed analytical formula for the nonlinear dispersion relation. In the second model, the full dependency of the Kerr nonlinearity on the electric field components is taken into account and no assumption is required on the amplitude of the nonlinear term. The disadvantage of this approach is that the field profiles must be computed numerically. Nevertheless analytical constraints are obtained to reduce the parameter space where the solutions of the nonlinear dispersion relations are sought.


title: A superconvergent HDG method for the Incompressible Navier-Stokes Equations on general polyhedral meshes

abstract: We present a superconvergent hybridizable discontinuous Galerkin (HDG) method for the steady-state incompressible Navier-Stokes equations on general polyhedral meshes. For arbitrary conforming polyhedral mesh, we use polynomials of degree k+1, k, k to approximate the velocity, velocity gradient and pressure, respectively. In contrast, we only use polynomials of degree k to approximate the numerical trace of the velocity on the interfaces. Since the numerical trace of the velocity field is the only globally coupled unknown, this scheme allows a very efficient implementation of the method. For the stationary case, and under the usual smallness condition for the source term, we prove that the method is well defined and that the global L2-norm of the error in each of the above-mentioned variables and the discrete H1-norm of the error in the velocity converge with the order of k+1 for k>=0. We also show that for k>=1, the global L2-norm of the error in velocity converges with the order of k+2. From the point of view of degrees of freedom of the globally coupled unknown: numerical trace, this method achieves optimal convergence for all the above-mentioned variables in L2-norm for k>=0, superconvergence for the velocity in the discrete H1-norm without postprocessing for k>=0, and superconvergence for the velocity in L2-norm without postprocessing for k>=1.


title: Elementary matrix reduction over locally stable rings

abstract: A commutative ring R is locally stable provided that for any $a,b\in R$ such that $aR+bR=R$, there exist some $y\in R$ such that $R/(a+by)R$ has stable range 1.For a Bezout ring $R$, we prove that $R$ is an elementary divisor ring if and only if $R$ is locally stable if and only if $R$ has neat range 1.


title: Learning Representations from Deep Networks Using Mode Synthesizers

abstract: Deep learning Networks play a crucial role in the evolution of a vast number of current machine learning models for solving a variety of real world non-trivial tasks. Such networks use big data which is generally unlabeled unsupervised and multi-layered requiring no form of supervision for training and learning data and has been used to successfully build automatic supervisory neural networks. However the question still remains how well the learned data represents interestingness, and their effectiveness i.e. efficiency in deep learning models or applications. If the output of a network of deep learning models can be beamed unto a scene of observables, we could learn the variational frequencies of these stacked networks in a parallel and distributive way.This paper seeks to discover and represent interesting patterns in an efficient and less complex way by incorporating the concept of Mode synthesizers in the deep learning process models


title: Optical Mass Flow Diagnostics in Herbig Ae/Be Stars

abstract: We examine a broad range of mass flow diagnostics in a large sample of Herbig Ae/Be stars (HAEBES) using high resolution optical spectra. The H-beta and He I 5876 angstrom lines show the highest incidence of P-Cygni (30%) and inverse P-Cygni (14%) morphologies, respectively. The Fe II 4924 angstrom line also shows a large incidence of P-Cygni profiles (11%). We find support for many of the conclusions reached in a study based on the analysis of the He I 10830 angstrom line in a large sample of HAEBES. Namely, HAEBES exhibit smaller fractions of both blue-shifted absorption (i.e. mass outflow) and red-shifted absorption (i.e. mass infall or accretion) than their lower mass cousins, the classical T Tauri stars (CTTSs). In particular, the optical data supports the conclusion that HAEBES displaying red-shifted absorption, in general, show maximum red-shifted absorption velocities that are smaller fractions of their stellar escape velocities than is found for CTTSs. This suggests that HAEBE accretion flows are originating deeper in the gravitational potentials of their stars than in CTTS systems. In addition, we find a lack of inner disk wind signatures in the blue-shifted absorption objects; only stellar wind signatures are clearly observed. These findings, along with the lack of detected magnetic fields around HAEBES, support the idea that large magnetospheres are not prevalent around HAEBES and that accretion flows are instead mediated by significantly smaller magnetospheres with relatively smaller truncation radii (e.g. 1-2 stellar radii). Red-shifted absorption is much more common around Herbig Ae stars than Be stars, suggesting that Herbig Be stars may accrete via a boundary layer rather than along magnetic field lines.


title: Spin from defects in two-dimensional quantum field theory

abstract: We build two-dimensional quantum field theories on spin surfaces starting from theories on oriented surfaces with networks of topological defect lines and junctions. The construction uses a combinatorial description of the spin structure in terms of a triangulation equipped with extra data. The amplitude for the spin surfaces is defined to be the amplitude for the underlying oriented surface together with a defect network dual to the triangulation. Independence of the triangulation and of the other choices follows if the line defect and junctions are obtained from a Delta-separable Frobenius algebra with involutive Nakayama automorphism in the monoidal category of topological defects. For rational conformal field theory we can give a more explicit description of the defect category, and we work out two examples related to free fermions in detail: the Ising model and the so(n) WZW model at level 1.


title: Efficient quasi-monoenergetic ion beams up to 18 MeV/nucleon via self-generated plasma fields in relativistic laser plasmas

abstract: Table-top laser-plasma ion accelerators seldom achieve narrow energy spreads, and never without serious compromises in efficiency, particle yield, etc. Using massive computer simulations, we identify a self-organizing scheme that exploits persisting self-generated plasma electric (~TV/m) and magnetic (~10${}^{4}$ Tesla) fields to reduce the ion energy spread after the laser exits the plasma - separating the ion acceleration from the energy spread reduction. Consistent with the scheme, we experimentally demonstrate aluminum and carbon ion beams with narrow spectral peaks at energies up to 310 MeV (11.5 MeV/nucleon) and 220 MeV (18.3 MeV/nucleon), respectively, with high conversion efficiency (~5%, i.e., 4J out of 80J laser). This is achieved with 0.12 PW high-contrast Gaussian laser pulses irradiating planar foils with optimal thicknesses of up to 250 nm that scale with laser intensity. When increasing the focused laser intensity fourfold (by reducing the focusing optic f/number twofold), the spectral-peak energy increases twofold. These results pave the way for next generation compact accelerators suitable for applications. For example, 400 MeV (33.3 MeV/nucleon) carbon-ion beam with narrow energy spread required for ion fast ignition could be generated using PW-class lasers.


title: Approximation of conformal mappings and novel applications to shape recognition of planar domains

abstract: Our goal is to provide a novel method of representing 2D shapes, where each shape will be assigned a unique fingerprint - a computable approximation to a conformal map of the given shape to a canonical shape in 2D or 3D space (see page 22 for a few examples). In this paper, we make the first significant step in this program where we address the case of simply, and doubly-connected planar domains. We prove uniform convergence of our approximation scheme to the appropriate conformal mapping. Along the way, we affirm a conjecture raised by Ken Stephenson in the 90's which predicts that the Riemann mapping can be approximated by a sequence of electrical networks. In fact, we first treat a more general case. Consider a planar annulus, i.e., a bounded, two-connected, Jordan domain, endowed with a sequence of triangulations exhausting it. We construct a corresponding sequence of maps which converge uniformly on compact subsets of the domain, to a conformal homeomorphism onto the interior of a Euclidean annulus bounded by two concentric circles. The resolution of Stephenson's Conjecture then follows by a limiting argument. Our methods involve harmonic mappings and boundary value problems: discrete and analytic in the plane. The scheme we constructed is programmable.


title: Connected components of irreducible maps and 1D quantum phases

abstract: We investigate elementary topological properties of sets of completely positive (CP) maps that arise in quantum Perron-Frobenius theory. We prove that the set of primitive CP maps of fixed Kraus rank is path-connected and we provide a complete classification of the connected components of irreducible CP maps at given Kraus rank and fixed peripheral spectrum in terms of a multiplicity index. These findings are then applied to analyse 1D quantum phases by studying equivalence classes of translational invariant Matrix Product States that correspond to the connected components of the respective CP maps. Our results extend the previously obtained picture in that they do not require blocking of physical sites, they lead to analytic paths and they allow to decompose into ergodic components and to study the breaking of translational symmetry.


title: On Soft Theorems And Form Factors In N=4 SYM Theory

abstract: Soft theorems for the form factors of 1/2-BPS and Konishi operator supermultiplets are derived at tree level in N=4 SYM theory. They have a form identical to the one in the amplitude case. For MHV sectors of stress tensor and Konishi supermultiplets loop corrections to soft theorems are considered at one loop level. They also appear to have universal form in soft limit. Possible generalization of the on-shell diagrams to the form factors based on leading soft behavior is suggested. Finally, we give some comments on inverse soft limit and integrability of form factors in the limit $q^2\to 0$


title: Splash: User-friendly Programming Interface for Parallelizing Stochastic Algorithms

abstract: Stochastic algorithms are efficient approaches to solving machine learning and optimization problems. In this paper, we propose a general framework called Splash for parallelizing stochastic algorithms on multi-node distributed systems. Splash consists of a programming interface and an execution engine. Using the programming interface, the user develops sequential stochastic algorithms without concerning any detail about distributed computing. The algorithm is then automatically parallelized by a communication-efficient execution engine. We provide theoretical justifications on the optimal rate of convergence for parallelizing stochastic gradient descent. Splash is built on top of Apache Spark. The real-data experiments on logistic regression, collaborative filtering and topic modeling verify that Splash yields order-of-magnitude speedup over single-thread stochastic algorithms and over state-of-the-art implementations on Spark.


title: Measuring the Luttinger liquid parameter with shot noise

abstract: We explore the low-frequency noise of interacting electrons in a one-dimensional structure (quantum wire or interaction-coupled edge states) with counterpropagating modes, assuming a single channel in each direction. The system is driven out of equilibrium by a quantum point contact (QPC) with an applied voltage, which induces a double-step energy distribution of incoming electrons on one side of the device. A second QPC serves to explore the statistics of outgoing electrons. We show that measurement of a low-frequency noise in such a setup allows one to extract the Luttinger liquid constant $K$ which is the key parameter characterizing an interacting 1D system. We evaluate the dependence of the zero-frequency noise on $K$ and on parameters of both QPCs (transparencies and voltages).


title: Double-jump stochastic volatility model for VIX: evidence from VVIX

abstract: The paper studies the continuous-time dynamics of VIX with stochastic volatility and jumps in VIX and volatility. Built on the general parametric affine model with stochastic volatility and jump in logarithm of VIX, we derive a linear relation between the stochastic volatility factor and VVIX index. We detect the existence of co-jump of VIX and VVIX and put forward a double-jump stochastic volatility model for VIX through its joint property with VVIX. With VVIX index as a proxy for the stochastic volatility, we use MCMC method to estimate the dynamics of VIX. Comparing nested models on VIX, we show the jump in VIX and the volatility factor is statistically significant. The jump intensity is also statedependent. We analyze the impact of jump factor on the VIX dynamics.


title: An Evolve-Then-Filter Regularized Reduced Order Model

abstract: In this paper, we propose a new evolve-then-filter reduced order model (EF-ROM). This is a regularized ROM (Reg-ROM), which aims at the numerical stabilization of proper orthogonal decomposition (POD) ROMs for convection-dominated flows. We also consider the Leray ROM (L-ROM). These two Reg-ROMs use explicit ROM spatial filtering to smooth (regularize) various terms in the ROMs. Two spatial filters are used: a POD projection onto a POD subspace (Proj) and a new POD differential filter (DF). The four Reg-ROM/filter combinations are tested in the numerical simulation of the three-dimensional flow past a circular cylinder at a Reynolds number $Re=1000$. Overall, the most accurate Reg-ROM/filter combination is EF-ROM-DF. Furthermore, the spatial filter has a higher impact on the Reg-ROM than the regularization used. Indeed, the DF generally yields better results than Proj for both the EF-ROM and L-ROM. Finally, the CPU times of the four Reg-ROM/filter combinations are orders of magnitude lower than the CPU time of the DNS.


title: Effects of anisotropy and disorder on the conductivity of Weyl semimetals

abstract: We study dc conductivity of a Weyl semimetal with uniaxial anisotropy (Fermi velocity ratio $\xi= v_\bot/v_\parallel\neq1$) considering the scattering of charge carriers by a wide class of impurity potentials, both short- and long-range. We obtain the ratio of transverse and longitudinal (with respect to the anisotropy axis) conductivities as a function of both $\xi$ and temperature. We find that the transverse and longitudinal conductivities exhibit different temperature dependence in the case of short-range disorder. For general long-range disorder, the temperature dependence ($\sim T^4$) of the conductivity turns out to be insensitive of the anisotropy in the limits of strong ($\xi\gg$ and $\ll1$) and weak ($\xi\approx1$) anisotropy.


title: The WZW term of the M5-brane and differential cohomotopy

abstract: We combine rational homotopy theory and higher Lie theory to describe the WZW term in the M5-brane sigma model. We observe that this term admits a natural interpretation as a twisted 7-cocycle on super-Minkowski spacetime with coefficients in the rational 4-sphere. This exhibits the WZW term as an element in twisted cohomology, with the twist given by the cocycle of the M2-brane. We consider integration of this rational situation to differential cohomology and differential cohomotopy.


title: Supernatural analogues of Beilinson monads

abstract: We use supernatural bundles to build GL-equivariant resolutions supported on the diagonal of P^n x P^n, in a way that extends Beilinson's resolution of the diagonal. We thus obtain results about supernatural bundles that largely parallel known results about exceptional collections. We apply this construction to Boij-S\"oderberg decompositions of cohomology tables of vector bundles, yielding a proof of concept for the idea that those positive rational decompositions should admit meaningful categorifications.


title: Simulation of Crab Waist Collisions In DA{\Phi}NE With KLOE-2 Interaction Region

abstract: After the successful completion of the SIDDHARTA experiment run with crab waist collisions, the electron- positron collider DA{\Phi}NE has started routine operations for the KLOE-2 detector. The new interaction region also exploits the crab waist collision scheme, but features certain complications including the experimental detector solenoid, compensating anti-solenoids, and tilted quadrupole magnets. We have performed simulations of the beam-beam collisions in the collider taking into account the real DA{\Phi}NE nonlinear lattice. In particular, we have evaluated the effect of crab waist sextupoles and beam-beam interactions on the DA{\Phi}NE dynamical aperture and energy acceptance, and estimated the luminosity that can be potentially achieved with and without crab waist sextupoles in the present working conditions. A numerical analysis has been performed in order to propose possible steps for further luminosity increase in DA{\Phi}NE such as a better working point choice, crab sextupole strength optimization, correction of the phase advance between the sextupoles and the interaction region. The proposed change of the e- ring working point was implemented and resulted in a significant performance increase.


title: Modulational Instability in the Whitham Equation with Surface Tension and Vorticity

abstract: We study modulational stability and instability in the Whitham equation, combining the dispersion relation of water waves and a nonlinearity of the shallow water equations, and modified to permit the effects of surface tension and constant vorticity. When the surface tension coefficient is large, we show that a periodic traveling wave of sufficiently small amplitude is unstable to long wavelength perturbations if the wave number is greater than a critical value, and stable otherwise, similarly to the Benjamin-Feir instability of gravity waves. In the case of weak surface tension, we find intervals of stable and unstable wave numbers, whose boundaries are associated with the extremum of the group velocity, the resonance between the first and second harmonics, the resonance between long and short waves, and a resonance between dispersion and the nonlinearity. For each constant vorticity we show that a periodic traveling wave of sufficiently small amplitude is unstable if the wave number is greater than a critical value, and stable otherwise. Moreover it can be made stable for a sufficiently large vorticity. The results agree with those based upon numerical computations or formal multiple-scale expansions for the physical problem.


title: Immersed Boundary Smooth Extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods

abstract: The Immersed Boundary method is a simple, efficient, and robust numerical scheme for solving PDE in general domains, yet it only achieves first-order spatial accuracy near embedded boundaries. In this paper, we introduce a new high-order numerical method which we call the Immersed Boundary Smooth Extension (IBSE) method. The IBSE method achieves high-order accuracy by smoothly extending the unknown solution of the PDE from a given smooth domain to a larger computational domain, enabling the use of simple Cartesian-grid discretizations (e.g. Fourier spectral methods). The method preserves much of the flexibility and robustness of the original IB method. In particular, it requires minimal geometric information to describe the boundary and relies only on convolution with regularized delta-functions to communicate information between the computational grid and the boundary. We present a fast algorithm for solving elliptic equations, which forms the basis for simple, high-order implicit-time methods for parabolic PDE and implicit-explicit methods for related nonlinear PDE. We apply the IBSE method to solve the Poisson, heat, Burgers', and Fitzhugh-Nagumo equations, and demonstrate fourth-order pointwise convergence for Dirichlet problems and third-order pointwise convergence for Neumann problems.


title: Numerical analysis of parasitic crossing compensation with wires in DA{\Phi}NE

abstract: Current bearing wire compensators were successfully used in the 2005-2006 run of the DA{\Phi}NE collider to mitigate the detrimental effects of parasitic beam-beam interactions. A marked improvement of the positron beam lifetime was observed in machine operation with the KLOE detector. In view of the possible application of wire beam-beam compensators for the High Luminosity LHC upgrade, we revisit the DA{\Phi}NE experiments. We use an improved model of the accelerator with the goal to validate the modern simulation tools and provide valuable input for the LHC upgrade project.


title: Benchmarking Machine Learning Technologies for Software Defect Detection

abstract: Machine Learning approaches are good in solving problems that have less information. In most cases, the software domain problems characterize as a process of learning that depend on the various circumstances and changes accordingly. A predictive model is constructed by using machine learning approaches and classified them into defective and non-defective modules. Machine learning techniques help developers to retrieve useful information after the classification and enable them to analyse data from different perspectives. Machine learning techniques are proven to be useful in terms of software bug prediction. This study used public available data sets of software modules and provides comparative performance analysis of different machine learning techniques for software bug prediction. Results showed most of the machine learning methods performed well on software bug datasets.


title: Spectral likelihood expansions for Bayesian inference

abstract: A spectral approach to Bayesian inference is presented. It pursues the emulation of the posterior probability density. The starting point is a series expansion of the likelihood function in terms of orthogonal polynomials. From this spectral likelihood expansion all statistical quantities of interest can be calculated semi-analytically. The posterior is formally represented as the product of a reference density and a linear combination of polynomial basis functions. Both the model evidence and the posterior moments are related to the expansion coefficients. This formulation avoids Markov chain Monte Carlo simulation and allows one to make use of linear least squares instead. The pros and cons of spectral Bayesian inference are discussed and demonstrated on the basis of simple applications from classical statistics and inverse modeling.


title: Simple commutative algebras in Deligne's categories Rep($S_t$)

abstract: We show that in the Deligne categories $\mathrm{Rep}(S_t)$ for $t$ a transcendental number, the only simple algebra objects are images of simple algebras in the category of representations of a symmetric group under a canonical induction functor. They come in families which interpolate the families of algebras of functions on the cosets of $H\times S_{n-k}$ in $S_n$, for a fixed subgroup $H$ of $S_k$.


title: Optimize Unsynchronized Garbage Collection in an SSD Array

abstract: Solid state disks (SSDs) have advanced to outperform traditional hard drives significantly in both random reads and writes. However, heavy random writes trigger fre- quent garbage collection and decrease the performance of SSDs. In an SSD array, garbage collection of individ- ual SSDs is not synchronized, leading to underutilization of some of the SSDs. We propose a software solution to tackle the unsyn- chronized garbage collection in an SSD array installed in a host bus adaptor (HBA), where individual SSDs are exposed to an operating system. We maintain a long I/O queue for each SSD and flush dirty pages intelligently to fill the long I/O queues so that we hide the performance imbalance among SSDs even when there are few parallel application writes. We further define a policy of select- ing dirty pages to flush and a policy of taking out stale flush requests to reduce the amount of data written to SSDs. We evaluate our solution in a real system. Experi- ments show that our solution fully utilizes all SSDs in an array under random write-heavy workloads. It improves I/O throughput by up to 62% under random workloads of mixed reads and writes when SSDs are under active garbage collection. It causes little extra data writeback and increases the cache hit rate.


title: Entropy, stability, and harmonic map flow

abstract: Inspired by work of Colding-Minicozzi on mean curvature flow, Zhang introduced a notion of entropy stability for harmonic map flow. We build further upon this work in several directions. First we prove the equivalence of entropy stability with a more computationally tractable $\mathcal F$-stability. Then, focusing on the case of spherical targets, we prove a general instability result for high-entropy solitons. Finally, we exploit results of Lin-Wang to observe long time existence and convergence results for maps into certain convex domains and how they relate to generic singularities of harmonic map flow.


title: Towards Resistance Sparsifiers

abstract: We study resistance sparsification of graphs, in which the goal is to find a sparse subgraph (with reweighted edges) that approximately preserves the effective resistances between every pair of nodes. We show that every dense regular expander admits a $(1+\epsilon)$-resistance sparsifier of size $\tilde O(n/\epsilon)$, and conjecture this bound holds for all graphs on $n$ nodes. In comparison, spectral sparsification is a strictly stronger notion and requires $\Omega(n/\epsilon^2)$ edges even on the complete graph. Our approach leads to the following structural question on graphs: Does every dense regular expander contain a sparse regular expander as a subgraph? Our main technical contribution, which may of independent interest, is a positive answer to this question in a certain setting of parameters. Combining this with a recent result of von Luxburg, Radl, and Hein~(JMLR, 2014) leads to the aforementioned resistance sparsifiers.


title: Dynamical History of the Local Group in $\Lambda$CDM

abstract: The positions and velocities of galaxies in the Local Group (LG) measure the gravitational field within it. This is mostly due to the Milky Way (MW) and Andromeda (M31). We constrain their masses using distance and radial velocity (RV) measurements of 32 LG galaxies. To do this, we follow the trajectories of many simulated particles starting on a pure Hubble flow at redshift 9. For each observed galaxy, we obtain a trajectory which today is at the same position. Its final velocity is the model prediction for the velocity of that galaxy. Unlike previous simulations based on spherical symmetry, ours are axisymmetric and include gravity from Centaurus A. We find the total LG mass is ${4.33^{+0.37}_{-0.32}\times{10}^{12}M_\odot}$, with $0.14 \pm 0.07$ of this being in the MW. We approximately account for IC 342, M81, the Great Attractor and the Large Magellanic Cloud. No plausible set of initial conditions yields a good match to the RVs of our sample of LG galaxies. Observed RVs systematically exceed those predicted by the best-fitting $\Lambda$CDM model, with a typical disagreement of ${45.1^{+7.0}_{-5.7}}$ km/s and a maximum of ${110 \pm 13}$ km/s for DDO 99. Interactions between LG dwarf galaxies can't easily explain this. One possibility is a past close flyby of the MW and M31. This arises in some modified gravity theories but not in $\Lambda$CDM. Gravitational slingshot encounters of material in the LG with either of these massive fast-moving galaxies could plausibly explain why some non-satellite LG galaxies are moving away from us even faster than a pure Hubble flow.


title: Fixed-point structure of low-dimensional relativistic fermion field theories: Universality classes and emergent symmetry

abstract: We investigate a class of relativistic fermion theories in 2<d<4 space-time dimensions with continuous chiral U(Nf)xU(Nf) symmetry. This includes a number of well-studied models, e.g., of Gross-Neveu and Thirring type, in a unified framework. Within the limit of pointlike interactions, the RG flow of couplings reveals a network of interacting fixed points, each of which defines a universality class. A subset of fixed points are "critical fixed points" with one RG relevant direction being candidates for critical points of second-order phase transitions. Identifying invariant hyperplanes of the RG flow and classifying their attractive/repulsive properties, we find evidence for emergent higher chiral symmetries as a function of Nf. For the case of the Thirring model, we discover a new critical flavor number that separates the RG stable large-Nf regime from an intermediate-Nf regime in which symmetry-breaking perturbations become RG relevant. This new critical flavor number has to be distinguished from the chiral-critical flavor number, below which the Thirring model is expected to allow spontaneous chiral symmetry breaking, and its existence offers a resolution to the discrepancy between previous results obtained in the continuum and the lattice Thirring models. Moreover, we find indications for a new feature of universality: details of the critical behavior can depend on additional "spectator symmetries" that remain intact across the phase transition. Implications for the physics of interacting fermions on the honeycomb lattice, for which our theory space provides a simple model, are given.


title: Indoor Location Estimation with Optical-based OFDM Communications

abstract: Visible Light Communication (VLC) using light emitting diodes (LEDs) has been gaining increasing attention in recent years as it is appealing for a wide range of applications such as indoor positioning. Orthogonal frequency division multiplexing (OFDM) has been applied to indoor wireless optical communications in order to mitigate the effect of multipath distortion of the optical channel as well as increasing data rate. In this paper, a novel OFDM VLC system is proposed which can be utilized for both communications and indoor positioning. A positioning algorithm based on power attenuation is used to estimate the receiver coordinates. We further calculate the positioning errors in all the locations of a room and compare them with those using single carrier modulation scheme, i.e., on-off keying (OOK) modulation. We demonstrate that OFDM positioning system outperforms its conventional counterpart. Finally, we investigate the impact of different system parameters on the positioning accuracy of the proposed OFDM VLC system.


title: The equivariant Cuntz semigroup

abstract: We introduce an equivariant version of the Cuntz semigroup, that takes an action of a compact group into account. The equivariant Cuntz semigroup is naturally a semimodule over the representation semiring of the given group. Moreover, this semimodule satisfies a number of additional structural properties. We show that the equivariant Cuntz semigroup, as a functor, is continuous and stable. Moreover, cocycle conjugate actions have isomorphic associated equivariant Cuntz semigroups. One of our main results is an analog of Julg's theorem: the equivariant Cuntz semigroup is canonically isomorphic to the Cuntz semigroup of the crossed product. We compute the induced semimodule structure on the crossed product, which in the abelian case is given by the dual action. As an application of our results, we show that freeness of a compact Lie group action on a compact Hausdorff space can be characterized in terms of a canonically defined map into the equivariant Cuntz semigroup, extending results of Atiyah and Segal for equivariant $K$-theory. Finally, we use the equivariant Cuntz semigroup to classify locally representable actions on direct limits of one-dimensional NCCW-complexes, generalizing work of Handelman and Rossmann.


title: Construction of the Lyapunov spectrum in a chaotic system displaying phase synchronization

abstract: We consider a three-dimensional chaotic system consisting of the suspension of Arnold's cat map coupled with a clock via a weak dissipative interaction. We show that the coupled system displays a synchronization phenomenon, in the sense that the relative phase between the suspension flow and the clock locks to a special value, thus making the motion fall onto a lower dimensional attractor. More specifically, we construct the attractive invariant manifold, of dimension smaller than three, using a convergent perturbative expansion. Moreover, we compute via convergent series the Lyapunov exponents, including notably the central one. The result generalizes a previous construction of the attractive invariant manifold in a similar but simpler model. The main novelty of the current construction relies in the computation of the Lyapunov spectrum, which consists of non-trivial analytic exponents. Some conjectures about a possible smoothening transition of the attractor as the coupling is increased are also discussed.


title: The influence of electron-doping on the ground state of (Sr{1-x}La{x})2IrO4

abstract: The evolution of the electronic properties of electron-doped (Sr{1-x}La{x})2IrO4 is experimentally explored as the doping limit of La is approached. As electrons are introduced, the electronic ground state transitions from a spin-orbit Mott phase into an electronically phase separated state, where long-range magnetic order vanishes beyond x = 0.02 and charge transport remains percolative up to the limit of La substitution (x~0.06). In particular, the electronic ground state remains inhomogeneous even beyond the collapse of the parent state's long-range antiferromagnetic order, while persistent short-range magnetism survives up to the highest La-substitution levels. Furthermore, as electrons are doped into Sr2IrO4, we observe the appearance of a low temperature magnetic glass-like state intermediate to the complete suppression of antiferromagnetic order. Universalities and differences in the electron-doped phase diagrams of single layer and bilayer Ruddlesden-Popper strontium iridates are discussed.


title: $Z'$ resonance and associated $Zh$ production at future Higgs boson factory: ILC and CLIC

abstract: We study the prospects of the B-L model with an additional $Z'$ boson to be a Higgs boson factory at high-energy and high-luminosity linear electron positron colliders, such as the ILC and CLIC, through the Higgs-strahlung process $e^{+}e^{-}\rightarrow (Z, Z') \to Zh$, including both the resonant and non-resonant effects. We evaluate the total cross section of $Zh$ and we calculate the total number of events for integrated luminosities of 500-2000\hspace{0.8mm}$fb^{-1}$ and center of mass energies between 500 and 3000\hspace{0.8mm}$GeV$. We find that the total number of expected $Zh$ events can reach $10^6$, which is a very optimistic scenario and it would be possible to perform precision measurements for both the $Z'$ and Higgs boson in future high-energy $e^+e^-$ colliders experiments.


title: Towards a Grand Unified Theory of Mathematics and Physics

abstract: Wigner's "unreasonable effectiveness of mathematics" in physics can be understood as a reflection of a deep and unexpected unity between the fundamental structures of mathematics and of physics. Some of the history of evidence for this is reviewed, emphasizing developments since Wigner's time and still poorly understood analogies between number theory and quantum field theory.


title: Ebb: A DSL for Physical Simulation on CPUs and GPUs

abstract: Designing programming environments for physical simulation is challenging because simulations rely on diverse algorithms and geometric domains. These challenges are compounded when we try to run efficiently on heterogeneous parallel architectures. We present Ebb, a domain-specific language (DSL) for simulation, that runs efficiently on both CPUs and GPUs. Unlike previous DSLs, Ebb uses a three-layer architecture to separate (1) simulation code, (2) definition of data structures for geometric domains, and (3) runtimes supporting parallel architectures. Different geometric domains are implemented as libraries that use a common, unified, relational data model. By structuring the simulation framework in this way, programmers implementing simulations can focus on the physics and algorithms for each simulation without worrying about their implementation on parallel computers. Because the geometric domain libraries are all implemented using a common runtime based on relations, new geometric domains can be added as needed, without specifying the details of memory management, mapping to different parallel architectures, or having to expand the runtime's interface. We evaluate Ebb by comparing it to several widely used simulations, demonstrating comparable performance to hand-written GPU code where available, and surpassing existing CPU performance optimizations by up to 9$\times$ when no GPU code exists.


title: Thermodynamic perturbation theory for associating fluids confined in a 1- dimensional pore

abstract: In this paper a new theory is developed for the self - assembly of associating molecules confined to a single spatial dimension, but allowed to explore all orientation angles. The interplay of the anisotropy of the pair potential and the low dimensional space, results in orientationally ordered associated clusters. This local order enhances association due to a decrease in orientational entropy. Unlike bulk 3D fluids which are orientationally homogeneous, association in 1D necessitates the self - consistent calculation of the orientational distribution function. To test the new theory, Monte Carlo simulations are performed and the theory is found to be accurate. The theory developed in this paper may be used as a tool to study hydrogen bonding of molecules in 1D zeolites as well as hydrogen bonding of water in carbon nanotubes.


title: Cosmic anisotropic doomsday in Bianchi type I universes

abstract: In order to investigate if the anisotropy of the spacetime may induce future singularities at a finite value of the cosmic time on the evolution of cosmological models, we study vacuum and non-vacuum Bianchi type I spacetimes exhibiting such future singularities. We show that in the case of Kasner vacuum cosmologies the spacetime may rip itself apart in a finite time, and only in the direction corresponding to the unique scale factor which diverges at this finite value of the cosmic time. The other two directional scale factors and the average scale factor do not diverge and tend to zero at this time, while the directional and average expansion rates also become infinite. Due to the absence of the matter content this anisotropic future singularity is induced by the shear scalar, which also blows up at this time. We call such a singularity "Vacuum Rip". For non-vacuum solutions we discuss fully anisotropic Bianchi type I spacetimes filled with a stiff fluid and ellipsoidal (axisymmetric) cosmological models filled with matter with isotropic and anisotropic barotropic pressure, and characterized by $\sigma/\theta=const$, where $\sigma$ and $\theta$ are the shear scalar and the expansion scalar respectively.


title: Moment Representations of the Exceptional $X_1$-Laguerre Orthogonal Polynomials

abstract: Exceptional orthogonal Laguerre polynomials can be viewed as an extension of the classical Laguerre polynomials per excluding polynomials of certain order(s) from being eigenfunctions for the corresponding exceptional differential operator. We are interested in the (so-called) Type I $X_1$-Laguerre polynomial sequence $\{L_n^\alpha\}_{n=1}^\infty$, $\text{deg} \,p_n = n$ and $\alpha>0$, where the constant polynomial is omitted. We derive two representations for the polynomials in terms of moments by using determinants. The first representation in terms of the canonical moments is rather cumbersome. We introduce adjusted moments and find a second, more elegant formula. We deduce a recursion formula for the moments and the adjusted ones. The adjusted moments are also expressed via a generating function. We observe a certain detachedness of the first two moments from the others.


title: Rigidity of Determinantal Point Processes with the Airy, the Bessel and the Gamma Kernel

abstract: A point process is said to be rigid if for any bounded domain in the phase space, the number of particles in the domain is almost surely determined by the restriction of the configuration to the complement of our bounded domain. The main result of this paper is that determinantal point processes with the Airy, the Bessel and the Gamma kernels are rigid. The proof follows the scheme of Ghosh [6], Ghosh and Peres [7]: the main step is the construction of a sequence of additive statistics with variance going to zero.


title: Too dynamic to fail. Empirical support for an autocatalytic model of Minsky's financial instability hypothesis

abstract: Solomon and Golo [1] have recently proposed an autocatalytic (self-reinforcing) feedback model which couples a macroscopic system parameter (the interest rate), a microscopic parameter that measures the distribution of the states of the individual agents (the number of firms in financial difficulty) and a peer-to-peer network effect (contagion across supply chain financing). In this model, each financial agent is characterized by its resilience to the interest rate. Above a certain rate the interest due on the firm's financial costs exceeds its earnings and the firm becomes susceptible to failure (ponzi). For the interest rate levels under a certain threshold level, the firm loans are smaller then its earnings and the firm becomes 'hedge.' In this paper, we fit the historical data (2002-2009) on interest rate data into our model, in order to predict the number of the ponzi firms. We compare the prediction with the data taken from a large panel of Italian firms over a period of 9 years. We then use trade credit linkages to discuss the connection between the ponzi density and the network percolation. We find that the 'top-down'-'bottom-up' positive feedback loop accounts for most of the Minsky crisis accelerator dynamics. The peer-to-peer ponzi companies contagion becomes significant only in the last stage of the crisis when the ponzi density is above a critical value. Moreover the ponzi contagion is limited only to the companies that were not dynamic enough to substitute their distressed clients with new ones. In this respect the data support a view in which the success of the economy depends on substituting the static 'supply-network' picture with an interacting dynamic agents one.


title: Resonant tunneling in a quantum oxide superlattice

abstract: Resonant tunnelling is a quantum mechanical process that has long been attracting both scientific and technological attention owing to its intriguing underlying physics and unique applications for high-speed electronics. The materials system exhibiting resonant tunnelling, however, has been largely limited to the conventional semiconductors, partially due to their excellent crystalline quality. Here we show that a deliberately designed transition metal oxide superlattice exhibits a resonant tunnelling behaviour with a clear negative differential resistance. The tunnelling occurred through an atomically thin, lanthanum {\delta}-doped SrTiO3 layer, and the negative differential resistance was realized on top of the bipolar resistance switching typically observed for perovskite oxide junctions. This combined process resulted in an extremely large resistance ratio (~10^5) between the high and low-resistance states. The unprecedentedly large control found in atomically thin {\delta}-doped oxide superlattices can open a door to novel oxide-based high-frequency logic devices.


title: Synchronization of ad hoc Clock Networks

abstract: We introduce a graph-theoretic approach to synchronizing clocks in an {\em ad hoc} network of $N$~timepieces. Clocks naturally drift away from being synchronized because of many physical factors. The manual way of clock synchronization suffers from an inherrent propagation of the so called "clock drift" due to "word-of-mouth effect." The current standard way of automated clock synchronization is either via radio band transmission of the global clock or via the software-based Network Time Protocol (NTP). Synchronization via radio band transmission suffers from the wave transmission delay, while the client-server-based NTP does not scale to increased number of clients as well as to unforeseen server overload conditions (e.g., flash crowd and time-of-day effects). Further, the trivial running time of NTP for synchronizing an $N$-node network, where each node is a clock and the NTP server follows a single-port communication model, is~$\bigO(N)$. We introduce in this paper a $\bigO(\log N)$ time for synchronizing the clocks in exchange for an increase of $\bigO(N)$ in space complexity, though through creative "tweaking," we later reduced the space requirement to~$\bigO(1)$. Our graph-theoretic protocol assumes that the network is $\K_N$, while the subset of clocks are in an embedded circulant graph $\C_{n<N}^q$ with $q$~jumps and clock information is communicated through circular shifts within the $\C_{n<N}^q$. All $N$~nodes communicate via a single-port duplex channel model. Theoretically, this synchronization protocol allows for $N(\log N)^{-1} - 1$ more synchronizations than the client-server-based one. Empirically through statistically replicated multi-agent-based microsimulation runs, our protocol allows at most 80\% of the clocks synchronized compared to the current protocol which only allows up to 30\% after some steady-state time.


title: Excitation, detection, and electrostatic manipulation of terahertz-frequency range plasmons in a two-dimensional electron system

abstract: Terahertz time domain spectroscopy employing free-space radiation has frequently been used to probe the elementary excitations of low-dimensional systems. The diffraction limit blocks its use for the in-plane study of individual laterally defined nanostructures, however. Here, we demonstrate a planar terahertz-frequency plasmonic circuit in which photoconductive material is monolithically integrated with a two-dimensional electron system. Plasmons with a broad spectral range (up to ~400 GHz) are excited by injecting picosecond-duration pulses, generated and detected by a photoconductive semiconductor, into a high mobility two-dimensional electron system. Using voltage modulation of a Schottky gate overlying the two-dimensional electron system, we form a tuneable plasmonic cavity, and observe electrostatic manipulation of the plasmon resonances. Our technique offers a direct route to access the picosecond dynamics of confined transport in a broad range of lateral nanostructures.


title: Asymptotic behavior of the Kohn-Sham exchange potential at a metal surface

abstract: The asymptotic structure of the Kohn-Sham exchange potential $v_x ({\bf r})$ in the classically forbidden region of a metal surface is investigated, together with that of the Slater exchange potential $V_x^S ({\bf r})$ and those of the approximate Krieger-Li-Iafrate $V_x^{KLI} ({\bf r})$ and Harbola-Sahni $W_x ({\bf r})$ exchange potentials. Particularly, the former is shown to have the form of $v_{x} (z \to \infty) = - \alpha_x / z$ with $\alpha_x$ a constant dependent only of bulk electron density. The same result in previous work is thus confirmed; in the meanwhile controversy raised recently gets resolved. The structure of the exchange hole $\rho_x ({\bf r}, {\bf r}')$ is examined, and the delocalization of it in the metal bulk when the electron is at large distance from the metal surface is demonstrated with analytical expressions. The asymptotic structures of $v_x ({\bf r})$, $V_x^S ({\bf r})$, $V_x^{KLI} ({\bf r})$, and $W_x ({\bf r})$ at a slab metal surface are also investigated. Particularly, $v_{x} (z \to \infty) = - 1/ z$ in the slab case. The distinction in this respect between the semi-infinite and the slab metal surfaces is elucidated.


title: On the set of catenary degrees of finitely generated cancellative commutative monoids

abstract: The catenary degree of an element $n$ of a cancellative commutative monoid $S$ is a nonnegative integer measuring the distance between the irreducible factorizations of $n$. The catenary degree of the monoid $S$, defined as the supremum over all catenary degrees occurring in $S$, has been heavily studied as an invariant of nonunique factorization. In this paper, we investigate the set $\mathsf C(S)$ of catenary degrees achieved by elements of $S$ as a factorization invariant, focusing on the case where $S$ in finitely generated (where $\mathsf C(S)$ is known to be finite). Answering an open question posed by Garc\'ia-S\'anchez, we provide a method to compute the smallest nonzero element of $\mathsf C(S)$ that parallels a well-known method of computing the maximum value. We also give several examples demonstrating certain extremal behavior for $\mathsf C(S)$, and present some open questions for further study.


title: Geometrogenesis under Quantum Graphity: problems with the ripening Universe

abstract: Quantum Graphity (QG) is a model of emergent geometry in which space is represented by a dynamical graph. The graph evolves under the action of a Hamiltonian from a high-energy pre-geometric state to a low-energy state in which geometry emerges as a coarse-grained effective property of space. Here we show the results of numerical modelling of the evolution of the QG Hamiltonian, a process we term "ripening" by analogy with crystallographic growth. We find that the model as originally presented favours a graph composed of small disjoint subgraphs. Such a disconnected space is a poor representation of our universe. A new term is introduced to the original QG Hamiltonian, which we call the hypervalence term. It is shown that the inclusion of a hypervalence term causes a connected lattice-like graph to be favoured over small isolated subgraphs.


title: DCLfix: A Recommendation System for Repairing Architectural Violations

abstract: Architectural erosion is a recurrent problem in software evolution. Despite this fact, the process is usually tackled in ad hoc ways, without adequate tool support at the architecture level. To address this shortcoming, this paper presents a recommendation system -- called DCLfix -- that provides refactoring guidelines for maintainers when tackling architectural erosion. In short, DCLfix suggests refactoring recommendations for violations detected after an architecture conformance process using DCL, an architectural constraint language


title: Ultracompact quantum splitter of degenerate photon pairs

abstract: Integrated sources of indistinguishable photons have attracted a lot of attention because of their applications in quantum communication and optical quantum computing. Here, we demonstrate an ultra-compact quantum splitter for degenerate single photons based on a monolithic chip incorporating Sagnac loop and a micro-ring resonator with a footprint of 0.011 mm2, generating and deterministically splitting indistinguishable photon pairs using time-reversed Hong-Ou-Mandel interference. The ring resonator provides enhanced photon generation rate, and the Sagnac loop ensures the photons travel through equal path lengths and interfere with the correct phase to enable the reversed HOM effect to take place. In the experiment, we observed a HOM dip visibility of 94.5 +- 3.3 %, indicating the photons generated by the degenerate single photon source are in a suitable state for further integration with other components for quantum applications, such as controlled-NOT gates.


title: Block Time Step Storage Scheme for Astrophysical N-body Simulations

abstract: Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used.


title: Tunable frequency-up/down conversion in gas-filled hollow-core photonic crystal fibers

abstract: Based on the interplay between photoionization and Raman effect in gas-filled photonic crystal fibers, we propose a new optical device to control frequency-conversion of ultrashort pulses. By tuning the input-pulse energy, the output spectrum can be either down-converted, up-converted, or even frequency-shift compensated. For low input energies, Raman effect is dominant and leads to a redshift that increases linearly during propagation. For larger pulse energies, photoionization starts to take over the frequency conversion process, and induces a strong blueshift. We have found also that the fiber-output pressure can provide an additional degree of freedom to control the spectrum shift.


title: Matrix elements and duality for type 2 unitary representations of the Lie superalgebra gl(m|n)

abstract: The characteristic identity formalism discussed in our recent articles is further utilized to derive matrix elements of type 2 unitary irreducible $gl(m|n)$ modules. In particular, we give matrix element formulae for all gl(m|n) generators, including the non-elementary generators, together with their phases on finite dimensional type 2 unitary irreducible representations. Remarkably, we find that the type 2 unitary matrix element equations coincide with the type 1 unitary matrix element equations for non-vanishing matrix elements up to a phase.


title: s.Baer and s.Rickart Modules

abstract: In this paper, we study module theoretic definitions of the Baer and related ring concepts. We say a module is s.Baer if the right annihilator of a nonempty subset of the module is generated by an idempotent in the ring. We show that s.Baer modules satisfy a number of closure properties. Under certain conditions, a torsion theory is established for the s.Baer modules, and we provide examples of s.Baer torsion modules and modules with a nonzero s.Baer radical. The other principal interest of this paper is to provide explicit connections between s.Baer modules and projective modules. Among other results, we show that every s.Baer module is an essential extension of a projective module. Additionally, we prove, with limited and natural assumptions, that in a generalized triangular matrix ring every s.Baer submodule of the ring is projective. As an application, we show that every prime ring with a minimal right ideal has the strong summand intersection property. Numerous examples are provided to illustrate, motivate, and delimit the theory.


title: Falconer distance problem, additive energy and Cartesian products

abstract: A celebrated result due to Wolff says if $E$ is a compact subset of ${\Bbb R}^2$, then the Lebesgue measure of the distance set $\Delta(E)=\{|x-y|: x,y \in E \}$ is positive if the Hausdorff dimension of $E$ is greater than $\frac{4}{3}$. In this paper we improve the $\frac{4}{3}$ barrier by a small exponent for Cartesian products. In higher dimensions, also in the context of Cartesian products, we reduce Erdogan's $\frac{d}{2}+\frac{1}{3}$ exponent to $\frac{d^2}{2d-1}$. The proof uses a combination of Fourier analysis and additive comibinatorics.


title: A comparison theorem for the law of large numbers in Banach spaces

abstract: Let $(\mathbf{B}, \|\cdot\|)$ be a real separable Banach space. Let $\{X, X_{n}; n \geq 1\}$ be a sequence of i.i.d. {\bf B}-valued random variables and set $S_{n} = \sum_{i=1}^{n}X_{i},~n \geq 1$. Let $\{a_{n}; n \geq 1\}$ and $\{b_{n}; n \geq 1\}$ be increasing sequences of positive real numbers such that $\lim_{n \rightarrow \infty} a_{n} = \infty$ and $\left\{b_{n}/a_{n};~ n \geq 1 \right\}$ is a nondecreasing sequence. In this paper, we provide a comparison theorem for the law of large numbers for i.i.d. {\bf B}-valued random variables. That is, we show that $\displaystyle \frac{S_{n}- n \mathbb{E}\left(XI\{\|X\| \leq b_{n} \} \right)}{b_{n}} \rightarrow 0$ almost surely (resp. in probability) for every {\bf B}-valued random variable $X$ with $\sum_{n=1}^{\infty} \mathbb{P}(\|X\| > b_{n}) < \infty$ (resp. $\lim_{n \rightarrow \infty}n\mathbb{P}(\|X\| > b_{n}) = 0$) if $S_{n}/a_{n} \rightarrow 0$ almost surely (resp. in probability) for every symmetric {\bf B}-valued random variable $X$ with $\sum_{n=1}^{\infty} \mathbb{P}(\|X\| > a_{n}) < \infty$ (resp. $\lim_{n \rightarrow \infty}n\mathbb{P}(\|X\| > a_{n}) = 0$). To establish this comparison theorem for the law of large numbers, we invoke two tools: 1) a comparison theorem for sums of independent {\bf B}-valued random variables and, 2) a symmetrization procedure for the law of large numbers for sums of independent {\bf B}-valued random variables. A few consequences of our main results are provided.


title: Degenerate Motions in Multicamera Cluster SLAM with Non-overlapping Fields of View

abstract: An analysis of the relative motion and point feature model configurations leading to solution degeneracy is presented, for the case of a Simultaneous Localization and Mapping system using multicamera clusters with non-overlapping fields-of-view. The SLAM optimization system seeks to minimize image space reprojection error and is formulated for a cluster containing any number of component cameras, observing any number of point features over two keyframes. The measurement Jacobian is transformed to expose a reduced-dimension representation such that the degeneracy of the system can be determined by the rank of a dense submatrix. A set of relative motions sufficient for degeneracy are identified for certain cluster configurations, independent of target model geometry. Furthermore, it is shown that increasing the number of cameras within the cluster and observing features across different cameras over the two keyframes reduces the size of the degenerate motion sets significantly.


title: The influence of initial solutions to exact solutions of the (2+1)-dimensional generalized Nizhnik-Novikov-Veselov equations

abstract: The (2+1)-dimensional generalized Nizhnik-Novikov-Veselov equations (GNNVEs) are investigated in order to search the influence of initial solution to exact solutions. The GNNVEs are converted into the combined equations of differently two bilinear forms by means of the homogeneous balance of undetermined coefficients method. Accordingly, the two class of exact N-soliton solutions and three wave solutions are obtained respectively by using the Hirota direct method combined with the simplified version of Hereman and the three wave method. The proposed method is also a standard and computable method, which can be generalized to deal with some nonlinear partial differential equations (NLPDEs).


title: On the Born-Oppenheimer approximation of diatomic molecular resonances

abstract: We give a new reduction of a general diatomic molecular Hamiltonian, without modifying it near the collision set of nuclei. The resulting effective Hamiltonian is the sum of a smooth semiclassical pseudodifferential operator (the semiclassical parameter being the inverse of the square-root of the nuclear mass), and a semibounded operator localised in the elliptic region corresponding to the nuclear collision set. We also study its behaviour on exponential weights, and give several applications where molecular resonances appear and can be well located.


title: Nodal length of Steklov eigenfunctions on real-analytic Riemannian surfaces

abstract: We prove sharp upper and lower bounds for the nodal length of Steklov eigenfunctions on real-analytic Riemannian surfaces with boundary. The argument involves frequency function methods for harmonic functions in the interior of the surface as well as the construction of exponentially accurate approximations for the Steklov eigenfunctions near the boundary.


title: Ultrafast Carrier Dynamics in the Large Magnetoresistance Material WTe$_{2}$

abstract: Ultrafast optical pump-probe spectroscopy is used to track carrier dynamics in the large magnetoresistance material WTe$_{2}$. Our experiments reveal a fast relaxation process occurring on a sub-picosecond time scale that is caused by electron-phonon thermalization, allowing us to extract the electron-phonon coupling constant. An additional slower relaxation process, occurring on a time scale of $\sim$5-15 picoseconds, is attributed to phonon-assisted electron-hole recombination. As the temperature decreases from 300 K, the timescale governing this process increases due to the reduction of the phonon population. However, below $\sim$50 K, an unusual decrease of the recombination time sets in, most likely due to a change in the electronic structure that has been linked to the large magnetoresistance observed in this material.


title: Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations

abstract: Due to the chaotic nature of the Solar System, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner planets, leading to close encounters and/or collisions can be initiated through a large increase in Mercury's eccentricity, with a currently assumed likelihood of ~1%. However, little is known at present about the robustness of this number. Here I report ensemble integrations of the full equations of motion of the eight planets and Pluto over 5 Gyr, including contributions from general relativity. The results show that different numerical algorithms lead to statistically different results for the evolution of Mercury's eccentricity (eM). For instance, starting at present initial conditions (eM ~= 0.21), Mercury's maximum eccentricity achieved over 5 Gyr is on average significantly higher in symplectic ensemble integrations using heliocentricthan Jacobi coordinates and stricter error control. In contrast, starting at a possible future configuration (eM ~= 0.53), Mercury's maximum eccentricity achieved over the subsequent 500 Myr is on average significantly lower using heliocentric than Jacobi coordinates. For example, the probability for eM to increase beyond 0.53 over 500 Myr is >90% (Jacobi) vs. only 40-55% (heliocentric). This poses a dilemma as the physical evolution of the real system - and its probabilistic behavior - cannot depend on the coordinate system or numerical algorithm chosen to describe it. Some tests of the numerical algorithms suggest that symplectic integrators using heliocentric coordinates underestimate the odds for destabilization of Mercury's orbit at high initial eM.


title: Analytic MMSE Bounds in Linear Dynamic Systems with Gaussian Mixture Noise Statistics

abstract: Using state-space representation, mobile object positioning problems can be described as dynamic systems, with the state representing the unknown location and the observations being the information gathered from the location sensors. For linear dynamic systems with Gaussian noise, the Kalman filter provides the Minimum Mean-Square Error (MMSE) state estimation by tracking the posterior. Hence, by approximating non-Gaussian noise distributions with Gaussian Mixtures (GM), a bank of Kalman filters or Gaussian Sum Filter (GSF), can provide the MMSE state estimation. However, the MMSE itself is not analytically tractable. Moreover, the general analytic bounds proposed in the literature are not tractable for GM noise statistics. Hence, in this work, we evaluate the MMSE of linear dynamic systems with GM noise statistics and propose its analytic lower and upper bounds. We provide two analytic upper bounds which are the Mean-Square Errors (MSE) of implementable filters, and we show that based on the shape of the GM noise distributions, the tighter upper bound can be selected. We also show that for highly multimodal GM noise distributions, the bounds and the MMSE converge. Simulation results support the validity of the proposed bounds and their behavior in limits.


title: Mesoscopic correlations in Tb2Ti2O7 spin liquid

abstract: We have studied the spin correlations with $\bf{k}$= ($\frac12$, $\frac12$, $\frac12$) propagation vector which appear below 0.4\, K in \tbti\ spin liquid by combining powder neutron diffraction and specific heat on Tb$_{2+x}$Ti$_{2-x}$O$_{7+y}$ samples with $x$=0, 0.01, -0.01. The $\bf{k}$= ($\frac12$, $\frac12$, $\frac12$) order clearly appears on all neutron patterns by subtracting a pattern at 1.2(1)\,K. Refining the subtracted patterns at 0.07\,K yields two possible spin structures, with spin-ice-like and monopole-like correlations respectively. Mesoscopic correlations involve Tb moments of 1 to 2 \mub\ ordered on a length scale of about 20 \AA. In addition, long range order involving a small spin component of 0.1 to 0.2 \mub\ is detected for the $x$= 0 and 0.01 samples showing a peak in the specific heat. Comparison with previous single crystals data suggests that the ($\frac12$, $\frac12$, $\frac12$) order settles in through nanometric spin textures with dominant spin ice character and correlated orientations, analogous to nanomagnetic twins.


title: Probing Parton Dynamics of QCD Matter with $\Omega$ and $\phi$ Production

abstract: We present measurements of $\Omega$ and $\phi$ production at mid-rapidity from Au+Au collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). Motivated by the coalescence formation mechanism for these strange hadrons, we study the ratios of $N(\Omega^{-}+\Omega^{+})/(2N(\phi))$. These ratios as a function of transverse momentum ($p_T$) fall on a consistent trend at high collision energies, but start to show deviations in peripheral collisions at $\sqrt{s_{NN}}$ = 19.6, 27 and 39 GeV, and in central collisions at 11.5 GeV in the intermediate $p_T$ region of 2.4-3.6 GeV/c. We further evaluate empirically the strange quark $p_T$ distributions at hadronization by studying the $\Omega/\phi$ ratios scaled by the number of constituent quarks. The NCQ-scaled $\Omega/\phi$ ratios show a suppression of strange quark production in central collisions at 11.5 GeV compared to $\sqrt{s_{NN}} >= 19.6$ GeV. The shapes of the presumably thermal strange quark distributions in 0-60% most central collisions at 7.7 GeV show significant deviations from those in 0-10% most central collisions at higher energies. These features suggest that there is likely a change of the underlying strange quark dynamics in the transition from quark-matter to hadronic matter at collision energies below 19.6 GeV.


title: Asymptotic Dynamics of Monopole Walls

abstract: We determine the asymptotic dynamics of the U(N) doubly periodic BPS monopole in Yang-Mills-Higgs theory, called a monopole wall, by exploring its Higgs curve using the Newton polytope and amoeba. In particular, we show that the monopole wall splits into subwalls when any of its moduli become large. The long-distance gauge and Higgs field interactions of these subwalls are abelian, allowing us to derive an asymptotic metric for the monopole wall moduli space.


title: E\"otv\"os Experiments with Supermassive Black Holes

abstract: By examining the locations of central black holes in two elliptical galaxies, M\,32 and M\,87, we derive constraints on the violation of the strong equivalence principle for purely gravitational objects, i.e. black holes, of less than about two-thirds, $\eta_N<0.68$ from the gravitational interaction of M\,87 with its neighbours in the Virgo cluster. Although M\,32 appears to be a good candidate for this technique, the high concentration of stars near its centre substantially weakens the constraints. On the other hand, if a central black hole is found in NGC 205 or one of the other satellite ellipticals of M\,31, substantially better constraints could be obtained. In all cases the constraints could improve dramatically with better astrometry.


title: Amplifying the Impact of Open Access: Wikipedia and the Diffusion of Science

abstract: With the rise of Wikipedia as a first-stop source for scientific knowledge, it is important to compare its representation of that knowledge to that of the academic literature. Here we identify the 250 most heavily used journals in each of 26 research fields (4,721 journals, 19.4M articles in total) indexed by the Scopus database, and test whether topic, academic status, and accessibility make articles from these journals more or less likely to be referenced on Wikipedia. We find that a journal's academic status (impact factor) and accessibility (open access policy) both strongly increase the probability of it being referenced on Wikipedia. Controlling for field and impact factor, the odds that an open access journal is referenced on the English Wikipedia are 47% higher compared to paywall journals. One of the implications of this study is that a major consequence of open access policies is to significantly amplify the diffusion of science, through an intermediary like Wikipedia, to a broad audience.


title: CRAFT: ClusteR-specific Assorted Feature selecTion

abstract: We present a framework for clustering with cluster-specific feature selection. The framework, CRAFT, is derived from asymptotic log posterior formulations of nonparametric MAP-based clustering models. CRAFT handles assorted data, i.e., both numeric and categorical data, and the underlying objective functions are intuitively appealing. The resulting algorithm is simple to implement and scales nicely, requires minimal parameter tuning, obviates the need to specify the number of clusters a priori, and compares favorably with other methods on real datasets.


title: CASSIS: The Cornell Atlas of Spitzer/Infrared Spectrograph Sources. II. High-resolution observations

abstract: The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the high-resolution spectra. The high-resolution observations represent approximately one third of all staring observations performed with the IRS instrument. While low-resolution observations are adapted to faint objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observations allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to minimize these effects through various improvements. A super sampled point-spread function was created in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline selects the best extraction method based on the spatial extent of the object. For unresolved sources, the optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extraction. We have developed several techniques for optimal extraction, including a differential method that eliminates low-level rogue pixels (even when no dedicated background observation was performed). The updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of mapping observations.

categories: astro-ph.IM astro-ph.GA astro-ph.HE astro-ph.SR


title: Joint community and anomaly tracking in dynamic networks

abstract: Most real-world networks exhibit community structure, a phenomenon characterized by existence of node clusters whose intra-edge connectivity is stronger than edge connectivities between nodes belonging to different clusters. In addition to facilitating a better understanding of network behavior, community detection finds many practical applications in diverse settings. Communities in online social networks are indicative of shared functional roles, or affiliation to a common socio-economic status, the knowledge of which is vital for targeted advertisement. In buyer-seller networks, community detection facilitates better product recommendations. Unfortunately, reliability of community assignments is hindered by anomalous user behavior often observed as unfair self-promotion, or "fake" highly-connected accounts created to promote fraud. The present paper advocates a novel approach for jointly tracking communities while detecting such anomalous nodes in time-varying networks. By postulating edge creation as the result of mutual community participation by node pairs, a dynamic factor model with anomalous memberships captured through a sparse outlier matrix is put forth. Efficient tracking algorithms suitable for both online and decentralized operation are developed. Experiments conducted on both synthetic and real network time series successfully unveil underlying communities and anomalous nodes.


title: Asymptotics and exact formulas for Zagier polynomials

abstract: In 1998 Don Zagier introduced the modified Bernoulli numbers $B_{n}^{*}$ and showed that they satisfy amusing variants of some properties of Bernoulli numbers. In particular, he studied the asymptotic behavior of $B_{2n}^{*}$, and also obtained an exact formula for them, the motivation for which came from the representation of $B_{2n}$ in terms of the Riemann zeta function $\zeta(2n)$. The modified Bernoulli numbers were recently generalized to Zagier polynomials $B_{n}^{*}(x)$. For $0<x<1$, an exact formula for $B_{2n}^{*}(x)$ involving infinite series of Bessel function of the second kind and Chebyshev polynomials, that yields Zagier's formula in a limiting case, is established here. Such series arise in diffraction theory. An analogous formula for $B_{2n+1}^{*}(x)$ is also presented. The $6$-periodicity of $B_{2n+1}^{*}$ is deduced as a limiting case of it. These formulas are reminiscent of the Fourier expansions of Bernoulli polynomials. Some new results, for example, the one yielding the derivative of the Bessel function of the first kind with respect to its order as the Fourier coefficient of a function involving Chebyshev polynomials, are obtained in the course of proving these exact formulas. The asymptotic behavior of Zagier polynomials is also derived from them. Finally, a Zagier-type exact formula is obtained for $B_{2n}^{*}\left(-\frac{3}{2}\right)+B_{2n}^{*}$.


title: Generalized Majorization-Minimization

abstract: Non-convex optimization is ubiquitous in machine learning. Majorization-Minimization (MM) is a powerful iterative procedure for optimizing non-convex functions that works by optimizing a sequence of bounds on the function. In MM, the bound at each iteration is required to \emph{touch} the objective function at the optimizer of the previous bound. We show that this touching constraint is unnecessary and overly restrictive. We generalize MM by relaxing this constraint, and propose a new optimization framework, named Generalized Majorization-Minimization (G-MM), that is more flexible. For instance, G-MM can incorporate application-specific biases into the optimization procedure without changing the objective function. We derive G-MM algorithms for several latent variable models and show empirically that they consistently outperform their MM counterparts in optimizing non-convex objectives. In particular, G-MM algorithms appear to be less sensitive to initialization.


title: A profusion of $1/2$ BPS Wilson loops in $\mathcal{N}=4$ Chern-Simons-matter theories

abstract: We initiate the study of $1/2$ BPS Wilson loops in $\mathcal{N}=4$ Chern-Simons-matter theories in three dimensions. We consider a circular or linear quiver with Chern-Simons levels $k$, $-k$ and $0$, and focus on loops preserving one of the two $SU(2)$ subgroups of the $R$-symmetry. In the cases with no vanishing Chern-Simons levels, we find a pair of Wilson loops for each pair of adjacent nodes on the quiver connected by a hypermultiplet (nodes connected by twisted hypermultiplets have Wilson loops preserving another set of supercharges). We expect this classical pairwise degeneracy to be lifted by quantum corrections. In the case with nodes with vanishing Chern-Simons terms connected by twisted hypermultiplets, we find that the usual $1/4$ BPS Wilson loops are automatically enlarged to $1/2$ BPS, as happens also in 3-dimensional Yang-Mills theory. When the nodes with vanishing Chern-Simons levels are connected by untwisted hypermultiplets, we do not find any Wilson loops coupling to those nodes which are classically invariant. Rather, we find several loops whose supersymmetry variation, while non zero, vanishes in any correlation function, so is weakly zero. We expect only one linear combination of those Wilson loops to remain BPS when quantum corrections are included. We analyze the M-theory duals of those Wilson loops and comment on their degeneracy. We also show that these Wilson loops are cohomologically equivalent to certain $1/4$ BPS Wilson loops whose expectation value can be evaluated by the appropriate localized matrix model.