title: Repulsive vacuum-induced forces on a magnetic particle

abstract: We study the possibility of obtaining a repulsive vacuum-induced force for a magnetic point particle near a surface. Considering the toy model of a particle with an electric-dipole transition and a large magnetic spin, we analyze the interplay between the repulsive magnetic-dipole and the attractive electric-dipole contributions to the total Casimir-Polder force. Particularly noting that the magnetic-dipole interaction is longer-ranged than the electric-dipole due to the difference in their respective characteristic transition frequencies, we find a regime where the repulsive magnetic contribution to the total force can potentially exceed the attractive electric part in magnitude for a sufficiently large spin. We analyze ways to further enhance the magnitude of the repulsive magnetic Casimir-Polder force for an excited particle, such as by preparing it in a "super-radiant" magnetic sub-level, and designing surface resonances close to the magnetic transition frequency.


title: Case Studies of Exocomets in the System of HD 10180

abstract: The aim of our study is to investigate the dynamics of possible comets in the HD 10180 system. This investigation is motivated by the discovery of exocomets in various systems, especially $\beta$ Pictoris, as well as in at least ten other systems. Detailed theoretical studies about the formation and evolution of star--planet systems indicate that exocomets should be quite common. Further observational results are expected in the foreseeable future, in part due to the availability of the Large Synoptic Survey Telescope. Nonetheless, the Solar System represents the best studied example for comets, thus serving as a prime motivation for investigating comets in HD 10180 as well. HD 10180 is strikingly similar to the Sun. This system contains six confirmed planets and (at least) two additional planets subject to final verification. In our studies, we consider comets of different inclinations and eccentricities and find an array of different outcomes such as encounters with planets, captures, and escapes. Comets with relatively large eccentricities are able to enter the inner region of the system facing early planetary encounters. Stable comets experience long-term evolution of orbital elements, as expected. We also tried to distinguish cometary families akin to our Solar System but no clear distinction between possible families was found. Generally, theoretical and observational studies of exoplanets have a large range of ramifications, involving the origin, structure and evolution of systems as well as the proliferation of water and prebiotic compounds to terrestrial planets, which will increase their chances of being habitable.


title: Invariant characterization of scalar third-order ODEs that admit the maximal point symmetry Lie algebra

abstract: The Cartan equivalence method is utilized to deduce an invariant characterization of the scalar third-order ordinary differential equation $u'"=f(x,u,u',u")$ which admits the maximal seven-dimensional point symmetry Lie algebra. The method provides auxiliary functions which can be used to efficiently obtain the point transformation that does the reduction to the simplest linear equation $\bar{u}'"=0$. Moreover, examples are given to illustrate the method.


title: Connected power domination in graphs

abstract: The study of power domination in graphs arises from the problem of placing a minimum number of measurement devices in an electrical network while monitoring the entire network. A power dominating set of a graph is a set of vertices from which every vertex in the graph can be observed, following a set of rules for power system monitoring. In this paper, we study the problem of finding a minimum power dominating set which is connected; the cardinality of such a set is called the connected power domination number of the graph. We show that the connected power domination number of a graph is NP-hard to compute in general, but can be computed in linear time in cactus graphs and block graphs. We also give various structural results about connected power domination, including a cut vertex decomposition and a characterization of the effects of various vertex and edge operations on the connected power domination number. Finally, we present novel integer programming formulations for power domination, connected power domination, and power propagation time, and give computational results.


title: The role of nonmagnetic d0 vs. d10 B-type cations on the magnetic exchange interactions in osmium double perovskites

abstract: Polycrystalline samples of double perovskites Ba2BOsO6 (B = Sc, Y, In) were synthesized by solid state reactions. They adopt the cubic double perovskite structures (space group, Fm-3m) with ordered B and Os arrangements. Ba2BOsO6 (B = Sc, Y, In) show antiferromagnetic transitions at 93 K, 69 K, and 28 K, respectively. The Weiss-temperatures are -590 K for Ba2ScOsO6, -571 K for Ba2YOsO6, and -155 K for Ba2InOsO6. Sc3+ and Y3+ have the open-shell d0 electronic configuration, while In3+ has the closed-shell d10. This indicates that a d0 B-type cation induces stronger overall magnetic exchange interactions in comparison to a d10. Comparison of Ba2BOsO6 (B = Sc, Y, In) to their Sr and Ca analogues shows that the structural distortions weaken the overall magnetic exchange interactions.


title: Noisy Natural Gradient as Variational Inference

abstract: Variational Bayesian neural nets combine the flexibility of deep learning with Bayesian uncertainty estimation. Unfortunately, there is a tradeoff between cheap but simple variational families (e.g.~fully factorized) or expensive and complicated inference procedures. We show that natural gradient ascent with adaptive weight noise implicitly fits a variational posterior to maximize the evidence lower bound (ELBO). This insight allows us to train full-covariance, fully factorized, or matrix-variate Gaussian variational posteriors using noisy versions of natural gradient, Adam, and K-FAC, respectively, making it possible to scale up to modern-size ConvNets. On standard regression benchmarks, our noisy K-FAC algorithm makes better predictions and matches Hamiltonian Monte Carlo's predictive variances better than existing methods. Its improved uncertainty estimates lead to more efficient exploration in active learning, and intrinsic motivation for reinforcement learning.


title: Utilizing hyperplanar substructures to perform efficient range-based WSN localization

abstract: A wireless sensor network (WSN) consists of multiple wireless sensor nodes that communicate each other to fulfill a particular task. In this paper, we emphasize on the networks whose deployments admit lower dimensional substructures, such as collinear groups in 2D, or coplanar groups in 3D. When these groups are given as a part of the input, we describe an algorithm to utilize this information to perform a low-cost localization. In emergency situations such as fire, earthquake etc. inside a building, wireless sensor networks might be very crucial to provide critical information and help the rescue teams to move very quickly by decreasing their burden of exploring the environment. Thus, it is very important to develop a system that provides information quickly and without consuming too much energy. We observe that in these type of environments, sensor nodes tend to form hyperplanar groups. A hyperplane is a subspace of one dimension less than its ambient space, and accordingly, a hyperplanar group of sensor nodes is a group of nodes that sit on the same hyperplane. When we consider a floor of a building, the nodes can be deployed on the corridors to form collinear groups, and when we consider the whole building, the sensor nodes will sit on the floors to form coplanar groups. Therefore, we study range-based WSN localization problem in 3D environments that induce hyperplanar groupings. First, we show that it is an \NP-hard problem to obtain hyperplanar groups, even if we are given the equations of the hyperplanes. Then, we describe an algorithm that assumes each sensor node is aware of its hyperplanar group, and performs a localization when the sensor nodes are deployed in the corridors on each floor of a building. In this case, grouping information allows us to localize a network, which cannot be localized by conventional range-based localization algorithms.


title: New Pulse Timing Measurements of the sdBV Star CS 1246

abstract: CS 1246 is a hot subdwarf B star discovered in 2009 to exhibit a single, large-amplitude radial pulsation. An O-C diagram constructed from this mode revealed reflex motion due to the presence of a low-mass M dwarf, as well as a long-term trend consistent with a decrease in the pulsational period. The orbital reflex motion was later confirmed with radial velocity measurements. Using eight years of data collected with the Skynet Robotic Telescope Network, we show that the pulsation amplitude of CS 1246 is decaying nonlinearly. We also present an updated O-C diagram, which might now indicate a positive $\dot P$ and a new $2.09 \pm 0.05$ yr oscillation consistent with orbital reflex motion of the entire inner sdB+dM binary, possibly due to the gravitational influence of a circumbinary planet with minimum mass $m\sin i = 3.3 \, \pm \, ^{2.1}_{1.2}$ $M_{\rm Jup}$. However, unlike the presence of the M dwarf, we hesitate to claim this object as a definitive detection since intrinsic variability of the pulsation phase could theoretically produce a similar effect.


title: Proposal for nanoscale cascaded plasmonic majority gates for non-Boolean computation

abstract: Surface-plasmon-polariton waves propagating at the interface between a metal and a dielectric, hold the key to future high-bandwidth, dense on-chip integrated logic circuits overcoming the diffraction limitation of photonics. While recent advances in plasmonic logic have witnessed the demonstration of basic and universal logic gates, these CMOS oriented digital logic gates cannot fully utilize the expressive power of this novel technology. Here, we aim at unraveling the true potential of plasmonics by exploiting an enhanced native functionality - the majority voter. Contrary to the state-of-the-art plasmonic logic devices, we use the phase of the wave instead of the intensity as the state or computational variable. We propose and demonstrate, via numerical simulations, a comprehensive scheme for building a nanoscale cascadable plasmonic majority logic gate along with a novel referencing scheme that can directly translate the information encoded in the amplitude and phase of the wave into electric field intensity at the output. Our MIM-based 3-input majority gate displays a highly improved overall area of only 0.636 {\mu}m$^2$ for a single-stage compared with previous works on plasmonic logic. The proposed device demonstrates non-Boolean computational capability and can find direct utility in highly parallel real-time signal processing applications like pattern recognition.


title: Generalized detailed balance relation for black hole horizon's fluctuations

abstract: In classical thermodynamics, irreversible processes are accomplished with an increase of entropy and a release of heat into the environment. In the case of black hole thermodynamics, instead, the increase of entropy is related with the absorption of energy from the outside of the black hole. Based on this analogy, the black hole is considered as a thermodynamic open system in the non-equilibrium regime during the absorption. Applying a novel result, we obtain a detailed balance relation for the macroscopic absorption of an astronomical object by an static spherically symmetric black hole. Because of the special characteristic of the thermodynamic system, it is necessary to include an additional term for the entropy produced by the absorption, allowing us to set a lower bound for the entropy production due to internal processes associated to the fluctuation.


