title: Additional Bending of Light in Sun's Vicinity by its Interior Index of Refraction

abstract: In the seventies, scientists observed discrepancies of the bending of light around the Sun based on Einstein's prediction of the curvature of star light due to the mass of the Sun. We claim that the interior electromagnetic properties of the Sun influence the curvature of the light path outside the Sun as well. In this paper, we investigate the additional deflection of light in the vacuum region surrounding the Sun by its electromagnetic parameters. Starting with Maxwell's equations, we show how the deflection of light passing the Sun depends on the electric permittivity and the magnetic permeability of the interior of the Sun. The electromagnetic field equations in Cartesian coordinates are transformed to the ones in an appropriately chosen Riemannian space. This coordinate transform is dictated by the introduction of a refractional potential. The geodetic lines with the shortest propagation time are constructed from this potential. As far as the deflection of light propagating along these geodetic lines is concerned, we show that the existence of a refractional potential influences the light path outside any object with a typical refractive index. Our results add new aspects to the bending of star light explained by general relativity. Some astrophysical observations, which cannot be explained by gravity in a satisfactory manner, are justified by the electromagnetic model. In particular, the frequency dependency of the light deflection is discussed. We show that the additional bending due to the refractive index is proportional to the third power of the inverse distance. The general relativity predicts that the bending due to the mass is proportional to the inverse distance.


title: Strong approximation in h-mass of rectifiable currents under homological constraint

abstract: Let h : R $\rightarrow$ R+ be a lower semi-continuous subbadditive and even function such that h(0) = 0 and h($\theta$) $\ge$ $\alpha$|$\theta$| for some $\alpha$ > 0. The h-mass of a k-polyhedral chain P =$\sum$j $\theta$j$\sigma$j in R n (0 $\le$ k $\le$ n) is defined as M h (P) := j h($\theta$j) H k ($\sigma$j). If T = $\tau$ (M, $\theta$, $\xi$) is a k-rectifiable chain, the definition extends to M h (T) := M h($\theta$) dH k. Given such a rectifiable flat chain T with M h (T) < $\infty$ and $\partial$T polyhedral, we prove that for every $\eta$ > 0, it decomposes as T = P + $\partial$V with P polyhedral, V rectifiable, M h (V) < $\eta$ and M h (P) < M h (T) + $\eta$. In short, we have a polyhedral chain P which strongly approximates T in h-mass and preserves the homological constraint $\partial$P = $\partial$T. These results are motivated by the study of approximations of M h by smoother functionals but they also provide explicit formulas for the lower semicontinuous envelope of T $\rightarrow$ M h (T) + I $\partial$S ($\partial$T) with respect to the topology of the flat norm.


title: L\'evy-Leblond fermions on the wormhole

abstract: We propose a simple model of entanglement generated by geometry, studying non-relativistic massive L\'evy-Leblond fermions in the $1+2$ geometry of a Bronnikov-Ellis wormhole. The model is equivalent to that of relativistic massless Dirac fermions in $1+3$ dimensions, where one spatial direction is flat. The effect of the wormhole is to generate quantum states that, far from the throat, are approximated by entangled particles on two flat, separated spacetime regions. An appealing feature of the model is that it has a condensed matter analogue, the regime of intermediate energies for two planes of bilayer graphene linked by a bilayer carbon nanotube. Therefore we expect that it might be possible to realize in the laboratory the entangled states studied here. We argue that generalisations of our solvable model which preserve the topology will have similar quantum behaviour.


title: Correlations of strongly interacting ultracold dipolar bosons in optical lattices

abstract: Strongly interacting dipolar bosons in optical lattices exhibit diverse quantum phases that are rich in physics. As the strength of the long-range boson-boson interaction increases, the system transitions across different phases: from a superfluid, through a Mott-insulator to a Tonks gas and, eventually, a crystal state. The signature of these phases and their transitions can be unequivocally identified by an experimentally detectable order parameter, recently described in arXiv:1708.07409. Herein, we calculate the momentum distributions and the normalized Glauber correlation functions of dipolar bosons in a one-dimensional optical lattice in order to characterize all their phases. To understand the behavior of the correlations across the phase transitions, we first investigate the eigenfunctions and eigenvalues of the one-body reduced density matrix as the function of the dipolar interaction strength. We then analyze the real- and momentum-space Glauber correlation functions, thereby gaining a spatially and momentum-resolved insight into the coherence properties of these quantum phases. We find an intriguing structure of non-local correlations that, independently of other system parameters, reveal the phase transitions of the system. In particular, spatial localization with synchronous momentum delocalization accompanies the formation of correlated islands in the density. Moreover, our study showcases that precise control of intersite correlations is possible through manipulation of the depth of the lattice, while intrasite correlations are influenced solely by changing the dipolar interaction strength.


title: MAP inference via Block-Coordinate Frank-Wolfe Algorithm

abstract: We present a new proximal bundle method for Maximum-A-Posteriori (MAP) inference in structured energy minimization problems. The method optimizes a Lagrangean relaxation of the original energy minimization problem using a multi plane block-coordinate Frank-Wolfe method that takes advantage of the specific structure of the Lagrangean decomposition. We show empirically that our method outperforms state-of-the-art Lagrangean decomposition based algorithms on some challenging Markov Random Field, multi-label discrete tomography and graph matching problems.


title: Constraints on the Thermal Contents of the X-ray Cavities of Cluster MS 0735.6+7421 with Sunyaev-Zel'dovich Effect Observations

abstract: Outbursts from active galactic nuclei (AGN) can inflate cavities in the intracluster medium (ICM) of galaxy clusters and are believed to play the primary role in offsetting radiative cooling in the ICM. However, the details of how the energy from AGN feedback thermalizes to heat the ICM is not well understood, partly due to the unknown composition and energetics of the cavities. The Sunyaev-Zel'dovich (SZ) effect, a measure of the integrated pressure along the line of sight, provides a means of measuring the thermal contents of the cavities, to discriminate between thermal, nonthermal, and other sources of pressure support. Here we report measurements of the SZ effect at 30 GHz toward the galaxy cluster MS 0735.6+7421 (MS0735), using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). MS0735 hosts the most energetic AGN outburst known and lobes of radio synchrotron emission coincident with a pair of giant X-ray cavities $\sim 200$ across. Our CARMA maps show a clear deficit in the SZ signal coincident with the X-ray identified cavities, when compared to a smooth X-ray derived pressure model. We find that the cavities have very little SZ-contributing material, suggesting that they are either supported by very diffuse thermal plasma with temperature in excess of hundreds of keV, or are not supported thermally. Our results represent the first detection (with $4.4 \sigma$ significance) of this phenomenon with the SZ effect.


title: Solvation in the Large Box Limit

abstract: In this paper, the authors study the limit of a sharp interface model for the solvation of charged molecules in an implicit solvent as the number of solute molecules and the size of the surrounding box tend to infinity. The energy is given by a combination of local terms accounting for the physical presence of the molecules in the solvent and a nonlocal electrical energy with or without an ionic effect. In the presence of an ionic effect, the authors prove a screening effect in the limit, i.e., the limit is completely localized and hence electrical long-range interactions of the molecules can be neglected. In the absence of the ionic effect, the authors show that the behavior of the energy depends on the scaling of the number of molecules with respect to the size of the surrounding box. All scaling regimes are identified and corresponding limit results proved. In regimes with many solute molecules this limit includes electrical interactions of $H^{-1}$-type between the molecules.


title: Generalized derivatives for the solution operator of the obstacle problem

abstract: We characterize generalized derivatives of the solution operator of the obstacle problem. This precise characterization requires the usage of the theory of so-called capacitary measures and the associated solution operators of relaxed Dirichlet problems. The generalized derivatives can be used to obtain a novel necessary optimality condition for the optimal control of the obstacle problem with control constraints. A comparison shows that this system is stronger than the known system of C-stationarity.


title: Geometric Lagrangian averaged Euler-Boussinesq and primitive equations

abstract: In this article we derive the equations for a rotating stratified fluid governed by inviscid Euler-Boussinesq and primitive equations that account for the effects of the perturbations upon the mean. Our method is based on the concept of geometric generalized Lagrangian mean recently introduced by Gilbert and Vanneste, combined with generalized Taylor and horizontal isotropy of fluctuations as turbulent closure hypotheses. The models we obtain arise as Euler-Poincar\'{e} equations and inherit from their parent systems conservation laws for energy and potential vorticity. They are structurally and geometrically similar to Euler-Boussinesq-$\alpha$ and primitive equations-$\alpha$ models, however feature a different regularizing second order operator.

categories: physics.flu-dyn math-ph math.MP physics.ao-ph


title: Support Vector Machine Application for Multiphase Flow Pattern Prediction

abstract: In this paper a data analytical approach featuring support vector machines (SVM) is employed to train a predictive model over an experimentaldataset, which consists of the most relevant studies for two-phase flow pattern prediction. The database for this study consists of flow patterns or flow regimes in gas-liquid two-phase flow. The term flow pattern refers to the geometrical configuration of the gas and liquid phases in the pipe. When gas and liquid flow simultaneously in a pipe, the two phases can distribute themselves in a variety of flow configurations. Gas-liquid two-phase flow occurs ubiquitously in various major industrial fields: petroleum, chemical, nuclear, and geothermal industries. The flow configurations differ from each other in the spatial distribution of the interface, resulting in different flow characteristics. Experimental results obtained by applying the presented methodology to different combinations of flow patterns demonstrate that the proposed approach is state-of-the-art alternatives by achieving 97% correct classification. The results suggest machine learning could be used as an effective tool for automatic detection and classification of gas-liquid flow patterns.


title: Nonuniform average sampling in multiply generated shift-invariant subspaces of mixed Lebesgue spaces

