title: Towers of solutions of qKZ equations and their applications to loop models

abstract: Cherednik's type A quantum affine Knizhnik-Zamolodchikov (qKZ) equations form a consistent system of linear $q$-difference equations for $V_n$-valued meromorphic functions on a complex $n$-torus, with $V_n$ a module over the GL${}_n$-type extended affine Hecke algebra $\mathcal{H}_n$. The family $(\mathcal{H}_n)_{n\geq 0}$ of extended affine Hecke algebras forms a tower of algebras, with the associated algebra morphisms $\mathcal{H}_n\rightarrow\mathcal{H}_{n+1}$ the Hecke algebra descends of arc insertion at the affine braid group level. In this paper we consider qKZ towers $(f^{(n)})_{n\geq 0}$ of solutions, which consist of twisted-symmetric polynomial solutions $f^{(n)}$ ($n\geq 0$) of the qKZ equations that are compatible with the tower structure on $(\mathcal{H}_n)_{n\geq 0}$. The compatibility is encoded by so-called braid recursion relations: $f^{(n+1)}(z_1,\ldots,z_{n},0)$ is required to coincide up to a quasi-constant factor with the push-forward of $f^{(n)}(z_1,\ldots,z_{n})$ by an intertwiner $\mu_{n}: V_{n}\rightarrow V_{n+1}$ of $\mathcal{H}_{n}$-modules, where $V_{n+1}$ is considered as an $\mathcal{H}_{n}$-module through the tower structure on $(\mathcal{H}_n)_{n\geq 0}$. We associate to the dense loop model on the half-infinite cylinder with nonzero loop weights a qKZ tower $(f^{(n)})_{n\geq 0}$ of solutions. The solutions $f^{(n)}$ are constructed from specialised dual non-symmetric Macdonald polynomials with specialised parameters using the Cherednik-Matsuo correspondence. In the special case that the extended affine Hecke algebra parameter is a third root of unity, $f^{(n)}$ coincides with the (suitably normalized) ground state of the inhomogeneous dense $O(1)$ loop model on the half-infinite cylinder with circumference $n$.


title: Probing the Need for Visual Context in Multimodal Machine Translation

abstract: Current work on multimodal machine translation (MMT) has suggested that the visual modality is either unnecessary or only marginally beneficial. We posit that this is a consequence of the very simple, short and repetitive sentences used in the only available dataset for the task (Multi30K), rendering the source text sufficient as context. In the general case, however, we believe that it is possible to combine visual and textual information in order to ground translations. In this paper we probe the contribution of the visual modality to state-of-the-art MMT models by conducting a systematic analysis where we partially deprive the models from source-side textual context. Our results show that under limited textual context, models are capable of leveraging the visual input to generate better translations. This contradicts the current belief that MMT models disregard the visual modality because of either the quality of the image features or the way they are integrated into the model.


title: $\kappa$-deformed phase spaces, Jordanian twists, Lorentz-Weyl algebra and dispersion relations

abstract: We consider $\kappa$-deformed relativistic quantum phase space and possible implementations of the Lorentz algebra. There are two ways of performing such implementations. One is a simple extension where the Poincar\'e algebra is unaltered, while the other is a general extension where the Poincar\'e algebra is deformed. As an example we fix the Jordanian twist and the corresponding realization of noncommutative coordinates, coproduct of momenta and addition of momenta. An extension with a one-parameter family of realizations of the Lorentz generators, dilatation and momenta closing the Poincar\'e-Weyl algebra is considered. The corresponding physical interpretation depends on the way the Lorentz algebra is implemented in phase space. We show how the spectrum of the relativistic hydrogen atom depends on the realization of the generators of the Poincar\'e-Weyl algebra.


title: A 68 uW 31 kS/s Fully-Capacitive Noise-Shaping SAR ADC with 102 dB SNDR

abstract: This paper presents a 17 bit analogue-to-digital converter that incorporates mismatch and quantisation noise-shaping techniques into an energy-saving 10 bit successive approximation quantiser to increase the dynamic range by another 42 dB. We propose a novel fully-capacitive topology which allows for high-speed asynchronous conversion together with a background calibration scheme to reduce the oversampling requirement by 10x compared to prior-art. A 0.18 um CMOS technology is used to demonstrate preliminary simulation results together with analytic measures that optimise parameter and topology selection. The proposed system is able to achieve a FoMS of 183 dB for a maximum signal bandwidth of 15.6 kHz while dissipating 68 uW from a 1.8 V supply. A peak SNDR of 102 dB is demonstrated for this rate with a 0.201 mm^2 area requirement.


title: Towards observing reverberating and superefficient pulsar wind nebulae

abstract: In a recent work, we numerically studied the radiative properties of the reverberation phase of pulsar wind nebulae (PWNe), i.e., when the reverse shock created by the supernova explosion travels back towards the pulsar, compressing the wind bubble. We focused on several well-characterized PWNe and used them as examples for introducing the concept of superefficiency. The latter is a period of the PWN evolution, happening within reverberation, where the luminosity in a given band exceeds the spin-down power at the time. Here, we explore a broad range of PWN models to study their reverberation and superefficiency phases in a systematic way. Armed with these models we consider two aspects: On the one hand, we analyze via Monte Carlo simulations how many Galactic PWNe are expected to be reverberating or in a superefficiency stage at any given time, providing the first such estimations. On the other hand, we focus on searching for observational signatures of such periods. We analyze archival observations and check for the existence of possible candidates for superefficient PWNe. We also provide predictions for the future evolution of the magnetar nebula J1834.9-0846 (which we consider to be starting its reverberation period) along the next 50 years. Using our simulations as input we study how sensitive current and future X-ray satellites (like eXTP or Athena) will be to observe such evolution, concluding that they will be able to track it in detail.


title: Im2Pencil: Controllable Pencil Illustration from Photographs

abstract: We propose a high-quality photo-to-pencil translation method with fine-grained control over the drawing style. This is a challenging task due to multiple stroke types (e.g., outline and shading), structural complexity of pencil shading (e.g., hatching), and the lack of aligned training data pairs. To address these challenges, we develop a two-branch model that learns separate filters for generating sketchy outlines and tonal shading from a collection of pencil drawings. We create training data pairs by extracting clean outlines and tonal illustrations from original pencil drawings using image filtering techniques, and we manually label the drawing styles. In addition, our model creates different pencil styles (e.g., line sketchiness and shading style) in a user-controllable manner. Experimental results on different types of pencil drawings show that the proposed algorithm performs favorably against existing methods in terms of quality, diversity and user evaluations.


title: Approximation of Fractional Local Times: Zero Energy and Derivatives

abstract: We consider empirical processes associated with high-frequency observations of a fractional Brownian motion (fBm) $X$ with Hurst parameter $H\in (0,1)$, and derive conditions under which these processes verify a (possibly uniform) law of large numbers, as well as a second order (possibly uniform) limit theorem. We devote specific emphasis to the `zero energy' case, corresponding to a kernel whose integral on the real line equals zero. Our asymptotic results are associated with explicit rates of convergence, and are expressed either in terms of the local time of $X$ or of its \blue{derivatives}: in particular, the full force of our finding applies to the `rough range' $0< H < 1/3$, on which the previous literature has been mostly silent. The {\color{black}use of the derivatives} of local times for studying the fluctuations of high-frequency observations of a fBm is new, and is the main technological breakthrough of the present paper. Our results are based on the use of Malliavin calculus and Fourier analysis, and extend and complete several findings in the literature, e.g. by Jeganathan (2004, 2006, 2008) and Podolskij and Rosenbaum (2018).


title: Addressing Temporal Variations in Qubit Quality Metrics for Parameterized Quantum Circuits

abstract: The public access to noisy intermediate-scale quantum (NISQ) computers facilitated by IBM, Rigetti, D-Wave, etc., has propelled the development of quantum applications that may offer quantum supremacy in the future large-scale quantum computers. Parameterized quantum circuits (PQC) have emerged as a major driver for the development of quantum routines that potentially improve the circuit's resilience to the noise. PQC's have been applied in both generative (e.g. generative adversarial network) and discriminative (e.g. quantum classifier) tasks in the field of quantum machine learning. PQC's have been also considered to realize high fidelity quantum gates with the available imperfect native gates of a target quantum hardware. Parameters of a PQC are determined through an iterative training process for a target noisy quantum hardware. However, temporal variations in qubit quality metrics affect the performance of a PQC. Therefore, the circuit that is trained without considering temporal variations exhibits poor fidelity over time. In this paper, we present training methodologies for PQC in a completely classical environment that can improve the fidelity of the trained PQC on a target NISQ hardware by as much as 42.5%.


title: Theory of correlated insulating behaviour and spin-triplet superconductivity in twisted double bilayer graphene

