title: A Two-State Picture of Water and the Funnel of Life

abstract: Here I show that experimental and simulation data on liquid water using vibrational (infrared and Raman) and X-ray (absorption and emission) spectroscopies, as well as recent data from X-ray scattering, are fully consistent with a two-state picture of water. At ambient conditions there are fluctuations between a dominating high-density liquid (HDL) and a low-density form (LDL). These are related to the two forms of amorphous ice at very low temperature, high-density amorphous (HDA) and low-density amorphous (LDA), which interconvert in a first-order-like transition. This transition line is assumed to continue into the so-called No-mans land as a liquid-liquid transition and terminate in a critical point with very large fluctuations between the two liquid forms. These fluctuations extend in a funnel-like region up to ambient temperatures and pressures and give water its unusual properties which are fundamental to life. With this picture we find simple, intuitive explanations of the anomalous properties of water, such as the density maximum at 4 C, why ice floats, and why the compressibility and heat capacity grow as the liquid is cooled. We summarize by noting that in this picture, water is not a complicated liquid, but two normal liquids with a complicated relationship.


title: Generalization ability of region proposal networks for multispectral person detection

abstract: Multispectral person detection aims at automatically localizing humans in images that consist of multiple spectral bands. Usually, the visual-optical (VIS) and the thermal infrared (IR) spectra are combined to achieve higher robustness for person detection especially in insufficiently illuminated scenes. This paper focuses on analyzing existing detection approaches for their generalization ability. Generalization is a key feature for machine learning based detection algorithms that are supposed to perform well across different datasets. Inspired by recent literature regarding person detection in the VIS spectrum, we perform a cross-validation study to empirically determine the most promising dataset to train a well-generalizing detector. Therefore, we pick one reference Deep Convolutional Neural Network (DCNN) architecture and three different multispectral datasets. The Region Proposal Network (RPN) originally introduced for object detection within the popular Faster R-CNN is chosen as a reference DCNN. The reason is that a stand-alone RPN is able to serve as a competitive detector for two-class problems such as person detection. Furthermore, current state-of-the-art approaches initially apply an RPN followed by individual classifiers. The three considered datasets are the KAIST Multispectral Pedestrian Benchmark including recently published improved annotations for training and testing, the Tokyo Multi-spectral Semantic Segmentation dataset, and the OSU Color-Thermal dataset including recently released annotations. The experimental results show that the KAIST Multispectral Pedestrian Benchmark with its improved annotations provides the best basis to train a DCNN with good generalization ability compared to the other two multispectral datasets. On average, this detection model achieves a log-average Miss Rate (MR) of 29.74 % evaluated on the reasonable test subsets of the three datasets.


title: Dark-ages Reionization and Galaxy Formation Simulation -- XIX: Predictions of infrared excess and cosmic star formation rate density from UV observations

abstract: We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX) -- $\beta$ relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the Meraxes semi-analytic model, and explore three different parametrisations for the dust optical depths, linked to star formation rate, dust-to-gas ratio and gas column density respectively. A Bayesian approach is employed to statistically calibrate model free parameters including star formation efficiency, mass loading factor, dust optical depths and reddening slope directly against UV luminosity functions and colour-magnitude relations at z ~ 4-7. The best-fit models show excellent agreement with the observations. We calculate IRX using energy balance arguments, and find that the large intrinsic scatter in the IRX -$\beta$ plane is driven by the specific star formation rate. Additionally, the difference among the three dust models suggests a factor of two systematic uncertainty in the dust-corrected star formation rate when using the Meurer IRX - $\beta$ relation at z > 4.


title: Sequences with almost Poissonian Pair Correlations

abstract: Although a generic uniformly distributed sequence has Poissonian pair correlations, only one explicit example has been found up to now. Additionally, it is even known that many classes of uniformly distributed sequences, like van der Corput sequences, Kronecker sequences and LS sequences, do not have Poissonian pair correlations. In this paper, we show that van der Corput sequences and the Kronecker sequence for the golden mean are as close to having Poissonian pair correlations as possible: they both have $\alpha$-pair correlations for all $0 < \alpha < 1$ but not for $\alpha = 1$ which corresponds to Poissonian pair correlations.


title: Code Design Principles for Ultra-Reliable Random Access with Preassigned Patterns

abstract: We study medium access control layer random access under the assumption that the receiver can perform successive interference cancellation, without feedback. During recent years, a number of protocols with impressive error performance have been suggested for this channel model. However, the random nature of these protocols causes an error floor which limits their usability when targeting ultra-reliable communications. In very recent works by Paolini et al. and Boyd et. al., it was shown that if each user employs predetermined combinatorial access patterns, this error floor disappears. In this paper, we develop code design criteria for deterministic random access protocols in the ultra-reliability region, and build codes based on these principles. The suggested design methods are supported by simulations.


title: Displaying Things in Common to Encourage Friendship Formation: A Large Randomized Field Experiment

abstract: Friendship formation is important to online social network sites and to society, but can suffer from informational friction. In this study, we demonstrate that social networks may effectively use an IT-facilitated intervention -- displaying things in common (TIC) between users (mutual hometown, interest, education, work, city) -- to encourage friendship formation. Displaying TIC updates an individual's belief about the shared similarity with another and reduces information friction that may be hard to overcome in offline communication. In collaboration with an online social network, we conduct a randomized field experiment that randomly varies the prominence of different things in common when a user is browsing a non-friend's profile. The dyad-level exogenous variation, orthogonal to any (un)observed structural factors in viewer-profile's network, allows us to cleanly isolate the role of preferences for TIC in driving network formation and homophily. We find that displaying TICs to viewers can significantly increase their probability of sending a friend requests and forming friendships, and is especially effective for pairs of people who have little in common. We also find that displaying TIC can improve friendship formation for a wide range of viewers with different demographics and is more effective when the TICs are more surprising to the viewer.


title: Authenticated teleportation with one-sided trust

abstract: We introduce a protocol for authenticated teleportation, which can be proven secure even when the receiver does not trust their measurement devices, and is experimentally accessible. We use the technique of self-testing from the device-independent approach to quantum information, where we can characterise quantum states and measurements from the exhibited classical correlations alone. First, we derive self-testing bounds for the Bell state and Pauli $\sigma_X, \sigma_Z$ measurements, that are robust enough to be implemented in the lab. Then, we use these to determine a lower bound on the fidelity of an untested entangled state to be used for teleportation. Finally, we apply our results to propose an experimentally feasible protocol for one-sided device-independent authenticated teleportation. This can be interpreted as a first practical authentication of a quantum channel, with additional one-sided device-independence.


title: Partial data inverse problems and simultaneous recovery of boundary and coefficients for semilinear elliptic equations

abstract: We study various partial data inverse boundary value problems for the semilinear elliptic equation $\Delta u+ a(x,u)=0$ in a domain in $\mathbb R^n$ by using the higher order linearization technique introduced in [LLS 19, FO19]. We show that the Dirichlet-to-Neumann map of the above equation determines the Taylor series of $a(x,z)$ at $z=0$ under general assumptions on $a(x,z)$. The determination of the Taylor series can be done in parallel with the detection of an unknown cavity inside the domain or an unknown part of the boundary of the domain. The method relies on the solution of the linearized partial data Calder\'on problem [FKSU09], and implies the solution of partial data problems for certain semilinear equations $\Delta u+ a(x,u) = 0$ also proved in [KU19]. The results that we prove are in contrast to the analogous inverse problems for the linear Schr\"odinger equation. There recovering an unknown cavity (or part of the boundary) and the potential simultaneously are long-standing open problems, and the solution to the Calder\'on problem with partial data is known only in special cases when $n \geq 3$.


title: Micrometeoroid Events in LISA Pathfinder

abstract: The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,\mu\textrm{Ns}$ to 230$\,\mu\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.


title: Maximum Entropy Method for Solving the Turbulent Channel Flow Problem

abstract: There are two components in this work that allow solutions of the turbulent channel problem: one is the Galilean-transformed Navier-Stokes equation which gives a theoretical expression for the Reynolds stress; and the second the maximum entropy principle which provides the spatial distribution of turbulent kinetic energy. The first concept transforms the momentum balance for a control volume moving at the local mean velocity, breaking the momentum exchange down to its basic components. The Reynolds stress gradient budget confirms this alternative interpretation of the turbulence momentum balance, as validated with DNS data. The second concept of maximum entropy principle states that turbulent kinetic energy in fully-developed flows will distribute itself until the maximum entropy is attained while conforming to the physical constraints. By equating the maximum entropy state with maximum allowable (viscous) dissipation at a given Reynolds number, along with other constraints, we arrive at function forms (inner and outer) for the turbulent kinetic energy. This allows us to compute the Reynolds stress, then integrate it to obtain the velocity profiles in channel flows. The results agree well with direct numerical simulation (DNS) data at Re = 400 and 1000.


title: Discrete Fractional Integration Operators Along the Primes

abstract: We prove that the discrete fractional integration operators along the primes \[ T^{\lambda}_{\mathbb{P}}f(x) := \sum_{p} \frac{f(x-p)}{p^{\lambda}} \cdot \log p \] are bounded $\ell^p\to \ell^{p'}$ whenever $ \frac{1}{p'} < \frac{1}{p} - (1-\lambda), \ p > 1.$ Here, the sum runs only over prime $p$.