abstract: In this paper, we study nonuniform average sampling problem in multiply generated shift-invariant subspaces of mixed Lebesgue spaces. We discuss two types of average sampled values: average sampled values $\{\left \langle f,\psi_{a}(\cdot-x_{j},\cdot-y_{k}) \right \rangle:j,k\in\mathbb{J}\}$ generated by single averaging function and average sampled values $\left\{\left \langle f,\psi_{x_{j},y_{k}}\right \rangle:j,k\in\mathbb{J}\right\}$ generated by multiple averaging functions. Two fast reconstruction algorithms for this two types of average sampled values are provided.


title: Multi-scale analysis of turbulence evolution in the density stratified intracluster medium

abstract: The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolution of turbulence in the intracluster medium (ICM) after major mergers, using the Omega500 non-radiative hydrodynamical cosmological simulations. By applying a novel wavelet analysis to study the radial dependence of the ICM turbulence spectrum, we discover that faster turbulence dissipation in the inner high density regions leads to the increasing turbulence amplitude with radius. We also find that the ICM turbulence at all radii decays in two phases after a major merger: an early fast decay phase followed by a slow secular decay phase. The buoyancy effects resulting from the ICM density stratification becomes increasingly important during turbulence decay, as revealed by a decreasing turbulence Froude number $Fr \sim \mathcal{O}(1)$. Our results indicate that the stronger density stratification and smaller eddy turn-over time are the likely causes of the faster turbulence dissipation rate in the inner regions of the cluster.


title: Unobservable False Data Injection Attacks against PMUs: Feasible Conditions and Multiplicative Attacks

abstract: This paper studies false data injection (FDI) attacks against phasor measurement units (PMUs). As compared to the conventional bad data detector (BDD), an enhanced BDD utilizing the effect of zero injection buses is proposed. Feasible conditions under which FDI attacks are unobservable to this enhanced BDD are discussed. In addition, a class of multiplicative FDI attacks that maintain the rank of the PMU measurement matrix is introduced. Simulation results on the IEEE RTS-24-bus system indicate that the these multiplicative unobservable attacks can avoid detection by both the enhanced BDD and a detector based on low-rank decomposition proposed in prior work.


title: Prospects of indirect searches for dark matter annihilations in the earth with ICAL@INO

abstract: We study the prospects of detecting muon events at the upcoming Iron CALorimeter (ICAL) detector to be built at the proposed India-based Neutrino Observatory (INO) facility due to neutrinos arising out of annihilation of Weakly Interactive Massive Particles (WIMP) in the centre of the earth. The atmospheric neutrinos coming from the direction of earth core presents an irreducible background. We consider 50kt $\times$ 10 years of ICAL running and WIMP masses between 10-100 GeV and present 90 \% C.L. exclusion sensitivity limits on $\sigma_{SI}$ which is the WIMP-nucleon Spin Independent (SI) interaction cross-section. The expected sensitivity limits calculated for ICAL for the WIMP annihilation in the earth are more stringent than the limits obtained by any other indirect detection experiment. For a WIMP mass of ~$52.14 \textup{ GeV}$, where the signal fluxes are enhanced due to resonance capture of WIMP in earth due to Fe nuclei, the sensitivity limits, assuming 100\% branching ratio for each channel, are : $\sigma_{SI} =1.02\times 10^{-44}~cm^2$ for the $\tau^{+} \tau^{-}$ channel and $\sigma_{SI} =5.36\times 10^{-44} ~cm^2$ for the $b ~\bar{b}$ channel.


title: Multilingual End-to-End Speech Recognition with A Single Transformer on Low-Resource Languages

abstract: Sequence-to-sequence attention-based models integrate an acoustic, pronunciation and language model into a single neural network, which make them very suitable for multilingual automatic speech recognition (ASR). In this paper, we are concerned with multilingual speech recognition on low-resource languages by a single Transformer, one of sequence-to-sequence attention-based models. Sub-words are employed as the multilingual modeling unit without using any pronunciation lexicon. First, we show that a single multilingual ASR Transformer performs well on low-resource languages despite of some language confusion. We then look at incorporating language information into the model by inserting the language symbol at the beginning or at the end of the original sub-words sequence under the condition of language information being known during training. Experiments on CALLHOME datasets demonstrate that the multilingual ASR Transformer with the language symbol at the end performs better and can obtain relatively 10.5\% average word error rate (WER) reduction compared to SHL-MLSTM with residual learning. We go on to show that, assuming the language information being known during training and testing, about relatively 12.4\% average WER reduction can be observed compared to SHL-MLSTM with residual learning through giving the language symbol as the sentence start token.


title: Long-Range Longitudinal Electric Wave in Vacuum Radiated by Electric Dipole: Part II

abstract: In this work we further advance theoretical investigation of radiation by the electric dipole under the assumption that wavelength is much smaller than charge separation distance of an electric dipole, which in turn is much smaller than a distance up to the point of observation. The electric dipole considered in this paper is the one with fixed charge positions, but oscillating charge magnitudes. Specifically, two cases were considered. In the first case phase delay between oscillations of the charge magnitudes was taken into account. In the second case, opposite current wave travel direction was considered. These results generalize the classical solution for electric dipole radiation, and they show that under the above assumptions the electric dipole emits both long-range longitudinal electric and transverse electromagnetic waves. Longitudinal and transverse electric fields emitted by the electric dipole with current wave traveling in both directions are displayed for specific values of the dipole system parameters. It was shown that under the standard assumptions that charge separation distance is much smaller than wavelength and non-zero phase delay between charge oscillations: a) our solution reduces to classical solution; b) transverse electromagnetic waves are dominant while longitudinal electric waves are negligible; c) total radiated power is proportional to the fourth degree of frequency and to the second degree of the charge separation distance. In case wavelength is much smaller than charge separation distance: a) the classical solution is invalid and it overestimates the total radiated power; b) longitudinal electric waves are dominant and transverse electromagnetic waves are negligible; c) total radiated power is proportional to the third degree of frequency and to the charge separation distance. All derived potential and force fields were summarized in the table.


title: Analysis of the image of pion-emitting sources in source center of mass frame

abstract: In this paper, we try a method to extract the image of pion-emitting source function in the center-of-mass frame of source (CMFS). We choose the identical pion pairs according to the difference of their energy and use these pion pairs to build the correlation function. The purpose is to reduce the effect of $\triangle E \triangle t$, thus the corresponding imaging result can tend to the real source function. We examine the effect of this method by comparing its results with real source functions extracted from models directly.


title: Two-grid discretizations and a local finite element scheme for a non-selfadjoint Stekloff eigenvalue problem

abstract: In this paper, for a new Stekloff eigenvalue problem which is non-selfadjoint and not $H^1$-elliptic, we establish and analyze two kinds of two-grid discretization scheme and a local finite element scheme. We present the error estimates of approximations of two-grid discretizations. We also prove a local error estimate which is suitable for the case that the local refined region contains singular points lying on the boundary of domain. Numerical experiments are reported finally to show the efficiency of our schemes.


title: Numerical methods for scattering problems from multi-layers with different periodicities

abstract: In this paper, we consider a numerical method to solve scattering problems with multi-periodic layers with different periodicities. The main tool applied in this paper is the Bloch transform. With this method, the problem is written into an equivalent coupled family of quasi-periodic problems. As the Bloch transform is only defined for one fixed period, the inhomogeneous layer with another period is simply treated as a non-periodic one. First, we approximate the refractive index by a periodic one where its period is an integer multiple of the fixed period, and it is decomposed by finite number of quasi-periodic functions. Then the coupled system is reduced into a simplified formulation. A convergent finite element method is proposed for the numerical solution, and the numerical method has been applied to several numerical experiments. At the end of this paper, relative errors of the numerical solutions will be shown to illustrate the convergence of the numerical algorithm.


title: Peculiar Spin Frequency and Radio Profile Evolution of PSR J1119$-$6127 Following Magnetar-like X-ray Bursts

abstract: We present the spin frequency and profile evolution of the radio pulsar J1119$-$6127 following magnetar-like X-ray bursts from the system in 2016 July. Using data from the Parkes radio telescope, we observe a smooth and fast spin-down process subsequent to the X-ray bursts resulting in a net change in the pulsar rotational frequency of $\Delta\nu\approx-4\times10^{-4}$\,Hz. During the transition, a net spin-down rate increase of $\Delta\dot\nu\approx-1\times10^{-10}$\,Hz\,s$^{-1}$ is observed, followed by a return of $\dot{\nu}$ to its original value. In addition, the radio pulsations disappeared after the X-ray bursts and reappeared about two weeks later with the flux density at 1.4\,GHz increased by a factor of five. The flux density then decreased and undershot the normal flux density followed by a slow recovery back to normal. The pulsar's integrated profile underwent dramatic and short-term changes in total intensity, polarization and position angle. Despite the complex evolution, we observe correlations between the spin-down rate, pulse profile shape and radio flux density. Strong single pulses have been detected after the X-ray bursts with their energy distributions evolving with time. The peculiar but smooth spin frequency evolution of PSR~J1119$-$6127 accompanied by systematic pulse profile and flux density changes are most likely to be a result of either reconfiguration of the surface magnetic fields or particle winds triggered by the X-ray bursts. The recovery of spin-down rate and pulse profile to normal provides us the best case to study the connection between high magnetic-field pulsars and magnetars.


title: Interplanetary radio type II and type IV bursts as indicators of propagating solar transients

abstract: Recent studies of interplanetary radio type II bursts and their source locations are reviewed. As these bursts are due to propagating shock waves, driven by coronal mass ejections, they can be followed to near-Earth distances and can be used to predict the arrival times of geo-effective disturbances. Radio type IV bursts, on the other hand, are usually due to moving magnetic structures in the low corona and trapped particles form at least part of the emission. The observed directivity of type IV emission may also be used for space weather purposes.


title: Perturbation determinants on Banach spaces and operator differentiability for Hirsch functional calculus

