title: On isometry groups of self-adjoint traceless and skew-symmetric matrices

abstract: This note is concerned with isometries on the spaces of self-adjoint traceless matrices. We compute the group of isometries with respect to any unitary similarity invariant norm. This completes and extends the result of Nagy on Schatten $p$-norm isometries. Furthermore, we point out that our proof techniques could be applied to obtain an old result concerning isometries on skew-symmetric matrices.


title: Embeddings for $\mathbb{A}$-weakly differentiable functions on domains

abstract: We prove that the critical embedding $\mathrm{W}^{\mathbb{A},1}(B)\hookrightarrow \mathrm{W}^{k-1,\frac{n}{n-1}}$ holds if and only if the $k$-homogeneous, linear differential operator $\mathbb{A}$ on $\mathbb{R}^n$ from $\mathbb{R}^N$ to $\mathbb{R}^m$ has finite dimensional null-space. Here $B$ is a ball in $\mathbb{R}^n$ and $\mathrm{W}^{\mathbb{A},1}(B)$ denotes the space of maps $u\in \mathrm{L}^1(B,\mathbb{R}^N)$ such that the vector valued distribution $\mathbb{A}u$ is an integrable map. The result was previously known only for several examples of $\mathbb{A}$. Our result contrasts the homogeneous embedding in full-space. Namely, Van Schaftingen proved that $\dot{\mathrm{W}}{^{\mathbb{A},1}}(\mathbb{R}^n)\hookrightarrow \dot{\mathrm{W}}{^{k-1,\frac{n}{n-1}}}$ if and only if $\mathbb{A}$ is elliptic and cancelling. We show that this condition is (strictly) implied by $\mathbb{A}$ having finite dimensional null-space.


title: Multidimensional nonlinear geometric optics for transport operators with applications to stable shock formation

abstract: In $n \geq 1$ spatial dimensions, we study the Cauchy problem for a quasilinear transport equation coupled to a quasilinear symmetric hyperbolic subsystem of a rather general type. For an open set (relative to a suitable Sobolev topology) of regular initial data that are close to the data of a simple plane wave, we give a sharp, constructive proof of shock formation in which the transport variable remains bounded but its first-order Cartesian coordinate partial derivatives blow up in finite time. Moreover, we prove that the singularity does not propagate into the symmetric hyperbolic variables: they and their first-order Cartesian coordinate partial derivatives remain bounded, even though they interact with the transport variable all the way up to its singularity. The formation of the singularity is tied to the finite-time degeneration, relative to the Cartesian coordinates, of a system of geometric coordinates adapted to the characteristics of the transport operator. Two crucial features of the proof are that relative to the geometric coordinates, all solution variables remain smooth, and that the finite-time degeneration coincides with the intersection of the characteristics. Compared to prior shock formation results in more than one spatial dimension, in which the blowup occurred in solutions to wave equations, the main new features of the present work are: i) we develop a theory of nonlinear geometric optics for transport operators, which is compatible with the coupling and which allows us to implement a quasilinear geometric vectorfield method, even though the regularity properties of the corresponding eikonal function are less favorable compared to the wave equation case and ii) we allow for a full quasilinear coupling, i.e., the principal coefficients in all equations are allowed to depend on all solution variables.


title: Normal subgroups generated by a single polynomial automorphism

abstract: We study criteria for deciding when the normal subgroup generated by a single polynomial automorphism of $\mathbb{A}^n$ is as large as possible, namely equal to the normal closure of the special linear group in the special automorphism group. In particular, we investigate $m$-triangular automorphisms, i.e. those that can be expressed as a product of affine automorphisms and $m$ triangular automorphisms. Over a field of characteristic zero, we show that every nontrivial $4$-triangular special automorphism generates the entire normal closure of the special linear group in the special tame subgroup, for any dimension $n \geq 2$. This generalizes a result of Furter and Lamy in dimension 2.


title: A Study of AI Population Dynamics with Million-agent Reinforcement Learning

abstract: We conduct an empirical study on discovering the ordered collective dynamics obtained by a population of intelligence agents, driven by million-agent reinforcement learning. Our intention is to put intelligent agents into a simulated natural context and verify if the principles developed in the real world could also be used in understanding an artificially-created intelligent population. To achieve this, we simulate a large-scale predator-prey world, where the laws of the world are designed by only the findings or logical equivalence that have been discovered in nature. We endow the agents with the intelligence based on deep reinforcement learning (DRL). In order to scale the population size up to millions agents, a large-scale DRL training platform with redesigned experience buffer is proposed. Our results show that the population dynamics of AI agents, driven only by each agent's individual self-interest, reveals an ordered pattern that is similar to the Lotka-Volterra model studied in population biology. We further discover the emergent behaviors of collective adaptations in studying how the agents' grouping behaviors will change with the environmental resources. Both of the two findings could be explained by the self-organization theory in nature.


title: Integrated Waveguide Brillouin Laser

abstract: The demand for high-performance chip-scale lasers has driven rapid growth in integrated photonics. The creation of such low-noise laser sources is critical for emerging on-chip applications, ranging from coherent optical communications, photonic microwave oscillators remote sensing and optical rotational sensors. While Brillouin lasers are a promising solution to these challenges, new strategies are needed to create robust, compact, low power and low cost Brillouin laser technologies through wafer-scale integration. To date, chip-scale Brillouin lasers have remained elusive due to the difficulties in realization of these lasers on a commercial integration platform. In this paper, we demonstrate, for the first time, monolithically integrated Brillouin lasers using a wafer-scale process based on an ultra-low loss Si3N4/SiO2 waveguide platform. Cascading of stimulated Brillouin lasing to 10 Stokes orders was observed in an integrated bus-coupled resonator with a loaded Q factor exceeding 28 million. We experimentally quantify the laser performance, including threshold, slope efficiency and cascading dynamics, and compare the results with theory. The large mode volume integrated resonator and gain medium supports a TE-only resonance and unique 2.72 GHz free spectral range, essential for high performance integrated Brillouin lasing. The laser is based on a non-acoustic guiding design that supplies a broad Brillouin gain bandwidth. Characteristics for high performance lasing are demonstrated due to large intra-cavity optical power and low lasing threshold power. Consistent laser performance is reported for multiple chips across multiple wafers. This design lends itself to wafer-scale integration of practical high-yield, highly coherent Brillouin lasers on a chip.


title: Quantum-Enhanced Sensing from Hyper-Entanglement

abstract: Hyperentanglement --- simultaneous entanglement between multiple degrees of freedom of two or more systems --- has been used to enhance quantum information tasks such as quantum communication and photonic quantum computing. Here we show that hyperentanglement can lead to increased quantum advantage in metrology, with contributions from the entanglement in each degree of freedom, allowing for Heisenberg scaling in the precision of parameter estimation. Our experiment employs photon pairs entangled in polarization and spatial degrees of freedom to estimate a small tilt angle of a mirror. Precision limits beyond shot noise are saturated through a simple binary measurement of the polarization state. The broad validity of the dynamics considered here implies that similar strategies based on hyperentanglement can offer improvement in a wide variety of metrological tasks.


title: Differentially Private Mixture of Generative Neural Networks

abstract: Generative models are used in a wide range of applications building on large amounts of contextually rich information. Due to possible privacy violations of the individuals whose data is used to train these models, however, publishing or sharing generative models is not always viable. In this paper, we present a novel technique for privately releasing generative models and entire high-dimensional datasets produced by these models. We model the generator distribution of the training data with a mixture of $k$ generative neural networks. These are trained together and collectively learn the generator distribution of a dataset. Data is divided into $k$ clusters, using a novel differentially private kernel $k$-means, then each cluster is given to separate generative neural networks, such as Restricted Boltzmann Machines or Variational Autoencoders, which are trained only on their own cluster using differentially private gradient descent. We evaluate our approach using the MNIST dataset, as well as call detail records and transit datasets, showing that it produces realistic synthetic samples, which can also be used to accurately compute arbitrary number of counting queries.


title: High thermoelectric performance of graphite nanofibers

abstract: Graphite nanofibers (GNFs) have been demonstrated to be a promising material for hydrogen storage and heat management in electronic devices. Here, by means of first-principles and transport simulations, we show that GNFs can also be an excellent material for thermoelectric applications thanks to the interlayer weak van der Waals interaction that induces low thermal conductance and a step-like shape in the electronic transmission with mini-gaps, which are necessary ingredients to achieve high thermoelectric performance. This study unveils that the platelet form of GNFs in which graphite layers are perpendicular to the fiber axis can exhibit outstanding thermoelectric properties with a figure of merit ZT reaching 3.55 in a 0.5 nm diameter fiber and 1.1 for a 1.1 nm diameter one. Interestingly, by introducing 14C isotope doping, ZT can even be enhanced up to more than 5, and more than 8 if we include the effect of finite phonon mean-free path, which demonstrates the amazing thermoelectric potential of GNFs.


title: Inverse square law isothermal property in relativistic charged static distributions

abstract: We analyse the impact of the inverse square law fall-off of the energy density in a charged isotropic spherically symmetric fluid. Initially we impose a linear barotropic equation of state $p=\alpha \rho$ but this leads to an intractable differential equation. Next we consider the neutral isothermal metric of Saslaw, Maharaj and Dadhich (1996) in an electric field and the usual inverse square law of energy density and pressure results thus preserving the equation of state. Additionally, we discard a linear equation of state and endeavour to find new classes of solutions with the inverse square law fall off of density. Certain prescribed forms of the spatial and temporal gravitational forms result in new exact solutions. An interesting result that emerges is that while isothermal fluid spheres are unbounded in the neutral case, this is not so when charge is involved. Indeed it was found that barotropic equations of state exist and hypersurfaces of vanishing pressure exist establishing a boundary in practically all models. One model was studied in depth and found to satisfy other elementary requirements for physical admissability such as a subluminal sound speed as well as gravitational surface redshifts smaller than 2. The Buchdahl (1959), Bohmer and Harko (2007) and Andreasson (2009) mass-radius bounds were also found to be satisfied. Graphical plots utilising constants selected from the boundary conditions established that the model displayed characteristics consistent with physically viable models.