title: Bounds on poloidal kinetic energy in plane layer convection

abstract: A numerical method is presented which conveniently computes upper bounds on heat transport and poloidal energy in plane layer convection for infinite and finite Prandtl numbers. The bounds obtained for the heat transport coincide with earlier results. These bounds imply upper bounds for the poloidal energy which follow directly from the definitions of dissipation and energy. The same constraints used for computing upper bounds on the heat transport lead to improved bounds for the poloidal energy.


title: Conflict-driven Hybrid Observer-based Anomaly Detection

abstract: This paper presents an anomaly detection method using a hybrid observer -- which consists of a discrete state observer and a continuous state observer. We focus our attention on anomalies caused by intelligent attacks, which may bypass existing anomaly detection methods because neither the event sequence nor the observed residuals appear to be anomalous. Based on the relation between the continuous and discrete variables, we define three conflict types and give the conditions under which the detection of the anomalies is guaranteed. We call this method conflict-driven anomaly detection. The effectiveness of this method is demonstrated mathematically and illustrated on a Train-Gate (TG) system.


title: Turbulence Generation from a stochastic wavelet model

abstract: This research presents a new turbulence generation method based on stochastic wavelets and tests its various properties in both homogeneous and inhomogeneous turbulence. Turbulence field can be generated with less basis compared to previous synthetic Fourier methods. Adaptive generation of inhomogeneous turbulence is achieved by scale reduction algorithm and lead to smaller computation cost. The generated turbulence shows good agreement with input data and theoretical results.


title: Software metadata: How much is enough?

abstract: Broad efforts are underway to capture metadata about research software and retain it across services; notable in this regard is the CodeMeta project. What metadata are important to have about (research) software? What metadata are useful for searching for codes? What would you like to learn about astronomy software? This BoF sought to gather information on metadata most desired by researchers and users of astro software and others interested in registering, indexing, capturing, and doing research on this software. Information from this BoF could conceivably result in changes to the Astrophysics Source Code Library (ASCL) or other resources for the benefit of the community or provide input into other projects concerned with software metadata.


title: Search for the flavor-changing neutral current interactions of the top quark and the Higgs boson which decays into a pair of b quarks at $\sqrt{s}=$ 13 TeV

abstract: A search for flavor-changing neutral currents (FCNC) in events with the top quark and the Higgs boson is presented. The Higgs boson decay to a pair of b quarks is considered. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ recorded by the CMS experiment at the LHC in proton-proton collisions at $\sqrt{s}=$ 13 TeV. Two channels are considered: single top quark FCNC production in association with the Higgs boson (pp $\to$ tH), and top quark pair production with FCNC decay of the top quark (t $\to$ qH). Final states with one isolated lepton and at least three reconstructed jets, among which at least two are associated with b quarks, are studied. No significant deviation is observed from the predicted background. Observed (expected) upper limits at 95% confidence level are set on the branching fractions of top quark decays, $\mathcal{B}$(t $\to$ uH) $<$ 0.47% (0.34%) and $\mathcal{B}$(t $\to$ cH) $<$ 0.47% (0.44%), assuming a single nonzero FCNC coupling.


title: Controllable Top-down Feature Transformer

abstract: We study the intrinsic transformation of feature maps across convolutional network layers with explicit top-down control. To this end, we develop top-down feature transformer (TFT), under controllable parameters, that are able to account for the hidden layer transformation while maintaining the overall consistency across layers. The learned generators capture the underlying feature transformation processes that are independent of particular training images. Our proposed TFT framework brings insights to and helps the understanding of, an important problem of studying the CNN internal feature representation and transformation under the top-down processes. In the case of spatial transformations, we demonstrate the significant advantage of TFT over existing data-driven approaches in building data-independent transformations. We also show that it can be adopted in other applications such as data augmentation and image style transfer.


title: FRW and domain walls in higher spin gravity

abstract: We present exact solutions to Vasiliev's bosonic higher spin gravity equations in four dimensions with positive and negative cosmological constant that admit an interpretation in terms of domain walls, quasi-instantons and Friedman-Robertson-Walker (FRW) backgrounds. Their isometry algebras are infinite dimensional higher-spin extensions of spacetime isometries generated by six Killing vectors. The solutions presented are obtained by using a method of holomorphic factorization in noncommutative twistor space and gauge functions. In interpreting the solutions in terms of Fronsdal-type fields in spacetime, a field-dependent higher spin transformation is required, which is implemented at leading order. To this order, the scalar field solves Klein-Gordon equation with conformal mass in (anti) de Sitter space. We interpret the FRW solution with de Sitter asymptotics in the context of inflationary cosmology and we expect that the domain wall and FRW solutions are associated with spontaneously broken scaling symmetries in their holographic description. We observe that the factorization method provides a convenient framework for setting up a perturbation theory around the exact solutions, and we propose that the nonlinear completion of particle excitations over FRW and domain wall solutions requires black hole-like states.


title: Differential Flatness of Quadrotor Dynamics Subject to Rotor Drag for Accurate Tracking of High-Speed Trajectories

abstract: In this paper, we prove that the dynamical model of a quadrotor subject to linear rotor drag effects is differentially flat in its position and heading. We use this property to compute feed-forward control terms directly from a reference trajectory to be tracked. The obtained feed-forward terms are then used in a cascaded, nonlinear feedback control law that enables accurate agile flight with quadrotors. Compared to state-of-the-art control methods, which treat the rotor drag as an unknown disturbance, our method reduces the trajectory tracking error significantly. Finally, we present a method based on a gradient-free optimization to identify the rotor drag coefficients, which are required to compute the feed-forward control terms. The new theoretical results are thoroughly validated trough extensive comparative experiments.


title: The oriented size Ramsey number of directed paths

abstract: An oriented graph is a directed graph with no bi-directed edges, i.e. if $xy$ is an edge then $yx$ is not an edge. The oriented size Ramsey number of an oriented graph $H$, denoted by $r(H)$, is the minimum $m$ for which there exists an oriented graph $G$ with $m$ edges, such that every $2$-colouring of $G$ contains a monochromatic copy of $H$. In this paper we prove that the oriented size Ramsey number of the directed paths on $n$ vertices satisfies $r(P_n) = \Omega(n^2 \log n)$. This improves a lower bound by Ben-Eliezer, Krivelevich and Sudakov. It also matches an upper bound by Buci\'c and the authors, thus establishing an asymptotically tight bound on $r(P_n)$. We also discuss how our methods can be used to improve the best known lower bound of the $k$-colour version of $r(P_n)$.


title: Effect of FePd alloy composition on the dynamics of artificial spin ice

abstract: Artificial spin ices (ASI) are arrays of single domain nano-magnetic islands, arranged in geometries that give rise to frustrated magnetostatic interactions. It is possible to reach their ground state via thermal annealing. We have made square ASI using different FePd alloys to vary the magnetization via co-sputtering. From a polarized state the samples were incrementally heated and we measured the vertex population as a function of temperature using magnetic force microscopy. For the higher magnetization FePd sample, we report an onset of dynamics at $T = 493$ K, with a rapid collapse into $>90\%$ ground state vertices. In contrast, the low magnetization sample started to fluctuate at lower temperatures, $T = 393$ K and over a wider temperature range but only reached a maximum of $25\%$ of ground state vertices. These results indicate that the interaction strength, dynamic temperature range and pathways can be finely tuned using a simple co-sputtering process. In addition we have compared our experimental values of the blocking temperature to those predicted using the simple N\'{e}el-Brown two-state model and find a large discrepancy which we attribute to activation volumes much smaller than the island volume.


title: Photospheric diagnostics of core helium burning in giant stars

abstract: Core helium burning primary red clump (RC) stars are evolved red giant stars which are excellent standard candles. As such, these stars are routinely used to map the Milky Way or determine the distance to other galaxies among other things. However distinguishing RC stars from their less evolved precursors, namely red giant branch (RGB) stars, is still a difficult challenge and has been deemed the domain of asteroseismology. In this work, we use a sample of 1,676 RGB and RC stars which have both single epoch infrared spectra from the APOGEE survey and asteroseismic parameters and classification to show that the spectra alone can be used to (1) predict asteroseismic parameters with precision high enough to (2) distinguish core helium burning RC from other giant stars with less than 2% contamination. This will not only allow for a clean selection of a large number of standard candles across our own and other galaxies from spectroscopic surveys, but also will remove one of the primary roadblocks for stellar evolution studies of mixing and mass loss in red giant stars.


title: MeV Pulsars: Modeling Spectra and Polarization

abstract: A sub-population of energetic rotation-powered pulsars show high fluxes of pulsed non-thermal hard X-ray emission. While this MeV pulsar population includes some radio-loud pulsars like the Crab, a significant number have no detected radio or GeV emission, a mystery since gamma- ray emission is a common characteristic of pulsars with high spin-down power. Their steeply rising hard X-ray spectral energy distributions (SEDs) suggest peaks at 0.1 - 1 MeV but they have not been detected above 200 keV. Several upcoming and planned telescopes may shed light on the MeV pulsars. The Neutron star Interior Composition ExploreR (NICER) will observe pulsars in the 0.2 - 12 keV band and may discover additional MeV pulsars. Planned telescopes, such as All-Sky Medium-Energy Gamma-Ray Observatory (AMEGO) and e-ASTROGAM, will detect emission above 0.2 MeV and polarization in the 0.2 - 10 MeV band. We present a model for the spectrum and polarization of MeV pulsars where the X-ray emission comes from electron- positron pairs radiating in the outer magnetosphere and current sheet. This model predicts that the peak of the SED increases with surface magnetic field strength if the pairs are produced in polar cap cascades. For small inclination angles, a range of viewing angles can miss both the radio pulse and the GeV pulse from particles accelerating near the current sheet. Characterizing the emission and geometry of MeV pulsars can thus provide clues to the source of pairs and acceleration in the magnetosphere.


