title: On monotonic determinacy and rewritability for recursive queries and views

abstract: A query Q is monotonically determined over a set of views if Q can be expressed as a monotonic function of the view image. In the case of relational algebra views and queries, monotonic determinacy coincides with rewritability as a union of conjunctive queries, and it is decidable in important special cases, such as for CQ views and queries. We investigate the situation for views and queries in the recursive query language Datalog. We give both positive and negative results about the ability to decide monotonic determinacy, and also about the co-incidence of monotonic determinacy with Datalog rewritability.


title: $SO(2)$ gauged Skyrmions in $4+1$ dimensions

abstract: We study the simplest $SO(2)$ gauged $O(5)$ Skyrme models in $4+1$ (flat) dimensions. In the gauge decoupled limit, the model supports topologically stable solitons (Skyrmions) and after gauging, the static energy of the solutions is bounded from below by a "baryon number". The studied model features both Maxwell and Maxwell--Chern-Simons dynamics. The considered configurations are subject to bi-azimuthal symmetry in the ${\mathbb R}^4$ subspace resulting in a two dimensional subsystem, as well as subject to an enhanced symmetry relating the two planes in the ${\mathbb R}^4$ subspace, which results in a one dimensional subsystem. Numerical solutions are constructed in both cases. In the purely magnetic case, fully bi-azimuthal solutions were given, while electrically charged and spinning solutions were constructed only in the radial (enhanced symmetric) case, both in the presence of a Chern-Simons term, and in its absence. We find that, in contrast with the analogous models in $2+1$ dimensions, the presence of the Chern-Simons term in the model under study here results only in quantitative effects.


title: Inhibition and Set-Shifting Tasks in Central Executive Function of Working Memory: An Event-Related Potential (ERP) Study

abstract: Understanding of neuro-dynamics of a complex higher cognitive process, Working Memory (WM) is challenging. In WM, information processing occurs through four subsystems: phonological loop, visual sketch pad, memory buffer and central executive function (CEF). CEF plays a principal role in WM. In this study, our objective was to understand the neurospatial correlates of CEF during inhibition and set-shifting processes. Thirty healthy educated subjects were selected. Event-Related Potential (ERP) related to visual inhibition and set-shifting task was collected using 32 channel EEG system. Activation of those ERPs components was analyzed using amplitudes of positive and negative peaks. Experiment was controlled using certain parametric constraints to judge behavior, based on average responses in order to establish relationship between ERP and local area of brain activation and represented using standardized low resolution brain electromagnetic tomography. The average score of correct responses was higher for inhibition task (87.5%) as compared to set-shifting task (59.5%). The peak amplitude of neuronal activity for inhibition task was lower compared to set-shifting task in fronto-parieto-central regions. Hence this proposed paradigm and technique can be used to measure inhibition and set-shifting neuronal processes in understanding pathological central executive functioning in patients with neuro-psychiatric disorders.


title: Quantum-Clustered Two-Photon Walks

abstract: We demonstrate a previously unknown two-photon effect in a discrete-time quantum walk. Two identical bosons with no mutual interactions nonetheless can remain clustered together as they walk on a lattice of directionally-reversible optical four-ports acting as Grover coins; both photons move in the same direction at each step due to a two-photon quantum interference phenomenon reminiscent of the Hong-Ou-Mandel effect. The clustered two-photon amplitude splits into two localized parts, one oscillating near the initial point, and the other moving ballistically without spatial spread, in soliton-like fashion. But the two photons are always clustered in the same part of the superposition, leading to potential applications for transport of entanglement and opportunities for novel two-photon interferometry experiments.


title: Dust evolution across the Horsehead Nebula

abstract: Micro-physical processes on interstellar dust surfaces are tightly connected to dust properties (i.e. dust composition, size and shape) and play a key role in numerous phenomena in the interstellar medium (ISM). The large disparity in physical conditions (i.e. density, gas temperature) in the ISM triggers an evolution of dust properties. The analysis of how dust evolves with the physical conditions is a stepping-stone towards a more thorough understanding of interstellar dust. The aim of this paper is to highlight dust evolution in the Horsehead Nebula PDR region. We use Spitzer/IRAC (3.6, 4.5, 5.8 and 8 {\mu}m), Spitzer/MIPS (24 {\mu}m) together with Herschel/PACS (70 and 160 {\mu}m) and Herschel/SPIRE (250, 350 and 500 {\mu}m) to map the spatial distribution of dust in the Horsehead over the entire emission spectral range. We model dust emission and scattering using the THEMIS interstellar dust model together with the 3D radiative transfer code SOC. We find that the nano-grains dust-to-gas ratio in the irradiated outer part of the Horsehead is 6 to 10 times lower than in the diffuse ISM. Their minimum size is 2 to 2.25 times larger than in the diffuse ISM and the power-law exponent of their size distribution, 1.1 to 1.4 times lower than in the diffuse ISM. Regarding the denser part of the Horsehead, it is necessary to use evolved grains (i.e. aggregates, with or without an ice mantle). It is not possible to explain the observations using grains from the diffuse medium. We therefore propose the following scenario to explain our results. In the outer part of the Horsehead, all the nano-grains have not yet had time to re-form completely through photo-fragmentation of aggregates and the smallest of the nano-grains that are sensitive to the radiation field are photo-destroyed. In the inner part of the Horsehead, grains most likely consist of multi-compositional, mantled aggregates.


title: A cow structural model for video analytics of cow health

abstract: In livestock farming, animal health directly influences productivity. For dairy cows, many health conditions can be evaluated by trained observers based on visual appearance and movement. However, to manually evaluate every cow in a commercial farm is expensive and impractical. This paper introduces a video-analytic system which automatically detects the cow structure from captured video sequences. A side-view cow structural model is designed to describe the spatial positions of the joints (keypoints) of the cow, and we develop a system using deep learning to automatically extract the structural model from videos. The proposed detection system can detect multiple cows in the same frame and provide robust performance under practical challenges like obstacles (fences) and poor illumination. Compared to other object detection methods, this system provides better detection results and successfully isolates the keypoints of each cow even when they are close to each other.


title: Direct determination of the Tomonaga-Luttinger parameter $K$ in quasi-one-dimensional spin systems

abstract: We test the analytical formula for the enhancement of the nuclear magnetic resonance rate $T_1^{-1}$ by the critical spin fluctuations, over the simple power-law dependence predicted for a purely one-dimensional spin system, recently derived in the random phase approximation [M. Dupont et al., Phys. Rev. B 98, 094403 (2018)]. This prediction is experimentally confirmed by excellent fits to the published temperature dependence of $T_1^{-1}$ data in the two representative spin compounds, (C$_7$H$_{10}$N)$_2$CuBr$_4$ (DIMPY) and BaCo$_2$V$_2$O$_8$, providing at the same time a direct and convenient experimental determination of the Tomonaga-Luttinger-liquid parameter $K$, very well in agreement with theoretical predictions.


title: Cascade EF-GAN: Progressive Facial Expression Editing with Local Focuses

abstract: Recent advances in Generative Adversarial Nets (GANs) have shown remarkable improvements for facial expression editing. However, current methods are still prone to generate artifacts and blurs around expression-intensive regions, and often introduce undesired overlapping artifacts while handling large-gap expression transformations such as transformation from furious to laughing. To address these limitations, we propose Cascade Expression Focal GAN (Cascade EF-GAN), a novel network that performs progressive facial expression editing with local expression focuses. The introduction of the local focus enables the Cascade EF-GAN to better preserve identity-related features and details around eyes, noses and mouths, which further helps reduce artifacts and blurs within the generated facial images. In addition, an innovative cascade transformation strategy is designed by dividing a large facial expression transformation into multiple small ones in cascade, which helps suppress overlapping artifacts and produce more realistic editing while dealing with large-gap expression transformations. Extensive experiments over two publicly available facial expression datasets show that our proposed Cascade EF-GAN achieves superior performance for facial expression editing.


title: Moments of the logarithmic derivative of characteristic polynomials from $SO(N)$ and $USp(2N)$

abstract: We study moments of the logarithmic derivative of characteristic polynomials of orthogonal and symplectic random matrices. In particular, we compute the asymptotics for large matrix size, $N$, of these moments evaluated at points which are approaching 1. This follows work of Bailey, Bettin, Blower, Conrey, Prokhorov, Rubinstein and Snaith where they compute these asymptotics in the case of unitary random matrices.


title: Optimal HDR and Depth from Dual Cameras

abstract: Dual camera systems have assisted in the proliferation of various applications, such as optical zoom, low-light imaging and High Dynamic Range (HDR) imaging. In this work, we explore an optimal method for capturing the scene HDR and disparity map using dual camera setups. Hasinoff et al. (2010) have developed a noise optimal framework for HDR capture from a single camera. We generalize this to the dual camera set-up for estimating both HDR and disparity map. It may seem that dual camera systems can capture HDR in a shorter time. However, disparity estimation is a necessary step, which requires overlap among the images captured by the two cameras. This may lead to an increase in the capture time. To address this conflicting requirement, we propose a novel framework to find the optimal exposure and ISO sequence by minimizing the capture time under the constraints of an upper bound on the disparity error and a lower bound on the per-exposure SNR. We show that the resulting optimization problem is non-convex in general and propose an appropriate initialization technique. To obtain the HDR and disparity map from the optimal capture sequence, we propose a pipeline which alternates between estimating the camera ICRFs and the scene disparity map. We demonstrate that our optimal capture sequence leads to better results than other possible capture sequences. Our results are also close to those obtained by capturing the full stereo stack spanning the entire dynamic range. Finally, we present for the first time a stereo HDR dataset consisting of dense ISO and exposure stack captured from a smartphone dual camera. The dataset consists of 6 scenes, with an average of 142 exposure-ISO image sequence per scene.