title: Visualizations for an Explainable Planning Agent

abstract: In this paper, we report on the visualization capabilities of an Explainable AI Planning (XAIP) agent that can support human in the loop decision making. Imposing transparency and explainability requirements on such agents is especially important in order to establish trust and common ground with the end-to-end automated planning system. Visualizing the agent's internal decision-making processes is a crucial step towards achieving this. This may include externalizing the "brain" of the agent -- starting from its sensory inputs, to progressively higher order decisions made by it in order to drive its planning components. We also show how the planner can bootstrap on the latest techniques in explainable planning to cast plan visualization as a plan explanation problem, and thus provide concise model-based visualization of its plans. We demonstrate these functionalities in the context of the automated planning components of a smart assistant in an instrumented meeting space.


title: Recurrent Saliency Transformation Network: Incorporating Multi-Stage Visual Cues for Small Organ Segmentation

abstract: We aim at segmenting small organs (e.g., the pancreas) from abdominal CT scans. As the target often occupies a relatively small region in the input image, deep neural networks can be easily confused by the complex and variable background. To alleviate this, researchers proposed a coarse-to-fine approach, which used prediction from the first (coarse) stage to indicate a smaller input region for the second (fine) stage. Despite its effectiveness, this algorithm dealt with two stages individually, which lacked optimizing a global energy function, and limited its ability to incorporate multi-stage visual cues. Missing contextual information led to unsatisfying convergence in iterations, and that the fine stage sometimes produced even lower segmentation accuracy than the coarse stage. This paper presents a Recurrent Saliency Transformation Network. The key innovation is a saliency transformation module, which repeatedly converts the segmentation probability map from the previous iteration as spatial weights and applies these weights to the current iteration. This brings us two-fold benefits. In training, it allows joint optimization over the deep networks dealing with different input scales. In testing, it propagates multi-stage visual information throughout iterations to improve segmentation accuracy. Experiments in the NIH pancreas segmentation dataset demonstrate the state-of-the-art accuracy, which outperforms the previous best by an average of over 2%. Much higher accuracies are also reported on several small organs in a larger dataset collected by ourselves. In addition, our approach enjoys better convergence properties, making it more efficient and reliable in practice.


title: Tracing the stellar component of low surface brightness Milky Way Dwarf Galaxies to their outskirts I: Sextans

abstract: We present results from deep and very spatially extended CTIO/DECam $g$ and $r$ photometry (reaching out to $\sim$ 2 mag below the oldest MSTO and covering $\sim$ 20 deg$^2$) around the Sextans dSph. We use this data-set to study the structural properties of Sextans overall stellar population and its different stellar evolutionary phases, as well as to search for signs of tidal disturbance from the MW, which would indicate departure from dynamical equilibrium. We perform the most accurate structural analysis to-date of Sextans' stellar components by applying Bayesian MCMC methods to the individual stars' positions. Surface density maps are built by decontaminating the sample through a matched filter analysis of the CMD, and then analysed for departures from axisymmetry. Sextans is found to be considerably less spatially extended than early studies suggested. No significant distortions or tidal disturbances are found down to a surface brightness of $\sim$ 31.8 mag/arcsec$^{-2}$ in V-band. We identify an overdensity in the central regions that may correspond to previously reported kinematic substructure(s). In agreement with previous findings, old $\&$ metal-poor stars such as BHB stars cover a much larger area than stars in other evolutionary phases, and bright BSs are less spatially extended than faint ones. However, the different spatial distribution of bright and faint BSs appears consistent with the general age/metallicity gradients found in Sextans' stellar component. This is compatible with BSs having formed by evolution of binaries and not necessarily due to the presence of a central disrupted globular cluster, as suggested in the literature. We provide structural parameters for the various populations analyzed and make publicly available the photometric catalogue of point-sources as well as a catalogue of literature spectroscopic measurements with updated membership probabilities.


title: Photonic Counterparts of Cooper Pairs

abstract: The microscopic theory of superconductivity raised the disruptive idea that electrons couple through the elusive exchange of virtual phonons, overcoming the strong Coulomb repulsion to form Cooper pairs. Light is also known to interact with atomic vibrations, as for example in the Raman effect. We show that photon pairs exchange virtual vibrations in transparent media, leading to an effective photon-photon interaction identical to that for electrons in BCS theory of supercondutivity, in spite of the fact that photons are bosons. In this scenario, photons may exchange energy without matching a quantum of vibration of the medium. As a result, pair correlations for photons scattered away from the Raman resonances are expected to be enhanced. Experimental demonstration of this effect is provided here by time correlated Raman measurements in different media. The experimental data confirm our theoretical interpretation of a photonic Cooper pairing, without the need for any fitting parameters.


title: Stellar Populations in a semi-analytic model I: bulges of Milky Way-like galaxies

abstract: We study the stellar populations of bulges of Milky Way-like (MW-like) galaxies with the aim of identifying the physical processes involved in the formation of the bulge of our Galaxy. We use the semi-analytic model of galaxy formation and evolution SAG adapted to this aim; this kind of models can trace the properties of galaxies and their components like stellar discs, bulges and halos, but resolution limits prevent them from reaching the scale of stellar populations (SPs). Properties of groups of stars formed during single star formation events are stored and tracked in the model and results are compared with observations of stars in the galactic bulge. MW-like galaxies are selected using two different criteria. One of them considers intrinsic photo-metric properties and the second is focused on the cosmological context of the local group of galaxies (LG). We compare our model results with spectroscopic and photometric stellar metallicity distributions. We find that 87% of stars in bulges of MWtype galaxies in our model are accreted and formed in starbursts during disc instability events. Mergers contribute to 13% of the mass budget of the bulge and are responsible for the low metallicity tail of the distribution. Abundance ratios of {\alpha} elements with respect to iron, [{\alpha}/Fe], are measured in SPs of model galaxies. The patterns found in the model for SPs with different origins help to explain the lack of a gradient of [{\alpha}/Fe] ratios in observed stars along the minor axis of the bulge.


title: Permanent spin currents in cavity-qubit systems

abstract: In a recent remarkable experiment [P. Roushan et al., Nature Physics 13, 146 (2017)], a spin current in an architecture of three superconducting qubits was produced during a few microseconds by creating synthetic magnetic fields. The life-time of the current was set by the typical dissipative mechanisms that occur in those systems. We propose a scheme for the generation of permanent currents, even in the presence of such imperfections, and scalable to larger system sizes. It relies on striking a subtle balance between multiple nonequilibrium drives and the dissipation mechanisms, in order to engineer and stimulate chiral excited states which can carry current.


title: Polishability of some groups of interval and circle diffeomorphisms

abstract: Let $M=I$ or $M=\mathbb{S}^1$ and let $k\geq 1$. We exhibit a new infinite class of Polish groups by showing that each group $\mathop{\rm Diff}_+^{k+AC}(M)$, consisting of those $C^k$ diffeomorphisms whose $k$-th derivative is absolutely continuous, admits a natural Polish group topology which refines the subspace topology inherited from $\mathop{\rm Diff}_+^k(M)$. By contrast, the group $\mathop{\rm Diff}_+^{1+BV}(M)$, consisting of $C^1$ diffeomorphisms whose derivative has bounded variation, admits no Polish group topology whatsoever.


title: Workflow Complexity for Collaborative Interactions: Where are the Metrics? -- A Challenge

abstract: In this paper, we introduce the problem of denoting and deriving the complexity of workflows (plans, schedules) in collaborative, planner-assisted settings where humans and agents are trying to jointly solve a task. The interactions -- and hence the workflows that connect the human and the agents -- may differ according to the domain and the kind of agents. We adapt insights from prior work in human-agent teaming and workflow analysis to suggest metrics for workflow complexity. The main motivation behind this work is to highlight metrics for human comprehensibility of plans and schedules. The planning community has seen its fair share of work on the synthesis of plans that take diversity into account -- what value do such plans hold if their generation is not guided at least in part by metrics that reflect the ease of engaging with and using those plans?


title: Partial Regularity for a Nonlinear Sigma Model with Gravitino in Higher Dimensions

abstract: We study the regularity problem of the nonlinear sigma model with gravitino fields in higher dimensions. After setting up the geometric model, we derive the Euler--Lagrange equations and consider the regularity of weak solutions defined in suitable Sobolev spaces. We show that any weak solution is actually smooth under some smallness assumption for certain Morrey norms. By assuming some higher integrability of the vector spinor, we can show a partial regularity result for stationary solutions, provided the gravitino is critical, which means that the corresponding supercurrent vanishes. Moreover, in dimension less than 6, partial regularity holds for stationary solutions with respect to general gravitino fields.


title: Phase and Power Control in the RF Magnetron Power Stations of Superconducting Accelerators

abstract: Phase and power control methods that satisfy the requirements of superconducting accelerators to magnetron RF sources were considered by a simplified kinetic model of a magnetron driven by a resonant injected signal. The model predicting and explaining stable, low noise operation of the tube below the threshold of self-excitation (the Hatrree voltage in free run mode) at a highest efficiency, a wide range of power control and a wide-band phase control was well verified in experiments demonstrating capabilities of the magnetron transmitters for powering of state of the art superconducting accelerators. Descriptions of the kinetic model, the experimental verification and a conceptual scheme of the highly-efficient magnetron RF transmitter for the accelerators are presented and discussed.


title: Control of Ultracold Photodissociation with Magnetic Fields

abstract: Photodissociation of a molecule produces a spatial distribution of photofragments determined by the molecular structure and the characteristics of the dissociating light. Performing this basic chemical reaction at ultracold temperatures allows its quantum mechanical features to dominate. In this regime, weak applied fields can be used to control the reaction. Here, we photodissociate ultracold diatomic strontium in magnetic fields below 10 G and observe striking changes in photofragment angular distributions. The observations are in excellent qualitative agreement with a multichannel quantum chemistry model that includes nonadiabatic effects and predicts strong mixing of partial waves in the photofragment energy continuum. The experiment is enabled by precise quantum-state control of the molecules.


title: Coordinates Adapted to Vector Fields: Canonical Coordinates

