title: Fock majorization in bosonic quantum channels with a passive environment

abstract: We introduce a class of quantum channels called passive-environment bosonic channels. These channels are relevant from a quantum thermodynamical viewpoint because they correspond to the energy-preserving linear coupling of a bosonic system with a bosonic environment that is in a passive state (no energy can be extracted from it by using a unitary transformation) followed by discarding the environment. The Fock-majorization relation defined in [New J. Phys. 18, 073047 (2016)] happens to be especially useful in this context as, unlike regular majorization, it connects the disorder of a state together with its energy. Our main result here is the preservation of Fock majorization across all passive-environment bosonic channels. This implies a similar preservation property for regular majorization over the set of passive states, and it also extends to passive-environment channels whose Stinespring dilation involves an active Gaussian unitary. Beyond bosonic systems, the introduced class of passive environment operations naturally generalizes thermal operations and is expected to provide new insights into the thermodynamics of quantum systems.


title: How to optimize preconditioners for the conjugate gradient method: a stochastic approach

abstract: The conjugate gradient method (CG) is typically used with a preconditioner which improves efficiency and robustness of the method. Many preconditioners include parameters and a proper choice of a preconditioner and its parameters is often not a trivial task. Although many convergence estimates exist which can be used for optimizing preconditioners, these estimates typically hold for all initial guess vectors, in other words, they reflect the worst convergence rate. To account for the mean convergence rate instead, in this paper, we follow a stochastic approach. It is based on trial runs with random initial guess vectors and leads to a functional which can be used to monitor convergence and to optimize preconditioner parameters in CG. Presented numerical experiments show that optimization of this new functional with respect to preconditioner parameters usually yields a better parameter value than optimization of the functional based on the spectral condition number.


title: On the coexistence of dipolar frustration and criticality in ferromagnets

abstract: In real magnets the tendency towards ferromagnetism, promoted by exchange coupling, is usually frustrated by dipolar interaction. As a result, the uniformly ordered phase is replaced by modulated (multi-domain) phases, characterized by different order parameters rather than the global magnetization. The transitions occurring within those modulated phases and towards the disordered phase are generally not of second-order type. Nevertheless, strong experimental evidence indicates that a standard critical behavior is recovered when comparatively small fields are applied that stabilize the uniform phase. The resulting power laws are observed with respect to a putative critical point that falls in the portion of the phase diagram occupied by modulated phases, in line with an avoided-criticality scenario. Here we propose a generalization of the scaling hypothesis for ferromagnets, which explains this observation assuming that the dipolar interaction acts as a relevant field, in the sense of renormalization group.


title: Computationally Efficient Estimation of the Spectral Gap of a Markov Chain

abstract: We consider the problem of estimating from sample paths the absolute spectral gap $\gamma_*$ of a reversible, irreducible and aperiodic Markov chain $(X_t)_{t \in \mathbb{N}}$ over a finite state space $\Omega$. We propose the ${\tt UCPI}$ (Upper Confidence Power Iteration) algorithm for this problem, a low-complexity algorithm which estimates the spectral gap in time ${\cal O}(n)$ and memory space ${\cal O}((\ln n)^2)$ given $n$ samples. This is in stark contrast with most known methods which require at least memory space ${\cal O}(|\Omega|)$, so that they cannot be applied to large state spaces. Furthermore, ${\tt UCPI}$ is amenable to parallel implementation.


title: Positive radial solutions for the Minkowski-curvature equation with Neumann boundary conditions

abstract: We analyze existence, multiplicity and oscillatory behavior of positive radial solutions to a class of quasilinear equations governed by the Lorentz-Minkowski mean curvature operator. The equation is set in a ball or an annulus of $\mathbb R^N$, is subject to homogeneous Neumann boundary conditions, and involves a nonlinear term on which we do not impose any growth condition at infinity. The main tool that we use is the shooting method for ODEs.


title: Planar Superconductor-Insulator-Superconductor Mixer Array Receivers for Wide Field of View Astronomical Observation

abstract: We present a conceptual framework of planar SIS mixer array receivers and the studies on the required techniques. This concept features membrane-based on-chip waveguide probes and a quasi-two-dimensional local-oscillator distribution waveguide network. This concept allows sophisticated functions, such as dual-polarization, balanced mixing and sideband separation, easily implemented with the SIS mixer array in the same planar circuit. We have developed a single-pixel prototype receiver by implementing the concept in the design. Initial measurement results show good evidences that support the feasibility of the concept.


title: Intensity mapping with SDSS/BOSS Lyman-alpha emission, quasars and their Lyman-alpha forest

abstract: We investigate the large-scale structure of Lyman-alpha emission intensity in the Universe at redshifts z=2-3.5 using cross-correlation techniques. Our Lya emission samples are spectra of BOSS Luminous Red Galaxies from Data Release 12 with the best fit model galaxies subtracted. We cross-correlate the residual flux in these spectra with BOSS quasars, and detect a positive signal on scales 1-15 Mpc/h. We identify and remove a source of contamination not previously accounted for, due to the effects of quasar clustering on cross-fibre light. Corrected, our quasar-Lya emission cross-correlation is 50 % lower than that seen by Croft et al. for DR10, but still significant. Because only 3% of space is within 15 Mpc/h of a quasar, the result does not fully explore the global large-scale structure of Lya emission. To do this, we cross-correlate with the Lya forest. We find no signal in this case. The 95% upper limit on the global Lya mean surface brightness from Lya emission-Lya forest cross-correlation is mu < 1.2x10^-22 erg/s/cm^2/A/arcsec^2 This null result rules out the scenario where the observed quasar-Lya emission cross-correlation is primarily due to the large scale structure of star forming galaxies, Taken in combination, our results suggest that Lya emitting galaxies contribute, but quasars dominate within 15 Mpc/h. A simple model for Lya emission from quasars based on hydrodynamic simulations reproduces both the observed forest-Lya emission and quasar-Lya emission signals. The latter is also consistent with extrapolation of observations of fluorescent emission from smaller scales r < 1 Mpc.


title: Positivity of hexagon perturbation theory

abstract: The hexagon-form-factor program was proposed as a way to compute three- and higher-point correlation functions in $\mathcal{N}=4$ super-symmetric Yang-Mills theory and in the dual AdS$_5\times$S$^5$ superstring theory, by exploiting the integrability of the theory in the 't Hooft limit. This approach is reminiscent of the asymptotic Bethe ansatz in that it applies to a large-volume expansion. Finite-volume corrections can be incorporated through L\"uscher-like formulae, though the systematics of this expansion is largely unexplored so far. Strikingly, finite-volume corrections may feature negative powers of the 't Hooft coupling $g$ in the small-$g$ expansion, potentially leading to a breakdown of the formalism. In this work we show that the finite-volume perturbation theory for the hexagon is positive and thereby compatible with the weak-coupling expansion for arbitrary $n$-point functions.


title: Intermediate-Mass Black Holes in Globular Cluster Systems

abstract: Using the Next Generation Very Large Array (ngVLA), we will make a comprehensive inventory of intermediate-mass black holes (IMBHs) in hundreds of globular cluster systems out to a distance of 25 Mpc. IMBHs have masses of about 100 to 100,000 solar masses. Finding them in globular clusters would validate a formation channel for seed black holes in the early universe and inform event predictions for gravitational wave facilities. Reaching a large number of globular clusters is key, as Fragione et al. (2018) predict that only a few percent will have retained their gravitational-wave fostering IMBHs.


title: Deep Lip Reading: a comparison of models and an online application

abstract: The goal of this paper is to develop state-of-the-art models for lip reading -- visual speech recognition. We develop three architectures and compare their accuracy and training times: (i) a recurrent model using LSTMs; (ii) a fully convolutional model; and (iii) the recently proposed transformer model. The recurrent and fully convolutional models are trained with a Connectionist Temporal Classification loss and use an explicit language model for decoding, the transformer is a sequence-to-sequence model. Our best performing model improves the state-of-the-art word error rate on the challenging BBC-Oxford Lip Reading Sentences 2 (LRS2) benchmark dataset by over 20 percent. As a further contribution we investigate the fully convolutional model when used for online (real time) lip reading of continuous speech, and show that it achieves high performance with low latency.


title: Fractal mesh refinement, rare events and type Ia supenova

abstract: Fractal mesh refinement enables mesh refinement regimes to the Gibson scale and beyond. Conventional multifractal subgrid models have cost considerations which also lead to local mesh refinement, but more importantly, the application of these models to compressible turbulent deflagration fronts is a topic for future research, while subgrid models tuned to this complex physics lack a multifractal search focus, and for cost reasons do not allow the large search volumes required to find the sought for rare event to trigger a DDT. Here we propose methods to resolve fine scales and locate rare possible DDT trigger events within large volumes while addressing multifractal issues at a feasible computational cost. We are motivated by the goal of confirming or assessing the hypothesis of a delayed detonation in type Ia supernova and to assessing the delay in this event if it is found to occur.


title: Classification with Fairness Constraints: A Meta-Algorithm with Provable Guarantees

abstract: Developing classification algorithms that are fair with respect to sensitive attributes of the data has become an important problem due to the growing deployment of classification algorithms in various social contexts. Several recent works have focused on fairness with respect to a specific metric, modeled the corresponding fair classification problem as a constrained optimization problem, and developed tailored algorithms to solve them. Despite this, there still remain important metrics for which we do not have fair classifiers and many of the aforementioned algorithms do not come with theoretical guarantees; perhaps because the resulting optimization problem is non-convex. The main contribution of this paper is a new meta-algorithm for classification that takes as input a large class of fairness constraints, with respect to multiple non-disjoint sensitive attributes, and which comes with provable guarantees. This is achieved by first developing a meta-algorithm for a large family of classification problems with convex constraints, and then showing that classification problems with general types of fairness constraints can be reduced to those in this family. We present empirical results that show that our algorithm can achieve near-perfect fairness with respect to various fairness metrics, and that the loss in accuracy due to the imposed fairness constraints is often small. Overall, this work unifies several prior works on fair classification, presents a practical algorithm with theoretical guarantees, and can handle fairness metrics that were previously not possible.