abstract: We consider a perturbation determinant for pairs of nonpositive (in a sense of Komatsu) operators on Banach space with nuclear difference and prove a generalization of the important formula for the logarithmic derivative of this determinant. To this end the Frechet differentiability of operator monotonic (negative complete Bernstein) functions of negative and nonpositive operators on Banach spaces is investigated. The results may be regarded as a contribution to the Hirsch functional calculus.


title: Strain-Gradient Plasticity as the $\Gamma$-Limit of a Nonlinear Dislocation Energy with Mixed Growth

abstract: In this paper a we derive by means of $\Gamma$-convergence a macroscopic strain-gradient plasticity from a semi-discrete model for dislocations in an infinite cylindrical crystal. In contrast to existing work, we consider an energy with subquadratic growth close to the dislocations. This allows to treat the stored elastic energy without the need to introduce an ad-hoc cut-off radius. As the main tool to prove a complementing compactness statement, we present a generalized version of the geometric rigidity result for fields with non-vanishing $\operatorname{curl}$. A main ingredient is a fine decomposition result for $L^1$-functions whose divergence is in certain critical Sobolev spaces.


title: Quantum Monte Carlo prerequisites for industrial catalysts. Accurately assessing H atom and H$_2$ physical adsorption energy on Pt(111)

abstract: The yardstick of new first-principles approaches to key points on reaction paths at metal surfaces is chemical accuracy compared to reliable experiment. By this we mean that such values as the activation barrier are required to within 1 kcal/mol. Quantum Monte Carlo (QMC) is a promising (albeit lengthy) first-principles method for this and we are now beyond the dawn of QMC benchmarks for these systems, since hydrogen dissociation on Cu(111) has been studied with quite adequate accuracy in two improving QMC studies \cite{hog0, kdd2}. Pt and Cu require the use of pseudo-potentials in these calculations and we show that those of Pt are less problematic than those for Cu, particularly for QMC work. This work determines physisorption energies for hydrogen atoms and molecules on Pt(111). These systems are used as asymptotes to determine reaction barriers. As such, they must be referred to clean Pt(111) surfaces and the isolated hydrogen atom or molecule. Previous work gave the activation barrier to hydrogen dissociation on Pt(111) using the bridging geometry. The method used is state-of-the-art ({\it ab initio}) QMC. This barrier agrees to better than chemical accuracy with a recent Specific Reaction Parameter (SRP-DFT) calculation, fitted to measurement \cite{iran}. The results give the dissociation barrier for hydrogen on Pt(111) as 5.4 kcal/mol (QMC and 6.2 by SRP-DFT) with a QMC standard error below 0.99 kcal/mol. This is encouraging for establishing less well-known benchmark values of industrial reaction barriers on Pt(111).


title: Minimizing Regret of Bandit Online Optimization in Unconstrained Action Spaces

abstract: We consider online convex optimization with a zero-order oracle feedback. In particular, the decision maker does not know the explicit representation of the time-varying cost functions, or their gradients. At each time step, she observes the value of the corresponding cost function evaluated at her chosen action (zero-order oracle). The objective is to minimize the regret, that is, the difference between the sum of the costs she accumulates and that of a static optimal action had she known the sequence of cost functions a priori. We present a novel algorithm to minimize regret in unconstrained action spaces. Our algorithm hinges on a classical idea of one-point estimation of the gradients of the cost functions based on their observed values. The algorithm is independent of problem parameters. Letting $T$ denote the number of queries of the zero-order oracle and $n$ the problem dimension, the regret rate achieved is $O(n^{2/3}T^{2/3})$. Moreover, we adapt the presented algorithm to the setting with two-point feedback and demonstrate that the adapted procedure achieves the theoretical lower bound on the regret of $(n^{1/2}T^{1/2})$.


title: Explicit estimates of sums related to the Nyman-Beurling criterion for the Riemann Hypothesis

abstract: We give an estimate for sums appearing in the Nyman-Beurling criterion for the Riemann Hypothesis. These sums contain the M\"obius function and are related to the imaginary part of the Estermann zeta function. The estimate is remarkably sharp in comparison to other sums containing the M\"obius function. The bound is smaller than the trivial bound - essentially the number of terms - by a fixed power of that number. The exponent is made explicit. The methods intensively use tools from the theory of continued fractions and from the theory of Fourier series.


title: Tri-comb spectroscopy

abstract: Multi-dimensional coherent spectroscopy (MDCS)1,2 is a powerful method for optical spectroscopy that has become an important tool for studying ultrafast dynamics in a wide range of systems. It is an optical analog of multidimensional nuclear magnetic resonance spectroscopy that enables the measurement of homogeneous linewidths in inhomogeneously broadened systems, many-body interactions, and coupling between excited resonances, all of which are not simultaneously accessible by any other linear or non-linear method. Current implementations of MDCS require a bulky apparatus and suffer from resolution and acquisition speed limitations that constrain their applications outside the laboratory3-5. Here we propose and demonstrate an approach to nonlinear coherent spectroscopy that utilizes three frequency combs with slightly different repetition rates. Unlike traditional nonlinear methods, tri-comb spectroscopy uses only a single photodetector and no mechanical moving elements to enable faster acquisition times, while also providing comb resolution. As a proof of concept, a multidimensional coherent spectrum with comb cross-diagonal resolution is generated using only 365 ms of data. These improvements make multidimensional coherent spectroscopy relevant for systems with narrow resonances (especially cold atomic and molecular systems). In addition the method has the potential to be field deployable for chemical sensing applications.


title: The entrainment matrix of a superfluid nucleon mixture at finite temperatures

abstract: It is considered a closed system of non-linear equations for the entrainment matrix of a non-relativistic mixture of superfluid nucleons at arbitrary temperatures below the onset of neutron superfluidity, which takes into account the essential dependence of the superfluid energy gap in the nucleon spectra on the velocities of superfluid flows. It is assumed that the protons condense into the isotropic $^1$S$_0$ state, and the neutrons are paired into the spin-triplet $^3$P$_2$ state. It is derived an analytic solution to the non-linear equations for the entrainment matrix under temperatures just below the critical value for the neutron superfluidity onset. In general case of an arbitrary temperature of the superfluid mixture the non-linear equations are solved numerically and fitted by simple formulas convenient for a practical use with an arbitrary set of the Landau parameters.


title: Application Of McDiarmid Inequality In Finite-Key-Length Decoy-State Quantum Key Distribution

abstract: In practical decoy-state quantum key distribution, the raw key length is finite. Thus, deviation of the estimated single photon yield and single photon error rate from their respective true values due to finite sample size can seriously lower the provably secure key rate $R$. Current method to obtain a lower bound of $R$ follows an indirect path by first bounding the yields and error rates both conditioned on the type of decoy used. These bounds are then used to deduce the single photon yield and error rate, which in turn are used to calculate a lower bound of the key rate $R$. Here I show how to directly compute a lower bound of $R$ via McDiarmid inequality in statistics. This method increases the provably secure key rate of realistic quantum channels by at least 30% when the raw key length is $\approx 10^5$ to $10^6$. More importantly, this is achieved by pure theoretical analysis without altering the experimental setup or the post-processing method. In a boarder context, this work introduces powerful concentration inequality techniques in statistics to tackle physics problem beyond straightforward statistical data analysis.


title: An extended framework of continuous-stage Runge-Kutta methods

abstract: We propose an extended framework for continuous-stage Runge-Kutta methods which enables us to treat more complicated cases especially for the case weighting on infinite intervals. By doing this, various types of weighted orthogonal polynomials (e.g., Jacobi polynomials, Laguerre polynomials, Hermite polynomials etc.) can be used in the construction of Runge-Kutta-type methods. Particularly, families of Runge-Kutta-type methods with geometric properties can be constructed in this new framework. As examples, some new symplectic integrators by using Legendre polynomials, Laguerre polynomials and Hermite polynomials are constructed.


title: Minimizing hysteresis in martensite phase transforming magnetocaloric Heusler alloys

abstract: The large magnetocaloric effect in Heusler alloys showing martensite phase transformation puts them forward as efficient materials for magnetic refrigeration. However, irreversibility of the magnetocaloric cooling cycle is a major challenge for real applications. This irreversibility is directly linked to the thermal hysteresis at the first-order martensite phase transition. Therefore, minimizing the hysteresis is essential in order to achieve reversibility. Here we show a large reduction in the thermal hysteresis at the martensite transition in the Ni$_{2}$Mn$_{1.4}$In$_{0.6}$ and Ni$_{1.8}$Co$_{0.2}$Mn$_{1.4}$In$_{0.6}$ Heusler alloys upon the application of hydrostatic pressure. Our pressure dependent X-ray diffraction study on Ni$_{2}$Mn$_{1.4}$In$_{0.6}$ reveals that with increasing pressure the lattice parameters of the two crystallographic phases (austenite and martensite) change in such a way that they increasingly satisfy the geometric compatibility (co-factor) condition. These results provide an opportunity to overcome the hysteresis problem and hence the irreversible behavior in Heusler materials using pressure as an external parameter.


title: The massive Feynman propagator on asymptotically Minkowski spacetimes II

abstract: We consider the massive Klein-Gordon equation on short-range asymptotically Minkowski spacetimes. Extending our results in [GW1], we show that the Klein-Gordon operator with Feynman type boundary conditions at infinite times is invertible and that its inverse, called the Feynman inverse, satisfies the microlocal conditions of Feynman parametrices in the sense of Duistermaat and H\"ormander. This supplements the recent work of Vasy with more explicit techniques.


title: Mixed-normal limit theorems for multiple Skorohod integrals in high-dimensions, with application to realized covariance

abstract: This paper develops mixed-normal approximations for probabilities that vectors of multiple Skorohod integrals belong to random convex polytopes when the dimensions of the vectors possibly diverge to infinity. We apply the developed theory to establish the asymptotic mixed normality of the realized covariance matrix of a high-dimensional continuous semimartingale observed at a high-frequency, where the dimension can be much larger than the sample size. We also present an application of this result to testing the residual sparsity of a high-dimensional continuous-time factor model.


