title: Defective fission correlation data from the 2E-2v method

abstract: The double-energy double-velocity ($2E$-$2v$) method allows assessing fission-fragment mass yields prior to and after prompt neutron emission with high resolution. It is, therefore, considered as a complementary technique to assess average prompt neutron multiplicity as a function of fragment properties. We have studied the intrinsic features of the $2E$-$2v$ method by means of event-wise generated fission-fragment data and found severe short-comings in the method itself as well as in some common practice of application. We find that the $2E$-$2v$ method leads to large deviations in the correlation between the prompt neutron multiplicity and pre-neutron mass, which deforms and exaggerates the so called `sawtooth' shape of $\bar{\nu}(A)$. We have identified the treatment of prompt neutron emission from the fragments as the origin of the problem. The intrinsic nature of this deficiency, risk to render $2E$-$2v$ experiments much less interesting. We suggest a method to correct the $2E$-$2v$ data, and recommend applying this method to previous data acquired in $2E$-$2v$ experiments, as well.


title: Force, Torque, Linear Momentum, and Angular Momentum in Classical Electrodynamics

abstract: The classical theory of electrodynamics is built upon Maxwell's equations and the concepts of electromagnetic (EM) field, force, energy, and momentum, which are intimately tied together by Poynting's theorem and by the Lorentz force law. Whereas Maxwell's equations relate the fields to their material sources, Poynting's theorem governs the flow of EM energy and its exchange between fields and material media, while the Lorentz law regulates the back-and-forth transfer of momentum between the media and the fields. An alternative force law, first proposed by Einstein and Laub, exists that is consistent with Maxwell's equations and complies with the conservation laws as well as with the requirements of special relativity. While the Lorentz law requires the introduction of hidden energy and hidden momentum in situations where an electric field acts on a magnetized medium, the Einstein-Laub (E-L) formulation of EM force and torque does not invoke hidden entities under such circumstances. Moreover, total force/torque exerted by EM fields on any given object turns out to be independent of whether the density of force/torque is evaluated using the law of Lorentz or that of Einstein and Laub. Hidden entities aside, the two formulations differ only in their predicted force and torque distributions inside matter. Such differences in distribution are occasionally measurable, and could serve as a guide in deciding which formulation, if either, corresponds to physical reality.


title: Analysis of Ultra High Energy Muons at the INO-ICAL Using Pair-Meter Technique

abstract: The proposed ICAL detector at INO is a large sized underground magnetized iron detector. ICAL is designed to reconstruct muon momentum using magnetic spectrometers. Energy measurement using magnets fail for muons in TeV range, since the angular deflection of the muon in the magnetic field is negligible and the muon tracks become nearly straight. A new technique for measuring the energy of muons in the TeV range is used by the CCFR neutrino detector, known as the Pair-Meter technique. This technique estimates muon energy from measurements of the energy deposited by the muon in many layers of an iron-calorimeter through e$^+$ and e$^-$ pair production. In this work we have performed Geant4 based preliminary analysis for iron plates and have demonstrated the observational feasibility of very high energy muons (1TeV-1000TeV) in a large mass underground detector operating as a pair-meter. This wide range of energy spectrum will be helpful for studying the cosmic rays in the Knee region and an understanding of the atmospheric neutrino flux for the present and future ultra high-energy atmospheric neutrino experiments.


title: Arbitrage and Geometry

abstract: This article introduces the notion of arbitrage for a situation involving a collection of investments and a payoff matrix describing the return to an investor of each investment under each of a set of possible scenarios. We explain the Arbitrage Theorem, discuss its geometric meaning, and show its equivalence to Farkas' Lemma. We then ask a seemingly innocent question: given a random payoff matrix, what is the probability of an arbitrage opportunity? This question leads to some interesting geometry involving hyperplane arrangements and related topics.


title: The magnetorotational and Tayler instabilities in the pulsar magnetosphere

abstract: The magnetospheres around neutron stars should be very particular because of their strong magnetic field and rapid rotation. A study of the pulsar magnetospheres is of a crucial importance since it is the key issue to understand how energy inflow to the exterior is produced. In this paper, I discuss magnetohydrodynamic processes in the pulsar magnetosphere. I consider in detail the properties of magnetohydrodynamic waves that can exist in the magnetosphere and their instabilities.These instabilities lead to formation of magnetic structures and can be responsible for a short-term variability ofthe pulsar emission.


title: A new look at distances and velocities of neutron stars

abstract: We take a fresh look at the determination of distances and velocities of neutron stars. The conversion of a parallax measurement into a distance, or distance probability distribution, has led to a debate quite similar to the one involving Cepheids, centering on the question whether priors can be used when discussing a single system. With the example of PSRJ0218+4232 we show that a prior is necessary to determine the probability distribution for the distance. The distance of this pulsar implies a gamma-ray luminosity larger than 10% of its spindown luminosity. For velocities the debate is whether a single Maxwellian describes the distribution for young pulsars. By limiting our discussion to accurate (VLBI) measurements we argue that a description with two Maxwellians, with distribution parameters sigma1=77 and sigma2=320 km/s, is significantly better. Corrections for galactic rotation, to derive velocities with respect to the local standards of rest, are insignificant.


title: The SPEDE spectrometer

abstract: The electron spectrometer, SPEDE, has been developed and will be employed in conjunction with the Miniball spectrometer at the HIE-ISOLDE facility, CERN. SPEDE allows for direct measurement of internal conversion electrons emitted in-flight, without employing magnetic fields to transport or momentum filter the electrons. Together with the Miniball spectrometer, it enables simultaneous observation of {\gamma} rays and conversion electrons in Coulomb-excitation experiments using radioactive ion beams.


title: ProCMotive: Bringing Programability and Connectivity into Isolated Vehicles

abstract: In recent years, numerous vehicular technologies, e.g., cruise control and steering assistant, have been proposed and deployed to improve the driving experience, passenger safety, and vehicle performance. Despite the existence of several novel vehicular applications in the literature, there still exists a significant gap between resources needed for a variety of vehicular (in particular, data-dominant, latency-sensitive, and computationally-heavy) applications and the capabilities of already-in-market vehicles. To address this gap, different smartphone-/Cloud-based approaches have been proposed that utilize the external computational/storage resources to enable new applications. However, their acceptance and application domain are still very limited due to programability, wireless connectivity, and performance limitations, along with several security/privacy concerns. In this paper, we present a novel architecture that can potentially enable rapid development of various vehicular applications while addressing shortcomings of smartphone-/Cloud-based approaches. The architecture is formed around a core component, called SmartCore, a privacy/security-friendly programmable dongle that brings general-purpose computational and storage resources to the vehicle and hosts in-vehicle applications. Based on the proposed architecture, we develop an application development framework for vehicles, that we call ProCMotive. ProCMotive enables developers to build customized vehicular applications along the Cloud-to-edge continuum, i.e., different functions of an application can be distributed across SmartCore, the user's personal devices, and the Cloud. To highlight potential benefits that the framework provides, we design and develop two different vehicular applications based on ProCMotive, namely, Amber Response and Insurance Monitor.


title: Vacuum dynamics in the Universe versus a rigid $\Lambda=$const

abstract: In this year, in which we celebrate 100 years of the cosmological term, $\Lambda$, in Einstein's gravitational field equations, we are still facing the crucial question whether $\Lambda$ is truly a fundamental constant or a mildly evolving dynamical variable. After many theoretical attempts to understand the meaning of $\Lambda$, and in view of the enhanced accuracy of the cosmological observations, it seems now mandatory that this issue should be first settled empirically before further theoretical speculations on its ultimate nature. In this work, we summarize the situation of some of these studies. Devoted analyses made recently show that the $\Lambda=$const. hypothesis, despite being the simplest, may well not be the most favored one. The overall fit to the cosmological observables $SNIa+BAO+H(z)+LSS+CMB$ singles out the class RVM of the "running" vacuum models, in which $\Lambda=\Lambda(H)$ is an affine power-law function of the Hubble rate. It turns out that the performance of the RVM as compared to the "concordance" $\Lambda$CDM model (with $\Lambda=$const.) is much better. The evidence in support of the RVM may reach $\sim 4\sigma$ c.l., and is bolstered with Akaike and Bayesian criteria providing strong evidence in favor of the RVM option. We also address the implications of this framework on the tension between the CMB and local measurements of the current Hubble parameter.


title: Low density, radiatively inefficient rotating-accretion flow onto a black hole

abstract: We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.


title: When Rational Sections Become Cyclic: Gauge Enhancement in F-theory via Mordell--Weil Torsion

abstract: We explore novel gauge enhancements from abelian to non-simply-connected gauge groups in F-theory. To this end we consider complex structure deformations of elliptic fibrations with a Mordell--Weil group of rank one and identify the conditions under which the generating section becomes torsional. For the specific case of Z2 torsion we construct the generic solution to these conditions and show that the associated F-theory compactification exhibits the global gauge group [SU(2) x SU(4)]/Z2 x SU(2). The subsolution with gauge group SU(2)/Z2 x SU(2), for which we provide a global resolution, is related by a further complex structure deformation to a genus-one fibration with a bisection whose Jacobian has a Z2 torsional section. While an analysis of the spectrum on the Jacobian fibration reveals an SU(2)/Z2 x Z2 gauge theory, reproducing this result from the bisection geometry raises some conceptual puzzles about F-theory on genus-one fibrations.


title: How the axial anomaly controls flavor mixing among mesons

abstract: It is well known that because of the axial anomaly in QCD, mesons with $J^P = 0^-$ are close to $SU(3)_{\mathrm{V}}$ eigenstates: the $eta^\prime(958)$ meson is largely a singlet, and the $\eta$ meson an octet. In contrast, states with $J^P = 1^-$ are flavor diagonal: \textit{e.g.}, the $\phi(1020)$ is almost pure $\bar{s} s$. Using effective Lagrangians, we show how this generalizes to states with higher spin, assuming that they can be classified according to the unbroken chiral symmetry of $G_{\mathrm{fl}} = SU(3)_{\mathrm{L}} \times SU(3)_{\mathrm{R}}$. We construct effective Lagrangians from terms invariant under $G_{\mathrm{fl}}$, and introduce the concept of \textit{hetero-} and \textit{homo}chiral multiplets. Because of the axial anomaly, only terms invariant under the $Z(3)_{\mathrm{A}}$ subgroup of the axial $U(1)_{\mathrm{A}}$ enter. For heterochiral multiplets, which begin with that including the $\eta$ and $\eta^\prime(958)$, there are $Z(3)_{\mathrm{A}}$ invariant terms with low mass dimension which cause states to mix according to $SU(3)_{\mathrm{V}}$ flavor. For homochiral multiplets, which begin with that including the $\phi(1020)$, there are no $Z(3)_{\mathrm{A}}$ invariant terms with low mass dimension, and so states are diagonal in flavor. In this way we predict the flavor mixing for the heterochiral multiplet with spin one, as well as for hetero- and homochiral multiplets with spin two and spin three.


title: Most pseudo-bulges can be formed at later stages of major mergers

abstract: Most giant spiral galaxies have pseudo or disk-like bulges that are considered to be the result of purely secular processes. This may challenge the hierarchical scenario predicting about one major merger per massive galaxy ($>$$3\times 10^{10} M_{\odot}$) since the last $\sim$ 9 billion years. Here we verify whether or not the association between pseudo-bulges and secular processes is irrevocable. Using GADGET2 N-body/SPH simulations, we have conducted a systematic study of remnants of major mergers which progenitors have been selected (1) to follow the gas richness-look back time relationship, and (2) with a representative distribution of orbits and spins in a cosmological frame. Analyzing the surface-mass density profile of both nearby galaxies and merger remnants with two components, we find that most of them show pseudo-bulges or bar dominated centers. Even if some orbits lead to classical bulges just after the fusion, the contamination by the additional gas that gradually accumulates to the center and forming stars later on, leads to remnants apparently dominated by pseudo-bulges. We also found that simple SPH simulations should be sufficient to form realistic spiral galaxies as remnants of ancient gas-rich mergers without need for specifically tuned feedback conditions. We then conclude that pseudo-bulges and bars in spiral galaxies are natural consequences of major mergers when they are realized in a cosmological context, i.e., with gas-rich progenitors as expected when selected in the distant Universe.


title: On symmetry inheritance of nonminimally coupled scalar fields

abstract: We present the first symmetry inheritance analysis of fields nonminimally coupled to gravity. In this work we are focused on the real scalar field $\phi$ with nonminimal coupling of the form $\xi\phi^2 R$. Possible cases of the symmetry noninheriting fields are constrained by the properties of the Ricci tensor and the scalar potential. Examples of such spacetimes can be found among those which are "dressed" with the stealth scalar field, a nontrivial scalar field configuration with the vanishing energy-momentum tensor. We classify the scalar field potentials which allow the symmetry noninheriting stealth field configurations on top of the exact solutions of the Einstein's gravitational field equation with the cosmological constant.


title: Internal dark matter structure of the most massive galaxy clusters

abstract: We investigate the evolution of the dark matter density profiles of the most massive galaxy clusters in the Universe. Using a `zoom-in' procedure on a large suite of cosmological simulations of total comoving volume of $3\,(h^{-1}\,\rm Gpc)^3$, we study the 25 most massive clusters in four redshift slices from $z\sim 1$ to the present. The minimum mass is $M_{500} > 5.5 \times 10^{14}$ M$_{\odot}$ at $z=1$. Each system has more than two million particles within $r_{500}$. Once scaled to the critical density at each redshift, the dark matter profiles within $r_{500}$ are strikingly similar from $z\sim1$ to the present day, exhibiting a low dispersion of 0.15 dex, and showing little evolution with redshift in the radial logarithmic slope and scatter. They have the running power law shape typical of the NFW-type profiles, and their inner structure, resolved to $3.8\,h^{-1}$ comoving kpc at $z=1$, shows no signs of converging to an asymptotic slope. Our results suggest that this type of profile is already in place at $z>1$ in the highest-mass haloes in the Universe, and that it remains exceptionally robust to merging activity.


