title: The bispectrum of polarized galactic foregrounds

abstract: Understanding the properties of the galactic emission at millimetre wavelengths is important for studies of the cosmic microwave background (CMB). In this work we explore the bispectrum, the harmonic equivalent of the three point function, from galactic dust and synchrotron emission. We investigate these effects across a broad range of frequencies using the synchrotron dominated S-band Polarization All Sky Survey (SPASS) maps at 2.3 GHz, the Planck satellite maps (30-857 GHz) and dust dominated Infrared Astronomical Satellite (IRAS) maps at 3 THz. We measure bispectra of total intensity fields, T, as well as the gradient, E, and curl modes, B, of the polarization field. We find that the synchrotron and galactic dust have strong temperature bispectra with significant contributions in the squeezed limit, which probes the correlations between two small scale modes with a large scale mode. Additionally, we find that the dust also has strong polarised bispectra that also peak in the squeezed configuration. We explore parity odd bispectra, such as BTT bispectra, and find strong parity odd bispectra for the galactic dust notably in BTT, BTE and BEE configurations. After masking bright sources, we find no evidence for polarised synchrotron bispectra and no evidence for cross bispectra between the dust and synchrotron emission. The strong foreground bispectra discussed here need to be carefully controlled to avoid biasing measurements of primordial non-Gaussianity. Finally we use these bispectra tools to test for residual foregrounds in the component separated Planck maps and find no evidence of residual foregrounds. These tools will be useful for characterizing residual foregrounds in component separated maps, particularly for experiments with less frequency coverage than the Planck satellite.


title: New physics in $b \to s e^+ e^-$?

abstract: At present, the measurements of some observables in $B \to K^* \mu^+\mu^-$ and $B_s^0 \to \phi \mu^+ \mu^-$ decays, and of $R_{K^{(*)}} \equiv {\cal B}(B \to K^{(*)} \mu^+ \mu^-)/{\cal B}(B \to K^{(*)} e^+ e^-)$, are in disagreement with the predictions of the standard model. While most of these discrepancies can be removed with the addition of new physics (NP) in $b \to s \mu^+ \mu^-$, a difference of $>\sim 1.7 \sigma$ still remains in the measurement of $R_{K^*}$ at small values of $q^2$, the dilepton invariant mass-squared. In the context of a global fit, this is not a problem. However, it does raise the question: if the true value of $R_{K^*}^{low}$ is near its measured value, what is required to explain it? In this paper, we show that, if one includes NP in $b \to s e^+ e^-$, one can generate values for $R_{K^*}^{low}$ that are within $\sim 1\sigma$ of its measured value. Using a model-independent, effective-field-theory approach, we construct many different possible NP scenarios. We also examine specific models containing leptoquarks or a $Z'$ gauge boson. Here, additional constraints from lepton-flavour-violating observables, $B_s^0$-${\bar B}_s^0$ mixing and neutrino trident production must be taken into account, but we still find a number of viable NP scenarios. For the various scenarios, we examine the predictions for $R_{K^{(*)}}$ in other $q^2$ bins, as well as for the observable $Q_5 \equiv P^{\prime\mu\mu}_5 -P^{\prime ee}_5$.


title: Total Hadronic Cross Sections via the Holographic Pomeron Exchange

abstract: The analysis on the hadron-hadron scattering at high energies in the framework of holographic QCD is presented. Combining the Pomeron exchange kernel and the hadron-Pomeron couplings which are obtained from the bottom-up AdS/QCD models, we calculate the total cross sections. Our calculation for the nucleon-nucleon scattering is consistent with the experimental data including the recent ones measured by the TOTEM collaboration at the LHC. Within the framework, one can consider other processes involving other hadrons. As examples to show this feature of the model, the resulting total cross sections of the pion-nucleon and pion-pion scattering are also presented.


title: Decentralized Poisson Multi-Bernoulli Filtering for Vehicle Tracking

abstract: A decentralized Poisson multi-Bernoulli filter is proposed to track multiple vehicles using multiple high-resolution sensors. Independent filters estimate the vehicles' presence, state, and shape using a Gaussian process extent model; a decentralized filter is realized through fusion of the filters posterior densities. An efficient implementation is achieved by parametric state representation, utilization of single hypothesis tracks, and fusion of vehicle information based on a fusion mapping. Numerical results demonstrate the performance.


title: Collinear matching for Sivers function at next-to-leading order

abstract: We evaluate the light-cone operator product expansion for the unpolarized transverse momentum dependent (TMD) operator in the background-field technique up twist-3 inclusively. The next-to-leading order (NLO) matching coefficient for the Sivers function is derived. The method, as well as many details of the calculation, are presented.


title: astroquery: An Astronomical Web-Querying Package in Python

abstract: astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.


title: Operator Entanglement in Interacting Integrable Quantum Systems: the Case of the Rule 54 Chain

