title: On Analytical and Topological Properties of Separatrices in 1-D Holomorphic Dynamical Systems and Complex-Time Newton Flows

abstract: Separatrices divide the phase space of some holomorphic dynamical systems into separate basins of attraction or 'stability regions' for distinct fixed points. 'Bundling' (high density) and mutual 'repulsion' of trajectories are often observed at separatrices in phase portraits, but their global mathematical characterisation is a difficult problem. For 1-D complex polynomial dynamical systems we prove the existence of a separatrix for each critical point at infinity via transformation to the Poincar\'e sphere. We show that introduction of complex time allows a significantly extended view with the study of corresponding Riemann surface solutions, their topology, geometry and their bifurcations/ramifications related to separatrices. We build a bridge to the Riemann $\xi$-function and present a polynomial approximation of its Newton flow solution manifold with precision depending on the polynomial degree.


title: An integrodifference model for vegetation patterns in semi-arid environments with seasonality

abstract: Vegetation patterns are a characteristic feature of semi-deserts occurring on all continents except Antarctica. In some semi-arid regions, the climate is characterised by seasonality, which yields a synchronisation of seed dispersal with the dry season or the beginning of the wet season. We reformulate the Klausmeier model, a reaction-advection-diffusion system that describes the plant-water dynamics in semi-arid environments, as an integrodifference model to account for the temporal separation of plant growth processes during the wet season and seed dispersal processes during the dry season. The model further accounts for nonlocal processes involved in the dispersal of seeds. Our analysis focusses on the onset of spatial patterns. The Klausmeier partial differential equations (PDE) model is linked to the integrodifference model in an appropriate limit, which yields a control parameter for the temporal separation of seed dispersal events. We find that the conditions for pattern onset in the integrodifference model are equivalent to those for the continuous PDE model and hence independent of the time between seed dispersal events. We thus conclude that in the context of seed dispersal, a PDE model provides a sufficiently accurate description, even if the environment is seasonal. This emphasises the validity of results that have previously been obtained for the PDE model. Further, we numerically investigate the effects of changes to seed dispersal behaviour on the onset of patterns. We find that long-range seed dispersal inhibits the formation of spatial patterns and that the seed dispersal kernel's decay at infinity is a significant regulator of patterning.


title: On a Babu\v{s}ka paradox for polyharmonic operators: spectral stability and boundary homogenization for intermediate problems

abstract: We analyse the spectral convergence of high order elliptic differential operators subject to singular domain perturbations and homogeneous boundary conditions of intermediate type. We identify sharp assumptions on the domain perturbations improving, in the case of polyharmonic operators of higher order, conditions known to be sharp in the case of fourth order operators. The optimality is proved by analysing in detail a boundary homogenization problem, which provides a smooth version of a polyharmonic Babu\v{s}ka paradox.


title: On the robustness and performance of entropy stable discontinuous collocation methods for the compressible Navier-Stokes equations

abstract: In computational fluid dynamics, the demand for increasingly multidisciplinary reliable simulations, for both analysis and design optimization purposes, requires transformational advances in individual components of future solvers. At the algorithmic level, hardware compatibility and efficiency are of paramount importance in determining viability at exascale and beyond. However, equally important (if not more so) is algorithmic robustness with minimal user intervention, which becomes progressively more challenging to achieve as problem size and physics complexity increase. We numerically show that low and high order entropy stable discontinuous spatial discretizations based on summation-by-part operators and simultaneous-approximation-terms technique provides an essential step toward a truly enabling technology in terms of reliability and robustness for both under-resolved turbulent flow simulations and flows with discontinuities.


title: Forecasting Human-Object Interaction: Joint Prediction of Motor Attention and Actions in First Person Video

abstract: We address the challenging task of anticipating human-object interaction in first person videos. Most existing methods ignore how the camera wearer interacts with the objects, or simply consider body motion as a separate modality. In contrast, we observe that the international hand movement reveals critical information about the future activity. Motivated by this, we adopt intentional hand movement as a future representation and propose a novel deep network that jointly models and predicts the egocentric hand motion, interaction hotspots and future action. Specifically, we consider the future hand motion as the motor attention, and model this attention using latent variables in our deep model. The predicted motor attention is further used to characterise the discriminative spatial-temporal visual features for predicting actions and interaction hotspots. We present extensive experiments demonstrating the benefit of the proposed joint model. Importantly, our model produces new state-of-the-art results for action anticipation on both EGTEA Gaze+ and the EPIC-Kitchens datasets. Our project page is available at https://aptx4869lm.github.io/ForecastingHOI/


title: An Investigation on Semismooth Newton based Augmented Lagrangian Method for Image Restoration

abstract: Augmented Lagrangian method (also called as method of multipliers) is an important and powerful optimization method for lots of smooth or nonsmooth variational problems in modern signal processing, imaging, optimal control and so on. However, one usually needs to solve the coupled and nonlinear system together and simultaneously, which is very challenging. In this paper, we proposed several semismooth Newton methods to solve the nonlinear subproblems arising in image restoration, which leads to several highly efficient and competitive algorithms for imaging processing. With the analysis of the metric subregularities of the corresponding functions, we give both the global convergence and local linear convergence rate for the proposed augmented Lagrangian methods with semismooth Newton solvers.


title: Intertwining and the Markov uniqueness problem on path spaces

abstract: There are two open problem on the analysis of continuous paths on a Riemannian manifold, the Markov uniqueness and the independence of the closure of the differential operator $d$ on its initial domain. The operator $d$ acts naturally on $BC^1$ functions, one is concerned with its extensions to the $L^2$ spaces. With a suitable choice of an initial domain we denote by $D^{2,1}$ its closure under the graph norm. For the Wiener space, the domain of $d$ can be classified, as a consequence its extension is unique whether the initial domain is smooth cylindrical or in $BC^1$ etc. This has not shown to be the same when the measure is the probability distribution of any smooth elliptic diffusion. In an earlier paper, we have shown that the closure of $BC^\infty$ functions agree with that of smooth cylindrical functions, leaving an undesirable gap. The Markov uniqueness is essentially concerned with the problem whether there exists a unique Markov process on the path space whose Markov generator agrees with the infinite-dimensional Laplacian on $C^\infty$ cylindrical functions. Here we reduce Markov uniqueness to whether the pull back of $D^{2,1}$ by the ito map is $ D^{2,1}$ (i.e. a surjection). We also propose a possible approach for tackle this problem.


title: Periodic Orbits of Active Particles induced by Hydrodynamic Monopoles

abstract: Terrestrial experiments on active particles, such as Volvox, involve gravitational forces, torques and accompanying monopolar fluid flows. Taking these into account, we analyse the dynamics of a pair of self-propelling, self-spinning active particles between widely separated parallel planes. Neglecting flow reflected by the planes, the dynamics of orientation and horizontal separation is symplectic, with a Hamiltonian exactly determining limit cycle oscillations. Near the bottom plane, gravitational torque damps and reflected flow excites this oscillator, sustaining a second limit cycle that can be perturbatively related to the first. Our work provides a theory for dancing Volvox and highlights the importance of monopolar flow in active matter.


title: Formulae for the derivatives of heat semigroups

abstract: We use a basic martingale method to show a differentiation formula for the derivatives $$d(P_tf)(x_0)(v_0)={1\over t} E f(x_t) \int_0^t \langle Y(x_s)(v_s),dB_t\rangle_{R^m}.$$ These are proved first on $R^n$, then on manifolds. Afterwards for solutions of heat equations on differential forms, and a second order formula.


title: On unbiased simulations of stochastic bridges conditioned on extrema

abstract: Stochastic bridges are commonly used to impute missing data with a lower sampling rate to generate data with a higher sampling rate, while preserving key properties of the dynamics involved in an unbiased way. While the generation of Brownian bridges and Ornstein-Uhlenbeck bridges is well understood, unbiased generation of such stochastic bridges subject to a given extremum has been less explored in the literature. After a review of known results, we compare two algorithms for generating Brownian bridges constrained to a given extremum, one of which generalises to other diffusions. We further apply this to generate unbiased Ornstein-Uhlenbeck bridges and unconstrained processes, both constrained to a given extremum, along with more tractable numerical approximations of these algorithms. Finally, we consider the case of drift, and applications to geometric Brownian motions.


title: Nanobenders: efficient piezoelectric actuators for widely tunable nanophotonics at CMOS-level voltages

abstract: Tuning and reconfiguring nanophotonic components is needed to realize systems incorporating many components. The electrostatic force can deform a structure and tune its optical response. Despite the success of electrostatic actuators, they suffer from trade-offs between tuning voltage, tuning range, and on-chip area. Piezoelectric actuation could resolve all these challenges. Standard materials possess piezoelectric coefficients on the order of ${0.01}~\text{nm/V}$, suggesting extremely small on-chip actuation using potentials on the order of one volt. Here we propose and demonstrate compact piezoelectric actuators, called nanobenders, that transduce tens of nanometers per volt. By leveraging the non-uniform electric field from submicron electrodes, we generate bending of a piezoelectric nanobeam. Combined with a sliced photonic crystal cavity to sense displacement, we show tuning of an optical resonance by $\sim 5~\text{nm/V}~({0.6}~\text{THz/V})$ and between $1520$ and $1560~\text{nm}$ ($\sim 400$ linewidths) with only $ {4}~\text{V}$. Finally, we consider other tunable nanophotonic components enabled by nanobenders.


title: The local information of equivariant sheaves and elliptic difference equations