title: Axion Isocurvature Perturbations in Low-Scale Models of Hybrid Inflation

abstract: The QCD axion solves the strong CP problem and represents an attractive particle candidate for cold dark matter (CDM). However, quantum fluctuations of the axion field during inflation easily result in large CDM isocurvature perturbations that are in conflict with observations of the cosmic microwave background (CMB). In this paper, we demonstrate how this problem can be solved in low-scale models of hybrid inflation that may emerge from supersymmetric grand unified theories. We consider both F-term hybrid inflation (FHI) and D-term hybrid inflation (DHI) in supergravity, explicitly taking into account the effect of hidden-sector supersymmetry breaking. We discuss the production of cosmic strings and show how the soft terms in the scalar potential readily allow to achieve the correct scalar spectral index. In both cases, we are able to identify large regions in parameter space that are consistent with all constraints. In particular, we find that evading the CDM isocurvature constraint always requires a small Yukawa or gauge coupling of O(10^-3) or smaller. This translates into upper bounds on the gravitino mass of O(10^5) GeV in FHI and O(10^9) GeV in DHI. Our results point to interesting scenarios in well-motivated parameter regions that will be tested in future axion and CMB experiments.


title: Multiple impurities and combined local density approximations in Site-Occupation Embedding Theory

abstract: Site-occupation embedding theory (SOET) is an in-principle-exact multi-determinantal extension of density-functional theory for model Hamiltonians. Various extensions of recent developments in SOET [Senjean et al., Phys. Rev. B 97, 235105 (2018)] are explored in this work. An important step forward is the generalization of the theory to multiple impurity sites. We also propose a new single-impurity density-functional approximation (DFA) where the density-functional impurity correlation energy of the two-level (2L) Hubbard system is combined with the Bethe ansatz local density approximation (BALDA) to the full correlation energy of the (infinite) Hubbard model. In order to test the new DFAs, the impurity-interacting wavefunction has been computed self-consistently with the density matrix renormalization group method (DMRG). Double occupation and per-site energy expressions have been derived and implemented in the one-dimensional case. A detailed analysis of the results is presented, with a particular focus on the errors induced either by the energy functionals solely or by the self-consistently converged densities. Among all the DFAs (including those previously proposed), the combined 2L-BALDA is the one that performs the best in all correlation and density regimes. Finally, extensions in new directions, like a partition-DFT-type reformulation of SOET, a projection-based SOET approach, or the combination of SOET with Green functions, are briefly discussed as a perspective.


title: Short-Range Neutrinoless Double Beta Decay Mechanisms

abstract: Neutrinoless double beta decay can significantly help to shed light on the issue of non-zero neutrino mass, as observation of this lepton number violating process would imply neutrinos are Majorana particles. However, the underlying interaction does not have to be as simple as the standard neutrino mass mechanism. The entire variety of neutrinoless double beta decay mechanisms can be approached effectively. In this work we focus on a theoretical description of short-range effective contributions to neutrinoless double beta decay, which are equivalent to 9-dimensional effective operators incorporating the appropriate field content. We give a detailed derivation of the nuclear matrix elements and phase space factors corresponding to individual terms of the effective Lagrangian. Using these, we provide general formulae for the neutrinoless double beta decay half-life and angular correlation of the outgoing electrons.


title: The Long-Term Secular Mass Accretion Rate of the Recurrent Nova T Pyxidis

abstract: We present Hubble Space Telescope ultraviolet spectroscopy of the recurrent nova T Pyxidis obtained more than 5 years after its 2011 outburst indicating that the system might not have yet reached its deep quiescent state. The ultraviolet data exhibit a 20% decline in the continuum flux from the pre-outburst deep quiescence state to the post-outburst near quiescent state. We suggest that a decline across each recurring nova eruption might help explain the proposed 2mag steady decline of the system since 1866. Using an improved version of our accretion disk model as well as International Ultraviolet Explorer ultraviolet and optical data, and the 4.8 kpc distance, we corroborate our previous findings that the quiescent mass accretion rate in T Pyx is of the order of 1e-6 Solar mass per year. Such a large mass accretion rate would imply that the mass of the white dwarf is increasing with time. However, with the just-release Gaia DR 2 distance of 3.3 kpc (after submission of the first version of this manuscript), we find a mass accretion of the order of 1e-7 Solar mass per year. Our results predict powerful soft X-ray or extreme ultraviolet emission from the hot inner region of the high accretion rate disk. Using constraining X-ray observations and assuming the accretion disk doesn't depart too much from the standard model, we are left with two possible scenarios. The disk either emits mainly extreme ultraviolet radiation which, at a distance of a few kpc, is completely absorbed by the interstellar medium, or the hot inner disk, emitting soft X-rays, is masked by the bulging disk seen at a higher inclination.


abstract: We treat the general theory of nonlinear ideals and extend as many notions as possible from the linear theory to the nonlinear theory. We define nonlinear ideals with special properties which associate new non-linear ideals to given ones and establish several properties and characterizations of them. Building upon the results of U. Matter we define a Lipschitz interpolative nonlinear ideal procedure between metric spaces and Banach spaces and establish this class of Lipschitz operators is an injective Banach nonlinear ideal and show several standard basic properties for such class. Extending the work of J. A. L\'{o}pez Molina and E. A. S\'{a}nchez P\'{e}rez we define a Lipschitz $\left(p,\theta, q, \nu\right)$-dominated operators for $1\leq p, q <\infty$; $0\leq \theta, \nu< 1$ and establish several characterizations. Afterwards we generalize a notion of Lipschitz interpolative nonlinear ideal procedure between arbitrary metric spaces and prove its a nonlinear ideal. Finally, we present certain basic counter examples of Lipschitz interpolative nonlinear ideal procedure between arbitrary metric spaces.


title: Computation of option greeks under hybrid stochastic volatility models via Malliavin calculus

abstract: This study introduces computation of option sensitivities (Greeks) using the Malliavin calculus under the assumption that the underlying asset and interest rate both evolve from a stochastic volatility model and a stochastic interest rate model, respectively. Therefore, it integrates the recent developments in the Malliavin calculus for the computation of Greeks: Delta, Vega, and Rho and it extends the method slightly. The main results show that Malliavin calculus allows a running Monte Carlo (MC) algorithm to present numerical implementations and to illustrate its effectiveness. The main advantage of this method is that once the algorithms are constructed, they can be used for numerous types of option, even if their payoff functions are not differentiable.


title: Distributed Optimization Strategy for Multi Area Economic Dispatch Based on Electro Search Optimization Algorithm

abstract: A new adopted evolutionary algorithm is presented in this paper to solve the non-smooth, non-convex and non-linear multi-area economic dispatch (MAED). MAED includes some areas which contains its own power generation and loads. By transmitting the power from the area with lower cost to the area with higher cost, the total cost function can be minimized greatly. The tie line capacity, multi-fuel generator and the prohibited operating zones are satisfied in this study. In addition, a new algorithm based on electro search optimization algorithm (ESOA) is proposed to solve the MAED optimization problem with considering all the constraints. In ESOA algorithm all probable moving states for individuals to get away from or move towards the worst or best solution needs to be considered. To evaluate the performance of the ESOA algorithm, the algorithm is applied to both the original economic dispatch with 40 generator systems and the multi-area economic dispatch with 3 different systems such as: 6 generators in 2 areas; and 40 generators in 4 areas. It can be concluded that, ESOA algorithm is more accurate and robust in comparison with other methods.


title: Probabilistic Trajectory Segmentation by Means of Hierarchical Dirichlet Process Switching Linear Dynamical Systems

abstract: Using movement primitive libraries is an effective means to enable robots to solve more complex tasks. In order to build these movement libraries, current algorithms require a prior segmentation of the demonstration trajectories. A promising approach is to model the trajectory as being generated by a set of Switching Linear Dynamical Systems and inferring a meaningful segmentation by inspecting the transition points characterized by the switching dynamics. With respect to the learning, a nonparametric Bayesian approach is employed utilizing a Gibbs sampler.


title: Protein Design by Algorithm

abstract: We review algorithms for protein design in general. Although these algorithms have a rich combinatorial, geometric, and mathematical structure, they are almost never covered in computer science classes. Furthermore, many of these algorithms admit provable guarantees of accuracy, soundness, complexity, completeness, optimality, and approximation bounds. The algorithms represent a delicate and beautiful balance between discrete and continuous computation and modeling, analogous to that which is seen in robotics, computational geometry, and other fields in computational science. Finally, computer scientists may be unaware of the almost direct impact of these algorithms for predicting and introducing molecular therapies that have gone in a short time from mathematics to algorithms to software to predictions to preclinical testing to clinical trials. Indeed, the overarching goal of these algorithms is to enable the development of new therapeutics that might be impossible or too expensive to discover using experimental methods. Thus the potential impact of these algorithms on individual, community, and global health has the potential to be quite significant.


title: The 21-SPONGE HI Absorption Line Survey II: The temperature of Galactic HI

