abstract: We investigate Bayesian and frequentist approaches to resonance searches using a toy model based on an ATLAS search for the Higgs boson in the diphoton channel. We draw pseudo-data from the background only model and background plus signal model at multiple luminosities, from $10^{-3}$/fb to $10^7$/fb. We chart the change in the Bayesian posterior of the background only model and the global p-value. We find that, as anticipated, the posterior converges to certainty about the model as luminosity increases. The p-value, on the other hand, randomly walks between 0 and 1 if the background only model is true, and otherwise converges to 0. After briefly commenting on the frequentist properties of the posterior, we make a direct comparison of the significances obtained in Bayesian and frequentist frameworks. We find that the well-known look-elsewhere effect reduces local significances by about 1$\sigma$. We furthermore find that significances from our Bayesian framework are typically about 1 to 2$\sigma$ smaller than the global significances, though the reduction depends on the prior, global significance and integrated luminosity. This suggests that even global significances could significantly overstate the evidence against the background only model. We checked that this effect --- the Bayes effect --- was robust with respect to fourteen choices of prior and investigated the Jeffreys-Lindley paradox for three of them.


title: Multi-wavelength properties of Type 1 and Type 2 AGN Host Galaxies in the Chandra-COSMOS Legacy Survey

abstract: We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass.


title: Ask Not What AI Can Do, But What AI Should Do: Towards a Framework of Task Delegability

abstract: While artificial intelligence (AI) holds promise for addressing societal challenges, issues of exactly which tasks to automate and to what extent to do so remain understudied. We approach this problem of task delegability from a human-centered perspective by developing a framework on human perception of task delegation to AI. We consider four high-level factors that can contribute to a delegation decision: motivation, difficulty, risk, and trust. To obtain an empirical understanding of human preferences in different tasks, we build a dataset of 100 tasks from academic papers, popular media portrayal of AI, and everyday life, and administer a survey based on our proposed framework. We find little preference for full AI control and a strong preference for machine-in-the-loop designs, in which humans play the leading role. Among the four factors, trust is the most correlated with human preferences of optimal human-machine delegation. This framework represents a first step towards characterizing human preferences of AI automation across tasks. We hope this work encourages future efforts towards understanding such individual attitudes; our goal is to inform the public and the AI research community rather than dictating any direction in technology development.


title: Millicharged Particles in Liquid Argon Neutrino Experiments

abstract: We investigate the potential of Liquid Argon (LAr) neutrino detectors to search for millicharged particles, a well-motivated extension of the standard model. Detectors located downstream of an intense proton beam that is striking a target may be exposed to a large flux of millicharged particles. Millicharged particles interact primarily through low momentum exchange producing electron recoil events near detector threshold. Recently, sub-MeV detection capabilities were demonstrated by the Fermilab ArgoNeuT detector, a small LAr detector which was exposed to the NuMI neutrino beam. Despite high background rates and its small size, we show that ArgoNeuT is capable of probing unexplored parameter space with its existing dataset. In particular, we show that the excellent spatial resolution in LAr detectors allows rejecting backgrounds by requiring two soft hits that are aligned with the upstream target. We further discuss the prospects of these types of searches in future larger LAr neutrino detectors such as the DUNE near detector.


title: Eigenstate thermalization and quantum chaos in the Holstein polaron model

abstract: The eigenstate thermalization hypothesis (ETH) is a successful theory that provides sufficient criteria for ergodicity in quantum many-body systems. Most studies were carried out for Hamiltonians relevant for ultracold quantum gases and single-component systems of spins, fermions, or bosons. The paradigmatic example for thermalization in solid-state physics are phonons serving as a bath for electrons. This situation is often viewed from an open-quantum system perspective. Here, we ask whether a minimal microscopic model for electron-phonon coupling is quantum chaotic and whether it obeys ETH, if viewed as a closed quantum system. Using exact diagonalization, we address this question in the framework of the Holstein polaron model. Even though the model describes only a single itinerant electron, whose coupling to dispersionless phonons is the only integrability-breaking term, we find that the spectral statistics and the structure of Hamiltonian eigenstates exhibit essential properties of the corresponding random-matrix ensemble. Moreover, we verify the ETH ansatz both for diagonal and offdiagonal matrix elements of typical phonon and electron observables, and show that the ratio of their variances equals the value predicted from random-matrix theory.


title: Spontaneous Josephson $\pi$ Junctions with Topological Superconductors

abstract: We study a junction between two time-reversal-invariant topological superconductors and show this system goes through a series of multiple transitions between a $0$-junction phase, where the free energy has its minimum for a superconducting phase difference of zero, and a $\pi$-junction phase, where the free energy has its minimum for a superconducting phase difference of $\pi$. These transitions occur in the absence of Coulomb blockade or magnetic impurities. Rather, they are driven by the spin orbit coupling in the junction, and can be probed, for example, by measuring the tunneling density of states or the critical current as a function of the junction's length or its Fermi velocity.


title: Insertion Transformer: Flexible Sequence Generation via Insertion Operations

abstract: We present the Insertion Transformer, an iterative, partially autoregressive model for sequence generation based on insertion operations. Unlike typical autoregressive models which rely on a fixed, often left-to-right ordering of the output, our approach accommodates arbitrary orderings by allowing for tokens to be inserted anywhere in the sequence during decoding. This flexibility confers a number of advantages: for instance, not only can our model be trained to follow specific orderings such as left-to-right generation or a binary tree traversal, but it can also be trained to maximize entropy over all valid insertions for robustness. In addition, our model seamlessly accommodates both fully autoregressive generation (one insertion at a time) and partially autoregressive generation (simultaneous insertions at multiple locations). We validate our approach by analyzing its performance on the WMT 2014 English-German machine translation task under various settings for training and decoding. We find that the Insertion Transformer outperforms many prior non-autoregressive approaches to translation at comparable or better levels of parallelism, and successfully recovers the performance of the original Transformer while requiring only logarithmically many iterations during decoding.


title: Amplitudes' Positivity, Weak Gravity Conjecture, and Modified Gravity

abstract: We derive new positivity bounds for scattering amplitudes in theories with a massless graviton in the spectrum in four spacetime dimensions, of relevance for the weak gravity conjecture and modified gravity theories. The bounds imply that extremal black holes are self-repulsive, $M/|Q|<1$ in suitable units, and that they are unstable to decay to smaller extremal black holes, providing an S-matrix proof of the weak gravity conjecture. We also present other applications of our bounds to the effective field theory of axions, $P(X)$ theories, weakly broken galileons, and curved spacetimes.


title: Invariant-equivariant representation learning for multi-class data

abstract: Representations learnt through deep neural networks tend to be highly informative, but opaque in terms of what information they learn to encode. We introduce an approach to probabilistic modelling that learns to represent data with two separate deep representations: an invariant representation that encodes the information of the class from which the data belongs, and an equivariant representation that encodes the symmetry transformation defining the particular data point within the class manifold (equivariant in the sense that the representation varies naturally with symmetry transformations). This approach is based primarily on the strategic routing of data through the two latent variables, and thus is conceptually transparent, easy to implement, and in-principle generally applicable to any data comprised of discrete classes of continuous distributions (e.g. objects in images, topics in language, individuals in behavioural data). We demonstrate qualitatively compelling representation learning and competitive quantitative performance, in both supervised and semi-supervised settings, versus comparable modelling approaches in the literature with little fine tuning.


title: High Performance Computing for gravitational lens modeling: single vs double precision on GPUs and CPUs

abstract: Strong gravitational lensing is a powerful probe of cosmology and the dark matter distribution. Efficient lensing software is already a necessity to fully use its potential and the performance demands will only increase with the upcoming generation of telescopes. In this paper, we study the possible impact of High Performance Computing techniques on a performance-critical part of the widely used lens modeling software LENSTOOL. We implement the algorithm once as a highly optimized CPU version and once with graphics card acceleration for a simple parametric lens model. In addition, we study the impact of finite machine precision on the lensing algorithm. While double precision is the default choice for scientific applications, we find that single precision can be sufficiently accurate for our purposes and lead to a big speedup. Therefore we develop and present a mixed precision algorithm which only uses double precision when necessary. We measure the performance of the different implementations and find that the use of High Performance Computing Techniques dramatically improves the code performance both on CPUs and GPUs. Compared to the current LENSTOOL implementation on 12 CPU cores, we obtain speedup factors of up to 170. We achieve this optimal performance by using our mixed precision algorithm on a high-end GPU which is common in modern supercomputers. We also show that these techniques reduce the energy consumption by up to 98%. Furthermore, we demonstrate that a highly competitive speedup can be reached with consumer GPUs. While they are an order of magnitude cheaper than the high-end graphics cards, they are rarely used for scientific computations due to their low double precision performance. Our mixed precision algorithm unlocks their full potential. The consumer GPU delivers a speedup which is only a factor of four lower than the best speedup achieved by a high-end GPU.


title: Skin Lesion Synthesis with Generative Adversarial Networks

abstract: Skin cancer is by far the most common type of cancer. Early detection is the key to increase the chances for successful treatment significantly. Currently, Deep Neural Networks are the state-of-the-art results on automated skin cancer classification. To push the results further, we need to address the lack of annotated data, which is expensive and require much effort from specialists. To bypass this problem, we propose using Generative Adversarial Networks for generating realistic synthetic skin lesion images. To the best of our knowledge, our results are the first to show visually-appealing synthetic images that comprise clinically-meaningful information.


title: Confirming the Existence of the strong CP Problem in Lattice QCD with the Gradient Flow