title: Determination of the scalar and vector polarizabilities of the cesium $6s \ ^2S_{1/2} \rightarrow 7s \ ^2S_{1/2}$ transition and implications for atomic parity non-conservation

abstract: Using recent high-precision measurements of electric dipole matrix elements of atomic cesium, we make an improved determination of the scalar ($\alpha$) and vector ($\beta$) polarizabilities of the cesium $6s \ ^2S_{1/2} \rightarrow 7s \ ^2S_{1/2} $ transition calculated through a sum-over-states method. We report values of $\alpha = -268.82 \ (30) \ a_0^3$ and $\beta = 27.139 \ (42) \ a_0^3$ with the highest precision to date. We find a discrepancy between our value of $\beta$ and the past preferred value, resulting in a significant shift in the value of the weak charge $Q_w$ of the cesium nucleus. Future work to resolve the differences in the polarizability will be critical for interpretation of parity non-conservation measurements in cesium, which have implications for physics beyond the Standard Model.


title: Gradient estimate for harmonic functions on K\"ahler manifolds

abstract: We prove a sharp integral gradient estimate for harmonic functions on noncompact K\"ahler manifolds. As application, we obtain a sharp estimate for the bottom of spectrum of the p-Laplacian and prove a splitting theorem for manifolds achieving this estimate.


title: Asymptotic behavior of age-structured and delayed Lotka-Volterra models

abstract: In this work we investigate some asymptotic properties of an age-structured Lotka-Volterra model, where a specific choice of the functional parameters allows us to formulate it as a delayed problem, for which we prove the existence of a unique coexistence equilibrium and characterize the existence of a periodic solution. We also exhibit a Lyapunov functional that enables us to reduce the attractive set to either the nontrivial equilibrium or to a periodic solution. We then prove the asymptotic stability of the nontrivial equilibrium where, depending on the existence of the periodic trajectory, we make explicit the basin of attraction of the equilibrium. Finally, we prove that these results can be extended to the initial PDE problem.


title: RACS: Rapid Analysis of ChIP-Seq data for contig based genomes

abstract: Background: Chromatin immunoprecipitation coupled to next generation sequencing (ChIP-Seq) is a widely used technique to investigate the function of chromatin-related proteins in a genome-wide manner. ChIP-Seq generates large quantities of data which can be difficult to process and analyse, particularly for organisms with contig based genomes. Contig-based genomes often have poor annotations for cis-elements, for example enhancers, that are important for gene expression. Poorly annotated genomes make a comprehensive analysis of ChIP-Seq data difficult and as such standardized analysis pipelines are lacking. Methods: We report a computational pipeline that utilizes traditional High-Performance Computing techniques and open source tools for processing and analysing data obtained from ChIP-Seq. We applied our computational pipeline "Rapid Analysis of ChIP-Seq data" (RACS) to ChIP-Seq data that was generated in the model organism Tetrahymena thermophila, an example of an organism with a genome that is available in contigs. Results: To test the performance and efficiency of RACs, we performed control ChIP-Seq experiments allowing us to rapidly eliminate false positives when analyzing our previously published data set. Our pipeline segregates the found read accumulations between genic and intergenic regions and is highly efficient for rapid downstream analyses. Conclusions: Altogether, the computational pipeline presented in this report is an efficient and highly reliable tool to analyze genome-wide ChIP-Seq data generated in model organisms with contig-based genomes. RACS is an open source computational pipeline available to download from: https://bitbucket.org/mjponce/racs --or-- https://gitrepos.scinet.utoronto.ca/public/?a=summary&p=RACS


title: Quantised Painlev\'e monodromy manifolds, Sklyanin and Calabi-Yau algebras

abstract: In this paper we study quantum del Pezzo surfaces belonging to a certain class. In particular we introduce the generalised Sklyanin-Painlev\'e algebra and characterise its PBW/PHS/Koszul properties. This algebra contains as limiting cases the generalised Sklyanin algebra, Etingof-Ginzburg and Etingof-Oblomkov-Rains quantum del Pezzo and the quantum monodromy manifolds of the Painlev\'e equations.


title: The Large High Altitude Air Shower Observatory (LHAASO) Science Book (2021 Edition)

abstract: Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment.

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


title: Lateral and axial resolution criteria in incoherent and coherent optics and holography, near- and far-field regimes

abstract: This work presents an overview of the spatial resolution criteria in classical optics, digital optics and holography. Although the classical Abbe and Rayleigh resolution criteria have been thoroughly discussed in the literature, there are still several issues which still need to be addressed, for example the axial resolution criterion for coherent and incoherent radiation, which is a crucial parameter of three-dimensional (3D) imaging, the resolution criteria in the Fresnel regime, and the lateral and axial resolution criteria in digital optics and holography. This work discusses these issues and provides a simple guide for which resolution criteria should be applied in each particular imaging scheme: coherent/incoherent, far- and near-field, lateral and axial resolution. Different resolution criteria such as two-point resolution and the resolution obtained from the image spectrum (diffraction pattern) are compared and demonstrated with simulated examples. Resolution criteria for spatial lateral and axial resolution are derived, and their application in imaging with coherent and incoherent (noncoherent) waves is considered. It is shown that for coherent light, the classical Abbe and Rayleigh resolution criteria do not provide an accurate estimation of the lateral and axial resolution. Lateral and axial resolution criteria based on an evaluation of the spectrum of the diffracted wave provide a more precise estimation of the resolution for coherent and incoherent light. It is also shown that resolution criteria derived in approximation of the far-field imaging regime can be applied for the near-field (Fresnel) regime.


title: Dynamic wake modulation induced by utility-scale wind turbine operation

abstract: Understanding wind turbine wake mixing and recovery is critical for improving the power generation and structural stability of downwind turbines in a wind farm. In the field, where incoming flow and turbine operation are constantly changing, wake recovery can be significantly influenced by dynamic wake modulation, which has not yet been explored. Here we present the first investigation of dynamic wake modulation in the near wake of a utility-scale turbine and quantify its relationship with changing conditions. This investigation is enabled using novel super-large-scale flow visualization with natural snowfall, providing unprecedented spatiotemporal resolution to resolve instantaneous changes of the wake envelope. These measurements reveal the significant influence of dynamic wake modulation on wake recovery. Further, our study uncovers the direct connection of dynamic wake modulation with operational parameters readily available to the turbine, paving the way for more precise wake prediction and control under field conditions for wind farm optimization.


title: Optimal Randomness in Swarm-Based Search

abstract: L\'{e}vy flights is a random walk where the step-lengths have a probability distribution that is heavy-tailed. It has been shown that L\'{e}vy flights can maximize the efficiency of resource searching in uncertain environments, and also movements of many foragers and wandering animals have been shown to follow a L\'{e}vy distribution. The reason mainly comes from that the L\'{e}vy distribution, has an infinite second moment, and hence is more likely to generate an offspring that is farther away from its parent. However, the investigation into the efficiency of other different heavy-tailed probability distributions in swarm-based searches is still insufficient up to now. For swarm-based search algorithms, randomness plays a significant role in both exploration and exploitation, or diversification and intensification. Therefore, it's necessary to discuss the optimal randomness in swarm-based search algorithms. In this study, CS is taken as a representative method of swarm-based optimization algorithms, and the results can be generalized to other swarm-based search algorithms. In this paper, four different types of commonly used heavy-tailed distributions, including Mittag-Leffler distribution, Pareto distribution, Cauchy distribution, and Weibull distribution, are considered to enhance the searching ability of CS. Then four novel CS algorithms are proposed and experiments are carried out on 20 benchmark functions to compare their searching performances. Finally, the proposed methods are used to system identification to demonstrate the effectiveness.


title: The Polarising Fragmentation Function and the Lambda polarisation in e+ e- processes

abstract: The surprising polarisation of Lambdas and other hyperons measured in many unpolarised hadronic processes, p N --> Lambda X, has been a long standing challenge for QCD phenomenology. One possible explanation was suggested, related to non perturbative properties of the quark hadronisation process, and encoded in the so-called Polarising Fragmentation Function (PFF). Recent Belle data have shown a non zero Lambda polarisation also in unpolarised e+ e- processes, e+ e- --> Lambda X and e+ e- --> Lambda h X. We consider the single inclusive case and the role of the PFFs. Adopting a simplified kinematics it is shown how they can originate a polarisation P_Lambda different from 0 and give explicit expressions for it in terms of the PFFs. Although the Belle data do not allow yet, in our approach, an extraction of the PFFs, some clear predictions are given, suggesting crucial measurements, and estimates of P_Lambda are computed, in qualitative agreement with the Belle data.


title: Robust temporal pumping in a magneto-mechanical topological insulator

abstract: The transport of energy through 1-dimensional (1D) waveguiding channels can be affected by sub-wavelength disorder, resulting in undesirable localization and backscattering phenomena. However, quantized disorder-resilient transport is observable in the edge currents of 2-dimensional (2D) topological band insulators with broken time-reversal symmetry. Topological pumps are able to reduce this higher-dimensional topological insulator phenomena to lower dimensionality by utilizing a pumping parameter (either space or time) as an artificial dimension. Here we demonstrate the first temporal topological pump that produces on-demand, robust transport of mechanical energy using a 1D magneto-mechanical metamaterial. We experimentally demonstrate that the system is uniquely resilient to defects occurring in both space and time Our findings open a new path towards exploration of higher-dimensional topological physics with time as a synthetic dimension.


