title: On the origin of GW190425

abstract: The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable "case BB" mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.


title: Renormalization to localization without a small parameter

abstract: We study the wave function localization properties in a d-dimensional model of randomly spaced particles with isotropic hopping potential depending solely on Euclidean interparticle distances. Due to the generality of this model usually called the Euclidean random matrix model, it arises naturally in various physical contexts such as studies of vibrational modes, artificial atomic systems, liquids and glasses, ultracold gases and photon localization phenomena. We generalize the known Burin-Levitov renormalization group approach, formulate universal conditions sufficient for localization in such models and inspect a striking equivalence of the wave function spatial decay between Euclidean random matrices and translation-invariant long-range lattice models with a diagonal disorder.


title: Turbulence in stratified atmospheres: implications for the intracluster medium

abstract: The gas motions in the intracluster medium (ICM) are governed by stratified turbulence. Stratified turbulence is fundamentally different from Kolmogorov (isotropic, homogeneous) turbulence; kinetic energy not only cascades from large to small scales, but it is also converted into buoyancy potential energy. To understand the density and velocity fluctuations in the ICM, we conduct high-resolution ($1024^2\times 1536$ grid points) hydrodynamical simulations of subsonic turbulence (with rms Mach number $\mathcal{M}\approx 0.25$) and different levels of stratification, quantified by the Richardson number $\mathrm{Ri}$, from $\mathrm{Ri}=0$ (no stratification) to $\mathrm{Ri}=13$ (strong stratification). We quantify the density, pressure and velocity fields for varying stratification because observational studies often use surface brightness fluctuations to infer the turbulent gas velocities of the ICM. We find that the standard deviation of the logarithmic density fluctuations ($\sigma_s$), where $s=\ln(\rho/\left<\rho(z)\right>)$, increases with $\mathrm{Ri}$. For weakly stratified subsonic turbulence ($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new $\sigma_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation, $\sigma_s^2=\ln(1+b^2\mathcal{M}^4+0.09\mathcal{M}^2\mathrm{Ri}H_P/H_S)$, where $b=1/3$--$1$ is the turbulence driving parameter, and $H_P$ and $H_S$ are the pressure and entropy scale heights respectively. We further find that the power spectrum of density fluctuations, $P(\rho_k/\left<\rho\right>)$, increases in magnitude with increasing $\mathrm{Ri}$, whereas the velocity power spectrum is invariant. Thus, the ratio between density and velocity power spectra strongly depends on $\mathrm{Ri}$. Pressure fluctuations, on the other hand, are independent of stratification and only depend on $\mathcal{M}$.

categories: astro-ph.GA astro-ph.CO physics.ao-ph physics.flu-dyn


title: Geometry -- dependent effects in Majorana nanowires

abstract: Starting from the Bogolubov - de Gennes theory describing the induced p-wave superconductivity in the Majorana wire of an arbitrary shape, we predict a number of intriguing phenomena such as the geometry-dependent phase battery (or a phi-Josephson junction with the spontaneous superconducting phase difference) and generation of additional quasiparticle modes at the Fermi level with the spatial position tuned by the external magnetic field direction. This tuning can be used to extend the capabilities of the braiding protocols in Majorana networks.


title: Quantum Paracrystalline Shear Modes of the Electron Liquid

abstract: Unlike classical fluids, a quantum Fermi liquid can support a long-lived and propagating shear sound wave at arbitrarily small wave vectors and frequencies, reminiscent of the transverse sound in crystals, despite lacking any form of long-range crystalline order. This mode is expected to be present in moderately interacting metals where the quasiparticle mass is renormalized to be more than twice the bare mass in two dimensions (2D), but it has remained undetected because it is hard to excite since it does not involve charge density fluctuations, in contrast to the conventional plasma mode. In this work we propose a strategy to excite and detect this unconventional mode in clean metallic channels. We show that the shear sound is responsible for the appearance of sharp dips in the ac conductance of narrow channels at resonant frequencies matching its dispersion. The liquid resonates while minimizing its dissipation in an analogous fashion to a sliding crystal. Ultra-clean 2D materials that can be tuned towards the Wigner crystallization transition such as silicon metal-oxide-semiconductor field-effect transistors, MgZnO/ZnO, p-GaAs, and AlAs quantum wells are promising platforms to experimentally discover the shear sound.


title: Cosmic recombination history in light of EDGES measurements of the cosmic dawn 21-cm signal

abstract: The recent EDGES measurements of the global 21-cm signal from the cosmic dawn suggest that the kinetic temperature of the inter-galactic medium (IGM) might be significantly lower compared to its expected value. The colder IGM directly affects the hydrogen recombination of the universe during the cosmic dawn and dark ages by enhancing the rate of recombinations. Here, we study and quantify, the impact of the colder IGM scenario on the recombination history of the universe in the context of DM-baryonic interaction model which is widely used to explain the depth of the EDGES 21-cm signal. We find that, in general, the hydrogen ionisation fraction gets suppressed during the dark ages and cosmic dawn and the suppression gradually increases at lower redshifts until X-ray heating turns on. However, accurate estimation of the ionisation fraction requires knowledge of the entire thermal history of the IGM, from the epoch of thermal decoupling of hydrogen gas and the CMBR to the cosmic dawn. It is possible that two separate scenarios which predict very similar HI differential temperature during the cosmic dawn and are consistent with the EDGES 21-cm signal might have very different IGM temperature during the dark ages. The evolutions of the ionisation fraction in these two scenarios are quite different. This prohibits us to accurately calculate the ionisation fraction during the cosmic dawn using the EDGES 21-cm signal alone. We find that the changes in the ionisation fraction w.r.t the standard scenario at redshift $z \sim 17 $ could be anything between $\sim 0 \%$ to $\sim 36 \%$. This uncertainty may be reduced if measurements of HI 21-cm differential temperature at multiple redshifts are simultaneously used.


title: Long-Period High-Amplitude Red Variables in the KELT Survey

abstract: We present a sample of 4,132 Mira-like variables (red variables with long periods and high amplitudes) in the Kilodegree Extremely Little Telescope (KELT) survey. Of these, 814 are new detections. We used 2MASS colors to identify candidate asymptotic giant branch (AGB) stars. We tested for photometric variability among the sample and used Lomb-Scargle to determine the periodicity of the variable sample. We selected variables with high amplitudes and strong periodic behavior using a Random Forest classifier. Of the sample of 4,132 Mira-like variables, we estimate that 70% are Miras, and 30% are semi-regular (SR) variables. We also adopt the method of using (W_{RP} - W_{K_s}) vs. (J - K_s) colors (Lebzelter et al. 2018) in distinguishing between O-rich and C-rich Miras and find it to be an improvement over 2MASS colors.


title: Temporal Interlacing Network

abstract: For a long time, the vision community tries to learn the spatio-temporal representation by combining convolutional neural network together with various temporal models, such as the families of Markov chain, optical flow, RNN and temporal convolution. However, these pipelines consume enormous computing resources due to the alternately learning process for spatial and temporal information. One natural question is whether we can embed the temporal information into the spatial one so the information in the two domains can be jointly learned once-only. In this work, we answer this question by presenting a simple yet powerful operator -- temporal interlacing network (TIN). Instead of learning the temporal features, TIN fuses the two kinds of information by interlacing spatial representations from the past to the future, and vice versa. A differentiable interlacing target can be learned to control the interlacing process. In this way, a heavy temporal model is replaced by a simple interlacing operator. We theoretically prove that with a learnable interlacing target, TIN performs equivalently to the regularized temporal convolution network (r-TCN), but gains 4% more accuracy with 6x less latency on 6 challenging benchmarks. These results push the state-of-the-art performances of video understanding by a considerable margin. Not surprising, the ensemble model of the proposed TIN won the $1^{st}$ place in the ICCV19 - Multi Moments in Time challenge. Code is made available to facilitate further research at https://github.com/deepcs233/TIN


title: Exceptional Collections for Mirrors of Invertible Polynomials

abstract: We prove the existence of a full exceptional collection for the derived category of equivariant matrix factorizations of an invertible polynomial with its maximal symmetry group. This proves a conjecture of Hirano--Ouchi. In the Gorenstein case, we also prove a stronger version of this conjecture due to Takahashi. Namely, that the full exceptional collection is strong.


title: The better half -- Asymmetric star-formation due to ram pressure in the EAGLE simulations

abstract: We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading half", and the opposite one "trailing half", (ii) perpendicular to the $\rm \texttt{radial}$ position of the satellite to the centre of the cluster, (iii) that maximizes the star-formation rate ($\rm SFR$) difference between the two halves. For (i), we find an enhancement of the $\rm SFR$, star formation efficiency ($\rm SFE$), and interstellar medium pressure in the leading half with respect to the trailing one and normal star-forming galaxies in the EAGLE simulation, and a clear overabundance of gas particles in their trailing. These results suggest that ram pressure (RP) is boosting the star formation by gas compression in the leading half, and transporting the gas to the trailing half. This effect is more pronounced in satellites of intermediate stellar masses $\rm 10^{9.5-10.5} M_{\odot}$, with gas masses above $\rm 10^{9} M_{\odot}$, and located within one virial radius or in the most massive clusters. In (iii) we find an alignment between the velocity and the vector perpendicular to the plane that maximizes the $\rm SFR$ difference between the two halves. It suggests that finding this plane in real galaxies can provide insights into the velocity direction.


title: Cech cohomology, homoclinic trajectories and robustness of non-saddle sets

abstract: In this paper we study flows $\varphi:M\times\mathbb{R}\longrightarrow M$ having an isolated non-saddle set. We see that the complexity of the region of influence of an isolated non-saddle set $K$ depends on the way in which $K$ sits on the phase space at the cohomological level. We construct flows in surfaces having isolated non-saddle sets with a prescribed structure for its region of influence. We also study parametrized families of smooth flows and continuations of non-saddle sets.


title: Study of the ionization efficiency for nuclear recoils in pure crystals

abstract: We study the basic integral equation in Lindhard's theory describing the energy given to atomic motion by nuclear recoils in a pure material when the atomic binding energy is taken into account. The numerical solution, which depends only on the slope of the velocity-proportional electronic stopping power and the binding energy, leads to an estimation of the ionization efficiency which is in good agreement with the available experimental measurements for Si and Ge. In this model, the quenching factor for nuclear recoils features a cut-off at an energy equal to twice the assumed binding energy. We argue that the model is a reasonable approximation for Ge even for energies close to the cutoff, while for Si is valid up to recoil energies greater than ~500 eV.


title: Regularity of optimal sets for some functional involving eigenvalues of an operator in divergence form

abstract: In this paper we consider minimizers of the functional \begin{equation*} \min \big\{ \lambda_1(\Omega)+\cdots+\lambda_k(\Omega) + \Lambda|\Omega|, \ : \ \Omega \subset D \text{ open} \big\} \end{equation*} where $D\subset\mathbb{R}^d$ is a bounded open set and where $0<\lambda_1(\Omega)\leq\cdots\leq\lambda_k(\Omega)$ are the first $k$ eigenvalues on $\Omega$ of an operator in divergence form with Dirichlet boundary condition and with H\"{o}lder continuous coefficients. We prove that the optimal sets $\Omega^\ast$ have finite perimeter and that their free boundary $\partial\Omega^\ast\cap D$ is composed of a regular part, which is locally the graph of a $C^{1,\alpha}$-regular function, and a singular part, which is empty if $d<d^\ast$, discrete if $d=d^\ast$ and of Hausdorff dimension at most $d-d^\ast$ if $d>d^\ast$, for some $d^\ast\in\{5,6,7\}$.


title: Urban Street Network Analysis in a Computational Notebook

abstract: Computational notebooks offer researchers, practitioners, students, and educators the ability to interactively conduct analytics and disseminate reproducible workflows that weave together code, visuals, and narratives. This article explores the potential of computational notebooks in urban analytics and planning, demonstrating their utility through a case study of OSMnx and its tutorials repository. OSMnx is a Python package for working with OpenStreetMap data and modeling, analyzing, and visualizing street networks anywhere in the world. Its official demos and tutorials are distributed as open-source Jupyter notebooks on GitHub. This article showcases this resource by documenting the repository and demonstrating OSMnx interactively through a synoptic tutorial adapted from the repository. It illustrates how to download urban data and model street networks for various study sites, compute network indicators, visualize street centrality, calculate routes, and work with other spatial data such as building footprints and points of interest. Computational notebooks help introduce methods to new users and help researchers reach broader audiences interested in learning from, adapting, and remixing their work. Due to their utility and versatility, the ongoing adoption of computational notebooks in urban planning, analytics, and related geocomputation disciplines should continue into the future.


title: On the embeddability of $[3]*K$

abstract: We relate the embeddability of the simplicial complex $[3]*K$ into $\mathbb{R}^{n+2}$ to that of $K$ into $\mathbb{R}^n$. In brief, the embeddability of $K$ into $\mathbb{R}^n$, in the metastable range $2n\geq 3(d+1)$, is equivalent to the embeddability of $[3]*K$ into $\mathbb{R}^{n+2}$. We show moreover than the van Kampen obstruction of $K$ vanishes if and only if the van Kampen obstruction of $[3]*K$ vanishes. It follows that for $d=2$, embeddabilty of $[3]*K$ is equivalent with the vanishing of the van Kampen obstruction for $K$, but not with the embeddability of $K$.


title: Robust Gaussian Joint Source-Channel Coding Under the Near-Zero Bandwidth Regime

abstract: Minimum power required to achieve a distortion-noise profile, i.e., a function indicating the maximum allowed distortion value for each noise level, is studied for the transmission of Gaussian sources over Gaussian channels under a regime of bandwidth approaching zero. A simple but instrumental lower bound to the minimum required power for a given profile is presented. For an upper bound, a dirty-paper based coding scheme is proposed and its power-distortion tradeoff is analyzed. Finally, upper and lower bounds to the minimum power is analyzed and compared for specific distortion-noise profiles, namely rational profiles with order one and two.


title: Compact groups with many elements of bounded order

abstract: L\'evai and Pyber proposed the following as a conjecture: Let $G$ be a profinite group such that the set of solutions of the equation $x^n=1$ has positive Haar measure. Then $G$ has an open subgroup $H$ and an element $t$ such that all elements of the coset $tH$ have order dividing $n$ (see Problem 14.53 of [The Kourovka Notebook, No. 19, 2019]). The validity of the conjecture has been proved in [Arch. Math. (Basel) 75 (2000) 1-7] for $n=2$. Here we study the conjecture for compact groups $G$ which are not necessarily profinite and $n=3$; we show that in the latter case the group $G$ contains an open normal $2$-Engel subgroup.


title: Trust in AutoML: Exploring Information Needs for Establishing Trust in Automated Machine Learning Systems

abstract: We explore trust in a relatively new area of data science: Automated Machine Learning (AutoML). In AutoML, AI methods are used to generate and optimize machine learning models by automatically engineering features, selecting models, and optimizing hyperparameters. In this paper, we seek to understand what kinds of information influence data scientists' trust in the models produced by AutoML? We operationalize trust as a willingness to deploy a model produced using automated methods. We report results from three studies -- qualitative interviews, a controlled experiment, and a card-sorting task -- to understand the information needs of data scientists for establishing trust in AutoML systems. We find that including transparency features in an AutoML tool increased user trust and understandability in the tool; and out of all proposed features, model performance metrics and visualizations are the most important information to data scientists when establishing their trust with an AutoML tool.


title: Variational-Correlations Approach to Quantum Many-body Problems

abstract: We investigate an approach for studying the ground state of a quantum many-body Hamiltonian that is based on treating the correlation functions as variational parameters. In this approach, the challenge set by the exponentially-large Hilbert space is circumvented by approximating the positivity of the density matrix, order-by-order, in a way that keeps track of a limited set of correlation functions. In particular, the density-matrix description is replaced by a correlation matrix whose dimension is kept linear in system size, to all orders of the approximation. Unlike the conventional variational principle which provides an upper bound on the ground-state energy, in this approach one obtains a lower bound instead. By treating several one-dimensional spin $1/2$ Hamiltonians, we demonstrate the ability of this approach to produce long-range correlations, and a ground-state energy that converges to the exact result. Possible extensions, including to higher-excited states are discussed.


title: A perturbation view of level-set methods for convex optimization

abstract: Level-set methods for convex optimization are predicated on the idea that certain problems can be parameterized so that their solutions can be recovered as the limiting process of a root-finding procedure. This idea emerges time and again across a range of algorithms for convex problems. Here we demonstrate that strong duality is a necessary condition for the level-set approach to succeed. In the absence of strong duality, the level-set method identifies $\epsilon$-infeasible points that do not converge to a feasible point as $\epsilon$ tends to zero. The level-set approach is also used as a proof technique for establishing sufficient conditions for strong duality that are different from Slater's constraint qualification.


title: Multiwavelength behaviour of the blazar 3C279: decade-long study from $\gamma$-ray to radio

abstract: We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.


title: The independence of Stone's Theorem from the Boolean Prime Ideal Theorem

abstract: We give a permutation model in which Stone's Theorem (every metric space is paracompact) is false and the Boolean Prime Ideal Theorem (every ideal in a Boolean algebra extends to a prime ideal) is true. The erring metric space in our model attains only rational distances and is not metacompact. Transfer theorems give the comparable independence in the Zermelo-Fraenkel setting, answering a question of Good, Tree and Watson.


title: QSO photometric redshifts from SDSS, WISE and GALEX colours

abstract: Machine learning techniques, specifically the k-nearest neighbour algorithm applied to optical band colours, have had some success in predicting photometric redshifts of quasi-stellar objects (QSOs): Although the mean of differences between the spectroscopic and photometric redshifts is close to zero, the distribution of these differences remains wide and distinctly non-Gaussian. As per our previous empirical estimate of photometric redshifts, we find that the predictions can be significantly improved by adding colours from other wavebands, namely the near-infrared and ultraviolet. Self-testing this, by using half of the 33 643 strong QSO sample to train the algorithm, results in a significantly narrower spread for the remaining half of the sample. Using the whole QSO sample to train the algorithm, the same set of magnitudes return a similar spread for a sample of radio sources (quasars). Although the matching coincidence is relatively low (739 of the 3663 sources having photometry in the relevant bands), this is still significantly larger than from the empirical method (2%) and thus may provide a method with which to obtain redshifts for the vast number of continuum radio sources expected to be detected with the next generation of large radio telescopes.


title: Tschirnhaus transformations after Hilbert

abstract: Let RD(n) denote the minimum d for which there exists a formula for the roots of the general degree n polynomial using only algebraic functions of d or fewer variables. In 1927, Hilbert sketched how the 27 lines on a cubic surface could be used to construct a 4-variable formula for the general degree 9 polynomial (implying $RD(9)\le 4$). In this paper, we turn Hilbert's sketch into a general method. We show this method produces best-to-date upper bounds on RD(n) for all n, improving earlier results of Hamilton, Sylvester, Segre and Brauer.


title: Higgs decay into a lepton pair and a photon revisited

abstract: We present new calculations of the differential decay rates for $H\to \ell^+\ell^- \gamma$ with $\ell=e$ or $\mu$ in the Standard Model. The branching fractions and forward-backward asymmetries, defined in terms of the flight direction of the photon relative to the lepton momenta, depend on the cuts on energies and invariant masses of the final state particles. For typical choices of these cuts we find the branching ratios $B(H\to e \bar e \gamma)=5.8\cdot 10^{-5}$ and $B(H\to \mu \bar \mu \gamma)=6.4\cdot 10^{-5}$ and the forward-backward asymmetries $\mathcal{A}^{(e)}_{\text{FB}}=0.343$ and $\mathcal{A}^{(\mu)}_{\text{FB}}=0.255$. We provide compact analytic expressions for the differential decay rates for the use in experimental analyses.


title: Population-based metaheuristics for Association Rule Text Mining

abstract: Nowadays, the majority of data on the Internet is held in an unstructured format, like websites and e-mails. The importance of analyzing these data has been growing day by day. Similar to data mining on structured data, text mining methods for handling unstructured data have also received increasing attention from the research community. The paper deals with the problem of Association Rule Text Mining. To solve the problem, the PSO-ARTM method was proposed, that consists of three steps: Text preprocessing, Association Rule Text Mining using population-based metaheuristics, and text postprocessing. The method was applied to a transaction database obtained from professional triathlon athletes' blogs and news posted on their websites. The obtained results reveal that the proposed method is suitable for Association Rule Text Mining and, therefore, offers a promising way for further development.


title: Large Field Excursions from Dimensional (De)construction

abstract: An inflation model based on dimensional (de)construction of a massive gauge theory is proposed. The inflaton in this model is the "zero-mode" of a component of the massive gauge field in the (de)constructed extra dimensions. The inflaton potential originates from the gauge invariant Stueckelberg potential. At low energy, the field range of the inflaton is enhanced by a factor $N^{\frac{d}{2}}$ compared to the field range of the original fields in the model, where $d$ is the number of the (de)constructed extra dimensions and $N$ is the number of the lattice points in each (de)constructed dimension. This enhancement of the field range is used to achieve a trans-Planckian inflaton field excursion. The extension of the mechanism ``excursions through KK modes'' to the case of (de)constructed extra dimensions is also studied. The burst of particle productions by this mechanism may have observable consequences in a region of the model parameter space.


title: On Scheduler Side-Channels in Dynamic-Priority Real-Time Systems

abstract: While the existence of scheduler side-channels has been demonstrated recently for fixed-priority real-time systems (RTS), there have been no similar explorations for dynamic-priority systems. The dynamic nature of such scheduling algorithms, e.g., EDF, poses a significant challenge in this regard. In this paper we demonstrate that side-channels exist in dynamic priority real-time systems. Using this side-channel, our proposed DyPS algorithm is able to effectively infer, with high precision, critical task information from the vantage point of an unprivileged (user space) task. Apart from demonstrating the effectiveness of DyPS, we also explore the various factors that impact such attack algorithms using a large number of synthetic task sets. We also compare against the state-of-the-art and demonstrate that our proposed DyPS algorithms outperform the ScheduLeak algorithms in attacking the EDF RTS.


title: Personality Traits and Drug Consumption. A Story Told by Data

abstract: This is a preprint version of the first book from the series: "Stories told by data". In this book a story is told about the psychological traits associated with drug consumption. The book includes: - A review of published works on the psychological profiles of drug users. - Analysis of a new original database with information on 1885 respondents and usage of 18 drugs. (Database is available online.) - An introductory description of the data mining and machine learning methods used for the analysis of this dataset. - The demonstration that the personality traits (five factor model, impulsivity, and sensation seeking), together with simple demographic data, give the possibility of predicting the risk of consumption of individual drugs with sensitivity and specificity above 70% for most drugs. - The analysis of correlations of use of different substances and the description of the groups of drugs with correlated use (correlation pleiades). - Proof of significant differences of personality profiles for users of different drugs. This is explicitly proved for benzodiazepines, ecstasy, and heroin. - Tables of personality profiles for users and non-users of 18 substances. The book is aimed at advanced undergraduates or first-year PhD students, as well as researchers and practitioners. No previous knowledge of machine learning, advanced data mining concepts or modern psychology of personality is assumed. For more detailed introduction into statistical methods we recommend several undergraduate textbooks. Familiarity with basic statistics and some experience in the use of probabilities would be helpful as well as some basic technical understanding of psychology.


title: On the comparison of nearby cycles via $b$-functions

abstract: In this article, we give a simple proof of the comparison of nearby and vanishing cycles in the sense of Riemann-Hilbert correspondence following the idea of Beilinson and Bernstein, without using the Kashiwara-Malgrange $V$-filtrations.


title: A Light Scalar Explanation of $(g-2)_{\mu}$ and the KOTO Anomaly

abstract: The KOTO experiment has recently performed a search for neutral Kaons decaying into neutral pions and a pair of neutrinos. Three events were observed in the KOTO signal region, with an expected background of about 0.05. Since no clear signal of systematic errors have been found, the excess of events in the decay $K_{L}\rightarrow\pi^0\nu\bar{\nu}$ is quite intriguing. One possibility to explain this anomaly would be the presence of a scalar $\phi$ with mass of the order of the pion mass and inducing decays $K_L \to \pi^0 \phi$ which mimic the observed signal.A scalar with mass of the order of the pion mass and a coupling to muons of the order of the Standard Model Higgs coupling could also explain the muon anomalous magnetic moment anomaly $(g-2)_{\mu}$. We built on these facts to show that a light singlet scalar with couplings to the leptons and quarks as the ones induced by mixing with Higgs states in two Higgs doublet models may lead to an explanation of both anomalies. More specifically, we show that this is the case in the so-called type-X models in which leptons and quarks couple to two different Higgs doublets, and for scalar masses that are in the range between 40 and 70 MeV. Due to the relatively large coupling to leptons required to fit $(g-2)_{\mu}$, the scalar lifetime accidentally falls into the sub-nanosecond range which is essential to evade the severe proton beam dump experiments and astrophysical constraints, though it becomes sensitive to constraints from electron beam dump experiments. The additional phenomenological properties of this model are discussed.


title: Large-area, high-NA Multi-level Diffractive Lens via inverse design

abstract: Flat lenses enable thinner, lighter, and simpler imaging systems. However, large-area and high-NA flat lenses have been elusive due to computational and fabrication challenges. Here, we applied inverse design to create a multi-level diffractive lens (MDL) with thickness <1.35{\mu}m, diameter of 4.13mm, NA=0.9 at wavelength of 850nm. Since the MDL is created in polymer, it can be cost-effectively replicated via imprint lithography.


title: Fermi surfaces of composite fermions

abstract: The fractional quantum Hall (FQH) effect was discovered in two-dimensional electron systems subject to a large perpendicular magnetic field nearly four decades ago. It helped launch the field of topological phases, and in addition, because of the quenching of the kinetic energy, gave new meaning to the phrase "correlated matter". Most FQH phases are gapped like insulators and superconductors; however, a small subset with even denominator fractional fillings nu of the Landau level, typified by nu = 1/2, are found to be gapless, with a Fermi surface akin to metals. We discuss our results, obtained numerically using the infinite Density Matrix Renormalization Group (iDMRG) scheme, on the effect of non-isotropic distortions with discrete N-fold rotational symmetry of the Fermi surface at zero magnetic field on the Fermi surface of the correlated nu = 1/2 state. We find that while the response for N = 2 (elliptical) distortions is significant (and in agreement with experimental observations with no adjustable parameters), it decreases very rapidly as N is increased. Other anomalies, like resilience to breaking the Fermi surface into disjoint pieces, are also found. This highlights the difference between Fermi surfaces formed from the kinetic energy, and those formed of purely potential energy terms in the Hamiltonian.


title: Effect of magnetic field on holographic insulator/superconductor phase transition in higher dimensional Gauss-Bonnet gravity

abstract: In this paper, we have investigated the effect of magnetic field numerically as well as analytically for holographic insulator/superconductor phase transition in higher dimensional Gauss-Bonnet gravity. First we have analysed the critical phenomena with magnetic field using two different numerical methods, namely, quasinormal modes method and the shooting method. Then we have carried out our calculation analytically using the St$\ddot{u}$rm-Liouville eigenvalue method. The methods show that marginally stable modes emerge at critical values of the chemical potential and the magnetic field satisfying the relation $\Lambda^2\equiv\mu^2-B$. We observe that the value of the chemical potential and hence the value of $\Lambda$ increases with higher values of the Gauss-Bonnet parameter and dimension of spacetime for a fixed mass of the scalar field. This clearly indicates that the phase transition from insulator to superconductor becomes difficult in the presence of the magnetic field for higher values of the Gauss-Bonnet parameter and dimension of spacetime. Our analytic results are in very good agreement with our numerical results.


title: Optimization of Radio Array Telescopes to Search for Fast RadioBursts

abstract: We present projected Fast Radio Burst detection rates from surveys carried out using a set of hypothetical close-packed array telescopes. The cost efficiency of such a survey falls at least as fast as the inverse square of the survey frequency. There is an optimum array element effective area in the range 0 to 25 $\rm{m^2}$. If the power law index of the FRB integrated source count versus fluence $\alpha = d ~ln R/d ~ln F > -1$ the most cost effective telescope layout uses individual dipole elements, which provides an all-sky field of view. If $\alpha <-1$ dish arrays are more cost effective.


title: Ground state stability in two spin glass models

abstract: An important but little-studied property of spin glasses is the stability of their ground states to changes in one or a finite number of couplings. It was shown in earlier work that, if multiple ground states are assumed to exist, then fluctuations in their energy differences --- and therefore the possibility of multiple ground states --- are closely related to the stability of their ground states. Here we examine the stability of ground states in two models, one of which is presumed to have a ground state structure that is qualitatively similar to other realistic short-range spin glasses in finite dimensions.


title: Activism by the AI Community: Analysing Recent Achievements and Future Prospects

abstract: The artificial intelligence community (AI) has recently engaged in activism in relation to their employers, other members of the community, and their governments in order to shape the societal and ethical implications of AI. It has achieved some notable successes, but prospects for further political organising and activism are uncertain. We survey activism by the AI community over the last six years; apply two analytical frameworks drawing upon the literature on epistemic communities, and worker organising and bargaining; and explore what they imply for the future prospects of the AI community. Success thus far has hinged on a coherent shared culture, and high bargaining power due to the high demand for a limited supply of AI talent. Both are crucial to the future of AI activism and worthy of sustained attention.


title: In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter

abstract: The dependence of the diffusion MRI signal on the diffusion time $t$ is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity $D(t)$ and mean kurtosis $K(t)$ in cortical gray matter and in 25 gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for $t=21.2$-100 ms, and is characterized by a power-law tail $\sim t^{-\vartheta}$ with dynamical exponent $\vartheta$. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using effective medium theory formalisms, we derive an exact relation between the power-law tails in $D(t)$ and $K(t)$. The estimated power-law dynamical exponent $\vartheta\simeq1/2$ in both $D(t)$ and $K(t)$ is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the $\vartheta=1/2$ exponent. In contrast, the K\"arger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from $K(t)$ falls below our measured $t$-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the $t^{-1/2}$ tail in $D(t)$ and $K(t)$ alltogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.


title: Multiblock variable influence on orthogonal projections (MB-VIOP) for enhanced interpretation of total, global, local and unique variations in OnPLS models

abstract: For multivariate data analysis involving only two input matrices, the previously published methods for variable influence on projection (e.g., VIPOPLS or VIPO2PLS) are widely used for variable selection purposes, including (i) variable importance assessment, (ii) dimensionality reduction of big data and (iii) interpretation enhancement of PLS, OPLS and O2PLS models. For multiblock analysis, the OnPLS models find relationships among multiple data matrices by calculating latent variables; however, a method for improving the interpretation of these latent variables by assessing the importance of the input variables was not available up to now. A method for variable selection in multiblock analysis, called multiblock variable influence on orthogonal projections (MB-VIOP) is explained in this paper. MB-VIOP is a model based variable selection method that uses the data matrices, the scores and the normalized loadings of an OnPLS model in order to sort the input variables of more than two data matrices according to their importance for both simplification and interpretation of the total multiblock model, and also of the unique, local and global model components separately. MB-VIOP has been tested using three multiblock datasets. MB-VIOP assesses the variable importance in any type of data. MB-VIOP connects the input variables of different data matrices according to their relevance for the interpretation of each latent variable, yielding enhanced interpretability for each OnPLS model component. Besides, MB-VIOP can deal with strong overlapping of types of variation, as well as with many data blocks with very different dimensionality. The ability of MB-VIOP for generating dimensionality reduced models with high interpretability makes this method ideal for big data mining, multi-omics data integration and any study that requires exploration and interpretation of large streams of data.


title: Forming Diverse Super-Earth Systems in Situ

abstract: Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multi-transiting systems like Kepler-11 to systems like our Solar System's terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super Earths, including the apparent dichotomy between single transiting and multiple transiting system. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler's transiting planets. Our simulations can account for the planets' distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for single than multi transiting planets; smaller eccentricities for larger planets; scatter in the mass-radius relation, including lower densities for planets with masses measured with TTVs than RVs; and similarity in planets' sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.


title: Fountains and storms: The role of AGN and mergers in disrupting the cool-core in the RomulusC simulation

abstract: The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem.


title: Study of excited $\Lambda_\mathrm{b}^0$ states decaying to $\Lambda_\mathrm{b}^0\pi^+\pi^-$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A study of excited $\Lambda_\mathrm{b}^0$ baryons is reported, based on a data sample collected in 2016-2018 with the CMS detector at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of up to 140 fb$^{-1}$. The existence of four excited $\Lambda_\mathrm{b}^0$ states: $\Lambda_\mathrm{b}$(5912)$^0$, $\Lambda_\mathrm{b}$(5920)$^0$, $\Lambda_\mathrm{b}$(6146)$^0$, and $\Lambda_\mathrm{b}$(6152)$^0$ in the $\Lambda_\mathrm{b}^0\pi^+\pi^-$ mass spectrum is confirmed, and their masses are measured. The $\Lambda_\mathrm{b}^0\pi^+\pi^-$ mass distribution exhibits a broad excess of events in the region of 6040-6100 MeV, whose origin cannot be discerned with the present data.


title: Strain-Controlled Magnetic Ordering in 2D Carbon Metamaterials

abstract: We use ab initio spin-polarized density functional theory to study the magnetic order in a Kagom\'e-like 2D metamaterial consisting of pristine or substitutionally doped phenalenyl radicals polymerized into a nanoporous, graphene-like structure. In this and in a larger class of related structures, the constituent polyaromatic hydrocarbon molecules can be considered as quantum dots that may carry a net magnetic moment. The structure of this porous system and the coupling between the quantum dots may be changed significantly by applying moderate strain, thus allowing to control the magnetic order and the underlying electronic structure.


title: CHAOS Challenge -- Combined (CT-MR) Healthy Abdominal Organ Segmentation

abstract: Segmentation of abdominal organs has been a comprehensive, yet unresolved, research field for many years. In the last decade, intensive developments in deep learning (DL) have introduced new state-of-the-art segmentation systems. In order to expand the knowledge on these topics, the CHAOS - Combined (CT-MR) Healthy Abdominal Organ Segmentation challenge has been organized in conjunction with IEEE International Symposium on Biomedical Imaging (ISBI), 2019, in Venice, Italy. CHAOS provides both abdominal CT and MR data from healthy subjects for single and multiple abdominal organ segmentation. Five different but complementary tasks have been designed to analyze the capabilities of current approaches from multiple perspectives. The results are investigated thoroughly, compared with manual annotations and interactive methods. The analysis shows that the performance of DL models for single modality (CT / MR) can show reliable volumetric analysis performance (DICE: 0.98 $\pm$ 0.00 / 0.95 $\pm$ 0.01) but the best MSSD performance remain limited (21.89 $\pm$ 13.94 / 20.85 $\pm$ 10.63 mm). The performances of participating models decrease significantly for cross-modality tasks for the liver (DICE: 0.88 $\pm$ 0.15 MSSD: 36.33 $\pm$ 21.97 mm) and all organs (DICE: 0.85 $\pm$ 0.21 MSSD: 33.17 $\pm$ 38.93 mm). Despite contrary examples on different applications, multi-tasking DL models designed to segment all organs seem to perform worse compared to organ-specific ones (performance drop around 5\%). Besides, such directions of further research for cross-modality segmentation would significantly support real-world clinical applications. Moreover, having more than 1500 participants, another important contribution of the paper is the analysis on shortcomings of challenge organizations such as the effects of multiple submissions and peeking phenomena.


title: Simpler Partial Derandomization of PPSZ for $k$-SAT

abstract: We give a simpler derandomization of the best known $k$-SAT algorithm PPSZ [FOCS'97, JACM'05] for $k$-SAT with \emph{sub-exponential} number of solutions. The existing derandomization uses a complicated construction of small sample space, while we only use \emph{hashing}. Our algorithm and theorem also have a nice byproduct: It outperforms the current fastest deterministic $k$-SAT algorithm when the formula has \emph{moderately exponential} number of solutions.


title: Neutral-Current Neutrino Scattering from the Deuteron

abstract: Neutral-current neutrino scattering from the deuteron leading to proton-neutron final states is considered using an approach that incorporates relativistic dynamics and consequently provides robust modeling at relatively high energies and momenta. In this work the focus is placed on the fully exclusive reaction where both the proton and neutron in the final state are assumed to be detected. Accordingly, the incident neutrino energy, the neutrino scattering angle and the scattered neutrino's energy can all be reconstructed. It is shown that for specific choices of kinematics the reaction proceeds mainly via scattering from the proton, while for other choices of kinematics it proceeds mainly from the neutron. Specific asymmetries are introduced to focus on these attributes. Measurements in both regions have the potential to yield valuable information on the nucleon's electroweak form factors at momentum transfers up to a (GeV/c)$^2$. In particular, the cross sections are shown to be very sensitive to the isoscalar axial-vector form factor, and sensitive but less so to the magnetic strangeness form factor. Comparisons with other reactions, specifically charge-changing neutrino reactions and both parity-conserving and -violating electron scattering, have the potential to provide new ways to test the Standard Model.


title: Adapting Grad-CAM for Embedding Networks

abstract: The gradient-weighted class activation mapping (Grad-CAM) method can faithfully highlight important regions in images for deep model prediction in image classification, image captioning and many other tasks. It uses the gradients in back-propagation as weights (grad-weights) to explain network decisions. However, applying Grad-CAM to embedding networks raises significant challenges because embedding networks are trained by millions of dynamically paired examples (e.g. triplets). To overcome these challenges, we propose an adaptation of the Grad-CAM method for embedding networks. First, we aggregate grad-weights from multiple training examples to improve the stability of Grad-CAM. Then, we develop an efficient weight-transfer method to explain decisions for any image without back-propagation. We extensively validate the method on the standard CUB200 dataset in which our method produces more accurate visual attention than the original Grad-CAM method. We also apply the method to a house price estimation application using images. The method produces convincing qualitative results, showcasing the practicality of our approach.


title: Multi-charged ion-water molecule collisions in a classical-trajectory time-dependent mean-field theory

abstract: A recently proposed classical-trajectory dynamical screening model for the description of multiple ionization and capture during ion-water molecule collisions is extended to incorporate dynamical screening on both the multi-center target potential and the projectile ion. Comparison with available experimental data for He$^{2+}$ + H$_2$O collisions at intermediate energies (10-150 keV/u) and Li$^{3+}$ + H$_2$O at higher energies (100-850 keV/u) demonstrates the importance of both screening mechanisms. The question of how to deal with the repartitioning of the capture flux into allowed capture channels is addressed. The model also provides insights for data on highly-charged projectile ions (C$^{6+}$, O$^{8+}$, Si$^{13+}$) in the MeV/u range where the question of saturation effects in net ionization was raised in the literature.


title: Self-consistent DFT+$U$+$V$ study of oxygen vacancies in SrTiO$_3$

abstract: Contradictory theoretical results for oxygen vacancies in SrTiO$_3$ (STO) were often related to the peculiar properties of STO, which is a $d^0$ transition metal oxide with mixed ionic-covalent bonding. Here, we apply, for the first time, density functional theory (DFT) within the extended Hubbard DFT+$U$+$V$ approach, including on-site as well as inter-site electronic interactions, to study oxygen-deficient STO with Hubbard $U$ and $V$ parameters computed self-consistently via density-functional perturbation theory. Our results demonstrate that the extended Hubbard functional is a promising approach to study defects in materials with electronic properties similar to STO. Indeed, DFT+$U$+$V$ provides a better description of stoichiometric STO compared to standard DFT or DFT+$U$, the band gap and crystal field splitting being in good agreement with experiments. In turn, also the description of the electronic properties of oxygen vacancies in STO is improved, with formation energies in excellent agreement with experiments as well as results obtained with the most frequently used hybrid functionals, however at a fraction of the computational cost. While our results do not fully resolve the contradictory findings reported in literature, our systematic approach leads to a deeper understanding of their origin, which stems from different cell sizes, STO phases, the exchange-correlation functional, and the treatment of structural relaxations and spin-polarization.


title: Neighborhood Structure Assisted Non-negative Matrix Factorization and its Application in Unsupervised Point-wise Anomaly Detection

abstract: Dimensionality reduction is considered as an important step for ensuring competitive performance in unsupervised learning such as anomaly detection. Non-negative matrix factorization (NMF) is a popular and widely used method to accomplish this goal. But NMF do not have the provision to include the neighborhood structure information and, as a result, may fail to provide satisfactory performance in presence of nonlinear manifold structure. To address that shortcoming, we propose to consider and incorporate the neighborhood structural similarity information within the NMF framework by modeling the data through a minimum spanning tree. We label the resulting method as the neighborhood structure assisted NMF. We further devise both offline and online algorithmic versions of the proposed method. Empirical comparisons using twenty benchmark datasets as well as an industrial dataset extracted from a hydropower plant demonstrate the superiority of the neighborhood structure assisted NMF and support our claim of merit. Looking closer into the formulation and properties of the neighborhood structure assisted NMF with other recent, enhanced versions of NMF reveals that inclusion of the neighborhood structure information using MST plays a key role in attaining the enhanced performance in anomaly detection.


title: Finite-size scaling analysis of protein droplet formation

abstract: The formation of biomolecular condensates inside cells often involve intrinsically disordered proteins (IDPs), and several of these IDPs are also capable of forming droplet-like dense assemblies on their own, through liquid-liquid phase separation. When modelings thermodynamic phase changes, it is well-known that finite-size scaling analysis can be a valuable tool. However, to our knowledge, this approach has not been applied before to the computationally challenging problem of modeling sequence-dependent biomolecular phase separation. Here, we implement finite-size scaling methods to investigate the phase behavior of two 10-bead sequences in a continuous hydrophobic/polar protein model. Combined with reversible explicit-chain Monte Carlo simulations of these sequences, finite-size scaling analysis turns out to be both feasible and rewarding, despite relying on theoretical results for asymptotically large systems. While both sequences form dense clusters at low temperature, this analysis shows that only one of them undergoes liquid-liquid phase separation. Furthermore, the transition temperature at which droplet formation sets in, is observed to converge slowly with system size, so that even for our largest systems the transition is shifted by about 8%. Using finite-size scaling analysis, this shift can be estimated and corrected for.


title: Siamese Graph Neural Networks for Data Integration

abstract: Data integration has been studied extensively for decades and approached from different angles. However, this domain still remains largely rule-driven and lacks universal automation. Recent development in machine learning and in particular deep learning has opened the way to more general and more efficient solutions to data integration problems. In this work, we propose a general approach to modeling and integrating entities from structured data, such as relational databases, as well as unstructured sources, such as free text from news articles. Our approach is designed to explicitly model and leverage relations between entities, thereby using all available information and preserving as much context as possible. This is achieved by combining siamese and graph neural networks to propagate information between connected entities and support high scalability. We evaluate our method on the task of integrating data about business entities, and we demonstrate that it outperforms standard rule-based systems, as well as other deep learning approaches that do not use graph-based representations.


title: Exploratory Projection to Latent Structure Models for use in Transcriptomic Analysis

abstract: In this paper, we ask if it is possible to increase the interpretability in multivariate analysis by aligning and projecting covariates onto comparative subspaces. We demonstrate our method as well as the interpretative power of PLS decomposed models and how robust interpretability can lead to quantitative insights. We discuss the statistical properties of the PLS weights, $p$-values associated with specific axes, as well as their alignment properties. The applicability of this approach within life science is also demonstrated by applying it to three use cases of publically available datasets. Further we present hierarchical pathway enrichment results stemming from aligned $p$-values, which are compared with results derived from enrichment analysis, as an external validation of our method. We find that the method can uncover known results from genomics for all of the studied use cases, i.e. microarray data from multiple sclerosis and diabetes patients as well as RNA sequencing data from breast cancer patients.


title: Machine learning and AI-based approaches for bioactive ligand discovery and GPCR-ligand recognition

abstract: In the last decade, machine learning and artificial intelligence applications have received a significant boost in performance and attention in both academic research and industry. The success behind most of the recent state-of-the-art methods can be attributed to the latest developments in deep learning. When applied to various scientific domains that are concerned with the processing of non-tabular data, for example, image or text, deep learning has been shown to outperform not only conventional machine learning but also highly specialized tools developed by domain experts. This review aims to summarize AI-based research for GPCR bioactive ligand discovery with a particular focus on the most recent achievements and research trends. To make this article accessible to a broad audience of computational scientists, we provide instructive explanations of the underlying methodology, including overviews of the most commonly used deep learning architectures and feature representations of molecular data. We highlight the latest AI-based research that has led to the successful discovery of GPCR bioactive ligands. However, an equal focus of this review is on the discussion of machine learning-based technology that has been applied to ligand discovery in general and has the potential to pave the way for successful GPCR bioactive ligand discovery in the future. This review concludes with a brief outlook highlighting the recent research trends in deep learning, such as active learning and semi-supervised learning, which have great potential for advancing bioactive ligand discovery.


title: Privacy Amplification of Iterative Algorithms via Contraction Coefficients

abstract: We investigate the framework of privacy amplification by iteration, recently proposed by Feldman et al., from an information-theoretic lens. We demonstrate that differential privacy guarantees of iterative mappings can be determined by a direct application of contraction coefficients derived from strong data processing inequalities for $f$-divergences. In particular, by generalizing the Dobrushin's contraction coefficient for total variation distance to an $f$-divergence known as $E_{\gamma}$-divergence, we derive tighter bounds on the differential privacy parameters of the projected noisy stochastic gradient descent algorithm with hidden intermediate updates.


title: Predictability limit of partially observed systems

abstract: Applications from finance to epidemiology and cyber-security require accurate forecasts of dynamic phenomena, which are often only partially observed. We demonstrate that a system's predictability degrades as a function of temporal sampling, regardless of the adopted forecasting model. We quantify the loss of predictability due to sampling, and show that it cannot be recovered by using external signals. We validate the generality of our theoretical findings in real-world partially observed systems representing infectious disease outbreaks, online discussions, and software development projects. On a variety of prediction tasks---forecasting new infections, the popularity of topics in online discussions, or interest in cryptocurrency projects---predictability irrecoverably decays as a function of sampling, unveiling fundamental predictability limits in partially observed systems.


title: Who voted for a No Deal Brexit? A Composition Model of Great Britains 2019 European Parliamentary Elections

abstract: The purpose of this paper is to use the votes cast at the 2019 European elections held in United Kingdom to re-visit the analysis conducted subsequent to its 2016 European Union referendum vote. This exercise provides a staging post on public opinion as the United Kingdom moves to leave the European Union during 2020. A composition data analysis in a seemingly unrelated regression framework is adopted that respects the compositional nature of the vote outcome; each outcome is a share that adds up to 100% and each outcome is related to the alternatives. Contemporary explanatory data for each counting area is sourced from the themes of socio-demographics, employment, life satisfaction and place. The study find that there are still strong and stark divisions in the United Kingdom, defined by age, qualifications, employment and place. The use of a compositional analysis approach produces challenges in regards to the interpretation of these models, but marginal plots are seen to aid the interpretation somewhat.


title: Correlations between Background Radiation inside a Multilayer Interleaving Structure, Geomagnetic Activity, and Cosmic Radiation: A Fourth Order Cumulant-based Correlation Analysis

abstract: In this work, we analyzed time-series of background radiation inside a multilayer interleaving structure, geomagnetic activity and cosmic-ray activity using the Pearson correlation coefficient and a new correlation measure based on the one-dimensional component of the fourth order cumulant. The new method is proposed based on the fact that the cumulant of a random process is zero if it is of Gaussian nature. The results show that this methodology is useful for detecting correlations between the analyzed variables.


title: Lower density selection schemes via small universal hitting sets with short remaining path length

abstract: Universal hitting sets are sets of words that are unavoidable: every long enough sequence is hit by the set (i.e., it contains a word from the set). There is a tight relationship between universal hitting sets and minimizers schemes, where minimizers schemes with low density (i.e., efficient schemes) correspond to universal hitting sets of small size. Local schemes are a generalization of minimizers schemes which can be used as replacement for minimizers scheme with the possibility of being much more efficient. We establish the link between efficient local schemes and the minimum length of a string that must be hit by a universal hitting set. We give bounds for the remaining path length of the Mykkeltveit universal hitting set. Additionally, we create a local scheme with the lowest known density that is only a log factor away from the theoretical lower bound.


title: Optical follow-up of gravitational wave triggers with DECam during the first two LIGO/VIRGO observing runs

abstract: Gravitational wave (GW) events detectable by LIGO and Virgo have several possible progenitors, including black hole mergers, neutron star mergers, black hole--neutron star mergers, supernovae, and cosmic string cusps. A subset of GW events are expected to produce electromagnetic (EM) emission that, once detected, will provide complementary information about their astrophysical context. To that end, the LIGO--Virgo Collaboration (LVC) sends GW candidate alerts to the astronomical community so that searches for their EM counterparts can be pursued. The DESGW group, consisting of members of the Dark Energy Survey (DES), the LVC, and other members of the astronomical community, uses the Dark Energy Camera (DECam) to perform a search and discovery program for optical signatures of LVC GW events. DESGW aims to use a sample of GW events as standard sirens for cosmology. Due to the short decay timescale of the expected EM counterparts and the need to quickly eliminate survey areas with no counterpart candidates, it is critical to complete the initial analysis of each night's images as quickly as possible. We discuss our search area determination, imaging pipeline, and candidate selection processes. We review results from the DESGW program during the first two LIGO--Virgo observing campaigns and introduce other science applications that our pipeline enables.


title: Positive Hankel operators, positive definite kernels and related topics

abstract: It is shown that a positive (bounded linear) operator on a Hilbert space with trivial kernel is unitarily equivalent to a Hankel operator that satisfies double positivity condition if and only if it is non-invertible and has simple spectrum (that is, if this operator admits a cyclic vector). More generally, for an arbitrary positive (bounded linear) operator A on a Hilbert space H with trivial kernel the collection V(A) of all linear isometries V from H into H such that AV is positive as well is investigated. In particular, operators A such that V(A) contains a pure isometry with a given deficiency index are characterized. Some applications to unbounded positive self-adjoint operators as well as to positive definite kernels are presented. In particular, positive definite matrix-type square roots of such kernels are studied and kernels that have a unique such root are characterized. The class of all positive definite kernels that have at least one such a square root is also investigated.


title: Measurements with silicon photomultipliers of dose-rate effects in the radiation damage of plastic scintillator tiles in the CMS hadron endcap calorimeter

abstract: Measurements are presented of the reduction of signal output due to radiation damage for plastic scintillator tiles used in the hadron endcap (HE) calorimeter of the CMS detector. The tiles were exposed to particles produced in proton-proton (pp) collisions at the CERN LHC with a center-of-mass energy of 13 TeV, corresponding to a delivered luminosity of 50 fb$^{-1}$. The measurements are based on readout channels of the HE that were instrumented with silicon photomultipliers, and are derived using data from several sources: a laser calibration system, a movable radioactive source, as well as hadrons and muons produced in pp collisions. Results from several irradiation campaigns using $^{60}$Co sources are also discussed. The damage is presented as a function of dose rate. Within the range of these measurements, for a fixed dose the damage increases with decreasing dose rate.


title: PARAFAC-Based Channel Estimation for Intelligent Reflective Surface Assisted MIMO System

abstract: Intelligent reflective surface (IRS) is an emergent technology for future wireless communications. It consists of a large 2D array of passive scattering elements that control the electromagnetic properties of radio-frequency waves so that the reflected signals add coherently at the intended receiver or destructively to reduce co-channel interference. The promised gains of IRS-assisted communications depend on the accuracy of the channel state information. In this paper, we propose two novel channel estimation methods for an IRS-assisted multiple-input multiple-output (MIMO) communication system. Assuming a structured time-domain pattern of pilots and IRS phase shifts, we show that the received signal follows a parallel factor (PARAFAC) tensor model that can be exploited to estimate the involved communication channels in closed-form or iteratively. Numerical results corroborate the effectiveness of the proposed channel estimation methods and highlight the involved tradeoffs.


title: Counterexamples to "The Blessings of Multiple Causes" by Wang and Blei

abstract: This note has been updated (April, 2020) to respond to "Towards Clarifying the Theory of the Deconfounder" by Yixin Wang, David M. Blei (arXiv:2003.04948). This original note, posted in January, 2020, is meant to complement our previous comment on "The Blessings of Multiple Causes" by Wang and Blei (2019). We provide a more succinct and transparent explanation of the fact that the deconfounder does not control for multi-cause confounding. The argument given in Wang and Blei (2019) makes two mistakes: (1) attempting to infer independence conditional on one variable from independence conditional on a different, unrelated variable, and (2) attempting to infer joint independence from pairwise independence. We give two simple counterexamples to the deconfounder claim.


title: Using Double Contractions to Derive the Structure of Slice-Wise Multiplications of Tensors with Applications to Semi-Blind MIMO OFDM

abstract: The slice-wise multiplication of two tensors is required in a variety of tensor decompositions (including PARAFAC2 and PARATUCK2) and is encountered in many applications, including the analysis of multidimensional biomedical data (EEG, MEG, etc.) or multi-carrier MIMO systems. In this paper, we propose a new tensor representation that is not based on a slice-wise (matrix) description, but can be represented by a double contraction of two tensors. Such a double contraction of two tensors can be efficiently calculated via generalized unfoldings. It leads to new tensor models of the investigated system that do not depend on the chosen unfolding and reveal the tensor structure of the data model (such that all possible unfoldings can be seen at the same time). As an example, we apply this new concept to the design of new receivers for multi-carrier MIMO systems in wireless communications. In particular, we consider MIMO OFDM systems with and without Khatri-Rao coding. The proposed receivers exploit the channel correlation between adjacent subcarriers, require the same amount of training symbols as traditional OFDM techniques, but have an improved performance in terms of the symbol error rate. Furthermore, we show that the spectral efficiency of the Khatri-Rao coded MIMO-OFDM can be increased by introducing "random coding" such that the "coding matrix" also contains useful information symbols. Considering this transmission technique, we derive a tensor model and two types of receivers for randomly coded MIMO-OFDM systems using the double contraction of two tensors.


title: Magic Cayley-Sudoku Tables

abstract: A Cayley-sudoku table of a finite group G is a Cayley table for G, the body of which is partitioned into uniformly sized rectangular blocks in such a way that each group element appears exactly once in each block. A Cayley-sudoku table is pandiagonal magic provided the blocks are square and the "sum" (using the group operation) of elements in every row, column, broken diagonal, and broken antidiagonal in each block gives the group identity. We provide sufficient conditions for a group to have a pandiagonal magic Cayley-sudoku table and give some examples.


title: Simplified 3D GCM modelling of the irradiated brown dwarf WD0137-349B

abstract: Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-IR emission excess at $\gtrsim$ 1.95 $\mu$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: We present a first attempt at simulating the atmosphere of a short period White dwarf - Brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the UV irradiation is required to more accurately capture the energy balance inside the Brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the Brown dwarf.


title: Single-layer graphene on epitaxial FeRh thin films

abstract: Graphene is a 2D material that displays excellent electronic transport properties with prospective applications in many fields. Inducing and controlling magnetism in the graphene layer, for instance by proximity of magnetic materials, may enable its utilization in spintronic devices. This paper presents fabrication and detailed characterization of single-layer graphene formed on the surface of epitaxial FeRh thin films. The magnetic state of the FeRh surface can be controlled by temperature, magnetic field or strain due to interconnected order parameters. Characterization of graphene layers by X-ray Photoemission and X-ray Absorption Spectroscopy, Low-Energy Ion Scattering, Scanning Tunneling Microscopy, and Low-Energy Electron Microscopy shows that graphene is single-layer, polycrystalline and covers more than 97% of the substrate. Graphene displays several preferential orientations on the FeRh(001) surface with unit vectors of graphene rotated by 30{\deg}, 15{\deg}, 11{\deg}, and 19{\deg} with respect to FeRh substrate unit vectors. In addition, the graphene layer is capable to protect the films from oxidation when exposed to air for several months. Therefore, it can be also used as a protective layer during fabrication of magnetic elements or as an atomically thin spacer, which enables incorporation of switchable magnetic layers within stacks of 2D materials in advanced devices.


title: Do model emission line galaxies live in filaments at z~1?

abstract: Current and future cosmological surveys are targeting star-forming galaxies at $z\sim 1$ with nebular emission lines. We use a state-of-the-art semi-analytical model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model ELGs are selected such that they can be compared directly with the DEEP2, VVDS, eBOSS-SGC and DESI surveys. The large scale environment of the ELGs is classified using velocity-shear-tensor and tidal-tensor algorithms. Half of the model ELGs live in filaments and about a third in sheets. Model ELGs which reside in knots have the largest satellite fractions. We find that the shape of the mean halo occupation distribution of model ELGs varies widely for different large scale environments. To interpret our results, we also study fixed number density samples of ELGs and galaxies selected using simpler criteria, with single cuts in stellar mass, star formation rate and [OII] luminosity. The fixed number density ELG selection produces samples that are close to L[OII] and SFR selected samples for densities above $10^{-4.2}h^{3}{\rm Mpc}^{-3}$. ELGs with an extra cut in stellar mass applied to fix their number density, present differences in sheets and knots with respect to the other samples. ELGs, SFR and L[OII] selected samples with equal number density have similar large scale bias but their clustering below separations of $1h^{-1}$Mpc is different.


title: Circllhist -- A Log-Linear Histogram Data Structure for IT Infrastructure Monitoring

abstract: The circllhist histogram is a fast and memory efficient data structure for summarizing large numbers of latency measurements. It is particularly suited for applications in IT infrastructure monitoring, and provides nano-second data insertion, full mergeability, accurate approximation of quantiles with a-priori bounds on the relative error. Open-source implementations are available for C/lua/python/Go/Java/JavaScript.


title: Set-Membership Filter for Discrete-Time Nonlinear Systems Using State Dependent Coefficient Parameterization

abstract: In this technical note, a recursive set-membership filtering algorithm for discrete-time nonlinear dynamical systems subject to unknown but bounded process and measurement noises is proposed. The nonlinear dynamics is represented in a pseudo-linear form using the state dependent coefficient (SDC) parameterization. Matrix Taylor expansions are utilized to expand the state dependent matrices about the state estimates. Upper bounds on the norms of remainders in the matrix Taylor expansions are calculated on-line using a non-adaptive random search algorithm at each time step. Utilizing these upper bounds and the ellipsoidal set description of the uncertainties, a two-step filter is derived that utilizes the `correction-prediction' structure of the standard Kalman Filter variants. At each time step, correction and prediction ellipsoids are constructed that contain the true state of the system by solving the corresponding semi-definite programs (SDPs). Finally, a simulation example is included to illustrate the effectiveness of the proposed approach.


title: Energy, work, entropy and heat balance in Marcus molecular junctions

abstract: We present a consistent theory of energy balance and conversion in a single-molecule junction with strong interactions between electrons on the molecular linker (dot) and phonons in the nuclear environment where the Marcus-type electron hopping processes predominate in the electron transport. It is shown that the environmental reorganization and relaxation that accompany electron hopping energy exchange between the electrodes and the nuclear (molecular and solvent) environment may bring a moderate local cooling of the latter in biased systems. The effect of a periodically driven dot level on the heat transport and power generated in the system is analyzed and energy conservation is demonstrated both within and beyond the quasistatic regime. Finally, a simple model of atomic scale engine based on a Marcus single-molecule junction with a driven electron level is suggested and discussed.


title: Media Forensics and DeepFakes: an overview

abstract: With the rapid progress of recent years, techniques that generate and manipulate multimedia content can now guarantee a very advanced level of realism. The boundary between real and synthetic media has become very thin. On the one hand, this opens the door to a series of exciting applications in different fields such as creative arts, advertising, film production, video games. On the other hand, it poses enormous security threats. Software packages freely available on the web allow any individual, without special skills, to create very realistic fake images and videos. So-called deepfakes can be used to manipulate public opinion during elections, commit fraud, discredit or blackmail people. Potential abuses are limited only by human imagination. Therefore, there is an urgent need for automated tools capable of detecting false multimedia content and avoiding the spread of dangerous false information. This review paper aims to present an analysis of the methods for visual media integrity verification, that is, the detection of manipulated images and videos. Special emphasis will be placed on the emerging phenomenon of deepfakes and, from the point of view of the forensic analyst, on modern data-driven forensic methods. The analysis will help to highlight the limits of current forensic tools, the most relevant issues, the upcoming challenges, and suggest future directions for research.


title: Emergence of Topologically Non-trivial Spin-polarized States in a Segmented Linear Chain

abstract: The synthesis of new materials with novel or useful properties is one of the most important drivers in the fields of condensed matter physics and materials science. Discoveries of this kind are especially significant when they point to promising future basic research and applications. Van der Waals bonded materials comprised of lower-dimensional building blocks have been shown to exhibit emergent properties when isolated in an atomically thin form1-8. Here, we report the discovery of a transition metal chalcogenide in a heretofore unknown segmented linear chain form, where basic building blocks each consisting of two hafnium atoms and nine tellurium atoms (Hf2Te9) are van der Waals bonded end-to-end. First-principle calculations based on density functional theory reveal striking crystal-symmetry-related features in the electronic structure of the segmented chain, including giant spin splitting and nontrivial topological phases of selected energy band states. Atomic-resolution scanning transmission electron microscopy reveals single segmented Hf2Te9 chains isolated within the hollow cores of carbon nanotubes, with a structure consistent with theoretical predictions. Van der Waals-bonded segmented linear chain transition metal chalcogenide materials could open up new opportunities in low-dimensional, gate-tunable, magnetic and topological crystalline systems.


title: Microscopic aspects of artificial ageing in Al-Mg-Si alloys

abstract: Al-Mg-Si alloys with total solute contents ranging from 0.8 to 1.4 wt.% were solutionised, quenched and then artificially aged (AA) at 180 {\deg}C, after which positron annihilation lifetime spectroscopy was applied to obtain information about precipitation and vacancy evolution during preceding ageing. Hardness and electrical resistivity measurements were carried out to complement these measurements. AA was carried out in four different heating media, which allowed for varying the heating rate from 2.4 K/s to 170 K/s. The main result of the study is that there is a competition between vacancy losses and precipitation. Any precipitation taking place during quenching or during heating to the AA temperature helps to prevent vacancies from going to sinks and allows them to assist in solute clustering. Higher solute content, slower heating to 180 {\deg}C and natural pre-ageing before AA were found to have a comparable effect.


title: A tail dependence-based MST and their topological indicators in modelling systemic risk in the European insurance sector

abstract: In the present work we analyse the dynamics of indirect connections between insurance companies that result from market price channels. In our analysis we assume that the stock quotations of insurance companies reflect market sentiments which constitute a very important systemic risk factor. Interlinkages between insurers and their dynamics have a direct impact on systemic risk contagion in the insurance sector. We propose herein a new hybrid approach to the analysis of interlinkages dynamics based on combining the copula-DCC-GARCH model and Minimum Spanning Trees (MST). Using the copula-DCC-GARCH model we determine the tail dependence coefficients. Then, for each analysed period we construct MST based on these coefficients. The dynamics is analysed by means of time series of selected topological indicators of the MSTs in the years 2005-2019. Our empirical results show the usefulness of the proposed approach to the analysis of systemic risk in the insurance sector. The times series obtained from the proposed hybrid approach reflect the phenomena occurring on the market. The analysed MST topological indicators can be considered as systemic risk predictors.


title: Chiral Ward identities for Dirac eigenmodes with staggered fermions

abstract: We study chiral properties of eigenvalue spectrum for staggered quarks. We present a new method to identify would-be zero modes and nonzero modes using their symmetry and chiral properties. Here, we review the traditional method with HYP improved staggered quarks, and extend it to a completely new method which uses the chiral Ward identities and leakage patterns to achieve the goal.


title: Temporally Varying Universal Gravitational "Constant" and Speed of Light in Energy Momentum Squared Gravity

abstract: Energy Momentum Squared Gravity (EMSG) is a cosmological model where the scale factor is non vanishing at all times and hence does not favor big bang cosmology. However, the profile of density in the radiation dominated universe shows that EMSG supports inflationary cosmology. Inflationary cosmological models are successful in providing convincing answers to major cosmological issues like horizon problem, flatness problem and small value of cosmological constant but hitherto no model of inflation has been observationally confirmed. Owing to this, Varying Speed of Light (VSL) were introduced which are a class of cosmological models which disfavor inflation and propose an alternative route to solve these cosmological issues by just allowing the speed of light (and Newtonian Gravitational constant) to vary. VSL theories were motivated to address the shortcomings of inflation but do not address the shortcomings related to the initial big bang singularity. In this spirit, we present here a novel cosmological model which is free from both the "initial big bang singularity" and "inflation" by incorporating a mutually varying speed of light $c(t)$ and Newtonian gravitational constant $G(t)$ in the framework of EMSG. We report that in EMSG, for a dust universe ($\omega=0$), cosmological models for a time varying $c(t)$ and $G(t)$ and constant $c$ and $G$ are indistinguishable, whereas for a radiation dominated universe ($\omega = 1/3$), a mutually varying $c(t)$ and $G(t)$ provides an exiting alternative to inflationary cosmology which is also free from initial big bang singularity. We further report that for an ansatz of scale factor representing a bouncing cosmological model, the VSL theory can be applied to a quadratic $T$ gravity model to get rid of "inflation" and "big bang singularity" and concurrently solve the above mentioned cosmological enigmas.


title: Harmonic Convolutional Networks based on Discrete Cosine Transform

abstract: Convolutional neural networks (CNNs) learn filters in order to capture local correlation patterns in feature space. We propose to learn these filters as combinations of preset spectral filters defined by the Discrete Cosine Transform (DCT). Our proposed DCT-based harmonic blocks replace conventional convolutional layers to produce partially or fully harmonic versions of new or existing CNN architectures. Using DCT energy compaction properties, we demonstrate how the harmonic networks can be efficiently compressed by truncating high-frequency information in harmonic blocks thanks to the redundancies in the spectral domain. We report extensive experimental validation demonstrating benefits of the introduction of harmonic blocks into state-of-the-art CNN models in image classification, object detection and semantic segmentation applications.


title: Persistent spin textures and currents in wurtzite nanowire-based quantum structures

abstract: We explore the spin and charge properties of electrons in wurtzite semiconductor nanowires where radial and axial confinement leads to tubular or ring-shaped quantum structures. Accounting for spin-orbit interaction induced by the wurtzite lattice as well as a radial potential gradient, we analytically derive the corresponding low-dimensional Hamiltonians. It is demonstrated that the resulting tubular spin-orbit Hamiltonian allows to construct spin states that are persistent in time and robust against disorder. We find that these special scenarios are characterized by distinctive features in the optical conductivity spectrum, which enable an unambiguous experimental verification. In both types of quantum structures, we discuss the dependence of the occurring persistent charge and spin currents on an axial magnetic field and Fermi energy which show clear fingerprints of the electronic subband structure. Here, the spin-preserving symmetries become manifest in the vanishing of certain spin current tensor components. Our analytic description relates the distinctive features of the optical conductivity and persistent currents to bandstructure characteristics which allows to deduce spin-orbit coefficients and other band parameters from measurements.


title: Light gauge boson interpretation for $(g-2)_\mu$ and the $K_L \rightarrow \pi^0 + \text{(invisible)}$ anomaly at the J-PARC KOTO experiment

abstract: We discuss a list of possible light gauge boson interpretations for the long-standing experimental anomaly in $(g-2)_\mu$ and also recent anomalous excess in $K_L \rightarrow \pi^0 + \text{(invisible)}$ events at the J-PARC KOTO experiment. We consider two models: $i$) $L_\mu - L_\tau$ gauge boson with heavy vector-like quarks and $ii$) $(L_\mu - L_\tau) + \epsilon (B_3 - L_\tau)$ gauge boson in the presence of right-handed neutrinos. When the light gauge boson has mass close to the neutral pion in order to satisfy the Grossman-Nir bound, the models successfully explain the anomalies simultaneously while satisfying all known experimental constraints. We extensively provide the future prospect of suggested models.


title: Homophily based on few attributes can impede structural balance

abstract: Two complementary mechanisms are thought to shape social groups: homophily between agents and structural balance in connected triads. Here we consider $N$ fully connected agents, where each agent has $G$ underlying attributes, and the similarity between agents in attribute space (i.e., homophily) is used to determine the link weight between them. To incorporate structural balance we use a triad-updating rule where only one attribute of one agent is changed intentionally in each update, but this also leads to accidental changes in link weights and even link polarities. The link weight dynamics in the limit of large $G$ is described by a Fokker-Planck equation from which the conditions for a phase transition to a fully balanced state with all links positive can be obtained. This "paradise state" of global cooperation is, however, difficult to achieve requiring $G > O(N^2)$ and $p>0.5$, where the parameter $p$ captures a willingness to consensus. Allowing edge weights to be a consequence of attributes naturally captures homophily and reveals that many real-world social systems would have a subcritical number of attributes necessary to achieve structural balance.


title: Unstable Entropy and Unstable Pressure for Partially Hyperbolic Endomorphisms

abstract: In this paper, unstable metric entropy, unstable topological entropy and unstable pressure for partially hyperbolic endomorphisms are introduced and investigated. A version of Shannon-McMillan-Breiman Theorem is established, and a variational principle is formulated, which gives a relationship between unstable metric entropy and unstable pressure (unstable topological entropy). As an application of the variational principle, some results on the $u$-equilibrium states are given.


title: Subdivided Phase Oracle for NISQ Search Algorithms

abstract: Because noisy, intermediate-scale quantum (NISQ) machines accumulate errors quickly, we need new approaches to designing NISQ-aware algorithms and assessing their performance. Algorithms with characteristics that appear less desirable under ideal circumstances, such as lower success probability, may in fact outperform their ideal counterparts on existing hardware. We propose an adaptation of Grover's algorithm, subdividing the phase flip into segments to replace a digital counter and complex phase flip decision logic. We applied this approach to obtaining the best solution of the MAX-CUT problem in sparse graphs, utilizing multi-control, Toffoli-like gates with residual phase shifts. We implemented this algorithm on IBM Q processors and succeeded in solving a 5-node MAX-CUT problem, demonstrating amplitude amplification on four qubits. This approach will be useful for a range of problems, and may shorten the time to reaching quantum advantage.


title: Inference for Network Structure and Dynamics from Time Series Data via Graph Neural Network

abstract: Network structures in various backgrounds play important roles in social, technological, and biological systems. However, the observable network structures in real cases are often incomplete or unavailable due to measurement errors or private protection issues. Therefore, inferring the complete network structure is useful for understanding complex systems. The existing studies have not fully solved the problem of inferring network structure with partial or no information about connections or nodes. In this paper, we tackle the problem by utilizing time series data generated by network dynamics. We regard the network inference problem based on dynamical time series data as a problem of minimizing errors for predicting future states and proposed a novel data-driven deep learning model called Gumbel Graph Network (GGN) to solve the two kinds of network inference problems: Network Reconstruction and Network Completion. For the network reconstruction problem, the GGN framework includes two modules: the dynamics learner and the network generator. For the network completion problem, GGN adds a new module called the States Learner to infer missing parts of the network. We carried out experiments on discrete and continuous time series data. The experiments show that our method can reconstruct up to 100% network structure on the network reconstruction task. While the model can also infer the unknown parts of the structure with up to 90% accuracy when some nodes are missing. And the accuracy decays with the increase of the fractions of missing nodes. Our framework may have wide application areas where the network structure is hard to obtained and the time series data is rich.