abstract: We calculate the electric dipole moment of the nucleon induced by the QCD theta term. We use the gradient flow to define the topological charge and use $N_f = 2+1$ flavors of dynamical quarks corresponding to pion masses of $700$, $570$, and $410$ MeV, and perform an extrapolation to the physical point based on chiral perturbation theory. We perform calculations at $3$ different lattice spacings in the range of $0.07~{\rm fm} < a < 0.11$ fm at a single value of the pion mass, to enable control on discretization effects. We also investigate finite size effects using $2$ different volumes. A novel technique is applied to improve the signal-to-noise ratio in the form factor calculations. The very mild discretization effects observed suggest a continuum-like behavior of the nucleon EDM towards the chiral limit. Under this assumption our results read $d_{n}=-0.00152(71)\ \bar\theta\ e~\text{fm}$ and $d_{p}=0.0011(10)\ \bar\theta\ e~\text{fm}$. Assuming the theta term is the only source of CP violation, the experimental bound on the neutron electric dipole moment limits $\left|\bar\theta\right| < 1.98\times 10^{-10}$ ($90\%$ CL). A first attempt at calculating the nucleon Schiff moment in the continuum resulted in $S_{p} = 0.50(59)\times 10^{-4}\ \bar\theta\ e~\text{fm}^3$ and $S_{n} = -0.10(43)\times 10^{-4}\ \bar\theta\ e~\text{fm}^3$.


title: Topological Control of Synchronization Patterns: Trading Symmetry for Stability

abstract: Symmetries are ubiquitous in network systems and have profound impacts on the observable dynamics. At the most fundamental level, many synchronization patterns are induced by underlying network symmetry, and a high degree of symmetry is believed to enhance the stability of identical synchronization. Yet, here we show that the synchronizability of almost any symmetry cluster in a network of identical nodes can be enhanced precisely by breaking its structural symmetry. This counterintuitive effect holds for generic node dynamics and arbitrary network structure and is, moreover, robust against noise and imperfections typical of real systems, which we demonstrate by implementing a state-of-the-art optoelectronic experiment. These results lead to new possibilities for the topological control of synchronization patterns, which we substantiate by presenting an algorithm that optimizes the structure of individual clusters under various constraints.


title: Practical Enclave Malware with Intel SGX

abstract: Modern CPU architectures offer strong isolation guarantees towards user applications in the form of enclaves. For instance, Intel's threat model for SGX assumes fully trusted enclaves, yet there is an ongoing debate on whether this threat model is realistic. In particular, it is unclear to what extent enclave malware could harm a system. In this work, we practically demonstrate the first enclave malware which fully and stealthily impersonates its host application. Together with poorly-deployed application isolation on personal computers, such malware can not only steal or encrypt documents for extortion, but also act on the user's behalf, e.g., sending phishing emails or mounting denial-of-service attacks. Our SGX-ROP attack uses new TSX-based memory-disclosure primitive and a write-anything-anywhere primitive to construct a code-reuse attack from within an enclave which is then inadvertently executed by the host application. With SGX-ROP, we bypass ASLR, stack canaries, and address sanitizer. We demonstrate that instead of protecting users from harm, SGX currently poses a security threat, facilitating so-called super-malware with ready-to-hit exploits. With our results, we seek to demystify the enclave malware threat and lay solid ground for future research on and defense against enclave malware.


title: Generalized Regular k-point Grid Generation On The Fly

abstract: In the DFT community, it is common practice to use regular k-point grids (Monkhorst-Pack, MP) for Brillioun zone integration. Recently Wisesa et. al.\cite{wisesa2016efficient} and Morgan et. al.\cite{MORGAN2018424} demonstrated that generalized regular (GR) grids offer advantages over traditional MP grids. GR grids have not been widely adopted because one must search through a large number of candidate grids. This work describes an algorithm that can quickly search over GR grids for those that have the most uniform distribution of points and the best symmetry reduction. The grids are ~60% more efficient, on average, than MP grids and can now be generated on the fly in seconds.


title: Work distributions on quantum fields

abstract: We study the work cost of processes in quantum fields without the need of projective measurements, which are always ill defined in quantum field theory. Inspired by interferometry schemes, we propose a work distribution that generalizes the two-point measurement scheme employed in quantum thermodynamics to the case of quantum fields and avoids the use of projective measurements. The distribution is calculated for local unitary processes performed on Kubo-Martin-Schwinger (thermal) states of scalar fields. Crooks theorem and the Jarzynski equality are shown to be satisfied for a family of spatio-temporally localized unitaries, and some features of the resulting distributions are studied as functions of temperature and the degree of localization of the unitary operation. We show how the work fluctuations become much larger than the average as the process becomes more localized in both time and space.


title: Production of $X(3872)$ Accompanied by a Pion in $B$ Meson Decay

abstract: If the $X(3872)$ is a weakly bound charm-meson molecule, it can be produced by the creation of $D^{*0} \bar{D}^0$ or $D^{0} \bar{D}^{*0}$ at short distances followed by the formation of the bound state from the charm-meson pairs. It can also be produced by the creation of $D^{*} \bar{D}^*$ at short distances followed by the rescattering of the charm mesons into $X \pi$. We use results of a previous isospin analysis of $B$ meson decays into $K D^{(*)} \bar D^{(*)}$ to estimate the short-distance amplitudes for creating $D^* \bar D^*$. We use an effective field theory for charm mesons and pions called XEFT to calculate the amplitudes for rescattering of $D^{*} \bar{D}^*$ into $X \pi$ with small relative momentum. The $X\pi$ invariant mass distribution is predicted to have a narrow peak near the $D^{*} \bar{D}^*$ threshold from a charm-meson triangle singularity. We estimate the branching fractions into the peak from the triangle singularity for the decays $B^0 \to K^+ X \pi^-$ and $B^+ \to K^0 X \pi^+$.


title: BMS4 Algebra, Its Stability and Deformations

abstract: We continue analysis of \cite{Parsa:2018kys} and study rigidity and stability of the BMS4 algebra and its centrally extended version. We construct and classify the family of algebras which appear as deformations of BMS4 and in general find the four-parameter family of algebras W(a,b;\bar{a},\bar{b}) as a result of the stabilization analysis, where BMS4=W(-1/2,-1/2;-1/2,-1/2). We then study the W(a,b;\bar{a},\bar{b}) algebra, its maximal finite subgroups and stability for different values of the four parameters. We prove stability of the W(a,b;\bar{a},\bar{b}) family of algebras for generic values of the parameters. For special cases of (a,b)=(\bar{a},\bar{b})=(0,0) and (a,b)=(0,-1), (\bar{a},\bar{b})=(0,0) the algebra can be deformed. In particular we show that centrally extended W(0,-1;0,0) algebra can be deformed to an algebra which has three copies of Virasoro as a subalgebra. We briefly discuss these deformed algebras as asymptotic symmetry algebras and the physical meaning of the stabilization and implications of our result.


title: Type Ia Supernovae are Excellent Standard Candles in the Near-Infrared

abstract: We analyze a set of 89 Type Ia supernovae (SN Ia) that have both optical and near-infrared (NIR) photometry to derive distances and construct low redshift ($z < 0.04$) Hubble diagrams. We construct mean light curve (LC) templates using a hierarchical Bayesian model. We explore both Gaussian process (GP) and template methods for fitting the LCs and estimating distances, while including peculiar velocity and photometric uncertainties. For the 56 SN Ia with both optical and NIR observations near maximum light, the GP method yields a NIR-only Hubble-diagram with a RMS of $0.117 \pm 0.014$ mag when referenced to the NIR maxima. For each NIR band, a comparable GP method RMS is obtained when referencing to NIR-max or B-max. Using NIR LC templates referenced to B-max yields a larger RMS value of $0.138 \pm 0.014$ mag. Fitting the corresponding optical data using standard LC fitters that use LC shape and color corrections yields larger RMS values of $0.179 \pm 0.018$ mag with SALT2 and $0.174 \pm 0.021$ mag with SNooPy. Applying our GP method to subsets of SN Ia NIR LCs at NIR maximum light, even without corrections for LC shape, color, or host-galaxy dust reddening, provides smaller RMS in the inferred distances, at the $\sim 2.3 - 4.1\sigma$ level, than standard optical methods that do correct for those effects. Our ongoing RAISIN program on the Hubble Space Telescope will exploit this promising infrared approach to limit systematic errors when measuring the expansion history of the universe to constrain dark energy.


title: The Cremona group and its subgroups

abstract: This survey deals with the Cremona group via its subgroups.


title: Survival and extinction of epidemics on random graphs with general degrees

abstract: In this paper, we establish the necessary and sufficient criterion for the contact process on Galton-Watson trees (resp. random graphs) to exhibit the phase of extinction (resp. short survival). We prove that the survival threshold $\lambda_1$ for a Galton-Watson tree is strictly positive if and only if its offspring distribution $\xi$ has an exponential tail, i.e., $\mathbb{E} e^{c\xi}<\infty$ for some $c>0$, settling a conjecture by Huang and Durrett [12]. On the random graph with degree distribution $\mu$, we show that if $\mu$ has an exponential tail, then for small enough $\lambda$ the contact process with the all-infected initial condition survives for $n^{1+o(1)}$-time w.h.p. (short survival), while for large enough $\lambda$ it runs over $e^{\Theta(n)}$-time w.h.p. (long survival). When $\mu$ is subexponential, we prove that the contact process w.h.p. displays long survival for any fixed $\lambda>0$.


title: FSNet: Compression of Deep Convolutional Neural Networks by Filter Summary