abstract: We study the singularities of algebraic difference equations on curves from the point of view of equivariant sheaves. We propose a definition for the formal local type of an equivariant sheaf at a point in the case of a reduced curve acted on by a group which is virtually the integers. We show that with this definition, equivariant sheaves can be glued from an "open cover". Precisely, we show that an equivariant sheaf can be uniquely recovered from the following data: the restriction to the complement of a point, the local type at the point itself, and an isomorphism between the two over the punctured neighborhood of said point. We study symmetric elliptic difference equations ("elliptic equations" from now on) from this point of view. We consider several natural notions for an algebraic version of symmetric elliptic difference equations, i.e. symmetric elliptic difference modules ("elliptic modules"). We show that different versions are not equivalent, but we detail how they are related: all the versions embed fully faithfully into the same category of equivariant sheaves. This implies that we can use the theory for equivariant sheaves to study singularities of elliptic equations as well. One reason to study elliptic equations is that they generalize, and degenerate to, (q-)difference equations (i.e. equivariant sheaves) and differential equations (i.e. D-modules) on the projective line. We discuss this from the elliptic module point of view, which requires studying elliptic modules on singular curves. We study the relation between elliptic modules on singular curves and their normalization. We show that for modules which are flat at the singular points there is an equivalence and we give examples showing that this cannot be improved upon.


title: Natural-mixing guided design of refractory high-entropy alloys with as-cast tensile ductility

abstract: Multi-principal-element metallic alloys have created a growing interest that is unprecedented in metallurgical history, in exploring the property limits of metals and the governing physical mechanisms. Refractory high-entropy alloys (RHEAs) have drawn particular attention due to their (i) high melting points and excellent softening-resistance, which are the two key requirements for high-temperature applications; and (ii) compositional space, which is immense even after considering cost and recyclability restrictions. However, RHEAs also exhibit intrinsic brittleness and oxidation-susceptibility, which remain as significant challenges for their processing and application. Here, utilizing natural-mixing characteristics amongst refractory elements, we designed a Ti38V15Nb23Hf24 RHEA that exhibits >20% tensile ductility already at the as-cast state, and physicochemical stability at high-temperatures. Exploring the underlying deformation mechanisms across multiple length-scales, we observe that a rare beta prime precipitation strengthening mechanism governs its intriguing mechanical response. These results also reveal the effectiveness of natural-mixing tendencies in expediting HEA discovery.


title: Intermolecular Cross-Correlations in the Dielectric Response of Glycerol

abstract: We suggest a way to disentangle self- from cross-correlation contributions in the dielectric spectra of glycerol. Recently it was demonstrated for monohydroxy alcohols that a detailed comparison of the dynamic susceptibilities of photon correlation and broadband dielectric spectroscopy allows to unambiguously disentangle a collective relaxation mode known as the Debye process, which could arises due to supramolecular structures, and the $\alpha$-relaxation, which proves to be identical in both methods. In the present paper, we apply the same idea and analysis to the paradigmatic glass former glycerol. For that purpose we present new light scattering data from photon correlation spectroscopy measurements and combine these with literature data to obtain a data set covering a dynamic range from $10^{-4}-10^{13}\,$Hz. Then we apply the above mentioned analysis by comparing this data set with a corresponding set of broadband dielectric data. Our finding is that even in a polyalcohol self- and cross-correlation contributions can approximately be disentangled in that way and that the emerging picture is very similar to that in monohydroxy alcohols. This is further supported by comparing the data with fast field cycling NMR measurements and dynamic shear relaxation data from the literature, and it turns out that, within the described approach, the $\alpha$-process appears very similar in all methods, while the pronounced differences observed in the spectral density are due to a different expression of the slow collective relaxational contribution. In the dielectric spectra the strength of this peak is reasonably well estimated by the Kirkwood correlation factor, which supports the view that it arises due to dynamic cross-correlations, which were previously often assumed to be negligible in dielectric measurements.


title: Mixed State Entanglement for Holographic Axion Model

abstract: We study the mixed state entanglement in a holographic axion model. We find that the holographic entanglement entropy (HEE), mutual information (MI) and entanglement of purification (EoP) exhibit very distinct behaviors with system parameters. The HEE exhibits universal monotonic behavior with system parameters, while the behaviors of MI and EoP relate to the specific system parameters and configurations. We find that MI and EoP can characterize mixed state entanglement better than HEE since they are less affected by thermal effects. Moreover, we argue that EoP is more suitable for describing mixed state entanglement than MI. Because the MI of large configurations are still dictated by the thermal entropy, while the EoP will never be controlled only by the thermal effects.


title: Modelling of Sickle Cell Anemia Patients Response to Hydroxyurea using Artificial Neural Networks

abstract: Hydroxyurea (HU) has been shown to be effective in alleviating the symptoms of Sickle Cell Anemia disease. While Hydroxyurea reduces the complications associated with Sickle Cell Anemia in some patients, others do not benefit from this drug and experience deleterious effects since it is also a chemotherapeutic agent. Therefore, to whom, should the administration of HU be considered as a viable option, is the main question asked by the responsible physician. We address this question by developing modeling techniques that can predict a patient's response to HU and therefore spare the non-responsive patients from the unnecessary effects of HU on the values of 22 parameters that can be obtained from blood samples in 122 patients. Using this data, we developed Deep Artificial Neural Network models that can predict with 92.6% accuracy, the final HbF value of a subject after undergoing HU therapy. Our current studies are focussing on forecasting a patient's HbF response, 30 days ahead of time.


title: Simple yet Effective Way for Improving the Performance of GAN

abstract: In adversarial learning, discriminator often fails to guide the generator successfully since it distinguishes between real and generated images using silly or non-robust features. To alleviate this problem, this brief presents a simple but effective way that improves the performance of generative adversarial network (GAN) without imposing the training overhead or modifying the network architectures of existing methods. The proposed method employs a novel cascading rejection (CR) module for discriminator, which extracts multiple non-overlapped features in an iterative manner using the vector rejection operation. Since the extracted diverse features prevent the discriminator from concentrating on non-meaningful features, the discriminator can guide the generator effectively to produce the images that are more similar to the real images. In addition, since the proposed CR module requires only a few simple vector operations, it can be readily applied to existing frameworks with marginal training overheads. Quantitative evaluations on various datasets including CIFAR-10, CelebA, CelebA-HQ, LSUN, and tiny-ImageNet confirm that the proposed method significantly improves the performance of GAN and conditional GAN in terms of Frechet inception distance (FID) indicating the diversity and visual appearance of the generated images.


title: Measurement of the Cosmic Microwave Background Polarization Lensing Power Spectrum from Two Years of POLARBEAR Data

abstract: We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.


title: Checking Chase Termination over Ontologies of Existential Rules with Equality

abstract: The chase is a sound and complete algorithm for conjunctive query answering over ontologies of existential rules with equality. To enable its effective use, we can apply acyclicity notions; that is, sufficient conditions that guarantee chase termination. Unfortunately, most of these notions have only been defined for existential rule sets without equality. A proposed solution to circumvent this issue is to treat equality as an ordinary predicate with an explicit axiomatisation. We empirically show that this solution is not efficient in practice and propose an alternative approach. More precisely, we show that, if the chase terminates for any equality axiomatisation of an ontology, then it terminates for the original ontology (which may contain equality). Therefore, one can apply existing acyclicity notions to check chase termination over an axiomatisation of an ontology and then use the original ontology for reasoning. We show that, in practice, doing so results in a more efficient reasoning procedure. Furthermore, we present equality model-faithful acyclicity, a general acyclicity notion that can be directly applied to ontologies with equality.


title: Formal Foundations of Continuous Graph Processing

abstract: With the growing need for online and iterative graph processing, software systems that continuously process large-scale graphs become widely deployed. With optimizations inherent as part of their design, these systems are complex, and have unique features beyond conventional graph processing. This paper describes CG Calculus, the first semantic foundation for continuous graph processing. The calculus captures the essential behavior of both the backend graph processing engine and the frontend application, with a focus on two essential features: temporal locality optimization (TLO) and incremental operation processing (IOP). A key design insight is that the operations continuously applied to the graph can be captured by a semantics defined over the operation stream flowing through the graph nodes. CG Calculus is a systematic study on the correctness of building continuous graph processing systems and applications. The most important result is result determinism: despite significant non-deterministic executions introduced by TLO and IOP, the results produced by CG Calculus are the same as conventional graph processing without TLO or IOP. The metatheory of CG Calculus is mechanized in Coq.


title: Light from an ion crystal: bunching or antibunching?

abstract: Photon statistics divides light sources into three different categories, characterized by bunched, antibunched or uncorrelated photon arrival times. Single atoms, ions, molecules, or solid state emitters display antibunching of photons, while classical thermal sources exhibit photon bunching. Here we demonstrate a light source in free space, where the photon statistics depends on the direction of observation, undergoing a continuous crossover between photon bunching and antibunching. We employ two trapped ions, observe their fluorescence under continuous laser light excitation, and record the spatially resolved autocorrelation function $g^{(2)}(\tau)$ with a movable Hanbury Brown and Twiss detector. Varying the detector position we find a minimum value for antibunching, $g^{(2)}(0) = 0.60(5)$ and a maximum of $g^{(2)}(0)=1.46(8)$ for bunching, demonstrating that this source radiates fundamentally different types of light alike. The observed variation of the autocorrelation function is understood in the Dicke model of heralded entangled states and the observed maximum and minimum values are modeled, taking independently measured experimental parameters into account.


title: Medium-band photometric reverberation mapping of AGNs at $0.1 < z < 0.8$. Techniques and sample

abstract: The most popular method of the broad-line region size estimation in active galactic nuclei (AGN) is the reverberation mapping based on measuring the time delay between the continuum flux and the flux in the emission lines. In our work, we apply the method of photometric reverberation mapping in mid-band filters, adapted for observations on the 1-m Zeiss-1000 telescope of Special Astrophysical Observatory of Russian Academy of Sciences, for the study of AGN with broad lines in the range of redshifts $0.1 < z < 0.8$. This paper provides a sample of 8 objects, describes the technique of observations and data processing for 2 studied objects to demonstrate the stability of the used method.


title: Manifestations of the electron-phonon interaction range in angle resolved photoemission spectra

