title: Fates of the oldest intermediate-mass stars - Primordial to Extremely Metal-Poor AGB and Super-AGB Stars: White Dwarf or Supernova progenitors

abstract: Abridged: Getting a better understanding of the evolution and nucleosynthetic yields of the most metal-poor stars (appr. Z<=10^-5) is critical because they are part of the big picture of the history of the primitive Universe. Yet many of the remaining unknowns of stellar evolution lie in the birth, life, and death of these objects. We review stellar evolution of intermediate-mass (IMS) Z<=10-5 models existing in the literature, with a focus on the problem of their final fates. The depth and efficiency of mixing episodes are critical to determine the mass limits for the formation of electron-capture supernovae, but our knowledge of these phenomena is not complete because they are strongly affected by the choice of input physics. We also consider the alternative SNI1/2 channel to form SNe out of the most metal-poor IMS. In this case, it is critical to understand the thermally-pulsing AGB evolution until the late stages. Efficient second dredge-up and, later, third dredge-up episodes could be able to pollute stellar envelopes enough for the stars to undergo thermal pulses in a way very similar to that of higher initial Z objects. Inefficient 2nd and/or 3rd dredge-up may leave an almost pristine envelope, unable to sustain strong stellar winds. This may allow the H-exhausted core to grow to M_Ch before the envelope is lost, and thus let the star explode as a SNI1/2. After reviewing the information available on these two possible channels for the formation of SNe, we discuss existing nucleosynthetic yields of stars of metallicity Z<=10^-5, and present an example of nucleosynthetic calculations for a thermally-pulsing Super-AGB star of Z=10^-5. We compare theoretical predictions with observations of the lowest [Fe/H] objects detected. The review closes by discussing current open questions as well as possible fruitful avenues for future research.


title: Quantum critical point of the Ising chain from boundary effects

abstract: We propose two easy-to-study observables in the quantum Ising chain with open boundary conditions. They measure the length at which boundaries affect the longitudinal or transverse magnetization. We show that their finite-size scaling behaviour encodes the position of the quantum critical point and the universal scaling exponent $\nu$. The applicability of proposed observables in small systems is also discussed. We expect that our results will be useful in quantum simulation of spin systems.


title: Dimensions of triangle sets

abstract: In this paper, we discuss some dimension results for triangle sets of compact sets in $\mathbb{R}^2$. In particular, we prove that for any compact set $F$ in $\mathbb{R}^2$, the triangle set $\Delta(F)$ satisfies \[ \dim_{\mathrm{A}} \Delta(F)\geq \frac{3}{2}\dim_{\mathrm{A}} F. \] If $\dim_{\mathrm{A}} F>1$ then we have \[ \dim_{\mathrm{A}} \Delta(F)\geq 1+\dim_{\mathrm{A}} F. \] If $\dim_{\mathrm{A}} F>4/3$ then we have the following better bound, \[ \dim_{\mathrm{A}} \Delta(F)\geq \min\left\{\frac{5}{2}\dim_{\mathrm{A}} F-1,3\right\}. \] Moreover, if $F$ satisfies a mild separation condition then the above result holds also for the box dimensions, namely, \[ \underline{\dim_{\mathrm{B}}} F\geq \frac{3}{2}\underline{\dim_{\mathrm{B}}} \Delta(F) \text{ and }\overline{\dim_{\mathrm{B}}} F\geq \frac{3}{2}\overline{\dim_{\mathrm{B}}} \Delta(F). \]


title: Selectivity correction in discrete-continuous models for the willingness to work as crowd-shippers and travel time tolerance

abstract: The objective of this study is to understand the different behavioral considerations that govern the choice of people to engage in a crowd-shipping market. Using novel data collected by the researchers in the US, we develop discrete-continuous models. A binary logit model has been used to estimate crowd-shippers' willingness to work, and an ordinary least-square regression model has been employed to calculate crowd-shippers' maximum tolerance for shipping and delivery times. A selectivity-bias term has been included in the model to correct for the conditional relationships of the crowd-shipper's willingness to work and their maximum travel time tolerance. The results show socio-demographic characteristics (e.g. age, gender, race, income, and education level), transporting freight experience, and number of social media usages significant influence the decision to participate in the crowd-shipping market. In addition, crowd-shippers pay expectations were found to be reasonable and concurrent with the literature on value-of-time. Findings from this research are helpful for crowd-shipping companies to identify and attract potential shippers. In addition, an understanding of crowd-shippers - their behaviors, perceptions, demographics, pay expectations, and in which contexts they are willing to divert from their route - are valuable to the development of business strategies such as matching criteria and compensation schemes for driver-partners.


title: Gyroscope-Aided Motion Deblurring with Deep Networks

abstract: We propose a deblurring method that incorporates gyroscope measurements into a convolutional neural network (CNN). With the help of such measurements, it can handle extremely strong and spatially-variant motion blur. At the same time, the image data is used to overcome the limitations of gyro-based blur estimation. To train our network, we also introduce a novel way of generating realistic training data using the gyroscope. The evaluation shows a clear improvement in visual quality over the state-of-the-art while achieving real-time performance. Furthermore, the method is shown to improve the performance of existing feature detectors and descriptors against the motion blur.


title: On GILP's group theoretic approach to Falconer's distance problem

abstract: In this paper, we follow and extend a group-theoretic method introduced by Greenleaf-Iosevich-Liu-Palsson (GILP) to study finite points configurations spanned by Borel sets in $\mathbb{R}^n,n\geq 2,n\in\mathbb{N}.$ We remove a technical continuity condition in a GILP's theorem in [GILP15]. This allows us to extend the Wolff-Erdogan dimension bound for distance sets to finite points configurations with $k$ points for $k\in\{2,\dots,n+1\}.$ At the end of this paper, we extend this group-theoretic method and illustrate a `Fourier free' approach to Falconer's distance set problem for the Lebesgue measure. We explain how to use tubular incidence estimates in distance set problems. Curiously, tubular incidence estimates are also related to the Kakeya problem.


title: Heterogeneous MacroTasking (HeMT) for Parallel Processing in the Public Cloud

abstract: Using tiny, equal-sized tasks (Homogeneous microTasking, HomT) has long been regarded an effective way of load balancing in parallel computing systems. When combined with nodes pulling in work upon becoming idle, HomT has the desirable property of automatically adapting its load distribution to the processing capacities of participating nodes - more powerful nodes finish their work sooner and, therefore, pull in additional work faster. As a result, HomT is deemed especially desirable in settings with heterogeneous (and possibly possessing dynamically changing) processing capacities. However, HomT does have additional scheduling and I/O overheads that might make this load balancing scheme costly in some scenarios. In this paper, we first analyze these advantages and disadvantages of HomT. We then propose an alternative load balancing scheme - Heterogeneous MacroTasking (HeMT) - wherein workload is intentionally partitioned according to nodes' processing capacity. Our goal is to study when HeMT is able to overcome the performance disadvantages of HomT. We implement a prototype of HeMT within the Apache Spark application framework with complementary enhancements to the Apache Mesos cluster manager. Spark's built-in scheduler, when parameterized appropriately, implements HomT. Our experimental results show that HeMT out-performs HomT when accurate workload-specific estimates of nodes' processing capacities are learned. As representative results, Spark with HeMT offers about 10% better average completion times for realistic data processing workloads over the default system.


title: Probing BSM and High-x Physics with SoLID at JLab

abstract: The program of parity violation with the proposed new SoLID spectrometer at JLab is presented. Physics topics include searched for physics beyond the Standard Model, studies of charge symmetry violation at the quark level, searched for quark-quark correlations, and a measurement of the ratio of up and down PDF's in the proton.


title: Finite index theorems for iterated Galois groups of unicritical polynomials

abstract: Let $K$ be the function field of a smooth, irreducible curve defined over $\overline{\mathbb{Q}}$. Let $f\in K[x]$ be of the form $f(x)=x^q+c$ where $q = p^{r}, r \ge 1,$ is a power of the prime number $p$, and let $\beta\in \overline{K}$. For all $n\in\mathbb{N}\cup\{\infty\}$, the Galois groups $G_n(\beta)=\mathop{\rm{Gal}}(K(f^{-n}(\beta))/K(\beta))$ embed into $[C_q]^n$, the $n$-fold wreath product of the cyclic group $C_q$. We show that if $f$ is not isotrivial, then $[[C_q]^\infty:G_\infty(\beta)]<\infty$ unless $\beta$ is postcritical or periodic. We are also able to prove that if $f_1(x)=x^q+c_1$ and $f_2(x)=x^q+c_2$ are two such distinct polynomials, then the fields $\bigcup_{n=1}^\infty K(f_1^{-n}(\beta))$ and $\bigcup_{n=1}^\infty K(f_2^{-n}(\beta))$ are disjoint over a finite extension of $K$.


title: Real-space probe for lattice quasiholes

abstract: We propose a real-space probe that is based on density measurements to extract distinct signatures of quasihole-like states of bosons experiencing a synthetic magnetic field in a two-dimensional lattice. We numerically show that certain ratios of the mean square radii of the particle cloud, obtainable through the density profile, approach the continuum values expected from Laughlin's ansatz wave functions quickly as the magnetic flux quanta per unit cell of the lattice decrease, even in a small lattice with few particles. This method could equally be used in both ultra-cold atomic and photonic systems.


title: Exact stationary solutions of the Kolmogorov-Feller equation in a bounded domain

abstract: We present the first detailed analysis of the statistical properties of jump processes bounded by a saturation function and driven by Poisson white noise, being a random sequence of delta pulses. The Kolmogorov-Feller equation for the probability density function (PDF) of such processes is derived and its stationary solutions are found analytically in the case of the symmetric uniform distribution of pulse sizes. Surprisingly, these solutions can exhibit very complex behavior arising from both the boundedness of pulses and processes. We show that all features of the stationary PDF (number of branches, their form, extreme values probability, etc.) are completely determined by the ratio of the saturation function width to the half-width of the pulse-size distribution. We verify all theoretical results by direct numerical simulations.


title: Ballot Permutations and Odd Order Permutations