title: Pressure dependence of the magic twist angle in graphene superlattices

abstract: The recently demonstrated unconventional superconductivity in twisted bilayer graphene (tBLG) opens the possibility for interesting applications of two-dimensional layers that involve correlated electron states. Here we explore the possibility of modifying electronic correlations by the application of uniaxial pressure on the weakly interacting layers, which results in increased interlayer coupling and a modification of the magic angle value and associated density of states. Our findings are based on first-principles calculations that accurately describe the height-dependent interlayer coupling through the combined use of Density Functional Theory and Maximally localized Wannier functions. We obtain the relationship between twist angle and external pressure for the magic angle flat bands of tBLG. This may provide a convenient method to tune electron correlations by controlling the length scale of the superlattice.


title: Giant fluctuations in logistic growth

abstract: We analyze the fluctuation of the number of individuals when two competing species, beginning with a few initial individuals, are submitted to a logistic growth. We show that when the total number of individuals reaches the carrying capacity, the number of each species is subject to giant fluctuations (variance $\sim$ mean 2) if the two species have similar growth rate. We show that the deterministic logistic equation can be used only when the growth rates are significantly different, otherwise such growth has to be investigated by stochastic processes tools. These results generalize to a wide class of growth law.


title: Complete Optimal Convex Approximations of Qubit States under $B_2$ Distance

abstract: We consider the optimal approximation of arbitrary qubit states with respect to an available states consisting the eigenstates of two of three Pauli matrices, the $B_2$-distance of an arbitrary target state. Both the analytical formulae of the $B_2$-distance and the corresponding complete optimal decompositions are obtained. The tradeoff relations for both the sum and the squared sum of the $B_2$-distances have been analytically and numerically investigated.


title: LASSO-Driven Inference in Time and Space

abstract: We consider the estimation and inference in a system of high-dimensional regression equations allowing for temporal and cross-sectional dependency in covariates and error processes, covering rather general forms of weak temporal dependence. A sequence of regressions with many regressors using LASSO (Least Absolute Shrinkage and Selection Operator) is applied for variable selection purpose, and an overall penalty level is carefully chosen by a block multiplier bootstrap procedure to account for multiplicity of the equations and dependencies in the data. Correspondingly, oracle properties with a jointly selected tuning parameter are derived. We further provide high-quality de-biased simultaneous inference on the many target parameters of the system. We provide bootstrap consistency results of the test procedure, which are based on a general Bahadur representation for the $Z$-estimators with dependent data. Simulations demonstrate good performance of the proposed inference procedure. Finally, we apply the method to quantify spillover effects of textual sentiment indices in a financial market and to test the connectedness among sectors.


title: Generalized quantum Rabi model with both one- and two-photon terms: A concise analytical study

abstract: A generalized quantum Rabi Hamiltonian with both one- and two-photon terms has emerged in the circuit quantum electrodynamics system for a decade. The usual parity symmetry is broken naturally in the simultaneous presence of both couplings, which complicates analytical treatments, even in the rotating wave approximations. In this paper, we propose an adiabatic approximation to this generic model by using Bogoliubov operators, and obtain a very concise analytical solution for both eigenvalues and eigenstates. Although the adiabatic approximation is only exact in the vanishing limit of the qubit frequency, the results for some physical observables nevertheless agree well with the numerical ones in a wide parameter regime. In the rotating-wave approximations, we also derive an analytical eigensolution. Two dominant Rabi frequencies are found in the Rabi oscillations of this generalized model. We also apply the present analytical theory to the vacuum Rabi splitting. It is found that some new phenomena emerge just because of the presence of the additional two-photon coupling term.


title: Multiple Instance Learning for Heterogeneous Images: Training a CNN for Histopathology

abstract: Multiple instance (MI) learning with a convolutional neural network enables end-to-end training in the presence of weak image-level labels. We propose a new method for aggregating predictions from smaller regions of the image into an image-level classification by using the quantile function. The quantile function provides a more complete description of the heterogeneity within each image, improving image-level classification. We also adapt image augmentation to the MI framework by randomly selecting cropped regions on which to apply MI aggregation during each epoch of training. This provides a mechanism to study the importance of MI learning. We validate our method on five different classification tasks for breast tumor histology and provide a visualization method for interpreting local image classifications that could lead to future insights into tumor heterogeneity.


title: Tilings, packings and expected Betti numbers in simplicial complexes

abstract: Let $K$ be a finite simplicial complex. We prove that the normalized expected Betti numbers of a random subcomplex in its $d$-th barycentric subdivision $\text{Sd}^d (K)$ converge to universal limits as $d$ grows to $+ \infty$. In codimension one, we use canonical filtrations of $\text{Sd}^d (K)$ to upper estimate these limits and get a monotony theorem which makes it possible to improve these estimates given any packing of disjoint simplices in $\text{Sd}^d (K)$. We then introduce a notion of tiling of simplicial complexes having the property that skeletons and barycentric subdivisions of tileable simplicial complexes are tileable. This enables us to tackle the problem: How many disjoint simplices can be packed in $\text{Sd}^d (K)$, $d \gg 0$?


title: Your 2 is My 1, Your 3 is My 9: Handling Arbitrary Miscalibrations in Ratings

abstract: Cardinal scores (numeric ratings) collected from people are well known to suffer from miscalibrations. A popular approach to address this issue is to assume simplistic models of miscalibration (such as linear biases) to de-bias the scores. This approach, however, often fares poorly because people's miscalibrations are typically far more complex and not well understood. In the absence of simplifying assumptions on the miscalibration, it is widely believed by the crowdsourcing community that the only useful information in the cardinal scores is the induced ranking. In this paper, inspired by the framework of Stein's shrinkage, empirical Bayes, and the classic two-envelope problem, we contest this widespread belief. Specifically, we consider cardinal scores with arbitrary (or even adversarially chosen) miscalibrations which are only required to be consistent with the induced ranking. We design estimators which despite making no assumptions on the miscalibration, strictly and uniformly outperform all possible estimators that rely on only the ranking. Our estimators are flexible in that they can be used as a plug-in for a variety of applications, and we provide a proof-of-concept for A/B testing and ranking. Our results thus provide novel insights in the eternal debate between cardinal and ordinal data.


title: Group Equivariant Capsule Networks

abstract: We present group equivariant capsule networks, a framework to introduce guaranteed equivariance and invariance properties to the capsule network idea. Our work can be divided into two contributions. First, we present a generic routing by agreement algorithm defined on elements of a group and prove that equivariance of output pose vectors, as well as invariance of output activations, hold under certain conditions. Second, we connect the resulting equivariant capsule networks with work from the field of group convolutional networks. Through this connection, we provide intuitions of how both methods relate and are able to combine the strengths of both approaches in one deep neural network architecture. The resulting framework allows sparse evaluation of the group convolution operator, provides control over specific equivariance and invariance properties, and can use routing by agreement instead of pooling operations. In addition, it is able to provide interpretable and equivariant representation vectors as output capsules, which disentangle evidence of object existence from its pose.


title: Seshadri constants of the anticanonical divisors of Fano manifolds with large index

abstract: Seshadri constants, introduced by Demailly, measure the local positivity of a nef divisor at a point. In this paper, we compute the Seshadri constants of the anticanonical divisors of Fano manifolds with coindex at most $3$ at a very general point. As a consequence, if $X$ is a nonsingular Fano threefold which is very general in its deformation family, then $\varepsilon(X,-K_X;x)\leq 1$ for all points $x\in X$ if and only if $\vert-K_X\vert$ is not base point free.


title: Evolution from Ferromagnetism to Antiferromagnetism in Yb(Rh1-xCox)2Si2

abstract: Yb(Rh1-xCox)2Si2 is a model system to address two challenging problems in the field of strongly correlated electron systems: The first is the intriguing competition between ferromagnetic (FM) and antiferromagnetic (AFM) order when approaching a magnetic quantum critical point (QCP). The second is the occurrence of magnetic order along a very hard crystalline electric field (CEF) direction, i.e. along the one with the smallest available magnetic moment. Here, we present a detailed study of the evolution of the magnetic order in this system from a FM state with moments along the very hard c direction at x = 0.27 towards the yet unknown magnetic state at x = 0. We first observe a transition towards an AFM canted state with decreasing x and then to a pure AFM state. This confirms that the QCP in YbRh2Si2 is AFM, but the phase diagram is very similar to those observed in some inherently FM systems like NbFe2 and CeRuPO, which suggests that the basic underlying instability might be FM. Despite the huge CEF anisotropy the ordered moment retains a component along the c-axis also in the AFM state. The huge CEF anisotropy in Yb(Rh1-xCox)2Si2 excludes that this hard-axis ordering originates from a competing exchange anisotropy as often proposed for other heavy-fermion systems. Instead, it points to an order-by-disorder based mechanism.


title: Particles-on-Demand for Kinetic Theory

abstract: A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes any restrictions on Mach number and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. In the new kinetic theory, discrete particles are determined by a rigorous limit process which avoids ad hoc assumptions about their velocities. Classical benchmarks for incompressible and compressible flows demonstrate that the proposed discrete-particles kinetic theory opens up an unprecedented wide domain of applications for computational fluid dynamics.


title: Exploring new physics contributions to $CP$ violation in $\tau^- \to K^-\pi^0\nu_{\tau} $