title: Interferometric view of the circumstellar envelopes of northern FU Orionis-type stars

abstract: FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars.


title: Modeling the effects of inhomogeneous aerosols on the hot Jupiter Kepler-7b's atmospheric circulation

abstract: Motivated by the observational evidence of inhomogeneous clouds in exoplanetary atmospheres, we investigate how proposed simple cloud distributions can affect atmospheric circulations and infrared emission. We simulated temperatures and winds for the hot Jupiter Kepler-7b using a three-dimensional atmospheric circulation model that included a simplified aerosol radiative transfer model. We prescribed fixed cloud distributions and scattering properties based on results previously inferred from Kepler-7b optical phase curves, including inhomogeneous aerosols centered along the western terminator and hypothetical cases in which aerosols additionally extended across much of the planet's night side. In all cases, a strong jet capable of advecting aerosols from a cooler nightside to dayside was found to persist, but only at the equator. Colder temperatures at mid- and polar-latitudes might permit aerosol to form on the dayside without the need for advection. By altering the deposition and redistribution of heat, aerosols along the western terminator produced an asymmetric heating that effectively shifts the hottest spot further east of the sub-stellar point than expected for a uniform distribution. The addition of opaque high clouds on the nightside can partly mitigate this enhanced shift by retaining heat that contributes to warming west of the hotspot. These expected differences in infrared phase curves could place constraints on proposed cloud distributions and their infrared opacities for brighter hot Jupiters.


title: Renormalization of tensor networks using graph independent local truncations

abstract: We introduce an efficient algorithm for reducing bond dimensions in an arbitrary tensor network without changing its geometry. The method is based on a novel, quantitative understanding of local correlations in a network. Together with a tensor network coarse-graining algorithm, it yields a proper renormalization group (RG) flow. Compared to existing methods, the advantages of our algorithm are its low computational cost, simplicity of implementation, and applicability to any network. We benchmark it by evaluating physical observables for the 2D classical Ising model and find accuracy comparable with the best existing tensor network methods. Because of its graph independence, our algorithm is an excellent candidate for implementation of real-space RG in higher dimensions. We discuss some of the details and the remaining challenges in 3D. Source code for our algorithm is freely available.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: Dynamical quantum phase transitions: a review

abstract: Quantum theory provides an extensive framework for the description of the equilibrium properties of quantum matter. Yet experiments in quantum simulators have now opened up a route towards generating quantum states beyond this equilibrium paradigm. While these states promise to show properties not constrained by equilibrium principles such as the equal a priori probability of the microcanonical ensemble, identifying general properties of nonequilibrium quantum dynamics remains a major challenge especially in view of the lack of conventional concepts such as free energies. The theory of dynamical quantum phase transitions attempts to identify such general principles by lifting the concept of phase transitions to coherent quantum real-time evolution. This review provides a pedagogical introduction to this field. Starting from the general setting of nonequilibrium dynamics in closed quantum many-body systems, we give the definition of dynamical quantum phase transitions as phase transitions in time with physical quantities becoming nonanalytic at critical times. We summarize the achieved theoretical advances as well as the first experimental observations, and furthermore provide an outlook onto major open questions as well as future directions of research.


title: The Baryonic Collapse Efficiency of Galaxy Groups in the RESOLVE and ECO Surveys

abstract: We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at group-integrated cold baryonic mass M_coldbary ~ 10^11 Msun. The SAM, however, has significantly fewer groups at the transition mass ~ 10^11 Msun and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ~2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M_halo ~ 10^11.4-12 Msun, which we label "nascent groups." Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.


title: The BMS conjecture

abstract: I explain an open conjecture by Braverman/Milatovic/Shubin (BMS) on the positivity of square integrable solutions $f$ of $(-\Delta+1)f\geq 0$ on a geodescially complete Riemannian manifold, and its connection to essential self-adjointness problems of covariant Schr\"odinger operators. The latter conjecture has remained open for more than 14 years now.


title: Primordial perturbations with pre-inflationary bounce

abstract: Based on the effective field theory (EFT) of nonsingular cosmologies, we build a stable model, without the ghost and gradient instabilities, of bounce inflation (inflation is preceded by a cosmological bounce). We perform a full simulation for the evolution of scalar perturbation, and find that the perturbation spectrum has a large-scale suppression (as expected), which is consistent with the power deficit of the cosmic microwave background (CMB) TT-spectrum at low multipoles, but unexpectedly, it also shows itself one marked lower valley, which actually provides a better fit to the dip at multipole $l\sim 20$. The depth of valley is relevant with the physics around the bounce scale, which is model-dependent.


title: Cosmological implications of Primordial Black Holes