abstract: A permutation $\pi$ is ballot if, for all $k$, the word $\pi_1\cdots \pi_k$ has at least as many ascents as it has descents. Let $b(n)$ denote the number of ballot permutations of order $n$, and let $p(n)$ denote the number of permutations which have odd order in the symmetric group $S_n$. Callan conjectured that $b(n)=p(n)$ for all $n$, which was proved by Bernardi, Duplantier, and Nadeau. We propose a refinement of Callan's original conjecture. Let $b(n,d)$ denote the number of ballot permutations with $d$ descents. Let $p(n,d)$ denote the number of odd order permutations with $M(\pi)=d$, where $M(\pi)$ is a certain statistic related to the cyclic descents of $\pi$. We conjecture that $b(n,d)=p(n,d)$ for all $n$ and $d$. We prove this stronger conjecture for the cases $d=1,\ 2,\ 3$, and $d=\lfloor(n-1)/2\rfloor$, and in each of these cases we establish formulas for $b(n,d)$ involving Eulerian numbers and Eulerian-Catalan numbers.


title: Bounds on Instantaneous Nonlocal Quantum Computation

abstract: Instantaneous nonlocal quantum computation refers to a process in which spacelike separated parties simulate a nonlocal quantum operation on their joint systems through the consumption of pre-shared entanglement. To prevent a violation of causality, this simulation succeeds up to local errors that can only be corrected after the parties communicate classically with one another. However, this communication is non-interactive, and it involves just the broadcasting of local measurement outcomes. We refer to this operational paradigm as local operations and broadcast communication (LOBC) to distinguish it from the standard local operations and (interactive) classical communication (LOCC). In this paper, we show that an arbitrary two-qubit gate can be implemented by LOBC with $\epsilon$-error using $O(\log(1/\epsilon))$ entangled bits (ebits). This offers an exponential improvement over the best known two-qubit protocols, whose ebit costs behave as $O(1/\epsilon)$. We also consider the family of binary controlled gates on dimensions $d_A\otimes d_B$. We find that any hermitian gate of this form can be implemented by LOBC using a single shared ebit. In sharp contrast, a lower bound of $\log d_B$ ebits is shown in the case of generic (i.e. non-hermitian) gates from this family, even when $d_A=2$. This demonstrates an unbounded gap between the entanglement costs of LOCC and LOBC gate implementation. Whereas previous lower bounds on the entanglement cost for instantaneous nonlocal computation restrict the minimum dimension of the needed entanglement, we bound its entanglement entropy. To our knowledge this is the first such lower bound of its kind.


title: Searches for gamma-ray lines and `pure WIMP' spectra from Dark Matter annihilations in dwarf galaxies with H.E.S.S

abstract: Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95\% confidence level of about $\langle \sigma v \rangle \lesssim 3 \times 10^{-25}$ cm$^3$ s$^{-1}$ are obtained in the mass range of 400 GeV to 1 TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.


title: Recent Results on Heavy Quark Production in High Energy Nucleus Collisions

abstract: The goal of the ultra-relativistic heavy ion program is to study Quantum Chromodynamics under finite temperature and density conditions. After a couple of decades of experiment, the focus at the top RHIC and the LHC energy has evolved to quantitative understanding of properties of the hot and dense medium, namely the strongly-coupled Quark Gluon Plasma (sQGP) created in these heavy-ion collisions, and to constrain transport parameters of the sQGP medium. Heavy quarks offer unique insights towards detailed understanding of the sQGP properties due to their large masses. Recent collider and detector advances have enabled precision measurements of heavy quark hadron production in heavy-ion collisions. In this talk, I will review recent results from heavy quark production measurements at RHIC and the LHC. These high quality data will offer stringent constraints on theoretical model calculations and help precision determination of QGP medium transport parameters. Finally, I look forward to a more prospective future of the heavy quark program with further improved detectors at both RHIC and the LHC.


title: A simple parameter-free and adaptive approach to optimization under a minimal local smoothness assumption

abstract: We study the problem of optimizing a function under a \emph{budgeted number of evaluations}. We only assume that the function is \emph{locally} smooth around one of its global optima. The difficulty of optimization is measured in terms of 1) the amount of \emph{noise} $b$ of the function evaluation and 2) the local smoothness, $d$, of the function. A smaller $d$ results in smaller optimization error. We come with a new, simple, and parameter-free approach. First, for all values of $b$ and $d$, this approach recovers at least the state-of-the-art regret guarantees. Second, our approach additionally obtains these results while being \textit{agnostic} to the values of both $b$ and $d$. This leads to the first algorithm that naturally adapts to an \textit{unknown} range of noise $b$ and leads to significant improvements in a moderate and low-noise regime. Third, our approach also obtains a remarkable improvement over the state-of-the-art SOO algorithm when the noise is very low which includes the case of optimization under deterministic feedback ($b=0$). There, under our minimal local smoothness assumption, this improvement is of exponential magnitude and holds for a class of functions that covers the vast majority of functions that practitioners optimize ($d=0$). We show that our algorithmic improvement is borne out in experiments as we empirically show faster convergence on common benchmarks.


title: Automated sequence and motion planning for robotic spatial extrusion of 3D trusses

abstract: While robotic spatial extrusion has demonstrated a new and efficient means to fabricate 3D truss structures in architectural scale, a major challenge remains in automatically planning extrusion sequence and robotic motion for trusses with unconstrained topologies. This paper presents the first attempt in the field to rigorously formulate the extrusion sequence and motion planning (SAMP) problem, using a CSP encoding. Furthermore, this research proposes a new hierarchical planning framework to solve the extrusion SAMP problems that usually have a long planning horizon and 3D configuration complexity. By decoupling sequence and motion planning, the planning framework is able to efficiently solve the extrusion sequence, end-effector poses, joint configurations, and transition trajectories for spatial trusses with nonstandard topologies. This paper also presents the first detailed computation data to reveal the runtime bottleneck on solving SAMP problems, which provides insight and comparing baseline for future algorithmic development. Together with the algorithmic results, this paper also presents an open-source and modularized software implementation called Choreo that is machine-agnostic. To demonstrate the power of this algorithmic framework, three case studies, including real fabrication and simulation results, are presented.


title: Searches for rare and non-Standard Model decays of the Higgs boson

abstract: Discovered in 2012, the Higgs boson has opened a new window on nature. The latest searches for its rare and non-Standard Model decays with the ATLAS detector at the LHC are presented. They represent a promising probe of the 1st and 2nd generation Yukawa couplings, and of new physics. Searches for rare exclusive decays of the Higgs boson to a meson and a photon are presented. Also presented are four searches for decays of the Higgs boson to pairs of beyond the Standard Model resonances, in various final states.


title: MiniBooNE, MINOS+ and IceCube data imply a baroque neutrino sector

abstract: The 4.8$\sigma$ anomaly in MiniBooNE data cannot be reconciled with MINOS+ and IceCube data within the vanilla framework of neutrino oscillations involving an eV-mass sterile neutrino. We show that an apparently consistent picture can be drawn if charged-current and neutral-current nonstandard neutrino interactions are at work in the 3+1 neutrino scheme. It appears that either the neutrino sector is more elaborate than usually envisioned, or one or more datasets needs revision.


title: A beamline for fundamental neutron physics at TRIUMF

abstract: This article describes the new primary proton beamline 1U at TRIUMF. The purpose of this beamline is to produce ultracold neutrons (UCN) for fundamental-physics experiments. It delivers up to 40 microA of 480 MeV protons from the TRIUMF cyclotron to a tungsten spallation target and uses a fast kicker to share the beam between the Center for Molecular and Materials Science and UCN. The beamline has been successfully commissioned and operated with a beam current up to 10 microA, facilitating first large-scale UCN production in Canada.


title: Probing the Strange Sea Quarks with Kaon SIDIS

abstract: It is well known that protons and neutrons are made from constituents, called quarks and gluons, which give substructure to these particles. The goal of this project is to make measurements of the spatial distributions and the momenta of the quarks that provide a three-dimensional map of quarks in the nuclear medium. This knowledge provides the basis of our understanding of nuclear matter in terms of the dynamics of their internal constituents. This abstract focuses on the study of the contribution of the sea quarks and in particular of the strange sea to the proton spin structure. This study is feasible with semi-inclusive deep inelastic scattering of electrons off proton and deuteron targets in Hall B at Jefferson Lab. To achieve the desired precision, a Ring Imaging CHerenkov (RICH) detector was built so pion, kaon and proton identification is well performed in the momentum range of 3 to 8 GeV/c. The experimental method and projected precision of the measurements of the parton distributions using Kaon SIDIS will be discussed and the status of the recently built Hybrid RICH detector will be presented.


title: Critical groups of van Lint-Schrijver Cyclotomic Strongly Regular Graphs

abstract: The \emph{critical} group of a finite connected graph is an abelian group defined by the Smith normal form of its Laplacian. Let $q$ be a power of a prime and $H$ be a multiplicative subgroup of $K=\mathbb{F}_{q}$. By $\mathrm{Cay}(K,H)$ we denote the Cayley graph on the additive group of $K$ with `connection' set $H$. A strongly regular graph of the form $\mathrm{Cay}(K,H)$ is called a \emph{cyclotomic strongly regular graph}. Let $p$ and $\ell >2$ be primes such that $p$ is primitive $\pmod{\ell}$. We compute the \emph{critical} groups of a family of \emph{cyclotomic strongly regular graphs} for which $q=p^{(\ell-1)t}$ (with $t\in \mathbb{N}$) and $H$ is the unique multiplicative subgroup of order $k=\frac{q-1}{\ell}$. These graphs were first discovered by van Lint and Schrijver in \cite{VS}.


title: The role of incidence angle in the laser ablation of planar ICF targets