title: Formation of interstellar SH$^+$ from vibrationally excited H$_2$: Quantum study of S$^+$ + H$_2$ $\rightleftarrows$ SH$^+$ + H reactions and inelastic collisions

abstract: The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3\Sigma^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H$_2$ reaction is endothermic by $\approx$ 9860 K, and requires at least two vibrational quanta on the H$_2$ molecule to yield significant reactivity. Quasi-classical calculations of the total formation rate constant for H$_2$($v=2$) are in very good agreement with the quantum results above 100K. Further quasi-classical calculations are then performed for $v=3$, 4, and 5 to cover all vibrationally excited H$_2$ levels significantly populated in dense photodissociation regions (PDR). The new calculated formation and destruction rate constants are two to six times larger than the previous ones and have been introduced in the Meudon PDR code to simulate the physical and illuminating conditions in the Orion bar prototypical PDR. New astrochemical models based on the new molecular data produce four times larger SH$^+$ column densities, in agreement with those inferred from recent ALMA observations of the Orion bar.


title: Uncertainty-Aware Data Aggregation for Deep Imitation Learning

abstract: Estimating statistical uncertainties allows autonomous agents to communicate their confidence during task execution and is important for applications in safety-critical domains such as autonomous driving. In this work, we present the uncertainty-aware imitation learning (UAIL) algorithm for improving end-to-end control systems via data aggregation. UAIL applies Monte Carlo Dropout to estimate uncertainty in the control output of end-to-end systems, using states where it is uncertain to selectively acquire new training data. In contrast to prior data aggregation algorithms that force human experts to visit sub-optimal states at random, UAIL can anticipate its own mistakes and switch control to the expert in order to prevent visiting a series of sub-optimal states. Our experimental results from simulated driving tasks demonstrate that our proposed uncertainty estimation method can be leveraged to reliably predict infractions. Our analysis shows that UAIL outperforms existing data aggregation algorithms on a series of benchmark tasks.


title: Robust Dense Mapping for Large-Scale Dynamic Environments