abstract: The possibility that a relevant fraction of the dark matter might be comprised of Primordial Black Holes (PBHs) has been seriously reconsidered after LIGO's detection of a $\sim 30 M_{\odot}$ binary black holes merger. Despite the strong interest in the model, there is a lack of studies on possible cosmological implications and effects on cosmological parameters inference. We investigate correlations with the other standard cosmological parameters using cosmic microwave background observations, finding significant degeneracies, especially with the tilt of the primordial power spectrum and the sound horizon at radiation drag. However, these degeneracies can be greatly reduced with the inclusion of small scale polarization data. We also explore if PBHs as dark matter in simple extensions of the standard $\Lambda$CDM cosmological model induces extra degeneracies, especially between the additional parameters and the PBH's ones. Finally, we present cosmic microwave background constraints on the fraction of dark matter in PBHs, not only for monochromatic PBH mass distributions but also for popular extended mass distributions. Our results show that extended mass distribution's constraints are tighter, but also that a considerable amount of constraining power comes from the high-$\ell$ polarization data. Moreover, we constrain the shape of such mass distributions in terms of the correspondent constraints on the PBH mass fraction.


title: The 6-girth-thickness of the complete graph

abstract: The $g$-girth-thickness $\theta(g,G)$ of a graph $G$ is the minimum number of planar subgraphs of girth at least $g$ whose union is $G$. In this paper, we determine the $6$-girth-thickness $\theta(6,K_n)$ of the complete graph $K_n$ in almost all cases. And also, we calculate by computer the missing value of $\theta(4,K_n)$.


title: Primordial Black Holes as Dark Matter: Converting Constraints from Monochromatic to Extended Mass Distributions

abstract: The model in which Primordial Black Holes (PBHs) constitute a non-negligible fraction of the dark matter has (re)gained popularity after the first detections of binary black hole mergers. Most of the observational constraints to date have been derived assuming a single mass for all the PBHs, although some more recent works tried to generalize constraints to the case of extended mass functions. Here we derive a general methodology to obtain constraints for any PBH Extended Mass Distribution (EMD) and any observables in the desired mass range. Starting from those obtained for a monochromatic distribution, we convert them into constraints for EMDs by using an equivalent, effective mass $M_{\rm eq}$ that depends on the specific observable. We highlight how limits of validity of the PBH modelling affect the EMD parameter space. Finally, we present converted constraints on the total abundance of PBH from microlensing, stellar distribution in ultra-faint dwarf galaxies and CMB accretion for Lognormal and Power Law mass distributions, finding that EMD constraints are generally stronger than monochromatic ones.


title: Uniform independence for Dehn twist automorphisms of a free group

abstract: McCarthy's Theorem for the mapping class group of a closed hyperbolic surface states that for any two mapping classes $\sigma,\tau \in \mathrm{Mod}(S)$ there is some power $N$ such that the group $\langle \sigma^N,\tau^N\rangle$ is either free of rank two or abelian, and gives a geometric criterion for the dichotomy. The analogous statement is false in linear groups, and unresolved for outer automorphisms of a free group. Several analogs are known for exponentially growing outer automorphisms satisfying various technical hypothesis. In this article we prove an analogous statement when $\sigma$ and $\tau$ are linearly growing outer automorphisms of $F_r$, and give a geometric criterion for the dichotomy. Further, Hamidi-Tehrani proved that for Dehn twists in the mapping class group this independence dichotomy is \emph{uniform}: $N=4$ suffices. In a similar style, we obtain an $N$ that depends only on the rank of the free group.


title: Massively parallel stochastic solution of the geophysical gravity problem

abstract: In this paper, we report the advantages of using a stochastic algorithm in the context of mineral exploration based on gravity measurements. This approach has the advantage over deterministic methods in that it allows one to find the solution of the Poisson equation in specified, isolated points without the need of meshing the computational domain and solving the Poisson equation over the entire domain. Moreover, the stochastic approach is embarrassingly parallelizable and therefore suitable for an implementation on multi-core compute clusters with or without GPUs. Benchmark tests are carried out that show that the stochastic approach can yield accurate results for both the gravitational potential and the gravitational acceleration and could hence provide an alternative to existing deterministic methods used in mineral exploration.


title: Learning Domain-Specific Word Embeddings from Sparse Cybersecurity Texts

abstract: Word embedding is a Natural Language Processing (NLP) technique that automatically maps words from a vocabulary to vectors of real numbers in an embedding space. It has been widely used in recent years to boost the performance of a vari-ety of NLP tasks such as Named Entity Recognition, Syntac-tic Parsing and Sentiment Analysis. Classic word embedding methods such as Word2Vec and GloVe work well when they are given a large text corpus. When the input texts are sparse as in many specialized domains (e.g., cybersecurity), these methods often fail to produce high-quality vectors. In this pa-per, we describe a novel method to train domain-specificword embeddings from sparse texts. In addition to domain texts, our method also leverages diverse types of domain knowledge such as domain vocabulary and semantic relations. Specifi-cally, we first propose a general framework to encode diverse types of domain knowledge as text annotations. Then we de-velop a novel Word Annotation Embedding (WAE) algorithm to incorporate diverse types of text annotations in word em-bedding. We have evaluated our method on two cybersecurity text corpora: a malware description corpus and a Common Vulnerability and Exposure (CVE) corpus. Our evaluation re-sults have demonstrated the effectiveness of our method in learning domain-specific word embeddings.


title: Stability of Spatial Smoothness and Cluster-Size Threshold Estimates in FMRI using AFNI

abstract: In a recent analysis of FMRI datasets [K Mueller et al, Front Hum Neurosci 11:345], the estimated spatial smoothness parameters and the statistical significance of clusters were found to depend strongly on the resampled voxel size (for the same data, over a range of 1 to 3 mm) in one popular FMRI analysis software package (SPM12). High sensitivity of thresholding results on such an arbitrary parameter as final spatial grid size is an undesirable feature in a processing pipeline. Here, we examine the stability of spatial smoothness and cluster-volume threshold estimates with respect to voxel resampling size in the AFNI software package's pipeline. A publicly available collection of resting-state and task FMRI datasets from 78 subjects was analyzed using standard processing steps in AFNI. We found that the spatial smoothness and cluster-volume thresholds are fairly stable over the voxel resampling size range of 1 to 3 mm, in contradistinction to the reported results from SPM12.


title: Unsupervised Contact Learning for Humanoid Estimation and Control

abstract: This work presents a method for contact state estimation using fuzzy clustering to learn contact probability for full, six-dimensional humanoid contacts. The data required for training is solely from proprioceptive sensors - endeffector contact wrench sensors and inertial measurement units (IMUs) - and the method is completely unsupervised. The resulting cluster means are used to efficiently compute the probability of contact in each of the six endeffector degrees of freedom (DoFs) independently. This clustering-based contact probability estimator is validated in a kinematics-based base state estimator in a simulation environment with realistic added sensor noise for locomotion over rough, low-friction terrain on which the robot is subject to foot slip and rotation. The proposed base state estimator which utilizes these six DoF contact probability estimates is shown to perform considerably better than that which determines kinematic contact constraints purely based on measured normal force.


title: On Two Theorems of Darboux

abstract: We provide precise formulations and proofs of two theorems from Darboux's lectures on orthogonal systems. These results provide local existence and uniqueness of solutions to certain types of first order PDE systems where each equation contains a single derivative for which it is solved: \[\frac{\partial u_i}{\partial x_j}(x)=f_{ij}(x,u(x)).\] The data prescribe values for the unknowns $u_i$ along certain hyperplanes through a given point $\bar x$. The first theorem applies to determined systems (the number of equations equals the number unknowns), and a unique, local solution is obtained via Picard iteration. While Darboux's statement of the theorem leaves unspecified "certaines conditions de continuit\'e," it is clear from his proof that he assumes Lipschitz continuity of the maps $f_{ij}$. On the other hand, he did not address the regularity of the data. We provide a precise formulation and proof of his first theorem. The second theorem is more involved and applies to overdetermined systems of the same general form. Under the appropriate integrability conditions, Darboux used his first theorem to treat the cases with two and three independent variables. We provide a proof for any number of independent variables. While the systems are rather special, they do appear in applications; e.g., the second theorem contains the classical Frobenius theorem on overdetermined systems as a special case. The key aspect of the proofs is that they apply to non-analytic situations. In an analytic setup the results are covered by the general Cartan-K\"ahler theorem.


title: Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation

abstract: We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time-evolution of the nuclear and electronic phase space variables as well as the Monodromy matrix, under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve-crossing in model two-level systems and show that in the quantum limit MQC-IVR is in good agreement with the exact quantum results, whereas in the classical limit the method yields results in keeping with mean-field approaches like the Linearized Semiclassical IVR. Finally, exploiting the ability of MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.


title: Computing bases of modular forms using the graded algebra structure

abstract: We develop a new algorithm to compute a basis for $M_k(\Gamma_0(N))$, the space of weight $k$ holomorphic modular forms on $\Gamma_0(N)$, in the case when the graded algebra of modular forms over $\Gamma_0(N)$ is generated at weight two. Our tests show that this algorithm significantly outperforms a commonly used algorithm which relies more heavily on modular symbols.


title: High-resolution Imaging of Transiting Extrasolar Planetary systems (HITEP). II. Lucky Imaging results from 2015 and 2016