abstract: The effect of the laser ray incidence angle on the mass ablation rate and ablation pressure of planar inertial confinement fusion (ICF) targets is explored using an idealized model. Polar direct drive (PDD) on the National Ignition Facility (NIF) requires the repointing of its 192 beams clustered within 50 degrees of the poles to minimize the imparted polar varying payload kinetic energy of the target. Due to this repointing, non-normal incidence angles of the beam centerlines are encountered in any PDD design. The formulation of a PDD scheme that minimizes non-uniformity is a significant challenge that requires an understanding of the induced differences in ablation including those of incidence angle. In this work, a modified version of the textbook model of laser ablation [Manheimer et al. Phys. Fluids 25, 1644 (1982)] is used to demonstrate that the mass ablation rate and ablation pressure scale with the 4/3 and 2/3 power of the cosine of the laser ray incidence angle for the planar case during the beginning portion of the laser ablation. This result is considered with an idealized 1-D two segment planar model using rays at different incidence angles for each segment. Within the model, the segments cannot have equal mass and velocity simultaneously without tailoring the ray's time-dependent intensity for each segment. However, with the correct segment-to-segment time-independent ray intensities, equal velocity and dynamic pressure can be achieved approximately without tailoring. Additionally, an analytic prediction for the conduction zone width as a function of incidence angle is provided. It is predicted that this width increases with incidence angle, resulting in a decrease in laser imprint, an effect which has been previously observed in experiments. These results, when generalized to spherical geometry, may provide insight into the implosion dynamics encountered in PDD.


title: plsRglm: Partial least squares linear and generalized linear regression for processing incomplete datasets by cross-validation and bootstrap techniques with R

abstract: The aim of the plsRglm package is to deal with complete and incomplete datasets through several new techniques or, at least, some which were not yet implemented in R. Indeed, not only does it make available the extension of the PLS regression to the generalized linear regression models, but also bootstrap techniques, leave-one-out and repeated $k$-fold cross-validation. In addition, graphical displays help the user to assess the significance of the predictors when using bootstrap techniques. Biplots (Fig. 4) can be used to delve into the relationship between individuals and variables.


title: Energy-momentum tensor of a ferromagnet

abstract: The energy-momentum tensor of a ferromagnet derived according to the standard prescription of Noether's theorem has a major flaw: the term originating from the spin Berry phase is gauge-dependent. As a consequence, some physical quantities computed from the tensor show unphysical behavior. For example, the presence of a spin-polarized current does not affect the energy of the domain wall in the commonly accepted gauge, which implies-incorrectly-the absence of the adiabatic spin torque. In other gauges, the spin torque shows unphysical glitches occurring when the plane of magnetization crosses the Dirac string associated with a magnetic monopole in spin space. We derive a gauge-invariant energy-momentum tensor that is free from these artifacts but requires the addition of an extra spatial dimension, with the ferromagnet living on its boundary. It can be obtained most directly from the Wess-Zumino action for spins, which relies on the same extra dimension.


title: How UX Practitioners Produce Findings in Usability Testing