abstract: Numerous angle resolved photoemission spectroscopy (ARPES) studies of a wide class of low-density metallic systems, ranging from doped transition metal oxides to quasi two-dimensional interfaces between insulators, exhibit phonon sidebands below the quasi-particle peak as a unique hallmark of polaronic correlations. Here, we single out properties of ARPES spectra that can provide a robust estimate of the effective range (screening length) of the electron-phonon interaction, regardless of the limited experimental resolution, dimensionality and particular features of the electronic structure, facilitating a general methodology for an analysis of a whole class of materials.


title: On Perfect Matchings in $k$-complexes

abstract: Keevash and Mycroft [\emph{Mem.~Amer.~Math.~Soc., 2015}] developed a geometric theory for hypergraph matchings and characterized the dense simplicial complexes that contain a perfect matching. Their proof uses the hypergraph regularity method and the hypergraph blow-up lemma recently developed by Keevash. In this note we give a new proof of their results, which avoids these complex tools. In particular, our proof uses the lattice-based absorbing method developed by the author and a recent probabilistic argument of Kohayakawa, Person and the author.


title: Dynamical quantum phase transition in a system of non-interacting bosons

abstract: We study a Bose-Einstein condensate at the low energy limit and show that their collective dynamics exhibit interesting quantum dynamical behavior. The system undergoes a dynamical quantum phase transition after a sudden quench into a properly distributed static force, provided its dispersion relation is linear. We corroborate the occurrence of the dynamical phase transition by calculating Fisher zeros of the Loschmidt amplitude and showing that they cross the real time axis in thermodynamic limit. A connection is established between the order of nonanalycity in the return rate function and the spectral density of the force. Furthermore, it is shown that a logarithmic or power law scaling holds at the critical times depending on the displacement spectrum. The scaling behaviors are studied for three different cases. Eventually, a scheme for the quantum simulation of such dynamical phase transition and its verification is proposed. We show that the behavior remains observable even by taking into account the environmental effects.


title: Sparsity through evolutionary pruning prevents neuronal networks from overfitting

abstract: Modern Machine learning techniques take advantage of the exponentially rising calculation power in new generation processor units. Thus, the number of parameters which are trained to resolve complex tasks was highly increased over the last decades. However, still the networks fail - in contrast to our brain - to develop general intelligence in the sense of being able to solve several complex tasks with only one network architecture. This could be the case because the brain is not a randomly initialized neural network, which has to be trained by simply investing a lot of calculation power, but has from birth some fixed hierarchical structure. To make progress in decoding the structural basis of biological neural networks we here chose a bottom-up approach, where we evolutionarily trained small neural networks in performing a maze task. This simple maze task requires dynamical decision making with delayed rewards. We were able to show that during the evolutionary optimization random severance of connections lead to better generalization performance of the networks compared to fully connected networks. We conclude that sparsity is a central property of neural networks and should be considered for modern Machine learning approaches.


title: Microscopy Image Restoration with Deep Wiener-Kolmogorov filters

abstract: Microscopy is a powerful visualization tool in biology, enabling the study of cells, tissues, and the fundamental biological processes; yet, the observed images typically suffer from blur and background noise. In this work, we propose a unifying framework of algorithms for Gaussian image deblurring and denoising. These algorithms are based on deep learning techniques for the design of learnable regularizers integrated into the Wiener-Kolmogorov filter. Our extensive experimentation line showcases that the proposed approach achieves a superior quality of image reconstruction and surpasses the solutions that rely either on deep learning or on optimization schemes alone. Augmented with the variance stabilizing transformation, the proposed reconstruction pipeline can also be successfully applied to the problem of Poisson image deblurring, surpassing the state-of-the-art methods. Moreover, several variants of the proposed framework demonstrate competitive performance at low computational complexity, which is of high importance for real-time imaging applications.


title: Against 'Particle Metaphysics' and 'Collapses' within the Definition of Quantum Entanglement

abstract: In this paper we argue against the orthodox definition of quantum entanglement which has been explicitly grounded on several "common sense" (metaphysical) presuppositions and presents today serious formal and conceptual drawbacks. This interpretation which some researchers in the field call "minimal", has ended up creating a narrative according to which QM talks about "small particles" represented by pure states (in general, superpositions) which --each time someone attempts to observe what is going on-- suddenly "collapse" to a single measurement outcome. After discussing the consequences of applying 'particle metaphysics' within the definition of entanglement we turn our attention to two recent approaches which might offer an interesting way out of this metaphysical conundrum. Both approaches concentrate their efforts in going beyond the notion of 'system'. While the first, called device-independent approach, proposes an operational anti-metaphysical scheme in which language plays an essential role; the second approach, takes an essentially metaphysical path which attempts to present a new non-classical representation grounded on intensive relations which, in turn, imposes the need to reconsider the definition and meaning of quantum entanglement.


title: Accurate Computations of Euler Products over Primes in Arithmetic Progressions

abstract: This note provides truncated formulae with explicit error terms to compute Euler products over primes in arithmetic progressions of rational fractions. It further provides such a formula for the product of terms of the shape $F(1/p, 1/p^s)$ when $F$ is a two-variable polynomial with coefficients in $\mathbb{C}$ and satisfying some restrictive conditions.


title: On 3-Hom-Lie-Rinehart algebras

abstract: We introduce the notion of 3-Hom-Lie-Rinehart algebra and systematically describe a cohomology complex by considering coefficient modules. Furthermore, we consider extensions of a 3-Hom-Lie-Rinehart algebra and characterize the first cohomology space in terms of the group of automorphisms of an $A$-split abelian extension and the equivalence classes of $A$-split abelian extensions. Finally, we study formal deformations of 3-Hom-Lie-Rinehart algebras.


title: C*-algebras generated by multiplication operators and composition operators by functions with self-similar branches

abstract: Let $K$ be a compact metric space and let $\varphi: K \to K$ be continuous. We study C*-algebra $\mathcal{MC}_\varphi$ generated by all multiplication operators by continuous functions on $K$ and a composition operator $C_\varphi$ induced by $\varphi$ on a certain $L^2$ space. Let $\gamma = (\gamma_1, \dots, \gamma_n)$ be a system of proper contractions on $K$. Suppose that $\gamma_1, \dots, \gamma_n$ are inverse branches of $\varphi$ and $K$ is self-similar. We consider the Hutchinson measure $\mu^H$ of $\gamma$ and the $L^2$ space $L^2(K, \mu^H)$. Then we show that the C*-algebra $\mathcal{MC}_\varphi$ is isomorphic to the C*-algebra $\mathcal{O}_\gamma (K)$ associated with $\gamma$ under some conditions.


title: PETALO: Time-of-Flight PET with liquid xenon

abstract: The fast scintillation decay time and the high scintillation yield of liquid xenon makes it an appropriate material for nuclear medicine. Moreover, being a continuous medium with a uniform response, liquid xenon allows one to avoid most of the geometrical distortions of conventional detectors based on scintillating crystals. In this paper, we describe how these properties have motivated the development of a novel concept for positron emission tomography scanners with Time-Of-Flight measurement, which uses liquid xenon as a scintillating material and silicon photomultipliers as sensors. Monte Carlo studies have indicated that this technology would provide a very good intrinsic time resolution, of around 70 ps. Moreover, being liquid xenon transparent to UV and blue wavelengths, both scintillation and Cherenkov light can be exploited. While the former can be used for energy measurements, the latter is a prompt signal (of a few picoseconds), which provides a very precise time measurement. Monte Carlo simulations point to a time resolution of 30-50 ps obtained using Cherenkov light. A first prototype is being built to demonstrate the high energy, spatial and time resolution of this concept, using a ring of around 30 cm of internal diameter and a depth of 3 cm instrumented with VUV--sensitive silicon photomultipliers.


title: DE/RM-MEDA: A New Hybrid Multi-Objective Generator

abstract: Under the condition of Karush-Kuhn-Tucker, the Pareto Set (PS) in the decision area of an m-objective optimization problem is a piecewise continuous (m-1)-D manifold. For illustrate the degree of convergence of the population, we employed the ratio of the sum of the first (m-1) largest eigenvalue of the population's covariance matrix of the sum of all eigenvalue. Based on this property, this paper proposes a new algorithm, called DE/RM-MEDA, which mix differential evolutionary (DE) and the estimation of distribution algorithm (EDA) to generate and adaptively adjusts the number of new solutions by the ratio. The proposed algorithm is experimented on nine tec09 problems. The comparison results between DE/RM-MEDA and the others algorithms, called NSGA-II-DE and RM-MEDA, show that the proposed algorithm perform better in terms of convergence and diversity metric.


title: Strong cosmic censorship for a scalar field in an Einstein-Maxwell-Gauss-Bonnet-de Sitter black hole

abstract: It has been shown that the Christodoulou version of the Strong Cosmic Censorship (SCC) conjecture can be violated for a scalar field in a near-extremal Reissner-Nordstrom-de Sitter black hole. In this paper, we investigate the effects of higher derivative corrections to the Einstein-Hilbert action on the validity of SCC, by considering a neutral massless scalar perturbation in 5- and 6-dimensional Einstein-Maxwell-Gauss-Bonnet-de Sitter black holes. Our numerical results exhibit that the higher derivative term plays a different role in the d = 5 case than it does in the d = 6 case. For d = 5, the SCC violation region increases as the strength of the higher derivative term increases. For d = 6, the SCC violation region first increases and then decreases as the higher derivative correction becomes stronger, and SCC can always be restored for a black hole with a fixed charge ratio when the higher derivative correction is strong enough. Finally, we find that the C2 version of SCC is respected in the d = 6 case, but can be violated in some near-extremal regime in the d = 5 case.


title: 6.7 GHz CH3OH absorption towards the N3 Galactic Center point-source