abstract: In a many-body quantum system, local operators in Heisenberg picture $O(t) = e^{i H t} O e^{-i H t}$ spread as time increases. Recent studies have attempted to find features of that spreading which could distinguish between chaotic and integrable dynamics. The operator entanglement - the entanglement entropy in operator space - is a natural candidate to provide such a distinction. Indeed, while it is believed that the operator entanglement grows linearly with time $t$ in chaotic systems, numerics suggests that it grows only logarithmically in integrable systems. That logarithmic growth has already been established for non-interacting fermions, however progress on interacting integrable systems has proved very difficult. Here, for the first time, a logarithmic upper bound is established rigorously for all local operators in such a system: the `Rule 54' qubit chain, a model of cellular automaton introduced in the 1990s [Bobenko et al., CMP 158, 127 (1993)], recently advertised as the simplest representative of interacting integrable systems. Physically, the logarithmic bound originates from the fact that the dynamics of the models is mapped onto the one of stable quasiparticles that scatter elastically; the possibility of generalizing this scenario to other interacting integrable systems is briefly discussed.

categories: cond-mat.stat-mech cond-mat.str-el math-ph math.MP quant-ph


title: Limits on electromagnetic counterparts of gravitational wave-detected binary black hole mergers

abstract: Unlike mergers of two compact objects containing a neutron star (NS), binary black hole (BBH) mergers are not accompanied by the production of tidally disrupted material, and hence lack the most direct source of accretion to power a jet and generate electromagnetic (EM) radiation. However, following a tentative detection by the Fermi GBM of a $\gamma$-ray counterpart to GW150914, several alternative ideas were proposed for driving a jet and producing EM radiation. {\em If} such jets were in fact produced, they would however lack the bright cocoon emission that makes jets from binary NSs bright also at large viewing angles. Here, via Monte Carlo simulations of a population of BBH mergers with properties consistent with those inferred from the existing LIGO/Virgo observations, and the angular emission characteristic of jets propagating into the interstellar medium, we derive limits on the allowed energetics and Lorentz factors of such jets from EM follow ups to GW-detected BBH merger events to date, and we make predictions which will help tighten these limits with broadband EM follow ups to events in future LIGO/Virgo runs. The condition that $\lesssim 1$ event out of 10 GW-detected BBH mergers be above the Fermi/GBM threshold imposes that any currently allowed emission model has to satisfy the condition $E_{\rm iso}/(10^{49} {\rm erg})(\theta_{jet}/20^\circ)\lesssim 1$.


title: Layer degradation triggers an abrupt structural transition in multiplex networks

abstract: Network robustness is a central point in network science, both from a theoretical and a practical point of view. In this paper, we show that layer degradation, understood as the continuous or discrete loss of links' weight, triggers a structural transition revealed by an abrupt change in the algebraic connectivity of the graph. Unlike traditional single layer networks, multiplex networks exist in two phases, one in which the system is protected from link failures in some of its layers and one in which all the system senses the failure happening in one single layer. We also give the exact critical value of the weight of the intra-layer links at which the transition occurs for continuous layer degradation and its relation to the value of the coupling between layers. This relation allows us to reveal the connection between the transition observed under layer degradation and the one observed under the variation of the coupling between layers.


title: Observations of Fast Radio Bursts at Frequencies down to 400 Megahertz

abstract: Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.


title: A Second Source of Repeating Fast Radio Bursts

abstract: The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.


title: Lepton Flavor Violation and Neutrino Masses from $A_5$ and CP in the Non-Universal MSSM

abstract: We analyze the phenomenological consequences of embedding a flavor symmetry based on the groups $A_5$ and CP in a supersymmetric framework. We concentrate on the leptonic sector, where two different residual symmetries are assumed to be conserved at LO for charged and neutral leptons. All possible realizations to generate neutrino masses at tree level are investigated. Sizable flavor violating effects in the charged lepton sector are unavoidable due to the non-universality of soft-breaking terms determined by the symmetry. We derive testable predictions for the neutrino spectrum, lepton mixing and flavor changing processes with non-trivial relations among observables.


title: Hopfions in lattice dimer model

abstract: In this paper, we study topological properties of 3D lattice dimer model. We demonstrate, that the dimer model on a bipartite lattice possesses topological defects, which are exactly characterized by Hopf invariant. We derive its explicit algebraic expression in terms of effective magnetic field of a dimer configuration. Thus, we solve the problem of topological classification of possible states in 3D lattice dimer model. Furthermore, since the lattice dimer model is known to be dual to spin ice, our work can be viewed as a proposal to search for hopfions in classical, as well as, artificial spin ice and related materials.


title: On monoids of ideals of orders in quadratic number fields

abstract: We determine the set of catenary degrees, the set of distances, and the unions of sets of lengths of the monoid of nonzero ideals and of the monoid of invertible ideals of orders in quadratic number fields.


title: Consistent Skyrme parametrizations constrained by GW170817

abstract: The high-density behavior of the stellar matter composed of nucleons and leptons under $\beta$~equilibrium and charge neutrality conditions is studied with the Skyrme parametrizations shown to be consistent (CSkP) with the nuclear matter, pure neutron matter, symmetry energy and its derivatives in a set of $11$ constraints [Dutra {\it et al.}, Phys. Rev. C 85, 035201 (2012)]. The predictions of these parametrizations on the tidal deformabilities related to the GW170817 event are also examined. The CSkP that produce massive neutron stars give a range of $11.86~\mbox{km} \leqslant R_{1.4} \leqslant 12.55~\mbox{km}$ for the canonical star radius, in agreement with other theoretical predictions. It is shown that the CSkP are compatible with the region of masses and radii obtained from the analysis of recent data from LIGO and Virgo Collaboration (LVC). A correlation between dimensionless tidal deformability and radius of the canonical star is found, namely, $\Lambda_{1.4} \approx 3.16\times10^{-6}R_{1.4}^{7.35}$, with results for the CSkP compatible with the recent range of $\Lambda_{1.4}=190_{-120}^{+390}$ from LVC. An analysis of the $\Lambda_1\times\Lambda_2$ graph shows that all the CSkP are compatible with the recent bounds obtained by LVC. Finally, the universal correlation between the moment of inertia and the deformability of a neutron star, named as the \mbox{$I$-Love} relation, is verified for the CSkP, that are also shown to be consistent with the prediction for the moment of inertia of the \mbox{PSR J0737-3039} primary component pulsar.


title: CrossNet: Latent Cross-Consistency for Unpaired Image Translation

abstract: Recent GAN-based architectures have been able to deliver impressive performance on the general task of image-to-image translation. In particular, it was shown that a wide variety of image translation operators may be learned from two image sets, containing images from two different domains, without establishing an explicit pairing between the images. This was made possible by introducing clever regularizers to overcome the under-constrained nature of the unpaired translation problem. In this work, we introduce a novel architecture for unpaired image translation, and explore several new regularizers enabled by it. Specifically, our architecture comprises a pair of GANs, as well as a pair of translators between their respective latent spaces. These cross-translators enable us to impose several regularizing constraints on the learnt image translation operator, collectively referred to as latent cross-consistency. Our results show that our proposed architecture and latent cross-consistency constraints are able to outperform the existing state-of-the-art on a variety of image translation tasks.


title: Statistical Models for the Number of Successful Cyber Intrusions

abstract: We propose several generalized linear models (GLMs) to predict the number of successful cyber intrusions (or "intrusions") into an organization's computer network, where the rate at which intrusions occur is a function of the following observable characteristics of the organization: (i) domain name server (DNS) traffic classified by their top-level domains (TLDs); (ii) the number of network security policy violations; and (iii) a set of predictors that we collectively call "cyber footprint" that is comprised of the number of hosts on the organization's network, the organization's similarity to educational institution behavior (SEIB), and its number of records on scholar.google.com (ROSG). In addition, we evaluate the number of intrusions to determine whether these events follow a Poisson or negative binomial (NB) probability distribution. We reveal that the NB GLM provides the best fit model for the observed count data, number of intrusions per organization, because the NB model allows the variance of the count data to exceed the mean. We also show that there are restricted and simpler NB regression models that omit selected predictors and improve the goodness-of-fit of the NB GLM for the observed data. With our model simulations, we identify certain TLDs in the DNS traffic as having significant impact on the number of intrusions. In addition, we use the models and regression results to conclude that the number of network security policy violations are consistently predictive of the number of intrusions.


title: Precise Micropatterning of a Porous Poly(Ionic Liquid) via Maskless Photolithography for High-Performance Non-Enzymatic H2O2 Sensing

abstract: Porous poly(ionic liquid)s (PILs) recently have been actively serving as a multifunctional, interdisciplinary materials platform in quite a few research areas, including separation, catalysis, actuator, sensor, and energy storage, just to name a few. In this context, the capability to photo-pattern PIL microstructures in a porous state on a substrate is still missing but is a crucial step for their real industrial usage. Here, we developed a method for in situ rapid patterning of porous PIL microstructures via a maskless photolithography approach coupled with a simple electrostatic complexation treatment. This breakthrough enables designs of miniaturized sensors. As exemplified in this work, upon loading Pt nanoparticles into porous PIL microstructures, the hybrid sensor showed outstanding performance, bearing both a high sensitivity and a wide detection range.


title: FEAST for differential eigenvalue problems

abstract: An operator analogue of the FEAST matrix eigensolver is developed to compute the discrete part of the spectrum of a differential operator in a region of interest in the complex plane. Unbounded search regions are handled with a novel rational filter for the right half-plane. If the differential operator is normal or self-adjoint, then the operator analogue preserves that structure and robustly computes eigenvalues to near machine precision accuracy. The algorithm is particularly adept at computing high-frequency modes of differential operators that possess self-adjoint structure with respect to weighted Hilbert spaces.


title: Transmission coefficients in compound-nucleus reaction theory

abstract: A recent article (D.A. Brown, et al., Phys. Rev. C98 024616 (2018)) rproposed a modification of the cross section formula used in practical calculations of compound nucleus reactions. We discuss the main concepts and approximations of statistical reaction theory and conclude that the standard practical implementations of the cross section formula remain the preferred choice.


title: A non-perturbative theory of effective Hamiltonians: example of moir\'e materials

abstract: We demonstrate that there exists a continuum Hamiltonian $H(\bf{r},\bf{p})$ that is formally the operator equivalent of the general tight-binding method, inheriting the associativity and Hermiticity of the latter operator. This provides a powerful and controlled method of obtaining effective Hamiltonians via Taylor expansion with respect to momentum and, optionally, deformation fields. In particular, for fundamentally non-perturbative defects, such as twist faults and partial dislocations, the method allows the deformation field to be retained to all orders, providing an efficient scheme for the generation of transparent and compact Hamiltonians for such defects. We apply the method to a survey of incommensurate physics in twist bilayers of graphene, graphdiyne, MoS$_2$, and phosphorene. For graphene we are able to reproduce the `reflected Dirac cones' of the $30^\circ$ quasi-crystalline bilayer found in a recent ARPES experiment, and show it is an example of a more general phenomena of coupling by the moir\'e momentum. We show that incommensurate physics is governed by the decay of the interlayer interaction on the scale of the single layer reciprocal lattices, and demonstrate that if this is slow incommensurate scattering effects lead to very rapid broadening of band manifolds as the twist angle is tuned through commensurate values.


title: All Optical Control of Beam Dynamics in a DLA