title: Color superconductivity and charge neutrality in Yukawa theory

abstract: It is generally believed that when Cooper pairing occurs between two different species of fermions, their Fermi surfaces become locked together so that the resultant state remains 'neutral', with equal number densities of the two species, even when subjected to a chemical potential that couples to the difference in number densities. This belief is based on mean-field calculations in models with a zero-range interaction, where the anomalous self-energy is independent of energy and momentum. Following up on an early report of a deviation from neutrality in a Dyson-Schwinger calculation of color-flavor-locked quark matter, we investigate the neutrality of a two-species condensate using a Yukawa model which has a finite-range interaction. In a mean field calculation we obtain the full energy-momentum dependence of the self energy and find that the energy dependence leads to a population imbalance in the Cooper-paired phase when it is stressed by a species-dependent chemical potential. This gives some support to the suggestion that the color-flavor-locked phase of quark matter might not be an insulator.


title: Deep Regionlets for Object Detection

abstract: In this paper, we propose a novel object detection framework named "Deep Regionlets" by establishing a bridge between deep neural networks and conventional detection schema for accurate generic object detection. Motivated by the abilities of regionlets for modeling object deformation and multiple aspect ratios, we incorporate regionlets into an end-to-end trainable deep learning framework. The deep regionlets framework consists of a region selection network and a deep regionlet learning module. Specifically, given a detection bounding box proposal, the region selection network provides guidance on where to select regions to learn the features from. The regionlet learning module focuses on local feature selection and transformation to alleviate local variations. To this end, we first realize non-rectangular region selection within the detection framework to accommodate variations in object appearance. Moreover, we design a "gating network" within the regionlet leaning module to enable soft regionlet selection and pooling. The Deep Regionlets framework is trained end-to-end without additional efforts. We perform ablation studies and conduct extensive experiments on the PASCAL VOC and Microsoft COCO datasets. The proposed framework outperforms state-of-the-art algorithms, such as RetinaNet and Mask R-CNN, even without additional segmentation labels.


title: Pulsar Emission Physics: The First Fifty Years

abstract: Over the last fifty years since the discovery of pulsars, our understanding of where and how pulsars emit the radiation we observe has undergone significant revision. The location and mechanisms of high-energy radiation are intimately tied to the sites of particle acceleration. The evolution of emission models has paralleled the development of increasingly more sensitive telescopes, especially at high energies. I will review the history of pulsar emission modeling, from the early days of gaps at the polar caps, to outer gaps and slot gaps in the outer magnetosphere, to the present era of global magnetosphere simulations that locate most acceleration and high- energy emission in the current sheets.


title: Statistics students' identification of inferential model elements within contexts of their own invention

abstract: Statistical thinking partially depends upon an iterative process by which essential features of a problem setting are identified and mapped onto an abstract model or archetype, and then translated back into the context of the original problem setting (Wild and Pfannkuch 1999). Assessment in introductory statistics often relies on tasks that present students with data in context and expects them to choose and describe an appropriate model. This study explores post-secondary student responses to an alternative task that prompts students to clearly identify a sample, population, statistic, and parameter using a context of their own invention. The data include free text narrative responses of a random sample of 500 students from a sample of more than 1600 introductory statistics students. Results suggest that students' responses often portrayed sample and population accurately. Portrayals of statistic and parameter were less reliable and were associated with descriptions of a wide variety of other concepts. Responses frequently attributed a variable of some kind to the statistic, or a study design detail to the parameter. Implications for instruction and research are discussed, including a call for emphasis on a modeling paradigm in introductory statistics.


title: The BAHAMAS project: the CMB--large-scale structure tension and the roles of massive neutrinos and galaxy formation

abstract: Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) if one accounts for internal tensions in the Planck CMB dataset (by allowing the lensing amplitude, A_Lens, to vary), invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV, can resolve the tension. This solution is fully consistent with separate constraints from the primary CMB and baryon acoustic oscillations.


title: Estimating the error variance in a high-dimensional linear model

abstract: The lasso has been studied extensively as a tool for estimating the coefficient vector in the high-dimensional linear model; however, considerably less is known about estimating the error variance in this context. In this paper, we propose the natural lasso estimator for the error variance, which maximizes a penalized likelihood objective. A key aspect of the natural lasso is that the likelihood is expressed in terms of the natural parameterization of the multiparameter exponential family of a Gaussian with unknown mean and variance. The result is a remarkably simple estimator of the error variance with provably good performance in terms of mean squared error. These theoretical results do not require placing any assumptions on the design matrix or the true regression coefficients. We also propose a companion estimator, called the organic lasso, which theoretically does not require tuning of the regularization parameter. Both estimators do well empirically compared to preexisting methods, especially in settings where successful recovery of the true support of the coefficient vector is hard. Finally, we show that existing methods can do well under fewer assumptions than previously known, thus providing a fuller story about the problem of estimating the error variance in high-dimensional linear models.


title: The flux homomorphism on closed hyperbolic surfaces and Anti-de Sitter three-dimensional geometry

abstract: Given a smooth spacelike surface $\Sigma$ of negative curvature in Anti-de Sitter space of dimension 3, invariant by a representation $\rho:\pi_1(S)\to\mathrm{PSL}_2\mathbb{R}\times\mathrm{PSL}_2\mathbb{R}$ where $S$ is a closed oriented surface of genus $\geq 2$, a canonical construction associates to $\Sigma$ a diffeomorphism $\phi_\Sigma$ of $S$. It turns out that $\phi_\Sigma$ is a symplectomorphism for the area forms of the two hyperbolic metrics $h$ and $h'$ on $S$ induced by the action of $\rho$ on $\mathbb{H}^2\times\mathbb{H}^2$. Using an algebraic construction related to the flux homomorphism, we give a new proof of the fact that $\phi_\Sigma$ is the composition of a Hamiltonian symplectomorphism of $(S,h)$ and the unique minimal Lagrangian diffeomorphism from $(S,h)$ to $(S,h')$.


title: Are optical and X-ray AGN mostly disjoint?

abstract: The relationship between the populations of optically and X-ray selected Active Galactic Nuclei (AGN) has been unclear due to divergent results from different studies. Arnold et al. (2009) claim that X-ray AGN are almost entirely disjoint from optical AGN, while the Swift-BAT 70-month hard X-ray survey reported that 553 of their 711 X-ray AGN are optical. In this work, we set out to understand this difference by cross-checking between these studies and examining their sampling and AGN-selection criteria. We also re-analyze the X-ray and optical AGN in 16 groups and clusters reported by Arnold et al. using our own optical spectrum fitting techniques. We find that 6 of the 8 X-ray AGN in the Arnold et al. sample are also optical AGN, contrary to Arnold et al.'s report that only 1 of the 8 X-ray AGN is also an optical AGN, thereby falsifies their conclusion that optical and X-ray AGN are nearly disjoint sets.


title: Avalanches and Criticality in Driven Magnetic Skyrmions

abstract: We show using numerical simulations that slowly driven skyrmions interacting with random pinning move via correlated jumps or avalanches. The avalanches exhibit power law distributions in their duration and size, and the average avalanche shape for different avalanche durations can be scaled to a universal function, in agreement with theoretical predictions for systems in a nonequilibrium critical state. A distinctive feature of skyrmions is the influence of the non-dissipative Magnus term. When we increase the ratio of the Magnus term to the damping term, a change in the universality class of the behavior occurs, the average avalanche shape becomes increasingly asymmetric, and individual avalanches exhibit motion in the direction perpendicular to their own density gradient.


title: Pieri rules for the Jack polynomials in superspace and the 6-vertex model