abstract: We present evidence of 6.7 GHz methanol (CH3OH)~and 4.8 GHz formaldehyde (H2CO) absorption towards the Galactic Center (GC) point-source `N3'. Both absorption features are unresolved and spatially aligned with N3. The 6.7 GHz CH3OH contains a single velocity component (centered at ~10 km/s) while the 4.8 GHz H2CO shows two velocity components (centered at ~-3 and +8 km/s). We find that the velocity of these absorption components are similar to that of emission lines from other molecules (e.g., SiO and H3CN) detected toward this compact-source (-13 to +25 km/s; `N3 cloud'). The detection of these absorption features is a firm indication that some of the molecular gas in the N3 molecular cloud is on the near-side of the continuum source. Analysis of the CH3OH absorption kinematics shows a relatively large velocity dispersion (3.8 km/s) for the size scale of this feature (<0.1'', <0.01 pc at the GC; Ludovici et al. 2016), when compared with other similar size GC clouds in the Larson linewidth-size relationship. Further, this linewidth is closer to velocity dispersion measurements for size scales of 1.3 pc, which is roughly the width of the N3 cloud (25''; 1.0 pc). We argue that this relatively broad linewidth, over a small cross-sectional area, is due to turbulence through the depth of the cloud, where the cloud has a presumed line-of-sight thickness of ~1 pc.


title: Statistical analysis of read-back signals in magnetic recording on granular media

abstract: The comprehensive simulation of magnetic recording, including the write and read-back process, on granular media becomes computationally expensive if the magnetization dynamics of each grain are explicitly computed. In addition, in heat-assisted magnetic recording, the writing of a single track becomes a random process since the temperature must be considered and thermal noise is involved. Further, varying grain structures of various granular media must also be taken into account to obtain correct statistics for the final read-back signal. Hence, it requires many repetitions of the write process to investigate the mean signal as well as the noise. This work presents a method that improves the statistical evaluation of the whole recording process. The idea is to avoid writing the magnetization to one of its binary states. Instead, we assign each grain its probability of occupying one of its stable states, which can be calculated in advance in terms of a switching probability phase diagram. In the read-back process, we combine the probabilities to calculate a mean signal and its variance. Afterwards, repetitions on different media lead to the final read-back signal. Using a recording example, we show that the statistical behavior of the evaluated signal-to-noise ratio can be significantly improved by applying this probability mapping method, while the computational effort remains low.


title: JKO estimates in linear and non-linear Fokker-Planck equations, and Keller-Segel: L p and Sobolev bounds

abstract: We analyze some parabolic PDEs with different drift terms which are gradient flows in the Wasserstein space and consider the corresponding discrete-in-time JKO scheme. We prove with optimal transport techniques how to control the L p and L $\infty$ norms of the iterated solutions in terms of the previous norms, essentially recovering well-known results obtained on the continuous-in-time equations. Then we pass to higher order results, and in particulat to some specific BV and Sobolev estimates, where the JKO scheme together with the so-called "five gradients inequality" allows to recover some inequalities that can be deduced from the Bakry-Emery theory for diffusion operators, but also to obtain some novel ones, in particular for the Keller-Segel chemiotaxis model. 1 Short introduction The goal of this paper is to present some estimates on evolution PDEs in the space of probability densities which share two important features: they include a linear diffusion term, and they are gradient flows in the Wasserstein space W2. These PDEs will be of the form $\partial$t$\rho$ -- $\Delta$$\rho$ -- $\nabla$ $\times$ ($\rho$$\nabla$u[$\rho$]) = 0, complemented with no-flux boundary conditions and an intial condition on $\rho$0. We will in particular concentrate on the Fokker-Plack case, where u[$\rho$] = V and V is a fixed function (with possible regularity assumptions) independent of $\rho$, on the case where u[$\rho$] = W * $\rho$ is obtained by convolution and models interaction between particles, and on the parabolic-elliptic Keller-Segel case where u[$\rho$] is related to $\rho$ via an elliptic equation. This last case models the evolution of a biological population $\rho$ subject to diffusion but attracted by the concentration of a chemo-attractant, a nutrient which is produced by the population itself, so that its distribution is ruled by a PDE where the density $\rho$ appears as a source term. Under the assumption that the production rate of this nutrient is much faster than the motion of the cells, we can assume that its distribution is ruled by a statical PDE with no explicit time-dependence, and gives rise to a system which is a gradient flow in the variable $\rho$ (the parabolic-parabolic case, where the time scale for the cells and for the nutrient are comparable, is also a gradient flow, in the product space W2 x L 2 , but we will not consider this case). Since we mainly concentrate on the case of bounded domains, in the Keller-Segel case the term u[$\rho$] cannot be expressed as a convoluton and requires ad-hoc computations. In all the paper, the estimates will be studied on a time-discretized version of these PDEs, consisting in the so-called JKO (Jordan-Kinderleherer-Otto) scheme, based on iterated optimization problems involving the Wasserstein distance W2. We will first present 0-order estimates, on the L p and L $\infty$ norms of the solution. This is just a translation into the JKO language of well-known properties of these equations. The main goal of this part is hence to


title: The Convex Information Bottleneck Lagrangian

abstract: The information bottleneck (IB) problem tackles the issue of obtaining relevant compressed representations $T$ of some random variable $X$ for the task of predicting $Y$. It is defined as a constrained optimization problem which maximizes the information the representation has about the task, $I(T;Y)$, while ensuring that a certain level of compression $r$ is achieved (i.e., $ I(X;T) \leq r$). For practical reasons, the problem is usually solved by maximizing the IB Lagrangian (i.e., $\mathcal{L}_{\text{IB}}(T;\beta) = I(T;Y) - \beta I(X;T)$) for many values of $\beta \in [0,1]$. Then, the curve of maximal $I(T;Y)$ for a given $I(X;T)$ is drawn and a representation with the desired predictability and compression is selected. It is known when $Y$ is a deterministic function of $X$, the IB curve cannot be explored and another Lagrangian has been proposed to tackle this problem: the squared IB Lagrangian: $\mathcal{L}_{\text{sq-IB}}(T;\beta_{\text{sq}})=I(T;Y)-\beta_{\text{sq}}I(X;T)^2$. In this paper, we (i) present a general family of Lagrangians which allow for the exploration of the IB curve in all scenarios; (ii) provide the exact one-to-one mapping between the Lagrange multiplier and the desired compression rate $r$ for known IB curve shapes; and (iii) show we can approximately obtain a specific compression level with the convex IB Lagrangian for both known and unknown IB curve shapes. This eliminates the burden of solving the optimization problem for many values of the Lagrange multiplier. That is, we prove that we can solve the original constrained problem with a single optimization.


title: Riesz bases of reproducing kernels in small Fock spaces

abstract: We give a complete characterization of Riesz bases of normalized reproducing kernels in the small Fock spaces $\mathcal{F}^2_{\varphi}$, the spaces of entire functions $f$ such that $f\mathrm{e}^{-\varphi} \in L^{2}(\mathbb{C})$, where $\varphi(z)= (\log^+|z|)^{\beta+1}$, $0< \beta \leq 1$.The first results in this direction are due to Borichev-Lyubarskii who showed that $\varphi$ with $\beta=1$ is the largest weight for which the corresponding Fock space admits Riesz bases of reproducing kernels. Later, such bases were characterized by Baranov-Dumont-Hartman-Kellay in the case when $\beta=1$. The present paper answers a question in Baranov et al. by extending their results for all parameters $\beta\in (0,1)$. Our results are analogous to those obtained for the case $\beta=1$ and those proved for Riesz bases of complex exponentials for the Paley-Wiener spaces. We also obtain a description of complete interpolating sequences in small Fock spaces with corresponding uniform norm.


title: ForestFit : An R package for modeling tree diameter distributions

abstract: Modeling the diameter distribution of trees in forest stands is a common forestry task that supports key biologically and economically relevant management decisions. The choice of model used to represent the diameter distribution and how to estimate its parameters has received much attention in the forestry literature; however, accessible software that facilitates comprehensive comparison of the myriad modeling approaches is not available. To this end, we developed an R package called ForestFit that simplifies estimation of common probability distributions used to model tree diameter distributions, including the two- and three-parameter Weibull distributions, Johnson's SB distribution, Birnbaum-Saunders distribution, and finite mixture distributions. Frequentist and Bayesian techniques are provided for individual tree diameter data, as well as grouped data. Additional functionality facilitates fitting growth curves to height-diameter data. The package also provides a set of functions for computing probability distributions and simulating random realizations from common finite mixture models.


title: A Model System of Mixed Ionized Gas Dynamics

abstract: The aim of this paper is to study a one dimensional model system of equations for ionized gas dynamics at high temperature where the gas is a mixture of two kinds of monatomic gas. In addition to the mass density, pressure, temperature and particle velocity, degrees of ionization of both gases are also involved. By assuming that the local thermal equilibrium is attained, Saha's ionization equations are added. Thus the equations are supplemented by the first and second law of thermodynamics, a single equation of state and, in addition, a set of thermodynamic equations. The equations constitute a strictly hyperbolic system, which guarantees that the initial value problem is well-posed locally in time for sufficiently smooth initial data. However the geometric properties of the system are rather complicated: in particular, we prove the existence of a region where convexity (genuine nonlinearity) fails for forward and backward characteristic fields. Also we study thermodynamic properties of shock waves by a detailed analysis of the Hugoniot locus, which is used in a mathematical study of existence and uniqueness of solutions to the shock tube problem.


title: Factorization and malleability of RSA modules, and counting points on elliptic curves modulo N

abstract: In this paper we address two different problems related with the factorization of an RSA module N. First we can show that factoring is equivalent in deterministic polynomial time to counting points on a pair of twisted Elliptic curves modulo N. Also we settle the malleability of factoring an RSA module, as described in [9], using the number of points of a single elliptic curve modulo N, and Coppersmith's algorithm.


title: Greedy Algorithms for Fair Division of Mixed Manna