abstract: We present 21-cm Spectral Line Observations of Neutral Gas with the VLA (21-SPONGE), a Karl G. Jansky Very Large Array (VLA) large project (~600 hours) for measuring the physical properties of Galactic neutral hydrogen (HI). 21-SPONGE is distinguished among previous Galactic HI studies as a result of: (1) exceptional optical depth sensitivity ($\sigma_{\tau} < 10^{-3}$ per $0.42\rm\,km\,s^{-1}$ channels over 57 lines of sight); (2) matching 21 cm emission spectra with highest-possible angular resolution (~4') from the Arecibo Observatory; (3) detailed comparisons with numerical simulations for assessing observational biases. We autonomously decompose 21 cm spectra and derive the physical properties (i.e., spin temperature, $T_s$, column density) of the cold neutral medium (CNM; $T_s<250\rm\,K$), thermally unstable medium (UNM; $250< T_s < 1000\rm\,K$) and warm neutral medium (WNM; $T_s > 1000\rm\,K$) simultaneously. We detect 50% of the total HI mass in absorption, the majority of which is CNM (56 +/- 10%, corresponding to 28% of the total HI mass). Although CNM is detected ubiquitously, the CNM fraction along most lines of sight is <50%. We find that 20% of the total HI mass is thermally unstable (41 +/- 10% of HI detected in absorption), with no significant variation with Galactic environment. Finally, although the WNM comprises 52% of the total HI mass, we detect little evidence for WNM absorption with $1000<T_s<4000\rm\,K$. Following spectral modeling, we detect a stacked residual absorption feature corresponding to WNM with $T_s\sim10^4\rm\,K$. We conclude that excitation in excess of collisions likely produces significantly higher WNM $T_s$ than predicted by steady-state models.


title: Theories of Class F and Anomalies

abstract: We consider the 6d (2,0) theory on a fibration by genus g curves, and dimensionally reduce along the fiber to 4d theories with duality defects. This generalizes class S theories, for which the fibration is trivial. The non-trivial fibration in the present setup implies that the gauge couplings of the 4d theory, which are encoded in the complex structures of the curve, vary and can undergo S-duality transformations. These monodromies occur around 2d loci in space-time, the duality defects, above which the fiber is singular. The key role that the fibration plays here motivates refering to this setup as theories of class F. In the simplest instance this gives rise to 4d N=4 Super-Yang-Mills with space-time dependent coupling that undergoes SL(2, Z) monodromies. We determine the anomaly polynomial for these theories by pushing forward the anomaly polynomial of the 6d (2,0) theory along the fiber. This gives rise to corrections to the anomaly polynomials of 4d N=4 SYM and theories of class S. For the torus case, this analysis is complemented with a field theoretic derivation of a U(1) anomaly in 4d N=4 SYM. The corresponding anomaly polynomial is tested against known expressions of anomalies for wrapped D3-branes with varying coupling, which are known field theoretically and from holography. Extensions of the construction to 4d N = 0 and 1, and 2d theories with varying coupling, are also discussed.


title: The Fourier transform of a projective group frame

abstract: Many tight frames of interest are constructed via their Gramian matrix (which determines the frame up to unitary equivalence). Given such a Gramian, it can be determined whether or not the tight frame is projective group frame, i.e., is the projective orbit of some group $G$ (which may not be unique). On the other hand, there is complete description of the projective group frames in terms of the irreducible projective representations of $G$. Here we consider the inverse problem of taking the Gramian of a projective group frame for a group $G$, and identifying the cocycle and constructing the frame explicitly as the projective group orbit of a vector $v$ (decomposed in terms of the irreducibles). The key idea is to recognise that the Gramian is a group matrix given by a vector $f\in\mathbb{C}^G$, and to take the Fourier transform of $f$ to obtain the components of $v$ as orthogonal projections. This requires the development of a theory of group matrices and the Fourier transform for projective representations. Of particular interest, we give a block diagonalisation of (projective) group matrices. This leads to a unique Fourier decomposition of the group matrices, and a further fine-scale decomposition into low rank group matrices.


title: Detecting Dead Weights and Units in Neural Networks

abstract: Deep Neural Networks are highly over-parameterized and the size of the neural networks can be reduced significantly after training without any decrease in performance. One can clearly see this phenomenon in a wide range of architectures trained for various problems. Weight/channel pruning, distillation, quantization, matrix factorization are some of the main methods one can use to remove the redundancy to come up with smaller and faster models. This work starts with a short informative chapter, where we motivate the pruning idea and provide the necessary notation. In the second chapter, we compare various saliency scores in the context of parameter pruning. Using the insights obtained from this comparison and stating the problems it brings we motivate why pruning units instead of the individual parameters might be a better idea. We propose some set of definitions to quantify and analyze units that don't learn and create any useful information. We propose an efficient way for detecting dead units and use it to select which units to prune. We get 5x model size reduction through unit-wise pruning on MNIST.


title: A New Semiclassical Picture of Vacuum Decay

abstract: We introduce a new picture of vacuum decay which, in contrast to existing semiclassical techniques, provides a real-time description and does not rely on classically-forbidden tunneling paths. Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations of vacuum fluctuations and evolving with the classical equations of motion. The decay rate obtained from an ensemble of simulations is in excellent agreement with existing techniques. Future applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states.


title: Formation of S0 galaxies through mergers. Morphological properties: tidal relics, lenses, ovals, and other inner components

abstract: Major mergers are popularly considered too destructive to produce the relaxed regular structures and the morphological inner components (ICs) usually observed in lenticular (S0) galaxies. We aim to test if major mergers can produce remnants with realistic S0 morphologies. We have selected a sample of relaxed discy remnants resulting from the dissipative merger simulations of the GalMer database and derived their properties mimicking the typical conditions of current observational data. We compare their global morphologies, visual components, and merger relics in mock photometric images with their real counterparts. Only $\sim$1-2 Gyr after the full merger, we find that: 1) many remnants (67 major and 29 minor events) present relaxed structures and typical S0 or E/S0 morphologies, for a wide variety of orbits and even in gas-poor cases. 2) Contrary to popular expectations, most of them do not exhibit any morphological traces of their past merger origin under typical observing conditions and at distances as nearby as 30 Mpc. 3) The merger relics are more persistent in minor mergers than in major ones for similar relaxing time periods. 4) No major-merger S0-like remnant develops a significant bar. 5) Nearly 58% of the major-merger S0 remnants host visually detectable ICs, such as embedded inner discs, rings, pseudo-rings, inner spirals, nuclear bars, and compact sources, very frequent in real S0s too. 6) All remnants contain a lens or oval, identically ubiquitous in local S0s. 7) These lenses and ovals do not come from bar dilution in major merger cases, but are associated with stellar halos or embedded inner discs instead (thick or thin). We conclude that the relaxed morphologies, lenses, ovals, and other ICs of real S0s do not necessarily come from internal secular evolution, gas infall or environmental mechanisms, as traditionally assumed, but they can result from major mergers as well.


title: Nonrelativistic String Theory and T-Duality

abstract: Nonrelativistic string theory in flat spacetime is described by a two-dimensional quantum field theory with a nonrelativistic global symmetry acting on the worldsheet fields. Nonrelativistic string theory is unitary, ultraviolet complete and has a string spectrum and spacetime S-matrix enjoying nonrelativistic symmetry. The worldsheet theory of nonrelativistic string theory is coupled to a curved spacetime background and to a Kalb-Ramond two-form and dilaton field. The appropriate spacetime geometry for nonrelativistic string theory is dubbed string Newton-Cartan geometry, which is distinct from Riemannian geometry. This defines the sigma model of nonrelativistic string theory describing strings propagating and interacting in curved background fields. We also implement T-duality transformations in the path integral of this sigma model and uncover the spacetime interpretation of T-duality. We show that T-duality along the longitudinal direction of the string Newton-Cartan geometry describes relativistic string theory on a Lorentzian geometry with a compact lightlike isometry, which is otherwise only defined by a subtle infinite boost limit. This relation provides a first principles definition of string theory in the discrete light cone quantization (DLCQ) in an arbitrary background, a quantization that appears in nonperturbative approaches to quantum field theory and string/M-theory, such as in Matrix theory. T-duality along a transverse direction of the string Newton-Cartan geometry equates nonrelativistic string theory in two distinct, T-dual backgrounds.


title: Bootstrapping two-loop Feynman integrals for planar N=4 sYM

abstract: We derive analytic results for the symbol of certain two-loop Feynman integrals relevant for seven- and eight-point two-loop scattering amplitudes in planar $\mathcal{N}=4$ super-Yang--Mills theory. We use a bootstrap inspired strategy, combined with a set of second-order partial differential equations that provide powerful constraints on the symbol ansatz. When the complete symbol alphabet is not available, we adopt a hybrid approach. Instead of the full function, we bootstrap a certain discontinuity for which the alphabet is known. Then we write a one-fold dispersion integral to recover the complete result. At six and seven points, we find that the individual Feynman integrals live in the same space of functions as the amplitude, which is described by the 9- and 42-letter cluster alphabets respectively. Starting at eight points however, the symbol alphabet of the MHV amplitude is insufficient for individual integrals. In particular, some of the integrals require algebraic letters involving four-mass box square-root singularities. We point out that these algebraic letters are relevant at the amplitude level directly starting with N$^2$MHV amplitudes even at one loop.


title: The progeny of a Cosmic Titan: a massive multi-component proto-supercluster in formation at z=2.45 in VUDS

abstract: [Abridged] We unveil the complex shape of a proto-supercluster at z~2.45 in the COSMOS field using the spectroscopic redshifts of the VIMOS Ultra-Deep Survey (VUDS), complemented by the zCOSMOS-Deep sample and high-quality photometric redshifts. The method relies on a 2D Voronoi tessellation in overlapping redshift slices that is converted into a 3D density field. The galaxy distribution is constructed using a statistical treatment of spectroscopic and photometric redshifts. We identify a proto-supercluster, dubbed "Hyperion" for its immense size and mass, which extends over a volume of ~60x60x150 comoving Mpc^3 and has an estimated total mass of ~4.8x 10^15 M_sun. This immensely complex structure contains at least 7 density peaks within 2.4 < z < 2.5, connected by filaments. Based on the peaks average matter density, we estimate their total mass, M_tot, and find a range of masses from ~0.1x10^14 M_sun to ~2.7x10^14 M_sun. By using spectroscopic members of each peak, we obtain the velocity dispersion of the galaxies in the peaks, and then their virial mass M_vir (under the strong assumption that they are virialised). The agreement between M_vir and M_tot is surprisingly good, considering that almost all the peaks are probably not yet virialised. According to the spherical collapse model, the peaks have already started or are about to start collapsing, and they are all predicted to be virialised by redshift z~0.8-1.6. We finally perform a careful comparison with the literature, given that smaller components of this proto-supercluster had previously been identified using heterogeneous galaxy samples. With VUDS, we obtain a panoramic view of this large structure, that encompasses, connects, and considerably expands in a homogeneous way on all previous detections of the various sub-components. This provides us the unique possibility to study a rich supercluster in formation.