abstract: The formation and dynamical history of hot Jupiters is currently debated, with wide stellar binaries having been suggested as a potential formation pathway. Additionally, contaminating light from both binary companions and unassociated stars can significantly bias the results of planet characterisation studies, but can be corrected for if the properties of the contaminating star are known. We search for binary companions to known transiting exoplanet host stars, in order to determine the multiplicity properties of hot Jupiter host stars. We also characterise unassociated stars along the line of sight, allowing photometric and spectroscopic observations of the planetary system to be corrected for contaminating light. We analyse lucky imaging observations of 97 Southern hemisphere exoplanet host stars, using the Two Colour Instrument on the Danish 1.54m telescope. For each detected companion star, we determine flux ratios relative to the planet host star in two passbands, and measure the relative position of the companion. The probability of each companion being physically associated was determined using our two-colour photometry. A catalogue of close companion stars is presented, including flux ratios, position measurements, and estimated companion star temperature. For companions that are potential binary companions, we review archival and catalogue data for further evidence. For WASP-77AB and WASP-85AB, we combine our data with historical measurements to determine the binary orbits, showing them to be moderately eccentric and inclined to the line of sight and planetary orbital axis. Combining our survey with the similar Friends of Hot Jupiters survey, we conclude that known hot Jupiter host stars show a deficit of high mass stellar companions compared to the field star population; however, this may be a result of the biases in detection and target selection by ground-based surveys.


title: On the observable spectrum of theories with a Brout-Englert-Higgs effect

abstract: The physical, observable spectrum in gauge theories is made up from gauge-invariant states. The Fr\"ohlich-Morchio-Strocchi mechanism allows in the standard model to map these states to the gauge-dependent elementary $W$, $Z$ and Higgs states. This is no longer necessarily the case in theories with a more general gauge group and Higgs sector. We classify and predict the physical spectrum for a wide range of such theories, with special emphasis on GUT-like cases, and show that discrepancies between the spectrum of elementary fields and physical particles frequently arise.


title: Particle-in-cell simulation of Buneman instability beyond quasilinear saturation

abstract: Spatio-temporal evolution of Buneman instability has been followed numerically till its quasilinear quenching and beyond, using an in-house developed electrostatic 1D particle-in-cell simulation code. For different initial drift velocities $k_{L}v_{0}/\omega_{pe} \approx 0.1 \, - \, 1$ and for a wide range of electron to ion mass ratios (m/M), growth rate obtained from simulation agrees well with the numerical solution of the fourth order dispersion relation. Quasi-linear saturation of Buneman instability occurs when ratio of electrostatic field energy density ($\sum\limits_{k} |E_{k}|^{2}/8{\pi}$) to initial electron drift kinetic energy density ($W_{0} = \frac{1}{2} n_{0}m v^{2}_{0}$) reaches up to a constant value, which as predicted by Hirose [Plasma Physics 20, 481(1978)], is independent of initial electron drift velocity but depends on electron to ion mass ratio m/M as $\sum\limits_{k} |E_{k}|^{2}/16{\pi}W_{0} \approx (m/M)^{1/3}$. This result stands verified in our simulations. Growth of the instability beyond the first saturation (quasilinear saturation ) till its final saturation [Ishihara et. al., PRL 44, 1404(1980)] follows an algebraic scaling with time. In contrast to the quasilinear saturation, the ratio of final saturated electrostatic field energy density to initial kinetic energy density, is relatively independent of electron to ion mass ratio and is found to depend only on the initial drift velocity. Beyond the final saturation, electron phase space holes coupled to large amplitude ion solitary waves, a state known as coupled hole-soliton , are seen in our simulations. The propagation characteristics ( amplitude - speed relation ) of these coherent modes is found to be consistent with the theory of Saeki et. al. [PRL 80, 1224(1998)].


title: Primordial non-Gaussianity and power asymmetry with quantum gravitational effects in loop quantum cosmology

abstract: Loop quantum cosmology (LQC) provides a resolution of the classical big bang singularity in the deep Planck era. The evolution, prior to the usual slow-roll inflation, naturally generates excited states at the onset of the slow-roll inflation. It is expected that these quantum gravitational effects could leave its fingerprints on the primordial perturbation spectrum and non-Gaussianity, and lead to some observational evidences in the cosmic microwave background (CMB). While the impact of the quantum effects on the primordial perturbation spectrum has been already studied and constrained by current data, in this paper we continue studying such effects on the non-Gaussianity of the primordial curvature perturbations. In this paper, we present detailed and analytical calculations of the non-Gaussianity and show explicitly that the corrections due to quantum effects are in the same magnitude of the slow-roll parameters in the observable scales and thus are well within current observational constraints. Despite this, we show that the non-Gaussianity in the squeezed limit can be enhanced at superhorizon scales and further, these effects may yield a large statistical anisotropy on the power spectrum through the Erickcek-Kamionkowski-Carroll mechanism.


title: Complexity of Scheduling Charging in the Smart Grid

abstract: In the smart grid, the intent is to use flexibility in demand, both to balance demand and supply as well as to resolve potential congestion. A first prominent example of such flexible demand is the charging of electric vehicles, which do not necessarily need to be charged as soon as they are plugged in. The problem of optimally scheduling the charging demand of electric vehicles within the constraints of the electricity infrastructure is called the charge scheduling problem. The models of the charging speed, horizon, and charging demand determine the computational complexity of the charge scheduling problem. For about 20 variants, we show, using a dynamic programming approach, that the problem is either in P or weakly NP-hard. We also show that about 10 variants of the problem are strongly NP-hard, presenting a potentially significant obstacle to their use in practical situations of scale.


title: First principles study of structural, magnetic and electronic properties of CrAs

abstract: We report ab initio density functional calculations of the structural and magnetic properties, and the electronic structure of CrAs. To simulate the observed pressure-driven experimental results, we perform our analysis for different volumes of the unit cell, showing that the structural, magnetic and electronic properties strongly depend on the size of the cell. We find that the calculated quantities are in good agreement with the experimental data, and we review our results in terms of the observed superconductivity.


title: Jacobians with prescribed eigenvectors

abstract: Let $\Omega\subset \mathbb{R}^n$ be open and let $\mathcal{R}$ be a partial frame on $\Omega$, that is a set of $m$ linearly independent vector fields prescribed on $\Omega$ ($m\leq n$). We consider the issue of describing the set of all maps $F:\Omega\to\mathbb{R}^n$ with the property that each of the given vector fields is an eigenvector of the Jacobian matrix of $F$. By introducing a coordinate independent definition of the Jacobian, we obtain an intrinsic formulation of the problem, which leads to an overdetermined PDE system, whose compatibility conditions can be expressed in an intrinsic, coordinate independent manner. To analyze this system we formulate and prove a generalization of the classical Frobenius integrability theorems. The size and structure of the solution set of this system depends on the properties of the partial frame, in particular, whether or not it is in involution. A particularly nice subclass of involutive partial frames, called rich, can be completely analyzed. Involutive, but non-rich case is somewhat harder to handle. We provide a complete answer in the case of $m=3$ and arbitrary $n$, as well as some general results for arbitrary $m$. The non-involutive case is far more challenging, and we only obtain a comprehensive analysis in the case $n=3$, $m=2$. Finally, we provide explicit examples illustrating the various possibilities. Our initial motivation for considering this problem comes from the geometric study of hyperbolic conservative systems in one spatial dimension.


title: Low Gilbert Damping Constant in Perpendicularly Magnetized W/CoFeB/MgO Films with High Thermal Stability

abstract: Perpendicular magnetic materials with low damping constant and high thermal stability have great potential for realizing high-density, non-volatile, and low-power consumption spintronic devices, which can sustain operation reliability for high processing temperatures. In this work, we study the Gilbert damping constant ({\alpha}) of perpendicularly magnetized W/CoFeB/MgO films with a high perpendicular magnetic anisotropy (PMA) and superb thermal stability. The {\alpha} of these PMA films annealed at different temperatures is determined via an all-optical Time-Resolved Magneto-Optical Kerr Effect method. We find that {\alpha} of these W/CoFeB/MgO PMA films decreases with increasing annealing temperature, reaches a minimum of {\alpha} = 0.016 at an annealing temperature of 350 {\deg}C, and then increases to 0.024 after post-annealing at 400 {\deg}C. The minimum {\alpha} observed at 350 {\deg}C is rationalized by two competing effects as the annealing temperature becomes higher: the enhanced crystallization of CoFeB and dead-layer growth occurring at the two interfaces of the CoFeB layer. We further demonstrate that {\alpha} of the 400 {\deg}C-annealed W/CoFeB/MgO film is comparable to that of a reference Ta/CoFeB/MgO PMA film annealed at 300 {\deg}C, justifying the enhanced thermal stability of the W-seeded CoFeB films.


title: WERd: Using Social Text Spelling Variants for Evaluating Dialectal Speech Recognition

abstract: We study the problem of evaluating automatic speech recognition (ASR) systems that target dialectal speech input. A major challenge in this case is that the orthography of dialects is typically not standardized. From an ASR evaluation perspective, this means that there is no clear gold standard for the expected output, and several possible outputs could be considered correct according to different human annotators, which makes standard word error rate (WER) inadequate as an evaluation metric. Such a situation is typical for machine translation (MT), and thus we borrow ideas from an MT evaluation metric, namely TERp, an extension of translation error rate which is closely-related to WER. In particular, in the process of comparing a hypothesis to a reference, we make use of spelling variants for words and phrases, which we mine from Twitter in an unsupervised fashion. Our experiments with evaluating ASR output for Egyptian Arabic, and further manual analysis, show that the resulting WERd (i.e., WER for dialects) metric, a variant of TERp, is more adequate than WER for evaluating dialectal ASR.