abstract: We present a novel method of compression of deep Convolutional Neural Networks (CNNs) by weight sharing through a new representation of convolutional filters. The proposed method reduces the number of parameters of each convolutional layer by learning a 1D vector termed Filter Summary (FS). The convolutional filters are located in FS as overlapping 1D segments, and nearby filters in FS share weights in their overlapping regions in a natural way. The resultant neural network based on such weight sharing scheme, termed Filter Summary CNNs or FSNet, has a FS in each convolution layer instead of a set of independent filters in the conventional convolution layer. FSNet has the same architecture as that of the baseline CNN to be compressed, and each convolution layer of FSNet has the same number of filters from FS as that of the basline CNN in the forward process. With compelling computational acceleration ratio, the parameter space of FSNet is much smaller than that of the baseline CNN. In addition, FSNet is quantization friendly. FSNet with weight quantization leads to even higher compression ratio without noticeable performance loss. We further propose Differentiable FSNet where the way filters share weights is learned in a differentiable and end-to-end manner. Experiments demonstrate the effectiveness of FSNet in compression of CNNs for computer vision tasks including image classification and object detection, and the effectiveness of DFSNet is evidenced by the task of Neural Architecture Search.


title: Amorphous carbons in view of multianalytical consideration: normal, expected, and new

abstract: The approach based on selected set of samples and selected set of analytical tools occurred quite efficient when applying to amorphous carbons thus leading to a transformation of the current representation of the issue based on particulars into that one based on a limited set of fixed commonalities. The approach first part implies a set of different origin solid samples. The second part concerns different analytical tools. The combining part means the application of each tool to the whole set of samples. In the current study two natural amorphous carbons shungite carbon and antraxolite, as well as two engineered products carbon blacks CB632 and CB624, all of the four belonging to the elitist highest-carbon-content species, were subjected to analytical study by using modern structural and compositional analytical techniques. The approach has allowed disclosing steady points that are common to the whole class of this carbon allotrope and that may lay the foundation of a consolidate representation of what are amorphous carbons.


title: Discovering Context Effects from Raw Choice Data

abstract: Many applications in preference learning assume that decisions come from the maximization of a stable utility function. Yet a large experimental literature shows that individual choices and judgements can be affected by "irrelevant" aspects of the context in which they are made. An important class of such contexts is the composition of the choice set. In this work, our goal is to discover such choice set effects from raw choice data. We introduce an extension of the Multinomial Logit (MNL) model, called the context dependent random utility model (CDM), which allows for a particular class of choice set effects. We show that the CDM can be thought of as a second-order approximation to a general choice system, can be inferred optimally using maximum likelihood and, importantly, is easily interpretable. We apply the CDM to both real and simulated choice data to perform principled exploratory analyses for the presence of choice set effects.


title: Constructing Selections Stepwise Over Cones of Simplicial Complexes

abstract: It is obtained a natural generalisation of Uspenskij's selection characterisation of paracompact $C$-spaces. The method developed to achieve this result is also applied to give a simplified proof of a similar characterisation of paracompact finite $C$-space obtained previously by Valov. Another application is a characterisation of finite-dimensional paracompact spaces which generalises both a remark done by Michael and a result obtained by the author.


title: Stratified $\beta$-numbers and traveling salesman in Carnot groups

abstract: We introduce a modified version of P. Jones's $\beta$-numbers for Carnot groups which we call {\it stratified $\beta$-numbers}. We show that an analogue of Jones's traveling salesman theorem on 1-rectifiability of sets holds for any Carnot group if we replace previous notions of $\beta$-numbers in Carnot groups with stratified $\beta$-numbers. As we generalize both directions of the traveling salesman theorem, we get a characterization of subsets of Carnot groups that lie on finite length rectifiable curves. Our proof expands upon previous analysis of the Hebisch-Sikora norm for Carnot groups. In particular, we find new estimates on the drift between almost parallel line segments that take advantage of the stratified $\beta$'s and also develop a Taylor expansion technique of the norm. We also give an example of a Carnot group for which a traveling salesman theorem based on the unmodified $\beta$-numbers must exhibit a gap between the necessary and sufficient directions.


title: Dynamically Defined Sequences with Small Discrepancy

abstract: We study the problem of constructing sequences $(x_n)_{n=1}^{\infty}$ on $[0,1]$ in such a way that $$ D_N^* = \sup_{0 \leq x \leq 1} \left| \frac{ \left\{1 \leq i \leq N: x_i \leq x \right\}}{N} - x \right|$$ is uniformly small. A result of Schmidt shows that necessarily $D_N^* \gtrsim (\log{N}) N^{-1}$ for infinitely many $N$ and there are several classical constructions attaining this growth. We describe a type of uniformly distributed sequence that seems to be completely novel: given $\left\{x_1, \dots, x_{N-1} \right\}$, we construct $x_N$ in a greedy manner $$ x_N = \arg\min_{\min_k |x-x_k| \geq N^{-10}} \sum_{k=1}^{N-1}{1-\log{(2\sin{(\pi |x-x_k|)})}}.$$ We prove that $D_N \lesssim (\log{N}) N^{-1/2}$ and conjecture that $D_N \lesssim (\log{N}) N^{-1}$. Numerical examples illustrate this conjecture in a very impressive manner. We also establish a discrepancy bound $D_N \lesssim (\log{N})^d N^{-1/2}$ for an analogous construction in higher dimensions and conjecture it to be $D_N \lesssim (\log{N})^d N^{-1}$.


title: The exit from a metastable state: concentration of the exit point distribution on the low energy saddle points

abstract: We consider the first exit point distribution from a bounded domain $\Omega$ of the stochastic process $(X_t)_{t\ge 0}$ solution to the overdamped Langevin dynamics $$d X_t = -\nabla f(X_t) d t + \sqrt{h} \ d B_t$$ starting from the quasi-stationary distribution in $\Omega$. In the small temperature regime ($h\to 0$) and under rather general assumptions on $f$ (in particular, $f$ may have several critical points in $\Omega$), it is proven that the support of the distribution of the first exit point concentrates on some points realizing the minimum of $f$ on $\partial \Omega$. The proof relies on tools to study tunnelling effects in semi-classical analysis. Extensions of the results to more general initial distributions than the quasi-stationary distribution are also presented.


title: Does the "Artificial Intelligence Clinician" learn optimal treatment strategies for sepsis in intensive care?

abstract: From 2017 to 2018 the number of scientific publications found via PubMed search using the keyword "Machine Learning" increased by 46% (4,317 to 6,307). The results of studies involving machine learning, artificial intelligence (AI), and big data have captured the attention of healthcare practitioners, healthcare managers, and the public at a time when Western medicine grapples with unmitigated cost increases and public demands for accountability. The complexity involved in healthcare applications of machine learning and the size of the associated data sets has afforded many researchers an uncontested opportunity to satisfy these demands with relatively little oversight. In a recent Nature Medicine article, "The Artificial Intelligence Clinician learns optimal treatment strategies for sepsis in intensive care," Komorowski and his coauthors propose methods to train an artificial intelligence clinician to treat sepsis patients with vasopressors and IV fluids. In this post, we will closely examine the claims laid out in this paper. In particular, we will study the individual treatment profiles suggested by their AI Clinician to gain insight into how their AI Clinician intends to treat patients on an individual level.


title: Scalable Holistic Linear Regression

abstract: We propose a new scalable algorithm for holistic linear regression building on Bertsimas & King (2016). Specifically, we develop new theory to model significance and multicollinearity as lazy constraints rather than checking the conditions iteratively. The resulting algorithm scales with the number of samples $n$ in the 10,000s, compared to the low 100s in the previous framework. Computational results on real and synthetic datasets show it greatly improves from previous algorithms in accuracy, false detection rate, computational time and scalability.


title: Learning Ontologies with Epistemic Reasoning: The EL Case

abstract: We investigate the problem of learning description logic ontologies from entailments via queries, using epistemic reasoning. We introduce a new learning model consisting of epistemic membership and example queries and show that polynomial learnability in this model coincides with polynomial learnability in Angluin's exact learning model with membership and equivalence queries. We then instantiate our learning framework to EL and show some complexity results for an epistemic extension of EL where epistemic operators can be applied over the axioms. Finally, we transfer known results for EL ontologies and its fragments to our learning model based on epistemic reasoning.


title: Faster Repetition-Aware Compressed Suffix Trees based on Block Trees

abstract: Suffix trees are a fundamental data structure in stringology, but their space usage, though linear, is an important problem for its applications. We design and implement a new compressed suffix tree targeted to highly repetitive texts, such as large genomic collections of the same species. Our suffix tree tree builds on Block Trees, a recent Lempel-Ziv-bounded data structure that captures the repetitiveness of its input. We use Block Trees to compress the topology of the suffix tree, and augment the Block Tree nodes with data that speeds up suffix tree navigation. Our compressed suffix tree is slightly larger than previous repetition-aware suffix trees based on grammars, but outperforms them in time, often by orders of magnitude. The component that represents the tree topology achieves a speed comparable to that of general-purpose compressed trees, while using 2.3--10 times less space, and might be of interest in other scenarios.


title: Milky Way Halo Vibrations and Incommensurate Stream Velocities

abstract: Collisionless dark matter galactic halos are expected to exhibit damped oscillations as a result of ongoing late time accretion. An n-body model of the cosmological assembly of a Milky Way-like halo is used to quantify the time dependence of its gravitational field. The simulation contains stellar streams whose incommensurate perpendicular velocities are found to have an approximately exponential distribution with a scale of 10-20\kms, depending on how the stars are selected, comparable to those reported for the Orphan stream. The fluctuations in the quadrupole moment of the dark matter halo are sufficient to largely explain the tangential velocities. If velocity measurements of a larger sample of Milky Way streams finds (or does not find) the expected distribution of transverse velocities it will lead to limits on the cross-section of self-interacting dark matter, in which kinetic viscosity can damp the oscillations more rapidly than the mixing processes of collisionless dark matter alone.


title: Exploring the Star Formation Histories of Galaxies in Different Environments from MaNGA Spectra