title: Dwarf carbon stars are likely metal-poor binaries and unlikely hosts to carbon planets

abstract: Dwarf carbon stars make up the largest fraction of carbon stars in the Galaxy with around 1200 candidates known to date primarily from the Sloan Digital Sky Survey. They either possess primordial carbon-enhancements, or are polluted by mass transfer from an evolved companion such that C/O is enhanced beyond unity. To directly test the binary hypothesis, a radial velocity monitoring survey has been carried out on 28 dwarf carbon stars, resulting in the detection of variations in 21 targets. Using Monte Carlo simulations, this detection fraction is found to be consistent with a 100% binary population and orbital periods on the order of hundreds of days. This result supports the post-mass transfer nature of dwarf carbon stars, and implies they are not likely hosts to carbon planets.


title: Stability of exomoons around the Kepler transiting circumbinary planets

abstract: The Kepler mission has detected a number of transiting circumbinary planets (CBPs). Although currently not detected, exomoons could be orbiting some of these CBPs, and they might be suitable for harboring life. A necessary condition for the existence of such exomoons is their long-term dynamical stability. Here, we investigate the stability of exomoons around the Kepler CBPs using numerical $N$-body integrations. We determine regions of stability and obtain stability maps in the (a_m,i_pm) plane, where a_m is the initial exolunar semimajor axis with respect to the CBP, and i_pm is the initial inclination of the orbit of the exomoon around the planet with respect to the orbit of the planet around the stellar binary. Ignoring any dependence on i_pm, for most Kepler CBPs the stability regions are well described by the location of the 1:1 mean motion commensurability of the binary orbit with the orbit of the moon around the CBP. This is related to a destabilizing effect of the binary compared to the case if the binary were replaced by a single body, and which is borne out by corresponding 3-body integrations. For high inclinations, the evolution is dominated by Lidov-Kozai oscillations, which can bring moons in dynamically stable orbits to close proximity within the CBP, triggering strong interactions such as tidal evolution, tidal disruption, or direct collisions. This suggests that there is a dearth of highly-inclined exomoons around the Kepler CBPs, whereas coplanar exomoons are dynamically allowed.


title: Binary Pulsar Distances and Velocities from Gaia Data Release 2

abstract: The second data release from the Gaia mission (Gaia DR2) includes, among its billion entries, as- trometric parameters for binary companions to a number of known pulsars, including white dwarf companions to millisecond pulsars and the non-degenerate components of so-called "black widow" and "redback" systems. We find 22 such counterparts in DR2, of which 12 have statistically sig- nificant measurements of parallax. These DR2 optical proper motions and parallaxes provide new measurements of the distances and transverse velocities of the associated pulsars. For the most part, the results agree with existing radio interferometric and pulsar timing-based astrometry, as well as other distance estimates based on photometry or associations, and for some pulsars they provide the best known distance and velocity estimates. In particular, two of these pulsars have no previous dis- tance measurement: PSR J1227$-$4853, for which Gaia measures a parallax of $0.62\pm 0.16$ mas, and PSR J1431$-$4715, with a Gaia parallax of $0.64\pm 0.16$ mas. Using the Gaia distance measurements, we find that dispersion measure-based distance estimates calculated using the NE2001 and YMW16 Galactic electron density models are on average slightly underestimated, which may be a selection effect due to the over-representation of pulsars at high Galactic latitudes in the present Gaia sample. While the Gaia DR2 results do not quite match the precision that can be achieved by dedicated pulsar timing or radio interferometry, taken together they constitute a small but important improvement to the pulsar distance scale, and the subset of millisecond pulsars with distances measured by Gaia may help improve the sensitivity of pulsar timing arrays to nanohertz gravitational waves.


title: The holographic entropy zoo

abstract: We study the holographic dual of a two parameter family of quantities known as the $\alpha$-$z$ divergences. Many familiar information theoretic quantities occur within this family, including the relative entropy, fidelity, and collision relative entropy. We find explicit bulk expressions for the boundary divergences to second order in a state perturbation whenever $\alpha$ is an integer and $z\geq0$, as well as when $z\in\{0,\infty\}$ and $\alpha\in \mathbb{R}$. Our results apply for perturbations around an arbitrary background state and in any dimension, under the assumption of the equality of bulk and boundary modular flows.


title: Instantaneous braids and Dehn twists in topologically ordered states

abstract: A defining feature of topologically ordered states of matter is the existence of locally indistinguishable states on spaces with non-trivial topology. These degenerate states form a representation of the mapping class group (MCG) of the space, which is generated by braids of defects or anyons, and by Dehn twists along non-contractible cycles. These operations can be viewed as fault-tolerant logical gates in the context of topological quantum error correcting codes and topological quantum computation. Here we show that braids and Dehn twists can in general be implemented by a constant depth quantum circuit, with a depth that is independent of code distance $d$ and system size. The circuit consists of a constant depth local quantum circuit (LQC) implementing a local geometry deformation of the quantum state, followed by a permutation on (relabelling of) the qubits. The permutation requires permuting qubits that are separated by a distance of order $d$; it can be implemented by collective classical motion of mobile qubits or as a constant depth circuit using long-range SWAP operations (with a range set by $d$) on immobile qubits. Applying these results to certain non-Abelian quantum error correcting codes demonstrates how universal logical gate sets can be implemented on encoded qubits using only constant depth unitary circuits.

categories: cond-mat.str-el hep-th math-ph math.MP math.QA quant-ph


title: Soft Bootstrap and Supersymmetry

abstract: The soft bootstrap is an on-shell method to constrain the landscape of effective field theories (EFTs) of massless particles via the consistency of the low-energy S-matrix. Given assumptions on the on-shell data (particle spectra, linear symmetries, and low-energy theorems), the soft bootstrap is an efficient algorithm for determining the possible consistency of an EFT with those properties. The implementation of the soft bootstrap uses the recently discovered method of soft subtracted recursion. We derive a precise criterion for the validity of these recursion relations and show that they fail exactly when the assumed symmetries can be trivially realized by independent operators in the effective action. We use this to show that the possible pure (real and complex) scalar, fermion, and vector exceptional EFTs are highly constrained. Next, we prove how the soft behavior of states in a supermultiplet must be related and illustrate the results in extended supergravity. We demonstrate the power of the soft bootstrap in two applications. First, for the N= 1 and N=2 CP^1 nonlinear sigma models, we show that on-shell constructibility establishes the emergence of accidental IR symmetries. This includes a new on-shell perspective on the interplay between N=2 supersymmetry, low-energy theorems, and electromagnetic duality. We also show that N=2 supersymmetry requires 3-point interactions with the photon that make the soft behavior of the scalar O(1) instead of vanishing, despite the underlying symmetric coset. Second, we study Galileon theories, including aspects of supersymmetrization, the possibility of a vector-scalar Galileon EFT, and the existence of higher-derivative corrections preserving the enhanced special Galileon symmetry. This is addressed by soft bootstrap and by application of double-copy/KLT relations applied to higher-derivative corrections of chiral perturbation theory.


title: The Global 21 cm Absorption from Cosmic Dawn with Inhomogeneous Gas Distribution

abstract: We make an analytical estimate of the maximum 21 cm absorption signal from the cosmic dawn, taking into account the inhomogeneity of gas distribution in the intergalactic medium (IGM) due to non-linear structure formation. The gas located near halos is over-dense but adiabatically heated, while the gas far from halos is under-dense and hence cooler. The cumulative effect of adiabatic heating and cooling from this gas inhomogeneity results in a reduction in the maximum global 21 cm absorption depth by about 40% as compared with homogeneous IGM model, assuming saturated coupling between the spin temperature of neutral hydrogen (HI) and the adiabatic gas kinetic temperature.


title: Uncertain fate of fair sampling in quantum annealing

abstract: Recently, it was demonstrated both theoretically and experimentally on the D-Wave quantum annealer that transverse-field quantum annealing does not find all ground states with equal probability. In particular, it was proposed that more complex driver Hamiltonians beyond transverse fields might mitigate this shortcoming. Here, we investigate the mechanisms of (un)fair sampling in quantum annealing. While higher-order terms can improve the sampling for selected small problems, we present multiple counterexamples where driver Hamiltonians that go beyond transverse fields do not remove the sampling bias. Using perturbation theory we explain why this is the case. In addition, we present large-scale quantum Monte Carlo simulations for spin glasses with known degeneracy in two space dimensions and demonstrate that the fair-sampling performance of quadratic driver terms is comparable to standard transverse-field drivers. Our results suggest that quantum annealing machines are not well suited for sampling applications, unless post-processing techniques to improve the sampling are applied.


title: The 2-Chu-Dialectica construction and the polycategory of multivariable adjunctions

abstract: Cheng, Gurski, and Riehl constructed a cyclic double multicategory of multivariable adjunctions. We show that the same information is carried by a double polycategory, in which opposite categories are polycategorical duals. Moreover, this double polycategory is a full substructure of a double Chu construction, whose objects are a sort of polarized category, and which is a natural home for 2-categorical dualities. We obtain the double Chu construction using a general "Chu-Dialectica" construction on polycategories, which includes both the Chu construction and the categorical Dialectica construction of de Paiva. The Chu and Dialectica constructions each impose additional hypotheses making the resulting polycategory representable (hence *-autonomous), but for different reasons; this leads to their apparent differences.


title: Stability Conditions for Cluster Synchronization in Networks of Heterogeneous Kuramoto Oscillators

abstract: In this paper we study cluster synchronization in networks of oscillators with heterogenous Kuramoto dynamics, where multiple groups of oscillators with identical phases coexist in a connected network. Cluster synchronization is at the basis of several biological and technological processes; yet the underlying mechanisms to enable cluster synchronization of Kuramoto oscillators have remained elusive. In this paper we derive quantitative conditions on the network weights, cluster configuration, and oscillators' natural frequency that ensure asymptotic stability of the cluster synchronization manifold; that is, the ability to recover the desired cluster synchronization configuration following a perturbation of the oscillators' states. Qualitatively, our results show that cluster synchronization is stable when the intra-cluster coupling is sufficiently stronger than the inter-cluster coupling, the natural frequencies of the oscillators in distinct clusters are sufficiently different, or, in the case of two clusters, when the intra-cluster dynamics is homogeneous. We illustrate and validate the effectiveness of our theoretical results via numerical studies.