abstract: We consider a multi-agent model for fair division of mixed manna (i.e. items for which agents can have positive, zero or negative utilities), in which agents have additive utilities for bundles of items. For this model, we give several general impossibility results and special possibility results for three common fairness concepts (i.e. EF1, EFX, EFX3) and one popular efficiency concept (i.e. PO). We also study how these interact with common welfare objectives such as the Nash, disutility Nash and egalitarian welfares. For example, we show that maximizing the Nash welfare with mixed manna (or minimizing the disutility Nash welfare) does not ensure an EF1 allocation whereas with goods and the Nash welfare it does. We also prove that an EFX3 allocation may not exist even with identical utilities. By comparison, with tertiary utilities, EFX and PO allocations, or EFX3 and PO allocations always exist. Also, with identical utilities, EFX and PO allocations always exist. For these cases, we give polynomial-time algorithms, returning such allocations and approximating further the Nash, disutility Nash and egalitarian welfares in special cases.


title: On the Problem of Infinite Spin in Total Collisions of the Planar N-Body Problem

abstract: For the planar $N$-body problem, we first introduce a class of moving frame suitable for orbits near central configurations, especially for total collision orbits, which is the main new ingredient of this paper. The moving frame allows us to reduce the degeneracy of the problem according to intrinsic symmetrical characteristic of the $N$-Body problem. First, we give a full answer to the infinite spin or $Painlev\acute{e}$-$Wintner$ problem in the case corresponding to nondegenerate central configurations. Then following some original ideas of C.L. Siegel, especially the idea of normal forms, and applying the theory of central manifolds, we give a partial answer to the problem in the case corresponding to degenerate central configurations. We completely answer the problem in the case corresponding to central configurations with one degree of degeneracy. Combining some results on the planar nonhyperbolic equilibrium point, we give a criterion for the case corresponding to central configurations with two degrees of degeneracy. We further answer the problem in the case corresponding to all known degenerate central configurations of four-body. Therefore, we solve the problem for almost every choice of the masses of the four-body problem. Finally, we give a measure of the set of initial conditions leading to total collisions.


title: Negative-hydrogen-ion production from a nanoporous 12CaO $\cdot$ 7Al2O3 (C12A7) electride surface

abstract: A high production rate of negative hydrogen ion was observed from a nanoporous C12A7 electride surface immersed in hydrogen/deuterium low-pressure plasmas. The target was biased at 10-170 V negatively and the target surface was bombarded by H3 + ions from the plasma. The production rate was 1-2 orders higher than a clean molybdenum surface. The measured H-energy spectrum indicates that the major production mechanism is desorption by sputtering. This material has potential to be used as production surface of cesium-free negative ion sources for accelerators, heating beams in nuclear fusion, and surface modification for industrial applications.


title: Korean-to-Chinese Machine Translation using Chinese Character as Pivot Clue

abstract: Korean-Chinese is a low resource language pair, but Korean and Chinese have a lot in common in terms of vocabulary. Sino-Korean words, which can be converted into corresponding Chinese characters, account for more than fifty of the entire Korean vocabulary. Motivated by this, we propose a simple linguistically motivated solution to improve the performance of the Korean-to-Chinese neural machine translation model by using their common vocabulary. We adopt Chinese characters as a translation pivot by converting Sino-Korean words in Korean sentences to Chinese characters and then train the machine translation model with the converted Korean sentences as source sentences. The experimental results on Korean-to-Chinese translation demonstrate that the models with the proposed method improve translation quality up to 1.5 BLEU points in comparison to the baseline models.


title: Efficient energy resolved quantum master equation for transport calculations in large strongly correlated systems

abstract: We introduce a systematic approximation for an efficient evaluation of Born--Markov master equations for steady state transport studies in open quantum systems out of equilibrium: the energy resolved master equation approach. The master equation is formulated in the eigenbasis of the open quantum system and build successively by including eigenstates with increasing grandcanonical energies. In order to quantify convergence of the approximate scheme we introduce quality factors to check preservation of trace, positivity and hermiticity. Furthermore, we discuss different types of master equations that go beyond the commonly used secular approximation in order to resolve coherences between quasi--degenerate states. For the discussion of complete positivity we introduce a canonical Redfield-Bloch master equation and compare it to a previously derived master equations in Lindblad form with and without using the secular approximation. The approximate scheme is benchmarked for a six orbital quantum system which shows destructive quantum interference under the application of a bias voltage. The energy resolved master equation approach presented here makes quantum transport calculations in many--body quantum systems numerically accessible also beyond six orbitals with a full Hilbert space of the order of $\sim 10^6$.


title: Event Recognition with Automatic Album Detection based on Sequential Processing, Neural Attention and Image Captioning

abstract: In this paper a new formulation of event recognition task is examined: it is required to predict event categories in a gallery of images, for which albums (groups of photos corresponding to a single event) are unknown. We propose the novel two-stage approach. At first, features are extracted in each photo using the pre-trained convolutional neural network. These features are classified individually. The scores of the classifier are used to group sequential photos into several clusters. Finally, the features of photos in each group are aggregated into a single descriptor using neural attention mechanism. This algorithm is optionally extended to improve the accuracy for classification of each image in an album. In contrast to conventional fine-tuning of convolutional neural networks (CNN) we proposed to use image captioning, i.e., generative model that converts images to textual descriptions. They are one-hot encoded and summarized into sparse feature vector suitable for learning of arbitrary classifier. Experimental study with Photo Event Collection and Multi-Label Curation of Flickr Events Dataset demonstrates that our approach is 9-20% more accurate than event recognition on single photos. Moreover, proposed method has 13-16% lower error rate than classification of groups of photos obtained with hierarchical clustering. It is experimentally shown that the image captions trained on Conceptual Captions dataset can be classified more accurately than the features from object detector, though they both are obviously not as rich as the CNN-based features. However, it is possible to combine our approach with conventional CNNs in an ensemble to provide the state-of-the-art results for several event datasets.


title: Response relation in Pb-Pb and p-Pb collisions at 5.02 TeV

abstract: We carry out simulations using a multi-phase transport (AMPT) model to describe the response relation between $v_2$ and $\varepsilon_2$ in Pb-Pb and p-Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV, respectively. To simulate such relation, two methods have been introduced in the calculation: one is the directed response (DR) method, which correlates the outgoing particles with the initial anisotropy directly, and the other one is the cumulants response (CR) method, which is constructed from a cumulants correlation between outgoing particles. Based on calculations of the DR and CR methods, the response relations as a function of the transverse momentum are both shown in Pb-Pb and p-Pb collisions. By comparing the DR and CR methods, we found that the linear response relations are almost identical in all the present collisions. Similar results of linear+cubic response relations are also shown in the higher multiplicity systems, and it has become a significant difference in the lower multiplicity systems, i.e., the peripheral Pb-Pb collisions and p-Pb collisions. Throughout the whole $p_{T}$-dependent simulations, the $\kappa_2$ in the linear response and in the linear+cubic response are almost identical, except for in the lower multiplicity systems by the DR method. If one implements a pseudorapidity gap by the CR calculation, the $p_{T}$-dependent and $\eta$-independent response relations are similarly shown in peripheral Pb-Pb systems and p-Pb systems, which may imply that a collective response exists in the most central p-Pb collisions. These collective behaviors are dominantly produced on the stage of the medium expansions.


title: A Stabilizing Control Algorithm for Asynchronous Parallel Quadratic Programming via Dual Decomposition

abstract: This paper proposes a control algorithm for stable implementation of asynchronous parallel quadratic programming (PQP) through dual decomposition technique. In general, distributed and parallel optimization requires synchronization of data at each iteration step due to the interdependency of data. The synchronization latency may incur a large amount of waiting time caused by an idle process during computation. We aim to mitigate this synchronization penalty in PQP problems by implementing asynchronous updates of the dual variable. The price to pay for adopting asynchronous computing algorithms is the unpredictability of the solution, resulting in a tradeoff between speedup and accuracy. In the worst case, the state of interest may become unstable owing to the stochastic behavior of asynchrony. We investigate the stability condition of asynchronous PQP problems by employing the switched system framework. A formal algorithm is provided to ensure the asymptotic stability of dual variables. Further, it is shown that the implementation of the proposed algorithm guarantees the uniqueness of optimal solutions, irrespective of asynchronous behavior. To verify the validity of the proposed methods, simulation results are presented.


title: On the analytical solution in non-inertial frame of R2BP

abstract: In this analytical study, we have presented a new type of solving procedure with aim to obtain the coordinates of small mass m, which moves around primary M_Sun, referred to non-inertial frame of restricted two-body problem (R2BP) with modified potential function (taking into account the components of variable velocity of central body M_Sun motion) instead of classical potential function for Kepler formulation of R2BP. Meanwhile, system of equations of motion has been successfully explored with respect to the existence of analytical way for presentation of the solution in polar coordinates with radial distance r = r(t). We have obtained analytical formula for function t = t(r) via appropriate elliptic integral. Having obtained the inversed dependence r = r(t), we can obtain the time-dependence for the polar angle as well. Also, we have pointed out how to express components of solution (including initial conditions) from cartesian to polar coordinates.


title: The Batchelor spectrum of passive scalar turbulence in stochastic fluid mechanics at fixed Reynolds number

abstract: In 1959, Batchelor predicted that the stationary statistics of passive scalars advected in fluids with small diffusivity $\kappa$ should display a $|k|^{-1}$ power spectrum along an inertial range contained in the viscous-convective range of the fluid model. This prediction has been extensively tested, both experimentally and numerically, and is a core prediction of passive scalar turbulence. In this article we provide a rigorous proof of a version of Batchelor's prediction in the $\kappa \to 0$ limit when the scalar is subjected to a spatially-smooth, white-in-time stochastic source and is advected by the 2D Navier-Stokes equations or 3D hyperviscous Navier-Stokes equations in $\mathbb{T}^d$ forced by sufficiently regular, nondegenerate stochastic forcing. Although our results hold for fluids at arbitrary Reynolds number, this value is fixed throughout. Our results rely on the quantitative understanding of Lagrangian chaos and passive scalar mixing established in our recent works. Additionally, in the $\kappa \to 0$ limit, we obtain statistically stationary, weak solutions in $H^{-\epsilon}$ to the stochastically-forced advection problem without diffusivity. These solutions are almost-surely not locally integrable distributions with non-vanishing average anomalous flux and satisfy the Batchelor spectrum at all sufficiently small scales. We also prove an Onsager-type criticality result which shows that no such dissipative, weak solutions with a little more regularity can exist.