abstract: We present Pieri rules for the Jack polynomials in superspace. The coefficients in the Pieri rules are, except for an extra determinant, products of quotients of linear factors in $\alpha$ (expressed, as in the usual Jack polynomial case, in terms of certain hook-lengths in a Ferrers' diagram). We show that, surprisingly, the extra determinant is related to the partition function of the 6-vertex model. We give, as a conjecture, the Pieri rules for the Macdonald polynomials in superspace.


title: Identified kaon production in Ar+Sc collisions at SPS energies

abstract: NA61/SHINE is a fixed target experiment at the CERN Super Proton Synchrotron. The main goals of the experiment are to discover the critical point of strongly interacting matter and to study the properties of the onset of deconfinement. In order to reach these goals, a study of hadron production properties is performed in nucleus-nucleus, proton-proton and proton-nucleus interactions as a function of collision energy and size of the colliding nuclei. In this talk, the newest preliminary results on kaon spectra produced in Ar+Sc collisions at three beam momenta (30A, 40A and 75A) will be shown. The distributions of transverse mass and rapidity will be compared with results of NA61/SHINE (p+p, Be+Be) and NA49 (Pb+Pb, C+C, Si+Si), as well as with available world data.


title: Dynamic quantum kagome ice

abstract: The search for two dimensional quantum spin liquids, exotic magnetic states with an entangled ground state remaining disordered down to zero temperature, has been a great challenge in frustrated magnetism during the last decades. Recently, fractionalized excitations, called spinons, emerging from these states, have been evidenced in kagome and triangular antiferromagnets. In contrast, quantum ferromagnetic spin liquids in two dimensions, namely quantum kagome ices, have been less investigated, yet their classical counterparts exhibit amazing properties, magnetic monopole crystals as well as magnetic fragmentation. Here we show that, by applying a magnetic field on the pyrochlore oxide Nd$_2$Zr$_2$O$_7$, which has been shown to develop three dimensional quantum magnetic fragmentation in zero field, we are able to reduce the dimension of the system and to create a dynamic kagome ice state. Our results open the way to the observation of the quantum kagome ice state which was recently investigated theoretically.


title: Localization of eigenfunctions via an effective potential

abstract: We consider the localization of eigenfunctions for the operator $L=-\mbox{div} A \nabla + V$ on a Lipschitz domain $\Omega$ and, more generally, on manifolds with and without boundary. In earlier work, two authors of the present paper demonstrated the remarkable ability of the landscape, defined as the solution to $Lu=1$, to predict the location of the localized eigenfunctions. Here, we explain and justify a new framework that reveals a richly detailed portrait of the eigenfunctions and eigenvalues. We show that the reciprocal of the landscape function, $1/u$, acts as an effective potential. Hence from the single measurement of $u$, we obtain, via $1/u$, explicit bounds on the exponential decay of the eigenfunctions of the system and estimates on the distribution of eigenvalues near the bottom of the spectrum. (This version strengthens and simplifies the results of the first one by replacing a global bi-Lipschitz hypothesis on the domain with a local bi-Lipschitz hypothesis. It improves on the second version by adding pictures and numerical examples. This version is identical to the third version; all that is changed is to correct some tex mistakes in symbols in this abstract. There are no changes to the paper itself.)


title: Angular momentum transport by heat-driven g-modes in slowly pulsating B stars

abstract: Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution code. Focusing first on a $4.21$ $M_\odot$ model, we simulate 1,000 years of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation the initially non-rotating surface layers are spun in the retrograde direction to $\approx30\%$ of the critical rate. However, the additional inclusion of magnetic stresses in our simulations, almost completely inhibits this spin-up. Expanding our simulations to cover the whole instability strip, we show that the same general behavior is seen in all SPB stars. After providing some caveats to contextualize our results, we hypothesize that the observed slower surface rotation of SPB stars (as compared to other B-type stars) may be the direct consequence of the angular momentum transport that our simulations demonstrate.


title: The Free-play Sandbox: a Methodology for the Evaluation of Social Robotics and a Dataset of Social Interactions

abstract: Evaluating human-robot social interactions in a rigorous manner is notoriously difficult: studies are either conducted in labs with constrained protocols to allow for robust measurements and a degree of replicability, but at the cost of ecological validity; or in the wild, which leads to superior experimental realism, but often with limited replicability and at the expense of rigorous interaction metrics. We introduce a novel interaction paradigm, designed to elicit rich and varied social interactions while having desirable scientific properties (replicability, clear metrics, possibility of either autonomous or Wizard-of-Oz robot behaviours). This paradigm focuses on child-robot interactions, and builds on a sandboxed free-play environment. We present the rationale and design of the interaction paradigm, its methodological and technical aspects (including the open-source implementation of the software platform), as well as two large open datasets acquired with this paradigm, and meant to act as experimental baselines for future research.


title: High-End State-of-Art Copper Bromide Vapor Lasers: Optical, Electric and Thermal Properties

abstract: Copper bromide vapor laser is the most successful sealed-off version of a plenty of copper vapor lasers launched, studied and produced for the last half of century. Being first reported in 1975, it became the leading and long-lasting R&D&T project followed by the Laboratory of metal vapor lasers in Institute of solid state physics of the Bulgarian Academy of Sciences ever since. This manuscript gathers major research & development results of the last studies unfolding unexplored and unknown optical, electric and thermal properties of CuBr laser held in the Laboratory of metal vapor lasers during foregoing 20+ years. Some of these new aftermaths are record-high power, efficiency and lifetime of CuBr laser; comprehensive studies of energy dissipation and dynamics of the electric field of CuBr laser; achievement of high-brightness diffraction-limited throughout-pulse emission of master-oscillator power-amplifier (MOPA) CuBr laser and high focusability of CuBr laser emission low-affected by beam profile intensity variations.


title: Tomographic Reconstruction using Global Statistical Prior

abstract: Recent research in tomographic reconstruction is motivated by the need to efficiently recover detailed anatomy from limited measurements. One of the ways to compensate for the increasingly sparse sets of measurements is to exploit the information from templates, i.e., prior data available in the form of already reconstructed, structurally similar images. Towards this, previous work has exploited using a set of global and patch based dictionary priors. In this paper, we propose a global prior to improve both the speed and quality of tomographic reconstruction within a Compressive Sensing framework. We choose a set of potential representative 2D images referred to as templates, to build an eigenspace; this is subsequently used to guide the iterative reconstruction of a similar slice from sparse acquisition data. Our experiments across a diverse range of datasets show that reconstruction using an appropriate global prior, apart from being faster, gives a much lower reconstruction error when compared to the state of the art.


title: Evidence for Merger-Driven Growth in Luminous, High-z, Obscured AGN in the CANDELS/COSMOS Field

abstract: While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN.


title: A method of immediate detection of objects with a near-zero apparent motion in series of CCD-frames

abstract: The paper deals with a computational method for detection of the solar system minor bodies (SSOs), whose inter-frame shifts in series of CCD-frames during the observation are commensurate with the errors in measuring their positions. These objects have velocities of apparent motion between CCD-frames not exceeding three RMS errors ($3\sigma$) of measurements of their positions. About 15\% of objects have a near-zero apparent motion in CCD-frames, including the objects beyond the Jupiter's orbit as well as the asteroids heading straight to the Earth. The proposed method for detection of the object's near-zero apparent motion in series of CCD-frames is based on the Fisher f-criterion instead of using the traditional decision rules that are based on the maximum likelihood criterion. We analyzed the quality indicators of detection of the object's near-zero apparent motion applying statistical and in situ modeling techniques in terms of the conditional probability of the true detection of objects with a near-zero apparent motion. The efficiency of method being implemented as a plugin for the Collection Light Technology (CoLiTec) software for automated asteroids and comets detection has been demonstrated. Among the objects discovered with this plugin, there was the sungrazing comet C/2012 S1 (ISON). Within 26 minutes of the observation, the comet's image has been moved by three pixels in a series of four CCD-frames (the velocity of its apparent motion at the moment of discovery was equal to 0.8 pixels per CCD-frame; the image size on the frame was about five pixels). Next verification in observations of asteroids with a near-zero apparent motion conducted with small telescopes has confirmed an efficiency of the method even in bad conditions (strong backlight from the full Moon). So, we recommend applying the proposed method for series of observations with four or more frames.


title: Compressive Phase Retrieval via Reweighted Amplitude Flow

abstract: The problem of reconstructing a sparse signal vector from magnitude-only measurements (a.k.a., compressive phase retrieval), emerges naturally in diverse applications, but it is NP-hard in general. Building on recent advances in nonconvex optimization, this paper puts forth a new algorithm that is termed compressive reweighted amplitude flow and abbreviated as CRAF, for compressive phase retrieval. Specifically, CRAF operates in two stages. The first stage seeks a sparse initial guess via a new spectral procedure. In the second stage, CRAF implements a few hard thresholding based iterations using reweighted gradients. When there are sufficient measurements, CRAF provably recovers the underlying signal vector exactly with high probability under suitable conditions. Moreover, its sample complexity coincides with that of the state-of-the-art procedures. Finally, substantial simulated tests showcase remarkable performance of the new spectral initialization, as well as improved exact recovery relative to competing alternatives.


title: High performance ultra-low-precision convolutions on mobile devices

abstract: Many applications of mobile deep learning, especially real-time computer vision workloads, are constrained by computation power. This is particularly true for workloads running on older consumer phones, where a typical device might be powered by a single- or dual-core ARMv7 CPU. We provide an open-source implementation and a comprehensive analysis of (to our knowledge) the state of the art ultra-low-precision (<4 bit precision) implementation of the core primitives required for modern deep learning workloads on ARMv7 devices, and demonstrate speedups of 4x-20x over our additional state-of-the-art float32 and int8 baselines.


title: Link mutations and Goeritz matrices

abstract: Extending theorems of J. E. Greene [Invent. Math. 192 (2013), 717-750] and A. S. Lipson [Enseign. Math. (2) 36 (1990), 93-114], we prove that the equivalence class of a classical link L under mutation is determined by Goeritz matrices associated to diagrams of L.


title: Emergence of a spectral gap in a class of random matrices associated with split graphs

abstract: Motivated by the intriguing behavior displayed in a dynamic network that models a population of extreme introverts and extroverts (XIE), we consider the spectral properties of ensembles of random split graph adjacency matrices. We discover that, in general, a gap emerges in the bulk spectrum between -1 and 0 that contains a single eigenvalue. An analytic expression for the bulk distribution is derived and verified with numerical analysis. We also examine their relation to chiral ensembles, which are associated with bipartite graphs.

categories: physics.soc-ph cond-mat.stat-mech math.SP nlin.AO


title: Direct numerical simulations of turbulent channel flow under transcritical conditions

abstract: Turbulent flows under transcritical conditions are present in regenerative cooling systems of rocker engines and extraction processes in chemical engineering. The turbulent flows and the corresponding heat transfer phenomena in these complex processes are still not well understood experimentally and numerically. The objective of this work is to investigate the turbulent flows under transcritical conditions using DNS of turbulent channel flows. A fully compressible solver is used in conjunction with a Peng-Robinson real-fluid equation of state to describe the transcritical flows. A channel flow with two isothermal walls is simulated with one heated and one cooled boundary layers. The grid resolution adopted in this study is slightly finer than that required for DNS of incompressible channel flows. The simulations are conducted using both fully (FC) and quasi-conservative (QC) schemes to assess their performance for transcritical wall-bounded flows. The instantaneous flows and the statistics are analyzed and compared with the canonical theories. It is found that results from both FC and QC schemes qualitatively agree well with noticeable difference near the top heated wall, where spurious oscillations in velocity can be observed. Using the DNS data, we then examine the usefulness of Townsend attached eddy hypothesis in the context of flows at transcritical conditions. It is shown that the streamwise energy spectrum exhibits the inverse wavenumber scaling and that the streamwise velocity structure function follows a logarithmic scaling, thus providing support to the attached eddy model at transcritical conditions.


title: Numerical analysis on mixing processes for transcritical real-fluid simulations

abstract: The accurate and robust simulation of transcritical real-fluid flows is crucial for many engineering applications. Diffused interface methods are frequently employed and several numerical schemes have been developed for simulating transcritical flows. These schemes can be categorized into two types, namely fully conservative and quasi-conservative schemes. An adaptive scheme which is a hybrid of the two is developed in this study. By considering several numerical test cases, it is shown that different schemes predict distinctly different mixing behaviors. It is shown that the mixing processes follow the isobaric-adiabatic and isobaric-isochoric mixing models for fully and quasi-conservative schemes, respectively, and the adaptive scheme yields a mixing behavior that spans both models. The distinct mixing behaviors are a consequence of numerical diffusion instead of physical diffusion and can be attributed to insufficient numerical spatial resolution. This work provides a better understanding on the interpretation of numerical simulation results and the mixing models that are commonly used to study transcritical flows.


title: Sparse learning of stochastic dynamic equations

abstract: With the rapid increase of available data for complex systems, there is great interest in the extraction of physically relevant information from massive datasets. Recently, a framework called Sparse Identification of Nonlinear Dynamics (SINDy) has been introduced to identify the governing equations of dynamical systems from simulation data. In this study, we extend SINDy to stochastic dynamical systems, which are frequently used to model biophysical processes. We prove the asymptotic correctness of stochastics SINDy in the infinite data limit, both in the original and projected variables. We discuss algorithms to solve the sparse regression problem arising from the practical implementation of SINDy, and show that cross validation is an essential tool to determine the right level of sparsity. We demonstrate the proposed methodology on two test systems, namely, the diffusion in a one-dimensional potential, and the projected dynamics of a two-dimensional diffusion process.


title: Identifying Data Noises, User Biases, and System Errors in Geo-tagged Twitter Messages (Tweets)

abstract: Many social media researchers and data scientists collected geo-tagged tweets to conduct spatial analysis or identify spatiotemporal patterns of filtered messages for specific topics or events. This paper provides a systematic view to illustrate the characteristics (data noises, user biases, and system errors) of geo-tagged tweets from the Twitter Streaming API. First, we found that a small percentage (1%) of active Twitter users can create a large portion (16%) of geo-tagged tweets. Second, there is a significant amount (57.3%) of geo-tagged tweets located outside the Twitter Streaming API's bounding box in San Diego. Third, we can detect spam, bot, cyborg tweets (data noises) by examining the "source" metadata field. The portion of data noises in geo-tagged tweets is significant (29.42% in San Diego, CA and 53.47% in Columbus, OH) in our case study. Finally, the majority of geo-tagged tweets are not created by the generic Twitter apps in Android or iPhone devices, but by other platforms, such as Instagram and Foursquare. We recommend a multi-step procedure to remove these noises for the future research projects utilizing geo-tagged tweets.


title: Filamentary Structure of Orion A

abstract: We present filamentary structure analysis based on the CARMA Orion Project through which we obtained high resolution (0.01pc) and high dynamic range (2~deg in DEC) images of the Orion A molecular cloud. We show a preliminary emission map of the $^{13}$CO emission with the addition of identified filaments. The data product resulted from 1500~hr of observing time at CARMA interferometer, and the Nobeyama 45~m telescope. We address basic filament properties along with calculated typical kinematical properties. Furthermore, we present a pilot study where we compare the observed [C{\sc II}] and [C{\sc I}] emission to $^{13}$CO emission towards 6 selected filaments in Orion~A.


title: Expectation of primordial gravity waves generated during inflation

abstract: The inflationary paradigm is extremely successful regarding predictions of temperature anisotropies in the CMB. However, inflation also makes predictions for a CMB B-mode polarization, which has not been detected. Moreover, the standard inflationary paradigm is unable to accommodate the evolution from the initial state, which is assumed to be symmetric, into a non-symmetric aftermath. In Phys. Rev. D 96, 101301(R), we show that the incorporation of an element capable of explaining such a transition drastically changes the prediction for the shape and size of the B-mode spectrum. In particular, employing a realistic objective collapse model in a well-defined semiclassical context, we find that, while predictions of temperature anisotropies are nor altered (with respect to standard predictions), the B-mode spectrum gets strongly suppressed--in accordance with observations. Here we present an in-depth discussion of that analysis, together with the details of the calculation.


title: Performance Improvement by Introducing Mobility in Wireless Communication Networks

abstract: Communication technology is a major contributor to our lifestyles. Improving the performance of communication system brings with various benefits to human beings. This report covers two typical such systems: wireless sensor networks and cellular networks. We study the influence of mobility in these networks. In wireless sensor networks, we consider the usage of mobile sinks to collect sensory data from static nodes. For the latter case, we design data collection protocols for a single mobile sink which aim at balancing the energy consumption among sensor nodes to improve the network lifetime. We also design an algorithm for multiple mobile sinks to collect urgent sensory information within the allowed latency. Regarding the mobility in cellular networks, we mean the service providers are mobile. Conventional approaches usually consider how to optimally deploy static base stations according to a certain metric, such as throughput. Base station densification is a solution, but it is not cost efficient. Introducing mobility is to let the base stations fly in the sky, resulting in less investment. We consider one of the key issues of using drones to serve cellular users: drone deployment. We have implemented all the proposed algorithms by either simulations or experiments, and the results have confirmed the effectiveness of these approaches.


title: The Chevalley-Weil Formula for Orbifold Curves

abstract: In the 1930s Chevalley and Weil gave a formula for decomposing the canonical representation on the space of differential forms of the Galois group of a ramified Galois cover of Riemann surfaces. In this article we prove an analogous Chevalley-Weil formula for ramified Galois covers of orbifold curves. We then specialize the formula to the case when the base orbifold curve is the (reduced) modular orbifold. As an application of this latter formula we decompose the canonical representations of modular curves of full, prime level and of Fermat curves of arbitrary exponent.


title: Real-time Video Processing in Web Applications

abstract: The OpenGL ES standard is implemented in modern desktop and mobile browsers through the WebGL API. This paper explores the potential for using OpenGL ES hardware acceleration for real time video processing in standard HTML5 applications. It analyses the WebGL performance across device types and compares it with the standard JavaScript and canvas performance.


title: Studies on the anomalous magnetic and electric dipole moments of the tau-neutrino in $pp$ collisions at the LHC

abstract: In this paper the production cross section $pp\rightarrow (\gamma, Z) \to \nu_\tau \bar \nu_\tau \gamma+X$ in $pp$ collisions at $\sqrt{s}=8, 13, 14, 33TeV$ is presented. Furthermore, we estimate bounds at the $95\%C. L.$ on the dipole moments of the tau-neutrino using integrated luminosity of ${\cal L}=20, 50, 100, 200, 500, 1000, 3000 fb^{-1}$ collected with the ATLAS detector at the LHC and we consider systematic uncertainties of $\delta_{sys}=0, 5, 10\%$. It is shown that the process under consideration is a good prospect for probing the dipole moments of the tau-neutrino at the LHC.


title: Dynamical Equilibration of Topological Properties

abstract: We study the dynamical process of equilibration of topological properties in quantum many-body systems undergoing a parameter quench between two topologically inequivalent Hamiltonians. This scenario is motivated by recent experiments on ultracold atomic gases, where a trivial initial state is prepared before the Hamiltonian is ramped into a topological insulator phase. While the many-body wave function must stay topologically trivial in the coherent post-quench dynamics, here we show how the topological properties of the single particle density matrix dynamically change and equilibrate in the presence of interactions. In this process, the single particle density matrix goes through a characteristic level crossing as a function of time, which plays an analogous role to the gap closing of a Hamiltonian in an equilibrium topological quantum phase transition. As an exact case study exemplifying this mechanism, we numerically solve the quench dynamics of an interacting one-dimensional topological insulator.


title: A Novel Model for Arbitration between Planning and Habitual Control Systems

abstract: It is well established that humans decision making and instrumental control uses multiple systems, some which use habitual action selection and some which require deliberate planning. Deliberate planning systems use predictions of action-outcomes using an internal model of the agent's environment, while habitual action selection systems learn to automate by repeating previously rewarded actions. Habitual control is computationally efficient but may be inflexible in changing environments. Conversely, deliberate planning may be computationally expensive, but flexible in dynamic environments. This paper proposes a general architecture comprising both control paradigms by introducing an arbitrator that controls which subsystem is used at any time. This system is implemented for a target-reaching task with a simulated two-joint robotic arm that comprises a supervised internal model and deep reinforcement learning. Through permutation of target-reaching conditions, we demonstrate that the proposed is capable of rapidly learning kinematics of the system without a priori knowledge, and is robust to (A) changing environmental reward and kinematics, and (B) occluded vision. The arbitrator model is compared to exclusive deliberate planning with the internal model and exclusive habitual control instances of the model. The results show how such a model can harness the benefits of both systems, using fast decisions in reliable circumstances while optimizing performance in changing environments. In addition, the proposed model learns very fast. Finally, the system which includes internal models is able to reach the target under the visual occlusion, while the pure habitual system is unable to operate sufficiently under such conditions.


title: Exploration of the High Entropy Alloy Space as a Constraint Satisfaction Problem

abstract: High Entropy Alloys (HEAs), Multi-principal Component Alloys (MCA), or Compositionally Complex Alloys (CCAs) are alloys that contain multiple principal alloying elements. While many HEAs have been shown to have unique properties, their discovery has been largely done through costly and time-consuming trial-and-error approaches, with only an infinitesimally small fraction of the entire possible composition space having been explored. In this work, the exploration of the HEA composition space is framed as a Continuous Constraint Satisfaction Problem (CCSP) and solved using a novel Constraint Satisfaction Algorithm (CSA) for the rapid and robust exploration of alloy thermodynamic spaces. The algorithm is used to discover regions in the HEA Composition-Temperature space that satisfy desired phase constitution requirements. The algorithm is demonstrated against a new (TCHEA1) CALPHAD HEA thermodynamic database. The database is first validated by comparing phase stability predictions against experiments and then the CSA is deployed and tested against design tasks consisting of identifying not only single phase solid solution regions in ternary, quaternary and quinary composition spaces but also the identification of regions that are likely to yield precipitation-strengthened HEAs.


title: A Macromodeling Approach to Efficiently Compute Scattering from Large Arrays of Complex Scatterers

abstract: Full-wave electromagnetic simulations of electrically large arrays of complex antennas and scatterers are challenging, as they consume large amount of memory and require long CPU times. This paper presents a new reduced-order modeling technique to compute scattering and radiation from large arrays of complex scatterers and antennas. In the proposed technique, each element of the array is replaced by an equivalent electric current distribution on a fictitious closed surface enclosing the element. This equivalent electric current density is derived using the equivalence theorem and it is related to the surface currents on the scatterer by the Stratton-Chu formulation. With the proposed approach, instead of directly solving for the unknown surface current density on the scatterers, we only need to solve for the unknowns on the equivalent surface. This approach leads to a reduction in the number of unknowns and better conditioning when it is applied to problems involving complex scatterers with multiscale features. Furthermore, the proposed approach is accelerated with the adaptive integral equation method to solve large problems. As illustrated in several practical examples, the proposed method yields speed up of up to 20 times and consumes up to 12 times less memory than the standard method of moments accelerated with the adaptive integral method.


title: Effects of non-Hermitian perturbations on Weyl Hamiltonians with arbitrary topological charges

abstract: We provide a systematic study of non-Hermitian topologically charged systems. Starting from a Hermitian Hamiltonian supporting Weyl points with arbitrary topological charge, adding a non-Hermitian perturbation transforms the Weyl points to one-dimensional exceptional contours. We analytical prove that the topological charge is preserved on the exceptional contours. In contrast to Hermitian systems, the addition of gain and loss allows for a new class of topological phase transition: when two oppositely charged exceptional contours touch, the topological charge can dissipate without opening a gap. These effects can be demonstrated in realistic photonics and acoustics systems.


title: Targeted Random Projection for Prediction from High-Dimensional Features

abstract: We consider the problem of computationally-efficient prediction from high dimensional and highly correlated predictors in challenging settings where accurate variable selection is effectively impossible. Direct application of penalization or Bayesian methods implemented with Markov chain Monte Carlo can be computationally daunting and unstable. Hence, some type of dimensionality reduction prior to statistical analysis is in order. Common solutions include application of screening algorithms to reduce the regressors, or dimension reduction using projections of the design matrix. The former approach can be highly sensitive to threshold choice in finite samples, while the later can have poor performance in very high-dimensional settings. We propose a TArgeted Random Projection (TARP) approach that combines positive aspects of both strategies to boost performance. In particular, we propose to use information from independent screening to order the inclusion probabilities of the features in the projection matrix used for dimension reduction, leading to data-informed sparsity. We provide theoretical support for a Bayesian predictive algorithm based on TARP, including both statistical and computational complexity guarantees. Examples for simulated and real data applications illustrate gains relative to a variety of competitors.


title: HyperPower: Power- and Memory-Constrained Hyper-Parameter Optimization for Neural Networks

abstract: While selecting the hyper-parameters of Neural Networks (NNs) has been so far treated as an art, the emergence of more complex, deeper architectures poses increasingly more challenges to designers and Machine Learning (ML) practitioners, especially when power and memory constraints need to be considered. In this work, we propose HyperPower, a framework that enables efficient Bayesian optimization and random search in the context of power- and memory-constrained hyper-parameter optimization for NNs running on a given hardware platform. HyperPower is the first work (i) to show that power consumption can be used as a low-cost, a priori known constraint, and (ii) to propose predictive models for the power and memory of NNs executing on GPUs. Thanks to HyperPower, the number of function evaluations and the best test error achieved by a constraint-unaware method are reached up to 112.99x and 30.12x faster, respectively, while never considering invalid configurations. HyperPower significantly speeds up the hyper-parameter optimization, achieving up to 57.20x more function evaluations compared to constraint-unaware methods for a given time interval, effectively yielding significant accuracy improvements by up to 67.6%.


title: On Colouring $(2P_2,H)$-Free and $(P_5,H)$-Free Graphs

abstract: The Colouring problem asks whether the vertices of a graph can be coloured with at most $k$ colours for a given integer $k$ in such a way that no two adjacent vertices receive the same colour. A graph is $(H_1,H_2)$-free if it has no induced subgraph isomorphic to $H_1$ or $H_2$. A connected graph $H_1$ is almost classified if Colouring on $(H_1,H_2)$-free graphs is known to be polynomial-time solvable or NP-complete for all but finitely many connected graphs $H_2$. We show that every connected graph $H_1$ apart from the claw $K_{1,3}$ and the $5$-vertex path $P_5$ is almost classified. We also prove a number of new hardness results for Colouring on $(2P_2,H)$-free graphs. This enables us to list all graphs $H$ for which the complexity of Colouring is open on $(2P_2,H)$-free graphs and all graphs $H$ for which the complexity of Colouring is open on $(P_5,H)$-free graphs. In fact we show that these two lists coincide. Moreover, we show that the complexities of Colouring for $(2P_2,H)$-free graphs and for $(P_5,H)$-free graphs are the same for all known cases.


title: Towards Methods for Model-Based Software Development

abstract: Software engineering is a young discipline. Despite efforts in recent years, some elements still require further development, research, and systematization. One of these elements are methods. They consist of a set of well-defined activities used in a software development to guide how the work should be performed to achieve the expected results. However, there is a lack of systematic knowledge that effectively guide how work should be done in a variety of areas. Despite the knowledge available in Software Engineering Body of Knowledge (SWEBOK), only the classical methods are comprised there, while emerging topics such as Agent-Oriented Software Engineering (AOSE), Aspect-Oriented Software Development (AOSD), and Model-Based Software Development (MBSE) are only briefly mentioned. We claim that it is also necessary to investigate methods for those emerging software engineering subtracks. In this direction, this paper presents a a speculative and preliminary emerging results on the establishment of methods for MBSE. We report some insights on the conception of methods for MBSE.


title: Quantifying how much sensory information in a neural code is relevant for behavior

abstract: Determining how much of the sensory information carried by a neural code contributes to behavioral performance is key to understand sensory function and neural information flow. However, there are as yet no analytical tools to compute this information that lies at the intersection between sensory coding and behavioral readout. Here we develop a novel measure, termed the information-theoretic intersection information $I_{II}(S;R;C)$, that quantifies how much of the sensory information carried by a neural response R is used for behavior during perceptual discrimination tasks. Building on the Partial Information Decomposition framework, we define $I_{II}(S;R;C)$ as the part of the mutual information between the stimulus S and the response R that also informs the consequent behavioral choice C. We compute $I_{II}(S;R;C)$ in the analysis of two experimental cortical datasets, to show how this measure can be used to compare quantitatively the contributions of spike timing and spike rates to task performance, and to identify brain areas or neural populations that specifically transform sensory information into choice.


abstract: In this paper, we introduce the concept of Continuous $\ast$-g-Frame in Hilbert $C^{\ast}$-Modules and we establish some results. We also discuss the stability problem for Continuous $\ast$-g-Frame.


title: Chiral asymmetry during the EWPT from CP-violating scattering off bubble walls

abstract: We compute a net electric current during a first order EWPT arising from the asymmetric propagation of fermion chiral modes due to a CP-violating interaction with the Higgs. The interaction is quantified in terms of a CP-violating phase in the bubble wall that separate both false and true vacuum phases. We comment on the possibility of this current to generate a seed magnetic field and its implications for primordial magne- togenesis in the early Universe.


title: Opinion Dynamics and Social Power Evolution: A Single-Timescale Model

abstract: This paper studies the evolution of self-appraisal and social power, for a group of individuals who discuss and form opinions. We consider a modification of the recently proposed DeGroot-Friedkin (DF) model, in which the opinion formation process takes place on the same timescale as the reflected appraisal process; we call this new model the single-timescale DF model. We provide a comprehensive analysis of the equilibria and convergence properties of the model for the settings of irreducible and reducible influence networks. For the setting of irreducible influence networks, the single-timescale DF model has the same behavior as the original DF model, that is, it predicts among other things that the social power ranking among individuals is asymptotically equal to their centrality ranking, that social power tends to accumulate at the top of the centrality ranking hierarchy, and that an autocratic (resp., democratic) power structure arises when the centrality scores are maximally nonuniform (resp., uniform). For the setting of reducible influence networks, the single-timescale DF model behaves differently from the original DF model in two ways. First, an individual, who corresponds to a reducible node in a reducible influence network, can keep all social power in the single-timescale DF model if the initial condition does so, whereas its social power asymptotically vanishes in the original DF model. Second, when the associated network has multiple sinks, the two models behave very differently: the original DF model has a single globally-attractive equilibrium, whereas any partition of social power among the sinks is allowable at equilibrium in the single-timescale DF model.


title: Beam Alignment and Tracking for Autonomous Vehicular Communication using IEEE 802.11ad-based Radar

abstract: Mobility scenarios involving short contact times pose a challenge for high bandwidth data transfer between autonomous vehicles and roadside base stations (BS). Millimeter wave bands are a viable solution as they offer enormous bandwidth in the 60GHz band with several Gbps data transfer rates. However, beamforming is used as a default mode in this band, which requires accurate and continuous alignment under relative motion. We propose a method in which an off-the-shelf IEEE 802.11ad WiFi router is configured to serve as the BS as well as a radar exploiting special structure of 802.11ad preamble. We embed the radar functionality within standards-compliant operations that do not modify the core structure of the frames beyond what is defined by the 802.11ad protocol. This not only reduces the beam training time, but also ensures scalability with increasing vehicular traffic because radar allows accurate ranging of up to 0.1m at distances up to 200m. We further analyze the ensuing cost-benefit trade-off between the time allotted to the proposed in-band radar and communication modes. Our results reveal 83% reduction on the overhead incurred during the beam training achieved for a specific simulated vehicular scenario over the classical 802.11ad operation.


title: On the structures of hive algebras and tensor product algebras for general linear groups of low rank

abstract: The tensor product algebra TA(n) for the complex general linear group GL(n), introduced by Howe et al., describes the decomposition of tensor products of irreducible polynomial representations of GL(n). Using the hive model for the Littlewood-Richardson coefficients, we provide a finite presentation of the algebra TA(n) for n=2, 3, 4 in terms of generators and relations, thereby giving a description of highest weight vectors of irreducible representations in the tensor products. We also compute the generating function of certain sums of Littlewood-Richardson coefficients.


title: Rado's conjecture and its Baire version

abstract: Rado's Conjecture is a compactness/reflection principle that says any nonspecial tree of height $\omega_1$ has a nonspecial subtree of size $\leq \aleph_1$. Though incompatible with Martin's Axiom, Rado's Conjecture turns out to have many interesting consequences that are consequences of forcing axioms. In this paper, we obtain consistency results concerning Rado's Conjecture and its Baire version. In particular, we show a fragment of PFA, that is the forcing axiom for \emph{Baire Indestructibly proper forcings}, is compatible with the Baire Rado's Conjecture. As a corollary, Baire Rado's Conjecture does not imply Rado's Conjecture. Then we discuss the strength and limitations of the Baire Rado's Conjecture regarding its interaction with simultaneous stationary reflection and some families of weak square principles. Finally we investigate the influence of the Rado's Conjecture on some polarized partition relations.


title: A generalized Schmidt subspace theorem for closed subschemes

abstract: We prove a generalized version of Schmidt's subspace theorem for closed subschemes in general position in terms of suitably defined Seshadri constants with respect to a fixed ample divisor. Our proof builds on previous work by Evertse and Ferretti, Corvaja and Zannier, and others, and uses standard techniques from algebraic geometry such as notions of positivity, blowing-ups and direct image sheaves. As an application, we recover a higher-dimensional Diophantine approximation theorem of K.F. Roth-type due to D. McKinnon and M. Roth with a significantly shortened proof, while simultaneously extending the scope of the use of Seshadri constants in this context in a natural way.


title: Variability of Deeply Embedded Protostars: A New Direction for Star Formation?

abstract: The formation of a star is a dynamic process fed by the gravitational collapse of a molecular cloud core. Theoretical models and observations suggest that the majority of this infalling material settles into a protoplanetary disk before reaching the (proto)star and therefore that disk accretion processes are responsible for the rate at which the (proto)star grows. There is no fundamental reason why infall and disk accretion need to be instantaneously identical. Indeed, even within the disk it might be anticipated that there are regions of strong and weak accretion. Together these facts suggest that (proto)stellar mass assembly should be both secular and stochastic and that the underlying physical processes leading to these time-variable accretion rates should manifest in observable time-dependent accretion luminosity variations.


title: Solitons and black hole in shift symmetric scalar-tensor gravity with cosmological constant

abstract: We demonstrate the existence of static, spherically symmetric globally regular, i.e. solitonic solutions of a shift-symmetric scalar-tensor gravity model with negative cosmological constant. The norm of the Noether current associated to the shift symmetry is finite in the full space-time. We also discuss the corresponding black hole solutions and demonstrate that the interplay between the scalar-tensor coupling and the cosmological constant leads to the existence of new branches of solutions. To linear order in the scalar-tensor coupling, the asymptotic space-time corresponds to an Anti-de Sitter space-time with a non-trivial scalar field on its conformal boundary. This allows the interpretation of our solutions in the context of the AdS/CFT correspondence. Finally, we demonstrate that - for physically relevant, small values of the scalar-tensor coupling - solutions with positive cosmological constant do not exist in our model.


title: Nature of the growth of plasma electrolyte oxidation coating on Aluminum

abstract: The localized dielectric breakdown had been considered as the driving force for growth of plasma electrolytic oxide (PEO) coatings for several decades. However, the TEM study here reveals the dielectric breakdown behavior has little contribution for coating thickening. The presented evidences show the nature of PEO coating growth in all three consecutive stages I-III is ionic migration behavior inside the amorphous alumina layer (AAL) at the coating/matrix interface. The evolution of morphological characterizations in the PEO process is attributed to the interfacial reactions of AAL/alkaline (stage I), AAL/discharge-on-surface (stage II) and AAL/discharge-in-pore (stage III).


title: Covers and partial transversals of Latin squares

abstract: We define a cover of a Latin square to be a set of entries that includes at least one representative of each row, column and symbol. A cover is minimal if it does not contain any smaller cover. A partial transversal is a set of entries that includes at most one representative of each row, column and symbol. A partial transversal is maximal if it is not contained in any larger partial transversal. We explore the relationship between covers and partial transversals. We prove the following: (1) The minimum size of a cover in a Latin square of order $n$ is $n+a$ if and only if the maximum size of a partial transversal is either $n-2a$ or $n-2a+1$. (2) A minimal cover in a Latin square of order $n$ has size at most $\mu_n=3(n+1/2-\sqrt{n+1/4})$. (3) There are infinitely many orders $n$ for which there exists a Latin square having a minimal cover of every size from $n$ to $\mu_n$. (4) Every Latin square of order $n$ has a minimal cover of a size which is asymptotically equal to $\mu_n$. (5) If $1\le k\le n/2$ and $n\ge5$ then there is a Latin square of order $n$ with a maximal partial transversal of size $n-k$. (6) For any $\epsilon>0$, asymptotically almost all Latin squares have no maximal partial transversal of size less than $n-n^{2/3+\epsilon}$.


title: Using the Linear Sigma Model with quarks to describe the QCD phase diagram and to locate the critical end point

abstract: We study the QCD phase diagram using the linear sigma model coupled to quarks. We compute the effective potential at finite temperature and quark chemical potential up to ring diagrams contribution. We show that, provided the values for the pseudo-critical temperature $T_c = 155$ MeV and critical baryon chemical potential $\mu_{Bc} \simeq 1$ GeV, together with the vacuum sigma and pion masses. The model couplings can be fixed and that these in turn help to locate the region where the crossover transition line becomes first order.


title: Cosmic star formation history revealed by AKARI and Hyper Suprime-Cam

abstract: Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshift; the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in 9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, we measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe much more precisely than previous work. In the AKARI NEP wide field, AKARI has obtained deep images in the mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage. To rectify the situation, we used the newly advent Subaru telescope's Hyper Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the AKARI NEP wide field. With this deep and wide optical data, we, for the first time, can use the entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from z=0 to z=2.2, all probed by the AKARI satellite.


title: CURE-TSR: Challenging Unreal and Real Environments for Traffic Sign Recognition

abstract: In this paper, we investigate the robustness of traffic sign recognition algorithms under challenging conditions. Existing datasets are limited in terms of their size and challenging condition coverage, which motivated us to generate the Challenging Unreal and Real Environments for Traffic Sign Recognition (CURE-TSR) dataset. It includes more than two million traffic sign images that are based on real-world and simulator data. We benchmark the performance of existing solutions in real-world scenarios and analyze the performance variation with respect to challenging conditions. We show that challenging conditions can decrease the performance of baseline methods significantly, especially if these challenging conditions result in loss or misplacement of spatial information. We also investigate the effect of data augmentation and show that utilization of simulator data along with real-world data enhance the average recognition performance in real-world scenarios. The dataset is publicly available at https://ghassanalregib.com/cure-tsr/.


title: Gravitoelectromagnetism in Teleparallel Equivalent of General Relativity: A New Alternative

abstract: Spaniol and Andrade introduced grvitoelectromagnetism in TEGR by considering superpotentials, times the determinant of tetrads, as the gravitoelectromagnetic fields. However, since this defined gravitoelectromagnetic field strength does not give rise to a complete set of Maxwell-like equations, we propose an alternative definition of the gravitoelectromagnetic field strength: instead of superpotentials, torsions are taken as the gravitoelectromagnetic field strengths. Based on this new proposal we are able to derive a complete set of Maxwell-like equations. We then apply them to obtain explicit expressions of the gravitoelectromagnetic fields both in Schwarzchilds spacetime and for gravitational waves.


title: Free-Boolean independence with amalgamation

abstract: In this paper, we develop the notion of free-Boolean independence in an amalgamation setting. We construct free-Boolean cumulants and show that the vanishing of mixed free-Boolean cumulants is equivalent to our free-Boolean independence with amalgamation. We also provide a characterization of free-Boolean independence by conditions in terms of mixed moments. In addition, we study free-Boolean independence over a $C^*$-algebra and prove a positivity property.


title: On Sub-Packetization and Access Number of Capacity-Achieving PIR Schemes for MDS Coded Non-Colluding Servers

abstract: Consider the problem of private information retrieval (PIR) over a distributed storage system where $M$ records are stored across $N$ servers by using an $[N,K]$ MDS code. For simplicity, this problem is usually referred as the coded PIR problem. In 2016, Banawan and Ulukus designed the first capacity-achieving coded PIR scheme with sub-packetization $KN^{M}$ and access number $MKN^{M}$, where capacity characterizes the minimal download size for retrieving per unit of data, and sub-packetization and access number are two metrics closely related to implementation complexity. In this paper, we focus on minimizing the sub-packetization and the access number for linear capacity-achieving coded PIR schemes. We first determine the lower bounds on sub-packetization and access number, which are $Kn^{M-1}$ and $MKn^{M-1}$, respectively, in the nontrivial cases (i.e. $N\!>\!K\!\geq\!1$ and $M\!>\!1$), where $n\!=\!N/{\rm gcd}(N,K)$. We then design a general linear capacity-achieving coded PIR scheme to simultaneously attain these two bounds, implying tightness of both bounds.


title: Deep Primal-Dual Reinforcement Learning: Accelerating Actor-Critic using Bellman Duality

abstract: We develop a parameterized Primal-Dual $\pi$ Learning method based on deep neural networks for Markov decision process with large state space and off-policy reinforcement learning. In contrast to the popular Q-learning and actor-critic methods that are based on successive approximations to the nonlinear Bellman equation, our method makes primal-dual updates to the policy and value functions utilizing the fundamental linear Bellman duality. Naive parametrization of the primal-dual $\pi$ learning method using deep neural networks would encounter two major challenges: (1) each update requires computing a probability distribution over the state space and is intractable; (2) the iterates are unstable since the parameterized Lagrangian function is no longer linear. We address these challenges by proposing a relaxed Lagrangian formulation with a regularization penalty using the advantage function. We show that the dual policy update step in our method is equivalent to the policy gradient update in the actor-critic method in some special case, while the value updates differ substantially. The main advantage of the primal-dual $\pi$ learning method lies in that the value and policy updates are closely coupled together using the Bellman duality and therefore more informative. Experiments on a simple cart-pole problem show that the algorithm significantly outperforms the one-step temporal-difference actor-critic method, which is the most relevant benchmark method to compare with. We believe that the primal-dual updates to the value and policy functions would expedite the learning process. The proposed methods might open a door to more efficient algorithms and sharper theoretical analysis.


title: The Inverse Problem for Nested Polygonal Relative Equilibria

abstract: We prove that for some potentials (including the Newtonian one, and the potential of Helmholtz vortices in the plane) relative equilibria consisting of two homothetic regular polygons of arbitrary size can only occur if the masses at each polygon are equal. The same result is true for many regular polygons as long as the ratio between the radii of the polygons are sufficient large. Moreover, under these hypotheses, the relative equilibrium always exist.


title: Asymptotic coverage probabilities of bootstrap percentile confidence intervals for constrained parameters

abstract: The asymptotic behaviour of the commonly used bootstrap percentile confidence interval is investigated when the parameters are subject to linear inequality constraints. We concentrate on the important one- and two-sample problems with data generated from general parametric distributions in the natural exponential family. The focus of this paper is on quantifying the coverage probabilities of the parametric bootstrap percentile confidence intervals, in particular their limiting behaviour near boundaries. We propose a local asymptotic framework to study this subtle coverage behaviour. Under this framework, we discover that when the true parameters are on, or close to, the restriction boundary, the asymptotic coverage probabilities can always exceed the nominal level in the one-sample case; however, they can be, remarkably, both under and over the nominal level in the two-sample case. Using illustrative examples, we show that the results provide theoretical justification and guidance on applying the bootstrap percentile method to constrained inference problems.


title: Improved dynamical constraints on the mass of the central black hole in NGC 404

abstract: We determine the dynamical black hole mass in NGC404 including modeling of the nuclear stellar populations. We combine \hst/STIS spectroscopy with WFC3 images to create a local color-\ml~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3$\sigma$ upper limit on the black hole mass of $1.5\times10^5$\Msun.


title: Engineering Entropy for the Inverse Design of Colloidal Crystals from Hard Shapes

abstract: Throughout the physical sciences, entropy stands out as a pivotal but enigmatic concept that, in materials design, often takes a backseat to energy. Here, we demonstrate how to precisely engineer entropy to achieve desired colloidal crystals. We demonstrate the inverse design of hard particles that assemble six different target colloidal crystals due solely to entropy maximization. Our approach efficiently samples $10^8$ particle shapes from 88- and 192-dimensional design spaces to discover thermodynamically optimal shapes. We design particle shapes that self assemble known crystals with optimized thermodynamic stability, as well as new crystal structures with no known atomic or other equivalent.

categories: cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech physics.comp-ph


title: Precise asymptotics for Fisher-KPP fronts

abstract: We consider the one-dimensional Fisher-KPP equation with step-like initial data. Nolen, Roquejoffre, and Ryzhik showed that the solution $u$ converges at long time to a traveling wave $\phi$ at a position $\tilde \sigma(t) = 2t - (3/2)\log t + \alpha_0- 3\sqrt{\pi}/\sqrt{t}$, with error $O(t^{\gamma-1})$ for any $\gamma>0$. With their methods, we find a refined shift $\sigma(t) = \tilde \sigma(t) + \mu_* (\log t)/t + \alpha_1/t$ such that in the frame moving with $\sigma$, the solution $u$ satisfies $u(t,x) = \phi (x) + \psi(x)/t + O(t^{\gamma-3/2})$ for a certain profile $\psi$ independent of initial data. The coefficient $\alpha_1$ depends on initial data, but $\mu_* = 9(5-6\log 2)/8$ is universal, and agrees with a finding of Berestycki, Brunet, and Derrida in a closely-related problem. Furthermore, we predict the asymptotic forms of $\sigma$ and $u$ to arbitrarily high order.


title: Relevance of Packing to Colloidal Self-Assembly

abstract: Since the 1920s, packing arguments have been used to rationalize crystal structures in systems ranging from atomic mixtures to colloidal crystals. Packing arguments have recently been applied to complex nanoparticle structures, where they often, but not always, work. We examine when, if ever, packing is a causal mechanism in hard particle approximations of colloidal crystals. We investigate three crystal structures comprised of their ideal packing shapes. We show that, contrary to expectations, the ordering mechanism cannot be packing, even when the thermodynamically self-assembled structure is the same as that of the densest packing. We also show that the best particle shapes for hard particle colloidal crystals in the infinite pressure limit are imperfect versions of the ideal packing shape.


title: Gathering in the Plane of Location-Aware Robots in the Presence of Spies

abstract: A set of mobile robots (represented as points) is distributed in the Cartesian plane. The collection contains an unknown subset of byzantine robots which are indistinguishable from the reliable ones. The reliable robots need to gather, i.e., arrive to a configuration in which at the same time, all of them occupy the same point on the plane. The robots are equipped with GPS devices and at the beginning of the gathering process they communicate the Cartesian coordinates of their respective positions to the central authority. On the basis of this information, without the knowledge of which robots are faulty, the central authority designs a trajectory for every robot. The central authority aims to provide the trajectories which result in the shortest possible gathering time of the healthy robots. The efficiency of a gathering strategy is measured by its competitive ratio, i.e., the maximal ratio between the time required for gathering achieved by the given trajectories and the optimal time required for gathering in the offline case, i.e., when the faulty robots are known to the central authority in advance. The role of the byzantine robots, controlled by the adversary, is to act so that the gathering is delayed and the resulting competitive ratio is maximized. The objective of our paper is to propose efficient algorithms when the central authority is aware of an upper bound on the number of byzantine robots. We give optimal algorithms for collections of robots known to contain at most one faulty robot. When the proportion of byzantine robots is known to be less than one half or one third, we provide algorithms with small constant competitive ratios. We also propose algorithms with bounded competitive ratio in the case where the proportion of faulty robots is arbitrary.


title: Sandwich Structures from Arbitrary Functions in Group Theory

abstract: Functions between groups with the property that all function con- jugates are inverse preserving are called sandwich morphisms. These maps pre- serve a structure within the group known as the sandwich structure. Sandwich structures are left distributive idempotent left involutary magmas. These pro- vide a generalisation of groups which we call a sandwich. This paper explores sandwiches and their relationship to groups.