title: Touching Data: A Discoverability-based Evaluation of a Visualization Interface for Tablet Computers

abstract: While a number of touch-based visualization systems have appeared in recent years, relatively little work has been done to evaluate these systems. The prevailing methods compare these systems to desktop-class applications or utilize traditional training-based usability studies. We argue that existing studies, while useful, fail to address a key aspect of mobile application usage - initial impression and discoverability-driven usability. Over the past few years, we have developed a tablet-based visualization system, Tangere, for analyzing tabular data in a multiple coordinated view configuration. This article describes a discoverability-based user study of Tangere in which the system is compared to a commercially available visualization system for tablets - Tableau's Vizable. The study highlights aspects of each system's design that resonate with study participants, and we reflect upon those findings to identify design principles for future tablet-based data visualization systems.


title: Hyperbolic model for Helmholtz equation with impedance boundary conditions

abstract: Solution of Helmholtz equation with impedance boundary condition on finite interval is equivalently reformulated as steady state of initial boundary value problem for first order hyperbolic system of partial differential equations. Particularly interesting property of the proposed hyperbolic model is that steady state is achieved in finite time. For large wavenumber the numerically challenging task for Helmholtz equation is achieving high accuracy with small number of nodal points. We successfully solved this problem by means of using well balanced scheme approach. Numerical tests demonstrate excellent computational potential of the proposed method: high accuracy is achieved for large wavenumber with small number of nodal points in space and time.


title: Minibatch Gibbs Sampling on Large Graphical Models

abstract: Gibbs sampling is the de facto Markov chain Monte Carlo method used for inference and learning on large scale graphical models. For complicated factor graphs with lots of factors, the performance of Gibbs sampling can be limited by the computational cost of executing a single update step of the Markov chain. This cost is proportional to the degree of the graph, the number of factors adjacent to each variable. In this paper, we show how this cost can be reduced by using minibatching: subsampling the factors to form an estimate of their sum. We introduce several minibatched variants of Gibbs, show that they can be made unbiased, prove bounds on their convergence rates, and show that under some conditions they can result in asymptotic single-update-run-time speedups over plain Gibbs sampling.


title: Coherent quantum dynamics launched by incoherent relaxation in a quantum circuit simulator of a light-harvesting complex

abstract: Engineering and harnessing coherent excitonic transport in organic nanostructures has recently been suggested as a promising way towards improving man-made light harvesting materials. However, realising and testing the dissipative system-environment models underlying these proposals is presently very challenging in supramolecular materials. A promising alternative is to use simpler and highly tunable quantum simulators built from programmable qubits, as recently achieved in a superconducting circuit by Potocnik et al. In this article, we simulate the real-time dynamics of an exciton coupled to a quantum bath as it moves through a network based on the quantum circuit of Potocnik et al. Using the numerically exact hierarchical equations of motion to capture the open quantum system dynamics, we find that an ultrafast but completely incoherent relaxation from a high-lying bright exciton into a doublet of closely spaced dark excitons can spontaneously generate electronic coherences and oscillatory real-space motion across the network (quantum beats). Importantly, we show that this behaviour also survives when the environmental noise is classically stochastic (effectively high temperature), as in present experiments, and also leads to a novel, transient violation of detailed balance in the population relaxation. These predictions highlight the possibilities of designing matched electronic and spectral noise structures for robust coherence generation that does not require coherent excitation or cold environments.


title: Integral geometry on manifolds with boundary and applications

abstract: We survey recent results on inverse problems for geodesic X-ray transforms and other linear and non-linear geometric inverse problems for Riemannian metrics, connections and Higgs fields defined on manifolds with boundary.


title: An absolute sodium abundance for a cloud-free 'hot Saturn' exoplanet

abstract: Broad absorption signatures from alkali metals, such as the sodium (Na I) and potassium (K I) resonance doublets, have long been predicted in the optical atmospheric spectra of cloud-free irradiated gas-giant exoplanets1,2,3. However, observations have only revealed the narrow cores of these features rather than the full pressure-broadened profiles4-6. Cloud and haze opacity at the day-night planetary terminator are considered responsible for obscuring the absorption-line wings, which hinders constraints on absolute atmospheric abundances7-9. Here we present an optical transmission spectrum for the 'hot-Saturn' WASP-96b obtained with the Very Large Telescope, which exhibits the complete pressure-broadened profile of the sodium absorption feature. The spectrum is in excellent agreement with cloud-free, solar-abundance models assuming chemical equilibrium. We are able to measure a precise, absolute sodium abundance of log\epsilon_Na=6.9+0.6-0.4, and use it as a proxy to the planet's atmospheric metallicity relative to the solar value (Z_p/Z_\star=2.3+8.9/--1.7). This result is consistent with the mass-metallicity trend observed for solar-system planets and exoplanets10-12.


title: Finite Groups with a Prescribed Number of Cyclic Subgroups

abstract: Marius T\u{a}rn\u{a}uceanu described the finite groups $G$ having $|G|-1$ cyclic subgroups. We describe the finite groups $G$ having $|G|-\Delta$ cyclic subgroups for $\Delta=2, 3, 4$ and $5$.


title: Solving $(k-1)$-Stable Instances of k-Terminal Cut with Isolating Cuts

abstract: The k-Terminal Cut problem, also known as the Multiway Cut problem, is defined on an edge-weighted graph with $k$ distinct vertices called "terminals." The goal is to remove a minimum weight collection of edges from the graph such that there is no path between any pair of terminals. The problem is NP-hard. Isolating cuts are minimum cuts that separate one terminal from the rest. The union of all the isolating cuts, except the largest, is a $(2-2/k)$-approximation to the optimal k-Terminal Cut. This is the only currently-known approximation algorithm for k-Terminal Cut which does not require solving a linear program. An instance of k-Terminal Cut is $\gamma$-stable if edges in the cut can be multiplied by up to $\gamma$ without changing the unique optimal solution. In this paper, we show that, in any $(k-1)$-stable instance of k-Terminal Cut, the source sets of the isolating cuts are the source sets of the unique optimal solution of that k-Terminal Cut instance. We conclude that the $(2-2/k)$-approximation algorithm returns the optimal solution on $(k-1)$-stable instances. Ours is the first result showing that this $(2-2/k)$-approximation is an exact optimization algorithm on a special class of graphs. We also show that our $(k-1)$-stability result is tight. We construct $(k-1-\epsilon)$-stable instances of the k-Terminal Cut problem which only have trivial isolating cuts: that is, the source set of the isolating cuts for each terminal is just the terminal itself. Thus, the $(2-2/k)$-approximation does not return an optimal solution.


title: Majorana-Hubbard model on the honeycomb lattice

abstract: Phase diagram of a Hubbard model for Majorana fermions on the honeycomb lattice is explored using a combination of field theory, renormalization group and mean-field arguments, as well as exact numerical diagonalization. Unlike the previously studied versions of the model we find that even weak interactions break symmetries and lead to interesting topological phases. We establish two topologically nontrivial phases at weak coupling, one gapped with chiral edge modes, and the other gapless with anti-chiral edge modes. At strong coupling a mapping onto a novel frustrated spin-${1\over 2}$ model suggests a highly entangled spin liquid ground state.


title: Optimizing the nonlinearity and dissipation of a SNAIL Parametric Amplifier for dynamic range

abstract: We present a new quantum-limited Josephson-junction-based 3-wave-mixing parametric amplifier, the SNAIL Parametric Amplifier (SPA), which uses an array of SNAILs (Superconducting Nonlinear Asymmetric Inductive eLements) as the source of tunable nonlinearity. We show how to engineer the nonlinearity over multiple orders of magnitude by varying the physical design of the device. As a function of design parameters, we systematically explore two important amplifier nonidealities that limit dynamic range: the phenomena of gain compression and intermodulation distortion, whose minimization are crucial for high-fidelity multi-qubit readout. Through a comparison with first-principles theory across multiple devices, we demonstrate how to optimize both the nonlinearity and the input-output port coupling of these SNAIL-based parametric amplifiers to achieve higher saturation power, without sacrificing any other desirable characteristics. The method elaborated in our work can be extended to improve all forms of parametrically induced mixing that can be employed for quantum information applications.


title: SPNets: Differentiable Fluid Dynamics for Deep Neural Networks

abstract: In this paper we introduce Smooth Particle Networks (SPNets), a framework for integrating fluid dynamics with deep networks. SPNets adds two new layers to the neural network toolbox: ConvSP and ConvSDF, which enable computing physical interactions with unordered particle sets. We use these lay- ers in combination with standard neural network layers to directly implement fluid dynamics inside a deep network, where the parameters of the network are the fluid parameters themselves (e.g., viscosity, cohesion, etc.). Because SPNets are imple- mented as a neural network, the resulting fluid dynamics are fully differentiable. We then show how this can be successfully used to learn fluid parameters from data, perform liquid control tasks, and learn policies to manipulate liquids.


title: Crime Event Embedding with Unsupervised Feature Selection

abstract: We present a novel event embedding algorithm for crime data that can jointly capture time, location, and the complex free-text component of each event. The embedding is achieved by regularized Restricted Boltzmann Machines (RBMs), and we introduce a new way to regularize by imposing a $\ell_1$ penalty on the conditional distributions of the observed variables of RBMs. This choice of regularization performs feature selection and it also leads to efficient computation since the gradient can be computed in a closed form. The feature selection forces embedding to be based on the most important keywords, which captures the common modus operandi (M. O.) in crime series. Using numerical experiments on a large-scale crime dataset, we show that our regularized RBMs can achieve better event embedding and the selected features are highly interpretable from human understanding.


title: Hidden structural order controls Li-ion transport in cation-disordered oxides for rechargeable lithium batteries