abstract: Usability testing has long been a core interest of HCI research and forms a key element of industry practice. Yet our knowledge of it harbours striking absences. There are few, if any detailed accounts of the contingent, material ways in which usability testing is actually practiced. Further, it is rare that industry practitioners' testing work is treated as indigenous and particular (instead subordinated as a `compromised' version). To service these problems, this paper presents an ethnomethodological study of usability testing practices in a design consultancy. It unpacks how findings are produced in and as the work of observers analysing the test as it unfolds between moderators taking participants through relevant tasks. The study nuances conventional views of usability findings as straightforwardly `there to be found' or `read off' by competent evaluators. It explores how evaluators / observers collaboratively work to locate relevant troubles in the test's unfolding. However, in the course of doing this work, potential candidate troubles may also routinely be dissipated and effectively `ignored' in one way or another. The implications of the study suggest refinements to current understandings of usability evaluations, and affirm the value to HCI in studying industry practitioners more deeply.


title: Learning Hash Codes via Hamming Distance Targets

abstract: We present a powerful new loss function and training scheme for learning binary hash codes with any differentiable model and similarity function. Our loss function improves over prior methods by using log likelihood loss on top of an accurate approximation for the probability that two inputs fall within a Hamming distance target. Our novel training scheme obtains a good estimate of the true gradient by better sampling inputs and evaluating loss terms between all pairs of inputs in each minibatch. To fully leverage the resulting hashes, we use multi-indexing. We demonstrate that these techniques provide large improvements to a similarity search tasks. We report the best results to date on competitive information retrieval tasks for ImageNet and SIFT 1M, improving MAP from 73% to 84% and reducing query cost by a factor of 2-8, respectively.


title: Power domination throttling

abstract: A power dominating set of a graph $G=(V,E)$ is a set $S\subset V$ that colors every vertex of $G$ according to the following rules: in the first timestep, every vertex in $N[S]$ becomes colored; in each subsequent timestep, every vertex which is the only non-colored neighbor of some colored vertex becomes colored. The power domination throttling number of $G$ is the minimum sum of the size of a power dominating set $S$ and the number of timesteps it takes $S$ to color the graph. In this paper, we determine the complexity of power domination throttling and give some tools for computing and bounding the power domination throttling number. Some of our results apply to very general variants of throttling and to other aspects of power domination.


title: The boundary value Minkowski problem for Weingarten curvatures

abstract: In this paper, we prove the existence of hypersurfaces in the Euclidean space with prescribed boundary and whose k-th Weingarten curvature equals a given function that depends on the normal of the hypersurface. The proof is based on the solvability of a fully nonlinear elliptic PDE. The required a priori estimates are established under the natural assumptions that the prescribed boundary is strictly convex and the prescribed function satisfies a Serrin type condition.


title: CNN-SVO: Improving the Mapping in Semi-Direct Visual Odometry Using Single-Image Depth Prediction

abstract: Reliable feature correspondence between frames is a critical step in visual odometry (VO) and visual simultaneous localization and mapping (V-SLAM) algorithms. In comparison with existing VO and V-SLAM algorithms, semi-direct visual odometry (SVO) has two main advantages that lead to state-of-the-art frame rate camera motion estimation: direct pixel correspondence and efficient implementation of probabilistic mapping method. This paper improves the SVO mapping by initializing the mean and the variance of the depth at a feature location according to the depth prediction from a single-image depth prediction network. By significantly reducing the depth uncertainty of the initialized map point (i.e., small variance centred about the depth prediction), the benefits are twofold: reliable feature correspondence between views and fast convergence to the true depth in order to create new map points. We evaluate our method with two outdoor datasets: KITTI dataset and Oxford Robotcar dataset. The experimental results indicate that the improved SVO mapping results in increased robustness and camera tracking accuracy.


title: Distributional Semantics Approach to Detect Intent in Twitter Conversations on Sexual Assaults

abstract: The recent surge in women reporting sexual assault and harassment (e.g., #metoo campaign) has highlighted a longstanding societal crisis. This injustice is partly due to a culture of discrediting women who report such crimes and also, rape myths (e.g., 'women lie about rape'). Social web can facilitate the further proliferation of deceptive beliefs and culture of rape myths through intentional messaging by malicious actors. This multidisciplinary study investigates Twitter posts related to sexual assaults and rape myths for characterizing the types of malicious intent, which leads to the beliefs on discrediting women and rape myths. Specifically, we first propose a novel malicious intent typology for social media using the guidance of social construction theory from policy literature that includes Accusational, Validational, or Sensational intent categories. We then present and evaluate a malicious intent classification model for a Twitter post using semantic features of the intent senses learned with the help of convolutional neural networks. Lastly, we analyze a Twitter dataset of four months using the intent classification model to study narrative contexts in which malicious intents are expressed and discuss their implications for gender violence policy design.


title: AI for Trustworthiness! Credible User Identification on Social Web for Disaster Response Agencies

abstract: Although social media provides a vibrant platform to discuss real-world events, the quantity of information generated can overwhelm decision making based on that information. By better understanding who is participating in information sharing, we can more effectively filter information as the event unfolds. Fine-grained understanding of credible sources can even help develop a trusted network of users for specific events or situations. Given the culture of relying on trusted actors for work practices in the humanitarian and disaster response domain, we propose to identify potential credible users as organizational and organizational-affiliated user accounts on social media in realtime for effective information collection and dissemination. Therefore, we examine social media using AI and Machine Learning methods during three types of humanitarian or disaster events and identify key actors responding to social media conversations as organization (business, group, or institution), organization-affiliated (individual with an organizational affiliation), and non-affiliated (individual without organizational affiliation) identities. We propose a credible user classification approach using a diverse set of social, activity, and descriptive representation features extracted from user profile metadata. Our extensive experiments showed a contrasting participation behavior of the user identities by their content practices, such as the use of higher authoritative content sharing by organization and organization-affiliated users. This study provides a direction for designing realtime credible content analytics systems for humanitarian and disaster response agencies.


title: Bayesian Policy Optimization for Model Uncertainty

abstract: Addressing uncertainty is critical for autonomous systems to robustly adapt to the real world. We formulate the problem of model uncertainty as a continuous Bayes-Adaptive Markov Decision Process (BAMDP), where an agent maintains a posterior distribution over latent model parameters given a history of observations and maximizes its expected long-term reward with respect to this belief distribution. Our algorithm, Bayesian Policy Optimization, builds on recent policy optimization algorithms to learn a universal policy that navigates the exploration-exploitation trade-off to maximize the Bayesian value function. To address challenges from discretizing the continuous latent parameter space, we propose a new policy network architecture that encodes the belief distribution independently from the observable state. Our method significantly outperforms algorithms that address model uncertainty without explicitly reasoning about belief distributions and is competitive with state-of-the-art Partially Observable Markov Decision Process solvers.


title: Convergence Rates for Empirical Estimation of Binary Classification Bounds

abstract: Bounding the best achievable error probability for binary classification problems is relevant to many applications including machine learning, signal processing, and information theory. Many bounds on the Bayes binary classification error rate depend on information divergences between the pair of class distributions. Recently, the Henze-Penrose (HP) divergence has been proposed for bounding classification error probability. We consider the problem of empirically estimating the HP-divergence from random samples. We derive a bound on the convergence rate for the Friedman-Rafsky (FR) estimator of the HP-divergence, which is related to a multivariate runs statistic for testing between two distributions. The FR estimator is derived from a multicolored Euclidean minimal spanning tree (MST) that spans the merged samples. We obtain a concentration inequality for the Friedman-Rafsky estimator of the Henze-Penrose divergence. We validate our results experimentally and illustrate their application to real datasets.


title: Approximating strange attractors and Lyapunov exponents of delay differential equations using Galerkin projections

abstract: Delay differential equations (DDEs) are infinite-dimensional systems, so even a scalar, unforced nonlinear DDE can exhibit chaos. Lyapunov exponents are indicators of chaos and can be computed by comparing the evolution of infinitesimally close trajectories. We convert DDEs into partial differential equations with nonlinear boundary conditions, then into ordinary differential equations (ODEs) using the Galerkin projection. The solution of the resulting ODEs approximates that of the original DDE system; for smooth solutions, the error decreases exponentially as the number of terms used in the Galerkin approximation increases. Examples demonstrate that the strange attractors and Lyapunov exponents of chaotic DDE solutions can be reliably approximated by a smaller number of ODEs using the proposed approach compared to the standard method-of-lines approach, leading to faster convergence and improved computational efficiency.


title: Fighting Against XSS Attacks: A Usability Evaluation of OWASP ESAPI Output Encoding

abstract: Cross Site Scripting (XSS) is one of the most critical vulnerabilities exist in web applications. XSS can be prevented by encoding untrusted data that are loaded into browser content of web applications. Security Application Programming Interfaces (APIs) such as OWASP ESAPI provide output encoding functionalities for programmers to use to protect their applications from XSS attacks. However, XSS still being ranked as one of the most critical vulnerabilities in web applications suggests that programmers are not effectively using those APIs to encode untrusted data. Therefore, we conducted an experimental study with 10 programmers where they attempted to fix XSS vulnerabilities of a web application using the output encoding functionality of OWASP ESAPI. Results revealed 3 types of mistakes that programmers made which resulted in them failing to fix the application by removing XSS vulnerabilities. We also identified 16 usability issues of OWASP ESAPI. We identified that some of these usability issues as the reason for mistakes that programmers made. Based on these results, we provided suggestions on how the usability of output encoding APIs should be improved to give a better experience to programmers.


title: Simultaneously Optimizing Weight and Quantizer of Ternary Neural Network using Truncated Gaussian Approximation

abstract: In the past years, Deep convolution neural network has achieved great success in many artificial intelligence applications. However, its enormous model size and massive computation cost have become the main obstacle for deployment of such powerful algorithm in the low power and resource-limited mobile systems. As the countermeasure to this problem, deep neural networks with ternarized weights (i.e. -1, 0, +1) have been widely explored to greatly reduce the model size and computational cost, with limited accuracy degradation. In this work, we propose a novel ternarized neural network training method which simultaneously optimizes both weights and quantizer during training, differentiating from prior works. Instead of fixed and uniform weight ternarization, we are the first to incorporate the thresholds of weight ternarization into a closed-form representation using the truncated Gaussian approximation, enabling simultaneous optimization of weights and quantizer through back-propagation training. With both of the first and last layer ternarized, the experiments on the ImageNet classification task show that our ternarized ResNet-18/34/50 only has 3.9/2.52/2.16% accuracy degradation in comparison to the full-precision counterparts.


title: Synchrotron radiation microtomography of brain hemisphere and spinal cord of a mouse model of multiple sclerosis revealed a correlation between capillary dilation and clinical score

abstract: Multiple sclerosis is a neurological disorder in which the myelin sheaths of axons are damaged by the immune response. We report here a three-dimensional structural analysis of brain and spinal cord tissues of a mouse model of multiple sclerosis, known as experimental autoimmune encephalomyelitis (EAE). EAE-induced mice were raised with or without administration of fingolimod, which is used in the treatment of multiple sclerosis. Brains and spinal cords dissected from the EAE mice were lyophilized so as to reconstitute the intrinsic contrast of tissue elements, such as axons, in X-ray images. Three-dimensional structures of the brain hemispheres and spinal cords of the EAE mice were visualized with synchrotron radiation microtomography. Microtomographic cross sections reconstructed from the X-ray images revealed dilation of capillary vessels and vacuolation in the spinal cord of the EAE mice. Vacuolation was also observed in the cerebellum, suggesting that the neuroinflammatory response progressed in the brain. The vessel networks and vacuolation lesions in the spinal cords were modelled by automatically tracing the three-dimensional image in order to analyze the tissue structures quantitatively. The results of the analysis indicated that the distribution of vacuolations was not uniform but three-dimensionally localized. The mean vessel diameter showed a linear correlation with the clinical score, indicating that vasodilation is relevant to paralysis severity in the disease model. We suggest that vasodilation and vacuolation are related with neurological symptoms of multiple sclerosis.


title: Application of Adaptive Optics for Illumination Stability in Precision Radial Velocity Measurements in Astronomical Spectroscopy

abstract: Adaptive optics (AO) have been used to correct wavefronts to achieve diffraction limited point spread functions in a broad range of optical applications, prominently ground-based astronomical telescopes operating in near infra-red. While most AO systems cannot provide diffraction-limited performance in the optical passband (400 nm - 900 nm), AO can improve image concentration, as well as both near and far field image stability, within an AO-fed spectrograph. Enhanced near and far field stability increase wavelength-scale stability in high dispersion spectrographs. In this work, we describe detailed modelling of the stability improvements achievable on extremely large telescopes. These improvements in performance may enable the mass measurement of Earth Twins by the precision radial velocity method, and the discovery of evidence of exobiotic activity in exoplanet atmospheres with the next generation of extremely large telescopes (ELTs). In this paper, we report on numerical simulations of the impact of AO on the performance of the GMT-Consortium Large Earth Finder (G-CLEF) instrument for the future Giant Magellan Telescope (GMT). The proximate cause of this study is to evaluate what improvements AO offer for exoplanet mass determination by the precision radial velocity (PRV) method and the discovery of biomarkers in exoplanet atmospheres. A modified AO system capable of achieving this improved stability even with changing conditions is proposed.


title: Large batch size training of neural networks with adversarial training and second-order information

abstract: The most straightforward method to accelerate Stochastic Gradient Descent (SGD) computation is to distribute the randomly selected batch of inputs over multiple processors. To keep the distributed processors fully utilized requires commensurately growing the batch size. However, large batch training often leads to poorer generalization. A recently proposed solution for this problem is to use adaptive batch sizes in SGD. In this case, one starts with a small number of processes and scales the processes as training progresses. Two major challenges with this approach are (i) that dynamically resizing the cluster can add non-trivial overhead, in part since it is currently not supported, and (ii) that the overall speed up is limited by the initial phase with smaller batches. In this work, we address both challenges by developing a new adaptive batch size framework, with autoscaling based on the Ray framework. This allows very efficient elastic scaling with negligible resizing overhead (0.32\% of time for ResNet18 ImageNet training). Furthermore, we propose a new adaptive batch size training scheme using second order methods and adversarial training. These enable increasing batch sizes earlier during training, which leads to better training time. We extensively evaluate our method on Cifar-10/100, SVHN, TinyImageNet, and ImageNet datasets, using multiple neural networks, including ResNets and smaller networks such as SqueezeNext. Our method exceeds the performance of existing solutions in terms of both accuracy and the number of SGD iterations (up to 1\% and $5\times$, respectively). Importantly, this is achieved without any additional hyper-parameter tuning to tailor our method in any of these experiments.


title: Extremely large telescopes for complex stellar populations around the Galactic centre

abstract: The Galactic centre and its surrounding space are important in studying galaxy-scale evolution, and stellar populations therein are expected to have imprints of the long-term evolution. Interstellar extinction, however, severely limits optical observations, thereby requiring infrared observations. In addition, many systems from those in the proximity of the central black hole to foreground objects in the disc overlap in each sightline, which complicates the interpretation of observations of a wide variety of objects. We discuss some important issues concerning the central regions, particularly the Galactic bulge and the Nuclear Stellar Disc (also known as the Central Molecular Zone). An obvious advantage of Extremely Large Telescopes (ELTs) is the deeper limiting magnitudes, but we emphasise the importance of the synergy between the data of deep ELTs and other observational data (e.g. astrometric measurements and the detection of interstellar absorption lines) in order to disentangle the complex stellar populations.


title: Principal bundles of open groupoids

abstract: Taking into account the correspondence between open groupoids and their quantales, we establish a bijective correspondence between the principal $G$-bundles whose left projection is an open surjection and the principal $\opens(G)$-locales for any open groupoid $G$ and its groupoid quantale $\opens(G)$, the latter requires a generalization of the supported modules for groupoid quantales.


title: On Pointwise Products of Elliptic Eigenfunctions

abstract: We consider eigenfunctions of Schr\"odinger operators on a $d-$dimensional bounded domain $\Omega$ (or a $d-$dimensional compact manifold $\Omega$) with Dirichlet conditions. These operators give rise to a sequence of eigenfunctions $(\phi_n)_{n \in \mathbb{N}}$. We study the subspace of all pointwise products $$ A_n = \mbox{span} \left\{ \phi_i(x) \phi_j(x): 1 \leq i,j \leq n\right\} \subseteq L^2(\Omega).$$ Clearly, that vector space has dimension $\mbox{dim}(A_n) = n(n+1)/2$. We prove that products $\phi_i \phi_j$ of eigenfunctions are simple in a certain sense: for any $\varepsilon > 0$, there exists a low-dimensional vector space $B_n$ that almost contains all products. More precisely, denoting the orthogonal projection $\Pi_{B_n}:L^2(\Omega) \rightarrow B_n$, we have $$ \forall~1 \leq i,j \leq n~ \qquad \|\phi_i\phi_j - \Pi_{B_n}( \phi_i \phi_j) \|_{L^2} \leq \varepsilon$$ and the size of the space $\mbox{dim}(B_n)$ is relatively small $$ \mbox{dim}(B_n) \lesssim \left( \frac{1}{\varepsilon} \max_{1 \leq i \leq n} \|\phi_i\|_{L^{\infty}} \right)^d n.$$ In the generic delocalized setting, this bound grows linearly up to logarithmic factors: pointwise products of eigenfunctions are low-rank. This has implications, among other things, for the validity of fast algorithms in electronic structure computations.


title: Impacts of the New Carbon Fusion Cross Sections on Type Ia Supernovae

abstract: Type Ia supernovae (SNe Ia) are thought to be thermonuclear explosion of white dwarfs (WDs). Their progenitors are not well understood. One popular scenario is the double-degenerate (DD) scenario, which attributes SNe Ia to WD-WD binary mergers. The fates of the WD mergers depend on the rate of $^{12}$C+$^{12}$C reaction. Recently, the $^{12}$C+$^{12}$C cross sections have been measured and the analysis of the data using the Trojan Horse Method suggested that the astrophysical reaction rate is larger than conventional rates at astrophysical temperatures due to possible resonances. The resonance contribution results in a decrease of the carbon burning ignition temperature. Therefore accretion induced collapse occurs more easily and increases the birthrate of Galactic neutron stars with the contribution of the DD scenario to the SNe Ia rate becoming even smaller.


title: Topology of complexity one quotients

abstract: We describe of the topology of the geometric quotients of 2n dimensional compact connected symplectic manifolds with n-1 dimensional torus actions. When the isotropy weights at each fixed point are in general position, the quotient is homeomorphic to a sphere.


title: Quantum-limited estimation of the axial separation of two incoherent point sources

abstract: Improving axial resolution is crucial for three-dimensional optical imaging systems. Here we present a scheme of axial superresolution for two incoherent point sources based on spatial mode demultiplexing. A radial mode sorter is used to losslessly decompose the optical fields into a radial mode basis set to extract the phase information associated with the axial positions of the point sources. We show theoretically and experimentally that, in the limit of a zero axial separation, our scheme allows for reaching the quantum Cram\'er-Rao lower bound and thus can be considered as one of the optimal measurement methods. Unlike other superresolution schemes, this scheme does not require neither activation of fluorophores nor sophisticated stabilization control. Moreover, it is applicable to the localization of a single point source in the axial direction. Our demonstration can be useful to a variety of applications such as far-field fluorescence microscopy.


title: Approximation of probability density functions for PDEs with random parameters using truncated series expansions

abstract: The probability density function (PDF) of a random variable associated with the solution of a partial differential equation (PDE) with random parameters is approximated using a truncated series expansion. The random PDE is solved using two stochastic finite element methods, Monte Carlo sampling and the stochastic Galerkin method with global polynomials. The random variable is a functional of the solution of the random PDE, such as the average over the physical domain. The truncated series are obtained considering a finite number of terms in the Gram-Charlier or Edgeworth series expansions. These expansions approximate the PDF of a random variable in terms of another PDF, and involve coefficients that are functions of the known cumulants of the random variable. To the best of our knowledge, their use in the framework of PDEs with random parameters has not yet been explored.


title: Density-Matrix Simulation of Logical Qubit using 3-qubit Quantum Error Correction Code

abstract: Fault-tolerant quantum computing demands many qubits with long lifetimes to conduct accurate quantum gate operations. However, external noise limits the computing time of physical qubits. Quantum error correction codes may extend such limits, but imperfect gate operations introduce errors to the correction procedure as well. The additional gate operations required due to the physical layout of qubits exacerbate the situation. Here, we use density-matrix simulations to investigate the performance change of logical qubits according to quantum error correction codes and qubit layouts and the expected performance of logical qubits with gate operation time and gate error rates. Considering current qubit technology, the small quantum error correction codes are chosen. Assuming 0.1% gate error probability, a logical qubit encoded by a 5-qubit quantum error correction code is expected to have a fidelity 0.25 higher than its physical counterpart.


title: Simulating the Chiral Magnetic Wave in a Box System

abstract: The chiral magnetic wave from the interplay between the chiral magnetic effect and the chiral separation effect is simulated in a box system with the periodic boundary condition based on the chiral kinetic equations of motion. Simulation results are compared with available limits from theoretical derivations, and effects of the temperature, the magnetic field, and the specific shear viscosity on the key properties of the chiral magnetic wave are discussed. Our study serves as a baseline for further simulations of chiral anomalies in relativistic heavy-ion collisions.


title: Unobstructedness of hyperk\"ahler twistor spaces

abstract: A family of irreducible holomorphic symplectic (ihs) manifolds over the complex projective line has unobstructed deformations if its period map is an embedding. This applies in particular to twistor spaces of ihs manifolds. Moreover, a family of ihs manifolds over a subspace of the period domain extends to a universal family over an open neighborhood in the period domain.


title: Reinforcement Learning with Perturbed Rewards

abstract: Recent studies have shown that reinforcement learning (RL) models are vulnerable in various noisy scenarios. For instance, the observed reward channel is often subject to noise in practice (e.g., when rewards are collected through sensors), and is therefore not credible. In addition, for applications such as robotics, a deep reinforcement learning (DRL) algorithm can be manipulated to produce arbitrary errors by receiving corrupted rewards. In this paper, we consider noisy RL problems with perturbed rewards, which can be approximated with a confusion matrix. We develop a robust RL framework that enables agents to learn in noisy environments where only perturbed rewards are observed. Our solution framework builds on existing RL/DRL algorithms and firstly addresses the biased noisy reward setting without any assumptions on the true distribution (e.g., zero-mean Gaussian noise as made in previous works). The core ideas of our solution include estimating a reward confusion matrix and defining a set of unbiased surrogate rewards. We prove the convergence and sample complexity of our approach. Extensive experiments on different DRL platforms show that trained policies based on our estimated surrogate reward can achieve higher expected rewards, and converge faster than existing baselines. For instance, the state-of-the-art PPO algorithm is able to obtain 84.6% and 80.8% improvements on average score for five Atari games, with error rates as 10% and 30% respectively.


title: Targeted Assembly and Synchronization of Self-Spinning Microgears

abstract: Self-assembly is the autonomous organization of components into patterns or structures: an essential ingredient of biology and a desired route to complex organization. At equilibrium, the structure is encoded through specific interactions, at an unfavorable entropic cost for the system. An alternative approach, widely used by Nature, uses energy input to bypass the entropy bottleneck and develop features otherwise impossible at equilibrium. Dissipative building blocks that inject energy locally were made available by recent advance in colloidal science but have not been used to control self-assembly. Here we show the robust formation of self-powered rotors and dynamical superstructures from active particles and harness non-equilibrium phoretic phenomena to tailor interactions and direct self-assembly. We use a photoactive component that consumes fuel, hematite, to devise phototactic microswimmers that form self-spinning microgears following spatiotemporal light patterns. The gears are coupled via their chemical clouds and constitute the elementary bricks of synchronized superstructures, which autonomously regulate their dynamics. The results are quantitatively rationalized on the basis of a stochastic description of diffusio-phoretic oscillators dynamically coupled by chemical gradients to form directional interactions. Our findings demonstrate that non-equilibrium phenomena can be harnessed to shape interactions and program hierarchical constructions. It lays the groundwork for the self-assembly of dynamical architectures and synchronized micro-machinery.


title: On the Betti numbers of Springer fibers for classical types

abstract: For a Weyl group $W$ of classical type, we present a formula to calculate the restriction of (graded) Springer representations of $W$ to a maximal parabolic subgroup $W'$ where the types of $W$ and $W'$ are in the same series. As a result, we obtain recursive formulas for the Betti numbers of Springer fibers for classical types.


title: Real-Time Planning with Multi-Fidelity Models for Agile Flights in Unknown Environments

abstract: Autonomous navigation through unknown environments is a challenging task that entails real-time localization, perception, planning, and control. UAVs with this capability have begun to emerge in the literature with advances in lightweight sensing and computing. Although the planning methodologies vary from platform to platform, many algorithms adopt a hierarchical planning architecture where a slow, low-fidelity global planner guides a fast, high-fidelity local planner. However, in unknown environments, this approach can lead to erratic or unstable behavior due to the interaction between the global planner, whose solution is changing constantly, and the local planner; a consequence of not capturing higher-order dynamics in the global plan. This work proposes a planning framework in which multi-fidelity models are used to reduce the discrepancy between the local and global planner. Our approach uses high-, medium-, and low-fidelity models to compose a path that captures higher-order dynamics while remaining computationally tractable. In addition, we address the interaction between a fast planner and a slower mapper by considering the sensor data not yet fused into the map during the collision check. This novel mapping and planning framework for agile flights is validated in simulation and hardware experiments, showing replanning times of 5-40 ms in cluttered environments.


title: Towards Online Learning from Corrective Demonstrations

abstract: Robots operating in real-world human environments will likely encounter task execution failures. To address this, we would like to allow co-present humans to refine the robot's task model as errors are encountered. Existing approaches to task model modification require reasoning over the entire dataset and model, limiting the rate of corrective updates. We introduce the State-Indexed Task Updates (SITU) algorithm to efficiently incorporate corrective demonstrations into an existing task model by iteratively making local updates that only require reasoning over a small subset of the model. In future work, we will evaluate this approach with a user study.


title: Heterogeneous Replica for Query on Cassandra

abstract: Cassandra is a popular structured storage system with high-performance, scalability and high availability, and is usually used to store data that has some sortable attributes. When deploying and configuring Cassandra, it is important to design a suitable schema of column families for accelerating the target queries. However, one schema is only suitable for a part of queries, and leaves other queries with high latency. In this paper, we propose a new replica mechanism, called heterogeneous replica, to reduce the query latency greatly while ensuring high write throughput and data recovery. With this replica mechanism, different replica has the same dataset while having different serialization on disk. By implementing the heterogeneous replica mechanism on Cassandra, we show that the read performance of Cassandra can be improved by two orders of magnitude with TPC-H data set.


title: Korselt Rational Bases of Prime Powers

abstract: Let $N$ be a positive integer, $\mathbb{A}$ be a subset of $\mathbb{Q}$ and $\alpha=\dfrac{\alpha_{1}}{\alpha_{2}}\in \mathbb{A}\setminus \{0,N\}$. $N$ is called an $\alpha$-Korselt number (equivalently $\alpha$ is said an $N$-Korselt base) if $\alpha_{2}p-\alpha_{1}$ divides $\alpha_{2}N-\alpha_{1}$ for every prime divisor $p$ of $N$. By the Korselt set of $N$ over $\mathbb{A}$, we mean the set $\mathbb{A}$-$\mathcal{KS}(N)$ of all $\alpha\in \mathbb{A}\setminus \{0,N\}$ such that $N$ is an $\alpha$-Korselt number. In this paper we determine explicitly for a given prime number $q$ and an integer $l\in \mathbb{N}\setminus \{0,1\}$, the set $\mathbb{Q}$-$\mathcal{KS}(q^{l})$ and we establish some connections between the $q^{l}$-Korselt bases in $\mathbb{Q}$ and others in $\mathbb{Z}$. The case of $\mathbb{A}=\mathbb{Q}\cap[-1,1[$ is studied where we prove that $(\mathbb{Q}\cap[-1,1[)$-$\mathcal{KS}(q^{l})$ is empty if and only if $l=2$. Moreover, we show that each nonzero rational $\alpha$ is an $N$-Korselt base for infinitely many numbers $N=q^{l}$ where $q$ is a prime number and $l\in\mathbb{N}$.


title: Loop corrections to dark matter direct detection in a pseudoscalar mediator dark matter model

abstract: If dark matter (DM) is a fermion and its interactions with the standard model particles are mediated by pseudoscalar particles, the tree-level amplitude for the DM-nucleon elastic scattering is suppressed by the momentum transfer in the non-relativistic limit. At the loop level, on the other hand, the spin-independent contribution to the cross section appears without such suppression. Thus, the loop corrections are essential to discuss the sensitivities of the direct detection experiments for the model prediction. The one-loop corrections were investigated in the previous works. However, the two-loop diagrams give the leading order contribution to the DM-gluon effective operator ($\bar{\chi} \chi G^{a}_{\mu \nu} G^{a \mu \nu} $) and have not been correctly evaluated yet. Moreover, some interaction terms which affect the scattering cross section were overlooked. In this paper, we show the cross section obtained by the improved analysis and discuss the region where the cross section becomes large.


title: Stabilizer Slicing: Coherent Error Cancellations in LDPC Codes

abstract: Coherent errors are a dominant noise process in many quantum computing architectures. Unlike stochastic errors, these errors can combine constructively and grow into highly detrimental overrotations. To combat this, we introduce a simple technique for suppressing systematic coherent errors in low-density parity-check (LDPC) stabilizer codes, which we call stabilizer slicing. The essential idea is to slice low-weight stabilizers into two equally-weighted Pauli operators and then apply them by rotating in opposite directions, causing their overrotations to interfere destructively on the logical subspace. With access to native gates generated by 3-body Hamiltonians, we can completely eliminate purely coherent overrotation errors, and for overrotation noise of 0.99 unitarity we achieve a 135-fold improvement in the logical error rate of Surface-17. For more conventional 2-body ion trap gates, we observe an 89-fold improvement for Bacon-Shor-13 with purely coherent errors which should be testable in near-term fault-tolerance experiments. This second scheme takes advantage of the prepared gauge degrees of freedom, and to our knowledge is the first example in which the state of the gauge directly affects the robustness of a code's memory. This work demonstrates that coherent noise is preferable to stochastic noise within certain code and gate implementations when the coherence is utilized effectively.


title: The $\mathbb{Q}$-Korselt Set of $\mathrm{pq}$

abstract: Let $N$ be a positive integer, $\mathbb{A}$ be a nonempty subset of $\mathbb{Q}$ and $\alpha=\dfrac{\alpha_{1}}{\alpha_{2}}\in \mathbb{A}\setminus \{0,N\}$. $\alpha$ is called an $N$-Korselt base (equivalently $N$ is said an $\alpha$-Korselt number) if $\alpha_{2}p-\alpha_{1}$ is a divisor of $\alpha_{2}N-\alpha_{1}$ for every prime $p$ dividing $N$. The set of all Korselt bases of $N$ in $\mathbb{A}$ is called the $\mathbb{A}$-Korselt set of $N$ and is simply denoted by $\mathbb{A}$-$\mathcal{KS}(N)$. Let $p$ and $q$ be two distinct prime numbers. In this paper, we study the $\mathbb{Q}$-Korselt bases of $pq$, where we give in detail how to provide $\mathbb{Q}$-$\mathcal{KS}(pq)$. Consequently, we finish the incomplete characterization of the Korselt set of $pq$ over $\mathbb{Z}$ given in [4], by supplying the set $\mathbb{Z}$-$\mathcal{KS}(pq)$ when $q <2p$.


title: Implementing the Lexicographic Maxmin Bargaining Solution

abstract: There has been much work on exhibiting mechanisms that implement various bargaining solutions, in particular, the Kalai-Smorodinsky solution \cite{moulin1984implementing} and the Nash Bargaining solution. Another well-known and axiomatically well-studied solution is the lexicographic maxmin solution. However, there is no mechanism known for its implementation. To fill this gap, we construct a mechanism that implements the lexicographic maxmin solution as the unique subgame perfect equilibrium outcome in the n-player setting. As is standard in the literature on implementation of bargaining solutions, we use the assumption that any player can grab the entire surplus. Our mechanism consists of a binary game tree, with each node corresponding to a subgame where the players are allowed to choose between two outcomes. We characterize novel combinatorial properties of the lexicographic maxmin solution which are crucial to the design of our mechanism.


title: Nondegenerate spheres in four dimensions

abstract: Non-degeneracy was first defined for hyperplanes by Elekes-T\'oth, and later extended to spheres by Apfelbaum-Sharir: given a set $P$ of $m$ points in $\mathbb{R}^d$ and some $\beta\in(0,1)$, a $(d-1)$-dimensional sphere (or a $(d-1)$-sphere) $S$ in $\mathbb{R}^d$ is called $\beta$-nondegenerate with respect to $P$ if $S$ does not contain a proper subsphere $S'$ such that $|S'\cap P|\geq \beta|S\cap P|$. Apfelbaum-Sharir found an upper bound for the number of incidences between points and nondegenerate spheres in three dimensions, which was recently used by Zahl to obtain the best known bound for the unit distance problem in three dimensions. In this paper, we show that the number of incidences between $m$ points and $n$ $\beta$-nondegenerate 3-spheres in $\mathbb{R}^4$ is $O_{\beta,\varepsilon}\left(m^{\frac{15}{19}+\varepsilon} n^{\frac{16}{19}}+mn^{\frac{2}{3}}\right)$. As a consequence, we obtain a bound of $O_{\varepsilon}(n^{2+4/11+\varepsilon})$ on the number of similar triangles formed by $n$ points in $\mathbb{R}^4$, an improvement over the previously best known bound $O(n^{2+2/5})$. While proving this, we find it convenient to work with a more general definition of nondegeneracy: a bipartite graph $G=(P,Q)$ is called $\beta$-nondegenerate if $|N(q_1)\cap N(q_2)|<\beta |N(q_1)|$ for any two distinct vertices $q_1,q_2\in Q$; here $N(q)$ denotes the set of neighbors of $q$ and $\beta$ is some positive constant less than 1. A $\beta$-nondegenerate graph can have up to $\Theta(|P||Q|)$ edges without any restriction, but must have much fewer edges if the graph is semi-algebraic or has bounded VC-dimension. We show that previous results for planes and spheres in three dimensions still hold under this new definition, and so does our new bound for spheres in four dimensions.


title: Topological insulators vs. topological Dirac semimetals in honeycomb compounds

abstract: Intriguing physical property of materials stems from their chemical constituent whereas the connection between them is often not clear. Here, we uncover a general chemical classification for the two quantum phases in the honeycomb ABX structure--topological insulator (TI) and topological Dirac semimetal (TDSM). First, we find among the 816 (existing as well as hypothetical) calculated compounds, 160 TI's (none were noted before), 96 TDSM's, 282 normal insulators (NI's), and 278 metals. Second, based on this classification, we have distilled a simple chemical regularity based on compound formulae for the selectivity between TI and TDSM: The ABX compounds that are TDSM have B atoms (part of the BX honeycomb layers) that come from the Periodic Table columns XI (Cu, Ag, Au) or XII (Zn, Cd, Hg), or Mg (group II), whereas the ABX compounds whose B atoms come from columns I (Li, Na, K, Rb, Cs) or II (Ca, Sr, Ba) are TI's. Third, focusing on the ABX Bismide compounds that are thermodynamically stable, we find a structural motif that delivers topological insulation and stability at the same time. This study opens the way to simultaneously design new topological materials based on the compositional rules indicated here.


title: A ${\mathbb Z}_2$-index of symmetry protected topological phases with time reversal symmetry for quantum spin chains

abstract: We introduce a ${\mathbb Z}_2$-index for time reversal invariant Hamiltonians with unique gapped ground state on quantum spin chains. We show this is an invariant of a $C^1$-classification of gapped Hamiltonians.

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


title: PhotoSafer: Content-Based and Context-Aware Private Photo Protection for Smartphones

abstract: Nowadays many people store photos in smartphones. Many of the photos contain sensitive, private information, such as a photocopy of driver's license and credit card. An arising privacy concern is with the unauthorized accesses to such private photos by installed apps. Coarse-grained access control systems such as the Android permission system offer all-or-nothing access to photos stored on smartphones, and users are unaware of the exact behavior of installed apps. Our analysis finds that 82% of the top 200 free apps in a popular Android app store have complete access to stored photos and network on a user's smartphone, which indicates possible private photo leakage. In addition, our user survey reveals that 87.5% of the 112 respondents are not aware that certain apps can access their photos without informing users, and all the respondents believe that the stored photos on their smartphones contain different types of private information. Hence, we propose PhotoSafer, a content-based, context-aware private photo protection system for Android phones. PhotoSafer can detect private photos based on photo content with a well-trained deep convolutional neural network, and control access to photos based on system status (e.g., screen locked or not) and app-running status (e.g., app in the background). Evaluations demonstrate that PhotoSafer can accurately identify private photos in real time. The efficacy and efficiency of the implemented prototype system show the potential for practical use.


title: Phenomenology in type-I minimally modified gravity

abstract: We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.


title: Privacy-Preserving Outsourcing of Large-Scale Nonlinear Programming to the Cloud

abstract: The increasing massive data generated by various sources has given birth to big data analytics. Solving large-scale nonlinear programming problems (NLPs) is one important big data analytics task that has applications in many domains such as transport and logistics. However, NLPs are usually too computationally expensive for resource-constrained users. Fortunately, cloud computing provides an alternative and economical service for resource-constrained users to outsource their computation tasks to the cloud. However, one major concern with outsourcing NLPs is the leakage of user's private information contained in NLP formulations and results. Although much work has been done on privacy-preserving outsourcing of computation tasks, little attention has been paid to NLPs. In this paper, we for the first time investigate secure outsourcing of general large-scale NLPs with nonlinear constraints. A secure and efficient transformation scheme at the user side is proposed to protect user's private information; at the cloud side, generalized reduced gradient method is applied to effectively solve the transformed large-scale NLPs. The proposed protocol is implemented on a cloud computing testbed. Experimental evaluations demonstrate that significant time can be saved for users and the proposed mechanism has the potential for practical use.


title: A Unified Framework for Clustering Constrained Data without Locality Property

abstract: In this paper, we consider a class of constrained clustering problems of points in $\mathbb{R}^{d}$, where $d$ could be rather high. A common feature of these problems is that their optimal clusterings no longer have the locality property (due to the additional constraints), which is a key property required by many algorithms for their unconstrained counterparts. To overcome the difficulty caused by the loss of locality, we present in this paper a unified framework, called {\em Peeling-and-Enclosing (PnE)}, to iteratively solve two variants of the constrained clustering problems, {\em constrained $k$-means clustering} ($k$-CMeans) and {\em constrained $k$-median clustering} ($k$-CMedian). Our framework is based on two standalone geometric techniques, called {\em Simplex Lemma} and {\em Weaker Simplex Lemma}, for $k$-CMeans and $k$-CMedian, respectively. The simplex lemma (or weaker simplex lemma) enables us to efficiently approximate the mean (or median) point of an unknown set of points by searching a small-size grid, independent of the dimensionality of the space, in a simplex (or the surrounding region of a simplex), and thus can be used to handle high dimensional data. If $k$ and $\frac{1}{\epsilon}$ are fixed numbers, our framework generates, in nearly linear time ({\em i.e.,} $O(n(\log n)^{k+1}d)$), $O((\log n)^{k})$ $k$-tuple candidates for the $k$ mean or median points, and one of them induces a $(1+\epsilon)$-approximation for $k$-CMeans or $k$-CMedian, where $n$ is the number of points. Combining this unified framework with a problem-specific selection algorithm (which determines the best $k$-tuple candidate), we obtain a $(1+\epsilon)$-approximation for each of the constrained clustering problems. We expect that our technique will be applicable to other constrained clustering problems without locality.


title: Electron doping of proposed quantum spin liquid kagom\'e Zn-Cu hydroxyl-halides produces localized states in the band gap

abstract: Carrier doping of quantum spin liquids is a long-proposed route to the emergence of high-temperature superconductivity. Electrochemical intercalation in kagome hydroxyl-halide materials shows that samples remain insulating across a wide range of electron counts. Here we demonstrate through first-principles density functional calculations corrected for self-interaction the mechanism by which electrons remain localized in various Zn-Cu hydroxyl-halides, independently of the chemical identity of the dopant - the formation of polaronic states with attendant lattice displacements and a dramatic narrowing of bandwidth upon electron addition. The same theoretical method applied to electron doping in cuprate Nd2CuO4 correctly produces a metallic state when the initially formed polaron dissolves into an extended state. Our general findings explain the insulating behavior in a wide range of doped quantum magnets and demonstrate that new quantum spin liquid host materials are needed to realize metallicity borne of a spin liquid.


title: Improved Parallel Rabin-Karp Algorithm Using Compute Unified Device Architecture

abstract: String matching algorithms are among one of the most widely used algorithms in computer science. Traditional string matching algorithms efficiency of underlaying string matching algorithm will greatly increase the efficiency of any application. In recent years, Graphics processing units are emerged as highly parallel processor. They out perform best of the central processing units in scientific computation power. By combining recent advancement in graphics processing units with string matching algorithms will allows to speed up process of string matching. In this paper we proposed modified parallel version of Rabin-Karp algorithm using graphics processing unit. Based on that, result of CPU as well as parallel GPU implementations are compared for evaluating effect of varying number of threads, cores, file size as well as pattern size.


title: Quantum optomechanics of a two-dimensional atomic array

abstract: We demonstrate that a two-dimensional (2D) atomic array can be used as a novel platform for quantum optomechanics. Such arrays feature both nearly-perfect reflectivity and ultra-light mass, leading to significantly-enhanced optomechanical phenomena. Considering the collective atom-array motion under continuous laser illumination, we study the nonlinear optical response of the array. We find that the spectrum of light scattered by the array develops multiple sidebands, corresponding to collective mechanical resonances, and exhibits nearly perfect quantum-noise squeezing. Possible extensions and applications for quantum nonlinear optomechanics are discussed.


title: Decentralized Accelerated Gradient Methods With Increasing Penalty Parameters

abstract: In this paper, we study the communication and (sub)gradient computation costs in distributed optimization and give a sharp complexity analysis for the proposed distributed accelerated gradient methods. We present two algorithms based on the framework of the accelerated penalty method with increasing penalty parameters. Our first algorithm is for smooth distributed optimization and it obtains the near optimal $O\left(\sqrt{\frac{L}{\epsilon(1-\sigma_2(W))}}\log\frac{1}{\epsilon}\right)$ communication complexity and the optimal $O\left(\sqrt{\frac{L}{\epsilon}}\right)$ gradient computation complexity for $L$-smooth convex problems, where $\sigma_2(W)$ denotes the second largest singular value of the weight matrix $W$ associated to the network and $\epsilon$ is the target accuracy. When the problem is $\mu$-strongly convex and $L$-smooth, our algorithm has the near optimal $O\left(\sqrt{\frac{L}{\mu(1-\sigma_2(W))}}\log^2\frac{1}{\epsilon}\right)$ complexity for communications and the optimal $O\left(\sqrt{\frac{L}{\mu}}\log\frac{1}{\epsilon}\right)$ complexity for gradient computations. Our communication complexities are only worse by a factor of $\left(\log\frac{1}{\epsilon}\right)$ than the lower bounds for the smooth distributed optimization. %As far as we know, our method is the first to achieve both communication and gradient computation lower bounds up to an extra logarithm factor for smooth distributed optimization. Our second algorithm is designed for non-smooth distributed optimization and it achieves both the optimal $O\left(\frac{1}{\epsilon\sqrt{1-\sigma_2(W)}}\right)$ communication complexity and $O\left(\frac{1}{\epsilon^2}\right)$ subgradient computation complexity, which match the communication and subgradient computation complexity lower bounds for non-smooth distributed optimization.


title: ChainQueen: A Real-Time Differentiable Physical Simulator for Soft Robotics

abstract: Physical simulators have been widely used in robot planning and control. Among them, differentiable simulators are particularly favored, as they can be incorporated into gradient-based optimization algorithms that are efficient in solving inverse problems such as optimal control and motion planning. Simulating deformable objects is, however, more challenging compared to rigid body dynamics. The underlying physical laws of deformable objects are more complex, and the resulting systems have orders of magnitude more degrees of freedom and therefore they are significantly more computationally expensive to simulate. Computing gradients with respect to physical design or controller parameters is typically even more computationally challenging. In this paper, we propose a real-time, differentiable hybrid Lagrangian-Eulerian physical simulator for deformable objects, ChainQueen, based on the Moving Least Squares Material Point Method (MLS-MPM). MLS-MPM can simulate deformable objects including contact and can be seamlessly incorporated into inference, control and co-design systems. We demonstrate that our simulator achieves high precision in both forward simulation and backward gradient computation. We have successfully employed it in a diverse set of control tasks for soft robots, including problems with nearly 3,000 decision variables.


title: A Fourier-Bessel method with a regularization strategy for the boundary value problems of the Helmholtz equation

abstract: This paper is concerned with the Fourier-Bessel method for the boundary value problems of the Helmholtz equation in a smooth simply connected domain. Based on the denseness of Fourier-Bessel functions, the problem can be approximated by determining the unknown coefficients in the linear combination. By the boundary conditions, an operator equation can be obtained. We derive a lower bound for the smallest singular value of the operator, and obtain a stability and convergence result for the regularized solution with a suitable choice of the regularization parameter. Numerical experiments are also presented to show the effectiveness of the proposed method.


title: CMOS-Integrated Diamond Nitrogen-Vacancy Quantum Sensor

abstract: The nitrogen vacancy (NV) center in diamond has emerged as a leading solid-state quantum sensor for applications including magnetometry, electrometry, thermometry, and chemical sensing. However, an outstanding challenge for practical applications is that existing NV-based sensing techniques require bulky and discrete instruments for spin control and detection. Here, we address this challenge by integrating NV based quantum sensing with complementary metal-oxide-semiconductor (CMOS) technology. Through tailored CMOS-integrated microwave generation and photodetection, this work dramatically reduces the instrumentation footprint for quantum magnetometry and thermometry. This hybrid diamond-CMOS integration enables an ultra-compact and scalable platform for quantum sensing and quantum information processing.


title: High-mobility indirect excitons in wide single quantum well

abstract: Indirect excitons (IXs) are bound pairs of electrons and holes confined in spatially separated layers. We present wide single quantum well (WSQW) heterostructures with high IX mobility, spectrally narrow IX emission, voltage-controllable IX energy, and long and voltage-controllable IX lifetime. This set of properties shows that WSQW heterostructures provide an advanced platform both for studying basic properties of IXs in low-disorder environments and for the development of high mobility excitonic devices.


title: Reducing Subspaces of de Branges-Rovnyak Spaces

abstract: For $b\in H^\infty_1$, the closed unit ball of $H^\infty$, the de Branges-Rovnyak spaces $\mathcal{H}(b)$ is a Hilbert space contractively contained in the Hardy space $H^2$ that is invariant by the backward shift operator $S^*$. We consider the reducing subspaces of the operator $S^{*2}|_{\mathcal{H}(b)}$. When $b$ is an inner function, $S^{*2}|_{\mathcal{H}(b)}$ is a truncated Toepltiz operator and its reducibility was characterized by Douglas and Foias using model theory. We use another approach to extend their result to the case where $b$ is extreme. We prove that if $b$ is extreme but not inner, then $S^{*2}|_{\mathcal{H}(b)}$ is reducible if and only if $b$ is even or odd, and describe the structure of reducing subspaces.


title: Super-Resolution Blind Channel-and-Signal Estimation for Massive MIMO with One-Dimensional Antenna Array

abstract: In this paper, we study blind channel-and-signal estimation by exploiting the burst-sparse structure of angular-domain propagation channels in massive MIMO systems. The state-of-the-art approach utilizes the structured channel sparsity by sampling the angular-domain channel representation with a uniform angle-sampling grid, a.k.a. virtual channel representation. However, this approach is only applicable to uniform linear arrays and may cause a substantial performance loss due to the mismatch between the virtual representation and the true angle information. To tackle these challenges, we propose a sparse channel representation with a super-resolution sampling grid and a hidden Markovian support. Based on this, we develop a novel approximate inference based blind estimation algorithm to estimate the channel and the user signals simultaneously, with emphasis on the adoption of the expectation-maximization method to learn the angle information. Furthermore, we demonstrate the low-complexity implementation of our algorithm, making use of factor graph and message passing principles to compute the marginal posteriors. Numerical results show that our proposed method significantly reduces the estimation error compared to the state-of-the-art approach under various settings, which verifies the efficiency and robustness of our method.


title: New Families of Large Band Gap 2D Topological Insulators in Ethynyl-Derivative Functionalized Compounds

abstract: The search for large band gap systems with dissipationless edge states is essential to developing materials that function under a wide range of temperatures. Two-dimensional (2D) topological insulators (TIs) have recently attracted significant attention due to their dissipationless transport, robust properties and excellent compatibility with device integration. However, a major barrier of 2D TIs is their small bulk band gap, which allows for applications only in extremely low temperatures. In this work, first principle calculations were used to analyze the geometric, electronic, and topological properties of PbC2X and BiC2X (X = H, Cl, F, Br, I) compounds. The band gap values are remarkably large, ranging from 0.79eV to 0.99eV. The nanoribbons of these compounds exhibited nontrivial topological order in the simulation, thus proving ethynyl-derivative functionalized Pb and Bi films to be new classes of giant band gap 2D TIs. In addition, these findings indicate that chemical functionalization with ethynyl-derivatives is an effective method to tune the band gap and preserve the nontrivial topological order. These novel materials that are applicable at both room temperature and high temperatures open the door to a new generation of electronics.


title: Feature Selection Approach with Missing Values Conducted for Statistical Learning: A Case Study of Entrepreneurship Survival Dataset

abstract: In this article, we investigate the features which enhanced discriminate the survival in the micro and small business (MSE) using the approach of data mining with feature selection. According to the complexity of the data set, we proposed a comparison of three data imputation methods such as mean imputation (MI), k-nearest neighbor (KNN) and expectation maximization (EM) using mutually the selection of variables technique, whereby t-test, then through the data mining process using logistic regression classification methods, naive Bayes algorithm, linear discriminant analysis and support vector machine hence comparing their respective performances. The experimental results will be spread in developing a model to predict the MSE survival, providing a better understanding in the topic once it is a significant part of the Brazilian' GPA and macroeconomy.


title: Analysis of evolutionary origins of genomic loci harboring 59,732 candidate human-specific regulatory sequences identifies genetic divergence patterns during evolution of Great Apes

abstract: Our view of the universe of genomic regions harboring various types of candidate human-specific regulatory sequences (HSRS) has been markedly expanded in recent years. To infer the evolutionary origins of loci harboring HSRS, analyses of conservations patterns of 59,732 loci in Modern Humans, Chimpanzee, Bonobo, Gorilla, Orangutan, Gibbon, and Rhesus genomes have been performed. Two major evolutionary pathways have been identified comprising thousands of sequences that were either inherited from extinct common ancestors (ECAs) or created de novo in humans after human/chimpanzee split. Thousands of HSRS appear inherited from ECAs yet bypassed genomes of our closest evolutionary relatives, presumably due to the incomplete lineage sorting and/or species-specific loss or regulatory DNA. The bypassing pattern is prominent for HSRS associated with development and functions of human brain. Common genomic loci that may contributed to speciation during evolution of Great Apes comprise 248 insertions sites of African Great Ape-specific retrovirus PtERV1 (45.9%; p = 1.03E-44) intersecting regions harboring 442 HSRS, which are enriched for HSRS associated with human-specific (HS) changes of gene expression in cerebral organoids. Among non-human primates (NHP), most significant fractions of candidate HSRS associated with HS expression changes in both excitatory neurons (347 loci; 67%) and radial glia (683 loci; 72%) are highly conserved in Gorilla genome. Modern Humans acquired unique combinations of regulatory sequences highly conserved in distinct species of six NHP separated by 30 million years of evolution. Concurrently, this unique mosaic of regulatory sequences inherited from ECAs was supplemented with 12,486 created de novo HSRS. These observations support the model of complex continuous speciation process during evolution of Great Apes that is not likely to occur as an instantaneous event.


title: Collective motion of an atom array under laser illumination

abstract: We develop a theoretical formalism for the study of light-induced motion of atoms trapped in a two-dimensional (2D) array, considering the effect of multiple scattering of light between the atoms. We find that the atomic motion can be described by a collective diffusion equation, wherein laser-induced dipole-dipole forces couple the motion of different atoms. This coupling leads to the formation of collective mechanical modes of the array atoms, whose spatial structure and stability depend on the parameters of the illuminating laser and the geometry of the 2D array. We demonstrate the application of our formalism for the analysis of light-induced heating of the 2D array. The presented approach should be useful for treating the optomechanical properties of recently proposed quantum optical platforms made of atomic arrays.


title: Improving Sentence Representations with Consensus Maximisation

abstract: Consensus maximisation learning can provide self-supervision when different views are available of the same data. The distributional hypothesis provides another form of useful self-supervision from adjacent sentences which are plentiful in large unlabelled corpora. Motivated by the observation that different learning architectures tend to emphasise different aspects of sentence meaning, we present a new self-supervised learning framework for learning sentence representations which minimises the disagreement between two views of the same sentence where one view encodes the sentence with a recurrent neural network (RNN), and the other view encodes the same sentence with a simple linear model. After learning, the individual views (networks) result in higher quality sentence representations than their single-view learnt counterparts (learnt using only the distributional hypothesis) as judged by performance on standard downstream tasks. An ensemble of both views provides even better generalisation on both supervised and unsupervised downstream tasks. Also, importantly the ensemble of views trained with consensus maximisation between the two different architectures performs better on downstream tasks than an analogous ensemble made from the single-view trained counterparts.


title: Lattice points in vector-dilated quadratic irrational polytopes

abstract: We study the Ehrhart theory of quadratic irrational polytopes that undergo vector dilations. That is, for a given polytope with vertices in $\mathbb{Q}(\sqrt{D})$, and a different dilation factor for each facet, we show that the leading term of the lattice-point count behaves similar to an Ehrhart polynomial, generalizing previous work of Borda on scalar dilations of quadratic irrational polytopes. As a result, a form of the Ehrhart-Macdonald reciprocity law is obtained for the leading term.


title: PDE Acceleration: A convergence rate analysis and applications to obstacle problems

abstract: This paper provides a rigorous convergence rate and complexity analysis for a recently introduced framework, called PDE acceleration, for solving problems in the calculus of variations, and explores applications to obstacle problems. PDE acceleration grew out of a variational interpretation of momentum methods, such as Nesterov's accelerated gradient method and Polyak's heavy ball method, that views acceleration methods as equations of motion for a generalized Lagrangian action. Its application to convex variational problems yields equations of motion in the form of a damped nonlinear wave equation rather than nonlinear diffusion arising from gradient descent. These accelerated PDE's can be efficiently solved with simple explicit finite difference schemes where acceleration is realized by an improvement in the CFL condition from $dt\sim dx^2$ for diffusion equations to $dt\sim dx$ for wave equations. In this paper, we prove a linear convergence rate for PDE acceleration for strongly convex problems, provide a complexity analysis of the discrete scheme, and show how to optimally select the damping parameter for linear problems. We then apply PDE acceleration to solve minimal surface obstacle problems, including double obstacles with forcing, and stochastic homogenization problems with obstacles, obtaining state of the art computational results.


title: Study of the polarization produced by the Zeeman effect in the solar Mg I b lines

abstract: The next generation of solar observatories aim to understand the magnetism of the solar chromosphere. Therefore, it is crucial to understand the polarimetric signatures of chromospheric spectral lines. For this purpose, we here examine the suitability of the three Fraunhofer Mg I b1, b2, and b4 lines at 5183.6, 5172.7, and 5167.3 A, respectively. We start by describing a simplified atomic model of only 6 levels and 3 line transitions for computing the atomic populations of the 3p-4s (multiplet number 2) levels involved in the Mg I b line transitions assuming non-local thermodynamic conditions and considering only the Zeeman effect using the field-free approximation. We test this simplified atom against more complex ones finding that, although there are differences in the computed profiles, they are small compared with the advantages provided by the simple atom in terms of speed and robustness. After comparing the three Mg I lines, we conclude that the most capable one is the b2 line as b1 forms at similar heights and always show weaker polarization signals while b4 is severely blended with photospheric lines. We also compare Mg I b2 with the K I D1 and Ca II 8542 A lines finding that the former is sensitive to the atmospheric parameters at heights that are in between those covered by the latter two lines. This makes Mg I b2 an excellent candidate for future multi-line observations that aim to seamlessly infer the thermal and magnetic properties of different features in the lower solar atmosphere.


title: Soviet pioneers of fractional calculus and its application: III. Timothy D. Shermergor

abstract: We give a brief outline of biography, review 3 of his works on the use of confluent hypergeometric function for description asymmetric relaxation spectrum and provide an overview of the Soviet works of the past (20th) that are in tune with this topic.