abstract: Dielectric laser acceleration draws upon nano-fabrication techniques to build photonic structures for high gradient electron acceleration. At the small spatial scales characteristic of these structures conventional accelerator techniques become ineffective at stabilizing the beam dynamics. Instead we propose a scheme to stabilize the motion by directly modulating the drive laser, in analogy to a radio-frequency-quadrupole. Here we present a design for a programmable `lattice' being built at UCLA's Pegasus laboratory. The accelerator accepts an unmodulated 3.5 MeV electron beam and then bunches and accelerates the beam by 1.5 MeV over a distance of 2 cm.


title: Two extensions of the Stone Duality to the category of zero-dimensional Hausdorff spaces

abstract: Extending the Stone Duality Theorem, we prove two duality theorems for the category ZHaus of zero-dimensional Hausdorff spaces and continuous maps. Both of them imply easily the Tarski Duality Theorem, as well as two new duality theorems for the category EDTych of extremally disconnected Tychonoff spaces and continuous maps. Also, we describe two categories which are dually equivalent to the category ZComp of zero-dimensional Hausdorff compactifications of zero-dimensional Hausdorff spaces and obtain as a corollary the Dwinger Theorem about zero-dimensional compactifications of a zero-dimensional Hausdorff space.


title: On the geometry of van Kampen diagrams of graph products of groups

abstract: In this article, we propose a geometric framework dedicated to the study of van Kampen diagrams of graph products of groups. As an application, we find information on the word and the conjugacy problems. The main new result of the paper deals with the computation of conjugacy length functions. More precisely, if $\Gamma$ is a finite graph and $\mathcal{G}= \{ G_u \mid u \in V(\Gamma) \}$ a collection of finitely generated groups indexed by the vertices of $\Gamma$, then $$\max\limits_{u \in V(\Gamma)} \mathrm{CLF}_{G_u}(n) \leq \mathrm{CLF}_{\Gamma \mathcal{G}}(n) \leq (D+1) \cdot n + \max\limits_{\Delta \subset \Gamma \ \text{complete}} \sum\limits_{u \in V(\Delta)} \mathrm{CLF}_{G_u}(n)$$ for every $n \geq 1$, where $D$ denotes the maximal diameter of a connected component of the opposite graph $\Gamma^{\mathrm{opp}}$. As a consequence, a graph product of groups with linear conjugacy length functions has linear conjugacy length function as well.


title: Index-energy estimates for Yang-Mills connections and Einstein metrics

abstract: We prove a conformally invariant estimate for the index of Schr\"odinger operators acting on vector bundles over four-manifolds, related to the classical Cwikel-Lieb-Rozenblum estimate. Applied to Yang-Mills connections we obtain a bound for the index in terms of its energy which is conformally invariant, and captures the sharp growth rate. Furthermore we derive an index estimate for Einstein metrics in terms of the topology and the Einstein-Hilbert energy. Lastly we derive conformally invariant estimates for the Betti numbers of an oriented four-manifold with positive scalar curvature.


title: Assessment of central serous chorioretinopathy (CSC) depicted on color fundus photographs using deep Learning

abstract: To investigate whether and to what extent central serous chorioretinopathy (CSC) depicted on color fundus photographs can be assessed using deep learning technology. We collected a total of 2,504 fundus images acquired on different subjects. We verified the CSC status of these images using their corresponding optical coherence tomography (OCT) images. A total of 1,329 images depicted CSC. These images were preprocessed and normalized. This resulting dataset was randomly split into three parts in the ratio of 8:1:1 respectively for training, validation, and testing purposes. We used the deep learning architecture termed InceptionV3 to train the classifier. We performed nonparametric receiver operating characteristic (ROC) analyses to assess the capability of the developed algorithm to identify CSC. The Kappa coefficient between the two raters was 0.48 (p < 0.001), while the Kappa coefficients between the computer and the two raters were 0.59 (p < 0.001) and 0.33 (p < 0.05).Our experiments showed that the computer algorithm based on deep learning can assess CSC depicted on color fundus photographs in a relatively reliable and consistent way.


title: Rare B, D and K decays, radiative and electroweak penguin decays, including constraints on $V_{td}/V_{ts}$ and $\epsilon'/\epsilon$: Summary of CKM 2018 working group 3

abstract: Rare $B$, $D$ and $K$ decays provide interesting probes of the Standard Model (SM), with a potential sensitivity to New Physics (NP) higher than other, more common, decays. Their experimental measurement is challenging, and their theoretical interpretation requires a precise knowledge of QCD at low energies, in the non-perturbative (hadronic) regime. The various $b$, $c$, $s$ decays provide different types of tests of the Standard Model, which can be performed in different experimental settings and which may exhibit correlated deviations in models of New Physics.


title: BoostNet: Bootstrapping detection of socialbots, and a case study from Guatemala

abstract: We present a method to reconstruct networks of socialbots given minimal input. Then we use Kernel Density Estimates of Botometer scores from 47,000 social networking accounts to find clusters of automated accounts, discovering over 5,000 socialbots. This statistical and data driven approach allows for inference of thresholds for socialbot detection, as illustrated in a case study we present from Guatemala.


title: All-Poly(ionic liquid) Membrane-derived Porous Carbon Membranes: Scalable Synthesis and Application for Photothermal Conversion in Seawater Desalination

abstract: Herein we firstly introduce a straightforward, scalable and technologically relevant strategy to manufacture charged porous polymer membranes (CPMs) in a controllable manner. The pore sizes and porous architectures of CPMs are well-controlled by rational choice of anions in poly(ionic liquid)s (PILs). Continuously, heteroatom-doped hierarchically porous carbon membrane (HCMs) can be readily fabricated via morphology-maintaining carbonization of as-prepared CPMs. These HCMs being as photothermal membranes exhibited excellent performance for solar seawater desalination, representing a promising strategy to construct advanced functional nanomaterials for portable water production technologies.


title: PointWise: An Unsupervised Point-wise Feature Learning Network

abstract: We present a novel approach to learning a point-wise, meaningful embedding for point-clouds in an unsupervised manner, through the use of neural-networks. The domain of point-cloud processing via neural-networks is rapidly evolving, with novel architectures and applications frequently emerging. Within this field of research, the availability and plethora of unlabeled point-clouds as well as their possible applications make finding ways of characterizing this type of data appealing. Though significant advancement was achieved in the realm of unsupervised learning, its adaptation to the point-cloud representation is not trivial. Previous research focuses on the embedding of entire point-clouds representing an object in a meaningful manner. We present a deep learning framework to learn point-wise description from a set of shapes without supervision. Our approach leverages self-supervision to define a relevant loss function to learn rich per-point features. We train a neural-network with objectives based on context derived directly from the raw data, with no added annotation. We use local structures of point-clouds to incorporate geometric information into each point's latent representation. In addition to using local geometric information, we encourage adjacent points to have similar representations and vice-versa, creating a smoother, more descriptive representation. We demonstrate the ability of our method to capture meaningful point-wise features through three applications. By clustering the learned embedding space, we perform unsupervised part-segmentation on point clouds. By calculating euclidean distance in the latent space we derive semantic point-analogies. Finally, by retrieving nearest-neighbors in our learned latent space we present meaningful point-correspondence within and among point-clouds.


title: Model Checking Clinical Decision Support Systems Using SMT

abstract: Individual clinical Knowledge Artifacts (KA) are designed to be used in Clinical Decision Support (CDS) systems at the point of care for delivery of safe, evidence-based care in modern healthcare systems. For formal authoring of a KA, syntax verification and validation is guaranteed by the grammar. However, there are no methods for semantic verification. Any semantic fallacy may lead to rejection of the outcomes by care providers. As a first step toward solving this problem, we present a framework for translating the logical segments of KAs into Satisfiability Modulo Theory (SMT) models. We present the effectiveness and efficiency of our work by automatically translating the logic fragment of publicly available KAs and verifying them using Z3 SMT solver.


title: de Sitter versus Anti de Sitter flows and the (super)gravity landscape

abstract: Generic solutions are studied in Einstein-scalar gravity in an ansatz that can interpolate between de Sitter and Anti-de Sitter regimes. The classification of regular solutions of \cite{exotic} is first extended to the dS regime. This implies, among others, the existence of cosmic clocks that reverse direction without passing through a curvature singularity. We then consider an ansatz for solutions that interpolate between the dS and AdS regimes. The structure of such more general solutions and their singularities are studied. It is shown that there are no regular solutions that interpolate between dS and AdS extrema for generic potentials. This is unlike the Centaur solutions that were shown to exist in two bulk dimensions. We also comment on the potential interplay with recent dS conjectures and the dS BF bounds.


title: Self-Supervised Deep Active Accelerated MRI

abstract: We propose to simultaneously learn to sample and reconstruct magnetic resonance images (MRI) to maximize the reconstruction quality given a limited sample budget, in a self-supervised setup. Unlike existing deep methods that focus only on reconstructing given data, thus being passive, we go beyond the current state of the art by considering both the data acquisition and the reconstruction process within a single deep-learning framework. As our network learns to acquire data, the network is active in nature. In order to do so, we simultaneously train two neural networks, one dedicated to reconstruction and the other to progressive sampling, each with an automatically generated supervision signal that links them together. The two supervision signals are created through Monte Carlo tree search (MCTS). MCTS returns a better sampling pattern than what the current sampling network can give and, thus, a better final reconstruction. The sampling network is trained to mimic the MCTS results using the previous sampling network, thus being enhanced. The reconstruction network is trained to give the highest reconstruction quality, given the MCTS sampling pattern. Through this framework, we are able to train the two networks without providing any direct supervision on sampling.


title: A cryostat to hold frozen-spin polarized HD targets in CLAS: HDice-II

abstract: The design, fabrication, operation, and performance of a helium-3/4 dilution refrigerator and superconducting magnet system for holding a frozen-spin polarized hydrogen deuteride target in the Jefferson Laboratory CLAS detector during photon beam running is reported. The device operates both vertically (for target loading) and horizontally (for target bombardment). The device proves capable of maintaining a base temperature of 50 mK and a holding field of 1 Tesla for extended periods. These characteristics enabled multi-month polarization lifetimes for frozen spin HD targets having proton polarization of up to 50% and deuteron up to 27%.


title: Walking through the Gaussian Primes

abstract: The Gaussian Moat problem asks whether one can walk to infinity in the Gaussian integers using the Gaussian primes as stepping stones and taking bounded length steps or not. In this paper, we have analyzed the Gaussian primes and also developed an algorithm to find the primes on the $\mathbb{R}^2$ plane which will help us to calculate the moat for higher value. We have also reduced a lot of computation with this algorithm to find the Gaussian prime though their distribution on the $\mathbb{R}^2$ plane is not so regular. A moat of value $\sqrt{26}$ is already an existing result done by Genther et.al. The focus of the problem is to show that primes are getting lesser as we are approaching infinity. We have shown this result with the help of our algorithm. We have calculated the moat and also calculated the time complexity of our algorithm and compared it with Genther-Wagon-Wick's algorithm. As a new ingredient, we have defined the notion of primality for the plane $\mathbb{R}^3$ and proposed a problem on it.


title: Magnetic dielectric- graphene- ferroelectric system as a promising non-volatile device for modern spintronics

abstract: The conductivity of the system magnetic dielectric (EuO) - graphene channel - ferroelectric substrate was considered. The magnetic dielectric locally transforms the band spectrum of graphene by inducing an energy gap in it and making it spin-asymmetric with respect to the free electrons. The range of spontaneous polarization (2- 5)mC/m2 that can be easily realized in thin films of proper and incipient ferroelectrics, was under examination. It was demonstrated, that if the Fermi level in the graphene channel belongs to energy intervals where the graphene band spectrum, modified by EuO, becomes sharply spin-asymmetric, such a device can be an ideal non-volatile spin filter. Controlling of the Fermi level (e.g. by temperature that changes ferroelectric polarization) can convert a spin filter to a spin valve.


title: Quantum Computing with sine-Gordon Qubits

abstract: A universal quantum computing scheme, with a universal set of logical gates, is proposed based on networks of 1D quantum systems. The encoding of information is in terms of universal features of gapped phases, for which effective field theories such as sine-Gordon field theory can be employed to describe a qubit. Primary logical gates are from twist, pump, glue, and shuffle operations that can be realized in principle by tuning parameters of the systems. Our scheme demonstrates the power of 1D quantum systems for robust quantum computing.


title: Submanifold geometry in symmetric spaces of noncompact type

abstract: In this survey article we provide an introduction to submanifold geometry in symmetric spaces of noncompact type. We focus on the construction of examples and the classification problems of homogeneous and isoparametric hypersurfaces, polar and hyperpolar actions, and homogeneous CPC submanifolds.


title: Reactor neutrino oscillations as constraints on Effective Field Theory

abstract: We study constraints on the Standard Model Effective Field Theory (SMEFT) from neutrino oscillations in short-baseline reactor experiments. We calculate the survival probability of reactors antineutrinos at the leading order in the SMEFT expansion, that is including linear effects of dimension-6 operators. It is shown that, at this order, reactor experiments alone cannot probe charged-current contact interactions between leptons and quarks that are of the (pseudo)vector (V$\pm$A) or pseudo-scalar type. We also note that flavor-diagonal (pseudo)vector coefficients do not have observable effects in oscillation experiments. In this we reach novel or different conclusions than prior analyses of non-standard neutrino interactions. On the other hand, reactor experiments offer a unique opportunity to probe tensor and scalar SMEFT operators that are off-diagonal in the lepton-flavor space. We derive constraints on the corresponding SMEFT parameters using the most recent data from the Daya Bay and RENO experiments.


title: De Sitter Space and Entanglement

abstract: We argue that the notion of entanglement in de Sitter space arises naturally from the non-trivial Lorentzian geometry of the spacetime manifold, which consists of two disconnected boundaries and a causally disconnected interior. In four bulk dimensions, we propose an holographic description of an inertial observer in terms of a thermofield double state in the tensor product of the two boundaries Hilbert spaces, whereby the Gibbons--Hawking formula arises as the holographic entanglement entropy between the past and future conformal infinities. When considering the bulk entanglement between the two causally disconnected Rindler wedges, we show that the corresponding entanglement entropy is given by one quarter of the area of the pair of codimension two minimal surfaces that define the set of fixed points of the dS$_4/\mathbb Z_q$ orbifold.


title: Music Artist Classification with Convolutional Recurrent Neural Networks

abstract: Previous attempts at music artist classification use frame level audio features which summarize frequency content within short intervals of time. Comparatively, more recent music information retrieval tasks take advantage of temporal structure in audio spectrograms using deep convolutional and recurrent models. This paper revisits artist classification with this new framework and empirically explores the impacts of incorporating temporal structure in the feature representation. To this end, an established classification architecture, a Convolutional Recurrent Neural Network (CRNN), is applied to the artist20 music artist identification dataset under a comprehensive set of conditions. These include audio clip length, which is a novel contribution in this work, and previously identified considerations such as dataset split and feature level. Our results improve upon baseline works, verify the influence of the producer effect on classification performance and demonstrate the trade-offs between audio length and training set size. The best performing model achieves an average F1 score of 0.937 across three independent trials which is a substantial improvement over the corresponding baseline under similar conditions. Additionally, to showcase the effectiveness of the CRNN's feature extraction capabilities, we visualize audio samples at the model's bottleneck layer demonstrating that learned representations segment into clusters belonging to their respective artists.


title: Thermal corrections to quantum friction and decoherence: a Closed-Time-Path approach to atom-surface interaction

abstract: In this paper we study the dissipative effects and decoherence induced on a particle moving at constant speed in front of a dielectric plate in quantum vacuum, developing a Closed-Time-Path (CTP) integral formulation in order to account for the corrections to these phenomena generated by finite temperatures. We compute the frictional force of the moving particle and find that it contains two different contributions: a pure quantum term due to quantum fluctuations (even present at vanishing temperatures) and a temperature dependent component generated by thermal fluctuations (bigger contribution the higher the temperature). We further estimate the decoherence timescale for the internal degree of freedom of the quantum particle. As expected, decoherence time is reduced by temperature, however, this feature is stronger for large velocities and for resonant situations. When the particle approaches relativistic speed, decoherence time becomes independent of temperature. The finite temperatures corrections to the force or even in the decoherence timescale could be used to track traces of quantum friction through the study of the velocity dependence since the solely evidence of this dependence provides an indirect testimony of the existence of a quantum frictional force


title: The Solar Benchmark: Rotational Modulation of the Sun Reconstructed from Archival Sunspot Records

abstract: We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online.


title: K2-291 b: A rocky super-Earth in a 2.2 day orbit

abstract: K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.


title: Identifying Sagittarius Stream Stars By Their APOGEE Chemical Abundance Signatures

abstract: The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey provides precise chemical abundances of 18 chemical elements for $\sim$ 176,000 red giant stars distributed over much of the Milky Way Galaxy (MW), and includes observations of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). The APOGEE chemical abundance patterns of Sgr have revealed that it is chemically distinct from the MW in most chemical elements. We employ a \emph{k}-means clustering algorithm to 6-dimensional chemical space defined by [(C+N)/Fe], [O/Fe], [Mg/Fe], [Al/Fe], [Mn/Fe], and [Ni/Fe] to identify 62 MW stars in the APOGEE sample that have Sgr-like chemical abundances. Of the 62 stars, 35 have \emph{Gaia} kinematics and positions consistent with those predicted by \emph{N}-body simulations of the Sgr stream, and are likely stars that have been stripped from Sgr during the last two pericenter passages ($<$ 2 Gyr ago). Another 20 of the 62 stars exhibit chemical abundances indistinguishable from the Sgr stream stars, but are on highly eccentric orbits with median $r_{\rm apo} \sim $ 25 kpc. These stars are likely the `accreted' halo population thought to be the result of a separate merger with the MW 8-11 Gyr ago. We also find one hypervelocity star candidate. We conclude that Sgr was enriched to [Fe/H] $\sim$ -0.2 before its most recent pericenter passage. If the `accreted halo' population is from one major accretion event, then this progenitor galaxy was enriched to at least [Fe/H] $\sim$ -0.6, and had a similar star formation history to Sgr before merging.