abstract: Crystal structures play a vital role in determining materials properties. In Li-ion cathodes, the crystal structure defines the dimensionality and connectivity of interstitial sites, thus determining Li-ion diffusion kinetics. While a perfect crystal has infinite structural coherence, a class of recently discovered high-capacity cathodes, Li-excess cation-disordered rocksalts, falls on the other end of the spectrum: Their cation sublattices are assumed to be randomly populated by Li and transition metal ions with zero configurational coherence based on conventional X-ray diffraction, such that the Li transport is purely determined by statistical effects. In contrast to this prevailing view, we reveal that cation short-range order, hidden in diffraction, is ubiquitous in these long-range disordered materials and controls the local and macroscopic environments for Li-ion transport. Our work not only discovers a crucial property that has previously been overlooked, but also provides new guidelines for designing and engineering disordered rocksalts cathode materials.


title: Arithmetic Circuits with Locally Low Algebraic Rank

abstract: In recent years, there has been a flurry of activity towards proving lower bounds for homogeneous depth-4 arithmetic circuits, which has brought us very close to statements that are known to imply $\textsf{VP} \neq \textsf{VNP}$. It is open if these techniques can go beyond homogeneity, and in this paper we make some progress in this direction by considering depth-4 circuits of low algebraic rank, which are a natural extension of homogeneous depth-4 circuits. A depth-4 circuit is a representation of an $N$-variate, degree-$n$ polynomial $P$ as \[ P = \sum_{i = 1}^T Q_{i1}\cdot Q_{i2}\cdot \cdots \cdot Q_{it} \; , \] where the $Q_{ij}$ are given by their monomial expansion. Homogeneity adds the constraint that for every $i \in [T]$, $\sum_{j} \operatorname{deg}(Q_{ij}) = n$. We study an extension, where, for every $i \in [T]$, the algebraic rank of the set $\{Q_{i1}, Q_{i2}, \ldots ,Q_{it}\}$ of polynomials is at most some parameter $k$. Already for $k = n$, these circuits are a generalization of the class of homogeneous depth-4 circuits, where in particular $t \leq n$ (and hence $k \leq n$). We study lower bounds and polynomial identity tests for such circuits and prove the following results. We show an $\exp{(\Omega(\sqrt{n}\log N))}$ lower bound for such circuits for an explicit $N$ variate degree $n$ polynomial family when $k \leq n$. We also show quasipolynomial hitting sets when the degree of each $Q_{ij}$ and the $k$ are at most $\operatorname{poly}(\log n)$. A key technical ingredient of the proofs, which may be of independent interest, is a result which states that over any field of characteristic zero, up to a translation, every polynomial in a set of polynomials can be written as a function of the polynomials in a transcendence basis of the set. We combine this with methods based on shifted partial derivatives to obtain our final results.


title: Unsupervised Training for 3D Morphable Model Regression

abstract: We present a method for training a regression network from image pixels to 3D morphable model coordinates using only unlabeled photographs. The training loss is based on features from a facial recognition network, computed on-the-fly by rendering the predicted faces with a differentiable renderer. To make training from features feasible and avoid network fooling effects, we introduce three objectives: a batch distribution loss that encourages the output distribution to match the distribution of the morphable model, a loopback loss that ensures the network can correctly reinterpret its own output, and a multi-view identity loss that compares the features of the predicted 3D face and the input photograph from multiple viewing angles. We train a regression network using these objectives, a set of unlabeled photographs, and the morphable model itself, and demonstrate state-of-the-art results.


title: K2-260 b: a hot Jupiter transiting an F star, and K2-261 b: a warm Saturn around a bright G star

abstract: We present the discovery and confirmation of two new transiting giant planets from the Kepler extended mission K2. K2-260 b is a hot Jupiter transiting a $V=12.7$ F6V star in K2 Field 13, with a mass and radius of $M_{\star}=1.39_{-0.06}^{+0.05} M_{\odot}$ and $R_{\star}=1.69 \pm 0.03 R_{\odot}$. The planet has an orbital period of $P=2.627$ days, and a mass and radius of $M_P=1.42^{+0.31}_{-0.32} M_J$ and $R_P=1.552^{+0.048}_{-0.057} R_J$. This is the first K2 hot Jupiter with a detected secondary eclipse in the Kepler bandpass, with a depth of $71 \pm 15$ ppm, which we use to estimate a geometric albedo of $A_g\sim0.2$. We also detected a candidate stellar companion at 0.6" from K2-260; we find that it is very likely physically associated with the system, in which case it would be an M5-6V star at a projected separation of $\sim400$ AU. K2-261 b is a warm Saturn transiting a bright ($V=10.5$) G7IV/V star in K2 Field 14. The host star is a metal-rich ([Fe/H]$=0.36 \pm 0.06$), mildly evolved $1.10_{-0.02}^{+0.01} M_{\odot}$ star with $R_{\star}=1.65 \pm 0.04 R_{\odot}$. Thanks to its location near the main sequence turn-off, we can measure a relatively precise age of $8.8_{-0.3}^{+0.4}$ Gyr. The planet has $P=11.633$ days, $M_P=0.223 \pm 0.031 M_J$, and $R_P=0.850^{+0.026}_{-0.022} R_J$, and its orbit is eccentric ($e=0.39 \pm 0.15$). Its brightness and relatively large transit depth make this one of the best known warm Saturns for follow-up observations to further characterize the planetary system.


title: The Limits of Post-Selection Generalization

abstract: While statistics and machine learning offers numerous methods for ensuring generalization, these methods often fail in the presence of adaptivity---the common practice in which the choice of analysis depends on previous interactions with the same dataset. A recent line of work has introduced powerful, general purpose algorithms that ensure post hoc generalization (also called robust or post-selection generalization), which says that, given the output of the algorithm, it is hard to find any statistic for which the data differs significantly from the population it came from. In this work we show several limitations on the power of algorithms satisfying post hoc generalization. First, we show a tight lower bound on the error of any algorithm that satisfies post hoc generalization and answers adaptively chosen statistical queries, showing a strong barrier to progress in post selection data analysis. Second, we show that post hoc generalization is not closed under composition, despite many examples of such algorithms exhibiting strong composition properties.


title: Polarization as a probe of thick dust disk in edge-on galaxies: application to NGC 891

abstract: Radiative transfer models were developed to understand the optical polarizations in edge-on galaxies, which are observed to occur even outside the geometrically thin dust disk, with a scale height of ~ 0.2 kpc. In order to reproduce the vertically extended polarization structure, we find it is essential to include a geometrically thick dust layer in the radiative transfer model, in addition to the commonly-known thin dust layer. The models include polarizations due to both dust scattering and dichroic extinction which is responsible for the observed interstellar polarization in the Milky Way. We also find that the polarization level is enhanced if the clumpiness of the interstellar medium, and the dichroic extinction by vertical magnetic fields in the outer regions of the dust lane are included in the radiative transfer model. The predicted degree of polarization outside the dust lane was found to be consistent with that (ranging from 1% to 4%) observed in NGC 891.


title: Behavior of H-FABP-fatty acid complex in a protein crystal simulation

abstract: Crystallographic data comes from a space-time average over all the unit cells within the crystal, so dynamic phenomena do not contribute significantly to the diffraction data. Many efforts have been made to reconstitute the movement of the macromolecules and explore the microstates that the confined proteins can adopt in the crystalline network. In this paper, we explored different strategies to simulate a heart fatty acid binding proteins (H-FABP) crystal starting from high resolution coordinates obtained at room temperature, describing in detail the procedure to study protein crystals (in particular H-FABP) by means of Molecular Dynamics simulations, and exploring the role of ethanol as a co-solute that can modify the stability of the protein and facilitate the interchange of fatty acids. Also, we introduced crystallographic restraints in our crystal models, according to experimental isotropic B-factors and analyzed the H-FABP crystal motions using Principal Component Analysis, isotropic and anisotropic B-factors. Our results suggest that restrained MD simulations based in experimental B-factors produce lower simulated B-factors than simulations without restraints, leading to more accurate predictions of the temperature factors. However, the systems without positional restraints represent a higher microscopic heterogeneity in the crystal.


title: Robust Approaches to Handling Complex Geometries with Galerkin Difference Methods

abstract: The Galerkin difference (GD) basis is a set of continuous, piecewise polynomials defined using a finite difference like grid of degrees of freedom. The one dimensional GD basis functions are naturally extended to multiple dimensions using the tensor product constructions to quadrilateral elements for discretizing partial differential equations. Here we propose two approaches to handling complex geometries using the GD basis within a discontinuous Galerkin finite element setting: (1) using non-conforming, curvilinear GD elements and (2) coupling affine GD elements with curvilinear simplicial elements. In both cases the (semidiscrete) discontinuous Galerkin method is provably energy stable even when variational crimes are committed and in both cases a weight-adjusted mass matrix is used, which ensures that only the reference mass matrix must be inverted. Additionally, we give sufficient conditions on the treatment of metric terms for the curvilinear, nonconforming GD elements to ensure that the scheme is both constant preserving and conservative. Numerical experiments confirm the stability results and demonstrate the accuracy of the coupled schemes.


title: Entangled end states with fractionalized spin projection in a time-reversal-invariant topological superconducting wire

abstract: We study the ground state and low-energy subgap excitations of a finite wire of a time-reversal-invariant topological superconductor (TRITOPS) with spin-orbit coupling. We solve the problem analytically for a long chain of a specific one-dimensional lattice model in the electron-hole symmetric configuration and numerically for other cases of the same model. We present results for the spin density of excitations in long chains with an odd number of particles. The total spin projection along the axis of the spin-orbit coupling $S_z= \pm 1/2$ is distributed with fractions $\pm 1/4$ localized at both ends, and shows even-odd alternation along the sites of the chain. We calculate the localization length of these excitations and find that it can be well approximated by a simple analytical expression. We show that the energy $E$ of the lowest subgap excitations of the finite chain defines tunneling and entanglement between end states.We discuss the effect of a Zeeman coupling $\Delta_Z$ on one of the ends of the chain only. For $\Delta_Z<E$, the energy difference of excitations with opposite spin orientation is $\Delta_Z/2$, consistent with a spin projection $\pm 1/4$. We argue that these physical features are not model dependent and can be experimentally observed in TRITOPS wires under appropriate conditions.