title: Supersymmetric localization, modularity and the Witten genus

abstract: Equivariant localization techniques give a rigorous interpretation of the Witten genus as an integral over the double loop space. This provides a geometric explanation for its modularity properties. It also reveals an interplay between the geometry of double loop spaces and complex analytic elliptic cohomology. In particular, we identify a candidate target for the elliptic Bismut-Chern character.


title: Persistent Magnitude

abstract: In this paper we introduce the persistent magnitude, a new numerical invariant of (sufficiently nice) graded persistence modules. It is a weighted and signed count of the bars of the persistence module, in which a bar of the form $[a,b)$ in degree $d$ is counted with weight $(e^{-a}-e^{-b})$ and sign $(-1)^d$. Persistent magnitude has good formal properties, such as additivity with respect to exact sequences and compatibility with tensor products, and has interpretations in terms of both the associated graded functor, and the Laplace transform. Our definition is inspired by Otter's notion of blurred magnitude homology: we show that the magnitude of a finite metric space is precisely the persistent magnitude of its blurred magnitude homology. Turning this result on its head, we obtain a strategy for turning existing persistent homology theories into new numerical invariants by applying the persistent magnitude. We explore this strategy in detail in the case of persistent homology of Morse functions, and in the case of Rips homology.


title: An End-to-End Framework for Cold Question Routing in Community Question Answering Services

abstract: Routing newly posted questions (a.k.a cold questions) to potential answerers with the suitable expertise in Community Question Answering sites (CQAs) is an important and challenging task. The existing methods either focus only on embedding the graph structural information and are less effective for newly posted questions, or adopt manually engineered feature vectors that are not as representative as the graph embedding methods. Therefore, we propose to address the challenge of leveraging heterogeneous graph and textual information for cold question routing by designing an end-to-end framework that jointly learns CQA node embeddings and finds best answerers for cold questions. We conducted extensive experiments to confirm the usefulness of incorporating the textual information from question tags and demonstrate that an end-2-end framework can achieve promising performances on routing newly posted questions asked by both existing users and newly registered users.


title: A Self-Adaptive Synthetic Over-Sampling Technique for Imbalanced Classification

abstract: Traditionally, in supervised machine learning, (a significant) part of the available data (usually 50% to 80%) is used for training and the rest for validation. In many problems, however, the data is highly imbalanced in regard to different classes or does not have good coverage of the feasible data space which, in turn, creates problems in validation and usage phase. In this paper, we propose a technique for synthesising feasible and likely data to help balance the classes as well as to boost the performance in terms of confusion matrix as well as overall. The idea, in a nutshell, is to synthesise data samples in close vicinity to the actual data samples specifically for the less represented (minority) classes. This has also implications to the so-called fairness of machine learning. In this paper, we propose a specific method for synthesising data in a way to balance the classes and boost the performance, especially of the minority classes. It is generic and can be applied to different base algorithms, e.g. support vector machine, k-nearest neighbour, deep networks, rule-based classifiers, decision trees, etc. The results demonstrated that: i) a significantly more balanced (and fair) classification results can be achieved; ii) that the overall performance as well as the performance per class measured by confusion matrix can be boosted. In addition, this approach can be very valuable for the cases when the number of actual available labelled data is small which itself is one of the problems of the contemporary machine learning.


title: Influence of Fermions on Vortices in SU(2)-QCD

abstract: Gauge fields control the dynamics of fermions, also a back reaction of fermions on the gauge field is expected. This back reaction is investigated within the vortex picture of the QCD vacuum. We show that the center vortex model reproduces the string tension of the full theory also with the presence of fermionic fields.


title: Fractional Hypocoercivity

abstract: This paper is devoted to kinetic equations without confinement. We investigate the large time behaviour induced by collision operators with fat tailed local equilibria. Such operators have an anomalous diffusion limit. In the appropriate scaling, the macroscopic equation involves a fractional diffusion operator so that the optimal decay rate is determined by a fractional Nash type inequality. At kinetic level we develop an $\mathrm L^2$-hypocoercivity approach and establish a rate of decay compatible with the fractional diffusion limit.


title: Recent advances in the description of reactions involving exotic nuclei

abstract: In this contribution to the proceedings of the International Nuclear Physics Conference 2019, I review recent developments made in reaction models used to analyse data measured at radioactive-ion beam facilities to study exotic nuclear structures. I focus in particular on reactions like elastic scattering and breakup, which are used to study halo nuclei. Although these peculiar nuclei challenge usual nuclear-structure models, some can now be computed ab initio. This brief review illustrates the progresses made in nuclear-reaction theory in the last few years to improve the description of the projectile within reaction models. I dedicate this contribution to the memory of Mahir Hussein, who has significantly contributed to this field and who passed away in May this year.


title: Metastability as a coexistence mechanism in a model for dryland vegetation patterns

abstract: Vegetation patterns are a ubiquitous feature of water-deprived ecosystems. Despite the competition for the same limiting resource, coexistence of several plant species is commonly observed. We propose a two-species reaction-diffusion model based on the single-species Klausmeier model, to analytically investigate the existence of states in which both species coexist. Ecologically, the study finds that coexistence is supported if there is a small difference in the plant species' average fitness, measured by the ratio of a species' capabilities to convert water into new biomass to its mortality rate. Mathematically, coexistence is not a stable solution of the system, but both spatially uniform and patterned coexistence states occur as metastable states. In this context, a metastable solution in which both species coexist corresponds to a long transient (exceeding $10^3$ years in dimensional parameters) to a stable one-species state. This behaviour is characterised by the small size of a positive eigenvalue which has the same order of magnitude as the average fitness difference between the two species. Two mechanisms causing the occurrence of metastable solutions are established: a spatially uniform unstable equilibrium and a stable one-species pattern which is unstable to the introduction of a competitor. We further discuss effects of asymmetric interspecific competition (e.g. shading) on the metastability property.


title: Estimation of the distance between two bodies inside an $n$-dimensional ball of unit volume

abstract: We consider the problem of estimating the distance between two bodies of volume $\varepsilon$ located inside a $n$-dimensional ball $U$ of unit volume for $n\to\infty$. Let $A$ be a closed set with a smooth boundary of the volume $\varepsilon$ ($0 \leq \varepsilon \leq 1/2$) inside a $n$-dimensional ball $U$ of unit volume that implements among all the sets of volume $\varepsilon$ is a set with the smallest possible free surface area, lying in one half-space with respect to a certain hyperplane that passes through the center of the ball. Then $A$ has the same free surface area as the set representing the intersection of a ball perpendicular to $U$ and the ball $U$ itself.


title: Architectural configurations, atlas granularity and functional connectivity with diagnostic value in Autism Spectrum Disorder

abstract: Currently, the diagnosis of Autism Spectrum Disorder (ASD) is dependent upon a subjective, time-consuming evaluation of behavioral tests by an expert clinician. Non-invasive functional MRI (fMRI) characterizes brain connectivity and may be used to inform diagnoses and democratize medicine. However, successful construction of deep learning models from fMRI requires addressing key choices about the model's architecture, including the number of layers and number of neurons per layer. Meanwhile, deriving functional connectivity (FC) features from fMRI requires choosing an atlas with an appropriate level of granularity. Once a model has been built, it is vital to determine which features are predictive of ASD and if similar features are learned across atlas granularity levels. To identify aptly suited architectural configurations, probability distributions of the configurations of high versus low performing models are compared. To determine the effect of atlas granularity, connectivity features are derived from atlases with 3 levels of granularity and important features are ranked with permutation feature importance. Results show the highest performing models use between 2-4 hidden layers and 16-64 neurons per layer, granularity dependent. Connectivity features identified as important across all 3 atlas granularity levels include FC to the supplementary motor gyrus and language association cortex, regions associated with deficits in social and sensory processing in ASD. Importantly, the cerebellum, often not included in functional analyses, is also identified as a region whose abnormal connectivity is highly predictive of ASD. Results of this study identify important regions to include in future studies of ASD, help assist in the selection of network architectures, and help identify appropriate levels of granularity to facilitate the development of accurate diagnostic models of ASD.


title: Women, politics and Twitter: Using machine learning to change the discourse

abstract: Including diverse voices in political decision-making strengthens our democratic institutions. Within the Canadian political system, there is gender inequality across all levels of elected government. Online abuse, such as hateful tweets, leveled at women engaged in politics contributes to this inequity, particularly tweets focusing on their gender. In this paper, we present ParityBOT: a Twitter bot which counters abusive tweets aimed at women in politics by sending supportive tweets about influential female leaders and facts about women in public life. ParityBOT is the first artificial intelligence-based intervention aimed at affecting online discourse for women in politics for the better. The goal of this project is to: $1$) raise awareness of issues relating to gender inequity in politics, and $2$) positively influence public discourse in politics. The main contribution of this paper is a scalable model to classify and respond to hateful tweets with quantitative and qualitative assessments. The ParityBOT abusive classification system was validated on public online harassment datasets. We conclude with analysis of the impact of ParityBOT, drawing from data gathered during interventions in both the $2019$ Alberta provincial and $2019$ Canadian federal elections.


title: ALMA Observations of the Extraordinary Carina Pillars: HH 901/902

abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm continuum and C$^{18}$O(2$-$1), N$_2$D$^{+}$(3$-$2), $^{13}$CS(5$-$4), and $^{12}$CO(2$-$1) line sensitive and high angular resolution ($\sim$0.3$''$) observations of the famous carina pillars and protostellar objects HH 901/902. Our observations reveal for the first time, the bipolar CO outflows and the dusty disks (plus envelopes) that are energizing the extended and irradiated HH objects far from the pillars. We find that the masses of the disks$+$envelopes are about 0.1 M$_\odot$ and of the bipolar outflows are between 10$^{-3}$ - 10$^{-4}$ M$_\odot$, which suggests that they could be low- or maybe intermediate- mass protostars. Moreover, we suggest that these young low-mass stars are likely embedded Class 0/I protostars with high-accretion rates. We also show the kinematics of the gas in the pillars together with their respective gas masses (0.1 -- 0.2 M$_\odot$). We estimate that the pillars will be photo-evaporated in 10$^4$ to 10$^5$ years by the massive and luminous stars located in the Trumpler 14 cluster. Finally, given the short photo-evaporated timescales and that the protostars in these pillars are still very embedded, we suggest that the disks inside of the pillars will be quickly affected by the radiation of the massive stars, forming proplyds, like those observed in Orion.


title: Scalarized Black Hole dynamics in Einstein dilaton Gauss-Bonnet Gravity

abstract: We report on a numerical investigation of the stability of scalarized black holes in Einstein dilaton Gauss-Bonnet (EdGB) gravity in the full dynamical theory, though restricted to spherical symmetry. We find evidence that for sufficiently small curvature-couplings the resulting scalarized black hole solutions are nonlinearly stable. For such small couplings, we show that an elliptic region forms inside these EdGB black hole spacetimes (prior to any curvature singularity), and give evidence that this region remains censored from asymptotic view. However, for coupling values "superextremal" relative to a given black hole mass, an elliptic region forms exterior to the horizon, implying the exterior Cauchy problem is ill-posed in this regime.


title: Deep Decomposition Learning for Inverse Imaging Problems

abstract: Deep learning is emerging as a new paradigm for solving inverse imaging problems. However, the deep learning methods often lack the assurance of traditional physics-based methods due to the lack of physical information considerations in neural network training and deploying. The appropriate supervision and explicit calibration by the information of the physic model can enhance the neural network learning and its practical performance. In this paper, inspired by the geometry that data can be decomposed by two components from the null-space of the forward operator and the range space of its pseudo-inverse, we train neural networks to learn the two components and therefore learn the decomposition, i.e. we explicitly reformulate the neural network layers as learning range-nullspace decomposition functions with reference to the layer inputs, instead of learning unreferenced functions. We empirically show that the proposed framework demonstrates superior performance over recent deep residual learning, unrolled learning and nullspace learning on tasks including compressive sensing medical imaging and natural image super-resolution. Our code is available at https://github.com/edongdongchen/DDN.


title: How complex is the cosmic web?

abstract: The growth of large-scale cosmic structure is a beautiful exemplification of how complexity can emerge in our Universe, starting from simple initial conditions and simple physical laws. Using {\enzo} cosmological numerical simulations, I applied tools from Information Theory (namely, "statistical complexity") to quantify the amount of complexity in the simulated cosmic volume, as a function of cosmic epoch and environment. This analysis can quantify how much difficult to predict, at least in a statistical sense, is the evolution of the thermal, kinetic and magnetic energy of the dominant component of ordinary matter in the Universe (the intragalactic medium plasma). The most complex environment in the simulated cosmic web is generally found to be the periphery of large-scale structures (e.g. galaxy clusters and filaments), where the complexity is on average $\sim 10-10^2$ times larger than in more rarefied regions, even if the latter dominate the volume-integrated complexity of the simulated Universe. If the energy evolution of gas in the cosmic web is measured on a $\approx 100 $ $\rm kpc/h$ resolution and over a $\approx 200$ $\rm Myr$ timescale, its total complexity is the range of $\sim 10^{16}-10^{17} \rm ~bits$, with little dependence on the assumed gas physics, cosmology or cosmic variance.


title: Making Learners (More) Monotone

abstract: Learning performance can show non-monotonic behavior. That is, more data does not necessarily lead to better models, even on average. We propose three algorithms that take a supervised learning model and make it perform more monotone. We prove consistency and monotonicity with high probability, and evaluate the algorithms on scenarios where non-monotone behaviour occurs. Our proposed algorithm $\text{MT}_{\text{HT}}$ makes less than $1\%$ non-monotone decisions on MNIST while staying competitive in terms of error rate compared to several baselines.


title: The $S^3_\bfw$ Sasaki Join Construction

abstract: The main purpose of this work is to generalize the $S^3_\bfw$ Sasaki join construction $M\star_\bfl S^3_\bfw$ described in \cite{BoTo14a} when the Sasakian structure on $M$ is regular, to the general case where the Sasakian structure is only quasi-regular. This gives one of the main results, Theorem 3.2, which describes an inductive procedure for constructing Sasakian metrics of constant scalar curvature. In the Gorenstein case ($c_1(\cald)=0$) we construct a polynomial whose coeffients are linear in the components of $\bfw$ and whose unique root in the interval $(1,\infty)$ completely determines the Sasaki-Einstein metric. In the more general case we apply our results to prove that there exists infinitely many smooth 7-manifolds each of which admit infinitely many inequivalent contact structures of Sasaki type admitting constant scalar curvature Sasaki metrics (see Corollary 6.15). We also discuss the relationship with a recent paper \cite{ApCa18} of Apostolov and Calderbank as well as the relation with K-stability.


title: An optimal regularity result for Kolmogorov equations and weak uniqueness for some critical SPDEs

abstract: We show uniqueness in law for the critical SPDE $$ dX_t = AX_t dt + (-A)^{1/2}F(X(t))dt + dW_t,\;\; X_0 =x \in H, $$ where $A$ $ : dom(A) \subset H \to H$ is a negative definite self-adjoint operator on a separable Hilbert space $H$ having $A^{-1}$ of trace class and $W$ is a cylindrical Wiener process on $H$. Here $F: H \to H $ can be continuous with at most linear growth (some functions $F$ which grow more than linearly can also be considered). This leads to new uniqueness results for generalized stochastic Burgers' equations and for three-dimensional stochastic Cahn-Hilliard type equations which have interesting applications. To get weak uniqueness we also establish a new optimal regularity result for the Kolmogorov equation $ \lambda u - Lu = f$ on $H$, where $\lambda >0$, $ f: H \to {\mathbb R}$ is Borel and bounded and $L$ is the Ornstein-Uhlenbeck operator related to the SPDE when $F=0$. In particular we show that the first derivative $Du : H \to H$ verifies $Du(x) \in \text{dom}((-A)^{1/2})$, for any $x \in H,$ and moreover $$ \sup_{x \in H} |(-A)^{1/2}Du (x)|_H = \| (-A)^{1/2}Du \|_{0} \le C \, \| f\|_{0}. $$


title: Gating Revisited: Deep Multi-layer RNNs That Can Be Trained

abstract: We propose a new STAckable Recurrent cell (STAR) for recurrent neural networks (RNNs), which has fewer parameters than widely used LSTM and GRU while being more robust against vanishing or exploding gradients. Stacking recurrent units into deep architectures suffers from two major limitations: (i) many recurrent cells (e.g., LSTMs) are costly in terms of parameters and computation resources; and (ii) deep RNNs are prone to vanishing or exploding gradients during training. We investigate the training of multi-layer RNNs and examine the magnitude of the gradients as they propagate through the network in the "vertical" direction. We show that, depending on the structure of the basic recurrent unit, the gradients are systematically attenuated or amplified. Based on our analysis we design a new type of gated cell that better preserves gradient magnitude. We validate our design on a large number of sequence modelling tasks and demonstrate that the proposed STAR cell allows to build and train deeper recurrent architectures, ultimately leading to improved performance while being computationally more efficient.


title: Failure Modes in Machine Learning Systems

abstract: In the last two years, more than 200 papers have been written on how machine learning (ML) systems can fail because of adversarial attacks on the algorithms and data; this number balloons if we were to incorporate papers covering non-adversarial failure modes. The spate of papers has made it difficult for ML practitioners, let alone engineers, lawyers, and policymakers, to keep up with the attacks against and defenses of ML systems. However, as these systems become more pervasive, the need to understand how they fail, whether by the hand of an adversary or due to the inherent design of a system, will only become more pressing. In order to equip software developers, security incident responders, lawyers, and policy makers with a common vernacular to talk about this problem, we developed a framework to classify failures into "Intentional failures" where the failure is caused by an active adversary attempting to subvert the system to attain her goals; and "Unintentional failures" where the failure is because an ML system produces an inherently unsafe outcome. After developing the initial version of the taxonomy last year, we worked with security and ML teams across Microsoft, 23 external partners, standards organization, and governments to understand how stakeholders would use our framework. Throughout the paper, we attempt to highlight how machine learning failure modes are meaningfully different from traditional software failures from a technology and policy perspective.


title: Estimating Planetary Mass with Deep Learning

abstract: While thousands of exoplanets have been confirmed, the known properties about individual discoveries remain sparse and depend on detection technique. To utilize more than a small section of the exoplanet dataset, tools need to be developed to estimate missing values based on the known measurements. Here, we demonstrate the use of a neural network that models the density of planets in a space of six properties that is then used to impute a probability distribution for missing values. Our results focus on planetary mass which neither the radial velocity nor transit techniques for planet identification can provide alone. The neural network can impute mass across the four orders of magnitude in the exoplanet archive, and return a distribution of masses for each planet that can inform about trends in the underlying dataset. The average error on this mass estimate from a radial velocity detection is a factor of 1.5 of the observed value, and 2.7 for a transit observation. The mass of Proxima Centauri b found by this method is $1.6^{\rm +0.46}_{\rm -0.36}$ M$_\oplus$, where the upper and lower bounds are derived from the root mean square deviation from the log mass probability distribution. The network can similarly impute the other potentially missing properties, and we use this to predict planet radius for radial velocity measurements, with an average error of a factor 1.4 of the observed value. The ability of neural networks to search for patterns in multidimensional data means that such techniques have the potential to greatly expand the use of the exoplanet catalogue.