abstract: The star formation history (SFH) of galaxies allow us to investigate when galaxies formed their stars and assembled their mass. We can constrain the SFH with high level of precision from galaxies with resolved stellar populations, since we are able to discriminate between stars of different ages from the spectrum they emit. However, the relative importance of secular evolution (nature) over nurture is not yet clear, and separating the effects of interaction-driven evolution in the observed galaxy properties is not trivial. The aim of this study is to use MaNGA (Mapping Nearby Galaxies at APO) Integral Field Unit (IFU) data, in combination with multi-wavelength data, to constrain the SFH of nearby isolated galaxies. We present here the new techniques we are developing to constrain the SFH with high level of precision from Spectral Energy Distribution (SED) fitting. This study is part of a China-Chile collaboration program where we are applying these new techniques to investigate how galaxies formed and evolve in different environments.


title: Data-driven network alignment

abstract: Biological network alignment (NA) aims to find a node mapping between species' molecular networks that uncovers similar network regions, thus allowing for transfer of functional knowledge between the aligned nodes. However, current NA methods do not end up aligning functionally related nodes. A likely reason is that they assume it is topologically similar nodes that are functionally related. However, we show that this assumption does not hold well. So, a paradigm shift is needed with how the NA problem is approached. We redefine NA as a data-driven framework, TARA (daTA-dRiven network Alignment), which attempts to learn the relationship between topological relatedness and functional relatedness without assuming that topological relatedness corresponds to topological similarity, like traditional NA methods do. TARA trains a classifier to predict whether two nodes from different networks are functionally related based on their network topological patterns. We find that TARA is able to make accurate predictions. TARA then takes each pair of nodes that are predicted as related to be part of an alignment. Like traditional NA methods, TARA uses this alignment for the across-species transfer of functional knowledge. Clearly, TARA as currently implemented uses topological but not protein sequence information for this task. We find that TARA outperforms existing state-of-the-art NA methods that also use topological information, WAVE and SANA, and even outperforms or complements a state-of-the-art NA method that uses both topological and sequence information, PrimAlign. Hence, adding sequence information to TARA, which is our future work, is likely to further improve its performance.


title: Large deviations and entropy production in viscous fluid flows

abstract: We study the motion of a particle in a random time-dependent vector field defined by the 2D Navier-Stokes system with a noise. Under suitable non-degeneracy hypotheses we prove that the empirical measures of the trajectories of the pair (velocity field, particle) satisfy the LDP with a good rate function. Moreover, we show that the law of a unique stationary solution restricted to the particle component possesses a positive smooth density with respect to the Lebesgue measure in any finite time. This allows one to define a natural concept of the entropy production, and to show that its time average is a bounded function of the trajectory. The proofs are based on a new criterion for the validity of the level-3 LDP for Markov processes and an application of a general result on the image of probability measures under smooth maps to the laws associated with the motion of the particle.


title: Unique Continuation Properties for solutions to the Camassa-Holm equation and other non-local equations

abstract: It is shown that if $\,u(x,t)\,$ is a solution of the initial value problem for the Camassa-Holm equation which vanishes in an open set $\,\Omega\subset \mathbb R\times [0,T]$, then $\,u(x,t)=0,\,(x,t)\in\mathbb R\times [0,T]$. This result also applies to solutions of the initial periodic boundary value problems associated to the Camassa-Holm equation. The argument of proof can be placed in a general setting to extend the above results to a class of non-linear non-local 1-dimensional models which includes the Degasperis-Procesi equation.


title: The intertwined roles of particle shape and surface roughness in controlling the shear strength of a granular material

abstract: Both the shape of individual particles and their surface properties contribute to the strength of a granular material under shear. Here we show the degree to which these two aspects can be intertwined. In experiments on assemblies of 3D printed, convex lens-shaped particles, we measure the stress-strain response under repeated compressive loading and find that the aggregate's shear strength falls rapidly when compared to other particle shapes. We probe the granular material at mm-scales with X-ray computed tomography and $\mu$m-scales with high-resolution surface metrology to look for the cause of the degradation. We find that wear due to accumulated deformation smooths out the lens surfaces in a controlled and systematic manner that correlates with a significant loss of shear strength observed for the assembly as a whole. The sensitivity of lenses to changes in surface properties contrasts with results for assemblies of 3D printed tetrahedra and spheres, which under the same load cycling are found to exhibit only minor degradation in strength. This case study provides insight into the relationship between particle shape, surface wear, and the overall material response, and suggests new strategies when designing a granular material with desired evolution of properties under repeated deformation.


title: Reactor Neutrinos: Toward Oscillations

abstract: I shall sketch the history of reactor neutrino physics over five decades since the Reines-Cowan proof of neutrino existence in the late 50s, till the advent of the present era of precision reactor neutrino oscillation experiments. There are three chapters of this story: i) Exploration of possibilities of the reactor neutrinos in the 60s and 70s; ii) Looking for oscillations under the streetlamp in the 80s and 90s; iii) Exploring oscillations in detail with known (almost) $\Delta m^2_{atm}$ and $\Delta m^2_{sol}$.}


title: A Light-Weight Authentication Scheme for Air Force Internet of Things

abstract: Internet of Things (IoT) is ubiquitous because of its broad applications and the advance in communication technologies. The capabilities of IoT also enable its important role in homeland security and tactical missions, including Reconnaissance, Intelligence, Surveillance, and Target Acquisition (RISTA). IoT security becomes the most critical issue before its extensive use in military operations. While the majority of research focuses on smart IoT devices, treatments for legacy dumb network-ready devices are lacking; moreover, IoT devices deployed in a hostile environment are often required to be dumb due to the strict hardware constraints, making them highly vulnerable to cyber attacks. To mitigate the problem, we propose a light-weight authentication scheme for dumb IoT devices, in a case study of the UAV-sensor collaborative RISTA missions. Our scheme utilizes the covert channels in the physical layer for authentications and does not request conventional key deployments, key generations which may cause security risks and large overhead that a dumb sensor cannot afford. Our scheme operates on the physical layer, and thus it is highly portable and generalizable to most commercial and military communication protocols. We demonstrate the viability of our scheme by building a prototype system and conducting experiments to emulate the behaviors of UAVs and sensors in real scenarios.


title: Machine learning and chord based feature engineering for genre prediction in popular Brazilian music

abstract: Music genre can be hard to describe: many factors are involved, such as style, music technique, and historical context. Some genres even have overlapping characteristics. Looking for a better understanding of how music genres are related to musical harmonic structures, we gathered data about the music chords for thousands of popular Brazilian songs. Here, 'popular' does not only refer to the genre named MPB (Brazilian Popular Music) but to nine different genres that were considered particular to the Brazilian case. The main goals of the present work are to extract and engineer harmonically related features from chords data and to use it to classify popular Brazilian music genres towards establishing a connection between harmonic relationships and Brazilian genres. We also emphasize the generalization of the method for obtaining the data, allowing for the replication and direct extension of this work. Our final model is a combination of multiple classification trees, also known as the random forest model. We found that features extracted from harmonic elements can satisfactorily predict music genre for the Brazilian case, as well as features obtained from the Spotify API. The variables considered in this work also give an intuition about how they relate to the genres.


title: FERAtt: Facial Expression Recognition with Attention Net

abstract: We present a new end-to-end network architecture for facial expression recognition with an attention model. It focuses attention in the human face and uses a Gaussian space representation for expression recognition. We devise this architecture based on two fundamental complementary components: (1) facial image correction and attention and (2) facial expression representation and classification. The first component uses an encoder-decoder style network and a convolutional feature extractor that are pixel-wise multiplied to obtain a feature attention map. The second component is responsible for obtaining an embedded representation and classification of the facial expression. We propose a loss function that creates a Gaussian structure on the representation space. To demonstrate the proposed method, we create two larger and more comprehensive synthetic datasets using the traditional BU3DFE and CK+ facial datasets. We compared results with the PreActResNet18 baseline. Our experiments on these datasets have shown the superiority of our approach in recognizing facial expressions.


title: Fast Sequence Segmentation using Log-Linear Models

abstract: Sequence segmentation is a well-studied problem, where given a sequence of elements, an integer K, and some measure of homogeneity, the task is to split the sequence into K contiguous segments that are maximally homogeneous. A classic approach to find the optimal solution is by using a dynamic program. Unfortunately, the execution time of this program is quadratic with respect to the length of the input sequence. This makes the algorithm slow for a sequence of non-trivial length. In this paper we study segmentations whose measure of goodness is based on log-linear models, a rich family that contains many of the standard distributions. We present a theoretical result allowing us to prune many suboptimal segmentations. Using this result, we modify the standard dynamic program for one-dimensional log-linear models, and by doing so reduce the computational time. We demonstrate empirically, that this approach can significantly reduce the computational burden of finding the optimal segmentation.


title: $(g,k)$-Fermat curves: an embedding of moduli spaces

abstract: A group $H \cong {\mathbb Z}_{k}^{2g}$, where $g,k \geq 2$ are integers, of conformal automorphisms of a closed Riemann surface $S$ is called a $(g,k)$-Fermat group if it acts freely with quotient $S/H$ of genus $g$. We study some properties of these type of objects, in particular, we observe that $S$ is non-hyperelliptic and, if $k=p^{r}$, where $p>84(g-1)$ is a prime integer and $r \geq 1$, then $H$ is the unique $(g,k)$-Fermat group of $S$. Let $\Gamma$ be a co-compact torsion free Fuchsian group such that $S/H={\mathbb H}^{2}/\Gamma$. If $\Gamma_{k}$ is its normal subgroup generated by its commutators and the $k$-powers of its elements, then there is a biholomorphism between $S$ and ${\mathbb H}^{2}/\Gamma_{k}$ congugating $H$ to $\Gamma/\Gamma_{k}$. The inclusion $\Gamma_{k} < \Gamma$ induces a natural holomorphic embedding $\Theta_{k}:{\mathcal T}(\Gamma) \hookrightarrow {\mathcal T}(\Gamma_{k})$ of the corresponding Teichm\"uller spaces. Such an embedding induces a holomorphic map, at the level of their moduli spaces, $\Phi_{k}:{\mathcal M}(\Gamma) \to {\mathcal M}(\Gamma_{k})$. As a consequence of the results on $(g,k)$-Fermat groups, we provide sufficient conditions for the injectivity of $\Phi_{k}$.