abstract: A general analysis of possible violation of CP in processes like $\tau \to K\pi \nu$, for unpolarized $\tau$ is presented. In this paper, we derive the new contributions to the effective Hamiltonian governs $\vert\Delta S \vert=1$ semileptonic tau decays in the framework of two Higgs doublet model with generic Yukawa structure and Leptoquarks models. Within these models, we list all operators, in the effective Hamiltonian and provide analytical expression for their corresponding Wilson coefficients. Moreover, we analyze the role of the different contributions, originating from the scalar, vecor and tensor hadronic currents, in generating direct CP asymmetry in the decay rate of $\tau^-\to K^-\pi^0\nu_\tau$. We show that non vanishing direct CP asymmetry in the decay rate of $\tau^-\to K^-\pi^0\nu_\tau$ can be generated due to the presence of both, the weak phase in the Wilson coefficient corresponding to the tensor operator and the strong phase difference resulting from the interference between the form factors expressing the matrix elements of the vector and tensor hadronic currents. After taking into account all relevant constraints, we find that the generated direct CP asymmetry is of order $10^{-8}$ which is several orders of magnitude larger than the standard model prediction. We show also that, in two Higgs doublet model with generic Yukawa structure , direct local or non integrated CP violation can be as large as $0.3$ % not far from experimental possibilities. This kind of asymmetry can be generated due to the interference between vector and scalar contributions with different weak phases which is not the case in the SM.


title: Estimating Achilles tendon healing progress with convolutional neural networks

abstract: Quantitative assessment of a treatment progress in the Achilles tendon healing process - one of the most common musculoskeletal disorder in modern medical practice - is typically a long and complex process: multiple MRI protocols need to be acquired and analysed by radiology experts. In this paper, we propose to significantly reduce the complexity of this assessment using a novel method based on a pre-trained convolutional neural network. We first train our neural network on over 500,000 2D axial cross-sections from over 3000 3D MRI studies to classify MRI images as belonging to a healthy or injured class, depending on the patient's condition. We then take the outputs of modified pre-trained network and apply linear regression on the PCA-reduced space of the features to assess treatment progress. Our method allows to reduce up to 5-fold the amount of data needed to be registered during the MRI scan without any information loss. Furthermore, we are able to predict the healing process phase with equal accuracy to human experts in 3 out of 6 main criteria. Finally, contrary to the current approaches to regeneration assessment that rely on radiologist subjective opinion, our method allows to objectively compare different treatments methods which can lead to improved diagnostics and patient's recovery.


title: A survey on fractional variational calculus

abstract: Main results and techniques of the fractional calculus of variations are surveyed. We consider variational problems containing Caputo derivatives and study them using both indirect and direct methods. In particular, we provide necessary optimality conditions of Euler-Lagrange type for the fundamental, higher-order, and isoperimetric problems, and compute approximated solutions based on truncated Gr\"{u}nwald--Letnikov approximations of Caputo derivatives.


title: Finite element methods for fourth order axisymmetric geometric evolution equations

abstract: Fourth order curvature driven interface evolution equations frequently appear in the natural sciences. Often axisymmetric geometries are of interest, and in this situation numerical computations are much more efficient. We will introduce and analyze several new finite element schemes for fourth order geometric evolution equations in an axisymmetric setting, and for selected schemes we will show existence, uniqueness and stability results. The presented schemes have very good mesh and stability properties, as will be demonstrated by several numerical examples.


title: A combinatorial approach to scattering diagrams

abstract: Scattering diagrams arose in the context of mirror symmetry, but a special class of scattering diagrams (the cluster scattering diagrams) were recently developed to prove key structural results on cluster algebras. We use the connection to cluster algebras to calculate the function attached to the limiting wall of a rank-2 cluster scattering diagram of affine type. In the skew-symmetric rank-2 affine case, this recovers a formula due to Reineke. In the same case, we show that the generating function for signed Narayana numbers appears in a role analogous to a cluster variable. In acyclic finite type, we construct cluster scattering diagrams of acyclic finite type from Cambrian fans and sortable elements, with a simple direct proof.


title: Optimal moment inequalities for order statistics from nonnegative random variables

abstract: We obtain the best possible upper bounds for the moments of a single order statistic from independent, non-negative random variables, in terms of the population mean. The main result covers the independent identically distributed case. Furthermore, the case of the sample minimum for merely independent (not necessarily identically distributed) random variables is treated in detail. Key-words and phrases: order statistics; optimal moment bounds; nonnegative random variables; sample minimum; reliability systems.


title: Introducing user-prescribed constraints in Markov chains for nonlinear dimensionality reduction

abstract: Stochastic kernel based dimensionality reduction approaches have become popular in the last decade. The central component of many of these methods is a symmetric kernel that quantifies the vicinity between pairs of data points and a kernel-induced Markov chain on the data. Typically, the Markov chain is fully specified by the kernel through row normalization. However, in many cases, it is desirable to impose user-specified stationary-state and dynamical constraints on the Markov chain. Unfortunately, no systematic framework exists to impose such user-defined constraints. Here, we introduce a path entropy maximization based approach to derive the transition probabilities of Markov chains using a kernel and additional user-specified constraints. We illustrate the usefulness of these Markov chains with examples.


title: On GW170817 and the Galactic Binary Neutron Star Population

abstract: GW170817/GRB170817A, a short gamma-ray burst arising from a low-mass compact object merger was the first multi-messenger discovery of a compact binary system outside the local galactic neighborhood. From gravitational-wave measurements, we know GW170817 has a wide range of plausible component masses, depending also on less well-constrained properties such as the spin and tidal deformability of the component stars. The kilonova light curve --- and hence the total ejecta mass from a given source --- depends on the relative contribution of dynamical ejecta and other sources such as disk winds. Electromagnetic observations and model fitting of the ejecta properties from the subsequent kilonova detection provided values of the ejecta mass from the merger. These values, when combined with the gravitational-wave measurement disfavors an equal-mass configuration, with the level of disagreement dependent on the assumed amount of ejecta mass of dynamical origin. Within the confines of our own galaxy, several binary neutron star systems along with measurements of their component masses have been made. If those distributions are indicative of a universal distribution, the joint measurement of the component masses of GW170817 represents an outlier. This tension is not easily resolvable from physical arguments, as the proposed pathways which form binary neutron stars do not often produce very asymmetrical pairs. Even accounting for the uncertainty associated with the total mass of the dynamical ejecta, this tension suggests that the distribution of binary neutron star masses in the galaxy is not indicative of those in other galaxies.


title: Bridging thermodynamics and metrology in non-equilibrium Quantum Thermometry

abstract: Single-qubit thermometry presents the simplest tool to measure the temperature of thermal baths with reduced invasivity. At thermal equilibrium, the temperature uncertainty is linked to the heat capacity of the qubit, however the best precision is achieved outside equilibrium condition. Here, we discuss a way to generalize this relation in a non-equilibrium regime, taking into account purely quantum effects such as coherence. We support our findings with an experimental photonic simulation.


title: Graph-Based Decoding for Event Sequencing and Coreference Resolution

abstract: Events in text documents are interrelated in complex ways. In this paper, we study two types of relation: Event Coreference and Event Sequencing. We show that the popular tree-like decoding structure for automated Event Coreference is not suitable for Event Sequencing. To this end, we propose a graph-based decoding algorithm that is applicable to both tasks. The new decoding algorithm supports flexible feature sets for both tasks. Empirically, our event coreference system has achieved state-of-the-art performance on the TAC-KBP 2015 event coreference task and our event sequencing system beats a strong temporal-based, oracle-informed baseline. We discuss the challenges of studying these event relations.


title: The Energy-Dependent Black-Disk Fraction in Proton-Proton Scattering

abstract: Recent work by Block et al on the energy-invariant edge in proton-proton scattering is interpreted by quantitatively introducing a energy-dependent black-disk fraction of the proton-proton interaction.


title: Order-book modelling and market making strategies

abstract: Market making is one of the most important aspects of algorithmic trading, and it has been studied quite extensively from a theoretical point of view. The practical implementation of so-called "optimal strategies" however suffers from the failure of most order book models to faithfully reproduce the behaviour of real market participants. This paper is twofold. First, some important statistical properties of order driven markets are identified, advocating against the use of purely Markovian order book models. Then, market making strategies are designed and their performances are compared, based on simulation as well as backtesting. We find that incorporating some simple non-Markovian features in the limit order book greatly improves the performances of market making strategies in a realistic context.


title: Combined feedback and sympathetic cooling of a mechanical oscillator coupled to ultracold atoms

abstract: A promising route to novel quantum technologies are hybrid quantum systems, which combine the advantages of several individual quantum systems. We have realized a hybrid atomic-mechanical experiment consisting of a SiN membrane oscillator cryogenically precooled to 500 mK and optically coupled to a cloud of laser cooled Rb atoms. Here, we demonstrate active feedback cooling of the oscillator to a minimum mode occupation of n = 16 corresponding to a mode temperature of T = 200 {\mu}K. Furthermore, we characterize in detail the coupling of the membrane to the atoms by means of sympathetic cooling. By simultaneously applying both cooling methods we demonstrate the possibility of preparing the oscillator near the motional ground state while it is coupled to the atoms. Realistic modifications of our setup will enable the creation of a ground state hybrid quantum system, which opens the door for coherent quantum state transfer, teleportation and entanglement as well as quantum enhanced sensing applications.


title: A new non-perturbative approach in quantum mechanics for time-independent Schr\"{o}dinger equations

abstract: A new non-perturbative approach is proposed to solve time-independent Schr\"{o}dinger equations in quantum mechanics and chromodynamics (QCD). It is based on the homotopy analysis method (HAM), which was developed by the author for highly nonlinear equations since 1992 and has been widely applied in many fields. Unlike perturbative methods, this HAM-based approach has nothing to do with small/large physical parameters. Besides, convergent series solution can be obtained even if the disturbance is far from the known status. A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be more than hundreds larger than the possible superior limit of the perturbative approach. This HAM-based approach could provide us rigorous theoretical results in quantum mechanics and chromodynamics (QCD), which can be directly compared with experimental data. Obviously, this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.


title: Improving Cytoarchitectonic Segmentation of Human Brain Areas with Self-supervised Siamese Networks

