title: Study of a degenerate reaction-diffusion system arising in particle dynamics with aggregation effects

abstract: In this paper we are interested in a degenerate parabolic system of reaction-diffusion equations arising in biology when studying cell adhesion at the protein level. In this modeling the unknown is the couple of the distribution laws of the freely diffusing proteins and of the fixed ones. Under sufficient conditions on the aggregation and unbinding probabilities, we prove the existence of solution of the considered system, as well as their positivity, boundedness and uniqueness. Moreover, we discuss the stability of the equilibrium solution. Finally, we show that the simpler and particular choice of an affine aggregation function and of a constant unbinding probability do not lead to pattern formation as expected in the application. These analytical results are also supported by some numerical simulation.


title: Statistical model of evolution of brain parcellation

abstract: We study the distribution of brain and cortical area sizes [parcellation units (PUs)] obtained for three species: mouse, macaque, and human. We find that the distribution of PU sizes is close to lognormal. We analyze the mathematical model of evolution of brain parcellation based on iterative fragmentation and specialization. In this model, each existing PU has a probability to be split that depends on PU size only. This model shows that the same evolutionary process may have led to brain parcellation in these three species. Our model suggests that region-to-region (macro) connectivity is given by the outer product form. We show that most experimental data on non-vanishing macaque cortex macroconnectivity (62% for area V1) can be explained by the outer product power-law form suggested by our model. We propose a multiplicative Hebbian learning rule for the macroconnectome that could yield the correct scaling of connection strengths between areas. We thus propose a universal evolutionary model that may have contributed to both brain parcellation and mesoscopic level connectivity in mammals.


title: Analysis of Photonic Quantum Nodes Based on Single Photon Raman Interaction

abstract: The long-standing goal of deterministically controlling a single photon using another was recently realized in various experimental settings. Among these, a particularly attractive demonstration relied on single-photon Raman interaction (SPRINT) in a three-level Lambda-system coupled to a single-mode waveguide. Beyond the ability to control the direction of propagation of one photon by the direction of another photon, this scheme has the potential to perform as a passive quantum memory and a universal quantum gate. Relying on interference, this all-optical, coherent scheme requires no additional control fields, and can therefore form the basis for scalable quantum networks composed of passive quantum nodes that interact with each other only with single photon pulses. Here we present an analytical and numerical study of SPRINT, and characterise its limitations and the parameters for optimal operation. Specifically, we study the effect of losses and the presence of multiple excited states. In both cases we discuss strategies for restoring the high fidelity of the device's operation.


title: A directed search for gravitational waves from Scorpius X-1 with initial LIGO

abstract: We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of 1.3e-24 and 8e-25 are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof of principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of ~1.4 for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.


title: Mean field solution of the Blume-Capel model under a random crystal field

abstract: In this work we investigate the Blume-Capel model with infinite-range ferromagnetic interactions and under the influence of a quenched disorder - a random crystal field. For a suitable choice of the random crystal field the model displays a wealth of multicritical behavior, continuous and first-order transition lines, as well as re-entrant behavior. The resulting phase diagrams show a variety of topologies as a function of the disorder parameter $\textit{p}$. A comparison with recent results on the Blume-Capel model in random crystal field is discussed.


title: How to monitor and mitigate stair-casing in l1 trend filtering

abstract: In this paper we study the estimation of changing trends in time-series using $\ell_1$ trend filtering. This method generalizes 1D Total Variation (TV) denoising for detection of step changes in means to detecting changes in trends, and it relies on a convex optimization problem for which there are very efficient numerical algorithms. It is known that TV denoising suffers from the so-called stair-case effect, which leads to detecting false change points. The objective of this paper is to show that $\ell_1$ trend filtering also suffers from a certain stair-case problem. The analysis is based on an interpretation of the dual variables of the optimization problem in the method as integrated random walk. We discuss consistency conditions for $\ell_1$ trend filtering, how to monitor their fulfillment, and how to modify the algorithm to avoid the stair-case false detection problem.


title: Rigidity results with applications to best constants and symmetry of Caffarelli-Kohn-Nirenberg and logarithmic Hardy inequalities

abstract: We take advantage of a rigidity result for the equation satisfied by an extremal function associated with a special case of the Caffarelli-Kohn-Nirenberg inequalities to get a symmetry result for a larger set of inequali-ties. The main ingredient is a reparametrization of the solutions to the Euler-Lagrange equations and estimates based on the rigidity result. The symmetry results cover a range of parameters which go well beyond the one that can be achieved by symmetrization methods or comparison techniques so far.


title: On modulus of continuity of differentiation operator on weighted Sobolev classes

abstract: In this paper we investigate the modulus of continuity of differential operator of order $k$, $k\in\mathbb{N}$, on the classes of functions defined on half-line that have positive non-increasing continuous majorants of functions and their higher derivatives.


title: Green Bank Telescope and Swift X-ray Telescope Observations of the Galactic Center Radio Magnetar SGR J1745-2900

abstract: We present results from eight months of Green Bank Telescope 8.7-GHz observations and nearly 18 months Swift X-ray telescope observations of the radio magnetar SGR J1745-2900, which is located 2.4" from Sgr A*. We tracked the magnetar's radio flux density, polarization properties, pulse profile evolution, rotation, and single-pulse behavior. We identified two main periods of activity in SGR J1745-2900. The first is characterized by approximately 5.5 months of relatively stable evolution in radio flux density, rotation, and profile shape, while in the second these properties varied substantially. Specifically, a third profile component emerged and the radio flux increased on average, but also became more variable. Bright single pulses are visible and are well described by a log-normal energy distribution at low to moderate energies, but with an excess at high energies. The 2-10 keV flux has decayed steadily since the initial X-ray outburst, in contrast with the radio flux. Our timing analysis includes Green Bank Telescope, Swift, and NuSTAR data. When we include the X-ray data in our analyses, we find that SGR J1745-2900 exhibited a level of timing noise unprecedented in a radio magnetar, though an analysis of only the radio data indicates timing noise at a level similar to that observed in other radio magnetars. We conclude that, while SGR J1745-2900 is similar to other radio magnetars in many regards, it differs by having experienced a period of relative stability in the radio band that now appears to have ended, while the X-ray properties have evolved independently.


title: Training and Operation of an Integrated Neuromorphic Network Based on Metal-Oxide Memristors

abstract: Despite all the progress of semiconductor integrated circuit technology, the extreme complexity of the human cerebral cortex makes the hardware implementation of neuromorphic networks with a comparable number of devices exceptionally challenging. One of the most prospective candidates to provide comparable complexity, while operating much faster and with manageable power dissipation, are so-called CrossNets based on hybrid CMOS/memristor circuits. In these circuits, the usual CMOS stack is augmented with one or several crossbar layers, with adjustable two-terminal memristors at each crosspoint. Recently, there was a significant progress in improvement of technology of fabrication of such memristive crossbars and their integration with CMOS circuits, including first demonstrations of their vertical integration. Separately, there have been several demonstrations of discrete memristors as artificial synapses for neuromorphic networks. Very recently such experiments were extended to crossbar arrays of phase-change memristive devices. The adjustment of such devices, however, requires an additional transistor at each crosspoint, and hence the prospects of their scaling are less impressive than those of metal-oxide memristors, whose nonlinear I-V curves enable transistor-free operation. Here we report the first experimental implementation of a transistor-free metal-oxide memristor crossbar with device variability lowered sufficiently to demonstrate a successful operation of a simple integrated neural network, a single layer-perceptron. The network could be taught in situ using a coarse-grain variety of the delta-rule algorithm to perform the perfect classification of 3x3-pixel black/white images into 3 classes. We believe that this demonstration is an important step towards the implementation of much larger and more complex memristive neuromorphic networks.


title: Exoflops in Two Dimensions

abstract: An exoflop occurs in the gauged linear $\sigma$-model by varying the Kahler form so that a subspace appears to shrink to a point and then reemerge "outside" the original manifold. This occurs for K3 surfaces where a rational curve is "flopped" from inside to outside the K3 surface. We see that whether a rational curve contracts to an orbifold phase or an exoflop depends on whether this curve is a line or conic. We study how the D-brane category of the smooth K3 surface is described by the exoflop and, in particular, find the location of a massless D-brane in the exoflop limit. We relate exoflops to noncommutative resolutions.


title: Novel experimental design for high pressure - high temperature electrical resistance measurements in a 'Paris-Edinburgh' large volume press

abstract: We present a novel experimental design for high sensitivity measurements of the electrical resistance of samples at high pressures (0-6GPa) and high temperatures (300-1000K) in a 'Paris-Edinburgh' type large volume press. Uniquely, the electrical measurements are carried out directly on a small sample, thus greatly increasing the sensitivity of the measurement. The sensitivity to even minor changes in electrical resistance can be used to clearly identify phase transitions in material samples. Electrical resistance measurements are relatively simple and rapid to execute and the efficacy of the present experimental design is demonstrated by measuring the electrical resistance of Pb, Sn and Bi across a wide domain of temperature-pressure phase space and employing it to identify the loci of phase transitions. Based on these results, the phase diagrams of these elements are reconstructed to high accuracy and found to be in excellent agreement with previous studies. In particular, by mapping the locations of several well-studied reference points in the phase diagram of Sn and Bi, it is demonstrated that a standard calibration exists for the temperature and pressure, thus eliminating the need for direct or indirect temperature and pressure measurements. The present technique will allow simple and accurate mapping of phase diagrams under extreme conditions and may be of particular importance in advancing studies of liquid state anomalies.


title: Classification and Reconstruction of High-Dimensional Signals from Low-Dimensional Features in the Presence of Side Information