title: Explicit Solutions for Optimal Resource Extraction Problems under Regime Switching L\'evy Models

abstract: This paper studies the problem of optimally extracting nonrenewable natural resources. Taking into account the fact that the market values of the main natural resources i.e. oil, natural gas, copper,..., etc, fluctuate randomly following global and seasonal macroeconomic parameters, the prices of natural resources are modeled using Markov switching L\'evy processes. We formulate this optimal extraction problem as an infinite-time horizon optimal control problem. We derive closed-form solutions for the value function as well as the optimal extraction policy. Numerical examples are presented to illustrate these results.


title: A Planetary Microlensing Event with an Unusually Red Source Star: MOA-2011-BLG-291

abstract: We present the analysis of planetary microlensing event MOA-2011-BLG-291, which has a mass ratio of $q=(3.8\pm0.7)\times10^{-4}$ and a source star that is redder (or brighter) than the bulge main sequence. This event is located at a low Galactic latitude in the survey area that is currently planned for NASA's WFIRST exoplanet microlensing survey. This unusual color for a microlensed source star implies that we cannot assume that the source star is in the Galactic bulge. The favored interpretation is that the source star is a lower main sequence star at a distance of $D_S=4.9\pm1.3\,$kpc in the Galactic disk. However, the source could also be a turn-off star on the far side of the bulge or a sub-giant in the far side of the Galactic disk if it experiences significantly more reddening than the bulge red clump stars. However, these possibilities have only a small effect on our mass estimates for the host star and planet. We find host star and planet masses of $M_{\rm host} =0.15^{+0.27}_{-0.10}M_\odot$ and $m_p=18^{+34}_{-12}M_\oplus$ from a Bayesian analysis with a standard Galactic model under the assumption that the planet hosting probability does not depend on the host mass or distance. However, if we attempt to measure the host and planet masses with host star brightness measurements from high angular resolution follow-up imaging, the implied masses will be sensitive to the host star distance. The WFIRST exoplanet microlensing survey is expected to use this method to determine the masses for many of the planetary systems that it discovers, so this issue has important design implications for the WFIRST exoplanet microlensing survey.


title: Quantum Resource Theories

abstract: Quantum resource theories (QRTs) offer a highly versatile and powerful framework for studying different phenomena in quantum physics. From quantum entanglement to quantum computation, resource theories can be used to quantify a desirable quantum effect, develop new protocols for its detection, and identify processes that optimize its use for a given application. Particularly, QRTs revolutionize the way we think about familiar properties of physical systems like entanglement, elevating them from just being interesting from a fundamental point of view to being useful in performing practical tasks. The basic methodology of a general QRT involves partitioning all quantum states into two groups, one consisting of free states and the other consisting of resource states. Accompanying the set of free states is a collection of free quantum operations arising from natural restrictions on physical systems, and that consists of all the physical processes allowed by the resource theory and which acts invariantly on the set of free states. The QRT then studies what information processing tasks become possible using the restricted operations. Despite the large degree of freedom in how one defines the free states and free operations, unexpected similarities emerge among different QRTs in terms of resource measures and resource convertibility. As a result, objects that appear quite distinct on the surface, such as entanglement and quantum reference frames, appear to have great similarity on a deeper structural level. In this article we review the general framework of a quantum resource theory, focusing on common structural features, operational tasks, and resource measures. To illustrate these concepts, an overview is provided on some of the more commonly studied QRTs in the literature.


title: Non-Negative Networks Against Adversarial Attacks

abstract: Adversarial attacks against neural networks are a problem of considerable importance, for which effective defenses are not yet readily available. We make progress toward this problem by showing that non-negative weight constraints can be used to improve resistance in specific scenarios. In particular, we show that they can provide an effective defense for binary classification problems with asymmetric cost, such as malware or spam detection. We also show the potential for non-negativity to be helpful to non-binary problems by applying it to image classification.


title: Nanoscale Plasmonic Slot Waveguides for Enhanced Raman Spectroscopy

abstract: We theoretically investigate several types of plasmonic slot waveguides for enhancing the measured signal in Raman spectroscopy, which is a consequence of electric field and Purcell factor enhancements, as well as an increase in light-matter interaction volume and the Raman signal collection efficiency. An intuitive methodology is presented for calculating the accumulated Raman enhancement factor of an ensemble of molecules in waveguide sensing, which exploits an analytical photon Green function expansion in terms of the waveguide normal modes, and we combine this with a quantum optics formalism of the molecule-waveguide interaction to model Raman scattering. We subsequently show how integrated plasmonic slot waveguides can attain significantly higher Raman enhancement factors: $\sim$5.3$\times$ compared to optofluidic fibers and $\sim$3.7$\times$ compared to planar integrated dielectric waveguides, with a device size and thus analyte volume of at least three-orders of magnitude less. We also provide a comprehensive comparison between the different types of plasmonic slot waveguides based on the important figures-of-merit, and determine the optimal approaches to maximize Raman enhancement.


title: Remarks on Fixed Point Assertions in Digital Topology

abstract: Several recent papers in digital topology have sought to obtain fixed point results by mimicking the use of tools from classical topology, such as complete metric spaces and homotopy invariant fixed point theory. We show that in many cases, researchers using these tools have derived conclusions that are incorrect or trivial.


title: Instrumental variables regression

abstract: IV regression in the context of a re-sampling is considered in the work. Comparatively, the contribution in the development is a structural identification in the IV model. The work also contains a multiplier-bootstrap justification.


title: Vibrational Coupling Modifies Spectral Diffusion in Core-Shell Colloidal Quantum Dots

abstract: We perform two-dimensional coherent spectroscopy on CdSe/CdZnS core/shell colloidal quantum dots at cryogenic temperatures. In the two-dimensonal spectra, sidebands due to electronic coupling with CdSe lattice LO-phonon modes are observed that have evolutions deviating from the exponential dephasing expected from Markovian spectral diffusion. Comparison to simulations provides further evidence that the LO-phonon mode coupling itself significantly modifies the exciton lineshapes and results in the underlying spectral diffusion on ultrafast timescales.


title: On the well-posedness of the inviscid multi-layer quasi-geostrophic equations

abstract: The inviscid multi-layer quasi-geostrophic equations are considered over an arbitrary bounded domain. The no-flux but non-homogeneous boundary conditions are imposed to accommodate the free fluctuations of the top and layer interfaces. Using the barotropic and baroclinic modes in the vertical direction, the elliptic system governing the streamfunctions and the potential vorticity is decomposed into a sequence of scalar elliptic boundary value problems, where the regularity theories from the two-dimensional case can be applied. With the initial potential vorticity being essentially bounded, the multi-layer quasi-equations are then shown to be globally well-posed, and the initial and boundary conditions are satisfied in the classical sense.


title: Formalizing Category Theory and Presheaf Models of Type Theory in Nuprl

abstract: This article is the first in a series of articles that explain the formalization of a constructive model of cubical type theory in Nuprl. In this document we discuss only the parts of the formalization that do not depend on the choice of base category. So, it spells out how we make the first steps of our formalization of cubical type theory.


title: Gradings on simple real Lie algebras

abstract: English: This work is a doctoral thesis in mathematics by compendium of four articles. Here we explain, using a language as simple as possible, the results achieved in those articles. The general objective is the classification of gradings on simple real Lie algebras. The text is written in both english and spanish. ----- Espa\~nol: Esta obra es una tesis doctoral en matem\'aticas por compendio de cuatro art\'iculos. Aqu\'i explicamos, utilizando un lenguaje lo m\'as sencillo posible, los resultados alcanzados en esos art\'iculos. El objetivo general es la clasificaci\'on de las graduaciones en las \'algebras de Lie reales simples. El texto est\'a escrito tanto en ingl\'es como en espa\~nol.


title: Stochastic WaveNet: A Generative Latent Variable Model for Sequential Data

abstract: How to model distribution of sequential data, including but not limited to speech and human motions, is an important ongoing research problem. It has been demonstrated that model capacity can be significantly enhanced by introducing stochastic latent variables in the hidden states of recurrent neural networks. Simultaneously, WaveNet, equipped with dilated convolutions, achieves astonishing empirical performance in natural speech generation task. In this paper, we combine the ideas from both stochastic latent variables and dilated convolutions, and propose a new architecture to model sequential data, termed as Stochastic WaveNet, where stochastic latent variables are injected into the WaveNet structure. We argue that Stochastic WaveNet enjoys powerful distribution modeling capacity and the advantage of parallel training from dilated convolutions. In order to efficiently infer the posterior distribution of the latent variables, a novel inference network structure is designed based on the characteristics of WaveNet architecture. State-of-the-art performances on benchmark datasets are obtained by Stochastic WaveNet on natural speech modeling and high quality human handwriting samples can be generated as well.


title: Adjoint for advection schemes on the sphere in ICON model

abstract: Among the most advanced and sophisticated methods for state analysis of an atmospheric system is the four dimensional variational data assimilation. The numerically challenging task of this approach is the development and application of the adjoint model components. For tracer transport in fluid dynamics accuracy of numerical advection schemes is vital. It is even more important for applications in space-time variational data assimilation with adjoint model version. We propose novel straightforward and efficient approach - artificial source term method - for adjoint advection solver development. It has several benefits compared to traditional adjoint model building technique. One of the attractive features of the new approach is that it reuses existing advection solver code, thus resulting into significant reduction of time needed for adjoint solver development. The stability, accuracy and convergence of the adjoint schemes are investigated. The method is implemented and evaluated for the linear advection equation on the sphere in the Icosahedral Nonhydrostatic Model (ICON). Adjoint solvers developed by the conventional and the methods are compared against each other on a collection of standard advection test cases and variational data assimilation test cases developed here. The advantages of the artificial source term method is especially obvious in case of monotonic advection equation solvers, granting, e.g. absence of oscillations and nonphysical negative concentrations.


title: Sparse Principal Component based High-Dimensional Mediation Analysis