abstract: Cytoarchitectonic parcellations of the human brain serve as anatomical references in multimodal atlas frameworks. They are based on analysis of cell-body stained histological sections and the identification of borders between brain areas. The de-facto standard involves a semi-automatic, reproducible border detection, but does not scale with high-throughput imaging in large series of sections at microscopical resolution. Automatic parcellation, however, is extremely challenging due to high variation in the data, and the need for a large field of view at microscopic resolution. The performance of a recently proposed Convolutional Neural Network model that addresses this problem especially suffers from the naturally limited amount of expert annotations for training. To circumvent this limitation, we propose to pre-train neural networks on a self-supervised auxiliary task, predicting the 3D distance between two patches sampled from the same brain. Compared to a random initialization, fine-tuning from these networks results in significantly better segmentations. We show that the self-supervised model has implicitly learned to distinguish several cortical brain areas -- a strong indicator that the proposed auxiliary task is appropriate for cytoarchitectonic mapping.


title: Stability and complexity of mixed discriminants

abstract: We show that the mixed discriminant of $n$ positive semidefinite $n \times n$ real symmetric matrices can be approximated within a relative error $\epsilon >0$ in quasi-polynomial $n^{O(\ln n -\ln \epsilon)}$ time, provided the distance of each matrix to the identity matrix in the operator norm does not exceed some absolute constant $\gamma_0 >0$. We deduce a similar result for the mixed discriminant of doubly stochastic $n$-tuples of matrices from the Marcus - Spielman - Srivastava bound on the roots of the mixed characteristic polynomial. Finally, we construct a quasi-polynomial algorithm for approximating the sum of $m$-th powers of principal minors of a matrix, provided the operator norm of the matrix is strictly less than 1. As is shown by Gurvits, for $m=2$ the problem is $\#P$-hard and covers the problem of computing the mixed discriminant of positive semidefinite matrices of rank 2.


title: DRE-Bot: A Hierarchical First Person Shooter Bot Using Multiple Sarsa({\lambda}) Reinforcement Learners

abstract: This paper describes an architecture for controlling non-player characters (NPC) in the First Person Shooter (FPS) game Unreal Tournament 2004. Specifically, the DRE-Bot architecture is made up of three reinforcement learners, Danger, Replenish and Explore, which use the tabular Sarsa({\lambda}) algorithm. This algorithm enables the NPC to learn through trial and error building up experience over time in an approach inspired by human learning. Experimentation is carried to measure the performance of DRE-Bot when competing against fixed strategy bots that ship with the game. The discount parameter, {\gamma}, and the trace parameter, {\lambda}, are also varied to see if their values have an effect on the performance.


title: Simplicial complexes of small codimension

abstract: We show that a Buchsbaum simplicial complex of small codimension must have large depth. More generally, we achieve a similar result for ${\rm CM}_t$ simplicial complexes, a notion generalizing Buchsbaum-ness, and we prove more precise results in the codimension 2 case. Along the paper, we show that the ${\rm CM}_t$ property is a topological invariant of a simplicial complex.


title: Holographic Automata for Ambient Immersive A. I. via Reservoir Computing

abstract: We prove the existence of a semilinear representation of Cellular Automata (CA) with the introduction of multiple convolution kernels. Examples of the technique are presented for rules akin to the "edge-of-chaos" including the Turing universal rule 110 for further utilization in the area of reservoir computing. We also examine the significance of their dual representation on a frequency or wavelength domain as a superposition of plane waves for distributed computing applications including a new proposal for a "Hologrid" that could be realized with present Wi-Fi,Li-Fi technologies.


title: Statistical mechanics of an elastically pinned membrane: Static profile and correlations

abstract: The relation between thermal fluctuations and the mechanical response of a free membrane has been explored in great detail, both theoretically and experimentally. However, understanding this relationship for membranes, locally pinned by proteins, is significantly more challenging. Given that the coupling of the membrane to the cell cytoskeleton, the extracellular matrix and to other internal structures is crucial for the regulation of a number of cellular processes, understanding the role of the pinning is of great interest. In this manuscript we consider a single protein (elastic spring of a finite rest length) pinning a membrane modelled in the Monge gauge. First, we determine the Green$'$s function for the system and complement this approach by the calculation of the mode coupling coefficients for the plane wave expansion, and the orthonormal fluctuation modes, in turn building a set of tools for numerical and analytic studies of a pinned membrane. Furthermore, we explore static correlations of the free and the pinned membrane, as well as the membrane shape, showing that all three are mutually interdependent and have an identical long-range behaviour characterised by the correlation length. Interestingly, the latter displays a non-monotonic behaviour as a function of membrane tension. Importantly, exploiting these relations allows for the experimental determination of the elastic parameters of the pinning. Last but not least, we calculate the interaction potential between two pinning sites and show that, even in the absence of the membrane deformation, the pinnings will be subject to an attractive force due to changes in membrane fluctuations.


title: Failure Analysis of the Interval-Passing Algorithm for Compressed Sensing

abstract: In this work, we perform a complete failure analysis of the interval-passing algorithm (IPA) for compressed sensing, an efficient iterative algorithm for reconstructing a $k$-sparse nonnegative $n$-dimensional real signal $\boldsymbol{x}$ from a small number of linear measurements $\boldsymbol{y}$. In particular, we show that the IPA fails to recover $\boldsymbol{x}$ from $\boldsymbol{y}$ if and only if it fails to recover a corresponding binary vector of the same support, and also that only positions of nonzero values in the measurement matrix are of importance to the success of recovery. Based on this observation, we introduce termatiko sets and show that the IPA fails to fully recover $\boldsymbol x$ if and only if the support of $\boldsymbol x$ contains a nonempty termatiko set, thus giving a complete (graph-theoretic) description of the failing sets of the IPA. Two heuristics to locate small-size termatiko sets are presented. For binary column-regular measurement matrices with no $4$-cycles, we provide a lower bound on the termatiko distance, defined as the smallest size of a nonempty termatiko set. For measurement matrices constructed from the parity-check matrices of array LDPC codes, upper bounds on the termatiko distance are provided for column-weight at most $7$, while for column-weight $3$, the exact termatiko distance and its corresponding multiplicity are provided. Next, we show that adding redundant rows to the measurement matrix does not create new termatiko sets, but rather potentially removes termatiko sets and thus improves performance. An algorithm is provided to efficiently search for such redundant rows. Finally, we present numerical results for different specific measurement matrices and also for protograph-based ensembles of measurement matrices, as well as simulation results of IPA performance, showing the influence of small-size termatiko sets.


title: A Chandra/LETGS Survey of Main Sequence Stars

abstract: We analyze the X-ray spectra of 19 main sequence stars observed by Chandra using its LETGS configuration. Emission measure (EM) distributions are computed based on emission line measurements, an analysis that also yields evaluations of coronal abundances. The use of newer atomic physics data results in significant changes compared to past published analyses. The stellar EM distributions correlate with surface X-ray flux (F_X) in a predictable way, regardless of spectral type. Thus, we provide EM distributions as a function of F_X, which can be used to estimate the EM distribution of any main sequence star with a measured broadband X-ray luminosity. Comparisons are made with solar EM distributions, both full-disk distributions and spatially resolved ones from active regions (ARs), flares, and the quiet Sun. For moderately active stars, the slopes and magnitudes of the EM distributions are in excellent agreement with those of solar ARs for log T<6.6, suggesting that such stars have surfaces completely filled with solar-like ARs. A stellar surface covered with solar X-class flares yields a reasonable approximation for the EM distributions of the most active stars. Unlike the EM distributions, coronal abundances are very spectral-type dependent, and we provide relations with surface temperature for both relative and absolute abundances. Finally, the coronal abundances of the exoplanet host star Tau Boo A (F7 V) are anomalous, and we propose that this is due to the presence of the exoplanet.


title: Comparing Fairness Criteria Based on Social Outcome

abstract: Fairness in algorithmic decision-making processes is attracting increasing concern. When an algorithm is applied to human-related decision-making an estimator solely optimizing its predictive power can learn biases on the existing data, which motivates us the notion of fairness in machine learning. while several different notions are studied in the literature, little studies are done on how these notions affect the individuals. We demonstrate such a comparison between several policies induced by well-known fairness criteria, including the color-blind (CB), the demographic parity (DP), and the equalized odds (EO). We show that the EO is the only criterion among them that removes group-level disparity. Empirical studies on the social welfare and disparity of these policies are conducted.


title: The $X(4140)$ and $X(4160)$ resonances in the $e^+e^-\to \gamma J/\psi \phi $ reaction

abstract: We investigate the $J/\psi \phi$ invariant mass distribution of the $e^+e^-\to \gamma J/\psi\phi$ reaction at a center of mass energy of $\sqrt{s}=4.6$~GeV measured by the BESIII collaboration, who concludes that no significant signals are observed for $e^+e^- \to \gamma X(4140)$ because of the low statistics. We show, however, that the three bump structures in the $J/\psi \phi$ invariant mass distribution, though given only by three data points, are compatible with the existence of the $X(4140)$ state, appearing as a peak, and a strong cusp structure at the $D^*_s\bar{D}^*_s$ threshold, resulting from the molecular nature of the $X(4160)$ state, which also provides a substantial contribution to the reaction. This is consistent with our previous analysis on the $B^+\to J/\psi\phi K^+$ decay measured by the LHCb collaboration. We strongly call for the measurement of this process with more statistics to further pin down the nature of the $X(4140)$ and $X(4160)$ resonances.


title: Earthquake Detection in 1-D Time Series Data with Feature Selection and Dictionary Learning