abstract: This paper offers a characterization of fundamental limits on the classification and reconstruction of high-dimensional signals from low-dimensional features, in the presence of side information. We consider a scenario where a decoder has access both to linear features of the signal of interest and to linear features of the side information signal; while the side information may be in a compressed form, the objective is recovery or classification of the primary signal, not the side information. The signal of interest and the side information are each assumed to have (distinct) latent discrete labels; conditioned on these two labels, the signal of interest and side information are drawn from a multivariate Gaussian distribution. With joint probabilities on the latent labels, the overall signal-(side information) representation is defined by a Gaussian mixture model. We then provide sharp sufficient and/or necessary conditions for these quantities to approach zero when the covariance matrices of the Gaussians are nearly low-rank. These conditions, which are reminiscent of the well-known Slepian-Wolf and Wyner-Ziv conditions, are a function of the number of linear features extracted from the signal of interest, the number of linear features extracted from the side information signal, and the geometry of these signals and their interplay. Moreover, on assuming that the signal of interest and the side information obey such an approximately low-rank model, we derive expansions of the reconstruction error as a function of the deviation from an exactly low-rank model; such expansions also allow identification of operational regimes where the impact of side information on signal reconstruction is most relevant. Our framework, which offers a principled mechanism to integrate side information in high-dimensional data problems, is also tested in the context of imaging applications.

categories: cs.IT cs.CV math.IT math.ST stat.ML stat.TH


title: Vortex Nucleation in a Dissipative Variant of the Nonlinear Schr\"odinger Equation under Rotation