title: Minimally Connected Hypergraphs

abstract: Graphs and hypergraphs are foundational structures in discrete mathematics. They have many practical applications, including the rapidly developing field of bioinformatics, and more generally, biomathematics. They are also a source of interesting algorithmic problems. In this paper, we define a \textit{construction process} for minimally connected $r$-uniform hypergraphs, which captures the intuitive notion of building a hypergraph piece-by-piece, and a numerical invariant called the \textit{tightness}, which is independent of the construction process used. Using these tools, we prove some fundamental properties of minimally connected hypergraphs. We also give bounds on their chromatic numbers and provide some results involving edge colorings. We show that every connected $r$-uniform hypergraph contains a minimally connected spanning subhypergraph and provide a polynomial-time algorithm for identifying such a subhypergraph.


title: Chaos and Scrambling in Quantum Small Worlds

abstract: Quantum small-worlds are quantum many-body systems that interpolate between completely ordered (nearest-neighbour, next-to-nearest-neighbour etc.) and completely random interactions. As such, they furnish a novel new laboratory to study quantum systems transitioning between regular and chaotic behaviour. In this article, we introduce the idea of a quantum small-world network by starting from a well understood integrable system, a spin-1 Heisenberg chain. We then inject a small number of long-range interactions into the spin chain and study its ability to scramble quantum information using two primary devices: the out-of-time-order correlator (OTOC) and the spectral form factor (SFF). We find that the system shows increasingly rapid scrambling as its interactions become progressively more random, with no evidence of quantum chaos.


title: Putting Fairness Principles into Practice: Challenges, Metrics, and Improvements

abstract: As more researchers have become aware of and passionate about algorithmic fairness, there has been an explosion in papers laying out new metrics, suggesting algorithms to address issues, and calling attention to issues in existing applications of machine learning. This research has greatly expanded our understanding of the concerns and challenges in deploying machine learning, but there has been much less work in seeing how the rubber meets the road. In this paper we provide a case-study on the application of fairness in machine learning research to a production classification system, and offer new insights in how to measure and address algorithmic fairness issues. We discuss open questions in implementing equality of opportunity and describe our fairness metric, conditional equality, that takes into account distributional differences. Further, we provide a new approach to improve on the fairness metric during model training and demonstrate its efficacy in improving performance for a real-world product


title: Special class of self-similar dynamics for Rayleigh-Taylor mixing with variable acceleration