abstract: Earthquakes can be detected by matching spatial patterns or phase properties from 1-D seismic waves. Current earthquake detection methods, such as waveform correlation and template matching, have difficulty detecting anomalous earthquakes that are not similar to other earthquakes. In recent years, machine-learning techniques for earthquake detection have been emerging as a new active research direction. In this paper, we develop a novel earthquake detection method based on dictionary learning. Our detection method first generates rich features via signal processing and statistical methods and further employs feature selection techniques to choose features that carry the most significant information. Based on these selected features, we build a dictionary for classifying earthquake events from non-earthquake events. To evaluate the performance of our dictionary-based detection methods, we test our method on a labquake dataset from Penn State University, which contains 3,357,566 time series data points with a 400 MHz sampling rate. 1,000 earthquake events are manually labeled in total, and the length of these earthquake events varies from 74 to 7151 data points. Through comparison to other detection methods, we show that our feature selection and dictionary learning incorporated earthquake detection method achieves an 80.1% prediction accuracy and outperforms the baseline methods in earthquake detection, including Template Matching (TM) and Support Vector Machine (SVM).


title: Interplay between electronic correlation and metal-ligand delocalization in the spectroscopy of transition metal compounds: case study on a series of planar Cu$^{2+}$ complexes

abstract: We present a comprehensive theoretical study of the physical phenomena that determine the relative energies of the three of the lowest electronic states of each of the square-planar copper complexes $\cucl$, $\cunh$ and $\cuwater$, and present a detailed analysis of the extent to which truncated configuration interaction (CI) and coupled cluster (CC) theories succeed in predicing the excitation energies. We find that ligand-metal charge transfer (CT) single excitations play a crucial role in the correct determination of the properties of these systems, even though the CT processes first occur at fourth order in perturbation theory, and propose a suitable choice of minimal active space for describing these systems with multi-reference theories. CCSD energy differences agree very well with near full CI values even though the $T_1$ diagnostics are large, which casts doubt on the usefulness of singles-amplitude based multi-reference diagnostics. CISD severely underestimates the excitation energies and the failure is a direct consequence of the size-inconsisency errors in CISD. Finally, we present reference values for the energy differences computed using explicitly correlated CCSD(T) and BCCD(T) theory.


title: Is Fickian Yet Non-Gaussian Diffusion Ubiquitous?

abstract: Recent studies unveiled the Fickian yet non-Gaussian (FNG) dynamics of many soft matter systems and suggested this phenomenon as a general characteristic of the diffusion in complex fluids. In particular, it was shown that the distribution of particle displacements in Fickian diffusion is not necessarily Gaussian, and thus the Einstein and Smoluchowski theory describing the Brownian motion of individual objects in a fluid would not be applicable. In this Letter, we investigate whether the FNG dynamics so far reported in gels, granular materials, biological and active matter systems, is also a distinctive feature of colloidal liquid crystals. To this end, we perform Brownian Dynamics simulations of oblate and prolate colloidal particles in the nematic phase. We detect a normal and Gaussian dynamics at short and long time scales, whereas, at intermediate time scales, a non-Fickian and non-Gaussian dynamics is found. Additionally, we revisit the nature of the decay of the self-van Hove correlation function, $G_s(r,t)$, which is here approximated with an ellipsoidal, rather than spherical, Gaussian distribution. The new expression that we propose is able to correctly assess the Gaussian dynamics in inherently anisotropic systems, like liquid crystals, where the standard Gaussian approximation of $G_s(r,t)$ would fail.


title: Learning to Shoot in First Person Shooter Games by Stabilizing Actions and Clustering Rewards for Reinforcement Learning

abstract: While reinforcement learning (RL) has been applied to turn-based board games for many years, more complex games involving decision-making in real-time are beginning to receive more attention. A challenge in such environments is that the time that elapses between deciding to take an action and receiving a reward based on its outcome can be longer than the interval between successive decisions. We explore this in the context of a non-player character (NPC) in a modern first-person shooter game. Such games take place in 3D environments where players, both human and computer-controlled, compete by engaging in combat and completing task objectives. We investigate the use of RL to enable NPCs to gather experience from game-play and improve their shooting skill over time from a reward signal based on the damage caused to opponents. We propose a new method for RL updates and reward calculations, in which the updates are carried out periodically, after each shooting encounter has ended, and a new weighted-reward mechanism is used which increases the reward applied to actions that lead to damaging the opponent in successive hits in what we term "hit clusters".


title: Exact solution for driven oscillations in plasmonic field-effect transistors

abstract: High-mobility field effect transistors can serve as resonant detectors of terahertz radiation due to excitation of plasmons in the channel. The modeling of these devices previously relied either on approximate techniques, or complex full-wave simulations. In this paper, we obtain an exact solution for driven electrical oscillations in plasmonic field-effect transistor with realistic contact geometry. The obtained solution highlights the importance of evanescent plasma waves excited near the contacts, which qualitatively modify the detector responsivity spectra. We derive the boundary condition on the ac floating electrodes of plasmonic FET which interpolates between open-circuit (Dyakonov-Shur) and short-circuit (clamped voltage) boundary conditions. In both limits, the FET photovoltage possesses resonant fringes, however, the absolute value of voltage is greater in the open-circuit regime.


title: Long monochromatic paths and cycles in 2-colored bipartite graphs

abstract: Gy\'arf\'as and Lehel and independently Faudree and Schelp proved that in any 2-coloring of the edges of $K_{n,n}$ there exists a monochromatic path on at least $2\lceil n/2\rceil$ vertices, and this is tight. We prove a stability version of this result which holds even if the host graph is not complete; that is, if $G$ is a balanced bipartite graph on $2n$ vertices with minimum degree at least $(3/4+o(1))n$, then in every 2-coloring of the edges of $G$, either there exists a monochromatic cycle on at least $(1+o(1))n$ vertices, or the coloring of $G$ is close to an extremal coloring -- in which case $G$ has a monochromatic path on at least $2\lceil n/2\rceil$ vertices and a monochromatic cycle on at least $2\lfloor n/2\rfloor$ vertices. Furthermore, we determine an asymptotically tight bound on the length of a longest monochromatic cycle in a 2-colored balanced bipartite graph on $2n$ vertices with minimum degree $\delta n$ for all $0\leq \delta\leq 1$.


title: First results from a hidden photon dark matter search in the meV sector using a plane-parabolic mirror system

abstract: We report on the first results from a new dish antenna search for hidden photon dark matter (HPDM) in the meV mass region. A double mirror system composed of a plane and a parabolic mirror is designed to convert HPDMs into photons focused on a receiver. In this phase 1 experiment we obtain an upper limit on the photon-HP kinetic mixing $\chi \lesssim 10^{-8}$ for the mass range of $0.67-0.92\,{\rm meV}$ using conventional mm-wave technology with a room-temperature receiver and a small-sized mirror system.


title: Adaptive Path Interpolation for Sparse Systems: Application to a Simple Censored Block Model

abstract: Recently a new adaptive path interpolation method has been developed as a simple and versatile scheme to calculate exactly the asymptotic mutual information of Bayesian inference problems defined on dense factor graphs. These include random linear and generalized estimation, sparse superposition codes, or low-rank matrix and tensor estimation. For all these systems, the adaptive interpolation method directly proves that the replica symmetric prediction is exact, in a simple and unified manner. When the underlying factor graph of the inference problem is sparse the replica prediction is considerably more complicated, and rigorous results are often lacking or obtained by rather complicated methods. In this work we show how to extend the adaptive path interpolation method to sparse systems. We concentrate on a Censored Block Model, where hidden variables are measured through a binary erasure channel, for which we fully prove the replica prediction.

categories: cs.IT cond-mat.dis-nn math-ph math.IT math.MP


title: Photon-assisted quantum state transfer and entanglement generation in spin chains

abstract: We propose a protocol for state transfer and entanglement generation between two distant spin qubits (sender and receiver) that have different energies. The two qubits are permanently coupled to a far off-resonant spin-chain, and the qubit of the sender is driven by an external field, which provides the energy required to bridge the energy gap between the sender and the receiver. State transfer and entanglement generation are achieved via virtual single-photon and multi-photon transitions to the eigenmodes of the channel.


title: Linearly Convergent Frank-Wolfe with Backtracking Line-Search

abstract: Structured constraints in Machine Learning have recently brought the Frank-Wolfe (FW) family of algorithms back in the spotlight. While the classical FW algorithm has poor local convergence properties, the Away-steps and Pairwise FW variants have emerged as improved variants with faster convergence. However, these improved variants suffer from two practical limitations: they require at each iteration to solve a 1-dimensional minimization problem to set the step-size and also require the Frank-Wolfe linear subproblems to be solved exactly. In this paper, we propose variants of Away-steps and Pairwise FW that lift both restrictions simultaneously. The proposed methods set the step-size based on a sufficient decrease condition, and do not require prior knowledge of the objective. Furthermore, they inherit all the favorable convergence properties of the exact line-search version, including linear convergence for strongly convex functions over polytopes. Benchmarks on different machine learning problems illustrate large performance gains of the proposed variants.


title: Unraveling conformal gravity amplitudes

abstract: Conformal supergravity amplitudes are obtained from the double-copy construction using gauge-theory amplitudes, and compared to direct calculations starting from conformal supergravity Lagrangians. We consider several different theories: minimal ${\cal N}=4$ conformal supergravity, non-minimal ${\cal N}=4$ Berkovits-Witten conformal supergravity, mass-deformed versions of these theories, as well as supersymmetry truncations thereof. Coupling the theories to a Yang-Mills sector is also considered. For all cases we give the gravity Lagrangians that the double copy implicitly generates. The two main results are: we determine a Lagrangian for the non-minimal Berkovits-Witten theory, and we uncover the double-copy prescription for the minimal ${\cal N}=4$ conformal supergravity.


title: Observability of planet-disc interactions in CO kinematics