title: Open data to evaluate academic researchers: an experiment with the Italian Scientific Habilitation

abstract: The need for scholarly open data is ever increasing. While there are large repositories of open access articles and free publication indexes, there are still a few examples of free citation networks and their coverage is partial. One of the results is that most of the evaluation processes based on citation counts rely on commercial citation databases. Things are changing under the pressure of the Initiative for Open Citations (I4OC), whose goal is to campaign for scholarly publishers to make their citations as totally open. This paper investigates the growth of open citations with an experiment on the Italian Scientific Habilitation, the National process for University Professor qualification which instead uses data from commercial indexes. We simulated the procedure by only using open data and explored similarities and differences with the official results. The outcomes of the experiment show that the amount of open citation data currently available is not yet enough for obtaining similar results.


title: Power Accretion in Social Systems

abstract: We consider a model of power distribution in a social system where a set of agents play a simple game on a graph: the probability of winning each round is proportional to the agent's current power, and the winner gets more power as a result. We show that, when the agents are distributed on simple 1D and 2D networks, inequality grows naturally up to a certain stationary value characterized by a clear division between a higher and a lower class of agents. High class agents are separated by one or several lower class agents which serve as a geometrical barrier preventing further flow of power between them. Moreover, we consider the effect of redistributive mechanisms, such as proportional (non-progressive) taxation. Sufficient taxation will induce a sharp transition towards a more equal society, and we argue that the critical taxation level is uniquely determined by the system geometry. Interestingly, we find that the roughness and Shannon entropy of the power distributions are a very useful complement to the standard measures of inequality, such as the Gini index and the Lorenz curve.


title: Range separation: The divide between local structures and field theories

abstract: This work presents parallel histories of the development of two modern theories of condensed matter: the theory of electron structure in quantum mechanics, and the theory of liquid structure in statistical mechanics. Comparison shows that key revelations in both are not only remarkably similar, but even follow along a common thread of controversy that marks progress from antiquity through to the present. This theme appears as a creative tension between two competing philosophies, that of short range structure (atomistic models) on the one hand, and long range structure (continuum or density functional models) on the other. The timeline and technical content are designed to build up a set of key relations as guideposts for using density functional theories together with atomistic simulation.


title: Motion Scaling Solutions for Improved Performance in High Delay Surgical Teleoperation

abstract: Robotic teleoperation brings great potential for advances within the field of surgery. The ability of a surgeon to reach patient remotely opens exciting opportunities. Early experience with telerobotic surgery has been interesting, but the clinical feasibility remains out of reach, largely due to the deleterious effects of communication delays. Teleoperation tasks are significantly impacted by unavoidable signal latency, which directly results in slower operations, less precision in movements, and increased human errors. Introducing significant changes to the surgical workflow, for example by introducing semi-automation or self-correction, present too significant a technological and ethical burden for commercial surgical robotic systems to adopt. In this paper, we present three simple and intuitive motion scaling solutions to combat teleoperated robotic systems under delay and help improve operator accuracy. Motion scaling offers potentially improved user performance and reduction in errors with minimal change to the underlying teleoperation architecture. To validate the use of motion scaling as a performance enhancer in telesurgery, we conducted a user study with 17 participants, and our results show that the proposed solutions do indeed reduce the error rate when operating under high delay.


title: Distance-based and RKHS-based Dependence Metrics in High Dimension

abstract: In this paper, we study distance covariance, Hilbert-Schmidt covariance (aka Hilbert-Schmidt independence criterion [Gretton et al. (2008)]) and related independence tests under the high dimensional scenario. We show that the sample distance/Hilbert-Schmidt covariance between two random vectors can be approximated by the sum of squared componentwise sample cross-covariances up to an asymptotically constant factor, which indicates that the distance/Hilbert-Schmidt covariance based test can only capture linear dependence in high dimension. As a consequence, the distance correlation based t-test developed by Szekely and Rizzo (2013) for independence is shown to have trivial limiting power when the two random vectors are nonlinearly dependent but component-wisely uncorrelated. This new and surprising phenomenon, which seems to be discovered for the first time, is further confirmed in our simulation study. As a remedy, we propose tests based on an aggregation of marginal sample distance/Hilbert-Schmidt covariances and show their superior power behavior against their joint counterparts in simulations. We further extend the distance correlation based t-test to those based on Hilbert-Schmidt covariance and marginal distance/Hilbert-Schmidt covariance. A novel unified approach is developed to analyze the studentized sample distance/Hilbert-Schmidt covariance as well as the studentized sample marginal distance covariance under both null and alternative hypothesis. Our theoretical and simulation results shed light on the limitation of distance/Hilbert-Schmidt covariance when used jointly in the high dimensional setting and suggest the aggregation of marginal distance/Hilbert-Schmidt covariance as a useful alternative.


title: A Note on the Paper "Optimality Conditions for Vector Optimization Problems with Difference of Convex Maps"

abstract: In this work, some counterexamples are given to refute some results reported in the paper by Guo and Li [8] (J Optim Theory Appl 162,(2014), 821-844). We correct the faulty in some of their theorems and we present alternative proofs. Moreover, we extend the definition of approximately pseudo-dissipative in the setting of metrizable topological vector spaces.


title: Automatic dimensionality selection for principal component analysis models with the ignorance score

abstract: Principal component analysis (PCA) is by far the most widespread tool for unsupervised learning with high-dimensional data sets. Its application is popularly studied for the purpose of exploratory data analysis and online process monitoring. Unfortunately, fine-tuning PCA models and particularly the number of components remains a challenging task. Today, this selection is often based on a combination of guiding principles, experience, and process understanding. Unlike the case of regression, where cross-validation of the prediction error is a widespread and trusted approach for model selection, there are no tools for PCA model selection which reach this level of acceptance. In this work, we address this challenge and evaluate the utility of the cross-validated ignorance score with both simulated and experimental data sets. Application of this method is based on the interpretation of PCA as a density model, as in probabilistic principal component analysis, and is shown to be a valuable tool to identify an optimal number of principal components.


title: The Transverse-Traceless Spin-2 Gravitational Wave Cannot Be A Standalone Observable Because It Is Acausal

abstract: We show, through an explicit calculation of the relevant Green's functions, that the transverse-traceless (TT) portion of the gravitational perturbations of Minkowski spacetime and of spatially flat cosmologies with a constant equation-of-state $w$ receive contributions from their isolated matter source(s) outside the past null cone of the observer. This implies the TT gravitational wave (GW) cannot be a standalone observable -- despite widespread (apparent) claims in the gravitational wave literature to the contrary. About a Minkowski background, all 4 of the gauge-invariant variables -- the two scalars, one vector and tensor -- play crucial roles to ensure the spatial tidal forces encoded within the gauge-invariant linearized Riemann tensor are causal. These gravitational tidal forces do not depend solely on the TT graviton but rather on the causal portion of its acceleration. However, in the far zone radiative limit, the flat spacetime `TT' graviton Green's function does reduce to the causal `tt' ones, which are the ones commonly used to compute gravitational waveforms. Similar remarks apply to the spin-1 photon; for instance, the electric field does not depend solely on the photon, but is the causal part of its velocity. As is known within the quantum theory of photons and linearized gravitons, there are obstacles to the construction of simultaneously gauge-invariant and Lorentz-covariant descriptions of these massless spin-1 and spin-2 states. Our results transparently demonstrate that the quantum operators associated with the helicity-1 photon and helicity-2 linear graviton both violate micro-causality: namely, they do not commute outside the light cone in flat and cosmological spacetimes.


title: On Generalizations of Some Distance Based Classifiers for HDLSS Data

abstract: In high dimension, low sample size (HDLSS) settings, classifiers based on Euclidean distances like the nearest neighbor classifier and the average distance classifier perform quite poorly if differences between locations of the underlying populations get masked by scale differences. To rectify this problem, several modifications of these classifiers have been proposed in the literature. However, existing methods are confined to location and scale differences only, and often fail to discriminate among populations differing outside of the first two moments. In this article, we propose some simple transformations of these classifiers resulting into improved performance even when the underlying populations have the same location and scale. We further propose a generalization of these classifiers based on the idea of grouping of variables. The high-dimensional behavior of the proposed classifiers is studied theoretically. Numerical experiments with a variety of simulated examples as well as an extensive analysis of real data sets exhibit advantages of the proposed methods.


title: Efficient Power Theft Detection for Residential Consumers Using Mean Shift Data Mining Knowledge Discovery Process

abstract: Energy theft constitutes an issue of great importance for electricity operators. The attempt to detect and reduce non-technical losses is a challenging task due to insufficient inspection methods. With the evolution of advanced metering infrastructure (AMI) in smart grids, a more complicated status quo in energy theft has emerged and many new technologies are being adopted to solve the problem. In order to identify illegal residential consumers, a computational method of analyzing and identifying electricity consumption patterns of consumers based on data mining techniques has been presented. Combining principal component analysis (PCA) with mean shift algorithm for different power theft scenarios, we can now cope with the power theft detection problem sufficiently. The overall research has shown encouraging results in residential consumers power theft detection that will help utilities to improve the reliability, security and operation of power network.


title: Probably the Best Itemsets