title: Upper bounds on the Holevo Cram\'er-Rao bound for multiparameter quantum parametric and semiparametric estimation

abstract: We formulate multiparameter quantum estimation in the parametric and semiparametric setting. While the Holevo Cram\'er-Rao bound (CRB) requires no substantial modifications in moving from the former to the latter, we generalize the Helstrom CRB appropriately. We show that the Holevo CRB cannot be greater than twice the generalized Helstrom CRB. We also present a tighter, intermediate, bound. Finally, we show that for parameters encoded in the first moments of a Gaussian state there always exists a Gaussian measurement that gives a classical Fisher information matrix that is one-half of the quantum Fisher information matrix.


title: Analysis of a model for banded vegetation patterns in semi-arid environments with nonlocal dispersal

abstract: Vegetation patterns are a characteristic feature of semi-arid regions. On hillsides these patterns occur as stripes running parallel to the contours. The Klausmeier model, a coupled reaction-advection-diffusion system, is a deliberately simple model describing the phenomenon. In this paper, we replace the diffusion term describing plant dispersal by a more realistic nonlocal convolution integral to account for the possibility of long-range dispersal of seeds. Our analysis focuses on the rainfall level at which there is a transition between uniform vegetation and pattern formation. We obtain results, valid to leading order in the large parameter comparing the rate of water flow downhill to the rate of plant dispersal, for a negative exponential dispersal kernel. Our results indicate that both a wider dispersal of seeds and an increase in dispersal rate inhibit the formation of patterns. Assuming an evolutionary trade-off between these two quantities, mathematically motivated by the limiting behaviour of the convolution term, allows us to make comparisons to existing results for the original reaction-advection-diffusion system. These comparisons show that the nonlocal model always predicts a larger parameter region supporting pattern formation. We then numerically extend the results to other dispersal kernels, showing that the tendency to form patterns depends on the type of decay of the kernel.


title: Subalgebras generated in degree two with minimal Hilbert function

abstract: What can be said about the subalgebras of the polynomial ring, with minimal or maximal Hilbert function? This question was discussed in a recent paper by M. Boij and A. Conca. In this paper we study the subalgebras generated in degree two with minimal Hilbert function. The problem to determine the generators of these algebras transfers into a combinatorial problem on counting maximal north-east lattice paths inside a shifted Ferrers diagram. We conjecture that the subalgebras generated in degree two with minimal Hilbert function are generated by an initial Lex or RevLex segment.


title: Spectrum of strange singly charmed baryons in the constituent quark model

abstract: Excited states masses of the strange singly charmed baryons are calculated using the non-relativistic approach of hypercentral Constituent Quark Model (hCQM). The hyper-Coulomb plus screened potential is used as a confinement potential with the first order correction. The spin-spin, spin-orbit and spin-tensor interaction terms are included perturbatively. Our calculated masses are allowed to construct the Regge trajectories in both $(n_r, M^2)$ and $(J, M^2)$ planes. The mass spectra and the Regge trajectories study predict the spin-parity of $\Xi{_c(2970)^{+/0}}$, $\Xi{_c(3080)^{+/0}}$, $\Xi{_c(3123)^+}$, $\Omega{_c}(3000)^0$ and $\Omega{_c}(3119)^0$ baryons. Moreover, the strong one pion decay rates of the isodoublet states of $\Xi{_c(2645)}$, $\Xi{_c(2790)}$ and $\Xi{_c(2815)}$ are analyzed in the framework of Heavy Hadron Chiral Perturbation Theory (HHChPT). Also, the ground state magnetic moments and the radiative decay rates based on the transition magnetic moments are calculated in the framework of constituent quark model.


title: Algebras of Non-Local Screenings and Diagonal Nichols Algebras

abstract: In a vertex algebra setting, we consider non-local screening operators associated to the basis of any non-integral lattice. We have previously shown that, under certain restrictions, these screening operators satisfy the relations of a quantum shuffle algebra or Nichols algebra associated to a diagonal braiding, which encodes the non-locality and non-integrality. In the present article, we take all finite-dimensional diagonal Nichols algebras, as classified by Heckenberger, and find all lattice realizations of the braiding that are compatible with reflections. Usually, the realizations are unique or come as one- or two-parameter families. Examples include realizations of Lie superalgebras. We then study the associated algebra of screenings with improved methods. Typically, for positive definite lattices we obtain the Nichols algebra, such as the positive part of the quantum group, and for negative definite lattices we obtain a certain extension of the Nichols algebra generalizing the infinite quantum group with a large center.


title: First hitting times to intermittent targets

abstract: In noisy environments such as the cell, many processes involve target sites that are often hidden or inactive, and thus not always available for reaction with diffusing entities. To understand reaction kinetics in these situations, we study the first hitting time statistics of a Brownian particle searching for a target site that switches stochastically between visible and hidden phases. At high crypticity, an unexpected rate limited power-law regime emerges for the first hitting time density, which markedly differs from the classic $t^{-3/2}$ scaling for steady targets. Our problem admits an asymptotic mapping onto a mixed, or Robin, boundary condition. Similar results are obtained with non-Markov targets and particles diffusing anomalously.


title: Efficient and Stochastic Multireference Perturbation Theory for Large Active Spaces within a Full Configuration Interaction Quantum Monte Carlo Framework

abstract: Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems, and was recently adapted for use as an active space solver and combined with orbital optimisation. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally-contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while CASPT2 seems less amenable to a stochastic reformulation, NEVPT2 is far more stable, requiring a similar number of walkers to converge the NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of NEVPT2 within a (24,24) active space in a biologically relevant system, and show that small numbers of walkers are sufficient for a faithful sampling of the NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.


title: Resampling-based Confidence Intervals for Model-free Robust Inference on Optimal Treatment Regimes

abstract: We propose a new procedure for inference on optimal treatment regimes in the model-free setting, which does not require to specify an outcome regression model. Existing model-free estimators for optimal treatment regimes are usually not suitable for the purpose of inference, because they either have nonstandard asymptotic distributions or do not necessarily guarantee consistent estimation of the parameter indexing the Bayes rule due to the use of surrogate loss. We first study a smoothed robust estimator that directly targets the parameter corresponding to the Bayes decision rule for optimal treatment regimes estimation. This estimator is shown to have an asymptotic normal distribution. Furthermore, we verify that a resampling procedure provides asymptotically accurate inference for both the parameter indexing the optimal treatment regime and the optimal value function. A new algorithm is developed to calculate the proposed estimator with substantially improved speed and stability. Numerical results demonstrate the satisfactory performance of the new methods.


title: Optical non-Hermitian para-Fermi oscillators

abstract: We present a proposal for the optical simulation of para-Fermi oscillators in arrays of coupled waveguides. We use a representation that arises as a deformation of the $su(2)$ algebra. This provides us with a set of chiral and a zero-energy-like normal modes. The latter is its own chiral pair and suggest the addition of controlled losses/gains following a pattern defined by parity. In these non-Hermitian para-Fermi oscillators, the analog of the zero-energy mode presents the largest effective loss/gain and it is possible to tune the system to show sequences of exceptional points and varying effective losses/gains. These arrays can be used for mode suppression or enhancement depending on the use of loss or gain, in that order. We compare our coupled mode theory predictions with finite element method simulations to good agreement.


title: Explicit Pieri Inclusions

abstract: By the Pieri rule, the tensor product of an exterior power and a finite-dimensional irreducible representation of a general linear group has a multiplicity-free decomposition. The embeddings of the constituents are called Pieri inclusions and were first studied by Weyman in his thesis and described explicitly by Olver. More recently, these maps have appeared in the work of Eisenbud, Fl{\o}stad, and Weyman and of Sam and Weyman to compute pure free resolutions for classical groups. In this paper, we give a new closed form, non-recursive description of Pieri inclusions. For partitions with a bounded number of distinct parts, the resulting algorithm has polynomial time complexity whereas the previously known algorithm has exponential time complexity.


title: Higher-Point Conformal Blocks in the Comb Channel

abstract: We compute $M$-point conformal blocks with scalar external and exchange operators in the so-called comb configuration for any $M$ in any dimension $d$. Our computation involves repeated use of the operator product expansion to increase the number of external fields. We check our results in several limits and compare with the expressions available in the literature when $M=5$ for any $d$, and also when $M$ is arbitrary while $d=1$.


title: Quantum differential equations and helices

abstract: These notes partly touch the topic of the talk given by the author at the XXXVIII Workshop on Geometric Methods in Physics, hold in June-July 2019 in Bia\l{}owie\.{z}a, Poland. They consist of a short and self-contained introduction to the isomonodromic approach to quantum cohomology, and Dubrovin's conjecture. An overview of recent results obtained in joint works with B. Dubrovin and D. Guzzetti (arXiv:1811.09235), and A. Varchenko (arXiv:1909.06582) is given.


title: Listing Conflicting Triples in Optimal Time

abstract: Different sources of information might tell different stories about the evolutionary history of a given set of species. This leads to (rooted) phylogenetic trees that "disagree" on triples of species, which we call "conflict triples". An important subtask of computing consensus trees which is interesting in its own regard is the enumeration of all conflicts exhibited by a pair of phylogenetic trees (on the same set of $n$ taxa). As it is possible that a significant part of the $n^3$ triples are in conflict, the trivial ${\Theta}(n^3)$-time algorithm that checks for each triple whether it constitutes a conflict, was considered optimal. It turns out, however, that we can do way better in the case that there are only few conflicts. In particular, we show that we can enumerate all d conflict triples between a pair of phylogenetic trees in $O(n + d)$ time. Since any deterministic algorithm has to spend ${\Theta}(n)$ time reading the input and ${\Theta}(d)$ time writing the output, no deterministic algorithm can solve this task faster than we do (up to constant factors).