abstract: In the present work, we motivate and explore the dynamics of a dissipative variant of the nonlinear Schr{\"o}dinger equation under the impact of external rotation. As in the well established Hamiltonian case, the rotation gives rise to the formation of vortices. We show, however, that the most unstable mode leading to this instability scales with an appropriate power of the chemical potential $\mu$ of the system, increasing proportionally to $\mu^{2/3}$. The precise form of the relevant formula, obtained through our asymptotic analysis, provides the most unstable mode as a function of the atomic density and the trap strength. We show how these unstable modes typically nucleate a large number of vortices in the periphery of the atomic cloud. However, through a pattern selection mechanism, prompted by symmetry-breaking, only few isolated vortices are pulled in sequentially from the periphery towards the bulk of the cloud resulting in highly symmetric stable vortex configurations with far fewer vortices than the original unstable mode. These results may be of relevance to the experimentally tractable realm of finite temperature atomic condensates.


title: Logical Entropy for Quantum States

abstract: The novel concept of quantum logical entropy is presented and analyzed. We prove several basic properties of this entropy with regard to density matrices. We hereby motivate a different approach for the assignment of quantum entropy to density matrices.


title: Radar In-Band and Out-of-Band Interference into LTE Macro and Small Cell Uplinks in the 3.5 GHz Band

abstract: National Telecommunications and Information Administration (NTIA) has proposed vast exclusions zones between radar and Worldwide Interoperability for Microwave Access (WiMAX) (WiMAX) systems which are also being considered as geographic separations between radars and 3.5 GHz Long Term Evolution (LTE) systems without investigating any changes induced by the distinct nature of LTE as opposed to WiMAX. This paper performs a detailed system-level analysis of the interference effects from shipborne radar systems into LTE systems. Even though the results reveal impacts of radar interference on LTE systems performance, they provide clear indications of conspicuously narrower exclusion zones for LTE vis-\`a-vis those for WiMAX and pave the way toward deploying LTE at 3.5 GHz within the coastline populous areas.


title: MT2 to the Rescue -- Searching for Sleptons in Compressed Spectra at the LHC

abstract: We propose a novel method for probing sleptons in compressed spectra at hadron colliders. The process under study is slepton pair production in R-parity conserving supersymmetry, where the slepton decays to a neutralino LSP of mass close to the slepton mass. In order to pass the trigger and obtain large missing energy, an energetic mono-jet is required. Both leptons need to be detected in order to suppress large standard model backgrounds with one charged lepton. We study variables that can be used to distinguish the signal from the remaining major backgrounds, which include tt, WW+jet, Z+jet, and single top production. We find that the dilepton MT2, bound by the mass difference, can be used as an upper bound to efficiently reduce the backgrounds. It is estimated that sleptons with masses up to about 150 GeV can be discovered at the 14 TeV LHC with 100/fb integrated luminosity.


title: Low-Mass Dielectron Production in pp, p-Pb and Pb-Pb Collisions with ALICE

abstract: The ALICE Collaboration measures the production of low-mass dielectrons in pp, p-Pb and Pb-Pb collisions at the LHC. The main detectors used in the analyses are the Inner Tracking System, Time Projection Chamber and Time-Of-Flight detector, all located around mid-rapidity. The production of virtual photons relative to the inclusive yield in pp collisions is determined by analyzing the dielectron excess with respect to the expected hadronic sources. The direct photon cross section is then calculated and found to be in agreement with NLO pQCD calculations. Results from the invariant mass analysis in p-Pb collisions show an overall agreement between data and hadronic cocktail. In Pb-Pb collisions, uncorrected background-subtracted yields have been extracted in two centrality classes. A feasibility study for LHC run 3 after the ALICE upgrade indicates the possibility for a future measurement of the early effective temperature.


title: Low-Rank Approximation and Completion of Positive Tensors

abstract: Unlike the matrix case, computing low-rank approximations of tensors is NP-hard and numerically ill-posed in general. Even the best rank-1 approximation of a tensor is NP-hard. In this paper, we use convex optimization to develop polynomial-time algorithms for low-rank approximation and completion of positive tensors. Our approach is to use algebraic topology to define a new (numerically well-posed) decomposition for positive tensors, which we show is equivalent to the standard tensor decomposition in important cases. Though computing this decomposition is a nonconvex optimization problem, we prove it can be exactly reformulated as a convex optimization problem. This allows us to construct polynomial-time randomized algorithms for computing this decomposition and for solving low-rank tensor approximation problems. Among the consequences is that best rank-1 approximations of positive tensors can be computed in polynomial time. Our framework is next extended to the tensor completion problem, where noisy entries of a tensor are observed and then used to estimate missing entries. We provide a polynomial-time algorithm that for specific cases requires a polynomial (in tensor order) number of measurements, in contrast to existing approaches that require an exponential number of measurements. These algorithms are extended to exploit sparsity in the tensor to reduce the number of measurements needed. We conclude by providing a novel interpretation of statistical regression problems with categorical variables as tensor completion problems, and numerical examples with synthetic data and data from a bioengineered metabolic network show the improved performance of our approach on this problem.


title: Low-frequency variability and heat transport in a low-order nonlinear coupled ocean-atmosphere model

abstract: We formulate and study a low-order nonlinear coupled ocean-atmosphere model with an emphasis on the impact of radiative and heat fluxes and of the frictional coupling between the two components. This model version extends a previous 24-variable version by adding a dynamical equation for the passive advection of temperature in the ocean, together with an energy balance model. The bifurcation analysis and the numerical integration of the model reveal the presence of low-frequency variability (LFV) concentrated on and near a long-periodic, attracting orbit. This orbit combines atmospheric and oceanic modes, and it arises for large values of the meridional gradient of radiative input and of frictional coupling. Chaotic behavior develops around this orbit as it loses its stability; this behavior is still dominated by the LFV on decadal and multi-decadal time scales that is typical of oceanic processes. Atmospheric diagnostics also reveals the presence of predominant low- and high-pressure zones, as well as of a subtropical jet; these features recall realistic climatological properties of the oceanic atmosphere. Finally, a predictability analysis is performed. Once the decadal-scale periodic orbits develop, the coupled system's short-term instabilities --- as measured by its Lyapunov exponents --- are drastically reduced, indicating the ocean's stabilizing role on the atmospheric dynamics. On decadal time scales, the recurrence of the solution in a certain region of the invariant subspace associated with slow modes displays some extended predictability, as reflected by the oscillatory behavior of the error for the atmospheric variables at long lead times.


title: Thermoelectric transport properties of a T-shaped double quantum dot system in the Coulomb blockade regime

abstract: We investigate the thermoelectric properties of a T-shaped double quantum dot system described by a generalized Anderson Hamiltonian. The system's electrical conduction (G) and the fundamental thermoelectric parameters such as the Seebeck coefficient ($S$) and the thermal conductivity ($\kappa$), along with the system's thermoelectric figure of merit (ZT) are numerically estimated based on a Green's function formalism that includes contributions up to the Hartree-Fock level. Our results account for finite onsite Coulomb interaction terms in both component quantum dots and discuss various ways leading to an enhanced thermoelectric figure of merit for the system. We demonstrate that the presence of Fano resonances in the Coulomb blockade regime is responsible for a strong violation of the Wiedemann-Franz law and a considerable enhancement of the system's figure of merit ($ZT$).


title: Material Recognition in the Wild with the Materials in Context Database

abstract: Recognizing materials in real-world images is a challenging task. Real-world materials have rich surface texture, geometry, lighting conditions, and clutter, which combine to make the problem particularly difficult. In this paper, we introduce a new, large-scale, open dataset of materials in the wild, the Materials in Context Database (MINC), and combine this dataset with deep learning to achieve material recognition and segmentation of images in the wild. MINC is an order of magnitude larger than previous material databases, while being more diverse and well-sampled across its 23 categories. Using MINC, we train convolutional neural networks (CNNs) for two tasks: classifying materials from patches, and simultaneous material recognition and segmentation in full images. For patch-based classification on MINC we found that the best performing CNN architectures can achieve 85.2% mean class accuracy. We convert these trained CNN classifiers into an efficient fully convolutional framework combined with a fully connected conditional random field (CRF) to predict the material at every pixel in an image, achieving 73.1% mean class accuracy. Our experiments demonstrate that having a large, well-sampled dataset such as MINC is crucial for real-world material recognition and segmentation.


title: Reconstruction of Randomly Sampled Sparse Signals Using an Adaptive Gradient Algorithm

abstract: Sparse signals can be recovered from a reduced set of samples by using compressive sensing algorithms. In common methods the signal is recovered in the sparse domain. A method for the reconstruction of sparse signal which reconstructs the remaining missing samples/measurements is recently proposed. The available samples are fixed, while the missing samples are considered as minimization variables. Recovery of missing samples/measurements is done using an adaptive gradient-based algorithm in the time domain. A new criterion for the parameter adaptation in this algorithm, based on the gradient direction angles, is proposed. It improves the algorithm computational efficiency. A theorem for the uniqueness of the recovered signal for given set of missing samples (reconstruction variables) is presented. The case when available samples are a random subset of a uniformly or nonuniformly sampled signal is considered in this paper. A recalculation procedure is used to reconstruct the nonuniformly sampled signal. The methods are illustrated on statistical examples.


title: Quipper: Concrete Resource Estimation in Quantum Algorithms

abstract: Despite the rich literature on quantum algorithms, there is a surprisingly small amount of coverage of their concrete logical design and implementation. Most resource estimation is done at the level of complexity analysis, but actual concrete numbers (of quantum gates, qubits, etc.) can differ by orders of magnitude. The line of work we present here is a formal framework to write, and reason about, quantum algorithms. Specifically, we designed a language, Quipper, with scalability in mind, and we are able to report actual resource counts for seven non-trivial algorithms found in the quantum computer science literature.


title: The Atacama Cosmology Telescope: Lensing of CMB Temperature and Polarization Derived from Cosmic Infrared Background Cross-Correlation

abstract: We present a measurement of the gravitational lensing of the Cosmic Microwave Background (CMB) temperature and polarization fields obtained by cross-correlating the reconstructed convergence signal from the first season of ACTPol data at 146 GHz with Cosmic Infrared Background (CIB) fluctuations measured using the Planck satellite. Using an overlap area of 206 square degrees, we detect gravitational lensing of the CMB polarization by large-scale structure at a statistical significance of 4.5 sigma. Combining both CMB temperature and polarization data gives a lensing detection at 9.1 sigma significance. A B-mode polarization lensing signal is present with a significance of 3.2 sigma. We also present the first measurement of CMB lensing--CIB correlation at small scales corresponding to l > 2000. Null tests and systematic checks show that our results are not significantly biased by astrophysical or instrumental systematic effects, including Galactic dust. Fitting our measurements to the best-fit lensing-CIB cross power spectrum measured in Planck data, scaled by an amplitude A, gives A=1.02 +0.12/-0.18 (stat.) +/-0.06(syst.), consistent with the Planck results.


title: The temperature distribution of horizontal branch stars: methods and first results

abstract: As part of a large project aimed at characterizing the ultraviolet (UV) properties of globular clusters, we present here a theoretical and observational analysis aimed at setting the framework for the determination of horizontal branch (HB) temperature distributions. Indeed this is a crucial information to understand the physical parameters shaping the HB morphology in globular clusters and to interpret the UV emission from unresolved stellar systems. We found that the use of zero age HB color-Teff relations is a robust way to derive effective temperatures of individual HB stars. We investigated the most suitable colors for temperature estimates, and the effect on the color-Teff relations of variations of the initial chemical composition, and of the evolution off the zero age horizontal branch. As a test case, we applied our color-Teff calibrations to the Galactic globular cluster M15. The photometry of M15 has been obtained with the Wide Field and Planetary Camera 2 on board the HST. The HB of M15 turned out to have a multimodal distribution, with a main component peaked at Teff~8,000 K and confined below Teff~10,000 K. The second component is peaked at Teff~14,000 K and extends up to Teff~20,000 K. The vast majority (~95%) of the HB stars in M15 is below 20,000 K, in agreement with the lack of a well populated extreme HB observed in other metal-poor globular clusters. We also verified that the temperatures derived with our analysis are consistent with spectroscopic estimates available in the literature.


title: On the Hamiltonicity of the $k$-regular graph game

abstract: We consider a game played on an initially empty graph where two players alternate drawing an edge between vertices subject to the condition that no degree can exceed $k$. We show that for $k=3$, either player can avoid a Hamilton cycle, and for $k\geq4$, either player can force the resulting graph to be Hamiltonian.


title: On the Unstable Directions and Lyapunov Exponents of Anosov Endomorphisms

abstract: Despite the invertible setting, Anosov endomorphisms may have infinitely many unstable directions. Here we prove, under transitivity assumption, that an Anosov endomorphism on a closed manifold $M,$ is either special (that is, every $x \in M$ has only one unstable direction) or for a typical point in $M$ there are infinitely many unstable directions. Other result of this work is the semi rigidity of the unstable Lyapunov exponent of a $C^{1+\alpha}$ codimension one Anosov endomorphism and $C^1$ close to a linear endomorphism of $\mathbb{T}^n$ for $(n \geq 2).$ In the appendix we give a proof for ergodicity of $C^{1+\alpha}, \alpha > 0,$ conservative Anosov endomorphism.


title: Batch Nonlinear Continuous-Time Trajectory Estimation as Exactly Sparse Gaussian Process Regression

abstract: In this paper, we revisit batch state estimation through the lens of Gaussian process (GP) regression. We consider continuous-discrete estimation problems wherein a trajectory is viewed as a one-dimensional GP, with time as the independent variable. Our continuous-time prior can be defined by any nonlinear, time-varying stochastic differential equation driven by white noise; this allows the possibility of smoothing our trajectory estimates using a variety of vehicle dynamics models (e.g., `constant-velocity'). We show that this class of prior results in an inverse kernel matrix (i.e., covariance matrix between all pairs of measurement times) that is exactly sparse (block-tridiagonal) and that this can be exploited to carry out GP regression (and interpolation) very efficiently. When the prior is based on a linear, time-varying stochastic differential equation and the measurement model is also linear, this GP approach is equivalent to classical, discrete-time smoothing (at the measurement times); when a nonlinearity is present, we iterate over the whole trajectory to maximize accuracy. We test the approach experimentally on a simultaneous trajectory estimation and mapping problem using a mobile robot dataset.


title: Paving over arbitrary MASAs in von Neumann algebras

abstract: We consider a paving property for a maximal abelian *-subalgebra (MASA) $A$ in a von Neumann algebra $M$, that we call so-paving, involving approximation in the so-topology, rather than in norm (as in classical Kadison-Singer paving). If $A$ is the range of a normal conditional expectation, then so-paving is equivalent to norm paving in the ultrapower inclusion $A^\omega\subset M^\omega$. We conjecture that any MASA in any von Neumann algebra satisfies so-paving. We use [MSS13] to check this for all MASAs in $\mathcal B(\ell^2 \mathbb N)$, all Cartan subalgebras in amenable von Neumann algebras and in group measure space II$_1$ factors arising from profinite actions. By [P13], the conjecture also holds true for singular MASAs in II$_1$ factors, and we obtain here an improved paving size $C\varepsilon^{-2}$, which we show to be sharp.


title: Hessian ideals of a homogeneous polynomial and generalized Tjurina algebras

abstract: Using the minors in Hessian matrices, we introduce new graded algebras associated to a homogeneous polynomial. When the associated projective hypersurface has isolated singularities, these algebras are related to some new local algebras associated to isolated hypersurface singularities, which generalize their Tjurina algebras. One consequence of our results is a new way to determine the number of weighted homogeneous singularities of such a hypersurface.


title: Rotational component spaces for infinite-type translation surfaces

abstract: Finite translation surfaces can be classified by the order of their singularities. When generalizing to infinite translation surfaces, however, the notion of order of a singularity is no longer well-defined and has to be replaced by new concepts. This article discusses the nature of two such concepts, recently introduced by Bowman and Valdez: linear approaches and rotational components. We show that there is a large flexibility in the spaces of rotational components and even more in the spaces of linear approaches. In particular, we prove that every finite topological space arises as space of rotational components. However, this space will still not contain enough information to describe an infinite translation surface. We showcase this through an uncountable family with the same space of rotational components but different spaces of linear approaches. Additionally, we study several known and new examples to illustrate the concept of linear approaches and rotational components.


title: Reversing Hydride Ion Formation in Quantum Information Experiments with Be$^+$

abstract: We demonstrate photodissociation of BeH$^+$ ions within a Coulomb crystal of thousands of $^9$Be$^+$ ions confined in a Penning trap. Because BeH$^+$ ions are created via exothermic reactions between trapped, laser-cooled Be$^+$($^2\text{P}_{3/2}$) and background H$_2$ within the vacuum chamber, they represent a major contaminant species responsible for infidelities in large-scale trapped-ion quantum information experiments. The rotational-state-insensitive dissociation scheme described here makes use of 157 nm photons to produce Be$^+$ and H as products, thereby restoring Be$^+$ ions without the need for reloading. This technique facilitates longer experiment runtimes at a given background H$_2$ pressure, and may be adapted for removal of MgH$^+$ and AlH$^+$ impurities.


title: Geometry and fluxes of SL(5) exceptional field theory

abstract: We use a geometric approach to construct a flux formulation for the SL(5) U-duality manifest exceptional field theory. The resulting formalism is well-suited for studying gauged supergravities with geometric and non-geometric fluxes. Here we describe all such fluxes for both M-theory and IIB supergravity including the Ramond-Ramond fields for compactifications to seven dimensions. We define the locally non-geometric "R-flux" and globally non-geometric "Q-flux" for M-theory and find a new locally non-geometric R-flux for the IIB theory. We show how these non-geometric fluxes can be understood geometrically and give some examples of how they fluxes can be generated by acting with dualities on solutions with geometric or field-strength flux.


title: Experimental Verification of Van Vleck Nature of Long-Range Ferromagnetic Order in Vanadium-Doped Three-Dimensional Topological Insulator Sb$_{2}$Te$_{3}$

abstract: We demonstrate by high resolution low temperature electron energy loss spectroscopy (EELS) measurements that the long range ferromagnetic (FM) order in vanadium (V)-doped topological insulator Sb$_2$Te$_3$ has the nature of van Vleck-type ferromagnetism. The positions and the relative amplitudes of two core-level peaks (L$_3$ and L$_2$) of the V EELS spectrum show unambiguous change when the sample is cooled from room temperature to T=10K. Magnetotransport and comparison of the measured and simulated EELS spectra confirm that these changes originate from onset of FM order. Crystal field analysis indicates that in V-doped Sb$_2$Te$_3$, partially filled core states contribute to the FM order. Since van Vleck magnetism is a result of summing over all states, this magnetization of core level verifies the van Vleck-type ferromagnetism in a direct manner.


title: ADM Analysis of Gravity Models within the Framework of Bimetric Variational Formalism

abstract: Bimetric variational formalism was recently employed to construct novel bimetric gravity models. In these models an affine connection is generated by an additional tensor field which is independent of the physical metric. In this work we demonstrate how the ADM decomposition can be applied to study such models and provide some technical intermediate details. Using ADM decomposition we are able to prove that a linear model is unstable as has previously been indicated by perturbative analysis. Moreover, we show that it is also very difficult if not impossible to construct a non-linear model which is ghost-free within the framework of bimetric variational formalism. However, we demonstrate that viable models are possible along similar lines of thought. To this end, we consider a set up in which the affine connection is a variation of the Levi-Civita one. As a proof of principle we construct a gravity model with a massless scalar field obtained this way.


title: Symplectic instanton bundles on P3 and 't Hooft instantons

abstract: We study the moduli space $I_{n,r}$ of rank-$2r$ symplectic instanton vector bundles on $\mathbb{P}^3$ with $r\ge2$ and second Chern class $n\ge r+1,\ n-r\equiv 1(\mathrm{mod}2)$. We introduce the notion of tame symplectic instantons by excluding a kind of pathological monads and show that the locus $I^*_{n,r}$ of tame symplectic instantons is irreducible and has the expected dimension, equal to $4n(r+1)-r(2r+1)$. The proof is inherently based on a relation between the spaces $I^*_{n,r}$ and the moduli spaces of 't Hooft instantons


title: Beating the Generator-Enumeration Bound for Solvable-Group Isomorphism

abstract: We consider the isomorphism problem for groups specified by their multiplication tables. Until recently, the best published bound for the worst-case was achieved by the n^(log_p n + O(1)) generator-enumeration algorithm. In previous work with Fabian Wagner, we showed an n^((1 / 2) log_p n + O(log n / log log n)) time algorithm for testing isomorphism of p-groups by building graphs with degree bounded by p + O(1) that represent composition series for the groups and applying Luks' algorithm for testing isomorphism of bounded degree graphs. In this work, we extend this improvement to the more general class of solvable groups to obtain an n^((1 / 2) log_p n + O(log n / log log n)) time algorithm. In the case of solvable groups, the composition factors can be large which prevents previous methods from outperforming the generator-enumeration algorithm. Using Hall's theory of Sylow bases, we define a new object that generalizes the notion of a composition series with small factors but exists even when the composition factors are large. By constructing graphs that represent these objects and running Luks' algorithm, we obtain our algorithm for solvable-group isomorphism. We also extend our algorithm to compute canonical forms of solvable groups while retaining the same complexity.


title: Computationally-efficient stochastic cluster dynamics method for modeling damage accumulation in irradiated materials

abstract: An improved version of a recently developed stochastic cluster dynamics (SCD) method {[}Marian, J. and Bulatov, V. V., {\it J. Nucl. Mater.} \textbf{415} (2014) 84-95{]} is introduced as an alternative to rate theory (RT) methods for solving coupled ordinary differential equation (ODE) systems for irradiation damage simulations. SCD circumvents by design the curse of dimensionality of the variable space that renders traditional ODE-based RT approaches inefficient when handling complex defect population comprised of multiple (more than two) defect species. Several improvements introduced here enable efficient and accurate simulations of irradiated materials up to realistic (high) damage doses characteristic of next-generation nuclear systems. The first improvement is a procedure for efficiently updating the defect reaction-network and event selection in the context of a dynamically expanding reaction-network. Next is a novel implementation of the $\tau$-leaping method that speeds up SCD simulations by advancing the state of the reaction network in large time increments when appropriate. Lastly, a volume rescaling procedure is introduced to control the computational complexity of the expanding reaction-network through occasional reductions of the defect population while maintaining accurate statistics. The enhanced SCD method is then applied to model defect cluster accumulation in iron thin films subjected to triple ion-beam ($\text{Fe}^{3+}$, $\text{He}^{+}$ and $ $$\text{H\ensuremath{{}^{+}}}$$ $) irradiations, for which standard RT or spatially-resolved kinetic Monte Carlo simulations are prohibitively expensive.


title: Kerr geometry in f(T) gravity

abstract: Null tetrads are shown to be a valuable tool in teleparallel theories of modified gravity. We use them to prove that Kerr geometry remains a solution for a wide family of f(T) theories of gravity.


title: Two-generator subgroups of right-angled Artin groups are quasi-isometrically embedded

abstract: We show that non-abelian two-generator subgroups of right-angled Artin groups are quasi-isometrically embedded free groups. This provides an alternate proof of a theorem of A. Baudisch: that all two-generator subgroups are free or free abelian. Additionally, it shows that they are quasi-isometrically embedded. Our theorem also gives a method for detecting groups that are not isomorphic to a subgroup of any RAAG. We present some counterexamples in subgroups with more than two generators.


title: Extremal shift rule for continuous-time zero-sum Markov games

abstract: In the paper we consider the controlled continuous-time Markov chain describing the interacting particles system with the finite number of types. The system is controlled by two players with the opposite purposes. The limiting game as the number of particles tends to infinity is a zero-sum differential game. Krasovskii--Subbotin extremal shift provides the optimal strategy in the limiting game. The main result of the paper is the near optimality of the Krasovskii--Subbotin extremal shift rule for the original Markov game.


title: Mapping of Secondary Virtual Networks onto Wireless Substrate based on Cognitive Radio: multi-objective formulation and analysis

abstract: There is a growing demand for wireless services with different requirements in a dense and heterogeneous wireless environment. Handling this complex ecosystem is becoming a challenging issue, and wireless virtualization emerges as an efficient solution. Although the inclusion of virtualization in wireless networks ensures a better use of resources, current approaches adopted for wireless virtualization can cause an underutilization of resources, since the resource allocated to a virtual wireless network is not shared with other one. This problem can be overcome by combining wireless virtualization with the cognitive radio (CR) technology and dynamic access spectrum (DSA) techniques. Thus, virtual wireless networks with different access priorities to resources (e.g., primary and secondary) can be deployed in an overlay form and share the same substrate wireless network, where the secondary networks use the resources opportunistically. However, challenges emerge in this new scenario, ranging from the mapping to the operation of these networks. This paper is the first to propose the cognitive radio virtual network environment and to formulate the problem of mapping the secondary virtual networks (SVNs) onto the wireless substrate network based on CR. To this end, a multi-objective formulation is designed. Moreover, an analysis of the metrics defined in the problem is performed in order to give the reader useful assistance in designing schemes to solve the problem of mapping SVNs onto substrate networks.


title: Adaptive False Discovery Rate Control for Heterogeneous Data

abstract: Efforts to develop more efficient multiple hypothesis testing procedures for false discovery rate (FDR) control have focused on incorporating an estimate of the proportion of true null hypotheses (such procedures are called adaptive) or exploiting heterogeneity across tests via some optimal weighting scheme. This paper combines these approaches using a weighted adaptive multiple decision function (WAMDF) framework. Optimal weights for a flexible random effects model are derived and a WAMDF that controls the FDR for arbitrary weighting schemes when test statistics are independent under the null hypotheses is given. Asymptotic and numerical assessment reveals that, under weak dependence, the proposed WAMDFs provide more efficient FDR control even if optimal weights are misspecified. The robustness and flexibility of the proposed methodology facilitates the development of more efficient, yet practical, FDR procedures for heterogeneous data. To illustrate, two different weighted adaptive FDR methods for heterogeneous sample sizes are developed and applied to data.


title: Explicit Multi-Matrix Topological Expansion for Quaternionic Random Matrices

abstract: We present an explicit formula for the expected value of a product of several independent symplectically invariant matrices in which the trace and real part function may be applied, possibly to different subexpressions. This takes the form of a topological expansion; however, each term has two topologies: one for the trace, and another for the real part. The traces and real parts can always be written in terms of index contraction, but in some cases, it is possible to write the expression as a product in which the two functions are applied to bracketed intervals in a legal bracket diagram. We present the conditions under which this may be done, and an algorithm to construct such an expression given the contracted indices when possible. The summands in the topological expansion are written in terms of matrix cumulants. We compute the matrix cumulants of quaternionic Ginibre, Gaussian symplectic, quaternionic Wishart, and Haar-distributed symplectic matrices, which allow direct computation of an expression constructed from several independent ensembles of any of these matrices.


title: Identification of Observables for Quark and Gluon Orbital Angular Momentum

abstract: A new debate has recently arisen on the subject of orbital angular momentum in QCD, in particular on its observability and on its partonic interpretation. Orbital momentum can be defined in QCD using two different decomposition schemes that yield a kinetic and a canonical definition, respectively. We argue that kinetic orbital angular momentum is intrinsically associated with twist three generalized parton distributions, and it is therefore more readily observable, while, due to parity constraints, canonical angular momentum, if defined as suggested in the literature in terms of generalized transverse momentum distributions, cannot be observed in scattering processes involving a single hadronic reaction plane.


title: Uniform stability of twisted constant scalar curvature K\"ahler metrics

abstract: We introduce a norm on the space of test configurations, which we call the minimum norm. We conjecture that uniform K-stability with respect to this norm is equivalent to the existence of a constant scalar curvature K\"ahler metric. This notion of uniform K-stability is analogous to coercivity of the Mabuchi functional. We characterise the triviality of test configurations, by showing that a test configuration has zero minimum norm if and only if it has zero $L^2$-norm, if and only if it is almost trivial. We prove that the existence of a twisted constant scalar curvature K\"ahler metric implies uniform twisted K-stability with respect to the minimum norm, when the twisting is ample. We give algebro-geometric proofs of uniform K-stability in the general type and Calabi-Yau cases, as well as in the Fano case under an alpha invariant condition. Our results hold for line bundles sufficiently close to the (anti)-canonical line bundle, and also in the twisted setting. We show that log K-stability implies twisted K-stability, and also that twisted K-semistability of a variety implies that the variety has mild singularities.


title: Libcint: An efficient general integral library for Gaussian basis functions

abstract: An efficient integral library Libcint was designed to automatically implement general integrals for Gaussian-type scalar and spinor basis functions. The library can handle arbitrary integral expressions on top of $\mathbf{p}$, $\mathbf{r}$ and $\sigma$ operators with one-electron overlap and nuclear attraction, two-electron Coulomb and Gaunt operators. Using a symbolic algebra tool, new integrals are derived and translated to C code programmatically. The generated integrals can be used in various types of molecular properties. In the present work, we computed the analytical gradients and NMR shielding constants at both non-relativistic and four-component relativistic Hartree-Fock level to demonstrate the capability of the integral library. Due to the use of kinetically balanced basis and gauge including atomic orbitals, the relativistic analytical gradients and shielding constants requires the integral library to handle the fifth-order electron repulsion integral derivatives. The generality of the integral library is achieved without losing efficiency. On the modern multi-CPU platform, Libcint can easily reach the overall throughput being many times of the I/O bandwidth. On a 20-core node, we are able to achieve an average output 7.9 GB/s for C$_{60}$ molecule with cc-pVTZ basis.


title: A review of "Mem-computing NP-complete problems in polynomial time using polynomial resources" (arXiv:1411.4798)

abstract: The reviewed paper describes an analog device that empirically solves small instances of the NP-complete Subset Sum Problem (SSP). The authors claim that this device can solve the SSP in polynomial time using polynomial space, in principle, and observe no exponential scaling in resource requirements. We point out that (a) the properties ascribed by the authors to their device are insufficient to solve NP-complete problems in poly-time, (b) runtime analysis offered does not cover the spectral measurement step, (c) the overall technique requires exponentially increasing resources when scaled up because of the spectral measurement step.


title: Evolution Differential Equations in Fr\'echet Space with Schauder Basis

abstract: We consider evolution differential equations in Fr\'echet spaces that possess unconditional Schauder basis and construct a version of the majorant functions method to obtain existence theorems for Cauchy problems. Applications to PDE and ODE have been considered.


title: Still Simpler Way of Introducing Interior-Point method for Linear Programming

abstract: Linear Programming is now included in Algorithm undergraduate and postgraduate courses for Computer Science majors. It is possible to teach interior-point methods directly with just minimal knowledge of Algebra and Matrices.


title: Time-Varying Potassium in High-Resolution Spectra of the Type Ia Supernova 2014J

abstract: We present a time series of the highest resolution spectra yet published for the nearby Type Ia supernova (SN) 2014J in M82. They were obtained at 11 epochs over 33 days around peak brightness with the Levy Spectrograph (resolution R~110,000) on the 2.4m Automated Planet Finder telescope at Lick Observatory. We identify multiple Na I D and K I absorption features, as well as absorption by Ca I H & K and several of the more common diffuse interstellar bands (DIBs). We see no evolution in any component of Na I D, Ca I, or in the DIBs, but do establish the dissipation/weakening of the two most blueshifted components of K I. We present several potential physical explanations, finding the most plausible to be photoionization of circumstellar material, and discuss the implications of our results with respect to the progenitor scenario of SN 2014J.


title: Series solutions of confluent Heun equations in terms of incomplete Gamma-functions

abstract: We present a simple systematic algorithm for construction of expansions of the solutions of ordinary differential equations with rational coefficients in terms of mathematical functions having indefinite integral representation. The approach employs an auxiliary equation involving only the derivatives of a solution of the equation under consideration. Using power-series expansions of the solutions of this auxiliary equation, we construct several expansions of the four confluent Heun equations' solutions in terms of the incomplete Gamma-functions. In the cases of single- and double-confluent Heun equations the coefficients of the expansions obey four-term recurrence relations, while for the bi- and tri-confluent Heun equations the recurrence relations in general involve five terms. Other expansions for which the expansion coefficients obey recurrence relations involving more terms are also possible. The particular cases when these relations reduce to ones involving less number of terms are identified. The conditions for deriving closed-form finite-sum solutions via right-hand side termination of the constructed series are discussed.


title: Confirmation of a Steep Luminosity Function for Lyman-alpha Emitters at z = 5.7: A Major Component of Reionization

abstract: We report the first direct and robust measurement of the faint-end slope of the Lyman-alpha emitter (LAE) luminosity function at z = 5.7. Candidate LAEs from a low-spectral-resolution blind search with IMACS on Magellan-Baade were targeted at higher resolution to distinguish high redshift LAEs from foreground galaxies. All but 2 of our 42 single-emission-line systems have flux F $< 2.0 \times 10^{-17}$ ergs s$^{-1}$ cm$^{-2}$, making these the faintest emission-lines observed for a z = 5.7 sample with known completeness, an essential property for determining the faint end slope of the LAE luminosity function. We find 13 LAEs as compared to 29 foreground galaxies, in very good agreement with the modeled foreground counts predicted in Dressler et al. (2011a) that had been used to estimate a faint-end slope of $\alpha$ = -2.0 for the LAE luminosity function. A 32% LAE fraction, LAE/(LAE+foreground), within the flux interval F = $2-20 \times 10^{-18}$ ergs s$^{-1}$ cm$^{-2}$, constrains the faint end slope of the luminosity function to -2.35 < $\alpha$ < -1.95 (1-$\sigma$). We show how this steep LF should provide, to the limit of our observations, more than 20% of the flux necessary to maintain ionization at z=5.7, with a factor-of-ten extrapolation in flux reaching more than 55%. This is in addition to a comparable contribution from Lyman Break Galaxies M$_{UV} \le$ -18. We suggest that this bodes well for a sufficient supply of Lyman continuum photons by similar, low-mass star forming galaxies within the reionization epoch at z $\approx$ 7, only 250 Myr earlier.


title: Equation-of-State Parameter for Reheating

abstract: Constraints to the parameters of inflation models are often derived assuming some plausible range for the number--e.g., $N_k=46$ to $N_k=60$--of $e$-folds of inflation that occurred between the time that our current observable Universe exited the horizon and the end of inflation. However, that number is, for any specific inflaton potential, related to an effective equation-of-state parameter $w_{\mathrm{re}}$ and temperature $T_{\mathrm{re}}$, for reheating. Although the physics of reheating is highly uncertain, there is a finite range of reasonable values for $w_{\mathrm{re}}$. Here we show that by restricting $w_{\mathrm{re}}$ to this range, more stringent constraints to inflation-model parameters can be derived than those obtained from the usual procedure. To do so, we focus in this work in particular on natural inflation and inflation with a Higgs-like potential, and on power law models as limiting cases of those. As one example, we show that the lower limit to the tensor-to-scalar ratio $r$, derived from current measurements of the scalar spectral index, is about 20%-25% higher (depending on the model) with this procedure than with the usual approach.


title: Maximum speed of hypervelocity stars ejected from binaries

abstract: The recent detection of hypervelocity stars (HVSs) as late-type B-stars and HVS candidate G/K-dwarfs raises the important question of their origin. In this Letter, we investigate the maximum possible velocities of such HVSs if they are produced from binaries which are disrupted via an asymmetric supernova explosion. We find that HVSs up to ~770 and ~1280 km/s are possible in the Galactic rest frame from this scenario for these two subclasses of HVSs, respectively. We conclude that whereas a binary origin cannot easily explain all of the observed velocities of B-type HVSs (in agreement with their proposed central massive black hole origin) it can indeed account for the far majority (if not all) of the recently detected G/K-dwarf HVS candidates.


title: Mock Quasar-Lyman-{\alpha} Forest Data-sets for the SDSS-III Baryon Oscillation Spectroscopic Survey

abstract: We describe mock data-sets generated to simulate the high-redshift quasar sample in Data Release 11 (DR11) of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The mock spectra contain Ly{\alpha} forest correlations useful for studying the 3D correlation function including Baryon Acoustic Oscillations (BAO). They also include astrophysical effects such as quasar continuum diversity and high-density absorbers, instrumental effects such as noise and spectral resolution, as well as imperfections introduced by the SDSS pipeline treatment of the raw data. The Ly{\alpha} forest BAO analysis of the BOSS collaboration, described in Delubac et al. 2014, has used these mock data-sets to develop and cross-check analysis procedures prior to performing the BAO analysis on real data, and for continued systematic cross checks. Tests presented here show that the simulations reproduce sufficiently well important characteristics of real spectra. These mock data-sets will be made available together with the data at the time of the Data Release 11.


title: 24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs

abstract: We have simulated, for the first time, the long term evolution of the Milky Way Galaxy using 51 billion particles on the Swiss Piz Daint supercomputer with our $N$-body gravitational tree-code Bonsai. Herein, we describe the scientific motivation and numerical algorithms. The Milky Way model was simulated for 6 billion years, during which the bar structure and spiral arms were fully formed. This improves upon previous simulations by using 1000 times more particles, and provides a wealth of new data that can be directly compared with observations. We also report the scalability on both the Swiss Piz Daint and the US ORNL Titan. On Piz Daint the parallel efficiency of Bonsai was above 95%. The highest performance was achieved with a 242 billion particle Milky Way model using 18600 GPUs on Titan, thereby reaching a sustained GPU and application performance of 33.49 Pflops and 24.77 Pflops respectively.


title: Performance characterization of a broadband vector Apodizing Phase Plate coronagraph

abstract: One of the main challenges for the direct imaging of planets around nearby stars is the suppression of the diffracted halo from the primary star. Coronagraphs are angular filters that suppress this diffracted halo. The Apodizing Phase Plate coronagraph modifies the pupil-plane phase with an anti-symmetric pattern to suppress diffraction over a 180 degree region from 2 to 7 {\lambda}/D and achieves a mean raw contrast of 10^-4 in this area, independent of the tip-tilt stability of the system. Current APP coronagraphs implemented using classical phase techniques are limited in bandwidth and suppression region geometry (i.e. only on 1 side of the star). In this paper, we show the vector-APP (vAPP) whose phase pattern is implemented by the orientation of patterned liquid crystals. Beam-splitting according to circular polarization states produces two, complementary PSFs with dark holes on either side. We have developed a prototype vAPP that consists of a stack of 3 twisting liquid crystal layers with a bandwidth of 500-900 nm. We characterize the properties of this device using reconstructions of the pupil-plane pattern, and of the ensuing PSF structures. By imaging the pupil between crossed and parallel polarizers we reconstruct the fast axis pattern, transmission, and retardance of the vAPP, and use this as input for a PSF model. This model includes aberrations of the laboratory set-up, and matches the measured PSF, which shows a raw contrast of 10^-3.8 between 2 and 7 {\lambda}/D in a 135 degree wedge. The vAPP coronagraph is relatively easy to manufacture and can be implemented together with a broadband quarter-wave plate and Wollaston prism in a pupil wheel in high-contrast imaging instruments. The manufacturing techniques permit the application of phase patterns with deeper contrasts inside the dark holes and enables unprecedented spectral bandwidths for phase-manipulation coronagraphy.


title: Cosmological simulations of galaxy formation with cosmic rays

abstract: We investigate the dynamical impact of cosmic rays in cosmological simulations of galaxy formation using adaptive-mesh refinement simulations of a $10^{12}$ solar mass halo. In agreement with previous work, a run with only our standard thermal energy feedback model results in a massive spheroid and unrealistically peaked rotation curves. However, the addition of a simple two-fluid model for cosmic rays drastically changes the morphology of the forming disk. We include an isotropic diffusive term and a source term tied to star formation due to (unresolved) supernova-driven shocks. Over a wide range of diffusion coefficients, the CRs generate thin, extended disks with a significantly more realistic (although still not flat) rotation curve. We find that the diffusion of CRs is key to this process, as they escape dense star forming clumps and drive outflows within the more diffuse ISM.


title: The physical nature of the cosmic accretion of baryons and dark matter into halos and their galaxies

abstract: The cosmic accretion of both dark matter and baryons into halos typically is measured using some evolving virial relation, but recent work suggests that most halo growth at late cosmic time (z < 2) is not physical but is rather the by-product of an evolving virial radius ("pseudo-evolution"). Using Omega25, a suite of cosmological simulations that incorporate both dark matter and gas dynamics with differing treatments of gas cooling, star formation, and thermal feedback, we systematically explore the physics that governs cosmic accretion into halos and their galaxies. Physically meaningful cosmic accretion of both dark matter and baryons occurs at z > 1 across our halo mass range: M_200m = 10^{11-14} M_sun. However, dark matter, because it is dissipationless, is deposited (in a time-average sense) at > R_200m(z) in a shell-like manner, such that dark-matter mass and density experience little-to-no physical growth at any radius within a halo at z < 1. In contrast, gas, because it is able to cool radiatively, experiences significant accretion at all radii, at a rate that roughly tracks the accretion rate at R_200m, at all redshifts. Infalling gas starts to decouple from dark matter at ~2 R_200m and continues to accrete to smaller radii until the onset of strong angular-momentum support at ~ 0.1 R_200m. Thus, while the growth of dark matter is subject to pseudo-evolution, the growth of baryons is not. The fact that the accretion rate of gas on galactic scales tracks the accretion rate near R_200m provides insight into the tight relations between the masses/sizes of galaxies and those of their host halos across cosmic time.


title: High-precision Photometric Redshifts from Spitzer/IRAC: Extreme [3.6]-[4.5] Colors Identify Galaxies in the Redshift Range z~6.6-6.9

abstract: One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.


title: The dynamical evolution of molecular clouds near the Galactic Centre - I. Orbital structure and evolutionary timeline

abstract: We recently proposed that the star-forming potential of dense molecular clouds in the Central Molecular Zone (CMZ, i.e. the central few 100 pc) of the Milky Way is linked to their orbital dynamics, potentially giving rise to an absolute-time sequence of star-forming clouds. In this paper, we present an orbital model for the gas stream(s) observed in the CMZ. The model is obtained by integrating orbits in the observed gravitational potential and represents a good fit to the distribution of dense gas, reproducing all of its key properties. The orbit is also consistent with observational constraints not included in the fitting process, such as the velocities of Sgr B2 and the Arches and Quintuplet clusters. It differs from previous models: (1) the orbit is open rather than closed due to the extended mass distribution in the CMZ, (2) its orbital velocity is twice as high as in previous models, and (3) Sgr A$^*$ coincides with the focus of the (eccentric) orbit rather than being offset. Our orbital solution supports the scenario in which the dust ridge between G0.253+0.016 ('the Brick') and Sgr B2 represents an absolute-time sequence of star-forming clouds, triggered by the tidal compression during their recent pericentre passage. We position the clouds on a common timeline and find that their pericentre passages occurred 0.30-0.74 Myr ago. Given their short free-fall times (0.3-0.4 Myr), the quiescent cloud G0.253+0.016 and the vigorously star-forming complex Sgr B2 are separated by a single free-fall time of evolution, implying that star formation proceeds rapidly once collapse has been initiated. We provide several quantitative predictions of our model and conclude with a discussion of the model in the Galactic context, highlighting its relation to large-scale gas accretion, the dynamics of the bar, the $x_2$ orbital family, and the origin of the Arches and Quintuplet clusters. (Abridged)


title: Gravitational Waves and the Scale of Inflation

abstract: We revisit alternative mechanisms of gravitational wave production during inflation and argue that they generically emit a non-negligible amount of scalar fluctuations. We find the scalar power is larger than the tensor power by a factor of order $1/\epsilon^2$. For an appreciable tensor contribution the associated scalar emission completely dominates the zero-point fluctuations of inflaton, resulting in a tensor-to-scalar ratio $r\sim \epsilon^2$. A more quantitative result can be obtained if one further assumes that gravitational waves are emitted by localized sub-horizon processes, giving $r_{\rm max} \simeq 0.3 \epsilon^2$. However, $\epsilon$ is generally time dependent, and this result for $r$ depends on its instantaneous value during the production of the sources, rather than just its average value, somewhat relaxing constraints from the tilt $n_s$. We calculate the scalar 3-point correlation function in the same class of models and show that non-Gaussianity cannot be made arbitrarily small, i.e. $f_{NL} \geq 1$, independently of the value of $r$. Possible exceptions in multi-field scenarios are discussed.


title: Algebraic Geometrization of the Kuramoto Model: Equilibria and Stability Analysis

abstract: Finding equilibria of the finite size Kuramoto model amounts to solving a nonlinear system of equations, which is an important yet challenging problem. We translate this into an algebraic geometry problem and use numerical methods to find all of the equilibria for various choices of coupling constants K, natural frequencies, and on different graphs. We note that for even modest sizes (N ~ 10-20), the number of equilibria is already more than 100,000. We analyze the stability of each computed equilibrium as well as the configuration of angles. Our exploration of the equilibrium landscape leads to unexpected and possibly surprising results including non-monotonicity in the number of equilibria, a predictable pattern in the indices of equilibria, counter-examples to popular conjectures, multi-stable equilibrium landscapes, scenarios with only unstable equilibria, and multiple distinct extrema in the stable equilibrium distribution as a function of the number of cycles in the graph.

categories: nlin.CD cond-mat.stat-mech math-ph math.AG math.DS math.MP


title: Dynamical determination of the gravitational degrees of freedom

abstract: $[n+1]$-dimensional ($n\geq 3$) smooth Einsteinian spaces of Euclidean and Lorentzian signature are considered. The base manifold $M$ is supposed to be smoothly foliated by a two-parameter family of codimension-two-surfaces which are orientable and compact without boundary in $M$. By applying a pair of nested $1+n$ and $1+[n-1]$ decompositions, the canonical form of the metric and the conformal structure of the foliating codimension-two-surfaces a gauge fixing, analogous to the one applied in arXiv:1409.4914, is introduced. In verifying that the true degrees of freedom of gravity may conveniently be represented by the conformal structure it is shown first that regardless whether the primary space is Riemannian or Lorentzian, in terms of the chosen geometrically distinguished new variables, the $1+n$ momentum constraint can be written as a first order symmetric hyperbolic system. It is also argued that in the generic case the Hamiltonian constraint can be solved as an algebraic equation. By combining the $1+n$ constraints with the part of the reduced system of the secondary $1+[n-1]$ decomposition that governs the evolution of the conformal structure---in the Riemannian case with an additional introduction of an imaginary `time'---a well-posed mixed hyperbolic-algebraic-hyperbolic system is formed. It is shown that regardless whether the primary space is of Riemannian or Lorentzian if a regular origin exists solutions to this system are also solutions to the full set of Einstein's equations. This, in particular, offers the possibility of developing a new method for solving Einstein's equations in the Riemannian case. The true degrees of freedom of gravity are also found to be subject of a nonlinear wave equation.


title: An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: Near-infrared morphologies and stellar sizes

abstract: We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$\mu$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $\mu_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a S\'ersic fit to the $H_{160}$ surface-brightness profiles we derive a median S\'ersic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$\mu$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the S\'ersic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.


title: Reality of the quantum state: Towards a stronger {\psi}-ontology theorem

abstract: The Pusey-Barrett-Rudolph (PBR) no-go theorem provides an argument for the reality of the quantum state by ruling out {\psi}-epistemic ontological theories, in which the quantum state is of a statistical nature. It applies under an assumption of preparation independence, the validity of which has been subject to debate. We propose two plausible and less restrictive alternatives: a weaker notion allowing for classical correlations, and an even weaker, physically motivated notion of independence, which merely prohibits the possibility of superluminal causal influences in the preparation process. The latter is a minimal requirement for enabling a reasonable treatment of subsystems in any theory. It is demonstrated by means of an explicit {\psi}-epistemic ontological model that the argument of PBR becomes invalid under the alternative notions of independence. As an intermediate step, we recover a result which is valid in the presence of classical correlations. Finally, we obtain a theorem which holds under the minimal requirement, approximating the result of PBR. For this, we consider experiments involving randomly sampled preparations and derive bounds on the degree of {\psi} epistemicity that is consistent with the quantum-mechanical predictions. The approximation is exact in the limit as the sample space of preparations becomes infinite.


title: Supernova deleptonization asymmetry: Impact on self-induced flavor conversion

abstract: During the accretion phase of a core-collapse supernova (SN), the deleptonization flux has recently been found to develop a global dipole pattern (LESA---Lepton Emission Self-sustained Asymmetry). The $\nu_e$ minus $\bar\nu_e$ flux essentially vanishes in one direction, potentially facilitating self-induced flavor conversion. On the other hand, below the stalled shock wave, self-induced flavor conversion is typically suppressed by multi-angle matter effects, preventing any impact of flavor conversion on SN explosion dynamics. In a schematic model of SN neutrino fluxes, we study the impact of modified $\bar\nu_e$-$\nu_e$ flux asymmetries on collective flavor conversion. In the parameter space consisting of matter density and effective neutrino density, the region of instability with regard to self-induced flavor conversion is much larger for a vanishing lepton number flux, yet this modification does not intersect a realistic SN profile. Therefore, it appears that, even in the presence of LESA, self-induced flavor conversion remains suppressed below the shock front.


title: High-Energy Limit of Quantum Electrodynamics beyond Sudakov Approximation

abstract: We study the high-energy behavior of the scattering amplitudes in quantum electrodynamics beyond the leading order of the small electron mass expansion in the leading logarithmic approximation. In contrast to the Sudakov logarithms, the mass-suppressed double-logarithmic radiative corrections are induced by a soft electron pair exchange and result in enhancement of the power-suppressed contribution. Possible applications of our result to the analysis of the high-energy processes in quantum chromodynamics is also discussed.


title: Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks

abstract: We study diffusive shock acceleration (DSA) of protons and electrons at nonrelativistic, high Mach number, quasiparallel, collisionless shocks by means of self-consistent 1D particle-in-cell simulations. For the first time, both species are found to develop power-law distributions with the universal spectral index $-4$ in momentum space, in agreement with the prediction of DSA. We find that scattering of both protons and electrons is mediated by right-handed circularly polarized waves excited by the current of energetic protons via non-resonant hybrid (Bell) instability. Protons are injected into DSA after a few gyro-cycles of shock drift acceleration (SDA), while electrons are first pre-heated via SDA, then energized via a hybrid acceleration process that involves both SDA and Fermi-like acceleration mediated by Bell waves, before eventual injection into DSA. Using the simulations we can measure the electron/proton ratio in accelerated particles, which is of paramount importance for explaining the cosmic ray fluxes measured at Earth and the multi-wavelength emission of astrophysical objects such as supernova remnants, radio supernovae, and galaxy clusters. We find the normalization of electron power-law is $\lesssim 10^{-2}$ that of the protons for strong nonrelativistic shocks.


title: Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell simulations of oblique pulsars

abstract: We present first-principles relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron-positron pairs. We inject secondary pairs at locations of primary energetic particles, whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere, with the Y-point and current sheet. Solutions with obliquities $\lt 40^{\circ}$ do not show pair production in the open field line region, because the local current density along magnetic field is below the Goldreich-Julian value. The bulk outflow in these solutions is charge separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved space-time or multipolar surface fields.


title: Nature or nurture of coplanar Tatooines: the aligned circumbinary Kuiper belt analogue around HD 131511

abstract: A key discovery of the Kepler mission is of the circumbinary planets known as "Tatooines", which appear to be well aligned with their host stars' orbits. Whether this alignment is due to initially coplanar circumbinary planet-forming discs (i.e. nature), or subsequent alignment of initially misaligned discs by warping the inner disc or torquing the binary (i.e. nurture), is not known. Tests of which scenario dominates may be possible by observing circumbinary Kuiper belt analogues ("debris discs"), which trace the plane of the primordial disc. Here, the 140 au diameter circumbinary debris disc around HD 131511 is shown to be aligned to within 10$^\circ$ of the plane of the near edge-on 0.2 au binary orbit. The stellar equator is also consistent with being in this plane. If the primordial disc was massive enough to pull the binary into alignment, this outcome should be common and distinguishing nature versus nurture will be difficult. However, if only the inner disc becomes aligned with the binary, the HD 131511 system was never significantly misaligned. Given an initial misalignment, the $\sim$Gyr main-sequence lifetime of the star allows secular perturbations to align the debris disc out to 100 au at the cost of an increased scale height. The observed debris disc scale height limits any misalignment to less than 25$^\circ$. With only a handful known, many more such systems need to be characterised to help test whether the alignment of circumbinary planets is nature or nurture.


title: Exo-zodi modelling for the Large Binocular Telescope Interferometer

abstract: Habitable zone dust levels are a key unknown that must be understood to ensure the success of future space missions to image Earth analogues around nearby stars. Current detection limits are several orders of magnitude above the level of the Solar System's Zodiacal cloud, so characterisation of the brightness distribution of exo-zodi down to much fainter levels is needed. To this end, the large Binocular Telescope Interferometer (LBTI) will detect thermal emission from habitable zone exo-zodi a few times brighter than Solar System levels. Here we present a modelling framework for interpreting LBTI observations, which yields dust levels from detections and upper limits that are then converted into predictions and upper limits for the scattered light surface brightness. We apply this model to the HOSTS survey sample of nearby stars; assuming a null depth uncertainty of 10$^{-4}$ the LBTI will be sensitive to dust a few times above the Solar System level around Sun-like stars, and to even lower dust levels for more massive stars.


title: Novel electronic behavior driving NdNiO3 metal-insulator transition

abstract: We present evidence that the metal-insulator transition (MIT) in a tensile strained NdNiO3 (NNO) film is facilitated by a redistribution of electronic density and neither requires Ni charge disproportionation nor symmetry change [1, 2]. Given epitaxial tensile strain in thin NNO films induces preferential occupancy of the $e_g$ $d_{x^2-y^2}$ orbital ($s_{3z^2-r^2}$) we propose the larger transfer integral of this orbital state with the O 2p mediates a redistribution of electronic density from the Ni atom. A decrease in Ni $d_{x^2-y^2}$ orbital occupation is directly observed by resonant inelastic x-ray scattering below the MIT temperature. Furthermore, an increase in Nd charge occupancy is measured by x-ray absorption at the Nd L3 edge. Both spin-orbit coupling and crystal field effects combine to break the degeneracy of the Nd 5d states shifting the energy of the Nd $e_g$ $d_{x^2-y^2}$ orbital towards the Fermi level allowing the A site to become an active acceptor during the MI transition. This work identifies the relocation of electrons from the Ni 3d to the Nd 5d orbitals across the MIT. We propose the insulating gap opens between the Ni 3d and O 2p resulting from Ni 3d electron localization mediated by charge loss. The transition seems neither purely Mott-Hubbard nor simple charge transfer.


title: General procedure for determining braiding and statistics of anyons using entanglement interferometry

abstract: Recently, it was argued that the braiding and statistics of anyons in a two-dimensional topological phase can be extracted by studying the quantum entanglement of the degenerate ground-states on the torus. This construction either required a lattice symmetry (such as $\pi/2$ rotation) or tacitly assumed that the `minimum entanglement states' (MESs) for two different bipartitions can be uniquely assigned quasiparticle labels. Here we describe a procedure to obtain the modular $\mathcal S$ matrix, which encodes the braiding statistics of anyons, which does not require making any of these assumptions. Our strategy is to compare MESs of three independent entanglement bipartitions of the torus, which leads to a unique modular $\mathcal S$. This procedure also puts strong constraints on the modular $\mathcal T$ and $\mathcal U$ matrices without requiring any symmetries, and in certain special cases, completely determines it. Our method applies equally to Abelian and non-Abelian topological phases.


title: Implications of the scalar tilt for the tensor-to-scalar ratio

abstract: We investigate the possible implications of the measured value of the scalar tilt $n_s$ for the tensor-to-scalar ratio $r$ in slow-roll, single-field inflationary models. The measured value of the tilt satisfies $n_s -1\sim 1/N_*$, where $N_* \sim 60$ is the number of $e$-folds for observationally relevant scales. If this is not a coincidence and the scaling holds for different values of $N$, it strongly suggests that either $r$ is as big as $10^{-1}$ (a possibility in tension with the recent data), or smaller than $10^{-2}$ and exponentially dependent on $n_s$. A large region of the ($n_s$,$r$) plane is not compatible with this scaling.


title: Head-to-head domain walls in one-dimensional nanostructures: an extended phase diagram ranging from strips to cylindrical wires

abstract: So far magnetic domain walls in one-dimensional structures have been described theoretically only in the cases of flat strips, or cylindrical structures with a compact cross-section, either square or disk. Here we describe an extended phase diagram unifying the two pictures, extensively covering the (width,thickness) space. It is derived on the basis of symmetry and phase-transition arguments, and micromagnetic simulations. A simple classification of all domain walls in two varieties is proposed on the basis of their topology: either with a combined transverse/vortex character, or of the Bloch-point type. The exact arrangement of magnetization within each variety results mostly from the need to decrease dipolar energy, giving rise to asymmetric and curling structures. Numerical evaluators are introduced to quantify curling, and scaling laws are derived analytically for some of the iso-energy lines of the phase diagram.


title: Fast Sublinear Sparse Representation using Shallow Tree Matching Pursuit

abstract: Sparse approximations using highly over-complete dictionaries is a state-of-the-art tool for many imaging applications including denoising, super-resolution, compressive sensing, light-field analysis, and object recognition. Unfortunately, the applicability of such methods is severely hampered by the computational burden of sparse approximation: these algorithms are linear or super-linear in both the data dimensionality and size of the dictionary. We propose a framework for learning the hierarchical structure of over-complete dictionaries that enables fast computation of sparse representations. Our method builds on tree-based strategies for nearest neighbor matching, and presents domain-specific enhancements that are highly efficient for the analysis of image patches. Contrary to most popular methods for building spatial data structures, out methods rely on shallow, balanced trees with relatively few layers. We show an extensive array of experiments on several applications such as image denoising/superresolution, compressive video/light-field sensing where we practically achieve 100-1000x speedup (with a less than 1dB loss in accuracy).


title: Near Optimal LP Rounding Algorithm for Correlation Clustering on Complete and Complete k-partite Graphs

abstract: We give new rounding schemes for the standard linear programming relaxation of the correlation clustering problem, achieving approximation factors almost matching the integrality gaps: - For complete graphs our appoximation is $2.06 - \varepsilon$ for a fixed constant $\varepsilon$, which almost matches the previously known integrality gap of $2$. - For complete $k$-partite graphs our approximation is $3$. We also show a matching integrality gap. - For complete graphs with edge weights satisfying triangle inequalities and probability constraints, our approximation is $1.5$, and we show an integrality gap of $1.2$. Our results improve a long line of work on approximation algorithms for correlation clustering in complete graphs, previously culminating in a ratio of $2.5$ for the complete case by Ailon, Charikar and Newman (JACM'08). In the weighted complete case satisfying triangle inequalities and probability constraints, the same authors give a $2$-approximation; for the bipartite case, Ailon, Avigdor-Elgrabli, Liberty and van Zuylen give a $4$-approximation (SICOMP'12).


title: Circumbinary Components of Contact Binaries

abstract: The period changes of contact binaries obtained by the analysis of eclipse minima timing are found mostly chaotic in nature. However, they are representable by a few cyclic changes superposed on a secular change. The cyclic changes are caused most probably by the third components revolving around the contact binaries. Some typical examples of the period changes of contact binaries are presented in the present contribution.


title: Robust Stability Assessment in the Presence of Load Dynamics Uncertainty

abstract: Dynamic response of loads has a significant effect on system stability and directly determines the stability margin of the operating point. Inherent uncertainty and natural variability of load models make the stability assessment especially difficult and may compromise the security of the system. We propose a novel mathematical "robust stability" criterion for the assessment of small-signal stability of operating points. Whenever the criterion is satisfied for a given operating point, it provides mathematical guarantees that the operating point will be stable with respect to small disturbances for any dynamic response of the loads. The criterion can be naturally used for identification of operating regions secure from the occurrence of Hopf bifurcation. Several possible applications of the criterion are discussed, most importantly the concept of Robust Stability Assessment (RSA) that could be integrated in dynamic security assessment packages and used in contingency screening and other planning and operational studies.


title: A Global Identifiability Condition for Consensus Networks with Tree Graphs

abstract: In this paper we present a sufficient condition that guarantees identifiability of linear network dynamic systems exhibiting continuous-time weighted consensus protocols with acyclic structure. Each edge of the underlying network graph $\mathcal G$ of the system is defined by a constant parameter, referred to as the weight of the edge, while each node is defined by a scalar state whose dynamics evolve as the weighted linear combination of its difference with the states of its neighboring nodes. Following the classical definitions of identifiability and indistinguishability, we first derive a condition that ensure the identifiability of the edge weights of $\mathcal G$ in terms of the associated transfer function. Using this characterization, we propose a sensor placement algorithm that guarantees identifiability of the edge weights. We describe our results using several illustrative examples.


title: Positive scalar curvature and product formulas for secondary index invariants

abstract: We introduce partial secondary invariants associated to complete Riemannian metrics which have uniformly positive scalar curvature outside a prescribed subset on a spin manifold. These can be used to distinguish such Riemannian metrics up to concordance relative to the prescribed subset. We exhibit a general external product formula for partial secondary invariants, from which we deduce product formulas for the higher rho-invariant of a metric with uniformly positive scalar curvature as well as for the higher relative index of two metrics with uniformly positive scalar curvature. Our methods yield a new conceptual proof of the secondary partitioned manifold index theorem and a refined version of the delocalized APS-index theorem of Piazza-Schick for the spinor Dirac operator in all dimensions. We establish a partitioned manifold index theorem for the higher relative index. We also show that secondary invariants are stable with respect to direct products with aspherical manifolds that have fundamental groups of finite asymptotic dimension. Moreover, we construct examples of complete metrics with uniformly positive scalar curvature on non-compact spin manifolds which can be distinguished up to concordance relative to subsets which are coarsely negligible in a certain sense. A technical novelty in this paper is that we use Yu's localization algebras in combination with the description of K-theory for graded C*-algebras due to Trout. This formalism allows direct definitions of all the invariants we consider in terms of the functional calculus of the Dirac operator and enables us to give concise proofs of the product formulas.


title: Practical variational tomography for critical 1D systems

abstract: We improve upon a recently introduced efficient quantum state reconstruction procedure targeted to states well-approximated by the multi-scale entanglement renormalization ansatz (MERA), e.g., ground states of critical models. We show how to numerically select a subset of experimentally accessible measurements which maximizes information extraction about renormalized particles, thus dramatically reducing the required number of physical measurements. We numerically estimate the number of measurements required to characterize the ground state of the critical 1D Ising (resp. XX) model and find that MERA tomography on 16-qubit (resp. 24-qubit) systems requires the same experimental effort than brute-force tomography on 8 qubits. We derive a bound computable from experimental data which certifies the distance between the experimental and reconstructed states.


title: Entanglement entropy in three dimensional gravity

abstract: The Ryu-Takayanagi and covariant Hubeny-Rangamani-Takayanagi proposals relate entanglement entropy in CFTs with holographic duals to the areas of minimal or extremal surfaces in the bulk geometry. We show how, in three dimensional pure gravity, the relevant regulated geodesic lengths can be obtained by writing a spacetime as a quotient of AdS3, with the problem reduced to a simple purely algebraic calculation. We explain how this works in both Lorentzian and Euclidean formalisms, before illustrating its use to obtain novel results in a number of examples, including rotating BTZ, the RP2 geon, and several wormhole geometries. This includes spatial and temporal dependence of single-interval entanglement entropy, despite these symmetries being broken only behind an event horizon. We also discuss considerations allowing HRT to be derived from analytic continuation of Euclidean computations in certain contexts, and a related class of complexified extremal surfaces.


title: Dissipation due to pure spin-current generated by spin pumping

abstract: Based on spin-dependent transport theory and thermodynamics, we develop a generalized theory of the Joule heating in the presence of a spin current. Along with the conventional Joule heating consisting of an electric current and electrochemical potential, it is found that the spin current and spin accumulation give an additional dissipation because the spin-dependent scatterings inside bulk and ferromagnetic/nonmagnetic interface lead to a change of entropy. The theory is applied to investigate the dissipation due to pure spin-current generated by spin pumping across a ferromagnetic/nonmagnetic/ferromagnetic multilayer. The dissipation arises from an interface because the spin pumping is a transfer of both the spin angular momentum and the energy from the ferromagnet to conduction electrons near the interface. It is found that the dissipation is proportional to the enhancement of the Gilbert damping constant by spin pumping.