abstract: One of the main current challenges in itemset mining is to discover a small set of high-quality itemsets. In this paper we propose a new and general approach for measuring the quality of itemsets. The method is solidly founded in Bayesian statistics and decreases monotonically, allowing for efficient discovery of all interesting itemsets. The measure is defined by connecting statistical models and collections of itemsets. This allows us to score individual itemsets with the probability of them occuring in random models built on the data. As a concrete example of this framework we use exponential models. This class of models possesses many desirable properties. Most importantly, Occam's razor in Bayesian model selection provides a defence for the pattern explosion. As general exponential models are infeasible in practice, we use decomposable models; a large sub-class for which the measure is solvable. For the actual computation of the score we sample models from the posterior distribution using an MCMC approach. Experimentation on our method demonstrates the measure works in practice and results in interpretable and insightful itemsets for both synthetic and real-world data.


title: The Effect of Time Series Distance Functions on Functional Climate Networks

abstract: Complex network theory provides an important tool for the analysis of complex systems such as the Earth's climate. In this context, functional climate networks can be constructed using a spatiotemporal climate dataset and a suitable time series distance function. The resulting coarse-grained view on climate variability consists of representing distinct areas on the globe (i.e., grid cells) by nodes and connecting pairs of nodes that present similar time series. One fundamental concern when constructing such a functional climate network is the definition of a metric that captures the mutual similarity between time series. Here we study systematically the effect of 29 time series distance functions on functional climate network construction based on global temperature data. We observe that the distance functions previously used in the literature commonly generate very similar networks while alternative ones result in rather distinct network structures and reveal different long-distance connection patterns. These patterns are highly important for the study of climate dynamics since they generally represent pathways for the long-distance transportation of energy and can be used to forecast climate variability on subseasonal to interannual or even decadal scales. Therefore, we propose the measures studied here as alternatives for the analysis of climate variability and to further exploit their complementary capability of capturing different aspects of the underlying dynamics that may help gaining a more holistic empirical understanding of the global climate system.


title: Generalization of Kuratowski problem in linear spaces

abstract: In this short paper, Kuratowski problem will be investigated in vector space. The highest number of distinct sets that can be generated from one convex set in linear space by repeatedly applying algebraic closure and complement in any order is 8. Keywords: Kuratowski problem, algebraic interior, algebraic closure


title: The mass distribution of Galactic double neutron stars

abstract: The conventional wisdom, dating back to 2012, is that the mass distribution of Galactic double neutron stars is well-fit by a Gaussian distribution with a mean of $1.33 M_\odot$ and a width of $0.09 M_\odot$. With the recent discovery of new Galactic double neutron stars and GW170817, the first neutron star merger event to be observed with gravitational waves, it is timely to revisit this model. In order to constrain the mass distribution of double neutron stars, we perform Bayesian inference using a sample of 17 Galactic double neutron stars effectively doubling the sample used in previous studies. We expand the space of models so that the recycled neutron star need not be drawn from the same distribution as the non-recycled companion. Moreover, we consider different functional forms including uniform, single-Gaussian, and two-Gaussian distributions. While there is insufficient data to draw firm conclusions, we find positive support (a Bayes factor of 9) for the hypothesis that recycled and non-recycled neutron stars have distinct mass distributions. The most probable model---preferred with a Bayes factor of 29 over the conventional model---is one in which the recycled neutron star mass is distributed according to a two-Gaussian distribution and the non-recycled neutron star mass is distributed uniformly. We show that precise component mass measurements of $\approx 20$ double neutron stars are required in order to determine with high confidence (a Bayes factor of 150) if recycled and non-recycled neutron stars come from a common distribution. Approximately $60$ are needed in order to establish the detailed shape of the distributions.

categories: astro-ph.HE astro-ph.IM astro-ph.SR gr-qc


title: Un lemme d'analyse dont use Ibn al-Haytham en gnomonique, dioptrique et cin\'ematique c\'eleste

abstract: In several places and for different purposes, Ibn al-Haytham has studied the variations of a well-known function. The context of these applications reveals some characteristics of Ibn al-Haytham's mathematical and physical thought. An in-depth study of his mathematical proofs of monotony show how close he was from an elaborate concept of the continuous.


title: Exponential decay of correlations in the two-dimensional random field Ising model at zero temperature

abstract: We study random field Ising model on $\mathbb Z^2$ where the external field is given by i.i.d.\ Gaussian variables with mean zero and positive variance. We show that at zero temperature the effect of boundary conditions on the magnetization in a finite box decays exponentially in the distance to the boundary.


title: Including off-diagonal anisotropies in anisotropic hydrodynamics

abstract: In this paper we present a method for efficiently including the effects of off-diagonal local rest frame momentum anisotropies in leading-order anisotropic hydrodynamics. The method relies on diagonalization of the space-like block of the anisotropy tensor and allows one to reduce the necessary moments of the distribution function in the off-diagonal case to a linear combination of diagonal-anisotropy integrals. Once reduced to diagonal-anisotropy integrals, the results can be computed efficiently using techniques described previously in the literature. We present a general framework for how to accomplish this and provide examples for off-diagonal anisotropy moments entering into the energy-momentum tensor and viscous update equations which emerge when performing anisotropic pressure matching.


title: Maximum Likelihood Detection in a Four-Dimensional Stokes-Space Receiver

abstract: The maximum likelihood detection rule for a four-dimensional direct-detection optical front-end is derived. The four dimensions are two intensities and two differential phases. Three different signal processing algorithms, composed of symbol-by-symbol, sequence and successive detection, are discussed. To remedy dealing with special functions in the detection rules, an approximation for high signal-to-noise ratios (SNRs) is provided. Simulation results show that, despite the simpler structure of the successive algorithm, the resulting performance loss, in comparison with the other two algorithms, is negligible. For example, for an 8-ring/8-ary phase constellation, the complexity of detection reduces by a factor of 8, while the performance, in terms of the symbol error rate, degrades by 0.5 dB. It is shown that the high-SNR approximation is very accurate, even at low SNRs. The achievable rates for different constellations are computed and compared by the Monte Carlo method. For example, for a 4-ring/8-ary phase constellation, the achievable rate is 10 bits per channel use at an SNR of 25 dB, while by using an 8-ring/8-ary phase constellation and an error correcting code of rate 5/6, this rate is achieved at an SNR of 20 dB.


title: Consistency models in distributed systems: A survey on definitions, disciplines, challenges and applications

abstract: The replication mechanism resolves some challenges with big data such as data durability, data access, and fault tolerance. Yet, replication itself gives birth to another challenge known as the consistency in distributed systems. Scalability and availability are the challenging criteria on which the replication is based upon in distributed systems which themselves require the consistency. Consistency in distributed computing systems has been employed in three different applicable fields, such as system architecture, distributed database, and distributed systems. Consistency models based on their applicability could be sorted from strong to weak. Our goal is to propose a novel viewpoint to different consistency models utilized in the distributed systems. This research proposes two different categories of consistency models. Initially, consistency models are categorized into three groups of data-centric, client-centric and hybrid models. Each of which is then grouped into three subcategories of traditional, extended, and novel consistency models. Consequently, the concepts and procedures are expressed in mathematical terms, which are introduced in order to present our models' behavior without implementation. Moreover, we have surveyed different aspects of challenges with respect to the consistency i.e., availability, scalability, security, fault tolerance, latency, violation, and staleness, out of which the two latter i.e. violation and staleness, play the most pivotal roles in terms of consistency and trade-off balancing. Finally, the contribution extent of each of the consistency models and the growing need for them in distributed systems are investigated.


title: A simple and efficient architecture for trainable activation functions

abstract: Learning automatically the best activation function for the task is an active topic in neural network research. At the moment, despite promising results, it is still difficult to determine a method for learning an activation function that is at the same time theoretically simple and easy to implement. Moreover, most of the methods proposed so far introduce new parameters or adopt different learning techniques. In this work we propose a simple method to obtain trained activation function which adds to the neural network local subnetworks with a small amount of neurons. Experiments show that this approach could lead to better result with respect to using a pre-defined activation function, without introducing a large amount of extra parameters that need to be learned.


title: A Microwave Oscillator Based on a Single Straintronic Magneto-tunneling Junction

abstract: There is growing interest in exploring nanomagnetic devices as potential replacements for electronic devices (e.g. transistors) in digital switching circuits and systems. A special class of nanomagnetic devices are switched with electrically generated mechanical strain leading to electrical control of magnetism. Straintronic magneto-tunneling junctions (s-MTJ) belong to this category. Their soft layers are composed of two-phase multiferroics comprising a magnetostrictive layer elastically coupled to a piezoelectric layer. Here, we show that a single straintronic magneto-tunneling junction with a passive resistor can act as a microwave oscillator whose traditional implementation would have required microwave operational amplifiers, capacitors and resistors. This reduces device footprint and cost, while improving device reliability. This is an analog application of magnetic devices where magnetic interactions (interaction between the shape anisotropy, strain anisotropy, dipolar coupling field and spin transfer torque in the soft layer of the s-MTJ) are exploited to implement an oscillator with reduced footprint.


title: Penalized linear regression with high-dimensional pairwise screening

abstract: In variable selection, most existing screening methods focus on marginal effects and ignore dependence between covariates. To improve the performance of selection, we incorporate pairwise effects in covariates for screening and penalization. We achieve this by studying the asymptotic distribution of the maximal absolute pairwise sample correlation among independent covariates. The novelty of the theory is in that the convergence is with respect to the dimensionality $p$, and is uniform with respect to the sample size $n$. Moreover, we obtain an upper bound for the maximal pairwise R squared when regressing the response onto two different covariates. Based on these extreme value results, we propose a screening procedure to detect covariates pairs that are potentially correlated and associated with the response. We further combine the pairwise screening with Sure Independence Screening and develop a new regularized variable selection procedure. Numerical studies show that our method is very competitive in terms of both prediction accuracy and variable selection accuracy.