title: Rotational Flows Generated by Microrobots Rotating Near Surfaces at Low Reynolds Number

abstract: In this study, we use numerical simulations to investigate the flow field induced by a single magnetic microrobot rotating with a constant angular speed about an axis perpendicular to an underlying surface. A parallel solver for steady Stokes flow equations based on the boundary-element method is used for simulating these flows. A simple transformation is introduced to extend the predictive capability of the solver to cases with small unsteadiness. Flows induced by four simple robot shapes are investigated: sphere, upright cylinder, horizontally-laid cylinder, and five-pointed star-shaped prism. Shapes with cross-sections that are axisymmetric about the rotation axis (sphere and upright cylinder) generate time-invariant flow fields, which could be useful for applications such as micromanipulation. Non-axisymmetric shapes (horizontally-laid cylinder and the star-shaped prism) induce significant unsteadiness inside the flow field, which could be desirable for applications such as micromixing. Furthermore, a slender horizontally-laid cylinder generates substantially three-dimensional flows, an added benefit for micromixing applications. The presence of nearby walls such as a bottom substrate or sidewalls has a retarding effect on the induced flows, which is quantified. Finally, we present the driving torque and power-consumption of these microrobots rotating in viscous liquids. The numerical modeling platform used in this work can enable future optimal microrobot designs for a given application requirement.


title: Extent of Fermi-surface reconstruction in the high-temperature superconductor HgBa$_2$CuO$_{4+\delta}$

abstract: High magnetic fields have revealed a surprisingly small Fermi-surface in underdoped cuprates, possibly resulting from Fermi-surface reconstruction due to an order parameter that breaks translational symmetry of the crystal lattice. A crucial issue concerns the doping extent of this state and its relationship to the principal pseudogap and superconducting phases. We employ pulsed magnetic field measurements on the cuprate HgBa$_2$CuO$_{4+\delta}$ to identify signatures of Fermi surface reconstruction from a sign change of the Hall effect and a peak in the temperature-dependent planar resistivity. We trace the termination of Fermi-surface reconstruction to two hole concentrations where the superconducting upper critical fields are found to be enhanced. One of these points is associated with the pseudogap end-point near optimal doping. These results connect the Fermi-surface reconstruction to both superconductivity and the pseudogap phenomena.


title: Long time behavior of the fractional Korteweg-de Vries equation with cubic nonlinearity

abstract: We prove global existence and modified scattering for the solutions of the Cauchy problem to the fractional Korteweg-de Vries equation with cubic nonlinearity for small, smooth and localized initial data.


title: A hint on the metal-free star formation rate density from 21cm-EDGES data

abstract: We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dot{\rho}_{*}^{III}$ required at high-redshifts ($z \gtrsim 16$) in order to reproduce both the amplitude and the redshift of the EDGES 21-cm global signal. Our model accounts for the Lyman Alpha (Ly$\alpha$), radio and X-ray backgrounds from both Pop III and metal-enriched Population II (Pop II) stars. For the latter, we use the star formation rate density estimates (and the Ly$\alpha$ background) from the {\it Delphi} semi-analytic model that has been shown to reproduce all key observables for galaxies at $z \gtrsim 5$; the radio and X-ray backgrounds are fixed using low-$z$ values. The constraints on the free parameters characterizing the properties of the Pop III stars are obtained using a Markov Chain Monte Carlo analysis. Our results yield a $\dot{\rho}_{*}^{III}$ that whilst increasing from $z \sim 21-16$ thereafter shows a sharp decline which is in excellent agreement with the results found by \citet{valiante2016} to simulate the growth of $z \sim 6 - 7$ quasars and their host galaxies, suggesting that the bulk of Pop III star formation occurs in the rarest and most massive metal-poor halos at $z \lesssim 20$. This allows Pop III stars to produce a rapidly growing Ly$\alpha$ background between $z \sim 21-15$. Further, Pop III stars are required to provide a radio background that is about $3-4$ orders of magnitude higher than that provided by Pop II stars although Pop II stars dominate the X-ray background.


title: Computing Robustly Forward Invariant Sets for Mixed-Monotone Systems

abstract: This work presents new tools for studying reachability and set invariance for continuous-time mixed-monotone dynamical systems subject to a disturbance input. The vector field of a mixed-monotone system is decomposable via a decomposition function into increasing and decreasing components, and this decomposition enables embedding the original dynamics in a higher-dimensional embedding system. While the original system is subject to an unknown disturbance input, the embedding system has no disturbances and its trajectories provide bounds for finite-time reachable sets of the original dynamics. Our main contribution is to show how one can efficiently identify robustly forward invariant and attractive sets for mixed-monotone systems by studying certain equilibria of this embedding system. We show also how this approach, when applied to the backward-time dynamics, establishes different robustly forward invariant sets for the original dynamics. Lastly, we present an independent result for computing decomposition functions for systems with polynomial dynamics. These tools and results are demonstrated through several examples and a case study.


title: Escaping Cannibalization? Correlation-Robust Pricing for a Unit-Demand Buyer

abstract: We consider a robust version of the revenue maximization problem, where a single seller wishes to sell $n$ items to a single unit-demand buyer. In this robust version, the seller knows the buyer's marginal value distribution for each item separately, but not the joint distribution, and prices the items to maximize revenue in the worst case over all compatible correlation structures. We devise a computationally efficient (polynomial in the support size of the marginals) algorithm that computes the worst-case joint distribution for any choice of item prices. And yet, in sharp contrast to the additive buyer case (Carroll, 2017), we show that it is NP-hard to approximate the optimal choice of prices to within any factor better than $n^{1/2-\epsilon}$. For the special case of marginal distributions that satisfy the monotone hazard rate property, we show how to guarantee a constant fraction of the optimal worst-case revenue using item pricing; this pricing equates revenue across all possible correlations and can be computed efficiently.


title: Path integral contour deformations for noisy observables

abstract: Monte Carlo studies of many quantum systems face exponentially severe signal-to-noise problems. We show that noise arising from complex phase fluctuations of observables can be reduced without introducing bias using path integral contour deformation techniques. A numerical study of contour deformations for correlation functions in Abelian gauge theory and complex scalar field theory demonstrates that variance can be reduced by orders of magnitude without modifying Monte Carlo sampling.


title: Protection of genomic information: a classical and a quantum approach

abstract: Splitting a literal genomic sequence into 4 binary files is enough to ensure confidentiality and integrity during storage and transfer of information. The binary files are resources for RSA or one-time-pad (OTP) cryptography protocols. It is speculated that representing nucleic bases as Bell states in a quantum view of a sequence would provide tools for genomic data protection when implemented in an authentic quantum computer, soon to come as a practical and readily available device.


title: Livelayer: A Semi-Automatic Software Program for Segmentation of Layers and Diabetic Macular Edema in Optical Coherence Tomography Images

abstract: Given the capacity of Optical Coherence Tomography (OCT) imaging to display symptoms of a wide variety of eye diseases and neurological disorders, the need for OCT image segmentation and the corresponding data interpretation is latterly felt more than ever before. In this paper, we wish to address this need by designing a semi-automatic software program for applying reliable segmentation of 8 different macular layers as well as outlining retinal pathologies such as diabetic macular edema. The software accommodates a novel graph-based semi-automatic method, called ''Livelayer'' which is designed for straightforward segmentation of retinal layers and fluids. This method is chiefly based on Dijkstra's Shortest Path (SPF) algorithm and the Live-wire function together with some preprocessing operations on the to-be-segmented images. The software is indeed suitable for obtaining detailed segmentation of layers, exact localization of clear or unclear fluid objects and the ground truth, demanding far less endeavor in comparison to a common manual segmentation method. It is also valuable as a tool for calculating the irregularity index in deformed OCT images. The amount of time (seconds) that Livelayer required for segmentation of ILM, IPL-INL, OPL-ONL was much less than that for the manual segmentation, 5s for the ILM (minimum) and 15.57s for the OPL-ONL (maximum). The unsigned errors (pixels) between the semi-automatically labeled and gold standard data was on average 2.7, 1.9, 2.1 for ILM, IPL-INL, OPL-ONL, respectively. The Bland-Altman plots indicated perfect concordance between the Livelayer and the manual algorithm and that they could be used interchangeably. The repeatability error was around one pixel for the OPL-ONL and < 1 for the other two. The dice scores for comparing the two algorithms and for obtaining the repeatability on segmentation of fluid objects were at acceptable levels.


title: Needmining: Identifying micro blog data containing customer needs

abstract: The design of new products and services starts with the identification of needs of potential customers or users. Many existing methods like observations, surveys, and experiments draw upon specific efforts to elicit unsatisfied needs from individuals. At the same time, a huge amount of user-generated content in micro blogs is freely accessible at no cost. While this information is already analyzed to monitor sentiments towards existing offerings, it has not yet been tapped for the elicitation of needs. In this paper, we lay an important foundation for this endeavor: we propose a Machine Learning approach to identify those posts that do express needs. Our evaluation of tweets in the e-mobility domain demonstrates that the small share of relevant tweets can be identified with remarkable precision or recall results. Applied to huge data sets, the developed method should enable scalable need elicitation support for innovation managers - across thousands of users, and thus augment the service design tool set available to him.