abstract: Empirical evidence of planets in gas-rich circumstellar discs is required to constrain giant planet formation theories. Here we study the kinematic patterns which arise from planet-disc interactions and their observability in CO rotational emission lines. We perform three-dimensional hydrodynamical simulations of single giant planets, and predict the emergent intensity field with radiative transfer. Pressure gradients at planet-carved gaps, spiral wakes and vortices bear strong kinematic counterparts. The iso-velocity contours in the CO(2-1) line centroids $v_\circ$ reveal large-scale perturbations, corresponding to abrupt transitions from below sub-Keplerian to super-Keplerian rotation along with radial and vertical flows. The increase in line optical depth at the edge of the gap also modulates $v_\circ$, but this is a mild effect compared to the dynamical imprint of the planet-disc interaction. The large-scale deviations from the Keplerian rotation thus allow the planets to be indirectly detected via the first moment maps of molecular gas tracers, at ALMA angular resolutions. The strength of these deviations depends on the mass of the perturber. This initial study paves the way to eventually determine the mass of the planet by comparison with more detailed models.


title: Parameter-Independent Strategies for pMDPs via POMDPs

abstract: Markov Decision Processes (MDPs) are a popular class of models suitable for solving control decision problems in probabilistic reactive systems. We consider parametric MDPs (pMDPs) that include parameters in some of the transition probabilities to account for stochastic uncertainties of the environment such as noise or input disturbances. We study pMDPs with reachability objectives where the parameter values are unknown and impossible to measure directly during execution, but there is a probability distribution known over the parameter values. We study for the first time computing parameter-independent strategies that are expectation optimal, i.e., optimize the expected reachability probability under the probability distribution over the parameters. We present an encoding of our problem to partially observable MDPs (POMDPs), i.e., a reduction of our problem to computing optimal strategies in POMDPs. We evaluate our method experimentally on several benchmarks: a motivating (repeated) learner model; a series of benchmarks of varying configurations of a robot moving on a grid; and a consensus protocol.


title: Stratification Trees for Adaptive Randomization in Randomized Controlled Trials

abstract: This paper proposes an adaptive randomization procedure for two-stage randomized controlled trials. The method uses data from a first-wave experiment in order to determine how to stratify in a second wave of the experiment, where the objective is to minimize the variance of an estimator for the average treatment effect (ATE). We consider selection from a class of stratified randomization procedures which we call stratification trees: these are procedures whose strata can be represented as decision trees, with differing treatment assignment probabilities across strata. By using the first wave to estimate a stratification tree, we simultaneously select which covariates to use for stratification, how to stratify over these covariates, as well as the assignment probabilities within these strata. Our main result shows that using this randomization procedure with an appropriate estimator results in an asymptotic variance which is minimal in the class of stratification trees. Moreover, the results we present are able to accommodate a large class of assignment mechanisms within strata, including stratified block randomization. In a simulation study, we find that our method, paired with an appropriate cross-validation procedure ,can improve on ad-hoc choices of stratification. We conclude by applying our method to the study in Karlan and Wood (2017), where we estimate stratification trees using the first wave of their experiment.


title: Positive powers of the Laplacian in the half-space under Dirichlet boundary conditions

abstract: We present explicit formulas for solutions to nonhomogeneous boundary value problems involving any positive power of the Laplacian in the half-space. For non-integer powers the operator becomes nonlocal and this requires a suitable extension of Dirichlet-type boundary conditions. A key ingredient in our proofs is a point inversion transformation which preserves harmonicity and allows us to use known results for the ball. We include uniqueness statements, regularity estimates, and describe the growth or decay of solutions at infinity and at the boundary.


title: What Is It Like Down There? Generating Dense Ground-Level Views and Image Features From Overhead Imagery Using Conditional Generative Adversarial Networks

abstract: This paper investigates conditional generative adversarial networks (cGANs) to overcome a fundamental limitation of using geotagged media for geographic discovery, namely its sparse and uneven spatial distribution. We train a cGAN to generate ground-level views of a location given overhead imagery. We show the "fake" ground-level images are natural looking and are structurally similar to the real images. More significantly, we show the generated images are representative of the locations and that the representations learned by the cGANs are informative. In particular, we show that dense feature maps generated using our framework are more effective for land-cover classification than approaches which spatially interpolate features extracted from sparse ground-level images. To our knowledge, ours is the first work to use cGANs to generate ground-level views given overhead imagery and to explore the benefits of the learned representations.


title: Detecting Speech Act Types in Developer Question/Answer Conversations During Bug Repair

abstract: This paper targets the problem of speech act detection in conversations about bug repair. We conduct a "Wizard of Oz" experiment with 30 professional programmers, in which the programmers fix bugs for two hours, and use a simulated virtual assistant for help. Then, we use an open coding manual annotation procedure to identify the speech act types in the conversations. Finally, we train and evaluate a supervised learning algorithm to automatically detect the speech act types in the conversations. In 30 two-hour conversations, we made 2459 annotations and uncovered 26 speech act types. Our automated detection achieved 69% precision and 50% recall. The key application of this work is to advance the state of the art for virtual assistants in software engineering. Virtual assistant technology is growing rapidly, though applications in software engineering are behind those in other areas, largely due to a lack of relevant data and experiments. This paper targets this problem in the area of developer Q/A conversations about bug repair.


title: Synthesizing simulation and field data of solar irradiance

abstract: Predicting the intensity and amount of sunlight as a function of location and time is an essential component in identifying promising locations for economical solar farming. Although weather models and irradiance data are relatively abundant, these have yet, to our knowledge, been hybridized on a continental scale. Rather, much of the emphasis in the literature has been on short-term localized forecasting. This is probably because the amount of data involved in a more global analysis is prohibitive with the canonical toolkit, via the Gaussian process (GP). Here we show how GP surrogate and discrepancy models can be combined to tractably and accurately predict solar irradiance on time-aggregated and daily scales with measurements at thousands of sites across the continental United States. Our results establish short term accuracy of bias-corrected weather-based simulation of irradiance, when realizations are available in real space-time (e.g., in future days), and provide accurate surrogates for smoothing in the more common situation where reliable weather data is not available (e.g., in future years).


title: Go with the Flow: Understanding Inflow Mechanisms in Galaxy Collisions

abstract: Dynamical interactions between colliding spiral galaxies strongly affect the state and distribution of their interstellar gas. Observations indicate that interactions funnel gas toward the nuclei, fueling bursts of star formation and nuclear activity. To date, most numerical simulations of galaxy mergers have assumed that the gaseous and stellar disks initially have the same distribution and size. However, observations of isolated disk galaxies show that this is seldom the case; in fact, most spirals have as much or more gas beyond their optical radii as they do within. Can gas in such extended disks be efficiently transported to the nuclei during interactions? To address this question, we examine the effect of various parameters on the transport of gas to the nuclei of interacting galaxies. In addition to the relative radii of the gaseous and stellar disks, these parameters include the pericentric separation, disk orientation, fractional gas mass, presence of a bulge, treatment of gas thermodynamics, and the spatial resolution of the numerical simulation. We found that gas accumulates in most of our simulated nuclei, but the efficiency of inflow is largely dependent upon the encounter geometry. Dissipation alone is not enough to produce inflows; an efficient mechanism for extracting angular momentum from the gas is necessary. Several different mechanisms are seen in these experiments. Aside from mode-driven inflows (such as, but not limited to, bars) and ram-pressure sweeping, both of which have been previously described and well studied, we supply the first quantitative study of an often-seen process: the formation of massive gas clumps in Jeans-unstable tidal shocks, and their subsequent delivery to the nuclei via dynamical friction.


title: Radial velocity survey of spatially resolved young, low-mass binaries

abstract: The identification and characterisation of low-mass binaries is of importance for a range of astrophysical investigations. Low-mass binaries in young ($\sim10 - 100 $ Myr) moving groups (YMGs) are of particular significance as they provide unique opportunities to calibrate stellar models and evaluate the ages and coevality of the groups themselves. Low-mass M-dwarfs have pre-main sequence life times on the order of $\sim100$ Myr and therefore are continually evolving along a mass-luminosity track throughout the YMG phase, providing ideal laboratories for precise isochronal dating, if a model-independent dynamical mass can be measured. AstraLux lucky imaging multiplicity surveys have recently identified hundreds of new YMG low-mass binaries, where a subsample of M-dwarf multiples have estimated orbital periods less than 50 years. We have conducted a radial velocity survey of a sample of 29 such targets to complement the astrometric data. This will allow enhanced orbital determinations and precise dynamical masses to be derived in a shorter timeframe than possible with astrometric monitoring alone, and allow for a more reliable isochronal analysis. Here we present radial velocity measurements for our sample over several epochs. We report the detection of the three-component spectroscopic multiple 2MASS J05301858-5358483, for which the C component is a new discovery, and forms a tight pair with the B component. Originally identified as a YMG member, we find that this system is a likely old field interloper, whose high chromospheric activity level is caused by tidal spin-up of the tight BC pair. Two other triple systems with a tight pair exist in the sample, but for the rest of the targets we find that additional tidally synchronized companions are highly unlikely, providing further evidence that their high chromospheric activity levels are generally signatures of youth.


title: Marginal Policy Gradients: A Unified Family of Estimators for Bounded Action Spaces with Applications

abstract: Many complex domains, such as robotics control and real-time strategy (RTS) games, require an agent to learn a continuous control. In the former, an agent learns a policy over $\mathbb{R}^d$ and in the latter, over a discrete set of actions each of which is parametrized by a continuous parameter. Such problems are naturally solved using policy based reinforcement learning (RL) methods, but unfortunately these often suffer from high variance leading to instability and slow convergence. Unnecessary variance is introduced whenever policies over bounded action spaces are modeled using distributions with unbounded support by applying a transformation $T$ to the sampled action before execution in the environment. Recently, the variance reduced clipped action policy gradient (CAPG) was introduced for actions in bounded intervals, but to date no variance reduced methods exist when the action is a direction, something often seen in RTS games. To this end we introduce the angular policy gradient (APG), a stochastic policy gradient method for directional control. With the marginal policy gradients family of estimators we present a unified analysis of the variance reduction properties of APG and CAPG; our results provide a stronger guarantee than existing analyses for CAPG. Experimental results on a popular RTS game and a navigation task show that the APG estimator offers a substantial improvement over the standard policy gradient.