title: Linear Dynamics & Control of Brain Networks

abstract: The brain is an intricately structured organ responsible for the rich emergent dynamics that support the complex cognitive functions we enjoy as humans. With around $10^{11}$ neurons and $10^{15}$ synapses, understanding how the human brain works has proven to be a daunting endeavor, requiring concerted collaboration across traditional disciplinary boundaries. In some cases, that collaboration has occurred between experimentalists and technicians, who offer new physical tools to measure and manipulate neural function. In other contexts, that collaboration has occurred between experimentalists and theorists, who offer new conceptual tools to explain existing data and inform new directions for empirical research. In this chapter, we offer an example of the latter. Specifically, we focus on the simple but powerful framework of linear systems theory as a useful tool both for capturing biophysically relevant parameters of neural activity and connectivity, and for analytical and numerical study. We begin with a brief overview of state-space representations and linearization of neural models for non-linear dynamical systems. We then derive core concepts in the theory of linear systems such as the impulse and controlled responses to external stimuli, achieving desired state transitions, controllability, and minimum energy control. Afterwards, we discuss recent advances in the application of linear systems theory to structural and functional brain data across multiple spatial and temporal scales, along with methodological considerations and limitations. We close with a brief discussion of open frontiers and our vision for the future.


title: Preserve or retreat? Willingness-to-pay for Coastline Protection in New South Wales

abstract: Coastal erosion is a global and pervasive phenomenon that predicates a need for a strategic approach to the future management of coastal values and assets (both built and natural), should we invest in protective structures like seawalls that aim to preserve specific coastal features, or allow natural coastline retreat to preserve sandy beaches and other coastal ecosystems. Determining the most suitable management approach in a specific context requires a better understanding of the full suite of economic values the populations holds for coastal assets, including non-market values. In this study, we characterise New South Wales residents willingness to pay to maintain sandy beaches (width and length). We use an innovative application of a Latent Class Binary Logit model to deal with Yea-sayers and Nay-sayers, as well as revealing the latent heterogeneity among sample members. We find that 65% of the population would be willing to pay some amount of levy, dependent on the policy setting. In most cases, there is no effect of degree of beach deterioration characterised as loss of width and length of sandy beaches of between 5% and 100% on respondents willingness to pay for a management levy. This suggests that respondents who agreed to pay a management levy were motivated to preserve sandy beaches in their current state irrespective of the severity of sand loss likely to occur as a result of coastal erosion. Willingness to pay also varies according to beach type (amongst Iconic, Main, Bay and Surf beaches) a finding that can assist with spatial prioritisation of coastal management. Not recognizing the presence of nay-sayers in the data or recognizing them but eliminating them from the estimation will result in biased WTP results and, consequently, biased policy propositions by coastal managers.


title: The asymptotically sharp geometric rigidity interpolation estimate in thin bi-Lipschitz domains

abstract: This work is part of a program of development of asymptotically sharp geometric rigidity estimates for thin domains. A thin domain in three dimensional Euclidean space is roughly a small neighborhood of regular enough two dimensional compact surface. We prove an asymptotically sharp geometric rigidity interpolation inequality for thin domains with little regularity. In contrast to that celebrated Friesecke-James-M\"uller rigidity estimate [\textit{Comm. Pure Appl. Math.,} 55(11):1461-1506, 2002] for plates, our estimate holds for any proper rotations $\BR\in SO(3).$ Namely, the estimate bounds the $L^p$ distance of the gradient of any $\By\in W^{1,p}(\Omega,\mathbb R^3)$ field from any constant proper rotation $\BR\in SO(3)$, in terms of the average $L^p$ distance (nonlinear strain) of the gradient $\nabla\By$ from the rotation group $SO(3)$, and the average $L^p$ distance of the field itself from the set of rigid motions corresponding to the rotation $\BR$. There are several remarkable facts about the estimate: 1. The constants in the estimate are sharp in terms of the domain thickness scaling for any thin domains with the required regularity. 2. In the special cases when the domain has positive or negative Gaussian curvature, the inequality reduces the problem of estimating the gradient $\nabla\By$ in terms of the prototypical nonlinear strain $\int_\Omega\mathrm{dist}^p(\nabla\By(x),SO(3))dx$ to the easier problem of estimating only the vector field $\By$ in terms of the nonlinear strain without any loss in the constant scalings as the Ans\"atze suggest. The later will be a geometric rigidity Korn-Poincar\'e type estimate. This passage is major progress in the thin domain rigidity problem. 3. For the borderline energy scaling (bending-to-stretching), the estimate implies improved strong compactness on the vector fields for free.


title: Subnanosecond Fluctuations in Low-Barrier Nanomagnets

abstract: Fast magnetic fluctuations due to thermal torques have useful technological functionality ranging from cryptography to probabilistic computing. The characteristic time of fluctuations in typical uniaxial anisotropy magnets studied so far is bounded from below by the well-known energy relaxation mechanism. This time scales as $\alpha^{-1}$, where $\alpha$ parameterizes the strength of dissipative processes. Here, we theoretically analyze the fluctuating dynamics in easy-plane and antiferromagnetically coupled nanomagnets. We find in such magnets, the dynamics are strongly influenced by fluctuating intrinsic fields, which give rise to an additional dephasing-type mechanism for washing out correlations. In particular, we establish two time scales for characterizing fluctuations (i) the average time for a nanomagnet to reverse|which for the experimentally relevant regime of low damping is governed primarily by dephasing and becomes independent of $\alpha$, (ii) the time scale for memory loss of a single nanomagnet|which scales as $\alpha^{-1/3}$ and is governed by a combination of energy dissipation and dephasing mechanism. For typical experimentally accessible values of intrinsic fields, the resultant thermal-fluctuation rate is increased by multiple orders of magnitude when compared with the bound set solely by the energy relaxation mechanism in uniaxial magnets. This could lead to higher operating speeds of emerging devices exploiting magnetic fluctuations.


title: Quasi-optimal and pressure robust discretizations of the Stokes equations by new augmented Lagrangian formulations

abstract: We approximate the solution of the stationary Stokes equations with various conforming and nonconforming inf-sup stable pairs of finite element spaces on simplicial meshes. Based on each pair, we design a discretization that is quasi-optimal and pressure robust, in the sense that the velocity $H^1$-error is proportional to the best $H^1$-error to the analytical velocity. This shows that such a property can be achieved without using conforming and divergence-free pairs. We bound also the pressure $L^2$-error, only in terms of the best approximation errors to the analytical velocity and the analytical pressure. Our construction can be summarized as follows. First, a linear operator acts on discrete velocity test functions, before the application of the load functional, and maps the discrete kernel into the analytical one. Second, in order to enforce consistency, we employ a new augmented Lagrangian formulation, inspired by Discontinuous Galerkin methods.


title: Routing in Networks on Chip with Multiplicative Circulant Topology

abstract: The development of multi-core processor systems is a demanded branch of science and technology. The appearance of processors with dozens and hundreds of cores poses to the developers the question of choosing the optimal topology capable to provide efficient routing in a network with a large number of nodes. In this paper, we consider the possibility of using multiplicative circulants as a topology for networks-on-chip. A specialized routing algorithm for networks with multiplicative circulant topology, taking into account topology features and having a high scalability, has been developed.


title: Prediction of glass formability from liquid properties

abstract: Glass formation is one of the most interesting phenomena in the condensed matter field. Considerable effort has gone into understanding and predicting the glass formability. However, the previous prediction requires the glass first made before the prediction can be performed. Here, we propose a new prediction formula using liquid properties only. Moreover, we demonstrated that the similarity between liquid and crystalline structure plays an important role in determine the glass formability. Previously, only the kinetics of nucleation and growth processes have been considered.


title: Bayesian Model Selection with Graph Structured Sparsity

abstract: We propose a general algorithmic framework for Bayesian model selection. A spike-and-slab Laplacian prior is introduced to model the underlying structural assumption. Using the notion of effective resistance, we derive an EM-type algorithm with closed-form iterations to efficiently explore possible candidates for Bayesian model selection. The deterministic nature of the proposed algorithm makes it more scalable to large-scale and high-dimensional data sets compared with existing stochastic search algorithms. When applied to sparse linear regression, our framework recovers the EMVS algorithm [Rockova and George, 2014] as a special case. We also discuss extensions of our framework using tools from graph algebra to incorporate complex Bayesian models such as biclustering and submatrix localization. Extensive simulation studies and real data applications are conducted to demonstrate the superior performance of our methods over its frequentist competitors such as $\ell_0$ or $\ell_1$ penalization.


title: PASTA: A Parallel Sparse Tensor Algorithm Benchmark Suite

abstract: Tensor methods have gained increasingly attention from various applications, including machine learning, quantum chemistry, healthcare analytics, social network analysis, data mining, and signal processing, to name a few. Sparse tensors and their algorithms become critical to further improve the performance of these methods and enhance the interpretability of their output. This work presents a sparse tensor algorithm benchmark suite (PASTA) for single- and multi-core CPUs. To the best of our knowledge, this is the first benchmark suite for sparse tensor world. PASTA targets on: 1) helping application users to evaluate different computer systems using its representative computational workloads; 2) providing insights to better utilize existed computer architecture and systems and inspiration for the future design. This benchmark suite is publicly released https://gitlab.com/tensorworld/pasta.


title: Reciprocal Space Imaging of Ionic Correlations in Intercalation Compounds