title: The Covering Path Problem on a Grid

abstract: This paper introduces the covering path problem on a grid (CPPG) which finds the cost-minimizing path connecting a subset of points in a grid such that each point that needs to be covered is within a predetermined distance of a point from the chosen subset. We leverage the geometric properties of the grid graph which captures the road network structure in many transportation problems, including our motivating setting of school bus routing. As defined in this paper, the CPPG is a bi-objective optimization problem comprised of one cost term related to path length and one cost term related to stop count. We develop a trade-off constraint which quantifies the trade-off between path length and stop count and provides a lower bound for the bi-objective optimization problem. We introduce simple construction techniques to provide feasible paths that match the lower bound within a constant factor. Importantly, this solution approach uses transformations of the general CPPG to either a discrete CPPG or continuous CPPG based on the value of the coverage radius. For both the discrete and continuous versions, we provide fast constant-factor approximations, thus solving the general CPPG.


title: Rogue Wave Morphology in Broadband Nonbreaking Seas

abstract: Here, we show through a statistical analysis that rogue waves in broadband non-breaking seas are spatially asymmetric. In addition to the top-down asymmetry due to nonlinear effects, we show that the two troughs adjacent to the rogue wave crest are generally of different depths, which is unlike the conventional picture of rogue waves with symmetric fore and aft troughs often obtained from model equations. The rogue-wave trough asymmetry is demonstrated for sea states 4 to 6 on Douglas Sea Scale. Considering the deepest trough leads to approximately 10% increase in the calculation of the mean rogue-wave height compared to previous results for rogue waves reported with symmetric troughs. This implies that estimates of rogue-wave trough-crest amplitudes based on model equations should be re-assessed upward for all realistic oceanic conditions.


title: Computing the Unique Information

abstract: Given a pair of predictor variables and a response variable, how much information do the predictors have about the response, and how is this information distributed between unique, redundant, and synergistic components? Recent work has proposed to quantify the unique component of the decomposition as the minimum value of the conditional mutual information over a constrained set of information channels. We present an efficient iterative divergence minimization algorithm to solve this optimization problem with convergence guarantees and evaluate its performance against other techniques.


title: Dynamics of sedimenting active Brownian particles

abstract: We investigate the stochastic dynamics of one sedimenting active Brownian particle in three dimensions under the influence of gravity and passive fluctuations in the translational and rotational motion. We present an analytical solution of the Fokker--Planck equation for the stochastic process which allows us to describe the dynamics of the active Brownian particle in three dimensions. We address the time evolution of the monopole, the polarization, and the steady-state solution. We also perform Brownian dynamics simulations and study the effect of the activity of the particles on their collective motion. These results qualitatively agree with our model. Finally, we compare our results with experiments [J. Palacci \emph{et al.}, Phys. Rev. Lett. \textbf{105}, 088304 (2010)] and find very good agreement.


title: A Radio Relic and a Search for the Central Black Hole in the Abell 2261 Brightest Cluster Galaxy

abstract: We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.


title: Combining the Ensemble and Franck-Condon Approaches for Spectral Shapes of Molecules in Solution

abstract: The correct treatment of vibronic effects is vital for the modeling of absorption spectra of solvated dyes, as many prominent spectral features can often be ascribed to vibronic transitions. Vibronic spectra can be computed within the Franck-Condon approximation for small dyes in solution using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions and provide only an approximate treatment of solvent polarization effects. Furthermore, temperature-dependent solvent broadening effects are often only accounted for using a broadening parameter chosen to match experimental spectra. On the other hand, ensemble approaches provide a straightforward way of accounting for solute-solvent interactions and temperature-dependent broadening effects by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. However, ensemble approaches do not explicitly account for vibronic effects and often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining the vibronic fine structure of an excitation obtained in the Franck-Condon picture at zero temperature with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is therefore treated classically through the sampling of configurations, with vibronic contributions included as a zero-temperature correction to each vertical transition. We test the proposed method on Nile Red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interaction, the approach yields a significant improvement over the ensemble approach, whereas for systems with weaker interactions, both the shape and the width of the spectra are in excellent agreement with experiment.


title: Density operator approach for Landau problem quantum Hamiltonians

abstract: In this work, the definition of the density operator on quantum states in Hilbert spaces and some of its aspects relevant in thermodynamics and information-theoretical entropy calculations are given. In this framework, a physical model describing an electron in a magnetic field is investigated. The so-called exotic Landau problem in noncommutative plane is also considered. Then, a model related to the fractional quantum Hall effect is revisited. Thanks to the completeness relations verified by the coherent states (CS) in these models, the thermodynamics is discussed by using the diagonal $P$-representation of thedensity operator. Specifically, the Q-Husimi distribution and the Wehrl entropy are determined.


title: Sparse-to-Dense: Depth Prediction from Sparse Depth Samples and a Single Image

abstract: We consider the problem of dense depth prediction from a sparse set of depth measurements and a single RGB image. Since depth estimation from monocular images alone is inherently ambiguous and unreliable, to attain a higher level of robustness and accuracy, we introduce additional sparse depth samples, which are either acquired with a low-resolution depth sensor or computed via visual Simultaneous Localization and Mapping (SLAM) algorithms. We propose the use of a single deep regression network to learn directly from the RGB-D raw data, and explore the impact of number of depth samples on prediction accuracy. Our experiments show that, compared to using only RGB images, the addition of 100 spatially random depth samples reduces the prediction root-mean-square error by 50% on the NYU-Depth-v2 indoor dataset. It also boosts the percentage of reliable prediction from 59% to 92% on the KITTI dataset. We demonstrate two applications of the proposed algorithm: a plug-in module in SLAM to convert sparse maps to dense maps, and super-resolution for LiDARs. Software and video demonstration are publicly available.


title: Big Data Systems Meet Machine Learning Challenges: Towards Big Data Science as a Service

abstract: Recently, we have been witnessing huge advancements in the scale of data we routinely generate and collect in pretty much everything we do, as well as our ability to exploit modern technologies to process, analyze and understand this data. The intersection of these trends is what is called, nowadays, as Big Data Science. Cloud computing represents a practical and cost-effective solution for supporting Big Data storage, processing and for sophisticated analytics applications. We analyze in details the building blocks of the software stack for supporting big data science as a commodity service for data scientists. We provide various insights about the latest ongoing developments and open challenges in this domain.


title: Sullivan constructions for transitive Lie algebroids - smooth case