title: Nonlinear friction in underdamped anharmonic stochastic oscillators

abstract: Stationary states of overdamped anharmonic stochastic oscillators driven by L\'evy noise are typically multimodal. The very same situation is recorded for an underdamped L\'evy noise driven motion in single-well potentials with linear friction. Within current manuscript we relax the assumption that the friction experienced by a particle is linear. Using computer simulations, we study underdamped motion in single-well potentials in the regime of nonlinear friction. We demonstrate that it is relatively easy to observe multimodality in the velocity distribution as it is determined by the friction itself and it is the same as the multimodality in the overdamped case with the analogous deterministic force. Contrary to the velocity marginal density, it is more difficult to produce multimodality in the position. Nevertheless, for fine-tuned nonlinear friction, the spatial multimodality can be recorded.


title: Damping of elementary excitations in one-dimensional dipolar Bose gases

abstract: In the presence of dipolar interactions the excitation spectrum of a Bose gas can acquire a local minimum. The corresponding quasiparticles are known as rotons. They are gaped and do not decay at zero temperature. Here we study the decay of rotons in one-dimensional Bose gases at low temperatures. It predominantly occurs due to the backscattering of thermal phonons on rotons. The resulting rate scales with the third power of temperature and is inversely proportional to the sixth power of the roton gap near the solidification phase transition. The hydrodynamic approach used here enables us to find the decay rate for quasiparticles at practically any momenta, with minimal assumptions on the exact form of the interparticle interactions. Our results are an essential prerequisite for the description of all the dissipative phenomena in dipolar gases and have direct experimental relevance.


title: High brightness CW electron beams from Superconducting RF photoemission gun

abstract: CW photoinjectors operating at high accelerating gradients promise to revolutionize many areas of science and applications. They can establish the basis for a new generation of monochromatic X-ray free electron lasers, high brightness hadron beams, or a new generation of microchip production. In this letter we report on the record-performing superconducting RF electron gun with $\textrm{CsK}_{2}\textrm{Sb}$ photocathode. The gun is generating high charge electron bunches (up to 10 nC/bunch) and low transverse emittances, while operating for months with a single photocathode. This achievement opens a new era in generating high-power beams with a very high average brightness.


title: On a class of elliptic free boundary problems with multiple solutions

abstract: We prove that a certain class of elliptic free boundary problems, which includes the Prandtl-Batchelor problem from fluid dynamics as a special case, has two distinct nontrivial solutions for large values of a parameter. The first solution is a global minimizer of the energy. The energy functional is nondifferentiable, so standard variational arguments cannot be used directly to obtain a second nontrivial solution. We obtain our second solution as the limit of mountain pass points of a sequence of $C^1$-functionals approximating the energy. We use careful estimates of the corresponding energy levels to show that this limit is neither trivial nor a minimizer.


title: The Changing Nature of Computational Science Software

abstract: How should software engineering be adapted for Computational Science (CS)? If we understood that, then we could better support software sustainability, verifiability, reproducibility, comprehension, and usability for CS community. For example, improving the maintainability of the CS code could lead to: (a) faster adaptation of scientific project simulations to new and efficient hardware (multi-core and heterogeneous systems); (b) better support for larger teams to co-ordinate (through integration with interdisciplinary teams); and (c) an extended capability to model complex phenomena. In order to better understand computational science, this paper uses quantitative evidence (from 59 CS projects in Github) to check 13 published beliefs about CS. These beliefs reflect on (a) the nature of scientific challenges; (b) the implications of limitations of computer hardware; and (c) the cultural environment of scientific software development. What we found was, using this new data from Github, only a minority of those older beliefs can be endorsed. More than half of the pre-existing beliefs are dubious, which leads us to conclude that the nature of CS software development is changing. Further, going forward, this has implications for (1) what kinds of tools we would propose to better support computational science and (2) research directions for both communities.


title: Efficient spin torques in antiferromagnetic CoO/Pt quantified by comparing field- and current- induced switching

abstract: We achieve current-induced switching in collinear insulating antiferromagnetic CoO/Pt, with fourfold in-plane magnetic anisotropy. This is measured electrically by spin Hall magnetoresistance and confirmed by the magnetic field-induced spin-flop transition of the CoO layer. By applying current pulses and magnetic fields, we quantify the efficiency of the acting current-induced torques and estimate a current-field equivalence ratio of $4x10^{-11} T A^{-1} m^2$. The N\'eel vector final state ($n \perp j$) is in line with a thermomagnetoelastic switching mechanism for a negative magnetoelastic constant of the CoO.


title: Social confinement and mesoscopic localization of epidemics on networks

abstract: Recommendations around epidemics tend to focus on individual behaviors, with much less efforts attempting to guide event cancellations and other collective behaviors since most models lack the higher-order structure necessary to describe large gatherings. Through a higher-order description of contagions on networks, we model the impact of a blanket cancellation of events larger than a critical size and find that epidemics can suddenly collapse when interventions operate over groups of individuals rather than at the level of individuals. We relate this phenomenon to the onset of mesoscopic localization, where contagions concentrate around dominant groups.


title: Realising a species-selective double well with multiple-radiofrequency-dressed potentials

abstract: Techniques to manipulate the individual constituents of an ultracold mixture are key to investigating impurity physics. In this work, we confine a mixture of the hyperfine ground states of Rb-87 in a double-well potential. The potential is produced by dressing the atoms with multiple radiofrequencies. The amplitude and phase of each frequency component of the dressing field are individually controlled to independently manipulate each species. Furthermore, we verify that our mixture of hyperfine states is collisionally stable, with no observable inelastic loss.


title: Learning distributed representations of graphs with Geo2DR

abstract: We present Geo2DR (Geometric to Distributed Representations), a GPU ready Python library for unsupervised learning on graph-structured data using discrete substructure patterns and neural language models. It contains efficient implementations of popular graph decomposition algorithms and neural language models in PyTorch which can be combined to learn representations of graphs using the distributive hypothesis. Furthermore, Geo2DR comes with general data processing and loading methods to bring substantial speed-up in the training of the neural language models. Through this we provide a modular set of tools and methods to quickly construct systems capable of learning distributed representations of graphs. This is useful for replication of existing methods, modification, or development of completely new methods. This paper serves to present the Geo2DR library and perform a comprehensive comparative analysis of existing methods re-implemented using Geo2DR across widely used graph classification benchmarks. Geo2DR displays a high reproducibility of results in published methods and interoperability with other libraries useful for distributive language modelling.


title: $N$-body simulations for parametrised modified gravity

abstract: We present $\texttt{MG-evolution}$, an $N$-body code simulating the cosmological structure formation for parametrised modifications of gravity. It is built from the combination of parametrised linear theory with a parametrisation of the deeply nonlinear cosmological regime extrapolated from modified spherical collapse computations that cover the range of known screening mechanisms. We test $\texttt {MG-evolution}$, which runs at the speed of conventional $\Lambda$CDM simulations, against a suit of existing exact model-specific codes, encompassing linearised and chameleon $f(R)$ gravity as well as the normal branch of the Dvali-Gabadadz-Porrati braneworld model, hence covering both large-field value and large-derivative screening effects. We compare the nonlinear power spectra produced by the parametrised and model-specific approaches over the full range of scales set by the box size and resolution of our simulations, $k=(0.05-2.5)$~h/Mpc, and for two redshift slices, $z=0$ and $z=1$. We find sub-percent to one-percent level recovery of all the power spectra generated with the model-specific codes for the full range of scales. $\texttt {MG-evolution}$ can be used for generalised and accurate tests of gravity and dark energy with the increasing wealth of high-precision cosmological survey data becoming available over the next decade.


title: On the Convergence of the Dynamic Inner PCA Algorithm

abstract: Dynamic inner principal component analysis (DiPCA) is a powerful method for the analysis of time-dependent multivariate data. DiPCA extracts dynamic latent variables that capture the most dominant temporal trends by solving a large-scale, dense, and nonconvex nonlinear program (NLP). A scalable decomposition algorithm has been recently proposed in the literature to solve these challenging NLPs. The decomposition algorithm performs well in practice but its convergence properties are not well understood. In this work, we show that this algorithm is a specialized variant of a coordinate maximization algorithm. This observation allows us to explain why the decomposition algorithm might work (or not) in practice and can guide improvements. We compare the performance of the decomposition strategies with that of the off-the-shelf solver Ipopt. The results show that decomposition is more scalable and, surprisingly, delivers higher quality solutions.


title: An agent-based model of multi-dimensional opinion dynamics and opinion alignment

abstract: It is known that individual opinions on different policy issues often align to a dominant ideological dimension (e.g. "left" vs. "right") and become increasingly polarized. We provide an agent-based model that reproduces these two stylized facts as emergent properties of an opinion dynamics in a multi-dimensional space of continuous opinions. The mechanisms for the change of agents' opinions in this multi-dimensional space are derived from cognitive dissonance theory and structural balance theory. We test assumptions from proximity voting and from directional voting regarding their ability to reproduce the expected emerging properties. We further study how the emotional involvement of agents, i.e. their individual resistance to change opinions, impacts the dynamics. We identify two regimes for the global and the individual alignment of opinions. If the affective involvement is high and shows a large variance across agents, this fosters the emergence of a dominant ideological dimension. Agents align their opinions along this dimension in opposite directions, i.e. create a state of polarization.