abstract: The intercalation of alkali ions into layered materials has played an essential role in battery technology since the development of the first lithium-ion electrodes. Coulomb repulsion between the intercalants leads to ordering of the intercalant sublattice, which hinders ionic diffusion and impacts battery performance. While conventional diffraction can identify the long-range order that can occur at discrete intercalant concentrations during the charging cycle, it cannot determine short-range order at other concentrations that also disrupt ionic mobility. In this article, we show that the use of real-space transforms of single crystal diffuse scattering, measured with high-energy synchrotron x-rays, allows a model-independent measurement of the temperature dependence of the length scale of ionic correlations along each of the crystallographic axes in a sodium-intercalated V$_2$O$_5$. The techniques described here provide a new way of probing the evolution of structural ordering in crystalline materials.


title: A Differentiable Augmented Lagrangian Method for Bilevel Nonlinear Optimization

abstract: Many problems in modern robotics can be addressed by modeling them as bilevel optimization problems. In this work, we leverage augmented Lagrangian methods and recent advances in automatic differentiation to develop a general-purpose nonlinear optimization solver that is well suited to bilevel optimization. We then demonstrate the validity and scalability of our algorithm with two representative robotic problems, namely robust control and parameter estimation for a system involving contact. We stress the general nature of the algorithm and its potential relevance to many other problems in robotics.


title: Active Area Coverage from Equilibrium

abstract: This paper develops a method for robots to integrate stability into actively seeking out informative measurements through coverage. We derive a controller using hybrid systems theory that allows us to consider safe equilibrium policies during active data collection. We show that our method is able to maintain Lyapunov attractiveness while still actively seeking out data. Using incremental sparse Gaussian processes, we define distributions which allow a robot to actively seek out informative measurements. We illustrate our methods for shape estimation using a cart double pendulum, dynamic model learning of a hovering quadrotor, and generating galloping gaits starting from stationary equilibrium by learning a dynamics model for the half-cheetah system from the Roboschool environment.


title: Exchangeable and Sampling Consistent Distributions on Rooted Binary Trees

abstract: We introduce a notion of finite sampling consistency for phylogenetic trees and show that the set of finitely sampling consistent and exchangeable distributions on n leaf phylogenetic trees is a polytope. We use this polytope to show that the set of all exchangeable and infinite sampling consistent distributions on 4 leaf phylogenetic trees is exactly Aldous' beta-splitting model and give a description of some of the vertices for the polytope of distributions on 5 leaves. We also introduce a new semialgebraic set of exchangeable and sampling consistent models we call the multinomial model and use it to characterize the set of exchangeable and sampling consistent distributions.


title: An Approach to Track Reading Progression Using Eye-Gaze Fixation Points

abstract: In this paper, we consider the problem of tracking the eye-gaze of individuals while they engage in reading. Particularly, we develop ways to accurately track the line being read by an individual using commercially available eye tracking devices. Such an approach will enable futuristic functionalities such as comprehension evaluation, interest level detection, and user-assisting applications like hands-free navigation and automatic scrolling. Existing commercial eye trackers provide an estimated location of the eye-gaze fixations every few milliseconds. However, this estimated data is found to be very noisy. As such, commercial eye-trackers are unable to accurately track lines while reading. In this paper we propose several statistical models to bridge the commercial gaze tracker outputs and eye-gaze patterns while reading. We then employ hidden Markov models to parametrize these statistical models and to accurately detect the line being read. The proposed approach is shown to yield an improvement of over 20% in line detection accuracy.


title: Teaching Physics at School and Colleges

abstract: Invited contribution to the inaugural issue of {\em Physics Educator.}


title: Similarity renormalization group evolution of hypernuclear Hamiltonians

abstract: Unitary transformations of a Hamiltonian generally induce interaction terms beyond the particle rank present in the untransformed Hamiltonian that have to be captured and included in a many-body calculation. In systems with strangeness such as hypernuclei, the three-body terms induced by the hyperon-nucleon interaction are strong, so their inclusion is crucial. We present in detail a procedure for computing hyperon-nucleon-nucleon interaction terms that are induced during a similarity renormalization group (SRG) flow. The SRG is carried out in a basis spanned by antisymmetric harmonic-oscillator states with respect to three-body Jacobi coordinates. We discuss basis construction, antisymmetrization, numerical evaluation of the flow equations, and separation of the genuine three-body terms. We then use the hypernuclear no-core shell model with SRG-evolved Hamiltonians, addressing the sensitivity of hypernuclear states and hyperon separation energies to changes in the nucleonic Hamiltonian by example of $^7_\Lambda$Li, $^9_\Lambda$Be, $^{11}_\Lambda$B, $^{13}_\Lambda$C, and the hyper-helium chain. We also present a survey of the hyper-hydrogen chain, exploring the structure of hypernuclei with extreme neutron-proton asymmetries.


title: A tracially AF algebra which is not $\mathcal Z$-absorbing

abstract: We show that there is a simple separable unital (non-nuclear) tracially AF algebra $A$ which does not absorb the Jiang-Su algebra $\mathcal Z$ tensorially, i.e., $A \ncong A\otimes\mathcal Z$.


title: Architecture Compression

abstract: In this paper we propose a novel approach to model compression termed Architecture Compression. Instead of operating on the weight or filter space of the network like classical model compression methods, our approach operates on the architecture space. A 1-D CNN encoder-decoder is trained to learn a mapping from discrete architecture space to a continuous embedding and back. Additionally, this embedding is jointly trained to regress accuracy and parameter count in order to incorporate information about the architecture's effectiveness on the dataset. During the compression phase, we first encode the network and then perform gradient descent in continuous space to optimize a compression objective function that maximizes accuracy and minimizes parameter count. The final continuous feature is then mapped to a discrete architecture using the decoder. We demonstrate the merits of this approach on visual recognition tasks such as CIFAR-10, CIFAR-100, Fashion-MNIST and SVHN and achieve a greater than 20x compression on CIFAR-10.


title: Censored Quantile Regression Forests

abstract: Random forests are powerful non-parametric regression method but are severely limited in their usage in the presence of randomly censored observations, and naively applied can exhibit poor predictive performance due to the incurred biases. Based on a local adaptive representation of random forests, we develop its regression adjustment for randomly censored regression quantile models. Regression adjustment is based on new estimating equations that adapt to censoring and lead to quantile score whenever the data do not exhibit censoring. The proposed procedure named censored quantile regression forest, allows us to estimate quantiles of time-to-event without any parametric modeling assumption. We establish its consistency under mild model specifications. Numerical studies showcase a clear advantage of the proposed procedure.


title: Sub-Eddington super-massive black hole activity in Fornax early type galaxies

abstract: We characterize the incidence and intensity of low-level super-massive black hole activity within the Fornax cluster, through X-ray observations of the nuclei of 29 quiescent early-type galaxies. Using the \textit{Chandra X-ray Telescope}, we target 17 galaxies from the HST Fornax Cluster Survey, down to a uniform (3$\sigma$) limiting X-ray luminosity threshold of $5\cdot10^{38}$ergs$^{-1}$, which we combine with deeper, archival observations for an additional 12 galaxies. A nuclear X-ray point-source is detected in 11 out of 29 targets. After accounting for the low mass X-ray binary contamination to the nuclear X-ray signal, the X-ray active fraction is measured at $26.6\% \pm 9.6\%$. The results from this analysis are compared to similar investigations targeting quiescent early types in the Virgo cluster, as well as the field. After correcting for the different mass distributions, the measured Fornax active fraction is less than the field fraction, at more than 3$\sigma$, confirming that the funneling of gas to the nuclear regions of cluster members is inhibited compared to those galaxies in the field. At the same time, we find no statistically significant difference between Fornax and Virgo galaxies, with only marginal evidence for a lower active fraction in Fornax (1 $\sigma$); if real, owing to Fornax's higher galaxy number density, this could indicate that galaxy-galaxy interactions are more effective at gas removal than galaxy-gas effects.


title: Continuity equation and vacuum regions in compressible flows

abstract: We investigate the creation and properties of eventual vacuum regions in the weak solutions of the continuity equation, in general, and in the weak solutions of compressible Navier--Stokes equations, in particular. The main results are based on the analysis of renormalized solutions to the continuity and pure transport equations and their inter-relations which are of independent interest.


title: Simulation of the Magnetization Dynamics of a Single Domain BiFeO$_3$ Thin Film

abstract: The switching dynamics of a single-domain BiFeO$_3$ thin films is investigated through combining the dynamics of polarization and Neel vector. The evolution of the ferroelectric polarization is described by the Landau-Khalatnikov (LK) equation, and the Landau-Lifshitz-Gilbert (LLG) equations for spins in two sublattices to model the time evolution of the antiferromagnetic order (Neel vector) in a G-type antiferromagnet. This work theoretically demonstrates that due to the rotation of the magnetic hard axis following the polarization reversal, the Neel vector can be switched by 180 degrees, while the weak magnetization can remain unchanged. The simulation results are consistent with the ab initio calculation, where the Neel vector rotates during polarization rotation, and also match our calculation of the dynamics of order parameter using Landau-Ginzburg theory. We also find that the switching time of the Neel vector is determined by the speed polarization switching and is predicted to be as short as 30 ps.


title: Systematization of Vulnerability Discovery Knowledge: Review Protocol

abstract: In this report, we describe the review protocol that will guide the systematic review of the literature in metrics-based discovery of vulnerabilities. The protocol have been developed in adherence with the guidelines for performing Systematic Literature Reviews in Software Engineering prescribed by Kitchenham and Charters.


title: Foliation-based approach to quantum gravity and applications to astrophysics

abstract: The recently proposed Holography-inspired approach to quantum gravity is reviewed and expanded. The approach is based on the foliation of the background spacetime and reduction of the offshell states to the physical states. Careful attention is paid to the boundary conditions. It is noted that the outstanding problems such as the cosmological constant problem and black hole information can be tackled from the common thread of the quantized gravity. One-loop renormalization of the coupling constants and the beta function analysis are illustrated. Active galactic nuclei and gravitational waves are discussed as the potential applications of the present quantization scheme to astrophysics.