abstract: We present a stereo-based dense mapping algorithm for large-scale dynamic urban environments. In contrast to other existing methods, we simultaneously reconstruct the static background, the moving objects, and the potentially moving but currently stationary objects separately, which is desirable for high-level mobile robotic tasks such as path planning in crowded environments. We use both instance-aware semantic segmentation and sparse scene flow to classify objects as either background, moving, or potentially moving, thereby ensuring that the system is able to model objects with the potential to transition from static to dynamic, such as parked cars. Given camera poses estimated from visual odometry, both the background and the (potentially) moving objects are reconstructed separately by fusing the depth maps computed from the stereo input. In addition to visual odometry, sparse scene flow is also used to estimate the 3D motions of the detected moving objects, in order to reconstruct them accurately. A map pruning technique is further developed to improve reconstruction accuracy and reduce memory consumption, leading to increased scalability. We evaluate our system thoroughly on the well-known KITTI dataset. Our system is capable of running on a PC at approximately 2.5Hz, with the primary bottleneck being the instance-aware semantic segmentation, which is a limitation we hope to address in future work. The source code is available from the project website (http://andreibarsan.github.io/dynslam).


title: Ab Initio Simulations of Light Nuclear Systems Using Eigenvector Continuation and Auxiliary Field Monte Carlo

abstract: In this work, we discuss a new method for calculation of extremal eigenvectors and eigenvalues in systems or regions of parameter space where direct calculation is problematic. This technique relies on the analytic continuation of the power series expansion for the eigenvector around a point in the complex plane. We start this document by introducing the background material relevant to understand the basics of quantum mechanics and quantum field theories on the lattice, how we perform our numerical simulations, and how this relates to the nuclear physics we probe. We then move to the mathematical formalism of the eigenvector continuation, which is rooted in analytic function theory and linear algebra. We then discuss how these techniques are implemented numerically, with a discussion about the computational costs. Finally, we discuss applications of this method to full, quantum many-body systems. These include neutron matter, the Bose-Hubbard model, the Lipkin model, and the Coulomb interaction in light nuclei with LO chiral forces. These systems cover two categories of interest to the field: systems with a substantial sign problem, or systems that exhibit quantum phase transitions.


title: Measurement of the atom-surface van der Waals interaction by transmission spectroscopy in a wedged nano-cell

abstract: We demonstrate a method for measuring atom-surface interactions using transmission spectroscopy of thermal vapors confined in a wedged nano-cell. The wedged shape of the cell allows complementary measurements of both the bulk atomic vapor and atoms close to surfaces experiencing strong van der Waals atom-surface interaction. These are used to tightly constrain the dipole-dipole collisional parameters of a theoretical model for transmission spectra that accounts for atom-surface interactions, cavity effects, collisions with the surface of the cell and atomic motion. We illustrate this method on a cesium vapor in a sapphire cell, find $C_3=1.3\pm0.1$\,kHz.$\mu$m$^3$ and demonstrate that even the weakest of the van der Waals atom-surface interaction coefficients - for ground-state alkali atom transitions - can be determined with a very good precision. This result paves the way towards a precise quantitative characterization of atom-surface interactions in a wide range of atom-based nano-devices.


title: Regularization along central convergence on second and third Wiener chaoses

abstract: Consider $F$ an element of the second Wiener chaos with variance one. In full generality, we show that, for every integer $p\ge 1$, there exists $\eta_p>0$ such that if $\kappa_4(F)<\eta_p$ then the Malliavin derivative of $F$ admits a negative moment of order $p$. This entails that any sequence of random variables in the second Wiener chaos converging in distribution to a non--degenerated Gaussian is getting more regular as its distribution is getting close to the normal law. This substantially generalizes some recent findings contained in \cite{hu2014convergence,hu2015density,nourdin2016fisher} where analogous statements were given with additional assumptions which we are able to remove here. Moreover, we provide a multivariate version of this Theorem. Our main contribution concerns the case of the third Wiener chaos which is notoriously more delicate as one cannot anymore decompose the random variables into a linear combination of i.i.d. random variables. We still prove that the same phenomenon of regularization along central convergence occurs. Unfortunately, we are not able to provide a statement as strong as the previous one, but we can show that the usual non--degeneracy estimates of the Malliavin derivative given by the Carbery-Wright inequality can be improved by a factor three. Our proof introduces new techniques such that a specific Malliavin gradient enabling us to encode the distribution of the Malliavin derivative by the spectrum of some Gaussian matrix. This allows us to revisit the fourth moment phenomenon in terms of the behavior of its spectral radius.


title: NoBLE for lattice trees and lattice animals

abstract: We study lattice trees (LTs) and animals (LAs) on the nearest-neighbor lattice $\mathbb{Z}^d$ in high dimensions. We prove that LTs and LAs display mean-field behavior above dimension $16$ and $17$, respectively. Such results have previously been obtained by Hara and Slade in sufficiently high dimensions. The dimension above which their results apply was not yet specified. We rely on the non-backtracking lace expansion (NoBLE) method that we have recently developed. The NoBLE makes use of an alternative lace expansion for LAs and LTs that perturbs around non-backtracking random walk rather than simple random walk, leading to smaller corrections. The NoBLE method then provides a careful computational analysis that improves the dimension above which the result applies. Universality arguments predict that the upper critical dimension, above which our results apply, is equal to $d_c=8$ for both models, as is known for sufficiently spread-out models by the results of Hara and Slade mentioned earlier. The main ingredients in this paper are (a) to derive a non-backtracking lace expansion for the LT and LA two-point functions; (b) to bound the non-backtracking lace-expansion coefficients, thus showing that our general NoBLE methodology can be applied; and (c) to obtain sharp numerical bounds on the coefficients. Our proof is complemented by a computer-assisted numerical analysis that verifies that the necessary bounds used in the NoBLE are satisfied.


title: Many-Body Expanded Full Configuration Interaction. II. Strongly Correlated Regime

abstract: In this second part of our series on the recently proposed many-body expanded full configuration interaction (MBE-FCI) method, we introduce the concept of multideterminantal expansion references. Through theoretical arguments and numerical validations, the use of this class of starting points is shown to result in a focussed compression of the MBE decomposition of the FCI energy, thus allowing chemical problems dominated by strong correlation to be addressed by the method. The general applicability and performance enhancements of MBE-FCI are verified for standard stress tests such as the bond dissociations in H$_2$O, N$_2$, C$_2$, and a linear H$_{10}$ chain. Furthermore, the benefits of employing a multideterminantal expansion reference in accelerating calculations of high accuracy are discussed, with an emphasis on calculations in extended basis sets. As an illustration of this latter quality of the MBE-FCI method, results for H$_2$O and C$_2$ in basis sets ranging from double- to pentuple-$\zeta$ quality are presented, demonstrating near-ideal parallel scaling on up to almost $25000$ processing units.


title: Gravitational Waves from Mirror World

abstract: In this paper we consider the properties of the 10 confirmed by the LIGO Collaboration gravitational wave signals from the black hole mergers. We want to explain non-observation of electromagnetic counterpart and higher then expected merging rates of these events, assuming the existence of their sources in the hidden mirror universe. Mirror matter, which interacts with our world only through gravity, is a candidate of dark matter and its density can exceed ordinary matter density five times. Since mirror world is considered to be colder, star formation there started earlier and mirror black holes had more time to pick up the mass and to create more binary systems within the LIGO reachable zone. In total, we estimate factor of 15 amplification of black holes merging rate in mirror world with respect to our world, which is consistent with the LIGO observations.


title: Cosmological Acceleration as a Consequence of Quantum de Sitter Symmetry

abstract: Physicists usually understand that physics cannot (and should not) derive that $c\approx 3\cdot 10^8m/s$ and $\hbar \approx 1.054\cdot 10^{-34}kg\cdot m^2/s$. At the same time they usually believe that physics should derive the value of the cosmological constant $\Lambda$ and that the solution of the dark energy problem depends on this value. However, background space in General Relativity (GR) is only a classical notion while on quantum level symmetry is defined by a Lie algebra of basic operators. We prove that the theory based on Poincare Lie algebra is a special degenerate case of the theories based on de Sitter (dS) or anti-de Sitter (AdS) Lie algebras in the formal limit $R\to\infty$ where R is the parameter of contraction from the latter algebras to the former one, and $R$ has nothing to do with the radius of background space. As a consequence, $R$ is necessarily finite, is fundamental to the same extent as $c$ and $\hbar$, and a question why $R$ is as is does not arise. Following our previous publications, we consider a system of two free bodies in dS quantum mechanics and show that in semiclassical approximation the cosmological dS acceleration is necessarily nonzero and is the same as in GR if the radius of dS space equals $R$ and $\Lambda=3/R^2$. This result follows from basic principles of quantum theory. It has nothing to do with existence or nonexistence of dark energy and therefore for explaining cosmological acceleration dark energy is not needed. The result is obtained without using the notion of dS background space (in particular, its metric and connection) but simply as a consequence of quantum mechanics based on the dS Lie algebra. Therefore, $\Lambda$ has a physical meaning only on classical level and the cosmological constant problem and the dark energy problem do not arise.


title: Variational training of neural network approximations of solution maps for physical models

abstract: A novel solve-training framework is proposed to train neural network in representing low dimensional solution maps of physical models. Solve-training framework uses the neural network as the ansatz of the solution map and train the network variationally via loss functions from the underlying physical models. Solve-training framework avoids expensive data preparation in the traditional supervised training procedure, which prepares labels for input data, and still achieves effective representation of the solution map adapted to the input data distribution. The efficiency of solve-training framework is demonstrated through obtaining solutions maps for linear and nonlinear elliptic equations, and maps from potentials to ground states of linear and nonlinear Schr\"odinger equations.


title: A Comparison of the Diffuser Method Versus the Defocus Method for Performing High-Precision Photometry with Small Telescope Systems

abstract: This paper compares the performance of two different high-precision, photometric measurement techniques for bright (<11 magnitude) stars using the small telescope systems that today's amateur astronomers typically use. One technique is based on recent work using a beam-shaping diffuser method (Stefansson et al., (2017).) The other is based on the widely used "defocusing" method. We also developed and used a statistical photometric performance model to better understand the error components of the measurements to identify and quantify any difference in performance between the two methods. AstroImageJ (Collins et al. (2017)) was used for the exoplanet image analysis to provide the measured values and exoplanet models described in this study. Both methods were used at the Mark Slade Remote Observatory (MSRO) to conduct in-transit exoplanet observations of exoplanets HAT-P-30b/WASP-51b, HAT-P-16b, and a partial of WASP-93b. Observations of exoplanets KELT-1b and K2-100b and other stars were also performed at the MSRO to further understand and characterize the performance of the diffuser method under various sky conditions. In addition, both in-transit and out-of-transit observations of exoplanets HAT-P-23b, HAT-P-33b, and HAT-P-34b were performed at the Conti Private Observatory. We found that for observing bright stars, the diffuser method outperformed the defocus method when using small telescopes with poor tracking. We also found the diffuser method noticeably reduced the scintillation noise compared with the defocus method and provided high-precision results in typical, average sky conditions through all lunar phases. For small telescopes using excellent auto-guiding techniques and effective calibration procedures, we found the defocus method was equal to or in some cases better than the diffuser method when observing with good-to-excellent sky conditions.


title: Fast Neural Network Approach for Direct Covariant Forces Prediction in Complex Multi-Element Extended Systems

abstract: Neural network force field (NNFF) is a method for performing regression on atomic structure-force relationships, bypassing expensive quantum mechanics calculation which prevents the execution of long ab-initio quality molecular dynamics simulations. However, most NNFF methods for complex multi-element atomic systems indirectly predict atomic force vectors by exploiting just atomic structure rotation-invariant features and the network-feature spatial derivatives which are computationally expensive. We develop a staggered NNFF architecture exploiting both rotation-invariant and covariant features separately to directly predict atomic force vectors without using spatial derivatives, thereby reducing expensive structural feature calculation by ~180-480x. This acceleration enables us to develop NNFF which directly predicts atomic forces in complex ternary and quaternary-element extended systems comprised of long polymer chains, amorphous oxide, and surface chemical reactions. The staggered rotation-invariant-covariant architecture described here can also directly predict complex covariant vector outputs from local physical structures in domains beyond computational material science.


title: Low redshift constraints on energy-momentum-powered gravity models

abstract: There has been recent interest in the cosmological consequences of energy-momentum-powered gravity models, in which the matter side of Einstein's equations is modified by the addition of a term proportional to some power, $n$, of the energy-momentum tensor, in addition to the canonical linear term. In this work we treat these models as phenomenological extensions of the standard $\Lambda$CDM, containing both matter and a cosmological constant. We also quantitatively constrain the additional model parameters using low redshift background cosmology data that are specifically from Type Ia supernovas and Hubble parameter measurements. We start by studying specific cases of these models with fixed values of $n,$ which lead to an analytic expression for the Friedmann equation; we discuss both their current constraints and how the models may be further constrained by future observations of Type Ia supernovas for WFIRST complemented by measurements of the redshift drift by the ELT. We then consider and constrain a more extended parameter space, allowing $n$ to be a free parameter and considering scenarios with and without a cosmological constant. These models do not solve the cosmological constant problem per se. Nonetheless these models can phenomenologically lead to a recent accelerating universe without a cosmological constant at the cost of having a preferred matter density of around $\Omega_M\sim0.4$ instead of the usual $\Omega_M\sim0.3$. Finally we also briefly constrain scenarios without a cosmological constant, where the single component has a constant equation of state which needs not be that of matter; we provide an illustrative comparison of this model with a more standard dynamical dark energy model with a constant equation of state.


title: Skin Lesion Classification Using CNNs with Patch-Based Attention and Diagnosis-Guided Loss Weighting

abstract: Objective: This work addresses two key problems of skin lesion classification. The first problem is the effective use of high-resolution images with pretrained standard architectures for image classification. The second problem is the high class imbalance encountered in real-world multi-class datasets. Methods: To use high-resolution images, we propose a novel patch-based attention architecture that provides global context between small, high-resolution patches. We modify three pretrained architectures and study the performance of patch-based attention. To counter class imbalance problems, we compare oversampling, balanced batch sampling, and class-specific loss weighting. Additionally, we propose a novel diagnosis-guided loss weighting method which takes the method used for ground-truth annotation into account. Results: Our patch-based attention mechanism outperforms previous methods and improves the mean sensitivity by 7%. Class balancing significantly improves the mean sensitivity and we show that our diagnosis-guided loss weighting method improves the mean sensitivity by 3% over normal loss balancing. Conclusion: The novel patch-based attention mechanism can be integrated into pretrained architectures and provides global context between local patches while outperforming other patch-based methods. Hence, pretrained architectures can be readily used with high-resolution images without downsampling. The new diagnosis-guided loss weighting method outperforms other methods and allows for effective training when facing class imbalance. Significance: The proposed methods improve automatic skin lesion classification. They can be extended to other clinical applications where high-resolution image data and class imbalance are relevant.


title: QuantumATK: An integrated platform of electronic and atomic-scale modelling tools

abstract: QuantumATK is an integrated set of atomic-scale modelling tools developed since 2003 by professional software engineers in collaboration with academic researchers. While different aspects and individual modules of the platform have been previously presented, the purpose of this paper is to give a general overview of the platform. The QuantumATK simulation engines enable electronic-structure calculations using density functional theory or tight-binding model Hamiltonians, and also offers bonded or reactive empirical force fields in many different parametrizations. Density functional theory is implemented using either a plane-wave basis or expansion of electronic states in a linear combination of atomic orbitals. The platform includes a long list of advanced modules, including Green's-function methods for electron transport simulations and surface calculations, first-principles electron-phonon and electron-photon couplings, simulation of atomic-scale heat transport, ion dynamics, spintronics, optical properties of materials, static polarization, and more. Seamless integration of the different simulation engines into a common platform allows for easy combination of different simulation methods into complex workflows. Besides giving a general overview and presenting a number of implementation details not previously published, we also present four different application examples. These are calculations of the phonon-limited mobility of Cu, Ag and Au, electron transport in a gated 2D device, multi-model simulation of lithium ion drift through a battery cathode in an external electric field, and electronic-structure calculations of the composition-dependent band gap of SiGe alloys.


title: Failure of the Brauer-Manin principle for a simply connected fourfold over a global function field, via orbifold Mordell

abstract: Almost one decade ago, Poonen constructed the first examples of algebraic varieties over global fields for which Skorobogatov's etale Brauer-Manin obstruction does not explain the failure of the Hasse principle. By now, several constructions are known, but they all share common geometric features such as large fundamental groups. In this paper, we construct simply connected fourfolds over global fields of positive characteristic for which the Brauer-Manin machinery fails. Contrary to earlier work in this direction, our construction does not rely on major conjectures. Instead, we establish a new diophantine result of independent interest: a Mordell-type theorem for Campana's "geometric orbifolds" over function fields of positive characteristic. Along the way, we also construct the first example of simply connected surface of general type over a global field with a non-empty, but non-Zariski dense set of rational points.


title: Collaborative and Privacy-Preserving Machine Teaching via Consensus Optimization

abstract: In this work, we define a collaborative and privacy-preserving machine teaching paradigm with multiple distributed teachers. We focus on consensus super teaching. It aims at organizing distributed teachers to jointly select a compact while informative training subset from data hosted by the teachers to make a learner learn better. The challenges arise from three perspectives. First, the state-of-the-art pool-based super teaching method applies mixed-integer non-linear programming (MINLP) which does not scale well to very large data sets. Second, it is desirable to restrict data access of the teachers to only their own data during the collaboration stage to mitigate privacy leaks. Finally, the teaching collaboration should be communication-efficient since large communication overheads can cause synchronization delays between teachers. To address these challenges, we formulate collaborative teaching as a consensus and privacy-preserving optimization process to minimize teaching risk. We theoretically demonstrate the necessity of collaboration between teachers for improving the learner's learning. Furthermore, we show that the proposed method enjoys a similar property as the Oracle property of adaptive Lasso. The empirical study illustrates that our teaching method can deliver significantly more accurate teaching results with high speed, while the non-collaborative MINLP-based super teaching becomes prohibitively expensive to compute.


title: Sparse multiresolution representations with adaptive kernels

abstract: Reproducing kernel Hilbert spaces (RKHSs) are key elements of many non-parametric tools successfully used in signal processing, statistics, and machine learning. In this work, we aim to address three issues of the classical RKHS based techniques. First, they require the RKHS to be known a priori, which is unrealistic in many applications. Furthermore, the choice of RKHS affects the shape and smoothness of the solution, thus impacting its performance. Second, RKHSs are ill-equipped to deal with heterogeneous degrees of smoothness, i.e., with functions that are smooth in some parts of their domain but vary rapidly in others. Finally, the computational complexity of evaluating the solution of these methods grows with the number of data points, rendering these techniques infeasible for many applications. Though kernel learning, local kernel adaptation, and sparsity have been used to address these issues, many of these approaches are computationally intensive or forgo optimality guarantees. We tackle these problems by leveraging a novel integral representation of functions in RKHSs that allows for arbitrary centers and different kernels at each center. To address the complexity issues, we then write the function estimation problem as a sparse functional program that explicitly minimizes the support of the representation leading to low complexity solutions. Despite their non-convexity and infinite dimensionality, we show these problems can be solved exactly and efficiently by leveraging duality, and we illustrate this new approach in simulated and real data.


title: Locally conformally symplectic reduction of the cotangent bundle

abstract: In a previous article, we introduced a reduction procedure for locally conformally symplectic manifolds at any regular value of the momentum mapping. We use this construction to prove an analogue of a well-known theorem in the symplectic setting about the reduction of cotangent bundles.


title: Carbon radio recombination lines from gigahertz to megahertz frequencies towards Orion A

abstract: Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $\mu$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $\mu$m-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15--230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion~A (M42). We complement these observations with SOFIA observations of the 158 $\mu$m-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C$^{+}$/C/CO interface. By modeling the emission of the $158$~$\mu$m-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure $>5\times10^{7}$ K cm$^{-3}$ and a radiation field $G_{0}\approx10^{5}$ close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158 $\mu$m-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C$^{+}$/C/CO interface.


title: Online and Offline Greedy Algorithms for Routing with Switching Costs

abstract: Motivated by the use of high speed circuit switches in large scale data centers, we consider the problem of circuit switch scheduling. In this problem we are given demands between pairs of servers and the goal is to schedule at every time step a matching between the servers while maximizing the total satisfied demand over time. The crux of this scheduling problem is that once one shifts from one matching to a different one a fixed delay $\delta$ is incurred during which no data can be transmitted. For the offline version of the problem we present a $(1-\frac{1}{e}-\epsilon)$ approximation ratio (for any constant $\epsilon >0$). Since the natural linear programming relaxation for the problem has an unbounded integrality gap, we adopt a hybrid approach that combines the combinatorial greedy with randomized rounding of a different suitable linear program. For the online version of the problem we present a (bi-criteria) $ ((e-1)/(2e-1)-\epsilon)$-competitive ratio (for any constant $\epsilon >0$ ) that exceeds time by an additive factor of $O(\frac{\delta}{\epsilon})$. We note that no uni-criteria online algorithm is possible. Surprisingly, we obtain the result by reducing the online version to the offline one.


title: Spin Berry phase in a helical edge state: Sz nonconservation and transport signatures

abstract: Topological protection of edge state in quantum spin Hall systems relies only on time-reversal symmetry. Hence, S z conservation on the edge can be relaxed which can have an interferometric manifestation in terms of spin Berry phase. Primarily it could lead to the generation of spin Berry phase arising from a closed loop dynamics of electrons. Our work provides a minimal framework to generate and detect these effects by employing both spin-unpolarized and spin-polarized leads. We show that spin-polarized leads could lead to resonances or anti-resonances in the two-terminal conductance of the interferometer. We further show that the positions of these anti-resonances (as a function of energy of the incident electron) get shifted owing to the presence of spin Berry phase. Finally, we present simulations of a device setup using KWANTpackage which put our theoretical predictions on firm footing.


title: W-Symmetries of Ito stochastic differential equations

abstract: We discuss W-symmetries of Ito stochastic differential equations, introduced in a recent paper by Gaeta and Spadaro [J. Math. Phys. 2017]. In particular, we discuss the general form of acceptable generators for continuous (Lie-point) W-symmetry, arguing they are related to the (linear) conformal group, and how W-symmetries can be used in the integration of Ito stochastic equations along Kozlov theory for standard (deterministic or random) symmetries. It turns out this requires, in general, to consider more general classes of stochastic equations than just Ito ones.


title: P3DFFT: a framework for parallel computations of Fourier transforms in three dimensions

abstract: Fourier and related transforms is a family of algorithms widely employed in diverse areas of computational science, notoriously difficult to scale on high-performance parallel computers with large number of processing elements (cores). This paper introduces a popular software package called P3DFFT implementing Fast Fourier Transforms (FFT) in three dimensions (3D) in a highly efficient and scalable way. It overcomes a well-known scalability bottleneck of 3D FFT implementations by using two-dimensional domain decomposition. Designed for portable performance, P3DFFT achieves excellent timings for a number of systems and problem sizes. On Cray XT5 system P3DFFT attains 45% efficiency in weak scaling from 128 to 65,536 computational cores. Library features include Fourier and Chebyshev transforms, Fortran and C interfaces, in- and out-of-place transforms, uneven data grids, single and double precision. P3DFFT is available as open source at http://code.google.com/p/p3dfft/. This paper discusses P3DFFT implementation and performance in a way that helps guide the user in making optimal choices for parameters of their runs.


title: A weighted setting for the stationary Navier Stokes equations under singular forcing

abstract: In two dimensions, we show existence of solutions to the stationary Navier Stokes equations on weighted spaces $\mathbf{H}^1_0(\omega,\Omega) \times L^2(\omega,\Omega)$, where the weight belongs to the Muckenhoupt class $A_2$. We show how this theory can be applied to obtain a priori error estimates for approximations of the solution to the Navier Stokes problem with singular sources.


title: Network Coding Gaps for Completion Times of Multiple Unicasts

abstract: We study network coding gaps for the problem of makespan minimization of multiple unicasts. In this problem distinct packets at different nodes in a network need to be delivered to a destination specific to each packet, as fast as possible. The network coding gap specifies how much coding packets together in a network can help compared to the more natural approach of routing. While makespan minimization using routing has been intensely studied for the multiple unicasts problem, no bounds on network coding gaps for this problem are known. We develop new techniques which allow us to upper bound the network coding gap for the makespan of $k$ unicasts, proving this gap is at most polylogarithmic in $k$. Complementing this result, we show there exist instances of $k$ unicasts for which this coding gap is polylogarithmic in $k$. Our results also hold for average completion time, and more generally any $\ell_p$ norm of completion times.


title: Twisted Dolbeault cohomology of nilpotent Lie algebras

abstract: It is well known that cohomology of any non-trivial 1-dimensional local system on a nilmanifold vanishes (this result is due to L. Alaniya). A complex nilmanifold is a quotient of a nilpotent Lie group equipped with a left-invariant complex structure by an action of a discrete, co-compact subgroup. We prove a Dolbeault version of Alaniya's theorem, showing that the Dolbeault cohomology of a nilpotent Lie algebra with coefficients in any non-trivial 1-dimensional local system vanishes. Note that the Dolbeault cohomology of the corresponding local system on the manifold is not necessarily zero. This implies that the twisted version of Console-Fino theorem is false (Console-Fino proved that the Dolbeault cohomology of a complex nilmanifold is equal to the Dolbeault cohomology of its Lie algebra). As an application, we give a new proof of a theorem due to H. Sawai, who obtained an explicit description of LCK nilmanifolds. An LCK structure on a manifold $M$ is a K\"ahler structure on its cover $\tilde M$ such that the deck transform map acts on $\tilde M$ by homotheties. We show that any complex nilmanifold admitting an LCK structure is Vaisman, and is obtained as a compact quotient of the product of a Heisenberg group and the real line.


title: Optimal pulse propagation in an inhomogeneously gas-filled hollow-core fiber

abstract: We study optical pulse propagation through a hollow-core fiber filled with a radially inhomogeneous cloud of cold atoms. A co-propagating control field establishes electromagnetically induced transparency. In analogy to a graded index fiber, the pulse experiences micro-lensing and the transmission spectrum becomes distorted. Based on a two-layer model of the complex index of refraction, we can analytically understand the cause of the aberration, which is corroborated by numerical simulations for a radial Gaussian-shaped function. With these insights, we show that the spectral distortions can be rectified by choosing an optimal detuning from one-photon resonance.


title: Continued fractions and Bessel functions

abstract: Elementary transformations of equations $A\psi=\lambda\psi$ are considered. The invertibility condition (Theorem 1) is established and similar transformations of Riccati equations in the case of second order differential operator $A$ are constructed (Theorem 2). Applications to continuous fractions for Bessel functions and Chebyshev polynomials are established. It is shown particularly that the elementary solutions of Bessel equations are related to a fixed point transformations of Riccati equations.


title: Higher order nonlocal operator method

abstract: We extend the nonlocal operator method to higher order scheme by using a higher order Taylor series expansion of the unknown field. Such a higher order scheme improves the original nonlocal operator method proposed by the authors in [A nonlocal operator method for solving partial differential equations], which can only achieve one-order convergence. The higher order nonlocal operator method obtains all partial derivatives with specified maximal order simultaneously without resorting to shape functions. The functional based on the nonlocal operators converts the construction of residual and stiffness matrix into a series of matrix multiplication on the nonlocal operator matrix. Several numerical examples solved by strong form or weak form are presented to show the capabilities of this method.


title: Decision Making with Machine Learning and ROC Curves

abstract: The Receiver Operating Characteristic (ROC) curve is a representation of the statistical information discovered in binary classification problems and is a key concept in machine learning and data science. This paper studies the statistical properties of ROC curves and its implication on model selection. We analyze the implications of different models of incentive heterogeneity and information asymmetry on the relation between human decisions and the ROC curves. Our theoretical discussion is illustrated in the context of a large data set of pregnancy outcomes and doctor diagnosis from the Pre-Pregnancy Checkups of reproductive age couples in Henan Province provided by the Chinese Ministry of Health.


title: Modeling the excess velocity of low-viscous Taylor droplets in square microchannels

abstract: Microscopic multiphase flows have gained broad interest due to their capability to transfer processes into new operational windows and achieving significant process intensification. However, the hydrodynamic behavior of Taylor droplets is not yet entirely understood. In this work, we introduce a model to determine the excess velocity of Taylor droplets in square microchannels. This velocity difference between the droplet and the total superficial velocity of the flow has a direct influence on the droplet residence time and is linked to the pressure drop. Since the droplet does not occupy the entire channel cross section, it enables the continuous phase to bypass the droplet through the corners. A consideration of the continuity equation generally relates the excess velocity to the mean flow velocity. We base the quantification of the bypass flow on a correlation for the droplet cap deformation from its static shape. The cap deformation reveals the forces of the flowing liquids exerted onto the interface and allows estimating the local driving pressure gradient for the bypass flow. The characterizing parameters are identified as the bypass length, the wall film thickness, the viscosity ratio between both phases and the Ca-number. The proposed model is adapted with a stochastic, metaheuristic optimization approach based on high-speed camera measurements. In addition, our model is successfully verified with published empirical data.


title: From GenderMag to InclusiveMag: An Inclusive Design Meta-Method

abstract: How can software practitioners assess whether their software supports diverse users? Although there are empirical processes that can be used to find "inclusivity bugs" piecemeal, what is often needed is a systematic inspection method to assess soft-ware's support for diverse populations. To help fill this gap, this paper introduces InclusiveMag, a generalization of GenderMag that can be used to generate systematic inclusiveness methods for a particular dimension of diversity. We then present a multi-case study covering eight diversity dimensions, of eight teams' experiences applying InclusiveMag to eight under-served populations and their "mainstream" counterparts.


title: Fixing Inclusivity Bugs for Information Processing Styles and Learning Styles

abstract: Most software systems today do not support cognitive diversity. Further, because of differences in problem-solving styles that cluster by gender, software that poorly supports cognitive diversity can also embed gender biases. To help software professionals fix gender bias "bugs" related to people's problem-solving styles for information processing and learning of new software we collected inclusivity fixes from three sources. The first two are empirical studies we conducted: a heuristics-driven user study and a field research industry study. The third is data that we obtained about a before/after user study of inclusivity bugs. The resulting seven potential inclusivity fixes show how to debug software to be more inclusive for diverse problem-solving styles.


title: Spin-Layer- and Spin-Valley-Locking in CVD-Grown AA'- and AB-Stacked Tungsten-Disulfide Bilayers

abstract: Valley-selective optical selection rules and a spin-valley locking in transition-metal dichalcogenide (TMDC) monolayers are at the heart of "valleytronic physics", which exploits the valley degree of freedom and has been a major research topic in recent years. In contrast, valleytronic properties of TMDC bilayers have not been in the focus so much by now. Here, we report on the valleytronic properties and optical characterization of bilayers of WS2 as a representative TMDC material. In particular, we study the influence of the relative layer alignment in TMDC homo-bilayer samples on their polarization-dependent optical properties. Therefore, CVD-grown WS2 bilayer samples have been prepared that favor either the inversion symmetric AA' stacking or AB stacking without inversion symmetry during synthesis. Subsequently, a detailed analysis of reflection contrast and photoluminescence spectra under different polarization conditions has been performed. We observe circular and linear dichroism of the photoluminescence that is more pronounced for the AB stacking configuration. Our experimental findings are supported by theoretical calculations showing that the observed dichroism can be linked to optical selection rules, that maintain the spin-valley locking in the AB-stacked WS2 bilayer, whereas a spin-layer-locking is present the inversion symmetric AA' bilayer instead. Furthermore, our theoretical calculations predict a small relative shift of the excitonic resonances in both stacking configurations, which is also experimentally observed.


title: Allostery and conformational changes upon binding as generic features of proteins: a high-dimension geometrical approach

abstract: A growing number of experimental evidence shows that it is general for a ligand binding protein to have a potential for allosteric regulation and for further evolution. In addition, such proteins generically change their conformation upon binding. O. Rivoire has recently proposed an evolutionary scenario that explains these properties as a generic byproduct of selection for exquisite discrimination between very similar ligands. The initial claim was supported by two classes of basic examples: continuous protein models with small numbers of degrees of freedom, on which the development of a conformational switch was established, and a 2-dimensional spin glass model supporting the rest of the statement. This work aimed to clarify the implication of the exquisite discrimination for smooth models with large number of degrees of freedom, the situation closer to real biological systems. With the help of differential geometry, jet-space analysis, and transversality theorems, it is shown that the claim holds true for any generic flexible system that can be described in terms of smooth manifolds. The result suggests that, indeed, evolutionary solutions to the exquisite discrimination problem, if exist, are located near a codimension-1 subspace of the appropriate genotypical space. This constraint, in turn, gives rise to a potential for the allosteric regulation of the discrimination via generic conformational changes upon binding.


title: Crossing, Modular Averages and $N \leftrightarrow k $ in WZW Models

abstract: We consider the construction of genus zero correlators of $SU(N)_k$ WZW models involving two Kac Moody primaries in the fundamental and two in the anti-fundamental representation from modular averaging of the contribution of the vacuum conformal block. In cases where we find the orbit of the vacuum conformal block to be finite, modular averaging reproduces the exact result for the correlators. In other cases, we perform the modular averaging numerically, the results are in agreement with the exact answers. We find a close relationship between the modular averaging sums of the theories related by level rank duality. We establish a one to one correspondence between elements of the orbits of the vacuum conformal blocks of dual theories. The contributions of paired terms to their respective correlators are simply related. One consequence of this is that the ratio between the OPE coefficients associated with dual correlators can be obtained analytically without performing the sums involved in the modular averagings. The pairing of terms in the modular averaging sums for dual theories suggests an interesting connection between level rank duality and semi-classical holographic computations of the correlators in the theories.


title: Conditions for stable equilibrium in Cournot duopoly models with tax evasion and time delay

abstract: We provide conditions for stable equilibrium in Cournot duopoly models with tax evasion and time delay. We prove that our conditions actually imply asymptotically stable equilibrium and delay independence. Conditions include the same marginal cost and equal probability for evading taxes. We give examples of cost and inverse demand functions satisfying the proposed conditions. Some economic interpretations of our results are also included.


title: Geodesic mappings and concircular vector fields

abstract: In the present paper we study geodesic mappings of special pseudo-Riemannian manifolds called $V_n(K)$-spaces. We prove that the set of solutions of the system of equations of geodesic mappings on $V_n(K)$-spaces $(K\neq0)$ forms a special Jordan algebra and the set of solutions generated by consircular fields is an ideal of this algebra. We show that pseudo-Riemannian manifolds admitting a concircular field of the basic type form the class of manifolds closed with respect to the geodesic mappings.


title: Refractive index matching (RIM) using double-binary liquid-liquid mixtures

abstract: Within the last decade microscopic multiphase flows have gained broad interest. An exact understanding of the underlying hydrodynamic interrelations is the key for successful reactor layout and reaction control. To examine the local hydrodynamic behavior, non-invasive optical measurements techniques like particle tracking velocimetry (PTV) or (micro)particle image velocimetry ((\textmu)PIV) are the method of choice, since they provide precise velocity measurement with excellent spatial resolution. Such optical approaches require refractive index matching (RIM) of the involved flow phases to prevent optical distortion due to light refraction and reflection at the interfaces. Established RIM approaches often provide a single one degree of freedom which is sufficient to match the RI of the flow phases solely. With that, the material properties ($Oh$ number) are fixed and the relevant dimensionless numbers ($Ca$, $Re$) may only be altered hydrodynamically or geometrically. To avoid expansive geometric scaling of the microchannels, we propose an approach using two binary mixtures (double-binary mixtures) to introduce an additional degree of freedom. The approach allows examining liquid-liquid two-phase flows at a distinct velocity while being able to change the material combination ($Oh$-Number). Therefore $Ca$ and $Re$ can be chosen individually and the RIM provides undisturbed optical access. We present 4 different binary mixtures to be used, e.g. with Taylor droplets. The relevant material parameters are successfully correlated to measurement data, which delivers a system of equations that determines the mass fractions and the velocities to address $Re$ and $Ca$ individually. A proof-of-principle for the proposed double binary mixture RIM-approach is successfully established using \textmu PIV raw images.


title: Stochastic Averaging of The Einstein Vacuum Equations on a Toroidal Manifold with Randomly Perturbed Radial Moduli: Stability Criteria and Induced 'Cosmological Constant' Terms

abstract: The Einstein vacuum equations on an (n+1)-dimensional toroidal manifold $\mathbb{M}^{n+1}=\mathbb{T}^{n}\times\mathbb{R}^{+}$ reduce to a system of n-dimensional nonlinear ODEs in terms of the set of toroidal radii $(a_{i})_{i=1}^{n}$ or the radial moduli fields $(\psi_{i})_{i=1}^{n}=(\log(a_{i}(t))_{i=1}^{n}$ of the n-torus $\mathbb{T}^{n}$. This geometry is also the basis of Kasner-Bianchi-type cosmologies. The equations are trivially satisfied for static solutions $\psi_{i}^{E}=\psi^{E}$ or radii $a_{i}^{E}=a^{E}$, describing an initially static toroidal 'micro-universe' or 'vacuum bubble'. It is Lyapunov stable to short-pulse deterministic perturbations, which have a sharp Gaussian profile: the perturbed radii rapidly converge to new attractors and therefore to new stable equilibria. These perturbations induce transitions between stable states. Introducing classical random fluctuations or perturbations, with a regulated covariance, the radial moduli become Gaussian random fields paramatrizing a 'toroidal random geometry'. The randomly perturbed Einstein equations are then interpreted as a stochastic n-dimensional nonlinear dynamical system. Non-vanishing 'cosmological constant' terms are retained within the averaged equations since they are nonlinear. This is analogous to averaging the Navier-Stokes equations in statistical turbulence theory, which yields an additional non-vanishing Reynolds term, since like the Einstein equations they are also of nonlinear hyperbolic type. The expectations of the randomly perturbed toric radii can be estimated from a cumulant expansion method. The initially static toroidal vacuum bubble undergoes eternal 'noise-induced' stochastic exponential growth or inflation. Random radial moduli fields within this scenario therefore act like a 'dark energy'. Finally, a class of random perturbations is considered for which this Einstein system is stable.


title: Overcoming the absorption limit in high-harmonic generation from crystals

abstract: Since the new millennium coherent extreme ultra-violet and soft x-ray radiation has revolutionized the understanding of dynamical physical, chemical and biological systems at the electron's natural timescale. Unfortunately, coherent laser-based upconversion of infrared photons to vacuum-ultraviolet and soft x-ray high-order harmonics in gaseous, liquid and solid targets is notoriously inefficient. In dense nonlinear media, the limiting factor is strong re-absorption of the generated high-energy photons. Here we overcome this limitation by allowing high-order harmonics generated from a periodic array of thin one-dimensional crystalline silicon ridge waveguides to propagate in the vacuum gaps between the ridges, thereby avoiding the high absorption loss of the bulk nonlinear material and resulting in a ~ 100-fold increase in propagation length. As the grating period is varied, each high-harmonic shows a different and marked modulation, indicating the onset of coherent addition which is otherwise suppressed in absorption-limited emission. By beating the absorption limit, our results pave the way for bright coherent short-wavelength sources and their implementation in nano-photonic devices.


title: LightTrack: A Generic Framework for Online Top-Down Human Pose Tracking

abstract: In this paper, we propose a novel effective light-weight framework, called LightTrack, for online human pose tracking. The proposed framework is designed to be generic for top-down pose tracking and is faster than existing online and offline methods. Single-person Pose Tracking (SPT) and Visual Object Tracking (VOT) are incorporated into one unified functioning entity, easily implemented by a replaceable single-person pose estimation module. Our framework unifies single-person pose tracking with multi-person identity association and sheds first light upon bridging keypoint tracking with object tracking. We also propose a Siamese Graph Convolution Network (SGCN) for human pose matching as a Re-ID module in our pose tracking system. In contrary to other Re-ID modules, we use a graphical representation of human joints for matching. The skeleton-based representation effectively captures human pose similarity and is computationally inexpensive. It is robust to sudden camera shift that introduces human drifting. To the best of our knowledge, this is the first paper to propose an online human pose tracking framework in a top-down fashion. The proposed framework is general enough to fit other pose estimators and candidate matching mechanisms. Our method outperforms other online methods while maintaining a much higher frame rate, and is very competitive with our offline state-of-the-art. We make the code publicly available at: https://github.com/Guanghan/lighttrack.


title: Tracking the Progression of Reading Through Eye-gaze Measurements

abstract: In this paper we consider the problem of tracking the progression of reading through eye-gaze measurements. Such an algorithm is novel and will ultimately help to develop a method of analyzing eye-gaze data which had been collected during reading activity in order to uncover crucial information regarding the individual's interest level and quality of experience while reading a passage of text or book. Additionally, such an approach will serve as a "visual signature" - a method of verifying if an individual has indeed given adequate attention to critical text-based information. Further, an accurate "reading-progression-tracker" has potential applications in educational institutions, e-readers and parenting solutions. Tracking the progression of reading remains a challenging problem due to the fact that eye-gaze movements are highly noisy and the eye-gaze is easily distracted in a limited space, like an e-book. In a prior work, we proposed an approach to analyze eye-gaze fixation points collected while reading a page of text in order to classify each measurement to a line of text; this approach did not consider tracking the progression of reading along the line of text. In this paper, we extend the capabilities of the previous algorithm in order to accurately track the progression of reading along each line. the proposed approach employs least squares batch estimation in order to estimate three states of the horizontal saccade: position, velocity and acceleration. First, the proposed approach is objectively evaluated on a simulated eye-gaze dataset. Then, the proposed algorithm is demonstrated on real data collected by a Gazepoint eye-tracker while the subject is reading several pages from an electronic book.


title: Relaxation regimes of the holographic electrons at charge neutrality after a local quench of chemical potential

abstract: In this work we study the relaxation of the system of strongly correlated electrons, at charge neutrality, when the chemical potential undergoes a local change. This setup is a model for for the X-ray absorbtion edge study in the half-filled graphene. We use holographic duality to describe the system as a classical Schwarzschild black hole in curved 4-dimensional AdS spacetime. Assuming the amplitude of the quench is small, we neglect the backreaction on the geometry. We numerically study the two relaxation regimes: the adiabatic relaxation when the quench is slow and the relaxation governed by the quasinormal modes of the system, when the quench is fast. We confirm the expectation that the scale of separation between the slow and fast regimes is set by the characteristic frequency of the quasinormal modes.


title: Toybox: A Suite of Environments for Experimental Evaluation of Deep Reinforcement Learning

abstract: Evaluation of deep reinforcement learning (RL) is inherently challenging. In particular, learned policies are largely opaque, and hypotheses about the behavior of deep RL agents are difficult to test in black-box environments. Considerable effort has gone into addressing opacity, but almost no effort has been devoted to producing high quality environments for experimental evaluation of agent behavior. We present TOYBOX, a new high-performance, open-source* subset of Atari environments re-designed for the experimental evaluation of deep RL. We show that TOYBOX enables a wide range of experiments and analyses that are impossible in other environments. *https://kdl-umass.github.io/Toybox/


title: Nanoparticles actively fragment armored droplets

abstract: Understanding the complexity of fragmentation processes is essential for regulating intercellular communication in mechanistic biology and developing novel bottom-up approaches in a large range of multiphase flow processes. In this context, self-fragmentation proceeds without any external mechanical energy input allowing one to create efficiently micro- and nanodroplets. Here we examine self-fragmentation in emulsion nanodroplets stabilized by solid particles with different surface features. Mesoscopic modelling and accelerated dynamics simulations allow us to overcome the limitations of atomistic simulations and offer detailed insight into the interplay between the evolution of the droplet shape and the particle finite-sized effects at the interface. We show that finite-sized nanoparticles play an active role in the necking breakup, behaving like nano-scale razors, and affect strongly the thermodynamic properties of the system. The role played by the particles during self-fragmentation might be of relevance to multifunctional biomaterial design and tuning of signaling pathways in mechanistic biology.


title: Mesoscale Organization and Dynamics in Binary Ionic Liquid Mixtures

abstract: The impact of mesoscale organization on dynamics and ion transport in binary ionic liquid mixtures is investigated by broadband dielectric spectroscopy, dynamic-mechanical spectroscopy, x-ray scattering, and molecular dynamics simulations. The mixtures are found to form distinct liquids with macroscopic properties that significantly deviate from weighted contributions of the neat components. For instance, it is shown that the mesoscale morphologies in ionic liquids can be tuned by mixing to enhance the static dielectric permittivity of the resulting liquid by as high as 100$\%$ relative to the neat ionic liquid components. This enhancement is attributed to the intricate role of interfacial dynamics associated with the changes in the mesoscopic aggregate morphologies in these systems. These results demonstrate the potential to design the physicochemical properties of ionic liquids through control of solvophobic aggregation


title: A SAT-based System for Consistent Query Answering

abstract: An inconsistent database is a database that violates one or more integrity constraints, such as functional dependencies. Consistent Query Answering is a rigorous and principled approach to the semantics of queries posed against inconsistent databases. The consistent answers to a query on an inconsistent database is the intersection of the answers to the query on every repair, i.e., on every consistent database that differs from the given inconsistent one in a minimal way. Computing the consistent answers of a fixed conjunctive query on a given inconsistent database can be a coNP-hard problem, even though every fixed conjunctive query is efficiently computable on a given consistent database. We designed, implemented, and evaluated CAvSAT, a SAT-based system for consistent query answering. CAvSAT leverages a set of natural reductions from the complement of consistent query answering to SAT and to Weighted MaxSAT. The system is capable of handling unions of conjunctive queries and arbitrary denial constraints, which include functional dependencies as a special case. We report results from experiments evaluating CAvSAT on both synthetic and real-world databases. These results provide evidence that a SAT-based approach can give rise to a comprehensive and scalable system for consistent query answering.


title: Experimental quantum thermodynamics with linear optics

abstract: The study of non-equilibrium physics from the perspective of the quantum limits of thermodynamics and fluctuation relations can be experimentally addressed with linear optical systems. We discuss recent experimental investigations in this scenario and present new proposed schemes and the potential advances they could bring to the field.


title: A Simple Proof of the Monotonicity of the Invariant Distribution for a Discrete Markov Chain

abstract: This note presents a simple proof of the monotonicity of the invariant distribution of a discrete Markov chain with a finite state space. This answers a question recently raised by David Siegmund.


title: Parallel-Plate Waveguides for Terahertz-Driven MeV Electron Bunch Compression

abstract: We demonstrate the electromagnetic performance of waveguides for femtosecond electron beam bunch manipulation and compression with strong-field terahertz (THz) pulses. The compressor structure is a dispersion-free exponentially-tapered parallel-plate waveguide (PPWG) that can focus single-cycle THz pulses along one dimension. We show test results of the tapered PPWG structure using electro-optic sampling (EOS) at the interaction region with peak fields of at least 300 kV/cm given 0.9 uJ of incoming THz energy. We also present a modified shorted design of the tapered PPWG for better beam manipulation and reduced magnetic field as an alternative to a dual-feed approach. As an example, we demonstrate that with 5 uJ of THz energy, the PPWG compresses a 2.5 MeV electron bunch by a compression factor of more than 4 achieving a bunch length of about 18 fs.


title: Memory Augmented Neural Network Adaptive Controllers: Performance and Stability

abstract: In this paper, we propose a novel control architecture, inspired from neuroscience, for adaptive control of continuous-time systems. The proposed architecture, in the setting of standard Neural Network (NN) based adaptive control, augments an external working memory to the NN. The controller, through a write operation, writes the hidden layer feature vector of the NN to the external working memory and can also update this information with the observed error in the output. Through a read operation, the controller retrieves information from the working memory to modify the final control signal. First, we consider a simpler estimation problem to theoretically study the effect of an external memory and prove that the estimation accuracy can be improved by incorporating memory. We then consider a model reference NN adaptive controller for linear systems with matched uncertainty to implement and illustrate our ideas. We prove that the resulting controller leads to a Uniformly Bounded (UB) stable closed loop system. Through extensive simulations and specific metrics, such as peak deviation and settling time, we show that memory augmentation improves learning significantly. Importantly, we also provide evidence for and insights on the mechanism by which this specific memory augmentation improves learning.


title: Quantization and Coorbit Spaces for Nilpotent Groups

abstract: We reconsider the quantization of symbols defined on the product between a nilpotent Lie algebra and its dual. To keep track of the non-commutative group background, the Lie algebra is endowed with the Baker-Campbell-Hausdorff product, making it via the exponential diffeomorphism a copy of its unique connected simply connected nilpotent Lie group. Using harmonic analysis tools, we emphasize the role of a Weyl system, of the associated Fourier-Wigner transformation and, at the level of symbols, of an important family of exponential functions. Such notions also serve to introduce a family of phase-space shifts. These are used to define and briefly study a new class of coorbit spaces of symbols and its relationship with coorbit spaces of vectors, defined via the Fourier-Wigner transform.


title: On the Characterization of $\tau_{(n)}$-Atoms

abstract: In 2011, Anderson and Frazier define the concept of $\tau_{(n)}$-factorization, where $\tau_{(n)}$ is a restriction of the modulo $n$ equivalence relation. These relations have been worked mostly for small values of $n$. However, it is sometimes difficult to extend findings to larger values of $n$. One of these problems is finding $\tau_{(n)}$-irreducible elements or $\tau_{(n)}$-atoms in order to characterize elements that have a $\tau_{(n)}$-factorization in $\tau_{(n)}$-atoms. The $\tau_{(n)}$-irreducible elements are well known for $n=0,1,2,3,4,5,6,8,10,12$. However, the problem of determining the $\tau_{(n)}$-atoms becomes much more difficult the larger $n$ is. In this work, we present an algorithm to construct families of $\tau_{(n)}$-atoms. It is shown that the algorithm terminates in finitely many steps when $n$ is the safe prime associated to a Sophie Germain prime.


title: Automatic Generation of Moment-Based Invariants for Prob-Solvable Loops

abstract: One of the main challenges in the analysis of probabilistic programs is to compute invariant properties that summarise loop behaviours. Automation of invariant generation is still at its infancy and most of the times targets only expected values of the program variables, which is insufficient to recover the full probabilistic program behaviour. We present a method to automatically generate moment-based invariants of a subclass of probabilistic programs, called Prob-Solvable loops, with polynomial assignments over random variables and parametrised distributions. We combine methods from symbolic summation and statistics to derive invariants as valid properties over higher-order moments, such as expected values or variances, of program variables. We successfully evaluated our work on several examples where full automation for computing higher-order moments and invariants over program variables was not yet possible.


title: The point-particle solution and the asymptotic flatness in 2+1 Horava gravity

abstract: We show that the solution corresponding to the gravitational field of a point particle at rest in 2+1 nonprojectable Horava is exactly the same as its analogous in 2+1 General Relativity. This solution is well known, it is a flat cone whose deficit angle is proportional to the mass of the particle. To establish the system we couple the Horava theory to a point particle with relativistic action. As a consequence of this result, we define the condition of asymptotic flatness exactly in the same way of 2+1 General Relativity. A remarkable feature of this condition is that the dominant mode is not fixed, but affected by the mass of the configuration. Anoter important coincidence with 2+1 General Relativity under asymptotic flatness is that the energy is the same (except for some coupling constants involved), the z = 1 term with the derivative of the lapse function does not contribute.


title: Berezin-Type Operators on the Cotangent Bundle of a Nilpotent Group

abstract: We define and study coherent states, a Berezin-Toeplitz quantization and covariant symbols on the product between a connected simply connected nilpotent Lie group and the dual of its Lie algebra. The starting point is a Weyl system codifying the natural Canonical Commutation Relations of the system. The formalism is meant to complement the quantization of the cotangent bundle by pseudo-differential operators, to which it is connected in an explicit way. Some extensions are indicated, concerning $\tau$-quantizations and variable magnetic fields.


title: Optimization Modulo the Theories of Signed Bit-Vectors and Floating-Point Numbers

abstract: Optimization Modulo Theories (OMT) is an important extension of SMT which allows for finding models that optimize given objective functions, typically consisting in linear-arithmetic or pseudo-Boolean terms. However, many SMT and OMT applications, in particular from SW and HW verification, require handling bit-precise representations of numbers, which in SMT are handled by means of the theory of Bit-Vectors (BV) for the integers and that of Floating-Point Numbers (FP) for the reals respectively. Whereas an approach for OMT with (unsigned) BV has been proposed by Nadel & Ryvchin, unfortunately we are not aware of any existing approach for OMT with FP. In this paper we fill this gap. We present a novel OMT approach, based on the novel concept of attractor and dynamic attractor, which extends the work of Nadel & Ryvchin to signed BV and, most importantly, to FP. We have implemented some OMT(BV) and OMT(FP) procedures on top of OptiMathSAT and tested the latter ones on modified problems from the SMT-LIB repository. The empirical results support the validity and feasibility of the novel approach.