abstract: Rayleigh-Taylor (RT) mixing has critical importance for a broad range of process in nature and technology, from supernovae and plasma fusion to oil recovery and nano-fabrication. In most instances, RT flows are driven by variable acceleration, whereas the bulk of existing studies have considered only constant and impulsive accelerations. By analyzing symmetries of RT dynamics for certain patterns of variable acceleration, we discover a special class of self-similar solutions and identify their scaling, correlations and spectra. We find that dynamics of RT mixing can vary from super-ballistic to subdiffusive depending on the acceleration and retain memory of deterministic conditions for any acceleration. These rich dynamic properties considerably impact the understanding and control of RT relevant phenomena in fluids, plasmas, materials, and reveal conditions at which turbulence can be realized in RT mixing.


title: Breaking the Radio - Gamma-Ray Connection in Arp 220

abstract: Recent analyses of the gamma-ray spectrum from the ultra-luminous infrared galaxy Arp 220 have revealed a discrepancy in the cosmic ray energy injection rates derived from the gamma-rays versus the radio emission. While the observed radio emission is consistent with the star formation rate inferred from infrared observations, a significantly higher cosmic ray population is necessary to accurately model the measured gamma-ray flux. To resolve this discrepancy between the radio and gamma-ray observations, we find that we must increase the cosmic ray energy injection rate and account for an infrared optical depth greater than unity. Raising the energy injection rate naturally raises the total gamma-ray flux but also raises the radio flux unless there is also an increase in the energy loss rate for cosmic ray leptons. A optically thick medium results in an increase in energy losses via inverse Compton for cosmic ray leptons and preserves agreement with submillimeter, millimeter, and infrared wavelength observations.


title: Gas Jet Morphology and the Very Rapidly Increasing Rotation Period of Comet 41P/Tuttle-Giacobini-Kresak

abstract: We present results from our 47-night imaging campaign of Comet 41P/Tuttle-Giacobini-Kresak conducted from Lowell Observatory between 2017 February 16 and July 2. Coma morphology revealed gas jets, whose appearance and motion as a function of time yielded the rotation period and other properties. All narrowband CN images exhibited either one or two jets; one jet appeared as a partial face-on spiral with clockwise rotation while the second jet evolved from a side-on corkscrew, through face-on, and finally corkscrew again, with only a slow evolution throughout the apparition due to progressive viewing geometry changes. A total of 78 period determinations were made over a 7-week interval, yielding a smooth and accelerating rotation period starting at 24 hr (March 21&22) and passing 48 hr on April 28. While this is by far the fastest rate of change ever measured for a comet nucleus, the torque required is readily within what can exist given likely properties of the nucleus. If the torque remained constant, we estimate that the nucleus could have stopped rotating and/or began to tumble as soon as only two months following perihelion, and will certainly reach this stage by early in the next apparition. Working backwards in time, Tuttle-Giacobini-Kresak would have been rotating near its rotational break-up velocity 3-4 orbits earlier, suggesting that its extreme 7-magnitude outburst observed in 2001 might have been caused by a partial fragmentation at that time, as might the pair of 1973 8-magnitude outbursts if there had been an earlier spin-down and spin-up cycle.


title: Optimal Power Flow with Step-Voltage Regulators in Multi-Phase Distribution Networks

abstract: This paper develops a branch-flow based optimal power flow (OPF) problem for multi-phase distribution networks that allows for tap selection of wye, closed-delta, and open-delta step-voltage regulators (SVRs). SVRs are assumed ideal and their taps are represented by continuous decision variables. To tackle the non-linearity, the branch-flow semidefinite programming framework of traditional OPF is expanded to accommodate SVR edges. Three types of non-convexity are addressed: (a) rank-1 constraints on non-SVR edges, (b) nonlinear equality constraints on SVR power flows and taps, and (c) trilinear equalities on SVR voltages and taps. Leveraging a practical phase-separation assumption on the SVR secondary voltage, novel McCormick relaxations are provided for (c) and certain rank-1 constraints of (a), while dropping the rest. A linear relaxation based on conservation of power is used in place of (b). Numerical simulations on standard distribution test feeders corroborate the merits of the proposed convex formulation.


title: Limiting Top-Higgs Interaction and Higgs-Boson Width from Multi-Top Productions

abstract: We demonstrate that the multi-top productions efficiently probe the CP-property of top-Higgs interaction and the Higgs-boson width at the LHC. The four top-quark production alone can exclude a purely CP-odd top-quark Yukawa coupling at the 13~TeV LHC with an integrated luminosity of $430~{\rm fb}^{-1}$, regardless the size of the Yukawa coupling. Combining the single Higgs-boson production, the $t\bar{t}H$ associated production and the four top-quark production, we show that the CP-phase of the top-quark Yukawa coupling and the Higgs-boson width can be stringently bounded at the LHC with integrated luminosities of $300~{\rm fb}^{-1}$ and $3000~{\rm fb}^{-1}$.


title: Geometry of Kantorovich Polytopes and Support of Optimizers for Repulsive Multi-Marginal Optimal Transport on Finite State Spaces

abstract: We consider symmetric multi-marginal Kantorovich optimal transport problems on finite state spaces with uniform-marginal constraint. These problems consist of minimizing a linear objective function over a high-dimensional polytope, here referred to as Kantorovich polytope. The presented results are of split nature, computational and theoretical. Within the computational part only small numbers of marginals $N$ and marginal sites $\ell$ are considered. This restriction allows us to computationally determine all extreme points of the Kantorovich polytope and investigate how many of them are in compliance with the in optimal transport typical Monge ansatz. Singling out the results for $\ell=3$ discretization points and pairwise symmetric cost functions enables us to visually compare Kantorovich's to Monge's ansatz space for a varying number of marginals. Finally we present a necessary support-condition for optimizers which is inspired by the insights the said model problem on three sites provided. This result is not limited to the case of $\ell=3$ sites and applies to symmetric pair-costs whose diagonal entries lie above a cost-specific threshold. In case $N$ and $\ell$ display certain relationships the discussed condition provides an optimizer in Monge-form and implies its uniqueness as a solution of the considered Kantorovich optimal transport problem.


abstract: We consider the long-standing problem of Rayleigh-Taylor instability with variable acceleration, and focus on the early-time dynamics of an interface separating incompressible ideal fluids of different densities subject to an acceleration being a power-law function of time for a spatially extended threedimensional flow periodic in the plane normal to the acceleration with symmetry group p6mm. By employing group theory and scaling analysis, we discover two distinct sub-regimes of the early time dynamics depending on the exponent of the acceleration power-law. The time-scale and the early-time dynamics are set by the acceleration for exponents greater than -2, and by the initial growth-rate (due to, e.g., initial conditions) for exponents smaller than -2. At the exponent value (-2) a transition occurs from one regime to the other with varying acceleration strength. For a broad range of the acceleration parameters, the instability growth-rate is explicitly found, the dependence of the dynamics on the initial conditions is investigated, and theory benchmarks are elaborated.


title: Single-Shot Readout Performance of Two Heterojunction-Bipolar-Transistor Amplification Circuits at Millikelvin Temperatures

abstract: High-fidelity single-shot readout of spin qubits requires distinguishing states much faster than the T1 time of the spin state. One approach to improving readout fidelity and bandwidth (BW) is cryogenic amplification, where the signal from the qubit is amplified before noise sources are introduced and room-temperature amplifiers can operate at lower gain and higher BW. We compare the performance of two cryogenic amplification circuits: a current-biased heterojunction bipolar transistor circuit (CB-HBT), and an AC-coupled HBT circuit (AC-HBT). Both circuits are mounted on the mixing-chamber stage of a dilution refrigerator and are connected to silicon metal oxide semiconductor (Si-MOS) quantum dot devices on a printed circuit board (PCB). The power dissipated by the CB-HBT ranges from 0.1 to 1 {\mu}W whereas the power of the AC-HBT ranges from 1 to 20 {\mu}W. Referred to the input, the noise spectral density is low for both circuits, in the 15 to 30 fA/$\sqrt{\textrm{Hz}}$ range. The charge sensitivity for the CB-HBT and AC-HBT is 330 {\mu}e/$\sqrt{\textrm{Hz}}$ and 400 {\mu}e/$\sqrt{\textrm{Hz}}$, respectively. For the single-shot readout performed, less than 10 {\mu}s is required for both circuits to achieve bit error rates below $10^{-3}$, which is a putative threshold for quantum error correction.


title: Real-Time Predictive Control Strategy Optimization

abstract: Traffic congestion has lead to an increasing emphasis on management measures for a more efficient utilization of existing infrastructure. In this context, this paper proposes a novel framework that integrates real-time optimization of control strategies (tolls, ramp metering rates, etc.) with guidance generation using predicted network states for Dynamic Traffic Assignment systems. The efficacy of the framework is demonstrated through a fixed demand dynamic toll optimization problem which is formulated as a non-linear program to minimize predicted network travel times. A scalable efficient genetic algorithm is applied to solve this problem that exploits parallel computing. Experiments using a closed-loop approach are conducted on a large scale road network in Singapore to investigate the performance of the proposed methodology. The results indicate significant improvements in network wide travel time of up to 9% with real-time computational performance.


title: On the fundamentals of Rayleigh-Taylor dynamics with variable acceleration