abstract: Let $M$ be a smooth manifold, smoothly triangulated by a simplicial complex $K$, and $\cA$ a transitive Lie algebroid on $M$. The Lie algebroid restriction of $\cA$ to a simplex $\Delta$ of $K$ is denoted by $\cA^{!!}_{\Delta}$. A piecewise smooth form of degree $p$ on $\cA$ is a family $\omega=(\omega_{\Delta})_{\Delta\in K}$ such that $\omega_{\Delta}\in \Omega^{p}(\cA^{!!}_{\Delta};\Delta)$ for each $\Delta\in K$, satisfying the compatibility condition concerning the restrictions of $\omega_{\Delta}$ to the faces of $\Delta$, that is, if $\Delta'$ is a face of $\Delta$, the restriction of the form $\omega_{\Delta}$ to the simplex $\Delta'$ coincides with the form $\omega_{\Delta'}$. The set $\Omega^{\ast}(\cA;K)$ of all piecewise smooth forms on $\cA$ is a cochain algebra. One has a natural morphism $$\Omega^{\ast}(\cA;M)\rightarrow \Omega^{\ast}(\cA;K)$$ of cochain algebras given by restriction of a smooth form defined on $\cA$ to a smooth form defined on $\cA^{!!}_{\Delta}$, for all simplices $\Delta$ of $K$. In this paper, we prove that, for triangulated compact manifolds, the cohomology of this construction is isomorphic to the Lie algebroid cohomology of $\cA$, in which the isomorphism is induced by the restriction map.


title: A Symbolic Approach to Safety LTL Synthesis

abstract: Temporal synthesis is the automated design of a system that interacts with an environment, using the declarative specification of the system's behavior. A popular language for providing such a specification is Linear Temporal Logic, or LTL. LTL synthesis in the general case has remained, however, a hard problem to solve in practice. Because of this, many works have focused on developing synthesis procedures for specific fragments of LTL, with an easier synthesis problem. In this work, we focus on Safety LTL, defined here to be the Until-free fragment of LTL in Negation Normal Form~(NNF), and shown to express a fragment of safe LTL formulas. The intrinsic motivation for this fragment is the observation that in many cases it is not enough to say that something "good" will eventually happen, we need to say by when it will happen. We show here that Safety LTL synthesis is significantly simpler algorithmically than LTL synthesis. We exploit this simplicity in two ways, first by describing an explicit approach based on a reduction to Horn-SAT, which can be solved in linear time in the size of the game graph, and then through an efficient symbolic construction, allowing a BDD-based symbolic approach which significantly outperforms extant LTL-synthesis tools.


title: Factorization method and general second order linear difference equation

abstract: This paper addresses an investigation on a factorization method for difference equations. It is proved that some classes of second order linear difference operators, acting in Hilbert spaces, can be factorized using a pair of mutually adjoint first order difference operators. These classes encompass equations of hypergeometic type describing classical orthogonal polynomials of a discrete variable.


title: Evidence for the Higgs boson decay to a bottom quark-antiquark pair

abstract: A search for the standard model (SM) Higgs boson ($\mathrm{H}$) decaying to $\mathrm{b}\overline{\mathrm{b}}$ when produced in association with an electroweak vector boson is reported for the following processes: $\mathrm{Z}(\nu\nu)\mathrm{H}$, $\mathrm{W}(\mu \nu)\mathrm{H}$, $\mathrm{W}(\mathrm{e} \nu)\mathrm{H}$, $\mathrm{Z}(\mu\mu)\mathrm{H}$, and $\mathrm{Z}(\mathrm{e}\mathrm{e})\mathrm{H}$. The search is performed in data samples corresponding to an integrated luminosity of 35.9 fb$^{-1}$ at $\sqrt{s}=$ 13 TeV recorded by the CMS experiment at the LHC during Run 2 in 2016. An excess of events is observed in data compared to the expectation in the absence of a $\mathrm{H} \to\mathrm{b}\overline{\mathrm{b}}$ signal. The significance of this excess is 3.3 standard deviations, where the expectation from SM Higgs boson production is 2.8. The signal strength corresponding to this excess, relative to that of the SM Higgs boson production, is 1.2 $\pm$ 0.4. When combined with the Run 1 measurement of the same processes, the signal significance is 3.8 standard deviations with 3.8 expected. The corresponding signal strength, relative to that of the SM Higgs boson, is 1.06 $^{+0.31}_{-0.29}$.


title: Decision making and uncertainty quantification for individualized treatments

abstract: Individualized treatment rules (ITR) can improve health outcomes by recognizing that patients may respond differently to treatment and assigning therapy with the most desirable predicted outcome for each individual. Flexible and efficient prediction models are desired as a basis for such ITRs to handle potentially complex interactions between patient factors and treatment. Modern Bayesian semiparametric and nonparametric regression models provide an attractive avenue in this regard as these allow natural posterior uncertainty quantification of patient specific treatment decisions as well as the population wide value of the prediction-based ITR. In addition, via the use of such models, inference is also available for the value of the Optimal ITR. We propose such an approach and implement it using Bayesian Additive Regression Trees (BART) as this model has been shown to perform well in fitting nonparametric regression functions to continuous and binary responses, even with many covariates. It is also computationally efficient for use in practice. With BART we investigate a treatment strategy which utilizes individualized predictions of patient outcomes from BART models. Posterior distributions of patient outcomes under each treatment are used to assign the treatment that maximizes the expected posterior utility. We also describe how to approximate such a treatment policy with a clinically interpretable ITR, and quantify its expected outcome. The proposed method performs very well in extensive simulation studies in comparison with several existing methods. We illustrate the usage of the proposed method to identify an individualized choice of conditioning regimen for patients undergoing hematopoietic cell transplantation and quantify the value of this method of choice in relation to the Optimal ITR as well as non-individualized treatment strategies.


title: Surfactant chain length and concentration influence on the interfacial tension of two immiscible model liquids a coarse grained approach

abstract: The interfacial tension between immiscible liquids is studied as a function of a model linear surfactant length and concentration using coarse grained, dissipative particle dynamics numerical simulations. The adsorption isotherms obtained from the simulations are found to be in agreement with Langmuir model. The reduction of the interfacial tension with increasing surfactant concentration is found to display some common characteristics for all the values of chain length modeled, with our predictions being in agreement with Szyszkowski equation. Lastly, the critical micelle concentration is predicted for all surfactant lengths, finding exponentially decaying behavior, in agreement with Kleven model. It is argued that these findings can be helpful guiding tools in the interpretation of available experiments and in the design of new ones with new surfactants and polymers.


title: Charge transfer in iridate-manganite superlattices

abstract: Charge transfer in superlattices consisting of SrIrO$_3$ and SrMnO$_3$ is investigated using density functional theory. Despite the nearly identical work function and non-polar interfaces between SrIrO$_3$ and SrMnO$_3$, rather large charge transfer was experimentally reported at the interface between them. Here, we report a microscopic model that captures the mechanism behind this phenomenon, providing a qualitative understanding of the experimental observation. This leads to unique strain dependence of such charge transfer in iridate-manganite superlattices. The predicted behavior is consistently verified by experiment with soft x-ray and optical spectroscopy. Our work thus demonstrates a new route to control electronic states in non-polar oxide heterostructures.


title: Reorientation-effect measurement of the first 2$^+$ state in $^{12}$C: confirmation of oblate deformation

abstract: A Coulomb-excitation reorientation-effect measurement using the TIGRESS $\gamma-$ray spectrometer at the TRIUMF/ISAC II facility has permitted the first determination of the $\langle 2^+_1\mid\mid \hat{E2} \mid\mid 2^+_1\rangle$ diagonal matrix element in $^{12}$C from particle$-\gamma$ coincidence data. Required state-of-the-art no-core shell model calculations of the nuclear polarizability for the ground and first-excited (2$^+_1$) states in $^{12}$C using chiral NN N$^4$LO500 and NN+3NF350 interactions have been performed. Consistent predictions show a larger polarizability than previously anticipated. The polarizability of the 2$^+_1$ state is introduced into the current and previous Coulomb-excitation reorientation-effect analysis of $^{12}$C. Spectroscopic quadrupole moments of $Q_{_S}(2_1^+)= +0.053(44)$ eb and $Q_{_S}(2_1^+)= +0.08(3)$ eb are determined, respectively, yielding a weighted average of $Q_{_S}(2_1^+)= +0.071(25)$ eb, in agreement with recent ab initio calculations. The present measurement confirms that the 2$^+_1$ state of $^{12}$C is oblate and emphasizes the important role played by the nuclear polarizability in Coulomb-excitation studies of light nuclei.


title: A Multimodal, Full-Surround Vehicular Testbed for Naturalistic Studies and Benchmarking: Design, Calibration and Deployment

abstract: Recent progress in autonomous and semi-autonomous driving has been made possible in part through an assortment of sensors that provide the intelligent agent with an enhanced perception of its surroundings. It has been clear for quite some while now that for intelligent vehicles to function effectively in all situations and conditions, a fusion of different sensor technologies is essential. Consequently, the availability of synchronized multi-sensory data streams are necessary to promote the development of fusion based algorithms for low, mid and high level semantic tasks. In this paper, we provide a comprehensive description of LISA-A: our heavily sensorized, full-surround testbed capable of providing high quality data from a slew of synchronized and calibrated sensors such as cameras, LIDARs, radars, and the IMU/GPS. The vehicle has recorded over 100 hours of real world data for a very diverse set of weather, traffic and daylight conditions. All captured data is accurately calibrated and synchronized using timestamps, and stored safely in high performance servers mounted inside the vehicle itself. Details on the testbed instrumentation, sensor layout, sensor outputs, calibration and synchronization are described in this paper.


title: Can the graviton have a large mass near black holes?

abstract: The mass of the graviton, if nonzero, is usually considered to be very small, e.g. of the Hubble scale, from several observational constraints. In this paper, we propose a gravity model where the graviton mass is very small in the usual weak gravity environments, below all the current graviton mass bounds, but becomes much larger in the strong gravity regime such as a black hole's vicinity. For black holes in this model, significant deviations from general relativity emerge very close to the black hole horizon and alter the black hole quasi-normal modes, which can be extracted from the ringdown waveform of black hole binary mergers. Also, the enhancement of the graviton mass near the horizon can result in echoes in the late time ringdown, which can be verified in the upcoming gravitational wave observations of higher sensitivity.


title: Hopf monoids and generalized permutahedra

abstract: Generalized permutahedra are a family of polytopes with a rich combinatorial structure and strong connections to optimization. We prove that they are the universal family of polyhedra with a certain Hopf algebraic structure. Their antipode is remarkably simple: the antipode of a polytope is the alternating sum of its faces. Our construction provides a unifying framework to organize numerous combinatorial structures, including graphs, matroids, posets, set partitions, linear graphs, hypergraphs, simplicial complexes, building sets, and simple graphs. We highlight three applications: 1. We obtain uniform proofs of numerous old and new results about the Hopf algebraic and combinatorial structures of these families. In particular, we give the optimal formula for the antipode of graphs, posets, matroids, hypergraphs, and building sets, and we answer questions of Humpert--Martin and Rota. 2. We show that the reciprocity theorems of Stanley and Billera--Jia--Reiner on chromatic polynomials of graphs, order polynomials of posets, and BJR-polynomials of matroids are instances of the same reciprocity theorem for generalized permutahedra. 3. We explain why the formulas for the multiplicative and compositional inverses of power series are governed by the face structure of permutahedra and associahedra, respectively, answering a question of Loday. Along the way, we offer a combinatorial user's guide to Hopf monoids.


title: Electrically controllable spin filtering based on superconducting helical states

abstract: The magnetoelectric effects in the surface states of the 3D TI are extremely strong due to the full spin-momentum locking. Here the microscopic theory of S/3D TI bilayer structures in terms of quasiclassical Green's functions is developed. On the basis of the developed formalism it is shown that the DOS in the S/TI bilayer manifests giant magnetoelectric behavior and, as a result, S/3D TI heterostructures can work as non-magnetic fully electrically controllable spin filters. It is shown that due to the full spin-momentum locking the amplitudes of the odd-frequency singlet and triplet components of the consensate wave function are equal. The same is valid for the even frequency singlet and triplet components. We unveil the connection between the odd-frequency pairing in S/3D TI heterostructures and magnetoelectric effects in the DOS.


title: An Empirical Dynamic Programming Algorithm for Continuous MDPs

abstract: We propose universal randomized function approximation-based empirical value iteration (EVI) algorithms for Markov decision processes. The `empirical' nature comes from each iteration being done empirically from samples available from simulations of the next state. This makes the Bellman operator a random operator. A parametric and a non-parametric method for function approximation using a parametric function space and the Reproducing Kernel Hilbert Space (RKHS) respectively are then combined with EVI. Both function spaces have the universal function approximation property. Basis functions are picked randomly. Convergence analysis is done using a random operator framework with techniques from the theory of stochastic dominance. Finite time sample complexity bounds are derived for both universal approximate dynamic programming algorithms. Numerical experiments support the versatility and effectiveness of this approach.


title: Dynamically Correlated Region in Sheared Colloidal Glasses Revealed by Neutron Scattering

abstract: The microscopic deformation mechanism of charged colloidal glasses with extended-range interactions under shear is investigated by in-situ small-angle neutron scattering, and a dynamically correlated region (DCR) is identified. This short-lived region provides the resistance to the configurational rearrangement imposed by the external deformation, as evidenced by the evolution of the size of DCR in the shear thinning regime and the quantitative agreement between the local stress sustained by DCR and the macroscopic stress from rheological measurements at low and mediate shear rates. This finding suggests that DCR is an important quantity for microscopically addressing the flow and deformation behavior of strongly interacting colloids.


title: Hijacking .NET to Defend PowerShell

abstract: With the rise of attacks using PowerShell in the recent months, there has not been a comprehensive solution for monitoring or prevention. Microsoft recently released the AMSI solution for PowerShell v5, however this can also be bypassed. This paper focuses on repurposing various stealthy runtime .NET hijacking techniques implemented for PowerShell attacks for defensive monitoring of PowerShell. It begins with a brief introduction to .NET and PowerShell, followed by a deeper explanation of various attacker techniques, which is explained from the perspective of the defender, including assembly modification, class and method injection, compiler profiling, and C based function hooking. Of the four attacker techniques that are repurposed for defensive real-time monitoring of PowerShell execution, intermediate language binary modification, JIT hooking, and machine code manipulation provide the best results for stealthy run-time interfaces for PowerShell scripting analysis.


title: Maximum oriented forcing number for complete graphs

abstract: The maximum oriented $k$-forcing number of a simple graph $G$, written $\MOF_k(G)$, is the maximum directed $k$-forcing number among all orientations of $G$. This invariant was recently introduced by Caro, Davila and Pepper in [CaroDavilaPepper], and in the current paper we study the special case where $G$ is the complete graph with order $n$, denoted $K_n$. While $\MOF_k(G)$ is an invariant for the underlying simple graph $G$, $\MOF_k(K_n)$ can also be interpreted as an interesting property for tournaments. Our main results further focus on the case when $k=1$. These include a lower bound on $\MOF(K_n)$ of roughly $\frac{3}{4}n$, and for $n\ge 2$, a lower bound of $n - \frac{2n}{\log_2(n)}$. Along the way, we also consider various lower bounds on the maximum oriented $k$-forcing number for the closely related complete $q$-partite graphs.


title: Electronic Structure Theory of Strained Two-Dimensional Materials with Hexagonal Symmetry

abstract: We derive electronic tight-binding Hamiltonians for strained graphene, hexagonal boron nitride and transition metal dichalcogenides based on Wannier transformation of {\it ab initio} density functional theory calculations. Our microscopic models include strain effects to leading order that respect the hexagonal crystal symmetry and local crystal configuration, and are beyond the central force approximation which assumes only pair-wise distance dependence. Based on these models, we also derive and analyze the effective low-energy Hamiltonians. Our {\it ab initio} approaches complement the symmetry group representation construction for such effective low-energy Hamiltonians and provide the values of the coefficients for each symmetry-allowed term. These models are relevant for the design of electronic device applications, since they provide the framework for describing the coupling of electrons to other degrees of freedom including phonons, spin and the electromagnetic field. The models can also serve as the basis for exploring the physics of many-body systems of interesting quantum phases.


title: Robust Optimization of Unconstrained Binary Quadratic Problems

abstract: In this paper we focus on the unconstrained binary quadratic optimization model, maximize x^t Qx, x binary, and consider the problem of identifying optimal solutions that are robust with respect to perturbations in the Q matrix.. We are motivated to find robust, or stable, solutions because of the uncertainty inherent in the big data origins of Q and limitations in computer numerical precision, particularly in a new class of quantum annealing computers. Experimental design techniques are used to generate a diverse subset of possible scenarios, from which robust solutions are identified. An illustrative example with practical application to business decision making is examined. The approach presented also generates a surface response equation which is used to estimate upper bounds in constant time for Q instantiations within the scenario extremes. In addition, a theoretical framework for the robustness of individual x_i variables is considered by examining the range of Q values over which the x_i are predetermined.


title: Symmetry and structure of carbon-nitrogen complexes in gallium arsenide from infrared spectroscopy and first-principles calculations

abstract: Molecular-like carbon-nitrogen complexes in GaAs are investigated both experimentally and theoretically. Two characteristic high-frequency stretching modes at \num{1973} and \SI{2060}{cm^{-1}}, detected by Fourier transform infrared absorption (FTIR) spectroscopy, appear in carbon- and nitrogen-implanted and annealed layers. From isotopic substitution it is deduced that the chemical composition of the underlying complexes is CN$_2$ and C$_2$N, respectively. Piezospectroscopic FTIR measurements reveal that both centers have tetragonal symmetry. For density functional theory (DFT) calculations linear entities are substituted for the As anion, with the axis oriented along the \hkl<100> direction, in accordance with the experimentally ascertained symmetry. The DFT calculations support the stability of linear N-C-N and C-C-N complexes in the GaAs host crystal in the charge states ranging from $+3$ to $-3$. The valence bonds of the complexes are analyzed using molecular-like orbitals from DFT. It turns out that internal bonds and bonds to the lattice are essentially independent of the charge state. The calculated vibrational mode frequencies are close to the experimental values and reproduce precisely the isotopic mass splitting from FTIR experiments. Finally, the formation energies show that under thermodynamic equilibrium CN$_2$ is more stable than C$_2$N.


title: Progress in computing parton distribution functions from the quasi-PDF approach

abstract: We discuss the current developments by the European Twisted Mass Collaboration in extracting parton distribution functions from the quasi-PDF approach. We concentrate on the non-perturbative renormalization prescription recently developed by us, using the RI$'$ scheme. We show results for the renormalization functions of matrix elements needed for the computation of quasi-PDFs, including the conversion to the $\overline{\rm MS}$ scheme, and for renormalized matrix elements. We discuss the systematic effects present in the $Z$-factors and the possible ways of addressing them in the future.


title: Critical random forests

abstract: Let $F(N,m)$ denote a random forest on a set of $N$ vertices, chosen uniformly from all forests with $m$ edges. Let $F(N,p)$ denote the forest obtained by conditioning the Erdos-Renyi graph $G(N,p)$ to be acyclic. We describe scaling limits for the largest components of $F(N,p)$ and $F(N,m)$, in the critical window $p=N^{-1}+O(N^{-4/3})$ or $m=N/2+O(N^{2/3})$. Aldous described a scaling limit for the largest components of $G(N,p)$ within the critical window in terms of the excursion lengths of a reflected Brownian motion with time-dependent drift. Our scaling limit for critical random forests is of a similar nature, but now based on a reflected diffusion whose drift depends on space as well as on time.


title: Eco-evolutionary dynamics and collective dispersal: implications for salmon metapopulation robustness

abstract: The spatial dispersal of individuals is known to play an important role in the dynamics of populations, and is central to metapopulation theory. At the same time, local adaptation to environmental conditions creates a geographic mosaic of evolutionary forces, where the combined drivers of selection and gene flow interact. Although the dispersal of individuals from donor to recipient populations provides connections within the metapopulation, promoting demographic and evolutionary rescue, it may also introduce maladapted individuals into habitats host to different environmental conditions, potentially lowering the fitness of the recipient population. Here we explore a model of the eco-evolutionary dynamics between two populations connected by dispersal, where the productivity of each is defined by a trait complex that is subject to local selection. Although general in nature, our model is inspired by salmon metapopulations, where dispersal between populations is defined in terms of the straying rate, which has been shown to be density-dependent. The results of our model reveal that increased straying between evolving populations leads to alternative stable states, which has large and nonlinear effects on two measures of metapopulation robustness: the portfolio effect and the time to recovery following an induced disturbance. We show that intermediate levels of straying result in large gains in robustness, and that increased habitat heterogeneity promotes robustness when straying rates are low, and erodes robustness when straying rates are high. Finally, we show that density-dependent straying promotes robustness, particularly when the aggregate biomass is low and straying is correspondingly high, which has important ramifications for the conservation of salmon metapopulations facing both natural and anthropogenic disturbances.


title: Characterization of the beam from the RFQ of the PIP-II Injector Test

abstract: A 2.1 MeV, 10 mA CW RFQ has been installed and commissioned at Fermilab's test accelerator known as PIP-II Injector Test. This report describes the measurements of the beam properties after acceleration in the RFQ, including the energy and emittance.


title: Blowup rate for mass critical rotational nonlinear Schr\"odinger equations

abstract: We consider the blowup rate for blowup solutions to $L^2$-critical, focusing NLS with a harmonic potential and a rotation term. Under a suitable spectral condition we prove that there holds the "$\log$-$\log$ law" when the initial data is slightly above the ground state. We also construct minimal mass blowup solutions near the ground state level with distinct blowup rates.


title: Bolt Detection Signal Analysis Method Based on ICEEMD

abstract: The construction quality of the bolt is directly related to the safety of the project, and as such, it must be tested. In this paper, the improved complete ensemble empirical mode decomposition (ICEEMD) method is introduced to the bolt detection signal analysis. The ICEEMD is used in order to decompose the anchor detection signal according to the approximate entropy of each intrinsic mode function (IMF). The noise of the IMFs is eliminated by the wavelet soft threshold de-noising technique. Based on the approximate entropy, and the wavelet de-noising principle, the ICEEMD-De anchor signal analysis method is proposed. From the analysis of the vibration analog signal, as well as the bolt detection signal, the result shows that the ICEEMD-De method is capable of correctly separating the different IMFs under noisy conditions, and also that the IMF can effectively identify the reflection signal of the end of the bolt.


title: Combining angular differential imaging and accurate polarimetry with SPHERE/IRDIS to characterize young giant exoplanets

abstract: Young giant exoplanets emit infrared radiation that can be linearly polarized up to several percent. This linear polarization can trace: 1) the presence of atmospheric cloud and haze layers, 2) spatial structure, e.g. cloud bands and rotational flattening, 3) the spin axis orientation and 4) particle sizes and cloud top pressure. We introduce a novel high-contrast imaging scheme that combines angular differential imaging (ADI) and accurate near-infrared polarimetry to characterize self-luminous giant exoplanets. We implemented this technique at VLT/SPHERE-IRDIS and developed the corresponding observing strategies, the polarization calibration and the data-reduction approaches. By combining ADI and polarimetry we can characterize planets that can be directly imaged with a very high signal-to-noise ratio. We use the IRDIS pupil-tracking mode and combine ADI and principal component analysis to reduce speckle noise. We take advantage of IRDIS' dual-beam polarimetric mode to eliminate differential effects that severely limit the polarimetric sensitivity (flat-fielding errors, differential aberrations and seeing), and thus further suppress speckle noise. To correct for instrumental polarization effects, we apply a detailed Mueller matrix model that describes the telescope and instrument and that has an absolute polarimetric accuracy $\leq0.1\%$. Using this technique we have observed the planets of HR 8799 and the (sub-stellar) companion PZ Tel B. Unfortunately, we do not detect a polarization signal in a first analysis. We estimate preliminary $1\sigma$ upper limits on the degree of linear polarization of $\sim1\%$ and $\sim0.1\%$ for the planets of HR 8799 and PZ Tel B, respectively. The achieved sub-percent sensitivity and accuracy show that our technique has great promise for characterizing exoplanets through direct-imaging polarimetry.


title: On Polyhedral Product Spaces over Polyhedral Joins

abstract: The construction of a simplicial complex given by polyhedral joins (introduced by Anton Ayzenberg), generalizes Bahri, Bendersky, Cohen and Gitler's $J$-construction and simplicial wedge construction. This article gives a cohomological decomposition of a polyhedral product over a polyhedral join for certain families of pairs of simplicial complexes. A formula for the Hilbert-Poincar\'{e} series is given, which generalizes Ayzenberg's formula for the moment-angle complex.


title: Topological analysis of nuclear pasta phases

abstract: In this paper the analysis of the result of numerical simulations of pasta phases using algebraic topology methods is presented. These considerations suggest that some phases can be further split into (sub)phases and therefore should be more refined in numerical simulations. The results presented in the paper can also be used to relate the Euler characteristic from numerical simulations to the geometry of the phases. The Betti numbers are used as they provide finer characterization of the phases. It is also shown that different boundary conditions give different outcomes.


title: A Paley-Wiener Type Theorem for Singular Measures on $\mathbb{T}$

abstract: For a fixed singular Borel probability measure $\mu$ on $\mathbb{T}$, we give several characterizations of when an entire function is the Fourier transform of some $f \in L^2(\mu)$. The first characterization is given in terms of criteria for sampling functions of the form $\hat{f}$ when $f \in L^2(\mu)$. The second characterization is given in terms of criteria for interpolation of bounded sequences on $\mathbb{N}_{0}$ by $\hat{f}$. Both characterizations use the construction of Fourier series for $f \in L^2(\mu)$ demonstrated in Herr and Weber via the Kaczmarz algorithm and classical results concerning the Cauchy transform of $\mu$.


title: Hamilton-Jacobi Reachability: A Brief Overview and Recent Advances

abstract: Hamilton-Jacobi (HJ) reachability analysis is an important formal verification method for guaranteeing performance and safety properties of dynamical systems; it has been applied to many small-scale systems in the past decade. Its advantages include compatibility with general nonlinear system dynamics, formal treatment of bounded disturbances, and the availability of well-developed numerical tools. The main challenge is addressing its exponential computational complexity with respect to the number of state variables. In this tutorial, we present an overview of basic HJ reachability theory and provide instructions for using the most recent numerical tools, including an efficient GPU-parallelized implementation of a Level Set Toolbox for computing reachable sets. In addition, we review some of the current work in high-dimensional HJ reachability to show how the dimensionality challenge can be alleviated via various general theoretical and application-specific insights.


title: Achieving Parsimony in Bayesian VARs with the Horseshoe Prior

abstract: In the context of a vector autoregression (VAR) model, or any multivariate regression model, the number of relevant predictors may be small relative to the information set available from which to build a prediction equation. It is well known that forecasts based off of (un-penalized) least squares estimates can overfit the data and lead to poor predictions. Since the Minnesota prior was proposed (Doan et al. (1984)), there have been many methods developed aiming at improving prediction performance. In this paper we propose the horseshoe prior (Carvalho et al. (2010), Carvalho et al. (2009)) in the context of a Bayesian VAR. The horseshoe prior is a unique shrinkage prior scheme in that it shrinks irrelevant signals rigorously to 0 while allowing large signals to remain large and practically unshrunk. In an empirical study, we show that the horseshoe prior competes favorably with shrinkage schemes commonly used in Bayesian VAR models as well as with a prior that imposes true sparsity in the coefficient vector. Additionally, we propose the use of particle Gibbs with backwards simulation (Lindsten et al. (2012), Andrieu et al. (2010)) for the estimation of the time-varying volatility parameters. We provide a detailed description of all MCMC methods used in the supplementary material that is available online.


title: Solving differential equations for Feynman integrals by expansions near singular points

abstract: We describe a strategy to solve differential equations for Feynman integrals by powers series expansions near singular points and to obtain high precision results for the corresponding master integrals. We consider Feynman integrals with two scales, i.e. nontrivially depending on one variable. The corresponding algorithm is oriented at situations where canonical form of the differential equations is impossible. We provide a computer implementation of our algorithm in a simple example of four-loop generalized sun-set integrals with three equal non-zero masses. Our code provides values of the master integrals at any given point on the real axis with a required accuracy and a given order of expansion in the regularization parameter $\epsilon$.


title: Massless on-shell box integral with arbitrary powers of propagators

abstract: The massless one-loop box integral with arbitrary indices in arbitrary space-time dimension $d$ is shown to reduce to a sum over three generalised hypergeometric functions. This result follows from the solution to the third order differential equation of hypergeometric type. To derive the differential equation, the Gr\"obner basis technique for integrals with noninteger powers of propagators was used. A complete set of recurrence relations from the Gr\"obner basis is presented. The first several terms in the $\varepsilon =(4-d)/2$ expansion of the result are given.


title: A posteriori multi-stage optimal trading under transaction costs and a diversification constraint

abstract: This paper presents a simple method for a posteriori (historical) multi-variate multi-stage optimal trading under transaction costs and a diversification constraint. Starting from a given amount of money in some currency, we analyze the stage-wise optimal allocation over a time horizon with potential investments in multiple currencies and various assets. Three variants are discussed, including unconstrained trading frequency, a fixed number of total admissable trades, and the waiting of a specific time-period after every executed trade until the next trade. The developed methods are based on efficient graph generation and consequent graph search, and are evaluated quantitatively on real-world data. The fundamental motivation of this work is preparatory labeling of financial time-series data for supervised machine learning.


title: Defining a Lingua Franca to Open the Black Box of a Na\"ive Bayes Recommender

abstract: Many AI systems have a black box nature that makes it difficult to understand how they make their recommendations. This can be unsettling, as the designer cannot be certain how the system will respond to novelty. To penetrate our Na\"ive Bayes recommender's black box, we first asked, what do we want to know from our system, and how can it be obtained? The answers led us to recursively define a common lexicon with the AI, a lingua franca, using the very items that the system ranks to create meta-symbols recognized by the system, and enabling us to understand the system's knowledge in plain terms and at different levels of abstraction. As one bonus, using its existing knowledge, the lingua franca can enable the system to extend recommendations to related, but entirely new areas, ameliorating the cold start problem. We also supplement the lingua franca with techniques for visualizing the system's knowledge state, develop metrics for evaluating the meaningfulness of terms in the lingua franca, and generalize the requirements for developing a similar lingua franca in other applications.