abstract: Given a finite collection of $C^1$ vector fields on a $C^2$ manifold which span the tangent space at every point, we consider the question of when there is locally a coordinate system in which these vector fields have a higher level of smoothness. For example, when is there a coordinate system in which the vector fields are smooth, or real analytic, or have Zygmund regularity of some finite order? We address this question in a quantitative way, which strengthens and generalizes previous works on the quantitative theory of sub-Riemannian (aka Carnot-Carath\'eodory) geometry due to Nagel, Stein, and Wainger, Tao and Wright, the second author, and others. Furthermore, we provide a diffeomorphism invariant version of these theories. This is the first part in a three part series of papers. In this paper, we study a particular coordinate system adapted to a collection of vector fields (sometimes called canonical coordinates) and present results related to the above questions which are not quite sharp; these results from the backbone of the series. The methods of this paper are based on techniques from ODEs. In the second paper, we use additional methods from PDEs to obtain the sharp results. In the third paper, we prove results concerning real analyticity and use methods from ODEs.


title: Mirror-assisted coherent backscattering from the Mollow sidebands

abstract: In front of a mirror, the radiation of weakly driven large disordered clouds presents an interference fringe in the backward direction, on top of an incoherent background. Although strongly driven atoms usually present little coherent scattering, we here show that the mirror-assisted version can produce high contrast fringes, for arbitrarily high saturation parameters. The contrast of the fringes oscillates with the Rabi frequency of the atomic transition and the distance between the mirror and the atoms, due to the coherent interference between the carrier and the Mollow sidebands of the saturated resonant fluorescence spectrum emitted by the atoms. The setup thus represents a powerful platform to study the spectral properties of ensembles of correlated scatterers.


title: State-Secrecy Codes for Networked Linear Systems

abstract: In this paper, we study the problem of remote state estimation, in the presence of a passive eavesdropper. An authorized user estimates the state of an unstable linear plant, based on the packets received from a sensor, while the packets may also be intercepted by the eavesdropper. Our goal is to design a coding scheme at the sensor, which encodes the state information, in order to impair the eavesdropper's estimation performance, while enabling the user to successfully decode the sent messages. We introduce a novel class of codes, termed State-Secrecy Codes, which use acknowledgment signals from the user and apply linear time-varying transformations to the current and previously received states. By exploiting the properties of the system's process noise, the channel physical model and the dynamics, these codes manage to be fast, efficient and, thus, suitable for real-time dynamical systems. We prove that under minimal conditions, State-Secrecy Codes achieve perfect secrecy, namely the eavesdropper's estimation error grows unbounded almost surely, while the user's estimation performance is optimal. These conditions only require that at least once, the user receives the corresponding packet while the eavesdropper fails to intercept it. Even one occurrence of this event renders the eavesdropper's error unbounded with asymptotically optimal rate of increase. State-Secrecy Codes are provided and studied for two cases, i) when direct state measurements are available, and ii) when we only have output measurements. The theoretical results are illustrated in simulations.


title: Negative-mass instability of the spin and motion of an atomic gas driven by optical cavity backaction

abstract: We realize a spin-orbit interaction between the collective spin precession and center-of-mass motion of a trapped ultracold atomic gas, mediated by spin- and position-dependent dispersive coupling to a driven optical cavity. The collective spin, precessing near its highest-energy state in an applied magnetic field, can be approximated as a negative-mass harmonic oscillator. When the Larmor precession and mechanical motion are nearly resonant, cavity mediated coupling leads to a negative-mass instability, driving exponential growth of a correlated mode of the hybrid system. We observe this growth imprinted on modulations of the cavity field and estimate the full covariance of the resulting two-mode state by observing its transient decay during subsequent free evolution.


title: Effective g factor of low-density two-dimensional holes in a Ge quantum well

abstract: We report measurements of the effective $g$ factor of low-density two-dimensional holes in a Ge quantum well. Using the temperature dependence of the Shubnikov-de Haas oscillations, we extract the effective $g$ factor in a magnetic field perpendicular to the sample surface. Very large values of the effective $g$ factor, ranging from $\sim13$ to $\sim28$, are observed in the density range of $1.4\times10^{10}$ cm$^{-2}$ to $1.4\times10^{11}$ cm$^{-2}$. When the magnetic field is oriented parallel to the sample surface, the effective $g$ factor is obtained from a protrusion in the magneto-resistance data that signifies full spin polarization. In the latter orientation, a small effective $g$ factor, $\sim1.3-1.4$, is measured in the density range of $1.5\times10^{10}$ cm$^{-2}$ to $2\times10^{10}$ cm$^{-2}$. This very strong anisotropy is consistent with theoretical predictions and previous measurements in other 2D hole systems, such as InGaAs and GaSb.


title: TASI Lectures on Collider Physics

abstract: These lectures provide an introduction to the physics of particle colliders. Topics covered include a quantitative examination of the design and operational parameters of Large Hadron Collider, kinematics and observables at colliders, such as rapidity and transverse mass, and properties of distributions, such as Jacobian peaks. In addition, the lectures provide a practical introduction to the decay modes and signatures of important composite and elementary and particles in the Standard Model, from pions to the Higgs boson. The aim of these lectures is provide a bridge between theoretical and experimental particle physics. Whenever possible, results are derived using intuitive arguments and estimates rather than precision calculations.


title: Modeling a Propagating Sawtooth Flare Ribbon Structure as a Tearing Mode in the Presence of Velocity Shear

abstract: On April 18, 2014 (SOL2014-04-18T13:03) an M-class flare was observed by IRIS. The associated flare ribbon contained a quasi-periodic sawtooth pattern that was observed to propagate along the ribbon, perpendicular to the IRIS spectral slit, with a phase velocity of $\sim 15$ km s$^{-1}$. This motion resulted in periodicities in both intensity and Doppler velocity along the slit. These periodicities were reported to be approximately $\pm0.5$ arcseconds in position and $\pm20$ km s$^{-1}$ in velocity and were measured to be $\sim180^{\circ}$ out of phase with one another. This quasi-periodic behavior has been attributed by others to bursty or patchy reconnection and slipping occurring during three-dimensional magnetic reconnection. While able to account for periodicities in both intensity and Doppler velocity these suggestions do not explicitly account for the phase velocity of the entire sawtooth structure, or for the relative phasing of the oscillations. Here we propose that the observations can be explained by a tearing mode instability occurring at a current sheet across which there is also a velocity shear. Using a linear model of this instability we reproduce the relative phase of the oscillations, as well as the phase velocity of the sawtooth structure. We suggest a geometry and local plasma parameters for the April 18 flare which would support our hypothesis. Under this proposal the combined spectral and spatial IRIS observations of this flare may provide the most compelling evidence to date of a tearing mode occurring in the solar magnetic field.


title: Axion dark matter in a $3-3-1$ model

abstract: Slightly extending a right-handed neutrino version of the $3-3-1$ model, we show that it is not only possible to solve the strong CP problem but also to give the total dark matter abundance reported by the Planck collaboration. Specifically, we consider the possibility of introducing a $3-3-1$ scalar singlet to implement a gravity stable Peccei-Quinn mechanism in this model. Remarkably, for allowed regions of the parameter space, the arising axions with masses $m_a\approx$ meV can both make up the total dark matter relic density through non-thermal production mechanisms and be very close to the region to be explored by the IAXO helioscope.


title: Dynamic tunable notch filters for the Antarctic Impulsive Transient Antenna (ANITA)

abstract: The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.


title: Squares, ascent paths, and chain conditions

abstract: With the help of various square principles, we obtain results concerning the consistency strength of several statements about trees containing ascent paths, special trees, and strong chain conditions. Building on a result that shows that Todor\v{c}evi\'{c}'s principle $\square(\kappa)$ implies an indexed version of $\square(\kappa,\lambda)$, we show that for all infinite, regular cardinals $\lambda<\kappa$, the principle $\square(\kappa)$ implies the existence of a $\kappa$-Aronszajn tree containing a $\lambda$-ascent path. We then provide a complete picture of the consistency strengths of statements relating the interactions of trees with ascent paths and special trees. As a part of this analysis, we construct a model of set theory in which $\aleph_2$-Aronszajn trees exist and all such trees contain $\aleph_0$-ascent paths. Finally, we use our techniques to show that the assumption that the $\kappa$-Knaster property is countably productive and the assumption that every $\kappa$-Knaster partial order is $\kappa$-stationarily layered both imply the failure of $\square(\kappa)$.


title: Petz recovery versus matrix reconstruction

abstract: The reconstruction of the state of a multipartite quantum mechanical system represents a fundamental task in quantum information science. At its most basic, it concerns a state of a bipartite quantum system whose subsystems are subjected to local operations. We compare two different methods for obtaining the original state from the state resulting from the action of these operations. The first method involves quantum operations called Petz recovery maps, acting locally on the two subsystems. The second method is called matrix (or state) reconstruction and involves local, linear maps which are not necessarily completely positive. Moreover, we compare the quantities on which the maps employed in the two methods depend. We show that any state which admits Petz recovery also admits state reconstruction. However, the latter is successful for a strictly larger set of states. We also compare these methods in the context of a finite spin chain. Here, the state of a finite spin chain is reconstructed from the reduced states of a few neighbouring spins. In this setting, state reconstruction is the same as the MPO (i.e. matrix product operator) reconstruction proposed by Baumgratz et al. in [Phys. Rev. Lett. 111, 020401 (2013)]. Finally, we generalize both these methods so that they employ long-range measurements instead of relying solely on short-range correlations embodied in such local reduced states. Long-range measurements enable the reconstruction of states which cannot be reconstructed from measurements of local few-body observables alone and hereby we improve existing methods for quantum state tomography of quantum many-body systems.


title: Bayesian Analysis of Hot Jupiter Radius Anomalies: Evidence for Ohmic Dissipation?

abstract: The cause of hot Jupiter radius inflation, where giant planets with $T_{\rm eq}$ $>1000$ K are significantly larger than expected, is an open question and the subject of many proposed explanations. Rather than examine these models individually, this work seeks to characterize the anomalous heating as a function of incident flux, $\epsilon(F)$, needed to inflate the population of planets to their observed sizes. We then compare that result to theoretical predictions for various models. We examine the population of about 300 giant planets with well-determined masses and radii and apply thermal evolution and Bayesian statistical models to infer the anomalous power as a function of incident flux that best reproduces the observed radii. First, we observe that the inflation of planets below about $M=0.5 \;\rm{M}_\rm{J}$ appears very different than their higher mass counterparts, perhaps as the result of mass loss or an inefficient heating mechanism. As such, we exclude planets below this threshold. Next, we show with strong significance that $\epsilon(F)$ increases with $T_{\rm{eq}}$ towards a maximum of $\sim 2.5\%$ at $T_{\rm{eq}} \approx 1500$ K, and then decreases as temperatures increase further, falling to $\sim0.2\%$ at $T_\rm{eff}= 2500$ K. This high-flux decrease in inflation efficiency was predicted by the Ohmic dissipation model of giant planet inflation but not other models. We also explicitly check the thermal tides model and find that it predicts far more variance in radii than is observed. Thus, our results provide evidence for the Ohmic dissipation model and a functional form for $\epsilon(F)$ that any future theories of hot Jupiter radii can be tested against.


title: Comparative study of electronic and magnetic properties of $M$Pc ($M$ = Fe, Co) molecules physisorbed on 2D MoS$_2$ and graphene

abstract: In this paper, we have done a comparative study of electronic and magnetic properties of iron phthalocyanine (FePc) and cobalt phthalocyanine (CoPc) molecules physisorbed on monolayer of MoS$_2$ and graphene by using density functional theory. Various different types of physisorption sites have been considered for both surfaces. Our calculations reveal that the $M$Pc molecules prefer the S-top position on MoS$_2$. However, on graphene, FePc molecule prefers the bridge position while CoPc molecule prefers the top position. The $M$Pc molecules are physisorbed strongly on the MoS$_2$ surface than the graphene ($\sim$ 2.5 eV higher physisorption energy). Analysis of magnetic properties indicates the presence of strong spin dipole moment opposite to the spin moment and hence a huge reduction of effective spin moment can be observed. Our calculations of magnetic anisotropy energies using both variational approach and $2^{nd}$ order perturbation approach indicate no significant changes after physisorption. In case of FePc, an out-of-plane easy axis and in case of CoPc, an in-plane easy axis can be seen. Calculations of work function indicate a reduction of MoS$_2$ work function $\sim$ 1 eV due to physisorption of $M$Pc molecules while it does not change significantly in case of graphene.


title: A trapped field of 17.7 T in a stack of high temperature superconducting tape

abstract: High temperature superconducting (HTS) tape can be cut and stacked to generate large magnetic fields at cryogenic temperatures after inducing persistent currents in the superconducting layers. A field of 17.7 T was trapped between two stacks of HTS tape at 8 K with no external mechanical reinforcement. 17.6 T could be sustained when warming the stack up to 14 K. A new type of hybrid stack was used consisting of a 12 mm square insert stack embedded inside a larger 34.4 mm diameter stack made from different tape. The magnetic field generated is the largest for any trapped field magnet reported and 30% greater than previously achieved in a stack of HTS tapes. Such stacks are being considered for superconducting motors as rotor field poles where the cryogenic penalty is justified by the increased power to weight ratio. The sample reported can be considered the strongest permanent magnet ever created.


title: Heavy Quarkonium ($1^{--}$) Meson-Hybrid Mixing from QCD Sum Rules

abstract: We use QCD Laplace sum-rules to explore mixing between conventional mesons and hybrids in the heavy quarkonium vector $J^{PC}\!=\!1^{--}$ channel. Our cross-correlator includes perturbation theory and contributions proportional to the four-dimensional and six-dimensional gluon condensates. We input experimentally determined charmonium and bottomonium hadron masses into both single and multi-resonance models in order to test them for conventional meson and hybrid components. In the charmonium sector we find evidence for meson-hybrid mixing in the $J/\psi$ and a $\approx4.3$ GeV resonance. In the bottomonium sector, we find that the $\Upsilon(1S)$, $\Upsilon(2S)$, $\Upsilon(3S)$, and $\Upsilon(4S)$ all exhibit mixing.


title: Data-Efficient Multirobot, Multitask Transfer Learning for Trajectory Tracking

abstract: Transfer learning has the potential to reduce the burden of data collection and to decrease the unavoidable risks of the training phase. In this letter, we introduce a multirobot, multitask transfer learning framework that allows a system to complete a task by learning from a few demonstrations of another task executed on another system. We focus on the trajectory tracking problem where each trajectory represents a different task, since many robotic tasks can be described as a trajectory tracking problem. The proposed multirobot transfer learning framework is based on a combined $\mathcal{L}_1$ adaptive control and an iterative learning control approach. The key idea is that the adaptive controller forces dynamically different systems to behave as a specified reference model. The proposed multitask transfer learning framework uses theoretical control results (e.g., the concept of vector relative degree) to learn a map from desired trajectories to the inputs that make the system track these trajectories with high accuracy. This map is used to calculate the inputs for a new, unseen trajectory. Experimental results using two different quadrotor platforms and six different trajectories show that, on average, the proposed framework reduces the first-iteration tracking error by 74% when information from tracking a different single trajectory on a different quadrotor is utilized.


title: Distortion of Magnetic Fields in a Starless Core II: 3D Magnetic Field Structure of FeSt 1-457

abstract: Three dimensional (3D) magnetic field information on molecular clouds and cores is important for revealing their kinematical stability (magnetic support) against gravity which is fundamental for studying the initial conditions of star formation. In the present study, the 3D magnetic field structure of the dense starless core FeSt 1-457 is determined based on the near-infrared polarimetric observations of the dichroic polarization of background stars and simple 3D modeling. With an obtained angle of line-of-sight magnetic inclination axis $\theta_{\rm inc}$ of $45^{\circ}\pm10^{\circ}$ and previously determined plane-of-sky magnetic field strength $B_{\rm pol}$ of $23.8\pm12.1$ $\mu{\rm G}$, the total magnetic field strength for FeSt 1-457 is derived to be $33.7\pm18.0$ $\mu{\rm G}$. The critical mass of FeSt 1-457, evaluated using both magnetic and thermal/turbulent support is ${M}_{\rm cr} = 3.70\pm0.92$ ${\rm M}_{\odot}$, which is identical to the observed core mass, $M_{\rm core}=3.55\pm0.75$ ${\rm M}_{\odot}$. We thus conclude that the stability of FeSt 1-457 is in a condition close to the critical state. Without infalling gas motion and no associated young stars, the core is regarded to be in the earliest stage of star formation, i.e., the stage just before the onset of dynamical collapse following the attainment of a supercritical condition. These properties would make FeSt 1-457 one of the best starless cores for future studies of the initial conditions of star formation.


title: Efficient Computation of Multiple Density-Based Clustering Hierarchies

abstract: HDBSCAN*, a state-of-the-art density-based hierarchical clustering method, produces a hierarchical organization of clusters in a dataset w.r.t. a parameter mpts. While the performance of HDBSCAN* is robust w.r.t. mpts in the sense that a small change in mpts typically leads to only a small or no change in the clustering structure, choosing a "good" mpts value can be challenging: depending on the data distribution, a high or low value for mpts may be more appropriate, and certain data clusters may reveal themselves at different values of mpts. To explore results for a range of mpts values, however, one has to run HDBSCAN* for each value in the range independently, which is computationally inefficient. In this paper, we propose an efficient approach to compute all HDBSCAN* hierarchies for a range of mpts values by replacing the graph used by HDBSCAN* with a much smaller graph that is guaranteed to contain the required information. An extensive experimental evaluation shows that with our approach one can obtain over one hundred hierarchies for the computational cost equivalent to running HDBSCAN* about 2 times.


title: Normalized Direction-preserving Adam

abstract: Adaptive optimization algorithms, such as Adam and RMSprop, have shown better optimization performance than stochastic gradient descent (SGD) in some scenarios. However, recent studies show that they often lead to worse generalization performance than SGD, especially for training deep neural networks (DNNs). In this work, we identify the reasons that Adam generalizes worse than SGD, and develop a variant of Adam to eliminate the generalization gap. The proposed method, normalized direction-preserving Adam (ND-Adam), enables more precise control of the direction and step size for updating weight vectors, leading to significantly improved generalization performance. Following a similar rationale, we further improve the generalization performance in classification tasks by regularizing the softmax logits. By bridging the gap between SGD and Adam, we also hope to shed light on why certain optimization algorithms generalize better than others.


title: Spin dependence of {\eta} meson production in proton-proton collisions close to threshold

abstract: Taking advantage of the high acceptance and axial symmetry of the WASA-at-COSY detector, and the high degree of the polarized proton beam of COSY, the reaction pp{\to} pp{\eta} has been measured close to threshold to explore the analyzing power Ay. The angular distribution of Ay is determined with the precision improved by more than one order of magnitude with respect to previous results allowing a first accurate comparison with theoretical predictions. The determined analyzing power is consistent with zero for an excess energy of Q = 15 MeV signaling s wave production with no evidence for higher partial waves. At Q = 72 MeV the data reveals strong interference of P s and P p partial waves and cancellation of (P p) and Ss*Sd contributions. These results rule out the presently available theoretical predictions for the production mechanism of the {\eta} meson.


title: Distinguishing between convergent evolution and violation of the molecular clock

abstract: We give a non-technical introduction to convergence-divergence models, a new modeling approach for phylogenetic data that allows for the usual divergence of species post speciation but also allows for species to converge, i.e. become more similar over time. By examining the $3$-taxon case in some detail we illustrate that phylogeneticists have been "spoiled" in the sense of not having to think about the structural parameters in their models by virtue of the strong assumption that evolution is treelike. We show that there are not always good statistical reasons to prefer the usual class of treelike models over more general convergence-divergence models. Specifically we show many $3$-taxon datasets can be equally well explained by supposing violation of the molecular clock due to change in the rate of evolution along different edges, or by keeping the assumption of a constant rate of evolution but instead assuming that evolution is not a purely divergent process. Given the abundance of evidence that evolution is not strictly treelike, our discussion is an illustration that as phylogeneticists we often need to think clearly about the structural form of the models we use.


title: Ignoring Distractors in the Absence of Labels: Optimal Linear Projection to Remove False Positives During Anomaly Detection

abstract: In the anomaly detection setting, the native feature embedding can be a crucial source of bias. We present a technique, Feature Omission using Context in Unsupervised Settings (FOCUS) to learn a feature mapping that is invariant to changes exemplified in training sets while retaining as much descriptive power as possible. While this method could apply to many unsupervised settings, we focus on applications in anomaly detection, where little task-labeled data is available. Our algorithm requires only non-anomalous sets of data, and does not require that the contexts in the training sets match the context of the test set. By maximizing within-set variance and minimizing between-set variance, we are able to identify and remove distracting features while retaining fidelity to the descriptiveness needed at test time. In the linear case, our formulation reduces to a generalized eigenvalue problem that can be solved quickly and applied to test sets outside the context of the training sets. This technique allows us to align technical definitions of anomaly detection with human definitions through appropriate mappings of the feature space. We demonstrate that this method is able to remove uninformative parts of the feature space for the anomaly detection setting.


title: A Computational Model of Afterimages based on Simultaneous and Successive Contrasts

abstract: Negative afterimage appears in our vision when we shift our gaze from an over stimulated original image to a new area with a uniform color. The colors of negative afterimages differ from the old stimulating colors in the original image when the color in the new area is either neutral or chromatic. The interaction between stimulating colors in the test and inducing field in the original image changes our color perception due to simultaneous contrast, and the interaction between changed colors perceived in the previously-viewed field and the color in the currently-viewed field also affects our perception of colors in negative afterimages due to successive contrast. Based on these observations we propose a computational model to estimate colors of negative afterimages in more general cases where the original stimulating color in the test field is chromatic, and the original stimulating color in the inducing field and the new stimulating color can be either neutral or chromatic. We validate our model with human experiments.


title: An Operational Semantic Basis for OpenMP Race Analysis

abstract: OpenMP is the de facto standard to exploit the on-node parallelism in new generation supercomputers.Despite its overall ease of use, even expert users are known to create OpenMP programs that harbor concurrency errors, of which one of the most insidious of errors are {\em data races}.OpenMP is also a rapidly evolving standard, which means that future data races may be introduced within unfamiliar contexts.A simple and rigorous operational semantics for OpenMP can help build reliable race checkers and ward off future errors through programmer education and better tooling.This paper's key contribution is a simple operational semantics for OpenMP, with primitive events matching those generated by today's popular OpenMP runtimes and tracing methods such as OMPT.This makes our operational semantics more than a theoretical document for intellectual edification; it can serve as a blueprint for OpenMP event capture and tool building.We back this statement by summarizing the workings of a new data race checker for OpenMP being built based on this semantics.The larger purpose served by our semantics is to serve the needs of the OpenMP community with regard to their contemplated extensions to OpenMP, as well as future tooling efforts.


title: On the Calder\'{o}n-Zygmund structure of Petermichl's kernel. Weighted inequalities

abstract: We show that Petermichl's dyadic operator $\mathcal{P}$ (S. Petermichl (2000), Dyadic shifts and a logarithmic estimate for Hankel operators with matrix symbol) is a Calder\'{o}n-Zygmund type operator on an adequate metric normal space of homogeneous type. As a consequence of a general result on spaces of homogeneous type, we get weighted boundedness of the maximal operator $\mathcal{P}^*$ of truncations of the singular integral. We show that dyadic $A_p$ weights are the good weights for the maximal operator $\mathcal{P}^*$ of the scale truncations of $\mathcal{P}$.


title: MOLTE: a Modular Optimal Learning Testing Environment

abstract: We address the relative paucity of empirical testing of learning algorithms (of any type) by introducing a new public-domain, Modular, Optimal Learning Testing Environment (MOLTE) for Bayesian ranking and selection problem, stochastic bandits or sequential experimental design problems. The Matlab-based simulator allows the comparison of a number of learning policies (represented as a series of .m modules) in the context of a wide range of problems (each represented in its own .m module) which makes it easy to add new algorithms and new test problems. State-of-the-art policies and various problem classes are provided in the package. The choice of problems and policies is guided through a spreadsheet-based interface. Different graphical metrics are included. MOLTE is designed to be compatible with parallel computing to scale up from local desktop to clusters and clouds. We offer MOLTE as an easy-to-use tool for the research community that will make it possible to perform much more comprehensive testing, spanning a broader selection of algorithms and test problems. We demonstrate the capabilities of MOLTE through a series of comparisons of policies on a starter library of test problems. We also address the problem of tuning and constructing priors that have been largely overlooked in optimal learning literature. We envision MOLTE as a modest spur to provide researchers an easy environment to study interesting questions involved in optimal learning.


title: Cross-section geometry effects in the subband structure and spin-related properties of a HgTe/CdTe nanowire

abstract: By means of a multiband effective mass Hamiltonian, a theoretical characterization of the effect of the geometrical features of the confinement profile --in particular, a longitudinal groove-- on the subbands dispersion and spin-related properties of a rectangular HgTe/CdTe nanowire is presented. Through an external electric field applied perpendicular to the wire, the interplay of the induced Rashba spin splitting and these geometrical features is investigated. It is found that by exploiting this interplay a rich complexity of the subband structure arises, permitting the generation and modulation of spin-polarized currents without magnetic fields.


title: Predicting Organic Reaction Outcomes with Weisfeiler-Lehman Network

abstract: The prediction of organic reaction outcomes is a fundamental problem in computational chemistry. Since a reaction may involve hundreds of atoms, fully exploring the space of possible transformations is intractable. The current solution utilizes reaction templates to limit the space, but it suffers from coverage and efficiency issues. In this paper, we propose a template-free approach to efficiently explore the space of product molecules by first pinpointing the reaction center -- the set of nodes and edges where graph edits occur. Since only a small number of atoms contribute to reaction center, we can directly enumerate candidate products. The generated candidates are scored by a Weisfeiler-Lehman Difference Network that models high-order interactions between changes occurring at nodes across the molecule. Our framework outperforms the top-performing template-based approach with a 10\% margin, while running orders of magnitude faster. Finally, we demonstrate that the model accuracy rivals the performance of domain experts.


title: AMORPH: A statistical program for characterizing amorphous materials by X-ray diffraction

abstract: AMORPH utilizes a new Bayesian statistical approach to interpreting X-ray diffraction results of samples with both crystalline and amorphous components. AMORPH fits X-ray diffraction patterns with a mixture of narrow and wide components, simultaneously inferring all of the model parameters and quantifying their uncertainties. The program simulates background patterns previously applied manually, providing reproducible results, and significantly reducing inter- and intra-user biases. This approach allows for the quantification of amorphous and crystalline materials and for the characterization of the amorphous component, including properties such as the centre of mass, width, skewness, and nongaussianity of the amorphous component. Results demonstrate the applicability of this program for calculating amorphous contents of volcanic materials and independently modeling their properties in compositionally variable materials.


title: Contact Geometry and Quantum Mechanics

abstract: We present a generally covariant approach to quantum mechanics in which generalized positions, momenta and time variables are treated as coordinates on a fundamental "phase-spacetime." We show that this covariant starting point makes quantization into a purely geometric flatness condition. This makes quantum mechanics purely geometric, and possibly even topological. Our approach is especially useful for time-dependent problems and systems subject to ambiguities in choices of clock or observer. As a byproduct, we give a derivation and generalization of the Wigner functions of standard quantum mechanics.

categories: hep-th gr-qc math-ph math.DG math.MP quant-ph


title: Using NLU in Context for Question Answering: Improving on Facebook's bAbI Tasks

abstract: For the next step in human to machine interaction, Artificial Intelligence (AI) should interact predominantly using natural language because, if it worked, it would be the fastest way to communicate. Facebook's toy tasks (bAbI) provide a useful benchmark to compare implementations for conversational AI. While the published experiments so far have been based on exploiting the distributional hypothesis with machine learning, our model exploits natural language understanding (NLU) with the decomposition of language based on Role and Reference Grammar (RRG) and the brain-based Patom theory. Our combinatorial system for conversational AI based on linguistics has many advantages: passing bAbI task tests without parsing or statistics while increasing scalability. Our model validates both the training and test data to find 'garbage' input and output (GIGO). It is not rules-based, nor does it use parts of speech, but instead relies on meaning. While Deep Learning is difficult to debug and fix, every step in our model can be understood and changed like any non-statistical computer program. Deep Learning's lack of explicable reasoning has raised opposition to AI, partly due to fear of the unknown. To support the goals of AI, we propose extended tasks to use human-level statements with tense, aspect and voice, and embedded clauses with junctures: and answers to be natural language generation (NLG) instead of keywords. While machine learning permits invalid training data to produce incorrect test responses, our system cannot because the context tracking would need to be intentionally broken. We believe no existing learning systems can currently solve these extended natural language tests. There appears to be a knowledge gap between NLP researchers and linguists, but ongoing competitive results such as these promise to narrow that gap.


title: Schmid's Formula for Higher Local Fields

abstract: In local class field theory, the Schmid-Witt symbol encodes interesting data about the ramification theory of $p$-extensi-ons of $K$ and can, for example, be used to compute the higher ramification groups of such extensions. In 1936, Schmid discovered an explicit formula for the Schmid-Witt symbol of Artin-Schreier extensions of local fields. Later, his formula was generalized to Artin-Schreier-Witt extensions, but still over a local field. In this paper we generalize Schmid's formula to compute the Artin-Schreier-Witt-Parshin symbol for Artin-Schreier-Witt extensions of two-dimensional local fields of positive characteristic.


title: Herbig Ae Young Star's Infrared Spectrum Identified By Hydrocarbon Pentagon-Hexagon Combined Molecules

abstract: Infrared spectrum (IR) of Herbig Ae young stars was reproduced and classified by hydrocarbon pentagon-hexagon combined molecules by the quantum chemical calculation. Observed IR list by B. Acke et al. was categorized to four classes. Among 53 Herbig Ae stars, 26 samples show featured IR pattern named Type-D, which shows common IR peaks at 6.2, 8.3, 9.2, 10.0, 11.3, 12.1, and 14.0 micrometer. Typical star is HD144432. Calculation on di-cation molecule (C12H8)2+ having hydrocarbon one pentagon and two hexagons shows best coincidence at 6.1, 8.2, 9.2, 9.9, 11.3, 12.2, and 14.1 micrometer. There are some variation in Type-D. Spectrum of HD37357 was explained by a mixture with di-cation (C12H8)2+ and tri-cation (C12H8)3+. Ubiquitously observed spectrum Type-B was observed in 12 samples of Acke's list. In case of HD85567, observed 16 peaks were precisely reproduced by a single molecule (C23H12)2+. There is a mixture case with Type-B and Type-D. Typical example was HD142527. In this study, we could identify hidden carrier molecules for all types of IR in Herbig Ae stars.


title: First Weak-lensing Results from "See Change": Quantifying Dark Matter in the Two Z>1.5 High-redshift Galaxy Clusters SPT-CL J2040-4451 and IDCS J1426+3508

abstract: We present a weak-lensing study of SPT-CLJ2040-4451 and IDCSJ1426+3508 at z=1.48 and 1.75, respectively. The two clusters were observed in our "See Change" program, a HST survey of 12 massive high-redshift clusters aimed at high-z supernova measurements and weak-lensing estimation of accurate cluster masses. We detect weak but significant galaxy shape distortions using IR images from the WFC3, which has not yet been used for weak-lensing studies. Both clusters appear to possess relaxed morphology in projected mass distribution, and their mass centroids agree nicely with those defined by both the galaxy luminosity and X-ray emission. Using an NFW profile, for which we assume that the mass is tightly correlated with the concentration parameter, we determine the masses of SPT-CL J2040-4451 and IDCS J1426+3508 to be M_{200}=8.6_{-1.4}^{+1.7}x10^14 M_sun and 2.2_{-0.7}^{+1.1}x10^14 M_sun, respectively. The weak-lensing mass of SPT-CLJ2040-4451 shows that the cluster is clearly a rare object. Adopting the central value, the expected abundance of such a massive cluster at z>1.48 is only ~0.07 in the parent 2500 sq. deg. survey. However, it is yet premature to claim that the presence of this cluster creates a serious tension with the current LCDM paradigm unless that tension will remain in future studies after marginalizing over many sources of uncertainties such as the accuracy of the mass function and the mass-concentration relation at the high mass end. The mass of IDCSJ1426+3508 is in excellent agreement with our previous ACS-based weak-lensing result while the much higher source density from our WFC3 imaging data makes the current statistical uncertainty ~40% smaller.


title: An efficient adaptive sparse grid collocation method through derivative estimation

abstract: For uncertainty propagation of highly complex and/or nonlinear problems, one must resort to sample-based non-intrusive approaches [1]. In such cases, minimizing the number of function evaluations required to evaluate the response surface is of paramount importance. Sparse grid approaches have proven effective in reducing the number of sample evaluations. For example, the discrete projection collocation method has the notable feature of exhibiting fast convergence rates when approximating smooth functions; however, it lacks the ability to accurately and efficiently track response functions that exhibit fluctuations, abrupt changes or discontinuities in very localized regions of the input domain. On the other hand, the piecewise linear collocation interpolation approach can track these localized variations in the response surface efficiently, but it converges slowly in the smooth regions. The proposed methodology, building on an existing work on adaptive hierarchical sparse grid collocation algorithm [2], is able to track localized behavior while also avoiding unnecessary function evaluations in smoother regions of the stochastic space by using a finite difference based one-dimensional derivative evaluation technique in all the dimensions. This derivative evaluation technique leads to faster convergence in the smoother regions than what is achieved in the existing collocation interpolation approaches. Illustrative examples show that this method is well suited to high-dimensional stochastic problems, and that stochastic elliptic problems with stochastic dimension as high as 100 can be dealt with effectively.


title: Two-Step Adsorption of $PtCl_6^{2-}$ Complexes at a Charged Langmuir Monolayer: Role of Hydration and Ion Correlations

abstract: Anion exchange at positively charged interfaces plays an important role in a variety of physical and chemical processes. However, the molecular scale details of these processes, especially with heavy and large anionic complexes, are not well-understood. We studied the adsorption of $PtCl_6^{2-}$ anionic complexes to floating DPTAP monolayers in the presence of excess $Cl^-$ as a function of the bulk chlorometalate concentration. In situ x-ray scattering and fluorescence measurements, which are element and depth sensitive, show that the chlorometalate ions only adsorb in the diffuse layer at lower concentrations, while they adsorb predominantly in the Stern layer at higher concentrations. The response of DPTAP molecules to the adsorbed ions is determined independently by grazing incidence x-ray diffraction, and supports this picture. Molecular dynamics simulations further elucidate the nanoscale structure of the interfacial complexes. The results suggest that ion hydration and ion-ion correlations play a key role in the competitive adsorption process.


title: Pairwise versus multiple network alignment

abstract: Biological network alignment (NA) aims to identify similar regions between molecular networks of different species. NA can be local or global. Just as the recent trend in the NA field, we also focus on global NA, which can be pairwise (PNA) and multiple (MNA). PNA produces aligned node pairs between two networks. MNA produces aligned node clusters between more than two networks. Recently, the focus has shifted from PNA to MNA, because MNA captures conserved regions between more networks than PNA (and MNA is thus considered to be more insightful), though at higher computational complexity. The issue is that, due to the different outputs of PNA and MNA, a PNA method is only compared to other PNA methods, and an MNA method is only compared to other MNA methods. Comparison of PNA against MNA must be done to evaluate whether MNA's higher complexity is justified by its higher accuracy. We introduce a framework that allows for this. We evaluate eight prominent PNA and MNA methods, on synthetic and real-world biological networks, using topological and functional alignment quality measures. We compare PNA against MNA in both a pairwise (native to PNA) and multiple (native to MNA) manner. PNA is expected to perform better under the pairwise evaluation framework. Indeed this is what we find. MNA is expected to perform better under the multiple evaluation framework. Shockingly, we find this not to always hold; PNA is often better than MNA in this framework, depending on the choice of evaluation test.


title: Self-Adaptive Tensor Network States with Multi-Site Correlators

abstract: We introduce the concept of self-adaptive tensor network states (SATNS) based on multi-site correlators. The SATNS ansatz gradually extends its variational space incorporating the most important next-order correlators into the ansatz for the wave function. The selection of these correlators is guided by entanglement-entropy measures from quantum information theory. By sequentially introducing variational parameters and adjusting them to the system under study, the SATNS ansatz achieves to keep their number significantly smaller than the total number of full-configuration interaction parameters. The SATNS ansatz is studied for manganocene in its lowest-energy sextet and doublet states, the latter of which is known to be difficult to describe. It is shown that the SATNS parametrization solves the convergence issues found for previous correlator-based tensor network states.


title: A group action on multivariate polynomials over finite fields

abstract: Let $\mathbb{F}_q$ be the finite field with $q$ elements, where $q$ is a power of a prime $p$. Recently, a particular action of the group $\mathrm{GL}_2(\mathbb F_q)$ on irreducible polynomials in $\mathbb F_q[x]$ has been introduced and many questions concerning the invariant polynomials have been discussed. In this paper, we give a natural extension of this action on the polynomial ring $\mathbb F_q[x_1, \ldots, x_n]$ and study the algebraic properties of the invariant elements.


title: Definable Combinatorics of Some Borel Equivalence Relations

abstract: If $X$ is a set, $E$ is an equivalence relation on $X$, and $n \in \omega$, then define $$[X]^n_E = \{(x_0, ..., x_{n - 1}) \in {}^nX : (\forall i,j)(i \neq j \Rightarrow \neg(x_i \ E \ x_j))\}.$$ For $n \in \omega$, a set $X$ has the $n$-J\'onsson property if and only if for every function $f : [X]^n_= \rightarrow X$, there exists some $Y \subseteq X$ with $X$ and $Y$ in bijection so that $f[[Y]^n_=] \neq X$. A set $X$ has the J\'onsson property if and only for every function $f : (\bigcup_{n \in \omega}[X]^n_=) \rightarrow X$, there exists some $Y \subseteq X$ with $X$ and $Y$ in bijection so that $f[\bigcup_{n \in \omega} [Y]^n_=] \neq X$. Let $n \in \omega$, $X$ be a Polish space, and $E$ be an equivalence relation on $X$. $E$ has the $n$-Mycielski property if and only if for all comeager $C \subseteq {}^nX$, there is some $\mathbf{\Delta_1^1}$ $A \subseteq X$ so that $E \leq_{\mathbf{\Delta_1^1}} E \upharpoonright A$ and $[A]^n_E \subseteq C$. The following equivalence relations will be considered: $E_0$ is defined on ${}^\omega2$ by $x \ E_0 \ y$ if and only if $(\exists n)(\forall k > n)(x(k) = y(k))$. $E_1$ is defined on ${}^\omega({}^\omega2)$ by $x \ E_1 \ y$ if and only if $(\exists n)(\forall k > n)(x(k) = y(k))$. $E_2$ is defined on ${}^\omega2$ by $x \ E_2 \ y$ if and only if $\sum\{\frac{1}{n + 1} : n \in x \ \triangle \ y\} < \infty$, where $\triangle$ denotes the symmetric difference. $E_3$ is defined on ${}^\omega({}^\omega2)$ by $x \ E_3 \ y$ if and only if $(\forall n)(x(n) \ E_0 \ y(n))$. Holshouser and Jackson have shown that $\mathbb{R}$ is J\'onsson under $\mathsf{AD}$. It will be shown that $E_0$ does not have the $3$-Mycielski property and that $E_1$, $E_2$, and $E_3$ do not have the $2$-Mycielski property. Under $\mathsf{ZF + AD}$, ${}^\omega 2 / E_0$ does not have the $3$-J\'onsson property.


title: Structural properties of edge-chromatic critical multigraphs

abstract: Appearing in different format, Gupta\,(1967), Goldberg\,(1973), Andersen\,(1977), and Seymour\,(1979) conjectured that if $G$ is an edge-$k$-critical graph with $k \ge \Delta +1$, then $|V(G)|$ is odd and, for every edge $e$, $E(G-e)$ is a union of disjoint near-perfect matchings, where $\Delta$ denotes the maximum degree of $G$. Tashkinov tree method shows that critical graphs contain a subgraph with two important properties named closed and elementary. Recently, efforts have been made in extending graphs beyond Tashkinov trees. However, these results can only keep one of the two essential properties. In this paper, we developed techniques to extend Tashkinov trees to larger subgraphs with both properties. Applying our result, we have improved almost all known results towards Goldberg's conjecture. In particular, we showed that Goldberg's conjecture holds for graph $G$ with $|V(G)| \le 39$ and $|\Delta(G)| \le 39$ and Jacobsen's equivalent conjecture holds for $m \le 39$ while the previous known bound is $23$.


title: REMOTEGATE: Incentive-Compatible Remote Configuration of Security Gateways

abstract: Imagine that a malicious hacker is trying to attack a server over the Internet and the server wants to block the attack packets as close to their point of origin as possible. However, the security gateway ahead of the source of attack is untrusted. How can the server block the attack packets through this gateway? In this paper, we introduce REMOTEGATE, a trustworthy mechanism for allowing any party (server) on the Internet to configure a security gateway owned by a second party, at a certain agreed upon reward that the former pays to the latter for its service. We take an interactive incentive-compatible approach, for the case when both the server and the gateway are rational, to devise a protocol that will allow the server to help the security gateway generate and deploy a policy rule that filters the attack packets before they reach the server. The server will reward the gateway only when the latter can successfully verify that it has generated and deployed the correct rule for the issue. This mechanism will enable an Internet-scale approach to improving security and privacy, backed by digital payment incentives.


title: Learning Unknown Markov Decision Processes: A Thompson Sampling Approach

abstract: We consider the problem of learning an unknown Markov Decision Process (MDP) that is weakly communicating in the infinite horizon setting. We propose a Thompson Sampling-based reinforcement learning algorithm with dynamic episodes (TSDE). At the beginning of each episode, the algorithm generates a sample from the posterior distribution over the unknown model parameters. It then follows the optimal stationary policy for the sampled model for the rest of the episode. The duration of each episode is dynamically determined by two stopping criteria. The first stopping criterion controls the growth rate of episode length. The second stopping criterion happens when the number of visits to any state-action pair is doubled. We establish $\tilde O(HS\sqrt{AT})$ bounds on expected regret under a Bayesian setting, where $S$ and $A$ are the sizes of the state and action spaces, $T$ is time, and $H$ is the bound of the span. This regret bound matches the best available bound for weakly communicating MDPs. Numerical results show it to perform better than existing algorithms for infinite horizon MDPs.


title: When Waiting is not an Option : Learning Options with a Deliberation Cost

abstract: Recent work has shown that temporally extended actions (options) can be learned fully end-to-end as opposed to being specified in advance. While the problem of "how" to learn options is increasingly well understood, the question of "what" good options should be has remained elusive. We formulate our answer to what "good" options should be in the bounded rationality framework (Simon, 1957) through the notion of deliberation cost. We then derive practical gradient-based learning algorithms to implement this objective. Our results in the Arcade Learning Environment (ALE) show increased performance and interpretability.


title: On the Korn interpolation and second inequalities in thin domains

abstract: We consider shells of non-constant thickness in three dimensional Euclidean space around surfaces which have bounded principal curvatures. We derive Korn's interpolation (or the so called first and a half (The inequality first introduced in [Gra.Har.1])) and second inequalities on that kind of domains for $\Bu\in H^1$ vector fields, imposing no boundary or normalization conditions on $\Bu.$ The constants in the estimates are asymptotically optimal in terms of the domain thickness $h,$ with the leading order constant having the scaling $h$ as $h\to 0.$ This is the first work that determines the asymptotics of the optimal constant in the classical Korn second inequality for shells in terms of the domain thickness in almost full generality, the inequality being fulfilled for practically all thin domains $\Omega\in\mathbb R^3$ and all vector fields $\Bu\in H^1(\Omega).$ Moreover, the Korn interpolation inequality is stronger than Korn's second inequality, and it reduces the problem of estimating the gradient $\nabla\Bu$ in terms of the symmetrized gradient $e(\Bu)$, in particular any linear geometric rigidity estimates for thin domains, to the easier problem of proving the corresponding Poincar\'e-like estimates on the field $\Bu$ itself.


title: Bloch oscillations and subwavelength focusing in stacked fishnet metamaterials

abstract: In this letter, we introduce stacked fishnet metamaterial for steering light in microwave region. We numerically demonstrate that optical Bloch oscillations and a focus of as small as one sixth of a wavelength can be achieved. The flexibility of varying geometrical parameters of the fishnet slabs provides an efficient way for tuning its local effective media parameters and opens the possibility for controlling light arbitrarily. The experiment verifies subwavelength-sized light focusing effect by scanning magnetic field at the surface of the sample directly.


title: Towards personalized human AI interaction - adapting the behavior of AI agents using neural signatures of subjective interest

abstract: Reinforcement Learning AI commonly uses reward/penalty signals that are objective and explicit in an environment -- e.g. game score, completion time, etc. -- in order to learn the optimal strategy for task performance. However, Human-AI interaction for such AI agents should include additional reinforcement that is implicit and subjective -- e.g. human preferences for certain AI behavior -- in order to adapt the AI behavior to idiosyncratic human preferences. Such adaptations would mirror naturally occurring processes that increase trust and comfort during social interactions. Here, we show how a hybrid brain-computer-interface (hBCI), which detects an individual's level of interest in objects/events in a virtual environment, can be used to adapt the behavior of a Deep Reinforcement Learning AI agent that is controlling a virtual autonomous vehicle. Specifically, we show that the AI learns a driving strategy that maintains a safe distance from a lead vehicle, and most novelly, preferentially slows the vehicle when the human passengers of the vehicle encounter objects of interest. This adaptation affords an additional 20\% viewing time for subjectively interesting objects. This is the first demonstration of how an hBCI can be used to provide implicit reinforcement to an AI agent in a way that incorporates user preferences into the control system.


title: Exfoliation energy, quasi-particle bandstructure, and excitonic properties of selenium and tellurium atomic chains

abstract: Effects that are not captured by the generalized-gradient density-functional theory play a prominent effect on the structural binding, and on the electronic and optical properties of reduced-dimensional and weakly-bound materials. Here, we report the exfoliation energy of selenium and tellurium atomic chains with non-empirical van der Waals corrections, and their electronic and optical properties with the GW and Bethe-Salpeter formalisms. The exfoliation energy is found to be within 0.547 to 0.719 eV/u.c. for the selenium atomic chain, and 0.737 to 0.926 eV/u.c. for the tellurium atomic chain (u.c. stands for unit cell), depending on the approximation for the van der Waals interaction and the numerical tool chosen. The GW electronic bandgap turned out to be 5.22--5.47 (4.44--4.59) eV for the Se (Te) atomic chains, with the lowest bound obtained with the Godby-Needs (GB), and the upper bound to the Hybertsen-Louie (HL) plasmon-pole models (PPMs). The binding energy of the ground-state excitonic state ranges between 2.69 to 2.72 eV for selenium chains within the HL and GB PPM, respectively, and turned out to be 2.35 eV for tellurium chains with both approximations. The ground state excitonic wave function is localized within 50 \AA{} along the axis for both types of atomic chains, and its energy lies within the visible spectrum: blue [2.50(GN)--2.78(HL) eV] for selenium, and yellow--green [2.09(GN)--2.28(HL) eV] for tellurium, which could be useful for LED applications in the visible spectrum.


title: Catalyst design using actively learned machine with non-ab initio input features towards CO2 reduction reactions

abstract: In conventional chemisorption model, the d-band center theory (augmented sometimes with the upper edge of d-band for imporved accuarcy) plays a central role in predicting adsorption energies and catalytic activity as a function of d-band center of the solid surfaces, but it requires density functional calculations that can be quite costly for large scale screening purposes of materials. In this work, we propose to use the d-band width of the muffin-tin orbital theory (to account for local coordination environment) plus electronegativity (to account for adsorbate renormalization) as a simple set of alternative descriptors for chemisorption, which do not demand the ab initio calculations. This pair of descriptors are then combined with machine learning methods, namely, artificial neural network (ANN) and kernel ridge regression (KRR), to allow large scale materials screenings. We show, for a toy set of 263 alloy systems, that the CO adsorption energy can be predicted with a remarkably small mean absolute deviation error of 0.05 eV, a significantly improved result as compared to 0.13 eV obtained with descriptors including costly d-band center calculations in literature. We achieved this high accuracy by utilizing an active learning algorithm, without which the accuracy was 0.18 eV otherwise. As a practical application of this machine, we identified Cu3Y@Cu as a highly active and cost-effective electrochemical CO2 reduction catalyst to produce CO with the overpotential 0.37 V lower than Au catalyst.


title: DeepVoting: A Robust and Explainable Deep Network for Semantic Part Detection under Partial Occlusion

abstract: In this paper, we study the task of detecting semantic parts of an object, e.g., a wheel of a car, under partial occlusion. We propose that all models should be trained without seeing occlusions while being able to transfer the learned knowledge to deal with occlusions. This setting alleviates the difficulty in collecting an exponentially large dataset to cover occlusion patterns and is more essential. In this scenario, the proposal-based deep networks, like RCNN-series, often produce unsatisfactory results, because both the proposal extraction and classification stages may be confused by the irrelevant occluders. To address this, [25] proposed a voting mechanism that combines multiple local visual cues to detect semantic parts. The semantic parts can still be detected even though some visual cues are missing due to occlusions. However, this method is manually-designed, thus is hard to be optimized in an end-to-end manner. In this paper, we present DeepVoting, which incorporates the robustness shown by [25] into a deep network, so that the whole pipeline can be jointly optimized. Specifically, it adds two layers after the intermediate features of a deep network, e.g., the pool-4 layer of VGGNet. The first layer extracts the evidence of local visual cues, and the second layer performs a voting mechanism by utilizing the spatial relationship between visual cues and semantic parts. We also propose an improved version DeepVoting+ by learning visual cues from context outside objects. In experiments, DeepVoting achieves significantly better performance than several baseline methods, including Faster-RCNN, for semantic part detection under occlusion. In addition, DeepVoting enjoys explainability as the detection results can be diagnosed via looking up the voting cues.


title: Rota-Baxter modules toward derived functors

abstract: In this paper we study Rota-Baxter modules with emphasis on the role played by the Rota-Baxter operators and resulting difference between Rota-Baxter modules and the usual modules over an algebra. We introduce the concepts of free, projective, injective and flat Rota-Baxter modules. We give the construction of free modules and show that there are enough projective, injective and flat Rota-Baxter modules to provide the corresponding resolutions for derived functor.


title: Autonomous Extracting a Hierarchical Structure of Tasks in Reinforcement Learning and Multi-task Reinforcement Learning

abstract: Reinforcement learning (RL), while often powerful, can suffer from slow learning speeds, particularly in high dimensional spaces. The autonomous decomposition of tasks and use of hierarchical methods hold the potential to significantly speed up learning in such domains. This paper proposes a novel practical method that can autonomously decompose tasks, by leveraging association rule mining, which discovers hidden relationship among entities in data mining. We introduce a novel method called ARM-HSTRL (Association Rule Mining to extract Hierarchical Structure of Tasks in Reinforcement Learning). It extracts temporal and structural relationships of sub-goals in RL, and multi-task RL. In particular,it finds sub-goals and relationship among them. It is shown the significant efficiency and performance of the proposed method in two main topics of RL.


title: Topological end states and Zak phase of rectangular armchair ribbon

abstract: We consider the end states of a half-filled rectangular armchair graphene ribbon (RAGR) in a staggered potential. Taking electron-electron interactions into account we find that, as the strength of the staggered potential varies, three types of couplings between the end states can occur: antiferromagnetic without or with spin splitting, and paramagnetic without spin-splitting. We find that a spin-splitting is present only in the staggered potential region $0<\Delta<\Delta_c$. The transition from the antiferromagnetic state at $\Delta=0$ to the paramagnetic state goes through an intermediate spin-split antiferromagnetic state, and this spin-splitting disappears suddenly at $\Delta_c$. For small and large values of $\Delta$ the end charge of a RAGR can be connected to the Zak phase of the periodic armchair graphene ribbon (PARG) with the same width, and it varies continuously as the strength of the potential changes.


title: Tuning the ferro- to para-electric transition temperature and dipole orientation of group-IV monochalcogenide monolayers

abstract: Coordination-related, two-dimensional (2D) structural phase transitions are a fascinating and novel facet of two-dimensional materials with structural degeneracies. Nevertheless, a unified theoretical account of these transitions remains absent, and the following points are established through {\em ab-initio} molecular dynamics and 2D discrete clock models here: Group-IV monochalcogenide (GeSe, SnSe, SnTe, ...) monolayers have four degenerate structural ground states, and a 2D phase transition from a three-fold coordinated onto a five-fold coordinated structure takes place at finite temperature. On unstrained samples, the 2D phase transition requires lattice parameters to freely evolve. A fundamental energy scale permits understanding this transition. The transition temperature $T_c$ and the orientation of the in-plane intrinsic electric dipole can be controlled by moderate uniaxial tensile strain, and a modified discrete clock model describes the transition on strained samples. These results establish a general underlying theoretical background to understand structural phase transitions in 2D materials and their effects on material properties.


title: Pisano period codes

abstract: The cyclic codes with parity check polynomial the reciprocal of the characteristic polynomial of the Fibonacci recurrence over a prime finite field are shown to have either one weight or two weights. When these codes are irreducible cyclic we obtain many counterexamples to the conjectural classification of two-weight irreducible cyclic codes of Schmidt and White (2002). When they are reducible and projective their duals are uniformly packed.


title: Acceleration of Histogram-Based Contrast Enhancement via Selective Downsampling

abstract: In this paper, we propose a general framework to accelerate the universal histogram-based image contrast enhancement (CE) algorithms. Both spatial and gray-level selective down-sampling of digital images are adopted to decrease computational cost, while the visual quality of enhanced images is still preserved and without apparent degradation. Mapping function calibration is novelly proposed to reconstruct the pixel mapping on the gray levels missed by downsampling. As two case studies, accelerations of histogram equalization (HE) and the state-of-the-art global CE algorithm, i.e., spatial mutual information and PageRank (SMIRANK), are presented detailedly. Both quantitative and qualitative assessment results have verified the effectiveness of our proposed CE acceleration framework. In typical tests, computational efficiencies of HE and SMIRANK have been speeded up by about 3.9 and 13.5 times, respectively.


title: Stochastic collocation approach with adaptive mesh refinement for parametric uncertainty analysis

abstract: Presence of a high-dimensional stochastic parameter space with discontinuities poses major computational challenges in analyzing and quantifying the effects of the uncertainties in a physical system. In this paper, we propose a stochastic collocation method with adaptive mesh refinement (SCAMR) to deal with high dimensional stochastic systems with discontinuities. Specifically, the proposed approach uses generalized polynomial chaos (gPC) expansion with Legendre polynomial basis and solves for the gPC coefficients using the least squares method. It also implements an adaptive mesh (element) refinement strategy which checks for abrupt variations in the output based on the second order gPC approximation error to track discontinuities or non-smoothness. In addition, the proposed method involves a criterion for checking possible dimensionality reduction and consequently, the decomposition of the full-dimensional problem to a number of lower-dimensional subproblems. Specifically, this criterion checks all the existing interactions between input dimensions of a specific problem based on the high-dimensional model representation (HDMR) method, and therefore automatically provides the subproblems which only involve interacting dimensions. The efficiency of the approach is demonstrated using both smooth and non-smooth function examples with input dimensions up to 300, and the approach is compared against other existing algorithms.


title: Two-weight codes and second order recurrences

abstract: Cyclic codes of dimension $2$ over a finite field are shown to have at most two nonzero weights. This extends a construction of Rao et al (2010) and disproves a conjecture of Schmidt-White (2002). We compute their weight distribution, and give a condition on the roots of their check polynomials for them to be MDS.


title: CO2 packing polymorphism under confinement in cylindrical nanopores

abstract: We investigate the effect of cylindrical nano-confinement on the phase behaviour of a rigid model of carbon dioxide using both molecular dynamics and well tempered metadynamics. To this aim we study a simplified pore model across a parameter space comprising pore diameter, CO2-pore wall potential and CO2 density. In order to systematically identify ordering events within the pore model we devise a generally applicable approach based on the analysis of the distribution of intermolecular orientations. Our simulations suggest that, while confinement in nano-pores inhibits the formation of known crystal structures, it induces a remarkable variety of ordered packings unrelated to their bulk counterparts, and favours the establishment of short range order in the fluid phase. We summarise our findings by proposing a qualitative phase diagram for this model.


title: The origin of regular Newtonian potential in infinite derivative gravity

abstract: It turns out that the infinite derivative gravity (IDG) is ghost-free and renormalizable when one chooses the exponential of an entire function. For this IDG case, the corresponding Newtonian potential generated from the delta function is non-singular at the origin. However, we will explicitly show that the source generating this non-singular potential is given not by the delta-function due to the point-like source of mass, but by the Gaussian mass distribution. This explains clearly why the IDG with the exponential of an entire function yields the finite potential at the origin.


title: Bayesian support for Evolution: detecting phylogenetic signal in a subset of the primate family

abstract: The theory of evolution states that the diversity of species can be explained by descent with modification. Therefore, all living beings are related through a common ancestor. This evolutionary process must have left traces in our molecular composition. In this work, we present a randomization procedure in order to determine if a group of 5 species of the primate family, namely, macaque, guereza, orangutan, chimpanzee and human, has retained these traces in its molecules. Firstly, we present the randomization methodology through two toy examples, which allow to understand its logic. We then carry out a DNA data analysis to assess if the group of primates contains phylogenetic information which links them in a joint evolutionary history. This is carried out by monitoring a Bayesian measure, called marginal likelihood, which we estimate by using nested sampling. We found that it would be unusual to get the relationship observed in the data among these primate species if they had not shared a common ancestor. The results are in total agreement with the theory of evolution.


title: Generalizing causal inferences from individuals in randomized trials to all trial-eligible individuals

abstract: We consider methods for causal inference in randomized trials nested within cohorts of trial-eligible individuals, including those who are not randomized. We show how baseline covariate data from the entire cohort, and treatment and outcome data only from randomized individuals, can be used to identify potential (counterfactual) outcome means and average treatment effects in the target population of all eligible individuals. We review identifiability conditions, propose estimators, and assess the estimators' finite-sample performance in simulation studies. As an illustration, we apply the estimators in a trial nested within a cohort of trial-eligible individuals to compare coronary artery bypass grafting surgery plus medical therapy vs. medical therapy alone for chronic coronary artery disease.


title: Simultaneous drop in mean free path and carrier density at the pseudogap onset in high-$T_c$ cuprates

abstract: High-temperature superconducting cuprates are distinguished by an enigmatic pseudogap which opens near optimal doping where the superconducting transition temperature is highest. Key questions concern its origin and whether it is essential in any way to superconductivity. Recent field-induced normal-state transport experiments on hole-doped cuprates have measured abrupt changes in the doping dependent Hall number and resistivity, consistent with a drop in carrier density from $1+p$ to $p$ holes per copper atom, on entering the pseudogap phase. In this work the change in resistivity is analyzed in terms of an antiferromagnetic-order-induced Fermi surface reconstruction model that has already successfully described the Hall number. In order for this model to describe the resistivity we find that the zero-temperature mean free path must also drop abruptly in proportion to the size of the Fermi surface. This suggests that intrapocket scattering underlies the observed upturn in resistivity in the pseudogap state.