abstract: Two monolayers of graphene twisted by a small `magic' angle exhibit nearly flat bands leading to correlated electronic states and superconductivity, whose precise nature including possible broken symmetries, remain under debate. Here we theoretically study a related but different system with reduced symmetry - twisted {\em double} bilayer graphene (TDBLG), consisting of {\em two} Bernal stacked bilayer graphene sheets, twisted with respect to one another. Unlike the monolayer case, we show that isolated flat bands only appear on application of a vertical displacement field $D$. We construct a phase diagram as a function of twist angle and $D$, incorporating interactions via a Hartree-Fock approximation. At half filling, ferromagnetic insulators are stabilized, typically with valley Chern number $C_v=2$. Ferromagnetic fluctuations in the metallic state are argued to lead to spin triplet superconductivity from pairing between electrons in opposite valleys. Response of these states to a magnetic field applied either perpendicular or parallel to the graphene sheets is obtained, and found to compare favorably with a recent experiment. We highlight a novel orbital effect arising from in-plane fields that can exceed the Zeeman effect and plays an important role in interpreting experiments.


title: Quantum variance on quaternion algebras, III

abstract: We determine the asymptotic quantum variance of microlocal lifts of Hecke--Maass cusp forms on the arithmetic compact hyperbolic surfaces attached to maximal orders in quaternion algebras. Our result extends those of Luo--Sarnak--Zhao concerning the non-compact modular surface. The results of this article's prequel (which involved the theta correspondence, Rallis inner product formula and equidistribution of translates of elementary theta functions) reduce the present task to some local problems over the reals involving the construction and analysis of microlocal lifts via integral operators on the group. We address these here using an analytic incarnation of the method of coadjoint orbits.


title: On approximate validation of models: A Kolmogorov-Smirnov based approach

abstract: Classical tests of fit typically reject a model for large enough real data samples. In contrast, often in statistical practice a model offers a good description of the data even though it is not the "true" random generator. We consider a more flexible approach based on contamination neighbourhoods around a model. Using trimming methods and the Kolmogorov metric we introduce a functional statistic measuring departures from a contaminated model and the associated estimator corresponding to its sample version. We show how this estimator allows testing of fit for the (slightly) contaminated model vs sensible deviations from it, with uniformly exponentially small type I and type II error probabilities. We also address the asymptotic behavior of the estimator showing that, under suitable regularity conditions, it asymptotically behaves as the supremum of a Gaussian process. As an application we explore methods of comparison between descriptive models based on the paradigm of model falseness. We also include some connections of our approach with the False-Discovery-Rate setting, showing competitive behavior when estimating the contamination level, although applicable in a wider framework.


title: Stochastic Gradient Descent with Polyak's Learning Rate

abstract: Stochastic gradient descent (SGD) for strongly convex functions converges at the rate $\bO(1/k)$. However, achieving good results in practice requires tuning the parameters (for example the learning rate) of the algorithm. In this paper we propose a generalization of the Polyak step size, used for subgradient methods, to Stochastic gradient descent. We prove a non-asymptotic convergence at the rate $\bO(1/k)$ with a rate constant which can be better than the corresponding rate constant for optimally scheduled SGD. We demonstrate that the method is effective in practice, and on convex optimization problems and on training deep neural networks, and compare to the theoretical rate.


title: Implicit Generation and Generalization in Energy-Based Models

abstract: Energy based models (EBMs) are appealing due to their generality and simplicity in likelihood modeling, but have been traditionally difficult to train. We present techniques to scale MCMC based EBM training on continuous neural networks, and we show its success on the high-dimensional data domains of ImageNet32x32, ImageNet128x128, CIFAR-10, and robotic hand trajectories, achieving better samples than other likelihood models and nearing the performance of contemporary GAN approaches, while covering all modes of the data. We highlight some unique capabilities of implicit generation such as compositionality and corrupt image reconstruction and inpainting. Finally, we show that EBMs are useful models across a wide variety of tasks, achieving state-of-the-art out-of-distribution classification, adversarially robust classification, state-of-the-art continual online class learning, and coherent long term predicted trajectory rollouts.


title: Efficient Inner Product Approximation in Hybrid Spaces

abstract: Many emerging use cases of data mining and machine learning operate on large datasets with data from heterogeneous sources, specifically with both sparse and dense components. For example, dense deep neural network embedding vectors are often used in conjunction with sparse textual features to provide high dimensional hybrid representation of documents. Efficient search in such hybrid spaces is very challenging as the techniques that perform well for sparse vectors have little overlap with those that work well for dense vectors. Popular techniques like Locality Sensitive Hashing (LSH) and its data-dependent variants also do not give good accuracy in high dimensional hybrid spaces. Even though hybrid scenarios are becoming more prevalent, currently there exist no efficient techniques in literature that are both fast and accurate. In this paper, we propose a technique that approximates the inner product computation in hybrid vectors, leading to substantial speedup in search while maintaining high accuracy. We also propose efficient data structures that exploit modern computer architectures, resulting in orders of magnitude faster search than the existing baselines. The performance of the proposed method is demonstrated on several datasets including a very large scale industrial dataset containing one billion vectors in a billion dimensional space, achieving over 10x speedup and higher accuracy against competitive baselines.


title: Tunable MEMS VCSEL on Silicon substrate

abstract: We present design, fabrication and characterization of a MEMS VCSEL which utilizes a silicon-on-insulator wafer for the microelectromechanical system and encapsulates the MEMS by direct InP wafer bonding, which improves the protection and control of the tuning element. This procedure enables a more robust fabrication, a larger free spectral range and facilitates bidirectional tuning of the MEMS element. The MEMS VCSEL device uses a high contrast grating mirror on a MEMS stage as the bottom mirror, a wafer bonded InP with quantum wells for amplification and a deposited dielectric DBR as the top mirror. A 40 nm tuning range and a mechanical resonance frequency in excess of 2 MHz are demonstrated.


title: The parallelism between galaxy clusters and early-type galaxies: I. The light and mass profiles

abstract: We have analyzed the parallelism between the properties of galaxy clusters and early-type galaxies (ETGs) by looking at the similarity between their light profiles. We find that the equivalent luminosity profiles of all these systems in the \vfilt\ band, once normalized to the effective radius \re\ and shifted in surface brightness, can be fitted by the S\'ersic's law \Sers\ and superposed with a small scatter ($\le0.3$ mag). By grouping objects in different classes of luminosity, the average profile of each class slightly deviates from the other only in the inner and outer regions (outside $0.1\leq r/R_e\leq 3$), but the range of values of $n$ remains ample for the members of each class, indicating that objects with similar luminosity have quite different shapes. The "Illustris" simulation reproduces quite well the luminosity profiles of ETGs, with the exception of in the inner and outer regions where feedback from supernovae and active galactic nuclei, wet and dry mergers, are at work. The total mass and luminosity of galaxy clusters as well as their light profiles are not well reproduced. By exploiting simulations we have followed the variation of the effective half-light and half-mass radius of ETGs up to $z=0.8$, noting that progenitors are not necessarily smaller in size than current objects. We have also analyzed the projected dark+baryonic and dark-only mass profiles discovering that after a normalization to the half-mass radius, they can be well superposed and fitted by the S\'ersic's law.


title: Using Local Experiences for Global Motion Planning

abstract: Sampling-based planners are effective in many real-world applications such as robotics manipulation, navigation, and even protein modeling. However, it is often challenging to generate a collision-free path in environments where key areas are hard to sample. In the absence of any prior information, sampling-based planners are forced to explore uniformly or heuristically, which can lead to degraded performance. One way to improve performance is to use prior knowledge of environments to adapt the sampling strategy to the problem at hand. In this work, we decompose the workspace into local primitives, memorizing local experiences by these primitives in the form of local samplers, and store them in a database. We synthesize an efficient global sampler by retrieving local experiences relevant to the given situation. Our method transfers knowledge effectively between diverse environments that share local primitives and speeds up the performance dramatically. Our results show, in terms of solution time, an improvement of multiple orders of magnitude in two traditionally challenging high-dimensional problems compared to state-of-the-art approaches.


title: Spectral singularities of a potential created by two coupled microring resonators

abstract: Two microring resonators, one with gain and one with loss, coupled to each other and to a bus waveguide, create an effective non-Hermitian potential for light propagating in the waveguide. Due to geometry, coupling for each microring resonator yields two counter-propagating modes with equal frequencies. We show that such a system enables implementation of many types of scattering peculiarities, which are either the second or fourth order. The spectral singularities separate parameter regions where the spectrum is either pure real or else comprises complex eigenvalues; hence, they represent the points of the phase transition. By modifying the gain-loss relation for the resonators such an optical scatterer can act as a laser, as a coherent perfect absorber, be unidirectionally reflectionless or transparent, and support bound states either growing or decaying in time. These characteristics are observed for a discrete series of the incident-radiation wavelengths.


title: Plane-wave scattering by an ellipsoid composed of an orthorhombic dielectric-magnetic material with arbitrarily oriented constitutive principal axes

abstract: The extended boundary condition method can be formulated to study plane-wave scattering by an ellipsoid composed of an orthorhombic dielectric-magnetic material whose relative permittivity dyadic is a scalar multiple of its relative permeability dyadic, when the constitutive principal axes are arbitrarily oriented with respect to the shape principal axes. Known vector spherical wavefunctions are used to represent the fields in the surrounding matter-free space. After deriving closed-form expressions for the vector spherical wavefunctions for the scattering material, the internal fields are represented as superpositions of those vector spherical wavefunctions. The unknown scattered-field coefficients are related to the known incident-field coefficients by a transition matrix. The total scattering and absorption efficiencies are highly affected by the orientation of the constitutive principal axes relative to the shape principal axes, and the effect of the orientational mismatch between the two sets of principal axes is more pronounced as the electrical size increases. The dependence of the total scattering efficiency, but not of the absorption efficiency, on the angle of rotation about a shape principal axis can be predicted qualitatively from the variation of a scalar function with respect to the angle of rotation. The total scattering and absorption efficiencies do not depend on the polarization state of incident plane wave when the scattering material is impedance-matched to free space. The polarization state of the incident plane wave has a more discernible effect on the total scattering and absorption efficiencies for ellipsoids compared to spheres.


title: MUSE observations of a changing-look AGN I: The re-appearance of the broad emission lines

abstract: Optical changing-look Active Galactic Nuclei (AGN) are a class of sources that change type within a short timescale of years or decades. This change is characterised by the appearance or disappearance of broad emission lines, often associated with dramatic AGN continuum flux changes that are orders of magnitude larger than those expected from typical AGN variability. In this work we study for the first time the host galaxy of a changing-look AGN, Mrk 590, using high spatial resolution optical and near-infrared observations. We discover that after ~ 10 yr absence, the optical broad emission lines of Mrk 590 have reappeared. The AGN optical continuum flux however, is still ~ 10 times lower than that observed during the most luminous state in the 1990s. The host galaxy shows a 4.5 kpc radius star-forming ring with knots of ionised and cold molecular gas emission. Extended ionised and warm molecular gas emission are detected in the nucleus, indicating that there is a reservoir of gas as close as 60 pc from the black hole. We observe a nuclear gas spiral between radii r ~ 0.5 - 2 kpc, which has been suggested as a dynamical mechanism able to drive the necessary gas to fuel AGN. We also discover blue-shifted and high velocity dispersion [O III] emission out to a radius of 1 kpc, tracing a nuclear gas outflow. The gas dynamics in Mrk 590 suggest a complex balance between gas inflow and outflow in the nucleus of the galaxy.


title: Trace decategorification of tensor product algebras

abstract: We show that in ADE type the trace of Webster's categorification of a tensor product of irreducibles for the quantum group is isomorphic to a tensor product of Weyl modules for the current algebra $\dot{U}(\mathfrak{g}[t])$. This extends a result of Beliakova, Habiro, Lauda, and Webster who showed that the trace of the categorified quantum group $\dot{\mathcal{U}}^*(\mathfrak{g})$ is isomorphic to $\dot{U}(\mathfrak{g}[t])$, and the trace of a cyclotomic quotient of $\dot{\mathcal{U}}^*(\mathfrak{g})$, which categorifies a single irreducible for the quantum group, is isomorphic to a Weyl module for $\dot{U}(\mathfrak{g}[t])$. We use a deformation argument based on Webster's technique of unfurling 2-representations.


title: Electric-field induced helium-helium resonances

abstract: The tunability of the helium-helium interaction through an external electric field is investigated. For a static external field, electric-field induced resonances and associated electric-field induced bound states are calculated for the $^4$He-$^4$He, $^3$He-$^4$He, and $^3$He-$^3$He systems. Qualitative agreement is found with the literature for the $^3$He-$^4$He and $^3$He-$^3$He systems [E. Nielsen, D. V. Fedorov, and A. S. Jensen, Phys. Rev. Lett. {\bf{82}}, 2844 (1999)]. The implications of the predicted electric-field induced resonances for $^4$He-$^4$He on the wave packet dynamics, initiated by intense laser pulses, are investigated. Our results are expected to guide next generation experiments.


title: Realization of a quantum autoencoder for lossless compression of quantum data

abstract: As a ubiquitous aspect of modern information technology, data compression has a wide range of applications. Therefore, a quantum autoencoder which can compress quantum information into a low-dimensional space is fundamentally important to achieve automatic data compression in the field of quantum information. Such a quantum autoencoder can be implemented through training the parameters of a quantum device using classical optimization algorithms. In this article, we analyze the condition of achieving a perfect quantum autoencoder and theoretically prove that a quantum autoencoder can losslessly compress high-dimensional quantum information into a low-dimensional space (also called latent space) if the number of maximum linearly independent vectors from input states is no more than the dimension of the latent space. Also, we experimentally realize a universal two-qubit unitary gate and design a quantum autoencoder device by applying machine learning method. Experimental results demonstrate that our quantum autoencoder is able to compress two two-qubit states into two one-qubit states. Besides compressing quantum information, the quantum autoencoder is used to experimentally discriminate two groups of nonorthogonal states.


title: Quantum non-Markovian collision models from colored-noise baths

abstract: A quantum collision model (CM), also known as repeated interactions model, can be built from the standard microscopic framework where a system S is coupled to a white-noise bosonic bath under the rotating wave approximation, which typically results in Markovian dynamics. Here, we discuss how to generalize the CM construction to the case of frequency-dependent system-bath coupling, which defines a class of colored-noise baths. This leads to an intrinsically non-Markovian CM, where each ancilla (bath subunit) collides repeatedly with S at different steps. We discuss the illustrative example of an atom in front a mirror in the regime of non-negligible retardation times.


title: LaserNet: An Efficient Probabilistic 3D Object Detector for Autonomous Driving

abstract: In this paper, we present LaserNet, a computationally efficient method for 3D object detection from LiDAR data for autonomous driving. The efficiency results from processing LiDAR data in the native range view of the sensor, where the input data is naturally compact. Operating in the range view involves well known challenges for learning, including occlusion and scale variation, but it also provides contextual information based on how the sensor data was captured. Our approach uses a fully convolutional network to predict a multimodal distribution over 3D boxes for each point and then it efficiently fuses these distributions to generate a prediction for each object. Experiments show that modeling each detection as a distribution rather than a single deterministic box leads to better overall detection performance. Benchmark results show that this approach has significantly lower runtime than other recent detectors and that it achieves state-of-the-art performance when compared on a large dataset that has enough data to overcome the challenges of training on the range view.


title: Resummed photon spectrum from dark matter annihilation for intermediate and narrow energy resolution

abstract: The annihilation cross section of weakly interacting TeV scale dark matter particles $\chi^0$ into photons is affected by large quantum corrections due to electroweak Sudakov logarithms and the Sommerfeld effect. We extend our previous work on the resummation of the semi-inclusive photon energy spectrum in $\chi^0\chi^0\to \gamma+X$ in the vicinity of the maximal photon energy $E_\gamma = m_\chi$ with NLL' accuracy from the case of narrow photon energy resolution $E^\gamma_{\rm res}$ of order $m_W^2/m_\chi$ to intermediate resolution of order $E^\gamma_{\rm res} \sim m_W$. We also provide details on the previous narrow resolution calculation. The two calculations, performed in different effective field theory set-ups for the wino dark matter model, are then shown to match well, providing an accurate representation up to energy resolutions of about 300 GeV.


title: On Stable and Unstable Behaviours of Certain Rotation Segments

abstract: In this paper, we study non-wandering homeomorphisms of the two torus in the identity homotopy class, whose rotation sets are non-trivial line segments from $(0,0)$ to some totally irrational vector $(\alpha,\beta)$. We show this rotation set is in fact a non-generic phenomenon for any $C^r$ diffeomorphisms, with $r \geq 1$. When such a rotation set does happen, assuming several natural conditions that are generically satisfied in the area-preserving world, we give a clearer description of its rotational behavior. More precisely, the dynamics admits bounded deviation along the direction $-(\alpha,\beta)$ in the lift, and the rotation set is locked inside an arbitrarily small cone with respect to small $C^0$-perturbations of the dynamics. On the other hand, for any non-wandering homeomorphism $f$ with this kind of rotation set, we also present a perturbation scheme in order for the rotation set to be eaten by rotation sets of nearby dynamics, in the sense that the later set has non-empty interior and contains the former one. These two flavors interplay and share the common goal of understanding the stability/instability properties of this kind of rotation set.


title: Omitted variable bias of Lasso-based inference methods: A finite sample analysis

abstract: We study the finite sample behavior of Lasso-based inference methods such as post double Lasso and debiased Lasso. We show that these methods can exhibit substantial omitted variable biases (OVBs) due to Lasso not selecting relevant controls. This phenomenon can occur even when the coefficients are sparse and the sample size is large and larger than the number of controls. Therefore, relying on the existing asymptotic inference theory can be problematic in empirical applications. We compare the Lasso-based inference methods to modern high-dimensional OLS-based methods and provide practical guidance.


title: Spectrophotometric redshifts for $\mathrm{z\sim1}$ galaxies and predictions for number densities with WFIRST and Euclid

abstract: We investigate the accuracy of 4000\AA/Balmer-break based redshifts by combining Hubble Space Telescope ({\it HST}) grism data with photometry. The grism spectra are from the Probing Evolution And Reionization Spectroscopically (PEARS) survey with {\it HST} using the G800L grism on the Advanced Camera for Surveys (ACS). The photometric data come from a compilation by the 3D-HST collaboration of imaging from multiple surveys (notably CANDELS and 3D-HST). We show evidence that spectrophotometric redshifts (SPZs) typically improve the accuracy of photometric redshifts by $\sim$17--60\%. Our SPZ method is a template fitting based routine which accounts for correlated data between neighboring points within grism spectra via the covariance matrix formalism and also accounts for galaxy morphology along the dispersion direction. We show that the robustness of the SPZ is directly related to the fidelity of the D4000 measurement. We also estimate the accuracy of continuum-based redshifts, i.e., for galaxies that do not contain strong emission lines, based on the grism data alone ($\sigma^{\rm NMAD}_{\Delta z/(1+z)}{\lesssim}0.06$). Given that future space-based observatories like WFIRST and Euclid will spend a significant fraction of time on slitless spectroscopic observations, we estimate number densities for objects with $\ | \mathrm{\Delta z/(1+z_{s})} \ | \leq 0.02$. We predict $\sim$700--4400 galaxies/degree$^2$ for galaxies with D4000$>$1.1 and $\ | \mathrm{\Delta z/(1+z_{s})} \ | \leq 0.02$ to a limiting depth of $i_{AB}$=24 mag. This is \emph{especially} important in the absence of an accompanying rich photometric dataset like the existing one for the CANDELS fields, where redshift accuracy from future surveys will rely only on the presence of a feature like the 4000\AA/Balmer breaks or the presence of emission lines within the grism spectra.


title: Strong approximation of stochastic differential equations driven by a time-changed Brownian motion with time-space-dependent coefficients

abstract: The rate of strong convergence is investigated for an approximation scheme for a class of stochastic differential equations driven by a time-changed Brownian motion, where the random time changes $(E_t)_{t\ge 0}$ considered include the inverses of stable and tempered stable subordinators as well as their mixtures. Unlike those in the work of Jum and Kobayashi (2016), the coefficients of the stochastic differential equations discussed in this paper depend on the regular time variable $t$ rather than the time change $E_t$. This alteration makes it difficult to apply the method used in that paper. To overcome this difficulty, we utilize a Gronwall-type inequality involving a stochastic driver to control the moment of the error process. Moreover, in order to guarantee that an ultimately derived error bound is finite, we establish a useful criterion for the existence of exponential moments of powers of the random time change.


title: Planewave scattering by an ellipsoid composed of an orthorhombic dielectric magnetic material

abstract: The extended boundary condition method can be used to study planewave scattering by an ellipsoid composed of an orthorhombic dielectric-magnetic material whose relative permittivity dyadic is a scalar multiple of its relative permeability dyadic. The scattered and internal field phasors can be expanded in terms of appropriate vector spherical wavefunctions with unknown expansion coefficients, whereas the incident field phasors can be similarly expanded but with known expansion coefficients. The scattered-field coefficients are related to the incident-field coefficients through a matrix. The scattering, absorption, and extinction efficiencies were calculated thereby in relation to the propagation direction and the polarization state of the incident plane wave, the constitutive-anisotropy parameters, and the nonsphericity parameters of the ellipsoid, when the eigenvectors of the real permittivity dyadic are aligned along the three semi-axes of the ellipsoid. As the electrical size of the ellipsoid increases, multiple lobes appear in the scattering pattern. The total scattering efficiency can be smaller than the absorption efficiency for some configurations of the incident plane wave but not necessarily for others. The nonsphericity of the object has a stronger influence on the total scattering efficiency than on the absorption efficiency. The forward-scattering efficiency increases monotonically with the electrical size for all configurations of the incident plane wave, and so does the backscattering efficiency for some configurations. For other configurations, the backscattering efficiency has an undulating behavior with increase in electrical size, and is highly affected by the shape and the constitutive anisotropy of the ellipsoid.


title: Accelerating Gradient Boosting Machine

abstract: Gradient Boosting Machine (GBM) is an extremely powerful supervised learning algorithm that is widely used in practice. GBM routinely features as a leading algorithm in machine learning competitions such as Kaggle and the KDDCup. In this work, we propose Accelerated Gradient Boosting Machine (AGBM) by incorporating Nesterov's acceleration techniques into the design of GBM. The difficulty in accelerating GBM lies in the fact that weak (inexact) learners are commonly used, and therefore the errors can accumulate in the momentum term. To overcome it, we design a "corrected pseudo residual" and fit best weak learner to this corrected pseudo residual, in order to perform the z-update. Thus, we are able to derive novel computational guarantees for AGBM. This is the first GBM type of algorithm with theoretically-justified accelerated convergence rate. Finally we demonstrate with a number of numerical experiments the effectiveness of AGBM over conventional GBM in obtaining a model with good training and/or testing data fidelity.


title: Stable loops and almost transverse surfaces

abstract: We show that the cone over a fibered face of a compact fibered hyperbolic 3-manifold is dual to the cone generated by the homology classes of finitely many curves called minimal stable loops living in the associated veering triangulation. We also present a new, more hands-on proof of Mosher's Transverse Surface Theorem.


title: Simplex space-time meshes in thermally coupled two-phase flow simulations of mold filling

abstract: The quality of plastic parts produced through injection molding depends on many factors. Especially during the filling stage, defects such as weld lines, burrs, or insufficient filling can occur. Numerical methods need to be employed to improve product quality by means of predicting and simulating the injection molding process. In the current work, a highly viscous incompressible non-isothermal two-phase flow is simulated, which takes place during the cavity filling. The injected melt exhibits a shear-thinning behavior, which is described by the Carreau-WLF model. Besides that, a novel discretization method is used in the context of 4D simplex space-time grids [2]. This method allows for local temporal refinement in the vicinity of, e.g., the evolving front of the melt [10]. Utilizing such an adaptive refinement can lead to locally improved numerical accuracy while maintaining the highest possible computational efficiency in the remaining of the domain. For demonstration purposes, a set of 2D and 3D benchmark cases, that involve the filling of various cavities with a distributor, are presented.


title: Anticanonical system of Fano fivefolds

abstract: We show that any Fano fivefold with canonical Gorenstein singularities has an effective anticanonical divisor. Moreover,if a general element of the anticanonical system is reduced, then it has canonical singularities. We also prove nonvanishing of anticanonical system in the case of log canonical Gorenstein singularities.


title: A Review of Sensing and Communication, Human Factors, and Controller Aspects for Information-Aware Connected and Automated Vehicles

abstract: Information-aware connected and automated vehicles (CAVs) have drawn great attention in recent years due to its potentially significant positive impacts on roadway safety and operational efficiency. In this paper, we conduct an in-depth review of three basic and key interrelated aspects of a CAV: sensing and communication technologies, human factors, and information-aware controller design. First, different vehicular sensing and communication technologies and their protocol stacks, to provide reliable information to the information-aware CAV controller, are thoroughly discussed. Diverse human factor issues, such as user comfort, preferences, and reliability, to design the CAV systems for mass adaptation are also discussed. Then, different layers of a CAV controller (route planning, driving mode execution, and driving model selection) considering human factors and information through connectivity are reviewed. In addition, critical challenges for the sensing and communication technologies, human factors, and information-aware controller are identified to support the design of a safe and efficient CAV system while considering user acceptance and comfort. Finally, promising future research directions of these three aspects are discussed to overcome existing challenges to realize a safe and operationally efficient CAV.


title: Multi-band non-equilibrium GW+EDMFT formalism for correlated insulators

abstract: We study the dynamics of charge-transfer insulators after a photo-excitation using the three-band Emery model which is relevant for the description of cuprate superconductors. We provide a detailed derivation of the nonequilibrium extension of the multi-band GW+EDMFT formalism and the corresponding downfolding procedure. The Peierls construction of the electron-light coupling is generalized to the multi-band case resulting in a gauge invariant combination of the Peierls intra-band acceleration and dipolar intra-band transitions. We apply the formalism to the study of momentum-dependent (inverse) photo-emission spectra and optical conductivities. The time-resolved spectral function shows a strong renormalization of the charge-transfer gap and a substantial broadening of some of the bands. While the upper Hubbard band exhibits a momentum-dependent broadening, an almost rigid band shift is observed for the ligand bands. The inverse photo-emission spectrum reveals that the inclusion of the non-local and inter-band charge fluctuations lead to a very fast relaxation of holes into the lower Hubbard band. Consistent with the changes in the spectral function, the optical conductivity shows a renormalization of the charge-transfer gap, which is proportional to the photo-doping. The details of the photo-induced changes strongly depend on the dipolar matrix elements, which calls for an ab-initio determination of these parameters.


title: The Astrophysical Multimessenger Observatory Network (AMON): Performance and Science Program

abstract: The Astrophysical Multimessenger Observatory Network (AMON) has been built with the purpose of enabling near real-time coincidence searches using data from leading multimessenger observatories and astronomical facilities. Its mission is to evoke discovery of multimessenger astrophysical sources, exploit these sources for purposes of astrophysics and fundamental physics, and explore multimessenger datasets for evidence of multimessenger source population AMON aims to promote the advancement of multimessenger astrophysics by allowing its participants to study the most energetic phenomena in the universe and to help answer some of the outstanding enigmas in astrophysics, fundamental physics, and cosmology. The main strength of AMON is its ability to combine and analyze sub-threshold data from different facilities. Such data cannot generally be used stand-alone to identify astrophysical sources. The analyses algorithms used by AMON can identify statistically significant coincidence candidates of multimessenger events, leading to the distribution of AMON alerts used by partner observatories for real-time follow-up that may identify and, potentially, confirm the reality of the multimessenger association. We present the science motivation, partner observatories, implementation and summary of the current status of the AMON project.


title: Multidimensional configurations in the primes with shifted prime steps

abstract: Let $\mathcal{P}$ denote the set of primes. For a fixed dimension $d$, Cook-Magyar-Titichetrakun, Tao-Ziegler and Fox-Zhao independently proved that any subset of positive relative density of $\mathcal{P}^d$ contains an arbitrary linear configuration. In this paper, we prove that there exists such configuration with the step being a shifted prime (prime minus $1$ or plus $1$).


title: Is Basketball a Game of Runs?

abstract: Basketball is often referred to as "a game of runs." We investigate the appropriateness of this claim using data from the full NBA 2016-17 season, comparing actual longest runs of scoring events to what long run theory predicts under the assumption that team "momentum" is not present. We provide several different variations of the analysis. Our results consistently indicate that the lengths of longest runs in NBA games are no longer than those that would occur naturally when scoring events are generated by a random process, rather than one that is influenced by "momentum".


title: An iterative staggered scheme for phase field brittle fracture propagation with stabilizing parameters

abstract: This paper concerns the analysis and implementation of a novel iterative staggered scheme for quasi-static brittle fracture propagation models, where the fracture evolution is tracked by a phase field variable. The model we consider is a two-field variational inequality system, with the phase field function and the elastic displacements of the solid material as independent variables. Using a penalization strategy, this variational inequality system is transformed into a variational equality system, which is the formulation we take as the starting point for our algorithmic developments. The proposed scheme involves a partitioning of this model into two subproblems; phase field and mechanics, with added stabilization terms to both subproblems for improved efficiency and robustness. We analyze the convergence of the proposed scheme using a fixed point argument, and find that under a natural condition, the elastic mechanical energy remains bounded, and, if the diffusive zone around crack surfaces is sufficiently thick, monotonic convergence is achieved. Finally, the proposed scheme is validated numerically with several bench-mark problems.


title: CRAFT: A multifunction online platform for speech prosody visualisation

abstract: There are many research tools which are also used for teaching the acoustic phonetics of speech rhythm and speech melody. But they were not purpose-designed for teaching-learning situations, and some have a steep learning curve. CRAFT (Creation and Recovery of Amplitude and Frequency Tracks) is custom-designed as a novel flexible online tool for visualisation and critical comparison of functions and transforms, with implementations of the Reaper, RAPT, PyRapt, YAAPT, YIN and PySWIPE F0 estimators, three Praat configurations, and two purpose-built estimators, PyAMDF, S0FT. Visualisations of amplitude and frequency envelope spectra, spectral edge detection of rhythm zones, and a parametrised spectrogram are included. A selection of audio clips from tone and intonation languages is provided for demonstration purposes. The main advantages of online tools are consistency (users have the same version and the same data selection), interoperability over different platforms, and ease of maintenance. The code is available on GitHub.


title: Nonsequenceable Steiner triple systems

abstract: A partial Steiner triple system is is $sequenceable$ if the points can be sequenced so that no proper segment can be partitioned into blocks. We show that, if $0 \leq a \leq (n-1)/3$, then there exists a nonsequenceable PSTS$(n)$ of size $\frac{1}{3}\binom{n}{2}-a$, for all $n \equiv 1 \pmod{6}$ except for $n=7$.


title: Adjoint-based SQP Method with Block-wise quasi-Newton Jacobian Updates for Nonlinear Optimal Control

abstract: Nonlinear model predictive control~(NMPC) generally requires the solution of a non-convex optimization problem at each sampling instant under strict timing constraints, based on a set of differential equations that can often be stiff and/or that may include implicit algebraic equations. This paper provides a local convergence analysis for the recently proposed adjoint-based sequential quadratic programming~(SQP) algorithm that is based on a block-structured variant of the two-sided rank-one~(TR1) quasi-Newton update formula to efficiently compute Jacobian matrix approximations in a sparsity preserving fashion. A particularly efficient algorithm implementation is proposed in case an implicit integration scheme is used for discretization of the optimal control problem, in which matrix factorization and matrix-matrix operations can be avoided entirely. The convergence analysis results as well as the computational performance of the proposed optimization algorithm are illustrated for two simulation case studies of nonlinear MPC.


title: How to glide in Schwarzschild spacetime

abstract: We investigate the motion of extended test objects in the Schwarzschild spacetime, particularly the radial fall of two point masses connected by a massless rod of a length given as a fixed, periodic function of time. We argue that such a model is inappropriate in the most interesting regimes of high and low oscillation frequencies.


title: Efficient parametric frequency conversions in lithium niobate nanophotonic chips

abstract: Chip-integrated nonlinear photonics holds the key for advanced optical information processing with superior performance and novel functionalities. Here, we present an optimally mode-matched, periodically poled lithium niobate nanowaveguide for efficient parametric frequency conversions on chip. Using a 4-mm nanowaveguide with subwavelength mode confinement, we demonstrate second harmonic generation with efficiency over $2200\%~W^{-1}cm^{-2}$, and broadband difference frequency generation with similar efficiency over a 4.5-THz spectral span. These allow us to generate correlated photon pairs over multiple frequency channels via spontaneous parametric down conversion, all in their fundamental spatial modes, with a coincidence to accidental ratio as high as 600. The high efficiency and dense integrability of the present chip devices may pave a viable route to scalable nonlinear applications in both classical and quantum domains.


title: Characterizations of equivariant Steiner and linear isometries operads

abstract: We study the indexing systems that correspond to equivariant Steiner and linear isometries operads. When $G$ is a finite abelian group, we prove that a $G$-indexing system is realized by a Steiner operad if and only if it is generated by cyclic $G$-orbits. When $G$ is a finite cyclic group, whose order is either a prime power or a product of two distinct primes greater than $3$, we prove that a $G$-indexing system is realized by a linear isometries operad if and only if it satisfies Blumberg and Hill's horn-filling condition. We also repackage the data in an indexing system as a certain kind of partial order. We call these posets transfer systems, and we develop basic tools for computing with them.


title: Stellar populations in hosts of giant radio galaxies and their neighbouring galaxies

abstract: Context: Giant radio galaxies (with projected linear size of radio structure larger than 0.7 Mpc) are very rare and unusual objects. Only $\sim$5% of extended radio sources reach such sizes. Understanding of the processes responsible for their large sizes is crucial to further our knowledge about the radio source's evolution. Aims: We check the hypothesis that giants become extremely large due to the specific history of their host galaxy formation, as well as in the context of the cluster or group of galaxies where they evolve. Therefore we study the star formation histories in their host galaxies and in galaxies located in their neighbourhood. Methods: We studied 41 giant-size radio galaxies as well as galaxies located within a radius of 5 Mpc around giants to verify whether the external conditions of the intergalactic medium somehow influence the internal evolution of galaxies in the group/cluster. We compared the results with a control sample of smaller-sized Fanaroff--Riley type II radio galaxies and their neighbouring galaxies. We fit stellar continua in all galaxy spectra using the spectral synthesis code STARLIGHT and provide statistical analysis of the results. Results: We find that hosts of giant radio galaxies have a larger amount of intermediate age stellar populations compared with smaller-sized FRII radio sources. The same result is also visible when we compare neighbouring galaxies located up to 1.5 Mpc around giants and FRIIs. This may be evidence that star formation in groups with giants was triggered due to global processes occurring in the ambient intergalactic medium. These processes may also contribute to mechanisms responsible for the extremely large sizes of giants.


title: Bourdieu, networks, and movements: Using the concepts of habitus, field and capital to understand a network analysis of gender differences in undergraduate physics

abstract: Current trends suggest that significant gender disparities exist within Science, Technology, Engineering, and Mathematics (STEM) education at university, with female students being underrepresented in physics, but more equally represented in life sciences (e.g., biology, medicine). To understand these trends, it is important to consider the context in which students make decisions about which university courses to enrol in. The current study seeks to investigate gender differences in STEM through a unique approach that combines network analysis of student enrolment data with an interpretive lens based on the sociological theory of Pierre Bourdieu. We generate a network of courses taken by around 9000 undergraduate physics students (from 2009 to 2014) to quantify Bourdieu's concept of field. We explore the properties of this network to investigate gender differences in transverse movements (between different academic fields) and vertical movements (changes in students' achievement rankings within a field). Our findings indicate that female students are more likely to make transverse movements into life science fields. We also find that university physics does a poor job in attracting high achieving students, and especially high achieving female students. Of the students who do choose to study physics, low achieving female students are less likely to continue than their male counterparts. The results and implications are discussed in the context of Bourdieu's theory, and previous research. We argue that in order to remove constraints on female student's study choices, the field of physics needs to provide a culture in which all students feel like they belong.


title: Observation of $C\!P$ violation in charm decays

abstract: A search for charge-parity ($C\!P$) violation in $D^0 \to K^- K^+$ and $D^0 \to \pi^- \pi^+$ decays is reported, using $pp$ collision data corresponding to an integrated luminosity of 6 $\mathrm{fb}^{-1}$ collected at a center-of-mass energy of 13 TeV with the LHCb detector. The flavor of the charm meson is inferred from the charge of the pion in $D^*(2010)^+ \to D^0 \pi^+$ decays or from the charge of the muon in $\overline{B} \to D^0 \mu^-\bar{\nu}_\mu X$ decays. The difference between the $C\!P$ asymmetries in $D^0 \to K^- K^+$ and $D^0 \to \pi^- \pi^+$ decays is measured to be $\Delta A_{C\!P} = [ -18.2 \pm 3.2\,(\rm stat.) \pm 0.9\,(\rm syst.) ] \times 10^{-4}$ for $\pi$-tagged and $\Delta A_{C\!P} = [ -9 \pm 8\,(\rm stat.) \pm 5\,(\rm syst.) ] \times 10^{-4} $ for $\mu$-tagged $D^0$ mesons. Combining these with previous LHCb results leads to $$\Delta A_{C\!P} = ( -15.4 \pm 2.9) \times 10^{-4},$$ where the uncertainty includes both statistical and systematic contributions. The measured value differs from zero by more than five standard deviations. This is the first observation of $C\!P$ violation in the decay of charm hadrons.


title: A stochastic Gronwall inequality and applications to moments, strong completeness, strong local Lipschitz continuity, and perturbations

abstract: There are numerous applications of the classical (deterministic) Gronwall inequality. Recently, Michael Scheutzow discovered a stochastic Gronwall inequality which provides upper bounds for $p$-th moments, $p\in(0,1)$, of the supremum of nonnegative scalar continuous processes which satisfy a linear integral inequality. In this article we complement this with upper bounds for $p$-th moments, $p\in[2,\infty)$, of the supremum of general It\^o processes which satisfy a suitable one-sided affine-linear growth condition. As example applications, we improve known results on strong local Lipschitz continuity in the starting point of solutions of stochastic differential equations (SDEs), on (exponential) moment estimates for SDEs, on strong completeness of SDEs, and on perturbation estimates for SDEs.


title: A new conservative/dissipative time integration scheme for nonlinear mechanical systems

abstract: We present a conservative/dissipative time integration scheme for nonlinear mechanical systems. Starting from a weak form, we derive algorithmic forces and velocities that guarantee the desired conservation/dissipation properties. Our approach relies on a collection of linearly constrained quadratic programs defining high order correction terms that modify, in the minimum possible way, the classical midpoint rule so as to guarantee the strict energy conservation/dissipation properties. The solution of these programs provides explicit formulas for the algorithmic forces and velocities which can be easily incorporated into existing implementations. Similarities and differences between our approach and well-established methods are discussed as well. The approach, suitable for reduced-order models, finite element models, or multibody systems, is tested and its capabilities are illustrated by means of several examples.


title: Relativistic fluid dynamics and its extensions as an effective field theory

abstract: We examine hydrodynamics from the perspective of an effective field theory. The microscopic scale in this case is the thermalization scale, and the macroscopic scale is the gradient, with thermal fluctuations playing the role of $\hbar$. We argue that this method can be applied both, to consistently include thermal fluctuations in the theory, and to extend hydrodynamics to systems whose microscopic structure is non-trivial. For the latter, we discuss the case of spin polarization and gauge theories.


title: A characterization of the U(Omega,m) sets of a hyperelliptic curve as Omega and m vary

abstract: In this article we consider a certain distinguished set $U(\Omega,m) \subseteq \{1,2,\ldots,2g+1,\infty\}$ that can be attached to a marked hyperelliptic curve of genus $g$ equipped with a small period matrix $\Omega$ for its polarized Jacobian. We show that as $\Omega$ and the marking $m$ vary, this set ranges over all possibilities prescribed by an argument of Poor.


title: Three Convolution Inequalities on the Real Line with Connections to Additive Combinatorics

abstract: We discuss three convolution inequalities that are connected to additive combinatorics. Cloninger and the second author showed that for nonnegative $f \in L^1(-1/4, 1/4)$, $$ \max_{-1/2 \leq t \leq 1/2} \int_{\mathbb{R}}{f(t-x) f(x) dx} \geq 1.28 \left( \int_{-1/4}^{1/4}{f(x) dx}\right)^2$$ which is related to $g-$Sidon sets (1.28 cannot be replaced by 1.52). We prove a dual statement, related to difference bases, and show that for $f \in L^1(\mathbb{R})$, $$ \min_{0 \leq t \leq 1}\int_{\mathbb{R}}{f(x) f(x+t) dx} \leq 0.42 \|f\|_{L^1}^2,$$ where the constant 1/2 is trivial, 0.42 cannot be replaced by 0.37. This suggests a natural conjecture about the asymptotic structure of $g-$difference bases. Finally, we show for all functions $f \in L^1(\mathbb{R}) \cap L^2(\mathbb{R})$, $$ \int_{-\frac{1}{2}}^{\frac{1}{2}}{ \int_{\mathbb{R}}{f(x) f(x+t) dx}dt} \leq 0.91 \|f\|_{L^1}\|f\|_{L^2}$$


title: Digital Memcomputing: from Logic to Dynamics to Topology

abstract: Digital memcomputing machines (DMMs) are a class of computational machines designed to solve combinatorial optimization problems. A practical realization of DMMs can be accomplished via electrical circuits of highly non-linear, point-dissipative dynamical systems engineered so that periodic orbits and chaos can be avoided. A given logic problem is first mapped into this type of dynamical system whose point attractors represent the solutions of the original problem. A DMM then finds the solution via a succession of elementary instantons whose role is to eliminate solitonic configurations of logical inconsistency ("logical defects") from the circuit. By employing a supersymmetric theory of dynamics, a DMM can be described by a cohomological field theory that allows for computation of certain topological matrix elements on instantons that have the mathematical meaning of intersection numbers on instantons. We discuss the "dynamical" meaning of these matrix elements, and argue that the number of elementary instantons needed to reach the solution cannot exceed the number of state variables of DMMs, which in turn can only grow at most polynomially with the size of the problem. These results shed further light on the relation between logic, dynamics and topology in digital memcomputing.


title: OGLE-2018-BLG-0022: A Nearby M-dwarf Binary

abstract: We report observations of the binary microlensing event OGLE-2018-BLG-0022, provided by the ROME/REA Survey, which indicate that the lens is a low-mass binary star consisting of M3 (0.375+/-0.020 Msun) and M7 (0.098+/-0.005 Msun) components. The lens is unusually close, at 0.998+/-0.047 kpc, compared with the majority of microlensing events, and despite its intrinsically low luminosity, it is likely that AO observations in the near future will be able to provide an independent confirmation of the lens masses.


title: Combination of multiple Deep Learning architectures for Offensive Language Detection in Tweets

abstract: This report contains the details regarding our submission to the OffensEval 2019 (SemEval 2019 - Task 6). The competition was based on the Offensive Language Identification Dataset. We first discuss the details of the classifier implemented and the type of input data used and pre-processing performed. We then move onto critically evaluating our performance. We have achieved a macro-average F1-score of 0.76, 0.68, 0.54, respectively for Task a, Task b, and Task c, which we believe reflects on the level of sophistication of the models implemented. Finally, we will be discussing the difficulties encountered and possible improvements for the future.


title: Finite element theory on curved domains with applications to DGFEMs

abstract: In this paper we provide key estimates used in the stability and error analysis of discontinuous Galerkin finite element methods (DGFEMs) on domains with curved boundaries. In particular, we review trace estimates, inverse estimates, discrete Poincar\'e--Friedrichs' inequalities, and optimal interpolation estimates in noninteger Hilbert-Sobolev norms, that are well known in the case of polytopal domains. We also prove curvature bounds for curved simplices, which does not seem to be present in the existing literature, even in the polytopal setting, since polytopal domains have piecewise zero curvature. We demonstrate the value of these estimates, by analysing the IPDG method for the Poisson problem, introduced by Douglas and Dupont [\emph{Computing Methods in Applied Sciences, Lecture Notes in Physics, vol 58. Springer, Berlin, Heidelberg}, pages 207--216. Springer, 1976], and by analysing a variant of the $hp$-DGFEM for the biharmonic problem introduced by Mozolevski and S\"{u}li [\emph{Computer Methods in Applied Mechanics and Engineering}, 196(13-16):1851--1863, 2007]. In both cases we prove stability estimates and optimal a priori error estimates. Numerical results are provided, validating the proven error estimates.


title: Notes on Markov embedding

abstract: The representation problem of finite-dimensional Markov matrices in Markov semigroups is revisited, with emphasis on concrete criteria for matrix subclasses of theoretical or practical relevance, such as equal-input, circulant, symmetric or doubly stochastic matrices. Here, we pay special attention to various algebraic properties of the embedding problem, and discuss the connection with the centraliser of a Markov matrix.


title: Virtual concordance and the generalized Alexander polynomial

abstract: We use the Bar-Natan Zh-correspondence to identify the generalized Alexander polynomial of a virtual knot with the Alexander polynomial of a two component welded link. We show that the Zh-map is functorial under concordance, and also that Satoh's Tube map (from welded links to ribbon knotted tori in $S^4$) is functorial under concordance. In addition, we extend classical results of Chen, Milnor, and Hillman on the lower central series of link groups to links in thickened surfaces. Our main result is that the generalized Alexander polynomial vanishes on any knot in a thickened surface which is concordant to a homologically trivial knot. In particular, this shows that it vanishes on virtually slice knots. We apply it to complete the calculation of the slice genus for virtual knots with four crossings and to determine non-sliceness for a number of 5-crossing and 6-crossing virtual knots.


title: Batch Policy Learning under Constraints

abstract: When learning policies for real-world domains, two important questions arise: (i) how to efficiently use pre-collected off-policy, non-optimal behavior data; and (ii) how to mediate among different competing objectives and constraints. We thus study the problem of batch policy learning under multiple constraints, and offer a systematic solution. We first propose a flexible meta-algorithm that admits any batch reinforcement learning and online learning procedure as subroutines. We then present a specific algorithmic instantiation and provide performance guarantees for the main objective and all constraints. To certify constraint satisfaction, we propose a new and simple method for off-policy policy evaluation (OPE) and derive PAC-style bounds. Our algorithm achieves strong empirical results in different domains, including in a challenging problem of simulated car driving subject to multiple constraints such as lane keeping and smooth driving. We also show experimentally that our OPE method outperforms other popular OPE techniques on a standalone basis, especially in a high-dimensional setting.


title: Russian Language Datasets in the Digitial Humanities Domain and Their Evaluation with Word Embeddings

abstract: In this paper, we present Russian language datasets in the digital humanities domain for the evaluation of word embedding techniques or similar language modeling and feature learning algorithms. The datasets are split into two task types, word intrusion and word analogy, and contain 31362 task units in total. The characteristics of the tasks and datasets are that they build upon small, domain-specific corpora, and that the datasets contain a high number of named entities. The datasets were created manually for two fantasy novel book series ("A Song of Ice and Fire" and "Harry Potter"). We provide baseline evaluations with popular word embedding models trained on the book corpora for the given tasks, both for the Russian and English language versions of the datasets. Finally, we compare and analyze the results and discuss specifics of Russian language with regards to the problem setting.


title: Gaussian wave packet transform based numerical scheme for the semi-classical Schr\"odinger equation with random inputs

abstract: In this work, we study the semi-classical limit of the Schr\"odinger equation with random inputs, and show that the semi-classical Schr\"odinger equation produces $O(\varepsilon)$ oscillations in the random variable space. With the Gaussian wave packet transform, the original Schr\"odinger equation is mapped to an ODE system for the wave packet parameters coupled with a PDE for the quantity $w$ in rescaled variables. Further, we show that the $w$ equation does not produce $\varepsilon$ dependent oscillations, and thus it is more amenable for numerical simulations. We propose multi-level sampling strategy in implementing the Gaussian wave packet transform, where in the most costly part, i.e. simulating the $w$ equation, it is sufficient to use $\varepsilon$ independent samples. We also provide extensive numerical tests as well as meaningful numerical experiments to justify the properties of the numerical algorithm, and hopefully shed light on possible future directions.


title: A parametric acceleration of multilevel Monte Carlo convergence for nonlinear variably saturated flow

abstract: We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that are modeled using the Richards' equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards' equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards' equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared to the original MC estimator.


title: Noisy Accelerated Power Method for Eigenproblems with Applications

abstract: This paper introduces an efficient algorithm for finding the dominant generalized eigenvectors of a pair of symmetric matrices. Combining tools from approximation theory and convex optimization, we develop a simple scalable algorithm with strong theoretical performance guarantees. More precisely, the algorithm retains the simplicity of the well-known power method but enjoys the asymptotic iteration complexity of the powerful Lanczos method. Unlike these classic techniques, our algorithm is designed to decompose the overall problem into a series of subproblems that only need to be solved approximately. The combination of good initializations, fast iterative solvers, and appropriate error control in solving the subproblems lead to a linear running time in the input sizes compared to the superlinear time for the traditional methods. The improved running time immediately offers acceleration for several applications. As an example, we demonstrate how the proposed algorithm can be used to accelerate canonical correlation analysis, which is a fundamental statistical tool for learning of a low-dimensional representation of high-dimensional objects. Numerical experiments on real-world data sets confirm that our approach yields significant improvements over the current state-of-the-art.


title: Phase transition in random contingency tables with non-uniform margins

abstract: For parameters $n,\delta,B,$ and $C$, let $X=(X_{k\ell})$ be the random uniform contingency table whose first $\lfloor n^{\delta} \rfloor $ rows and columns have margin $\lfloor BCn \rfloor$ and the last $n$ rows and columns have margin $\lfloor Cn \rfloor$. For every $0<\delta<1$, we establish a sharp phase transition of the limiting distribution of each entry of $X$ at the critical value $B_{c}=1+\sqrt{1+1/C}$. In particular, for $1/2<\delta<1$, we show that the distribution of each entry converges to a geometric distribution in total variation distance, whose mean depends sensitively on whether $B<B_{c}$ or $B>B_{c}$. Our main result shows that $\mathbb{E}[X_{11}]$ is uniformly bounded for $B<B_{c}$, but has sharp asymptotic $C(B-B_{c}) n^{1-\delta}$ for $B>B_{c}$. We also establish a strong law of large numbers for the row sums in top right and top left blocks.


title: The Unexpected Spectrum of the Innermost Ejecta of the Red Hypergiant VY CMa

abstract: HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm VO as a circumstellar molecule. In this letter we discuss the spectra of the features, their motions and ages, and the identification of the molecular emission. The strength of the atomic and molecular features in the small clumps present an astrophysical problem for the excitation process. We show that the clumps must have a nearly clear line of sight to the star's radiation.


title: Reconciling dark matter, $R_{K^{(*)}}$ anomalies and $(g-2)_{\mu}$ in an ${L_{\mu}-L_{\tau}}$ scenario

abstract: We propose an anomaly free unified scenario by invocation of an extra local ${\rm U(1)}_{L_{\mu}-L_{\tau}}$ gauge symmetry. This scenario simultaneously resolves the $R_{K^{(*)}}$ anomalies, the dark matter puzzle and the long-standing discrepancy in muon's anomalous magnetic moment. A complex scalar ($\eta$) having nonzero ${L_{\mu}-L_{\tau}}$ charge has been introduced to break this new U(1) symmetry spontaneously. Moreover, for the purpose of studying dark matter phenomenology and $R_{K^{(*)}}$ anomalies in a correlated manner, we introduce an inert ${\rm SU(2)}_L$ scalar doublet ($\Phi$), a $\mathbb{Z}_2$-odd real singlet scalar ($S$) and a $\mathbb{Z}_2$-odd coloured fermion ($\chi$) which transforms vectorially under the ${\rm U(1)}_{L_{\mu}-L_{\tau}}$ symmetry. This extra gauge symmetry provides a new gauge boson $Z_{\mu\tau}$ which not only gives additional contribution to both $b\to s \ell\ell$ transition and $(g-2)_{\mu}$ but also provides a crucial annihilation channel for dark matter candidate $\rho_1$ of the present scenario. This $\rho_1$ is an admixture of CP-even neutral component of $\Phi$ and $S$. Our analysis shows that the low mass dark matter regime ($M_{\rho_1}\lesssim 60$ GeV) is still allowed by the experiments like XENON1T, LHC (via Higgs invisible branching) and Fermi-LAT, making the dark matter phenomenology drastically different from the standard Inert Doublet and the Scalar Singlet models. Furthermore, the present model is also fairly consistent with the observed branching ratio of $B\to X_s\gamma$ in $3\sigma$ range and is quite capable of explaining neutrino masses and mixings via Type-I seesaw mechanism if we add three right handed neutrinos in the particle spectrum. Finally, we use the latest ATLAS data of non-observation of a resonant $\ell^+\ell^-$ signal at the 13 TeV LHC to constrain the mass-coupling plane of $Z_{\mu\tau}$.


title: Affordance Learning In Direct Perception for Autonomous Driving

abstract: Recent development in autonomous driving involves high-level computer vision and detailed road scene understanding. Today, most autonomous vehicles are using mediated perception approach for path planning and control, which highly rely on high-definition 3D maps and real time sensors. Recent research efforts aim to substitute the massive HD maps with coarse road attributes. In this paper, we follow the direct perception based method to train a deep neural network for affordance learning in autonomous driving. Our goal in this work is to develop the affordance learning model based on freely available Google Street View panoramas and Open Street Map road vector attributes. Driving scene understanding can be achieved by learning affordances from the images captured by car-mounted cameras. Such scene understanding by learning affordances may be useful for corroborating base maps such as HD maps so that the required data storage space is minimized and available for processing in real time. We compare capability in road attribute identification between human volunteers and our model by experimental evaluation. Our results indicate that this method could act as a cheaper way for training data collection in autonomous driving. The cross validation results also indicate the effectiveness of our model.


title: Statistical Methods for Replicability Assessment

abstract: Large-scale replication studies like the Reproducibility Project: Psychology (RP:P) provide invaluable systematic data on scientific replicability, but most analyses and interpretations of the data fail to agree on the definition of "replicability" and disentangle the inexorable consequences of known selection bias from competing explanations. We discuss three concrete definitions of replicability based on (1) whether published findings about the signs of effects are mostly correct, (2) how effective replication studies are in reproducing whatever true effect size was present in the original experiment, and (3) whether true effect sizes tend to diminish in replication. We apply techniques from multiple testing and post-selection inference to develop new methods that answer these questions while explicitly accounting for selection bias. Our analyses suggest that the RP:P dataset is largely consistent with publication bias due to selection of significant effects. The methods in this paper make no distributional assumptions about the true effect sizes.


title: Forecasting of the Thermosphere via Assimilation of Electron Density and Temperature Data

abstract: The paper presents experiments of driving a physics-based thermosphere model by assimilating electron density (Ne) and temperature (Tn) data using the ensemble adjustment Kalman filter (EAKF) technique. This study not only helps to gauge the accuracy of the assimilation, to explain the inherent model bias, and to understand the limitations of the framework, but it also establishes EAKF as a viable technique in the presence of realistic data assimilation scenarios to forecast the highly dynamical thermosphere. The results from perfect model scenarios show that data assimilation changes and, more often than not, improves the model state. Data from Swarm-A, Swarm-C, CHAMP, and GRACE-A are used to validate the resulting analysis states. The independent validation results show that the Ne-guided thermosphere state does not outperform the model state without data assimilation along the considered orbits. This may be due to the limited number of bonafide Ne profiles available for the thermosphere specification tasks in the experiments. More importantly, the results show that the Ne-guided thermosphere state does not deteriorate much in performance during geomagnetic storm time. The results reveal a few challenges of using Ne profiles in a hypothetical operational data assimilation exercise. The experiment with assimilating Tn shows more promise over Ne in terms of estimating mass density along the orbits of both CHAMP and GRACE-A satellites. The results show that the improvement gained in the overall forecasted thermosphere state is better during solar minimum compared to that of solar maximum. These results also provide insights into the biases inherent in the physics-based model. The systematic biases that the paper highlight could be an indication that the specification of plasma-neutral interactions in the model needs further adjustments.


title: A note on higher order perturbative corrections to squirming speed in weakly viscoelastic fluids

abstract: Many microorganisms swim in fluids with complex rheological properties. Although much is now understood about motion of these swimmers in Newtonian fluids, the understanding is still developing in non-Newtonian fluids --- this understanding is crucial for various biomimetic and biomedical applications. Here we study a common model for microswimmers, the squirmer model, in two common viscoelastic fluid models, the Giesekus fluid model and fluids of differential type (grade three), at zero Reynolds number. Through this article we address a recent commentary that discussed suitable values of parameters in these model and pointed at higher order viscoelastic effects on the squirming motion.


title: Coherent transition radiation from the geomagnetically-induced current in cosmic-ray air showers: Implications for the anomalous events observed by ANITA

abstract: We show that coherent transition radiation from the electrically-neutral transverse geomagnetic current (CTR- GM) in a cosmic-ray air shower provides a natural, standard model, explanation to the recent ``anomalous'' events observed by the ANITA detector. We demonstrate that for zenith angles less than roughly 70 degrees, combined with high surface elevation, the inclusion of CTR-GM can significantly alter the emitted electric field from a cosmic-ray air shower. CTR-GM therefore has to be included in radio emission models to provide a full description of the radio emission from a high-energy cosmic-ray air shower traversing a dielectric boundary.


title: Plasma 2020 - Intracluster Medium Plasmas

abstract: Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.


title: Byzantine Fault Tolerant Distributed Linear Regression

abstract: This paper considers the problem of Byzantine fault tolerance in distributed linear regression in a multi-agent system. However, the proposed algorithms are given for a more general class of distributed optimization problems, of which distributed linear regression is a special case. The system comprises of a server and multiple agents, where each agent is holding a certain number of data points and responses that satisfy a linear relationship (could be noisy). The objective of the server is to determine this relationship, given that some of the agents in the system (up to a known number) are Byzantine faulty (aka. actively adversarial). We show that the server can achieve this objective, in a deterministic manner, by robustifying the original distributed gradient descent method using norm based filters, namely 'norm filtering' and 'norm-cap filtering', incurring an additional log-linear computation cost in each iteration. The proposed algorithms improve upon the existing methods on three levels: i) no assumptions are required on the probability distribution of data points, ii) system can be partially asynchronous, and iii) the computational overhead (in order to handle Byzantine faulty agents) is log-linear in number of agents and linear in dimension of data points. The proposed algorithms differ from each other in the assumptions made for their correctness, and the gradient filter they use.


title: Astro2020 Science White Paper: Triggered High-Priority Observations of Dynamic Solar System Phenomena

abstract: Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.


title: Stability and Error Analysis for Optimization and Generalized Equations

abstract: Stability and error analysis remain challenging for problems that lack regularity properties near solutions, are subject to large perturbations, and might be infinite dimensional. We consider nonconvex optimization and generalized equations defined on metric spaces and develop bounds on solution errors using the truncated Hausdorff distance applied to graphs and epigraphs of the underlying set-valued mappings and functions. In the process, we extend the calculus of such distances to cover compositions and other constructions that arise in nonconvex problems. The results are applied to constrained problems with feasible sets that might have empty interiors, solution of KKT systems, and optimality conditions for difference-of-convex functions and composite functions.


title: Using Ego-Clusters to Measure Network Effects at LinkedIn

abstract: A network effect is said to take place when a new feature not only impacts the people who receive it, but also other users of the platform, like their connections or the people who follow them. This very common phenomenon violates the fundamental assumption underpinning nearly all enterprise experimentation systems, the stable unit treatment value assumption (SUTVA). When this assumption is broken, a typical experimentation platform, which relies on Bernoulli randomization for assignment and two-sample t-test for assessment of significance, will not only fail to account for the network effect, but potentially give highly biased results. This paper outlines a simple and scalable solution to measuring network effects, using ego-network randomization, where a cluster is comprised of an "ego" (a focal individual), and her "alters" (the individuals she is immediately connected to). Our approach aims at maintaining representativity of clusters, avoiding strong modeling assumption, and significantly increasing power compared to traditional cluster-based randomization. In particular, it does not require product-specific experiment design, or high levels of investment from engineering teams, and does not require any changes to experimentation and analysis platforms, as it only requires assigning treatment an individual level. Each user either has the feature or does not, and no complex manipulation of interactions between users is needed. It focuses on measuring the one-out network effect (i.e the effect of my immediate connection's treatment on me), and gives reasonable estimates at a very low setup cost, allowing us to run such experiments dozens of times a year.


title: Distributed Vector Representations of Folksong Motifs

abstract: This article presents a distributed vector representation model for learning folksong motifs. A skip-gram version of word2vec with negative sampling is used to represent high quality embeddings. Motifs from the Essen Folksong collection are compared based on their cosine similarity. A new evaluation method for testing the quality of the embeddings based on a melodic similarity task is presented to show how the vector space can represent complex contextual features, and how it can be utilized for the study of folksong variation.


title: Studying Magnetic Fields in Star Formation and the Turbulent Interstellar Medium

abstract: Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.


title: Cohomologie des fibr\'es en droites sur SL3 /B en caract\'eristique positive : deux filtrations et cons\'equences

abstract: In this thesis, I will prove the existence of two filtrations of the cohomology of line bundles on SL_3/B. The first one is a two-step filtration that exists for $H^1(\mu)$ and $H^2(\mu)$ if $\mu$ is in the Griffith region. The second one exists for all $H^i(\mu)$, which is similar to the p-filtration that has been considered by Jens Carsten Jantzen.


title: Characterizing black hole metrics in quadratic gravity

abstract: The recent discovery of non-Schwarzschild black hole spacetimes has opened new directions of research in higher-derivative gravitational theories. However, despite intense analytical and numerical efforts, the link with the linearized theory is still poorly understood. In this work we address this point for the Einstein-Weyl Lagrangian, whose weak field limit is characterized by the standard massless graviton and a spin-2 ghost. We show that the strength of the Yukawa term at infinity determines the thermal properties of the black hole and the structure of the singularity near $r=0$. Moreover, inspired by recent results in the Asymptotic Safety scenario we investigate the consequences of an imaginary ghost mass. In this case we find a countable set of solutions all characterized by spatial oscillations of typical wavelength determined by the mass of the spin-2 field.


title: Nucleon Femtography from Exclusive Reactions

abstract: Major breakthroughs over the last two decades have led us to access information on how the nucleon's mass, spin and mechanical properties are generated from its quark and gluon degrees of freedom. On one side, a theoretical framework has been developed which enables the extraction of 3D parton distributions from deeply virtual exclusive scattering experiments. On the other hand, the so called gravitomagnetic form factors parameterizing the QCD energy momentum tensor of the nucleon have been connected to the Mellin moments of the 3D parton distributions. Current efforts in both experiment and theory are being directed at using information from electron scattering experiments to map out and eventually visualize these 3D distributions as well as the mass, spin and mechanical properties of the nucleon and of spin 0, 1/2, and 1 nuclei. A new science of nucleon femtography is emerging where the 3D structure of the nucleon will be studied merging information from current and forthcoming data from various facilities including the future EIC and the new influx of data science, imaging, and visualization.


title: Chaos in Three-dimensional Higher Spin Gravity

abstract: We discuss $SL(N,\mathbb{C})$ Chern-Simons higher spin gauge theories in Euclidean AdS$_3$. With appropriate boundary term, we derive the higher spin generalization of Schwarzian on-shell action. We investigate gravitationally dressed bi-locals, and we study the soft higher spin mode expansion to obtain soft mode eigenfunctions. We also derive the spin-$s$ eigenfunction from $\mathcal{W}$-Ward identity and a recursion relation. Using the on-shell action, we evaluate the contributions of the soft higher spin modes to the out-of-time-order correlators, and the corresponding Lyapunov exponent of $SL(N)$ higher spin gravity is found to be ${2\pi \over \beta}(N-1)$.


title: Large-Scale Online Experimentation with Quantile Metrics

abstract: Online experimentation (or A/B testing) has been widely adopted in industry as the gold standard for measuring product impacts. Despite the wide adoption, few literatures discuss A/B testing with quantile metrics. Quantile metrics, such as 90th percentile page load time, are crucial to A/B testing as many key performance metrics including site speed and service latency are defined as quantiles. However, with LinkedIn's data size, quantile metric A/B testing is extremely challenging because there is no statistically valid and scalable variance estimator for the quantile of dependent samples: the bootstrap estimator is statistically valid, but takes days to compute; the standard asymptotic variance estimate is scalable but results in order-of-magnitude underestimation. In this paper, we present a statistically valid and scalable methodology for A/B testing with quantiles that is fully generalizable to other A/B testing platforms. It achieves over 500 times speed up compared to bootstrap and has only $2\%$ chance to differ from bootstrap estimates. Beyond methodology, we also share the implementation of a data pipeline using this methodology and insights on pipeline optimization.


title: A Systematic and Meta-analysis Survey of Whale Optimization Algorithm

abstract: Whale Optimization Algorithm (WOA) is a nature-inspired meta-heuristic optimization algorithm, which was proposed by Mirjalili and Lewis in 2016. This algorithm has shown its ability to solve many problems. Comprehensive surveys have been conducted about some other nature-inspired algorithms, such as ABC, PSO, etc.Nonetheless, no survey search work has been conducted on WOA. Therefore, in this paper, a systematic and meta analysis survey of WOA is conducted to help researchers to use it in different areas or hybridize it with other common algorithms. Thus, WOA is presented in depth in terms of algorithmic backgrounds, its characteristics, limitations, modifications, hybridizations, and applications. Next, WOA performances are presented to solve different problems. Then, the statistical results of WOA modifications and hybridizations are established and compared with the most common optimization algorithms and WOA. The survey's results indicate that WOA performs better than other common algorithms in terms of convergence speed and balancing between exploration and exploitation. WOA modifications and hybridizations also perform well compared to WOA. In addition, our investigation paves a way to present a new technique by hybridizing both WOA and BAT algorithms. The BAT algorithm is used for the exploration phase, whereas the WOA algorithm is used for the exploitation phase. Finally, statistical results obtained from WOA-BAT are very competitive and better than WOA in 16 benchmarks functions. WOA-BAT also outperforms well in 13 functions from CEC2005 and 7 functions from CEC2019.


title: Generalized Framework for Nonlinear Acceleration

abstract: Nonlinear acceleration algorithms improve the performance of iterative methods, such as gradient descent, using the information contained in past iterates. However, their efficiency is still not entirely understood even in the quadratic case. In this paper, we clarify the convergence analysis by giving general properties that share several classes of nonlinear acceleration: Anderson acceleration (and variants), quasi-Newton methods (such as Broyden Type-I or Type-II, SR1, DFP, and BFGS) and Krylov methods (Conjugate Gradient, MINRES, GMRES). In particular, we propose a generic family of algorithms that contains all the previous methods and prove its optimal rate of convergence when minimizing quadratic functions. We also propose multi-secants updates for the quasi-Newton methods listed above. We provide a Matlab code implementing the algorithm.