title: A categorical framework for glider representations

abstract: Fragment and glider representations (introduced by F. Caenepeel, S. Nawal, and F. Van Oystaeyen) form a generalization of filtered modules over a filtered ring. Given a $\Gamma$-filtered ring $FR$ and a subset $\Lambda \subseteq \Gamma$, we provide a category $\operatorname{Glid}_\Lambda FR$ of glider representations, and show that it is a complete and cocomplete deflation quasi-abelian category. We discuss its derived category, and its subcategories of natural gliders and Noetherian gliders. If $R$ is a bialgebra over a field $k$ and $FR$ is a filtration by bialgebras, we show that $\operatorname{Glid}_\Lambda FR$ is a monoidal category which is derived equivalent to the category of representations of a semi-Hopf category (in the sense of E. Batista, S. Caenepeel, and J. Vercruysse). We show that the monoidal category of glider representations associated to the one-step filtration $k \cdot 1 \subseteq R$ of a bialgebra $R$ is sufficient to recover the bialgebra $R$ by recovering the usual fiber functor from $\operatorname{Glid}_\Lambda FR.$ When applied to a group algebra $kG$, this shows that the monoidal category $\operatorname{Glid}_\Lambda F(kG)$ alone is sufficient to distinguish even isocategorical groups.


title: Effects of Gravitational Chern-Simons during Axion-SU(2) Inflation

abstract: In this paper, we examine the viability of inflation models with a spectator axion field coupled to both gravitational and SU(2) gauge fields via Chern-Simons couplings. Requiring phenomenological success of the axion-SU(2) sector constrains the coupling strength of the gravitational Chern-Simons term. We find that the impact of this term on the production and propagation of gravitational waves can be as large as fifty percent enhancement for the helicity that is not sourced by the gauge field, if the cut-off scale is as low as $\Lambda$ = 20H. The effect becomes smaller for a larger value of $\Lambda$, while the impact on the helicity sourced by the gauge field is negligible regardless of $\Lambda$.


title: Three Short Period Jupiters from TESS

abstract: We report the confirmation and mass determination of three hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HIP 65Ab (TOI-129, TIC-201248411) is an ultra-short-period Jupiter orbiting a bright (V=11.1 mag) K4-dwarf every 0.98 days. It is a massive 3.213 +/- 0.078 Mjup planet in a grazing transit configuration with an impact parameter of b = 1.17 +0.10/-0.08. As a result the radius is poorly constrained, 2.03 +0.61/-0.49 Rjup. The planet's distance to its host star is less than twice the separation at which it would be destroyed by Roche lobe overflow. It is expected to spiral into HIP 65A on a timescale ranging from 80 Myr to a few gigayears, assuming a reduced tidal dissipation quality factor of Qs' = 10^7 - 10^9. We performed a full phase-curve analysis of the TESS data and detected both illumination- and ellipsoidal variations as well as Doppler boosting. HIP 65A is part of a binary stellar system, with HIP 65B separated by 269 AU (3.95 arcsec on sky). TOI-157b (TIC 140691463) is a typical hot Jupiter with a mass of 1.18 +/- 0.13 Mjup and a radius of 1.29 +/- 0.02 Rjup. It has a period of 2.08 days, which corresponds to a separation of just 0.03 AU. This makes TOI-157 an interesting system, as the host star is an evolved G9 sub-giant star (V=12.7). TOI-169b (TIC 183120439) is a bloated Jupiter orbiting a V=12.4 G-type star. It has a mass of 0.79 +/- 0.06 Mjup and a radius of 1.09 +0.08/-0.05 Rjup. Despite having the longest orbital period (P = 2.26 days) of the three planets, TOI-169b receives the most irradiation and is situated on the edge of the Neptune desert. All three host stars are metal rich with [Fe/H] ranging from 0.18 - 0.24.


title: Polarization as a tuning parameter for Floquet engineering: magnetism in the honeycomb, square, and triangular Mott insulators

abstract: Magnetic exchange couplings can be tuned by coupling to periodic light, where the frequency and amplitude are typically varied: a process known as Floquet engineering. The polarization of the light is also important, and in this paper, we show how different polarizations, including several types of unpolarized light, can tune the exchange couplings in distinct ways. Using unpolarized light, for example, it is possible to tune the material without breaking either time-reversal or any lattice symmetries. To illustrate these effects generically, we consider single-band Hubbard models at half-filling on the honeycomb, square and triangular lattices. We derive the effective Heisenberg spin models to fourth order in perturbation theory for arbitrary fixed polarizations, and several types of unpolarized light that preserve time-reversal and lattice symmetries. Coupling these models to periodic light tunes first, second and third neighbor exchange couplings, as well as ring exchange terms on the square and triangular lattices. Circularly polarized light induces chiral fields for the honeycomb and triangular lattices, which favors non-coplanar magnetism and potential chiral spin liquids. We discuss how to maximize the enhancement of the couplings without inducing heating.


title: Unified simulations of planetary formation and atmospheric evolution: Effects of pebble accretion, giant impacts, and stellar irradiations on super-Earth formation

abstract: A substantial number of super-Earths have been discovered, and atmospheres of transiting super-Earths have also been observed by transmission spectroscopy. Several lines of observational evidence indicate that most super-Earths do not possess massive H$_2$/He atmospheres. However, accretion and retention of less massive atmospheres on super-Earths challenge planet formation theory. We consider the following three mechanisms: (i) envelope heating by pebble accretion, (ii) mass loss during giant impacts, and (iii) atmospheric loss by stellar X-ray and EUV photoevaporation. We investigate whether these mechanisms influence the amount of the atmospheres that form around super-Earths. We develop a code combining an N-body simulation of pebble-driven planetary formation and an atmospheric evolution simulation. We demonstrate that the observed orbital properties of super-Earths are well reproduced by the results of our simulations. However, (i) heating by pebble accretion ceases prior to disk dispersal, (ii) the frequency of giant impact events is too low to sculpt massive atmospheres, and (iii) many super-Earths having H$_2$/He atmospheres of $\gtrsim 10$ wt% survive against stellar irradiations for 1 Gyr. Therefore, it is likely that other mechanisms such as suppression of gas accretion are required to explain less massive atmospheres ($\lesssim 10$ wt%) of super-Earths.


title: Fertility Monotonicity and Average Complexity of the Stack-Sorting Map

abstract: Let $\mathcal D_n$ denote the average number of iterations of West's stack-sorting map $s$ that are needed to sort a permutation in $S_n$ into the identity permutation $123\cdots n$. We prove that \[0.62433\approx\lambda\leq\liminf_{n\to\infty}\frac{\mathcal D_n}{n}\leq\limsup_{n\to\infty}\frac{\mathcal D_n}{n}\leq \frac{3}{5}(7-8\log 2)\approx 0.87289,\] where $\lambda$ is the Golomb-Dickman constant. Our lower bound improves upon West's lower bound of $0.23$, and our upper bound is the first improvement upon the trivial upper bound of $1$. We then show that fertilities of permutations increase monotonically upon iterations of $s$. More precisely, we prove that $|s^{-1}(\sigma)|\leq|s^{-1}(s(\sigma))|$ for all $\sigma\in S_n$, where equality holds if and only if $\sigma=123\cdots n$. This is the first theorem that manifests a law-of-diminishing-returns philosophy for the stack-sorting map that B\'ona has proposed. Along the way, we note some connections between the stack-sorting map and the right and left weak orders on $S_n$.


title: Singlet night in Feynman-ville: one-loop matching of a real scalar

abstract: A complete one-loop matching calculation for real singlet scalar extensions of the Standard Model to the Standard Model effective field theory (SMEFT) of dimension-six operators is presented. We compare our analytic results obtained by using Feynman diagrams to the expressions derived in the literature by a combination of the universal one-loop effective action (UOLEA) approach and Feynman calculus. After identifying contributions that have been overlooked in the existing calculations, we find that the pure diagrammatic approach and the mixed method lead to identical results. We highlight some of the subtleties involved in computing one-loop matching corrections in SMEFT.


title: Suppressing dissipation in a Floquet-Hubbard system