abstract: Rayleigh-Taylor instability (RTI) has critical importance for a broad range of processes in nature and technology, from supernovae to plasma fusion. In most instances RTI is driven by variable acceleration whereas the bulk of existing studies have considered constant acceleration. This work focuses on RTI driven by acceleration with power-law time-dependence, and applies group theory to solve the classical problem. For early time dynamics, we find dependence of RTI growth-rate on acceleration parameters and initial conditions. For late time dynamics, we directly link interface dynamics to interfacial shear, find continuous family of regular asymptotic solutions, and discover invariance properties of nonlinear RTI. Our results reveal the interfacial and multi-scale character of RTI with variable acceleration. The former is exhibited in structure of flow fields with intense fluid motion near the interface and effectively no motion in the bulk; the latter follows from the invariance properties of nonlinear dynamics defined by the interplay of two macroscopic length-scales - the wavelength and the amplitude. Our theory resolves the long-standing problem of RTI nonlinear dynamics, achieves excellent agreement with observations, and elaborates diagnostic benchmarks for future experiments and simulations.


title: Rejuvenated galaxies with very old bulges at the origin of the bending of the main sequence and of the "green valley"

abstract: We investigate the nature of star-forming galaxies with reduced specific star formation rate (sSFR) and high stellar masses, those `green valley' objects that seemingly cause a reported bending, or flattening, of the star-forming main sequence. The fact that such objects host large bulges recently led some to suggest that the internal formation of bulges was a late event that induced the sSFRs of massive galaxies to drop in a slow downfall, and thus the main sequence to bend. We have studied in detail a sample of 10 galaxies at $0.45<z<1$ with secure SFR from Herschel, deep Keck optical spectroscopy, and HST imaging from CANDELS allowing us to perform multi-wavelength bulge to disc decomposition, and to derive star formation histories for the separated bulge and disc components. We find that the bulges hosted in these systems below main sequence are virtually all maximally old, with ages approaching the age of the Universe at the time of observation, while discs are young ($\langle$ T$_{50}\rangle \sim 1.5$ Gyr). We conclude that, at least based on our sample, the bending of the main sequence is, for a major part, due to rejuvenation, and we disfavour mechanisms that postulate the internal formation of bulges at late times. The very old stellar ages of our bulges suggest a number density of Early Type Galaxies at $z=1-3$ higher than actually observed. If confirmed, this might represent one of the first direct validations of hierarchical assembly of bulges at high redshifts.


title: The Adaptive LQ Regulator

abstract: The optimal adaptive control of a linear system in a signal-plus-noise setting with infinite horizon LQ regulator cost is studied. The class of partially observed linear systems for which the certainty equivalence property holds is identified. It is also shown that a linear system is adaptively stabilizable if and only if it is uniformly stabilizable, and the class of partially observed linear systems for which the certainty equivalence value function is a supersolution of the Bellman equation is identified.


title: Stability and flow fields structure for interfacial dynamics with interfacial mass flux

abstract: We analyze from a far field the evolution of an interface that separates ideal incompressible fluids of different densities and has an interfacial mass flux. We develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations, and analyze the interplay of interface stability with flow fields structure, by directly linking rigorous mathematical attributes to physical observables. New mechanisms are identified of the interface stabilization and destabilization. We find that interfacial dynamics is stable when it conserves the fluxes of mass, momentum and energy. The stabilization is due to inertial effects causing small oscillations of the interface velocity. In the classic Landau dynamics, the postulate of perfect constancy of the interface velocity leads to the development of the Landau-Darrieus instability. This destabilization is also associated with the imbalance of the perturbed energy at the interface, in full consistency with the classic results. We identify extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including formal properties of fundamental solutions and qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.


title: On geometric complexity theory: Multiplicity obstructions are stronger than occurrence obstructions

abstract: Geometric Complexity Theory as initiated by Mulmuley and Sohoni in two papers (SIAM J Comput 2001, 2008) aims to separate algebraic complexity classes via representation theoretic multiplicities in coordinate rings of specific group varieties. The papers also conjecture that the vanishing behavior of these multiplicities would be sufficient to separate complexity classes (so-called occurrence obstructions). The existence of such strong occurrence obstructions has been recently disproven in 2016 in two successive papers, Ikenmeyer-Panova (Adv. Math.) and B\"urgisser-Ikenmeyer-Panova (J. AMS). This raises the question whether separating group varieties via representation theoretic multiplicities is stronger than separating them via occurrences. This paper provides for the first time a setting where separating with multiplicities can be achieved, while the separation with occurrences is provably impossible. Our setting is surprisingly simple and natural: We study the variety of products of homogeneous linear forms (the so-called Chow variety) and the variety of polynomials of bounded border Waring rank (i.e. a higher secant variety of the Veronese variety). As a side result we prove a slight generalization of Hermite's reciprocity theorem, which proves Foulkes' conjecture for a new infinite family of cases.


title: Definable coaisles over rings of weak global dimension at most one

abstract: In the setting of the unbounded derived category D(R) of a ring R of weak global dimension at most one we consider t-structures with a definable coaisle. The t-structures among these which are stable (that is, the t-structures which consist of a pair of triangulated subcategories) are precisely the ones associated to a smashing localization of the derived category. In this way, our present results generalize those of [B\v{S}17] to the non-stable case. As in the stable case [B\v{S}17], we confine for the most part to the commutative setting, and give a full classification of definable coaisles in the local case, that is, over valuation domains. It turns out that unlike in the stable case of smashing subcategories, the definable coaisles do not always arise from homological ring epimorphisms. We also consider a non-stable version of the telescope conjecture for t-structures and give a ring-theoretic characterization of the commutative rings of weak global dimension at most one for which it is satisfied.


title: On Rayleigh-Taylor interfacial mixing

abstract: The Rayleigh-Taylor instability develops when fluids are accelerated counter to their density gradients; intense interfacial fluid mixing ensues with time. The Rayleigh-Taylor mixing controls a broad range of processes in fluids, plasmas, materials, at astrophysical and at atomic scales. In this perspective paper we briefly review theoretical and experimental approaches, and apply theory and experiment to investigate order and disorder in Rayleigh-Taylor flows. The theory finds that properties of heterogeneous anisotropic accelerated Rayleigh-Taylor mixing depart from those of homogeneous isotropic inertial turbulence, including strong correlations, weak fluctuations, and sensitivity to deterministic conditions. The experiment unambiguously observes the heterogeneity and anisotropy of Rayleigh-Taylor mixing at very high Reynolds numbers, and the stabilizing effects of acceleration and accelerated shear on the interfacial dynamics. The theory and the experiment agree with one another and evince that Rayleigh- Taylor mixing may exhibit order and laminarize, similarly to other accelerated flows, thus opening new perspectives for research of complex processes in nature and technology.


title: Experimental Evidence about "A factorisation algorithm in adiabatic quantum computation" by T. D. Kieu

abstract: Computations show that the logic about a quantum factoring algorithm does not hold in reality on a D-Wave quantum computer. We demonstrate this for the integers 15 = 3 x 5, 91 = 7 x 13 and 899 = 29 x 31. The likely cause is the D-Wave hardware that does not accept input terms that are a number, i.e., only terms that contain a Boolean variable can be an input. Without terms that are numbers, the relative magnitude of the coefficients in the factoring algorithm is too great to differentiate values.


title: Impact of ICM disturbances on the mean pressure profile of galaxy clusters: a prospective study of the NIKA2 SZ large program with MUSIC synthetic clusters

abstract: $Context$. The mean pressure profile of the galaxy cluster population plays an essential role in cosmological analyses. An accurate characterization of the shape, intrinsic scatter, and redshift evolution of this profile is necessary to estimate some of the biases and systematic effects that currently prevent cosmological analyses based on thermal Sunyaev-Zel'dovich (tSZ) surveys from obtaining precise and unbiased cosmological constraints. This is one of the main goals of the ongoing NIKA2 SZ large program, which aims at mapping the tSZ signal of a representative cluster sample selected from the $Planck$ and ACT catalogs and spans a redshift range $0.5 < z < 0.9$. $Aims$. To estimate the impact of intracluster medium (ICM) disturbances that can be detected by NIKA2 on the mean pressure profile of galaxy clusters, we realized a study based on a synthetic cluster sample that is similar to that of the NIKA2 SZ large program. $Methods$. To reach this goal we employed the hydrodynamical simulation $Marenostrum~MUltidark~SImulations~of~galaxy~Clusters$ (MUSIC). We simulated realistic NIKA2 and $Planck$ tSZ observations, which were jointly analyzed to estimate the ICM pressure profile of each cluster. A comparison of the deprojected profiles with the true radial profiles directly extracted from the MUSIC simulation allowed us to validate the NIKA2 tSZ pipeline and to study the impact of ICM disturbances on the characterization of the ICM pressure distribution even at high redshift. After normalizing each profile by the integrated quantities estimated under the hydrostatic equilibrium hypothesis, we evaluated the mean pressure profile of the twin sample and show that it is compatible with that extracted directly from the MUSIC simulation in the scale range that can be recovered by NIKA2. We studied the impact of cluster dynamical state on both its shape and associated scatter. (abridged)


title: Covariant polarized radiative transfer on cosmological scales for investigating large-scale magnetic field structures

abstract: Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covariant. It accounts for cosmological and relativistic effects in a self-consistent manner and explicitly treats Faraday rotation, as well as Faraday conversion, emission, and absorption processes. The formulation is derived from the first principles of conservation of phase-space volume and photon number. Without loss of generality, we consider a flat Friedmann-Robertson-Walker (FRW) space-time metric and construct the corresponding polarized radiative transfer equations. We propose an all-sky CPRT calculation algorithm, based on a ray-tracing method, which allows cosmological simulation results to be incorporated and, thereby, model templates of polarization maps to be constructed. Such maps will be crucial in our interpretation of polarized data, such as those to be collected by the Square Kilometer Array (SKA). We describe several tests which are used for verifying the code and demonstrate applications in the study of the polarization signatures in different distributions of electron number density and magnetic fields. We present a pencil-beam CPRT calculation and an all-sky calculation, using a simulated galaxy cluster or a model magnetized universe obtained from GCMHD+ simulations as the respective input structures. The implications on large-scale magnetic field studies are discussed; remarks on the standard methods using rotation measure are highlighted.


title: Revisiting the dynamics of planets in binaries: evolutionary timescales and the effect of early stellar evolution

abstract: The discovery of planets in binaries is one the most interesting outcomes of planetary research. With the growing number of discoveries has also grown the interest on describing their formation, long-term evolution and potential habitability. In this work we revisit the dynamics of planets in S-type binary systems. For that purpose we develop explicit formulas for the secularized octupolar Hamiltonian, coupled with general relativistic corrections and non-conservative interactions. We implemented those formulas in an open-source package \texttt{SecDev3B}, that can be used to reproduce our results or test improved versions of the models. In order to test it, we study the long-term dynamical evolution of S-type binary planets during the pre-main-sequence phase of stellar evolution. During that phase, stellar radius significantly changes in timescales similar to secular timescales. We hypothesize that when close-encounters between the planet and its host star happens (e.g. via Lidov-Kozai effect), particularities in the secular formalism plus changes in stellar radius may alter significantly the dynamical evolution. We study the well-known binary planet HD 80606b and found that an octupolar expansion of the conservative Hamiltonian is required to properly predict its dynamical evolution. We also apply the dynamical model, enriched with results coming from stellar evolutionary models, to demonstrate that in S-type systems around low-mass stars, with relative high inclinations ($\itot\ge 60^\circ$), moderate eccentricities ($0.2\le e\le 0.4$) and planets located around 1 AU, the evolution of stellar radius during the first few hundreds of Myr, alters significantly the timescales of dynamical evolution.


title: Platoon Forming Algorithms for Intelligent Street Intersections

abstract: We study intersection access control for autonomous vehicles. Platoon forming algorithms, which aim to organize individual vehicles in platoons, are very promising. To create those platoons, we slow down vehicles before the actual arrival at the intersection in such a way that each vehicle can traverse the intersection at high speed. This increases the capacity of the intersection significantly, offering huge potential savings with respect to travel time compared to nowadays traffic. We propose several new platoon forming algorithms and provide an approximate mean delay analysis for our algorithms. A comparison between the current day practice at intersections (through a case study in SUMO) and our proposed algorithms is provided. Simulation results for fairness are obtained as well, showing that platoon forming algorithms with a low mean delay sometimes are relatively unfair, indicating a potential need for balancing mean delay and fairness.


title: Towards Personalized Management of Type B Aortic Dissection Using STENT: a STandard cta database with annotation of the ENtire aorta and True-false lumen

abstract: Type B Aortic Dissection(TBAD) is a rare aortic disease with a high 5-year mortality.Personalized and precise management of TBAD has been increasingly desired in clinic which requires the geometric parameters of TBAD specific to the patient be measured accurately.This remains to be a challenging task for vascular surgeons as manual measurement is highly subjective and imprecise. To solve this problem,we introduce STENT-a STandard cta database with annotation of the ENtire aorta and True-false lumen. The database contains 274 CT angiography (CTA) scans from 274 unique TBAD patients and is split into a training set(254 cases including 210 preoperative and 44 postoperative scans ) and a test set(20 cases).Based on STENT,we develop a series of methods including automated TBAD segmentation and automated measurement of TBAD parameters that facilitate personalized and precise management of the disease. In this work, the database and the proposed methods are thoroughly introduced and evaluated and the results of our study shows the feasibility and effectiveness of our approach to easing the decision-making process for vascular surgeons during personalized TBAD management.


title: Agent-Based Modelling Approach for Distributed Decision Support in an IoT Network

abstract: An increasing number of emerging applications, e.g., internet of things, vehicular communications, augmented reality, and the growing complexity due to the interoperability requirements of these systems, lead to the need to change the tools used for the modeling and analysis of those networks. Agent-Based Modeling (ABM) as a bottom-up modeling approach considers a network of autonomous agents interacting with each other, and therefore represents an ideal framework to comprehend the interactions of heterogeneous nodes in a complex environment. Here, we investigate the suitability of ABM to model the communication aspects of a road traffic management system, as an example of an Internet of Things (IoT) network. We model, analyze and compare various Medium Access Control (MAC) layer protocols for two different scenarios, namely uncoordinated and coordinated. Besides, we model the scheduling mechanisms for the coordinated scenario as a high level MAC protocol by using three different approaches: Centralized Decision Maker, DESYNC and decentralized learning MAC (L-MAC). The results clearly show the importance of coordination between multiple decision makers in order to improve the accuracy of information and spectrum utilization of the system.


title: Spectroscopy of Extremal (and Near-Extremal) Kerr Black Holes

abstract: We investigate linear, spin-field perturbations of Kerr black holes in the extremal limit throughout the complex-frequency domain. We calculate quasi-normal modes of extremal Kerr, as well as of near-extremal Kerr, via a novel approach: using the method of Mano, Suzuki and Takasugi (MST). We also show how, in the extremal limit, a branch cut is formed at the superradiant-bound frequency, $\omega_{SR}$, via a simultaneous accumulation of quasi-normal modes and totally-reflected modes. For real frequencies, we calculate the superradiant amplification factor, which yields the amount of rotational energy that can be extracted from a black hole. In the extremal limit, this factor is the largest and it displays a discontinuity at $\omega_{SR}$ for some modes. Finally, we find no exponentially-growing modes nor branch points on the upper-frequency plane in extremal Kerr after a numerical investigation, thus providing evidence of the mode-stability of this space-time away from the horizon.


title: Human few-shot learning of compositional instructions

abstract: People learn in fast and flexible ways that have not been emulated by machines. Once a person learns a new verb "dax," he or she can effortlessly understand how to "dax twice," "walk and dax," or "dax vigorously." There have been striking recent improvements in machine learning for natural language processing, yet the best algorithms require vast amounts of experience and struggle to generalize new concepts in compositional ways. To better understand these distinctively human abilities, we study the compositional skills of people through language-like instruction learning tasks. Our results show that people can learn and use novel functional concepts from very few examples (few-shot learning), successfully applying familiar functions to novel inputs. People can also compose concepts in complex ways that go beyond the provided demonstrations. Two additional experiments examined the assumptions and inductive biases that people make when solving these tasks, revealing three biases: mutual exclusivity, one-to-one mappings, and iconic concatenation. We discuss the implications for cognitive modeling and the potential for building machines with more human-like language learning capabilities.


title: Optimal Needle Diameter, Shape, and Path in Autonomous Suturing

abstract: Needle shape, diameter, and path are critical parameters that directly affect suture depth and tissue trauma in autonomous suturing. This paper presents an optimization-based approach to specify these parameters. Given clinical suturing guidelines, a kinematic model of needle-tissue interaction was developed to quantify suture parameters and constraints. The model was further used to formulate constant curvature needle path planning as a nonlinear optimization problem. The optimization results were confirmed experimentally with the Raven II surgical system. The proposed needle path planning algorithm guarantees minimal tissue trauma and complies with a wide range of suturing requirements.


title: The integral Cauchy problem for generalized Boussinesq equations with general leading parts

abstract: In this paper, the integral initial value problems for Boussinesq type equations are studied. The equation include the general differential operators. The existence, uniqueness and regularity properties of solution of these problems are obtained. By choosing differential operators including in the equation, the regularity properties of the Cauchy problem for different type of Boussinesg equations are studied.


title: A filtration on the higher Chow group of zero cycles on an abelian variety

abstract: In this paper we extend Gazaki's results on the Chow groups of abelian varieties to the higher Chow groups. We introduce a Gazaki type filtration on the higher Chow group of zero-cycles on an abelian variety, whose graded quotients are connected to the Somekawa type $K$-group. Via the \'{e}tale cycle map, we will compare this filtration with a filtration on the \'{e}tale cohomology induced by the Hochschild-Serre spectral sequence. As an application over local fields, we obtain an estimate of the kernel of the reciprocity map.


title: On the Alperin-McKay conjecture for simple groups of type $\mathrm{A}$

abstract: In this paper characters of the normaliser of $d$-split Levi subgroups in $\mathrm {SL}_n(q)$ and $\mathrm {SU}_n(q)$ are parametrized with a particular focus on the Clifford theory between the Levi subgroup and its normalizer.These results are applied to verify the Alperin-McKay conjecture for primes $\ell$ with $\ell\nmid 6(q^2-1)$ and the Alperin weight conjecture for $\ell$-blocks of those quasi-simple groups with abelian defect. The inductive Alperin-McKay condition and inductive Alperin weight condition by the second author are verified for certain blocks of $\mathrm {SL}_n(q)$ and $\mathrm {SU}_n(q)$.


title: Interpretable machine learning: definitions, methods, and applications

abstract: Machine-learning models have demonstrated great success in learning complex patterns that enable them to make predictions about unobserved data. In addition to using models for prediction, the ability to interpret what a model has learned is receiving an increasing amount of attention. However, this increased focus has led to considerable confusion about the notion of interpretability. In particular, it is unclear how the wide array of proposed interpretation methods are related, and what common concepts can be used to evaluate them. We aim to address these concerns by defining interpretability in the context of machine learning and introducing the Predictive, Descriptive, Relevant (PDR) framework for discussing interpretations. The PDR framework provides three overarching desiderata for evaluation: predictive accuracy, descriptive accuracy and relevancy, with relevancy judged relative to a human audience. Moreover, to help manage the deluge of interpretation methods, we introduce a categorization of existing techniques into model-based and post-hoc categories, with sub-groups including sparsity, modularity and simulatability. To demonstrate how practitioners can use the PDR framework to evaluate and understand interpretations, we provide numerous real-world examples. These examples highlight the often under-appreciated role played by human audiences in discussions of interpretability. Finally, based on our framework, we discuss limitations of existing methods and directions for future work. We hope that this work will provide a common vocabulary that will make it easier for both practitioners and researchers to discuss and choose from the full range of interpretation methods.


title: Magnetized interstellar molecular clouds: II. The Large-Scale Structure and Dynamics of Filamentary Molecular Clouds

abstract: We perform ideal MHD high resolution AMR simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ~4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fiber/web-like sub-structures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to 1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shape configuration which are strikingly similar to the Infrared dark cloud SDC13 leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (M_A ~ 1) is crucial for maintaining a long and slender filamentary cloud for a long period of time ~0.5 million years.


title: Phonon and electronic structures and resistance of layered electride Ca2N: DFT calculations

abstract: The phonon and electronic properties, the Eliashberg function and the temperature dependence of resistance of electride Ca2N are investigated by the DFT-LDA plane-wave method. The phonon dispersion, the partial phonon density of states and the atomic eigenvectors of zero-center phonons are studied. The electronic band dispersion and partial density of states conclude that Ca2N is a metal and the Ca 3p, 4s and N 2p orbitals are hybridized. For the analysis of an electron - phonon interaction (EPI) and its contribution to resistance the Eliashberg function was calculated and a temperature dependence of resistance caused EPI was found. The present results are in good agreement with experiment data.


title: Artificial Intelligence Assists Discovery of Reaction Coordinates and Mechanisms from Molecular Dynamics Simulations

abstract: Exascale computing holds great opportunities for molecular dynamics (MD) simulations. However, to take full advantage of the new possibilities, we must learn how to focus computational power on the discovery of complex molecular mechanisms, and how to extract them from enormous amounts of data. Both aspects still rely heavily on human experts, which becomes a serious bottleneck when a large number of parallel simulations have to be orchestrated to take full advantage of the available computing power. Here, we use artificial intelligence (AI) both to guide the sampling and to extract the relevant mechanistic information. We combine advanced sampling schemes with statistical inference, artificial neural networks, and deep learning to discover molecular mechanisms from MD simulations. Our framework adaptively and autonomously initializes simulations and learns the sampled mechanism, and is thus suitable for massively parallel computing architectures. We propose practical solutions to make the neural networks interpretable, as illustrated in applications to molecular systems.


title: AET vs. AED: Unsupervised Representation Learning by Auto-Encoding Transformations rather than Data

abstract: The success of deep neural networks often relies on a large amount of labeled examples, which can be difficult to obtain in many real scenarios. To address this challenge, unsupervised methods are strongly preferred for training neural networks without using any labeled data. In this paper, we present a novel paradigm of unsupervised representation learning by Auto-Encoding Transformation (AET) in contrast to the conventional Auto-Encoding Data (AED) approach. Given a randomly sampled transformation, AET seeks to predict it merely from the encoded features as accurately as possible at the output end. The idea is the following: as long as the unsupervised features successfully encode the essential information about the visual structures of original and transformed images, the transformation can be well predicted. We will show that this AET paradigm allows us to instantiate a large variety of transformations, from parameterized, to non-parameterized and GAN-induced ones. Our experiments show that AET greatly improves over existing unsupervised approaches, setting new state-of-the-art performances being greatly closer to the upper bounds by their fully supervised counterparts on CIFAR-10, ImageNet and Places datasets.


title: Exact Power of the Rank-Sum Test for a Continuous Variable

abstract: Accurate power calculations are essential in small studies containing expensive experimental units or high-stakes exposures. Herein, exact power of the Wilcoxon Mann-Whitney rank-sum test of a continuous variable is formulated using a Monte Carlo approach and defining P(X < Y) = p as a measure of effect size, where X and Y denote random observations from two distributions hypothesized to be equal under the null. Effect size p fosters productive communications because researchers understand p = 0.5 is analogous to a fair coin toss, and p near 0 or 1 represents a large effect. This approach is feasible even without background data. Simulations were conducted comparing the exact power approach to existing approaches by Rosner & Glynn (2009), Shieh et al. (2006), Noether (1987), and O'Brien-Castelloe (2006). Approximations by Noether and O'Brien-Castelloe are shown to be inaccurate for small sample sizes. The Rosner & Glynn and Shieh et al. approaches performed well in many small sample scenarios, though both are restricted to location-shift alternatives and neither approach is theoretically justified for small samples. The exact method is recommended and available in the R package wmwpow. KEYWORDS: Mann-Whitney test, Monte Carlo simulation, non-parametric, power analysis, Wilcoxon rank-sum test


title: Precision Annealing Monte Carlo Methods for Statistical Data Assimilation: Metropolis-Hastings Procedures

abstract: Statistical Data Assimilation (SDA) is the transfer of information from field or laboratory observations to a user selected model of the dynamical system producing those observations. The data is noisy and the model has errors; the information transfer addresses properties of the conditional probability distribution of the states of the model conditioned on the observations. The quantities of interest in SDA are the conditional expected values of functions of the model state, and these require the approximate evaluation of high dimensional integrals. We introduce a conditional probability distribution and use the Laplace method with annealing to identify the maxima of the conditional probability distribution. The annealing method slowly increases the precision term of the model as it enters the Laplace method. In this paper, we extend the idea of precision annealing (PA) to Monte Carlo calculations of conditional expected values using Metropolis-Hastings methods.


title: An Ensemble Method for Interval-Censored Time-to-Event Data

abstract: Interval-censored data analysis is important in biomedical statistics for any type of time-to-event response where the time of response is not known exactly, but rather only known to occur between two assessment times. Many clinical trials and longitudinal studies generate interval-censored data; one common example occurs in medical studies that entail periodic follow-up. In this paper we propose a survival forest method for interval-censored data based on the conditional inference framework. We describe how this framework can be adapted to the situation of interval-censored data. We show that the tuning parameters have a non-negligible effect on the survival forest performance and guidance is provided on how to tune the parameters in a data-dependent way to improve the overall performance of the method. Using Monte Carlo simulations we find that the proposed survival forest is at least as effective as a survival tree method when the underlying model has a tree structure, performs similarly to an interval-censored Cox proportional hazards model fit when the true relationship is linear, and outperforms the survival tree method and Cox model when the true relationship is nonlinear. We illustrate the application of the method on a tooth emergence data set.


title: Probing the transport properties of Quark-Gluon Plasma via heavy-flavor Boltzmann and Langevin dynamics

abstract: The heavy quark propagation behavior inside the quark-gluon plasma (QGP), is usually described in terms of the Boltzmann dynamics, which can be reduced to the Langevin approach by assuming a small momentum transfer for the scattering processes between heavy quarks and the QGP constituents. In this work, the temperature and energy dependence of the transport coefficients are calculated in the framework of both Boltzmann and Langevin dynamics. The derived transport coefficients are found to be systematically larger in the Boltzmann approach as compared with the Langevin, in particular in the high temperature and high energy region. Within each of the two theoretical frameworks, we simulate the charm quark production and the subsequent evolution processes in relativistic heavy-ion collisions. We find that the total in-medium energy loss is larger from the Langevin dynamics, resulting in a smaller (larger) $R_{\rm AA}$ at high (low) $p_{\rm T}$, for both the charm quark and heavy-flavor mesons. Meanwhile, the Boltzmann model is found to induce larger $v_{\rm 2}$, in particular at moderate $p_{\rm T}$, as well as stronger broadening behavior for the azimuthal distributions. By comparing the model calculations with available experimental measurements for D-mesons, we find that the Langevin approach is more favored by the $R_{\rm AA}$ data while the Boltzmann approach is more favored favor by the $v_{\rm 2}$ data. A simultaneous description of both observables appear challenging for both models.


title: Equilibrium time-correlation functions of the long-range interacting Fermi-Pasta-Ulam model

abstract: We present a numerical study of dynamical correlations (structure factors) of the long-range generalization of the Fermi-Pasta-Ulam oscillator chain, where the strength of the interaction between two lattice sites decays as a power $\alpha$ of the inverse of their distance. The structure factors at finite energy density display distinct peaks, corresponding to long-wavelength propagating modes, whose dispersion relation is compatible with the predictions of the linear theory. We demonstrate that dynamical scaling holds, with a dynamical exponent $z$ that depends weakly on $\alpha$ in the range $1<\alpha<3$. The lineshapes have a non-trivial functional form and appear somehow independent of $\alpha$. Within the accessible time and size ranges, we also find that the short-range limit is hardly attained even for relatively large values of $\alpha$.