abstract: Causal mediation analysis aims to quantify the intermediate effect of a mediator on the causal pathway from treatment to outcome. With multiple mediators, which are potentially causally dependent, the possible decomposition of pathway effects grows exponentially with the number of mediators. Huang and Pan (2016) introduced a principal component analysis (PCA) based approach to address this challenge, in which the transformed mediators are conditionally independent given the orthogonality of the PCs. However, the transformed mediator PCs, which are linear combinations of original mediators, are difficult to interpret. In this study, we propose a sparse high-dimensional mediation analysis approach by adopting the sparse PCA method introduced by Zou and others (2006) to the mediation setting. We apply the approach to a task-based functional magnetic resonance imaging study, and show that our proposed method is able to detect biologically meaningful results related to the identified mediator.


title: Generalized Dynamic Programming Principle and Sparse Mean-Field Control Problems

abstract: In this paper we study optimal control problems in Wasserstein spaces, which are suitable to describe macroscopic dynamics of multi-particle systems. The dynamics is described by a parametrized continuity equation, in which the Eulerian velocity field is affine w.r.t. some variables. Our aim is to minimize a cost functional which includes a control norm, thus enforcing a \emph{control sparsity} constraint. More precisely, we consider a nonlocal restriction on the total amount of control that can be used depending on the overall state of the evolving mass. We treat in details two main cases: an instantaneous constraint on the control applied to the evolving mass and a cumulative constraint, which depends also on the amount of control used in previous times. For both constraints, we prove the existence of optimal trajectories for general cost functions and that the value function is viscosity solution of a suitable Hamilton-Jacobi-Bellmann equation. Finally, we discuss an abstract Dynamic Programming Principle, providing further applications in the Appendix.


title: Rheoscopic Fluids in a Post-Kalliroscope World

abstract: In rheoscopic flow visualization, the working fluid is seeded with small reflective flakes that align preferentially in the flow due to their anisotropy. This leads to directed light scattering, which can be exploited to distinguish qualitatively different regions of the flow. For the past four decades, the gold standard in rheoscopic flow visualization has been Kalliroscope, a commercial product consisting of crystalline guanine particles. Recently, however, worldwide production of crystalline guanine has dropped precipitously, leading the Kalliroscope Corporation to halt production in 2014. Here, we present a short survey of alternative rheoscopic flow visualization techniques and introduce an inexpensive rheoscopic fluid based on stearic acid crystals extracted from shaving cream, which has a performance similar to, and in certain respects superior to, Kalliroscope.


title: Machine learning for prediction of extreme statistics in modulation instability

abstract: A central area of research in nonlinear science is the study of instabilities that drive the emergence of extreme events. Unfortunately, experimental techniques for measuring such phenomena often provide only partial characterization. For example, real-time studies of instabilities in nonlinear fibre optics frequently use only spectral data, precluding detailed predictions about the associated temporal properties. Here, we show how Machine Learning can overcome this limitation by predicting statistics for the maximum intensity of temporal peaks in modulation instability based only on spectral measurements. Specifically, we train a neural network based Machine Learning model to correlate spectral and temporal properties of optical fibre modulation instability using data from numerical simulations, and we then use this model to predict the temporal probability distribution based on high-dynamic range spectral data from experiments. These results open novel perspectives in all systems exhibiting chaos and instability where direct time-domain observations are difficult.

categories: physics.comp-ph cs.LG physics.optics stat.ML


title: Fairness Under Composition

abstract: Algorithmic fairness, and in particular the fairness of scoring and classification algorithms, has become a topic of increasing social concern and has recently witnessed an explosion of research in theoretical computer science, machine learning, statistics, the social sciences, and law. Much of the literature considers the case of a single classifier (or scoring function) used once, in isolation. In this work, we initiate the study of the fairness properties of systems composed of algorithms that are fair in isolation; that is, we study fairness under composition. We identify pitfalls of naive composition and give general constructions for fair composition, demonstrating both that classifiers that are fair in isolation do not necessarily compose into fair systems and also that seemingly unfair components may be carefully combined to construct fair systems. We focus primarily on the individual fairness setting proposed in [Dwork, Hardt, Pitassi, Reingold, Zemel, 2011], but also extend our results to a large class of group fairness definitions popular in the recent literature, exhibiting several cases in which group fairness definitions give misleading signals under composition.


title: On the Relationship between Data Efficiency and Error for Uncertainty Sampling

abstract: While active learning offers potential cost savings, the actual data efficiency---the reduction in amount of labeled data needed to obtain the same error rate---observed in practice is mixed. This paper poses a basic question: when is active learning actually helpful? We provide an answer for logistic regression with the popular active learning algorithm, uncertainty sampling. Empirically, on 21 datasets from OpenML, we find a strong inverse correlation between data efficiency and the error rate of the final classifier. Theoretically, we show that for a variant of uncertainty sampling, the asymptotic data efficiency is within a constant factor of the inverse error rate of the limiting classifier.


title: Supervised Fuzzy Partitioning

abstract: Centroid-based methods including k-means and fuzzy c-means are known as effective and easy-to-implement approaches to clustering purposes in many applications. However, these algorithms cannot be directly applied to supervised tasks. This paper thus presents a generative model extending the centroid-based clustering approach to be applicable to classification and regression tasks. Given an arbitrary loss function, the proposed approach, termed Supervised Fuzzy Partitioning (SFP), incorporates labels information into its objective function through a surrogate term penalizing the empirical risk. Entropy-based regularization is also employed to fuzzify the partition and to weight features, enabling the method to capture more complex patterns, identify significant features, and yield better performance facing high-dimensional data. An iterative algorithm based on block coordinate descent scheme is formulated to efficiently find a local optimum. Extensive classification experiments on synthetic, real-world, and high-dimensional datasets demonstrate that the predictive performance of SFP is competitive with state-of-the-art algorithms such as SVM and random forest. SFP has a major advantage over such methods, in that it not only leads to a flexible, nonlinear model but also can exploit any convex loss function in the training phase without compromising computational efficiency.


title: Impact of Channel Models on the End-to-End Performance of mmWave Cellular Networks

abstract: Communication at mmWave frequencies is one of the major innovations of the fifth generation of cellular networks, because of the potential multi-gigabit data rate given by the large amounts of available bandwidth. The mmWave channel, however, makes reliable communications particularly challenging, given the harsh propagation environment and the sensitivity to blockage. Therefore, proper modeling of the mmWave channel is fundamental for accurate results in system simulations of mmWave cellular networks. Nonetheless, complex models, such as the 3GPP channel model for frequencies above 6 GHz, may introduce a significant overhead in terms of computational complexity. In this paper we investigate the trade offs related to the accuracy and the simplicity of the channel model in end-to-end network simulations, and the impact on the performance evaluation of transport protocols.


title: Tight Bound of Incremental Cover Trees for Dynamic Diversification

abstract: Dynamic diversification---finding a set of data points with maximum diversity from a time-dependent sample pool---is an important task in recommender systems, web search, database search, and notification services, to avoid showing users duplicate or very similar items. The incremental cover tree (ICT) with high computational efficiency and flexibility has been applied to this task, and shown good performance. Specifically, it was empirically observed that ICT typically provides a set with its diversity only marginally ($\sim 1/ 1.2$ times) worse than the greedy max-min (GMM) algorithm, the state-of-the-art method for static diversification with its performance bound optimal for any polynomial time algorithm. Nevertheless, the known performance bound for ICT is 4 times worse than this optimal bound. With this paper, we aim to fill this very gap between theory and empirical observations. For achieving this, we first analyze variants of ICT methods, and derive tighter performance bounds. We then investigate the gap between the obtained bound and empirical observations by using specially designed artificial data for which the optimal diversity is known. Finally, we analyze the tightness of the bound, and show that the bound cannot be further improved, i.e., this paper provides the tightest possible bound for ICT methods. In addition, we demonstrate a new use of dynamic diversification for generative image samplers, where prototypes are incrementally collected from a stream of artificial images generated by an image sampler.


title: An operator splitting scheme for the fractional kinetic Fokker-Planck equation

abstract: In this paper, we develop an operator splitting scheme for the fractional kinetic Fokker-Planck equation (FKFPE). The scheme consists of two phases: a fractional diffusion phase and a kinetic transport phase. The first phase is solved exactly using the convolution operator while the second one is solved approximately using a variational scheme that minimizes an energy functional with respect to a certain Kantorovich optimal transport cost functional. We prove the convergence of the scheme to a weak solution to FKFPE. As a by-product of our analysis, we also establish a variational formulation for a kinetic transport equation that is relevant in the second phase. Finally, we discuss some extensions of our analysis to more complex systems.


title: Characterizing $d-$dimensional quantum channels by means of quantum process tomography

abstract: In this work we propose a simple optical architecture, based on phase-only programmable spatial light modulators, in order to characterize general processes on photonic spatial quantum systems in a $d>2$ Hilbert space. We demonstrate the full reconstruction of typical noises affecting quantum computing, as amplitude shifts, phase shifts, and depolarizing channel in dimension $d=5$. We have also reconstructed simulated atmospheric turbulences affecting a free-space transmission of qudits in dimension $d=4$. In each case, quantum process tomography (QPT) was performed in order to obtain the matrix $\chi$ that fully describe the corresponding quantum channel, $\mathcal{E}$. Fidelities between the states experimentally obtained after go through the channel and the expected ones are above $97\%$.


title: Categorical notions of fibration

abstract: Fibrations over a category $B$, introduced to category theory by Grothendieck, encode pseudo-functors $B^{op} \rightsquigarrow {\bf Cat}$, while the special case of discrete fibrations encode presheaves $B^{op} \to {\bf Set}$. A two-sided discrete variation encodes functors $B^{op} \times A \to {\bf Set}$, which are also known as profunctors from $A$ to $B$. By work of Street, all of these fibration notions can be defined internally to an arbitrary 2-category or bicategory. While the two-sided discrete fibrations model profunctors internally to ${\bf Cat}$, unexpectedly, the dual two-sided codiscrete cofibrations are necessary to model $\cal V$-profunctors internally to $\cal V$-$\bf Cat$.