abstract: The concept of `Floquet engineering' relies on an external periodic drive to realise novel, effectively static Hamiltonians. This technique is being explored in experimental platforms across physics, including ultracold atoms, laser-driven electron systems, nuclear magnetic resonance, and trapped ions. The key challenge in Floquet engineering is to avoid the uncontrolled absorption of photons from the drive, especially in interacting systems in which the excitation spectrum becomes effectively dense. The resulting dissipative coupling to higher-lying modes, such as the excited bands of an optical lattice, has been explored in recent experimental and theoretical works, but the demonstration of a broadly applicable method to mitigate this effect is lacking. Here, we show how two-path quantum interference, applied to strongly-correlated fermions in a driven optical lattice, suppresses dissipative coupling to higher bands and increases the lifetime of double occupancies and spin-correlations by two to three orders of magnitude. Interference is achieved by introducing a weak second modulation at twice the fundamental driving frequency with a definite relative phase. This technique is shown to suppress dissipation in both weakly and strongly interacting regimes of a driven Hubbard system, opening an avenue to realising low-temperature phases of matter in interacting Floquet systems.


title: Geodesic Distance Field-based Curved Layer Volume Decomposition for Multi-Axis Support-free Printing

abstract: This paper presents a new curved layer volume decomposition method for multi-axis support-free printing of freeform solid parts. Given a solid model to be printed that is represented as a tetrahedral mesh, we first establish a geodesic distance field embedded on the mesh, whose value at any vertex is the geodesic distance to the base of the model. Next, the model is naturally decomposed into curved layers by interpolating a number of iso-geodesic distance surfaces (IGDSs). These IGDSs morph from bottom-up in an intrinsic and smooth way owing to the nature of geodesics, which will be used as the curved printing layers that are friendly to multi-axis printing. In addition, to cater to the collision-free requirement and to improve the printing efficiency, we also propose a printing sequence optimization algorithm for determining the printing order of the IGDSs, which helps reduce the air-move path length. Ample experiments in both computer simulation and physical printing are performed, and the experimental results confirm the advantages of our method.


title: Some new results about a conjecture by Brian Alspach

abstract: In this paper we consider the following conjecture, proposed by Brian Alspach, concerning partial sums in finite cyclic groups: given a subset $A$ of $\mathbb{Z}_n\setminus \{0\}$ of size $k$ such that $\sum_{z\in A} z\not= 0$, it is possible to find an ordering $(a_1,\ldots,a_k)$ of the elements of $A$ such that the partial sums $s_i=\sum_{j=1}^i a_j$, $i=1,\ldots,k$, are nonzero and pairwise distinct. This conjecture is known to be true for subsets of size $k\leq 11$ in cyclic groups of prime order. Here, we extend such result to any torsion-free abelian group and, as a consequence, we provide an asymptotic result in $\mathbb{Z}_n$. We also consider a related conjecture, originally proposed by Ronald Graham: given a subset $A$ of $\mathbb{Z}_p\setminus\{0\}$, where $p$ is a prime, there exists an ordering of the elements of $A$ such that the partial sums are all distinct. Working with the methods developed by Hicks, Ollis and Schmitt, based on the Alon's combinatorial Nullstellensatz, we prove the validity of such conjecture for subsets $A$ of size $12$.


title: Statistical deprojection of intervelocities, interdistances and masses in the Isolated Galaxy Pair Catalog

abstract: In order to study the internal dynamics of actual galaxy pairs, we need to derive the probability distribution function (PDF) of true 3D (orbital) intervelocities and interdistances between pair members from their observed (projected) values, and of the pair masses from Kepler's third law. Our Isolated Galaxy Pair Catalog (IGPC) of 13114 pairs \cite{Nottale2018a} is used here for this research. The algorithms of statistical deprojection elaborated in \cite{Nottale2018b} are applied to these observational data. We derive the orbital velocity PDFs for the whole catalog and for several selected subsamples. The interdistance PDF is deprojected and compared to analytical profiles which are expected from semi-theoretical arguments. The PDF of deprojected pair orbital velocities is characterized by the existence of a main probability peak around $\approx 150$ km.s$^{-1}$ for all subsamples of the IGPC as well as for the UGC pair catalog \cite{Chamaraux2016}. The interdistance PDFs of both the projected and deprojected data are described at large distances by the same power law with exponent $\approx -2$. The whole distributions, including their cores, are fairly fitted by King profiles. The mass deprojection yields a mass/luminosity ratio for the pairs of $M/L=(30 \pm 5)$ in Solar units. The orbital velocity probability peak is observed at the same value, $\approx 150$ km/s, as the main exoplanet velocity peak, which points toward a possible universality of Keplerian structures, whatever the scale. The pair $M/L$ ratio is just 5 times the standard ratio for luminous matter, which does not require the existence of non-baryonic dark matter in these systems.


title: Data Driven Flavour Model

abstract: A bottom-up approach has been adopted to identify a flavour model that agrees with present experimental measurements. The charged fermion mass hierarchies suggest that only the top Yukawa term should be present at the renormalisable level. Similarly, describing the lightness of the active neutrinos through the type-I Seesaw mechanism, right-handed neutrino mass terms should also be present at the renormalisable level. The flavour symmetry of the Lagrangian including the fermionic kinetic terms and only the top Yukawa is then a combination of $U(2)$ and $U(3)$ factors. Once considering the Majorana neutrino terms, the associated symmetry is $O(3)$. Lighter charged fermion and active neutrino masses and quark and lepton mixings arise considering specific spurion fields {\it \`a la} Minimal Flavour Violation. The associated phenomenology is investigated and the model turns out to have almost the same flavour protection as the Minimal Flavour Violation in both quark and lepton sectors. Promoting the spurions to dynamical fields, the associated scalar potential is also studied and a minimum is identified such that fermion masses and mixings are correctly reproduced. Very precise predictions for the Majorana phases follow from the minimisation of the scalar potential and thus the neutrinoless-double-beta decay may represent a smoking gun for the model.


title: An empirical "high-confidence" candidate zone for $Fermi$ BL Lacertae objects

abstract: In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.


title: Strong biases in retrieved atmospheric composition caused by strong day-night chemical heterogeneities

abstract: Most planets currently amenable to transit spectroscopy are close enough to their host star to exhibit a relatively strong day to night temperature gradient. For hot planets, this leads to cause a chemical composition dichotomy between the two hemispheres. In the extreme case of ultra hot jupiters, some species, such as molecular hydrogen and water, are strongly dissociated on the day-side while others, such as carbon monoxide, are not. However, most current retrieval algorithm rely on 1D forward models that are unable to model this effect. We thus investigate how the 3D structure of the atmosphere biases the abundances retrieved using commonly used algorithms. We study the case of Wasp-121b as a prototypical ultra hot Jupiter. We use the simulations of this planet performed with the SPARC/MIT global climate model (GCM) and generate transmission spectra that fully account for the 3D structure of the atmosphere with Pytmopsh3R. These spectra are then analyzed using the \taurex retrieval code. We find that such ultra hot jupiter's transmission spectra exhibit muted H$_2$O features that originate in the night-side where the temperature, hence the scale-height, is smaller than on the day-side. However, the spectral features of molecules present on the day-side are boosted by both its high temperature and low mean molecular weight. As a result, the retrieved parameters are strongly biased compared to the ground truth. In particular the [CO]/[H$_2$O] is overestimated by one to three orders of magnitude. This must be kept in mind when using such retrieval analysis to infer the C/O ratio of a planet's atmosphere. We also discuss whether indicators can allow us to infer the 3D structure of an observed atmosphere. Finally we show that HST/WFC3 transmission data of Wasp-121b are compatible with the day-night thermal and compositional dichotomy predicted by models.


title: L-GALAXIES 2020: Spatially resolved cold gas phases, star formation and chemical enrichment in galactic discs

abstract: We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of discs is stored in 12 concentric rings. The new model retains the success of its predecessor in reproducing the observed evolution of the galaxy population, in particular, stellar mass functions and passive fractions over the redshift range 0<=z<=3 and mass range 8<=log(M_*/Msun)<=12, the black hole-bulge mass relation at z=0, galaxy morphology as a function of stellar mass and the mass-metallicity relations of both stellar and gas components. In addition, its detailed modelling of the radial structure of discs allows qualitatively new comparisons with observation, most notably with the relative sizes and masses of the stellar, atomic and molecular components in discs. Good agreement is found with recent data. Comparison of results obtained for simulations differing in mass resolution by more than two orders of magnitude shows that all important distributions are numerically well converged even for this more detailed model. An examination of metallicity and surface-density gradients in the stars and gas indicates that our new model, with star formation, chemical enrichment and feedback calculated self-consistently on local disc scales, reproduces some but not all of the trends seen in recent many-galaxy IFU surveys.


title: No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

abstract: We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.


title: Constraining the Source of the High-velocity Ejecta in Type Ia SN 2019ein

abstract: We present multiwavelength photometric and spectroscopic observations of SN 2019ein, a high-velocity Type Ia supernova (SN Ia) discovered in the nearby galaxy NGC 5353 with a two-day nondetection limit. SN 2019ein exhibited some of the highest measured expansion velocities of any SN Ia, with a Si II absorption minimum blueshifted by 24,000 km s$^{-1}$ at 14 days before peak brightness. More unusually, we observed the emission components of the P Cygni profiles to be blueshifted upward of 10,000 km s$^{-1}$ before B-band maximum light. This blueshift, among the highest in a sample of 28 other Type Ia supernovae, is greatest at our earliest spectroscopic epoch and subsequently decreases toward maximum light. We discuss possible progenitor systems and explosion mechanisms that could explain these extreme absorption and emission velocities. Radio observations beginning 14 days before B-band maximum light yield nondetections at the position of SN 2019ein, which rules out symbiotic progenitor systems, most models of fast optically thick accretion winds, and optically thin shells of mass $\lesssim 10^{-6}$ M$_\odot$ at radii $< 100$ AU. Comparing our spectra to models and observations of other high-velocity SNe Ia, we find that SN 2019ein is well fit by a delayed-detonation explosion. We propose that the high emission velocities may be the result of abundance enhancements due to ejecta mixing in an asymmetric explosion, or optical depth effects in the photosphere of the ejecta at early times. These findings may provide evidence for common explosion mechanisms and ejecta geometries among high-velocity SNe Ia.


title: Inflationary Twin Models

abstract: We introduce the concept of inflationary twin models, that is, a class of generalized (kinetic or "K-inflation") field theories which lead exactly to the same cosmological evolution during the inflationary period as a given standard scalar field theory of inflation. The twin concept permits to introduce generalized, K-inflation theories in a controlled manner, maintaining some inflationary predictions unaltered. Further, this concept allows to extend analytical tools like the slow-roll expansion to the realm of K-inflation theories, facilitating their investigation. Twins of a standard scalar field model of inflation still may result in different results for some inflationary observables, because of their nontrivial scalar speed of sound. This implies that non-standard twins may lead to completely viable models of inflation, even if their standard twin is ruled out by observations. We provide some explicit examples of this possibility, including the physical case of a dilaton--Dirac-Born-Infeld theory.


title: Topology of superconductors beyond mean-field theory

abstract: The study of topological superconductivity is largely based on the analysis of mean-field Hamiltonians that violate particle number conservation and have only short-range interactions. Although this approach has been very successful, it is not clear that it captures the topological properties of real superconductors, which are described by number-conserving Hamiltonians with long-range interactions. To address this issue, we study topological superconductivity directly in the number-conserving setting. We focus on a diagnostic for topological superconductivity that compares the fermion parity $\mathcal{P}$ of the ground state of a system in a ring geometry and in the presence of zero vs. $\Phi_{\text{sc}}=\frac{h}{2e} \equiv \pi$ flux of an external magnetic field. A version of this diagnostic exists in any dimension and provides a $\mathbb{Z}_2$ invariant $\nu=\mathcal{P}_0\mathcal{P}_{\pi}$ for topological superconductivity. In this paper we prove that the mean-field approximation correctly predicts the value of $\nu$ for a large family of number-conserving models of spinless superconductors. Our result applies directly to the cases of greatest physical interest, including $p$-wave and $p_x+ip_y$ superconductors in one and two dimensions, and gives strong evidence for the validity of the mean-field approximation in the study of (at least some aspects of) topological superconductivity.

categories: cond-mat.supr-con hep-th math-ph math.MP


title: Collisionless kinetic theory for parametrically pumped magnons

abstract: We discuss collisionless kinetic equations describing the non-equilibrium dynamics of magnons in a ferromagnet exposed to an oscillating microwave field. Previously, this problem has been treated within the so-called "S-theory" where the collision integral in the kinetic equation for the magnon distribution is either neglected or taken into account phenomenologically via an effective relaxation time. However, the possibility of magnon condensation has not been included in S-theory. Moreover, the momentum integrations appearing in the magnon self-energies are usually decoupled by retaining only the term where the loop momentum is equal to the external momentum. In this work we critically examine the accuracy of these approximations and develop the proper extensions of S-theory. We show that these extensions can significantly modify the time evolution of the magnon distribution.


title: A Link Between Ram Pressure Stripping and Active Galactic Nuclei

abstract: The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within "jellyfish galaxies," galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. We find unambiguous morphological evidence for ram pressure stripping. For lower mass galaxies (with stellar masses < 10^9.5 solar masses) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be the found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which then produces heating and outflows due to AGN feedback. AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.


title: Dynamical mass inference of galaxy clusters with neural flows

abstract: We present an algorithm for inferring the dynamical mass of galaxy clusters directly from their respective phase-space distributions, i.e. the observed line-of-sight velocities and projected distances of galaxies from the cluster centre. Our method employs normalizing flows, a deep neural network capable of learning arbitrary high-dimensional probability distributions, and inherently accounts, to an adequate extent, for the presence of interloper galaxies which are not bounded to a given cluster, the primary contaminant of dynamical mass measurements. We validate and showcase the performance of our neural flow approach to robustly infer the dynamical mass of clusters from a realistic mock cluster catalogue. A key aspect of our novel algorithm is that it yields the probability density function of the mass of a particular cluster, thereby providing a principled way of quantifying uncertainties, in contrast to conventional machine learning approaches. The neural network mass predictions, when applied to a contaminated catalogue with interlopers, have a mean overall logarithmic residual scatter of 0.028 dex, with a log-normal scatter of 0.126 dex, which goes down to 0.089 dex for clusters in the intermediate to high mass range. This is an improvement by nearly a factor of four relative to the classical cluster mass scaling relation with the velocity dispersion, and outperforms recently proposed machine learning approaches. We also apply our neural flow mass estimator to a compilation of galaxy observations of some well-studied clusters with robust dynamical mass estimates, further substantiating the efficacy of our algorithm.


title: Leakage reduction in fast superconducting qubit gates via optimal control

abstract: Reaching high speed, high fidelity qubit operations requires precise control over the shape of the underlying pulses. For weakly anharmonic systems, such as superconducting transmon qubits, short gates lead to leakage to states outside of the computational subspace. Control pulses designed with open-loop optimal control may reduce such leakage. However, model inaccuracies can severely limit the usability of such pulses. We implemented a closed-loop optimization that simultaneously adapts all control parameters based on measurements of a cost function built from Clifford gates. By parameterizing pulses with a piecewise-constant representation that matches the capabilities of the control hardware we create a $4.16~\rm{ns}$ single-qubit pulse with $99.76\,\%$ fidelity and $0.044\,\%$ leakage. This is a seven-fold reduction of the leakage rate of the best DRAG pulse we have calibrated at such short durations on the same system.


title: The case for jets in cataclysmic variables

abstract: For decades cataclysmic variables (CVs) were thought to be one of the few classes of accreting compact objects to not launch jets, and have consequently been used to constrain jet launching models. However, recent theoretical and observational advances indicate that CVs do in fact launch jets. Specifically, it was demonstrated that their accretion-outflow cycle is analogous to that of their higher mass cousins -- the X-ray Binaries (XRBs). Subsequent observations of the CV SS Cygni confirmed this and have consistently shown radio flaring equivalent to that in the XRBs that marks a transient jet. Based on this finding and the emission properties, several studies have concluded that the radio emission is most likely from a transient jet. Observations of other CVs, while not conclusive, are consistent with this interpretation. However, the issue is not yet settled. Later observations have raised a number of questions about this model, as well as about potential alternative radio emission mechanisms. CVs are non-relativistic and many have well-determined distances; these properties would make them ideal candidates with which to address many of our outstanding questions about fundamental jet physics. Here we review the case for jets in CVs, discuss the outstanding questions and issues, and outline the future work necessary to conclusively answer the question of whether CVs launch jets.


title: Emergent QCD$_3$ Quantum Phase Transitions of Fractional Chern Insulators

abstract: Motivated by the recent work of QED$_3$-Chern-Simons quantum critical points of fractional Chern insulators (Phys. Rev. X \textbf{8}, 031015, (2018)), we study its non-Abelian generalizations, namely QCD$_3$-Chern-Simons quantum phase transitions of fractional Chern insulators. These phase transitions are described by Dirac fermions interacting with non-Abelian Chern-Simons gauge fields ($U(N)$, $SU(N)$, $USp(N)$, etc.). Utilizing the level-rank duality of Chern-Simons gauge theory and non-Abelian parton constructions, we discuss two types of QCD$_3$ quantum phase transitions. The first type happens between two Abelian states in different Jain sequences, as opposed to the QED3 transitions between Abelian states in the same Jain sequence. A good example is the transition between $\sigma^{xy}=1/3$ state and $\sigma^{xy}=-1$ state, which has $N_f=2$ Dirac fermions interacting with a $U(2)$ Chern-Simons gauge field. The second type is naturally involving non-Abelian states. For the sake of experimental feasibility, we focus on transitions of Pfaffian-like states, including the Moore-Read Pfaffian, anti-Pfaffian, particle-hole Pfaffian, etc. These quantum phase transitions could be realized in experimental systems such as fractional Chern insulators in graphene heterostructures.


title: Post-Estimation Smoothing: A Simple Baseline for Learning with Side Information

abstract: Observational data are often accompanied by natural structural indices, such as time stamps or geographic locations, which are meaningful to prediction tasks but are often discarded. We leverage semantically meaningful indexing data while ensuring robustness to potentially uninformative or misleading indices. We propose a post-estimation smoothing operator as a fast and effective method for incorporating structural index data into prediction. Because the smoothing step is separate from the original predictor, it applies to a broad class of machine learning tasks, with no need to retrain models. Our theoretical analysis details simple conditions under which post-estimation smoothing will improve accuracy over that of the original predictor. Our experiments on large scale spatial and temporal datasets highlight the speed and accuracy of post-estimation smoothing in practice. Together, these results illuminate a novel way to consider and incorporate the natural structure of index variables in machine learning.


title: Reconfiguration of Colourings and Dominating Sets in Graphs: a Survey

abstract: We survey results concerning reconfigurations of colourings and dominating sets in graphs. The vertices of the $k$-colouring graph $\mathcal{C}_{k}(G)$ of a graph $G$ correspond to the proper $k$-colourings of a graph $G$, with two $k$-colourings being adjacent whenever they differ in the colour of exactly one vertex. Similarly, the vertices of the $k$-edge-colouring graph $\mathcal{EC}_{k}(G)$ of $g$ are the proper $k$-edge-colourings of $G$, where two $k$-edge-colourings are adjacent if one can be obtained from the other by switching two colours along an edge-Kempe chain, i.e., a maximal two-coloured alternating path or cycle of edges. The vertices of the $k$-dominating graph $\mathcal{D}_{k}(G)$ are the (not necessarily minimal) dominating sets of $G$ of cardinality $k$ or less, two dominating sets being adjacent in $\mathcal{D}_{k}(G)$ if one can be obtained from the other by adding or deleting one vertex. On the other hand, when we restrict the dominating sets to be minimum dominating sets, for example, we obtain different types of domination reconfiguration graphs, depending on whether vertices are exchanged along edges or not. We consider these and related types of colouring and domination reconfiguration graphs. Conjectures, questions and open problems are stated within the relevant sections.


title: Integrable Matrix Models in Discrete Space-Time

abstract: We introduce a class of integrable dynamical systems of interacting classical matrix-valued fields propagating on a discrete space-time lattice, realized as many-body circuits built from elementary symplectic two-body maps. The models provide an efficient integrable Trotterization of non-relativistic $\sigma$-models with complex Grassmannian manifolds as target spaces, including, as special cases, the higher-rank analogues of the Landau-Lifshitz field theory on complex projective spaces. As an application, we study transport of Noether charges in canonical local equilibrium states. We find a clear signature of superdiffusive behavior in the Kardar-Parisi-Zhang universality class, irrespectively of the chosen underlying global unitary symmetry group and the quotient structure of the compact phase space, providing a strong indication of superuniversal physics.

categories: cond-mat.stat-mech hep-th math-ph math.MP


title: Optimal market making with persistent order flow

abstract: \noindent We address the issue of market making on electronic markets when taking into account the clustering and long memory properties of market order flows. We consider a market model with one market maker and order flows driven by general Hawkes processes. We formulate the market maker's objective as a stochastic control problem. We characterize an optimal control by proving existence and uniqueness of a viscosity solution to the associated Hamilton-Jacobi-Bellman equation. Finally we propose a fully consistent numerical method allowing to implement this optimal strategy in practice.


title: Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

abstract: In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.


title: On the Bloch-Kato conjecture for GSp(4)

abstract: We prove an explicit reciprocity law for the Euler system attached to the spin motive of a genus 2 Siegel modular form. As consequences, we obtain one inclusion of the Iwasawa Main Conjecture for such motives, and the Bloch--Kato conjecture in analytic rank 0 for their critical twists.


title: W2S: Microscopy Data with Joint Denoising and Super-Resolution for Widefield to SIM Mapping

abstract: In fluorescence microscopy live-cell imaging, there is a critical trade-off between the signal-to-noise ratio and spatial resolution on one side, and the integrity of the biological sample on the other side. To obtain clean high-resolution (HR) images, one can either use microscopy techniques, such as structured-illumination microscopy (SIM), or apply denoising and super-resolution (SR) algorithms. However, the former option requires multiple shots that can damage the samples, and although efficient deep learning based algorithms exist for the latter option, no benchmark exists to evaluate these algorithms on the joint denoising and SR (JDSR) tasks. To study JDSR on microscopy data, we propose such a novel JDSR dataset, Widefield2SIM (W2S), acquired using a conventional fluorescence widefield and SIM imaging. W2S includes 144,000 real fluorescence microscopy images, resulting in a total of 360 sets of images. A set is comprised of noisy low-resolution (LR) widefield images with different noise levels, a noise-free LR image, and a corresponding high-quality HR SIM image. W2S allows us to benchmark the combinations of 6 denoising methods and 6 SR methods. We show that state-of-the-art SR networks perform very poorly on noisy inputs. Our evaluation also reveals that applying the best denoiser in terms of reconstruction error followed by the best SR method does not necessarily yield the best final result. Both quantitative and qualitative results show that SR networks are sensitive to noise and the sequential application of denoising and SR algorithms is sub-optimal. Lastly, we demonstrate that SR networks retrained end-to-end for JDSR outperform any combination of state-of-the-art deep denoising and SR networks


title: Tube-based Robust Model Predictive Control for a Distributed Parameter System Modeled as a Polytopic LPV (extended version)

abstract: Distributed parameter systems (DPS) are formulated as partial differential equations (PDE). Especially, under time-varying boundary conditions, PDE introduce force coupling. In the case of the flexible stacker crane (STC), nonlinear coupling is introduced. Accordingly, online trajectory planning and tracking can be addressed using a nonlinear model predictive control (NMPC). However, due to the high computational demands of a NMPC, this paper discusses a possibility of embedding nonlinearities inside a linear parameter varying (LPV) system and thus make a use of a numerically low-demanding linear MPC. The resulting mismatches are treated as parametric and additive uncertainties in the context of robust tube-based MPC (TMPC). For the proposed approach, most of the computations are carried out offline. Only a simple convex quadratic program (QP) is conducted online. Additionally a soft-constrained extension was briefly proposed. Simulation results are used to illustrate the good performance, closed-loop stability and recursive feasibility of the proposed approach despite uncertainties.


title: Multichannel asymmetric transmission through a dimer defect with saturable inter-site nonlinearity

abstract: We consider the asymmetric transmission properties of a Discrete Nonlinear Schr{\"o}dinger type dimer with a saturable nonlinear intersite coupling between the dimer sites, in addition to a cubic onsite nonlinearity and asymmetric linear onsite potentials. In contrast to previously studied cases with pure onsite nonlinearities, the transmission coefficient for stationary transmission is shown to be a multivalued function of the transmitted intensity, in regimes of low saturability and small or moderate transmitted intensity. The corresponding backward transfer map is analyzed analytically and numerically, and shown to have either one or three distinct solutions for saturable coupling, and zero or two solutions for the purely cubic nonlinear coupling. As saturation strength is increased, the multi-solution regimes disappear through bifurcations and the single-valued regime of the onsite model is recovered. The existence of multiple solution branches yields novel mechanisms for asymmetric left/right stationary transmission: in addition to shifts of the positions of transmission peaks, peaks for transmission in one direction may correspond to nonexistence of stationary solutions propagating in the opposite direction, at the corresponding branch and transmitted intensity. Moreover, one of these transmission channels behaves as a nearly perfect mirror for incoming signals. The linear stability of the stationary solutions is analyzed, and instabilities are typically observed, and illustrated by direct numerical simulations, in regimes with large transmission coefficient. Intersite nonlinearities are found to prevent the formation of a localized dimer mode in the instability-induced dynamics. Finally, the partial reflection and transmission of a Gaussian excitation is studied, and the asymmetric transmission properties are compared to previously studied onsite models.


title: Internally heated porous convection: an idealised model for Enceladus' hydrothermal activity

abstract: Recent planetary data and geophysical modelling suggest that hydrothermal activity is ongoing under the ice crust of Enceladus, one of Saturn's moons. According to these models, hydrothermal flow in the porous, rocky core of the satellite is driven by tidal deformation that induces dissipation and volumetric internal heating. Despite the effort in the modelling of Enceladus' interior, systematic understanding---and even basic scaling laws---of internally-heated porous convection and hydrothermal activity are still lacking. In this article, using an idealised model of an internally-heated porous medium, we explore numerically and theoretically the flows that develop close and far from the onset of convection. In particular, we quantify heat-transport efficiency by convective flows as well as the typical extent and intensity of heat-flux anomalies created at the top of the porous layer. With our idealised model, we derive simple and general laws governing the temperature and hydrothermal velocity that can be driven in the oceans of icy moons. In the future, these laws could help better constraining models of the interior of Enceladus and other icy satellites.


title: Observation of strong two-electron--one-photon transitions in few-electron ion

abstract: We resonantly excite the $K$ series of O$^{5+}$ and O$^{6+}$ up to principal quantum number $n=11$ with monochromatic x rays, producing $K$-shell holes, and observe their relaxation by soft-x-ray emission. Some photoabsorption resonances of O$^{5+}$ reveal strong two-electron--one-photon (TEOP) transitions. We find that for the $[(1s\,2s)_1\,5p_{3/2}]_{3/2;1/2}$ states, TEOP relaxation is by far stronger than the radiative decay and competes with the usually much faster Auger decay path. This enhanced TEOP decay arises from a strong correlation with the near-degenerate upper states $[(1s\,2p_{3/2})_1\,4s]_{3/2;1/2}$ of a Li-like satellite blend of the He-like $K\alpha$ transition. Even in three-electron systems, TEOP transitions can play a dominant role, and the present results should guide further research on the ubiquitous and abundant many-electron ions where electronic energy degeneracies are far more common and configuration mixing is stronger.


title: Hypersensitive tunable Josephson escape sensor for gigahertz astronomy

abstract: Sensitive photon detection in the gigahertz band constitutes the cornerstone to study different phenomena in astronomy, such as radio burst sources, galaxy formation, cosmic microwave background, axions, comets, gigahertz-peaked spectrum radio sources and supermassive black holes. Nowadays, state of the art detectors for astrophysics are mainly based on transition edge sensors and kinetic inductance detectors. Overall, most sensible nanobolometers so far are superconducting detectors showing a noise equivalent power (NEP) as low as 2x10-20 W/Hz1/2. Yet, fast thermometry at the nanoscale was demonstrated as well with Josephson junctions through switching current measurements. In general, detection performance are set by the fabrication process and limited by used materials. Here, we conceive and demonstrate an innovative tunable Josephson escape sensor (JES) based on the precise current control of the temperature dependence of a fully superconducting one-dimensional nanowire Josephson junction. The JES might be at the core of future hypersensitive in situ-tunable bolometers or single-photon detectors working in the gigahertz regime. Operated as a bolometer the JES points to a thermal fluctuation noise (TFN) NEP_TFN 1x10-25 W/Hz1/2, which as a calorimeter bounds the frequency resolution above 2 GHz, and resolving power below 40 at 50 GHz, as deduced from the experimental data. Beyond the obvious applications in advanced ground-based and space telescopes for gigahertz astronomy, the JES might represent a breakthrough in several fields of quantum technologies ranging from subTHz communications and quantum computing to cryptography and quantum key distribution.


title: Rank two free groups and integer points on real cubic surfaces

abstract: Counting integer points on the Markoff cubic is closely related to questions in hyperbolic geometry. In a previous work with Igor Rivin we investigated the regularity of the geodesic length function for a punctured torus. Here we extend this work to the three holed sphere and related orbifolds.


title: Statistical Inference for High Dimensional Panel Functional Time Series

abstract: In this paper we develop statistical inference tools for high dimensional functional time series. We introduce a new concept of physical dependent processes in the space of square integrable functions, which adopts the idea of basis decomposition of functional data in these spaces, and derive Gaussian and multiplier bootstrap approximations for sums of high dimensional functional time series. These results have numerous important statistical consequences. Exemplarily, we consider the development of joint simultaneous confidence bands for the mean functions and the construction of tests for the hypotheses that the mean functions in the spatial dimension are parallel. The results are illustrated by means of a small simulation study and in the analysis of Canadian temperature data.


title: Fibre inflation and precision CMB data

abstract: Generic features of models of inflation obtained from string compactifications are the correlations between the model parameters and the postinflationary evolution of the universe. Thus, the postinflationary evolution depends on the inflationary model parameters and accurate inflationary predictions require that this be incorporated in the evolution of the primordial spectrum. The fibre inflation model is a promising model of inflation constructed in type IIB string theory. This model has two interesting features in its postinflationary evolution. The reheating temperature of the model is directly correlated with the model parameters. The model also necessarily predicts some dark radiation, which can be sizable for certain choices of discrete parameters in the model. We analyze this model in detail using publicly available codes - ModeChord and CosmoMC with the latest Planck+BICEP2/Keck array data to constrain the model parameters and $N_{\rm pivot}$ (the number of $e$-foldings between horizon exit of the CMB pivot mode and the end of inflation). We also carry out the same analysis using the publicly available code Cobaya. We find the results of both the analysis to be in agreement. Our analysis provides the basic methods necessary to extract precise inflationary prediction in string models incorporating correlations between model parameters and postinflationary evolution.


title: LiDAR guided Small obstacle Segmentation

abstract: Detecting small obstacles on the road is critical for autonomous driving. In this paper, we present a method to reliably detect such obstacles through a multi-modal framework of sparse LiDAR(VLP-16) and Monocular vision. LiDAR is employed to provide additional context in the form of confidence maps to monocular segmentation networks. We show significant performance gains when the context is fed as an additional input to monocular semantic segmentation frameworks. We further present a new semantic segmentation dataset to the community, comprising of over 3000 image frames with corresponding LiDAR observations. The images come with pixel-wise annotations of three classes off-road, road, and small obstacle. We stress that precise calibration between LiDAR and camera is crucial for this task and thus propose a novel Hausdorff distance based calibration refinement method over extrinsic parameters. As a first benchmark over this dataset, we report our results with 73% instance detection up to a distance of 50 meters on challenging scenarios. Qualitatively by showcasing accurate segmentation of obstacles less than 15 cms at 50m depth and quantitatively through favourable comparisons vis a vis prior art, we vindicate the method's efficacy. Our project-page and Dataset is hosted at https://small-obstacle-dataset.github.io/


title: Semi-adsorption-controlled growth window for half Heusler FeVSb epitaxial films

abstract: The electronic, magnetic, thermoelectric, and topological properties of Heusler compounds (composition $XYZ$ or $X_2 YZ$) are highly sensitive to stoichiometry and defects. Here we establish the existence and experimentally map the bounds of a \textit{semi} adsorption-controlled growth window for semiconducting half Heusler FeVSb films, grown by molecular beam epitaxy (MBE). We show that due to the high volatility of Sb, the Sb stoichiometry is self-limiting for a finite range of growth temperatures and Sb fluxes, similar to the growth of III-V semiconductors such as GaSb and GaAs. Films grown within this window are nearly structurally indistinguishable by X-ray diffraction (XRD) and reflection high energy electron diffraction (RHEED). The highest electron mobility and lowest background carrier density are obtained towards the Sb-rich bound of the window, suggesting that Sb-vacancies may be a common defect. Similar \textit{semi} adsorption-controlled bounds are expected for other ternary intermetallics that contain a volatile species $Z=$\{Sb, As, Bi\}, e.g., CoTiSb, LuPtSb, GdPtBi, and NiMnSb. However, outstanding challenges remain in controlling the remaining Fe/V ($X/Y$) transition metal stoichiometry.


title: Skyrmion Dynamics and Transverse Mobility: Skyrmion Hall Angle Reversal on 2D Periodic Substrates with dc and Biharmonic ac Drives

abstract: We numerically examine the dynamics of a skyrmion interacting with a two-dimensional periodic substrate under dc and biharmonic ac drives. We show that the Magnus force of the skyrmion produces circular orbits that can resonate with the ac drive and the periodicity of the substrate to create quantized motion both parallel and perpendicular to the dc drive. The skyrmion Hall angle exhibits a series of increasing and/or decreasing steps along with strongly fluctuating regimes. In the phase locked regimes, the skyrmion Hall angle is constant and the skyrmion motion consists of periodic orbits encircling an integer number of obstacles per every or every other ac drive cycle. We also observe phases in which the skyrmion moves at $90^\circ$ with respect to the driving direction even in the presence of damping, a phenomenon called absolute transverse mobility that can exhibit reentrance as a function of dc drive. When the biharmonic ac drives have different amplitudes, in the two directions we find regimes in which the skyrmion Hall angle shows a sign reversal from positive to negative, as well as a reentrant pinning effect in which the skyrmion is mobile at low drives but becomes pinned at higher drives. These behaviors arise due to the combination of the Magnus force with the periodic motion of the skyrmions, which produce Shapiro steps, directional locking, and ratchet effects.


title: AskCI Server: Collaborative and version controlled knowledge base

abstract: AskCI server is a collaborative, open source documentation server that uses GitHub for automation and version control of shared knowledge. A programmatic application programming interface, friendly user interface, and organization of concepts into questions makes it versatile as a support or collaborative knowledge base.


title: Vortices of Electro-osmotic Flow in Heterogeneous Porous Media

abstract: Traditional models of electrokinetic transport in porous media are based on homogenized material properties, which neglect any macroscopic effects of microscopic fluctuations. This perspective is taken not only for convenience, but also motivated by the expectation of irrotational electro-osmotic flow, proportional to the electric field, for uniformly charged surfaces (or constant zeta potential) in the limit of thin double layers. Here, we show that the inherent heterogeneity of porous media generally leads to macroscopic vortex patterns, which have important implications for convective transport and mixing. These vortical flows originate due to competition between pressure-driven and electro-osmotic flows, and their size are characterized by the correlation length of heterogeneity in permeability or surface charge. The appearance of vortices is controlled by a single dimensionless control parameter, defined as the ratio of a typical electro-osmotic velocity to the total mean velocity.


title: A sharp inequality for the variance with respect to the Ewens Sampling Formula

abstract: The variance of a linear statistic defined on the symmetric group endowed with the Ewens probability is examined. Despite the dependence of the summands, it can be bounded from above by a constant multiple of the sum of variances. We find the exact value of this constant. The analysis of the appearing quadratic forms and eigenvalue search is built upon the exponential matrices and discrete Hahn's polynomials.


title: LHC at 10: the physics legacy

abstract: Ten years into the LHC physics programme, I reflect on the vast amount of extraordinary results that have emerged so far.


title: Existence and uniqueness theorems for pointwise slant immersions in complex space forms

abstract: An isometric immersion $f: M^{n} \rightarrow \tilde M^{m}$ from an $n$-dimensional Riemannian manifold $M^{n}$ into an almost Hermitian manifold $\tilde M^{m}$ of complex dimension $m$ is called pointwise slant if its Wirtinger angles define a function defined on $M$. In this paper we establish the existence and uniqueness theorems for pointwise slant immersions of Riemannian manifolds $M^{n}$ into a complex space form $\tilde M^{n}(c)$ of constant holomorphic sectional curvature $c$.


title: Brighter-fatter effect in near-infrared detectors -- III. Fourier-domain treatment of flat field correlations and application to WFIRST

abstract: Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. WFIRST will fly a focal plane containing 18 Teledyne H4RG-10 near infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work has shown that correlation functions of flat field images are effective tools for disentangling linear and non-linear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical non-linearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the non-linear IPC. The magnitude of the BFE is such that the effective area of a pixel is increased by $(3.54\pm0.03)\times 10^{-7}$ for every electron deposited in a neighboring pixel. We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113%. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors.


title: Power and Sample Size for Marginal Structural Models

abstract: Marginal structural models fit via inverse probability of treatment weighting are commonly used to control for confounding when estimating causal effects from observational data. When planning a study that will be analyzed with marginal structural modeling, determining the required sample size for a given level of statistical power is challenging because of the effect of weighting on the variance of the estimated causal means. This paper considers the utility of the design effect to quantify the effect of weighting on the precision of causal estimates. The design effect is defined as the ratio of the variance of the causal mean estimator divided by the variance of a naive estimator if, counter to fact, no confounding had been present and weights were not needed. A simple, closed-form approximation of the design effect is derived that is outcome invariant and can be estimated during the study design phase. Once the design effect is approximated for each treatment group, sample size calculations are conducted as for a randomized trial, but with variances inflated by the design effects to account for weighting. Simulations demonstrate the accuracy of the design effect approximation, and practical considerations are discussed.