title: Minimum variance estimation of statistical anisotropy via galaxy survey

abstract: We consider the benefits of measuring cosmic statistical anisotropy from redshift-space correlators of the galaxy number density fluctuation and the peculiar velocity field without adopting the plane-parallel (PP) approximation. Since the correlators are decomposed using the general tripolar spherical harmonic (TripoSH) basis, we can deal with wide-angle contributions untreatable by the PP approximation, and at the same time, target anisotropic signatures can be cleanly extracted. We, for the first time, compute the covariance of the TripoSH decomposition coefficient and the Fisher matrix to forecast the detectability of statistical anisotropy. The resultant expression of the covariance is free from nontrivial mixings between each multipole moment caused by the PP approximation and hence the detectability is fully optimized. Compared with the analysis under the PP approximation, the superiority in detectability is always confirmed, and it is highlighted, especially in the cases that the shot noise level is large and that target statistical anisotropy has a blue-tilted shape in Fourier space. The application of the TripoSH-based analysis to forthcoming all-sky survey data could result in constraints on anisotropy comparable to or tighter than the current cosmic microwave background ones.


title: Minimal complete sets for two pseudoscalar meson photoproduction

abstract: For photoproduction reactions with final states consisting of two pseudoscalar mesons and a spin-1/2 baryon, 8 complex amplitudes need to be determined uniquely. A modified version of Moravcsik's theorem is employed for these reactions, resulting in slightly over-complete sets of polarization observables that are able to determine the amplitudes uniquely. Further steps were taken to reduce the found sets to minimal complete sets. As a final result, multiple minimal complete sets without any remaining ambiguities are presented for the first time. These sets consist of $2N=16$ observables, containing only one triple polarization observable.


title: Charge radii of the nucleon from its flavor dependent Dirac form factors

abstract: We have determined the proton and the neutron charge radii from a global analysis of the proton and the neutron elastic form factors, after first performing a flavor decomposition of these form factors under charge symmetry in the light cone frame formulation. We then extracted the transverse mean-square radii of the flavor dependent quark distributions. In turn, these are related in a model-independent way to the proton and neutron charge radii but allow us to take into account motion effects of the recoiling nucleon for data at finite but high momentum transfer. In the proton case we find $\langle r_p \rangle = 0.852 \pm0.002_{\rm (stat.)} \pm0.009_{\rm (syst.)}~({\rm fm})$, consistent with the proton charge radius obtained from muonic hydrogen spectroscopy \cite{pohl:2010,antog2013}. The current method improves on the precision of the $\langle r_p \rangle$ extraction based on the form factor measurements. Furthermore, we find no discrepancy in the $\langle r_p \rangle$ determination among the different electron scattering measurements, all of which, utilizing the current method of extraction, result in a value that is consistent with the smallest $\langle r_p \rangle$ extraction from the electron scattering measurements \cite{Xiong:2019umf}. Concerning the neutron case, past results relied solely on the neutron-electron scattering length measurements, which suffer from an underestimation of underlying systematic uncertainties inherent to the extraction technique. Utilizing the present method we have performed the first extraction of the neutron charge radius based on nucleon form factor data, and we find $\langle r_n^2 \rangle = -0.122 \pm0.004_{\rm (stat.)} \pm0.010_{\rm (syst.)}~({\rm fm}^2)$.


title: Controlling rotation in the molecular-frame with an optical centrifuge

abstract: We computationally demonstrate a new method for coherently controlling the rotation-axis direction in asymmetric top molecules with an optical centrifuge. Appropriately chosen electric-field strengths and the centrifuge's acceleration rate allow to generate a nearly arbitrary rotational wavepacket. For D$_2$S and 2H-imidazole (C$_3$H$_4$N$_2$) we created wavepackets at large values of the rotational quantum number $J$ with the desired projections of the total angular momentum onto two of the molecules' principal axes of inertia. One application of the new method is three-dimensional alignment with a molecular axis aligned along the laser's wave vector, which is important for the three-dimensional imaging of molecules yet not accessible in standard approaches. The simultaneous orientation of the angular momentum in the laboratory frame and in the molecular frame could also be used in robust control of scattering experiments.


title: Towards forecast techniques for business analysts of large commercial data sets using matrix factorization methods

abstract: This research article suggests that there are significant benefits in exposing demand planners to forecasting methods using matrix completion techniques. This study aims to contribute to a better understanding of the field of forecasting with multivariate time series prediction by focusing on the dimension of large commercial data sets with hierarchies. This research highlights that there has neither been sufficient academic research in this sub-field nor dissemination among practitioners in the business sector. This study seeks to innovate by presenting a matrix completion method for short-term demand forecast of time series data on relevant commercial problems. Albeit computing intensive, this method outperforms the state of the art while remaining accessible to business users. The object of research is matrix completion for time series in a big data context within the industry. The subject of the research is forecasting product demand using techniques for multivariate hierarchical time series prediction that are both precise and accessible to non-technical business experts. Apart from a methodological innovation, this research seeks to introduce practitioners to novel methods for hierarchical multivariate time series prediction. The research outcome is of interest for organizations requiring precise forecasts yet lacking the appropriate human capital to develop them.


title: Modified Gravitation Theory (MOG) and the aLIGO GW190521 Gravitational Wave Event

abstract: A consequence of adopting a modified gravitational theory (MOG) for the aLIGO GW190521 gravitational wave detection involving binary black hole sources is to fit the aLIGO strain and chirp data with lower mass, compact coalescing binary systems such as neutron star-neutron star (NS-NS), black hole - neutron star (BH-NS), and black hole-black hole (BH-BH) systems. In MOG BH - BH component masses can be smaller than the component masses $m_1=85M_\odot$ and $m_2=66M_\odot$ inferred from the aLIGO GW190521 gravitational wave event. This reduces the mass of the final remnant mass $M_f=150M_\odot$ and allows the primary, secondary and final remnant masses of the black holes to be formed by conventional stellar collapse models.


title: Baryogenesis through baryon capture by black holes

abstract: A novel mechanism of cosmological baryogenesis through baryon capture by primordial black holes is suggested. In contrast to the conventional scenarios it does not demand non-conservation of baryonic number in particle physics and can proceed in thermal equilibrium. For implementation of this mechanism a heavy superweakly interacting particle a with non-zero baryon number is necessary.


title: Integrated Benchmarking and Design for Reproducible and Accessible Evaluation of Robotic Agents

abstract: As robotics matures and increases in complexity, it is more necessary than ever that robot autonomy research be reproducible. Compared to other sciences, there are specific challenges to benchmarking autonomy, such as the complexity of the software stacks, the variability of the hardware and the reliance on data-driven techniques, amongst others. In this paper, we describe a new concept for reproducible robotics research that integrates development and benchmarking, so that reproducibility is obtained "by design" from the beginning of the research/development processes. We first provide the overall conceptual objectives to achieve this goal and then a concrete instance that we have built: the DUCKIENet. One of the central components of this setup is the Duckietown Autolab, a remotely accessible standardized setup that is itself also relatively low-cost and reproducible. When evaluating agents, careful definition of interfaces allows users to choose among local versus remote evaluation using simulation, logs, or remote automated hardware setups. We validate the system by analyzing the repeatability of experiments conducted using the infrastructure and show that there is low variance across different robot hardware and across different remote labs.


abstract: Evidence for Higgs boson decay to a pair of muons is presented. This result combines searches in four exclusive categories targeting the production of the Higgs boson via gluon fusion, via vector boson fusion, in association with a vector boson, and in association with a top quark-antiquark pair. The analysis is performed using proton-proton collision data at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 137 fb$^{-1}$, recorded by the CMS experiment at the CERN LHC. An excess of events over the background expectation is observed in data with a significance of 3.0 standard deviations, where the expectation for the standard model (SM) Higgs boson with mass of 125.38 GeV is 2.5. The combination of this result with that from data recorded at $\sqrt{s} =$ 7 and 8 TeV, corresponding to integrated luminosities of 5.1 and 19.7 fb$^{-1}$, respectively, increases both the expected and observed significances by 1%. The measured signal strength, relative to the SM prediction, is 1.19 $^{+0.40}_{-0.39}$ (stat) $^{+0.15}_{-0.14}$ (syst). This result constitutes the first evidence for the decay of the Higgs boson to second generation fermions and is the most precise measurement of the Higgs boson coupling to muons reported to date.


title: A Random Antenna Subset Selection Jamming Method against Multistatic Radar System

abstract: Multistatic radar system (MSRS) is considered an effective scheme to suppress mainlobe jamming, since it has higher spatial resolution enabling jamming cancellation from spatial domain. To develop electronic countermeasures against MSRS, a random array subset selection (RASS)jamming method is proposed in this paper. In the RASS jammer, elements of the array antenna are activated randomly, leading to stable mainlobe and random sidelobes, different from the traditional jammer that applies the complete antenna array enjoying constant mainlobe and sidelobes. We study the covariance matrix of jamming signals received by radars, and derive its rank, revealing that the covariance matrix is of full rank. We also calculate the output jamming to signal and noise ratio (JSNR) after the subspace-based jamming suppression methods used in MSRS under the proposed jamming method, which demonstrates that the full rank property invalidates such suppression methods. Numerical results verify our analytical deduction and exhibit the improved countermeasure performance of our proposed RASS jamming method compared to the traditional one.


title: Plant Diseases recognition on images using Convolutional Neural Networks: A Systematic Review

abstract: Plant diseases are considered one of the main factors influencing food production and minimize losses in production, and it is essential that crop diseases have fast detection and recognition. The recent expansion of deep learning methods has found its application in plant disease detection, offering a robust tool with highly accurate results. In this context, this work presents a systematic review of the literature that aims to identify the state of the art of the use of convolutional neural networks(CNN) in the process of identification and classification of plant diseases, delimiting trends, and indicating gaps. In this sense, we present 121 papers selected in the last ten years with different approaches to treat aspects related to disease detection, characteristics of the data set, the crops and pathogens investigated. From the results of the systematic review, it is possible to understand the innovative trends regarding the use of CNNs in the identification of plant diseases and to identify the gaps that need the attention of the research community.


title: Effect of Coronal Loop Structure on Wave Heating by Phase Mixing

abstract: The mechanism behind coronal heating still elude direct observation and modelling of viable theoretical processes and the subsequent effect on coronal structures is one of the key tools available to assess possible heating mechanisms. Wave-heating via phase-mixing of Magnetohydrodynamics (MHD) transverse waves has been proposed as a possible way to convert magnetic energy into thermal energy but increasingly, MHD models suggest this is not a sufficiently efficient mechanism. We model heating by phase-mixing of transverse MHD waves in various configurations, to investigate whether certain circumstances can enhance the heating sufficiently to sustain the million degree solar corona and to assess the impact of the propagation and phase-mixing of transverse MHD waves on the structure of the boundary shell of coronal loops. We use 3D MHD simulations of a pre-existing density enhancement in magnetised medium and a boundary driver to trigger the propagation of transverse waves with the same power spectrum as measured by the COmP. We consider different density structures, boundary conditions at the non-drive footpoint, characteristics of the driver, and different forms of magnetic resistivity. We find that different initial density structures affect the evolution of the boundary shell and some driver configurations enhance the heating generated from the dissipation of the MHD waves. In particular, drivers coherent on a larger spatial scale and higher dissipation coefficients generate significant heating, although it is still insufficient. We conclude that while phase-mixing of transverse MHD waves is unlikely to sustain the thermal structure of the corona, there are configurations that allow for an enhanced efficiency of this mechanism. We provide possible signatures to identify the presence of such configurations, such as the location of where the heating is deposited along the coronal loop


title: The molecular nature of some exotic hadrons

abstract: The exciting discovery by LHCb of the $P_c(4312)^+$ and $P_c(4450)^+$ pentaquarks, or the suggestion of a tetraquark nature for the $Z_c(3900)$ state seen at BESIII and Belle, have triggered a lot of activity in the field of hadron physics, with new experiments planned for searching other exotic mesons and baryons, and many theoretical developments trying to disentangle the true multiquark nature from their possible molecular origin. After a brief review of the present status of these searches, this paper focusses on recently seen or yet to be discovered exotic heavy baryons that may emerge from a conveniently unitarized meson-baryon interaction model in coupled channels. In particular, we will show how interferences between the different coupled-channel amplitudes of the model may reveal the existence of a $N^*$ resonance around 2 GeV having a meson-baryon quasi-bound state nature. We also discuss the possible interpretation of some of the $\Omega_c$ states recently discovered at LHCb as being hadron molecules. The model also predicts the existence of doubly-charmed quasibound meson-baryon $\Xi_{cc}$ states, which would be excited states of the ground-state $\Xi_{cc}(3621)$ MeV, whose mass has only been recently established. Extensions of these results to the bottom sector will also be presented.


title: The Magellanic Corona and the formation of the Magellanic Stream

abstract: The dominant gaseous structure in the Galactic halo is the Magellanic Stream, an extended network of neutral and ionized filaments surrounding the Large and Small Magellanic Clouds (LMC/SMC), the two most massive satellite galaxies of the Milky Way. Recent observations indicate that the Clouds are on their first passage around our Galaxy, the Stream is made up of gas stripped from both the LMC and the SMC, and the majority of this gas is ionized. While it has long been suspected that tidal forces and ram-pressure stripping contributed to the Stream's formation, a full understanding of its origins has defied modelers for decades. Several recent developments, including the discovery of dwarf galaxies associated with the Magellanic Group, the high mass of the LMC, the detection of highly ionized gas toward stars in the LMC and the predictions of cosmological simulations all support the existence of a halo of warm ionized gas around the LMC at a temperature of $\sim5\times10^{5}\;\mathrm{K}$. Here we show that by including this "Magellanic Corona" in hydrodynamic simulations of the Magellanic Clouds falling onto the Galaxy, we can simultaneously reproduce the Stream and its Leading Arm. Our simulations explain the Stream's filamentary structure, spatial extent, radial velocity gradient, and total ionized gas mass. We predict that the Magellanic Corona will be unambiguously observable via high-ionization absorption lines in the ultraviolet spectra of background quasars lying near the LMC. This prediction is directly testable with the Cosmic Origins Spectrograph on the Hubble Space Telescope.


title: Viscous shock solutions to the stochastic Burgers equation

abstract: We define a notion of a viscous shock solution of the stochastic Burgers equation that connects "top" and "bottom" spatially stationary solutions of the same equation. Such shocks generally travel in space, but we show that they admit time-invariant measures when viewed in their own reference frames. Under such a measure, the viscous shock is a deterministic function of the bottom and top solutions and the shock location. However, the measure of the bottom and top solutions must be tilted to account for the change of reference frame. We also show a convergence result to these stationary shock solutions from solutions initially connecting two constants, as time goes to infinity.


title: A simple solver for the two-fluid plasma model based on PseudoSpectral Time-Domain algorithm

abstract: We present a solver of 3D two-fluid plasma model for the simulation of short-pulse laser interactions with plasma. This solver resolves the equations of the two-fluid plasma model with ideal gas closure. We also include the Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral Time-Domain (PSTD) method to solve Maxwell's curl equations. We use a Strang splitting to integrate Euler equations with source term: while Euler equations are solved with a composite scheme mixing Lax-Friedrichs and Lax-Wendroff schemes, the source term is integrated with a fourth-order Runge-Kutta scheme. This two-fluid plasma model solver is simple to implement because it only relies on finite difference schemes and Fast Fourier Transforms. It does not require spatially staggered grids. The solver was tested against several well-known problems of plasma physics. Numerical simulations gave results in excellent agreement with analytical solutions or with previous results from the literature.


title: Infinitely many embedded eigenvalues for the Neumann-Poincar\'e operator in 3D

abstract: This article constructs a surface whose Neumann-Poincar\'e (NP) integral operator has infinitely many eigenvalues embedded in its essential spectrum. The surface is a sphere perturbed by smoothly attaching a conical singularity, which imparts essential spectrum. Rotational symmetry allows a decomposition of the operator into Fourier components. Eigenvalues of infinitely many Fourier components are constructed so that they lie within the essential spectrum of other Fourier components and thus within the essential spectrum of the full NP operator. The proof requires the perturbation to be sufficiently small, with controlled curvature, and the conical singularity to be sufficiently flat.


title: A Generalized Online Algorithm for Translation and Scale Invariant Prediction with Expert Advice

abstract: In this work, we aim to create a completely online algorithmic framework for prediction with expert advice that is translation-free and scale-free of the expert losses. Our goal is to create a generalized algorithm that is suitable for use in a wide variety of applications. For this purpose, we study the expected regret of our algorithm against a generic competition class in the sequential prediction by expert advice problem, where the expected regret measures the difference between the losses of our prediction algorithm and the losses of the 'best' expert selection strategy in the competition. We design our algorithm using the universal prediction perspective to compete against a specified class of expert selection strategies, which is not necessarily a fixed expert selection. The class of expert selection strategies that we want to compete against is purely determined by the specific application at hand and is left generic, which makes our generalized algorithm suitable for use in many different problems. We show that no preliminary knowledge about the loss sequence is required by our algorithm and its performance bounds, which are second order, expressed in terms of sums of squared losses. Our regret bounds are stable under arbitrary scalings and translations of the losses.


title: Variance Reduced EXTRA and DIGing and Their Optimal Acceleration for Strongly Convex Decentralized Optimization

abstract: We study stochastic decentralized optimization for the problem of training machine learning models with large-scale distributed data. We extend the widely used EXTRA and DIGing methods with variance reduction (VR), and propose two methods: VR-EXTRA and VR-DIGing. The proposed VR-EXTRA requires the time of $O((\kappa_s+n)\log\frac{1}{\epsilon})$ stochastic gradient evaluations and $O((\kappa_b+\kappa_c)\log\frac{1}{\epsilon})$ communication rounds to reach precision $\epsilon$, which are the best complexities among the non-accelerated gradient-type methods, where $\kappa_s$ and $\kappa_b$ are the stochastic condition number and batch condition number for strongly convex and smooth problems, respectively, $\kappa_c$ is the condition number of the communication network, and $n$ is the sample size on each distributed node. The proposed VR-DIGing has a little higher communication cost of $O((\kappa_b+\kappa_c^2)\log\frac{1}{\epsilon})$. Our stochastic gradient computation complexities are the same as the ones of single-machine VR methods, such as SAG, SAGA, and SVRG, and our communication complexities keep the same as those of EXTRA and DIGing, respectively. To further speed up the convergence, we also propose the accelerated VR-EXTRA and VR-DIGing with both the optimal $O((\sqrt{n\kappa_s}+n)\log\frac{1}{\epsilon})$ stochastic gradient computation complexity and $O(\sqrt{\kappa_b\kappa_c}\log\frac{1}{\epsilon})$ communication complexity. Our stochastic gradient computation complexity is also the same as the ones of single-machine accelerated VR methods, such as Katyusha, and our communication complexity keeps the same as those of accelerated full batch decentralized methods, such as MSDA.


title: Assessing Game Balance with AlphaZero: Exploring Alternative Rule Sets in Chess

abstract: It is non-trivial to design engaging and balanced sets of game rules. Modern chess has evolved over centuries, but without a similar recourse to history, the consequences of rule changes to game dynamics are difficult to predict. AlphaZero provides an alternative in silico means of game balance assessment. It is a system that can learn near-optimal strategies for any rule set from scratch, without any human supervision, by continually learning from its own experience. In this study we use AlphaZero to creatively explore and design new chess variants. There is growing interest in chess variants like Fischer Random Chess, because of classical chess's voluminous opening theory, the high percentage of draws in professional play, and the non-negligible number of games that end while both players are still in their home preparation. We compare nine other variants that involve atomic changes to the rules of chess. The changes allow for novel strategic and tactical patterns to emerge, while keeping the games close to the original. By learning near-optimal strategies for each variant with AlphaZero, we determine what games between strong human players might look like if these variants were adopted. Qualitatively, several variants are very dynamic. An analytic comparison show that pieces are valued differently between variants, and that some variants are more decisive than classical chess. Our findings demonstrate the rich possibilities that lie beyond the rules of modern chess.


title: Distribution networks achieve uniform perfusion through geometric self-organization

abstract: A generic flow distribution network typically does not deliver its load at a uniform rate across a service area, instead oversupplying regions near the nutrient source while leaving downstream regions undersupplied. In this work we demonstrate how a local adaptive rule coupling tissue growth with nutrient density results in a flow network that self-organizes to deliver nutrients uniformly. This geometric adaptive rule can be generalized and imported to mechanics-based adaptive models to address the effects spatial gradients in nutrients or growth factors in tissues.


title: Micromagnetic instabilities in spin-transfer switching of perpendicular magnetic tunnel junctions

abstract: Micromagnetic instabilities and non-uniform magnetization states play a significant role in spin transfer induced switching of nanometer scale magnetic elements. Here we model domain wall mediated switching dynamics in perpendicularly magnetized magnetic tunnel junction nanopillars. We show that domain wall surface tension always leads to magnetization oscillations and instabilities associated with the disk shape of the junction. A collective coordinate model is developed that captures aspects of these instabilities and illustrates their physical origin. Model results are compared to those of micromagnetic simulations. The switching dynamics is found to be very sensitive to the domain wall position and phase, which characterizes the angle of the magnetization in the disk plane. This sensitivity is reduced in the presence of spin torques and the spin current needed to displace a domain wall can be far less than the threshold current for switching from a uniformly magnetized state. A prediction of this model is conductance oscillations of increasing frequency during the switching process.


title: Filter pairs and natural extensions of logics

abstract: We adjust the notion of finitary filter pair, which was coined for creating and analyzing finitary logics, in such a way that we can treat logics of cardinality $\kappa$, where $\kappa$ is a regular cardinal. The corresponding new notion is called $\kappa$-filter pair. A filter pair can be seen as a presentation of a logic, and we ask what different $\kappa$-filter pairs give rise to a fixed logic of cardinality $\kappa$. To make the question well-defined we restrict to a subcollection of filter pairs and establish a bijection from that collection to the set of natural extensions of that logic by a set of variables of cardinality $\kappa$. Along the way we use $\kappa$-filter pairs to construct natural extensions for a given logic, work out the relationships between this construction and several others proposed in the literature, and show that the collection of natural extensions forms a complete lattice. In an optional section we introduce and motivate the concept of a general filter pair.


title: A Suzaku sample of unabsorbed narrow-line and broad-line Seyfert 1 galaxies: I. X-ray spectral properties

abstract: A sample of narrow-line (NLS1) and broad-line Seyfert 1 (BLS1) galaxies observed with Suzaku is presented. The final sample consists of 22 NLS1s and 47 BLS1s, for a total of 69 AGN that are all at low redshift (z<0.5) and exhibit low host galaxy column densities (<10^22 cm^-2). The average spectrum for each object is fit with a toy model to characterise important parameters, including the photon index, soft excess, Compton hump (or hard excess), narrow iron line strength, luminosity and X-ray Eddington ratio (L_x/L_Edd). We confirm previous findings that NLS1s have steeper power laws and higher X-ray Eddington ratios, but also find that NLS1 galaxies have stronger soft and hard excesses than their BLS1 counterparts. Studying the correlations between parameters shows that the soft and hard excesses are correlated for NLS1 galaxies, while no such correlation is observed for BLS1s. Performing a principal component analysis (PCA) on the measured X-ray parameters shows that while the X-ray Eddington ratio is the main source of variations within our sample (PC1), variations in the soft and hard excesses form the second principal component (PC2) and it is dominated by the NLS1s. The correlation between the soft and hard excess in NLS1 galaxies may suggest a common origin for the two components, such as a blurred reflection model. The presented Suzaku sample of Seyfert 1 galaxies is a useful tool for analysis of the X-ray properties of AGN, and for the study of the soft and hard excesses observed in AGN.


title: Some new series for $1/\pi$ motivated by congruences

abstract: In this paper, we deduce a family of six new series for $1/\pi$; for example, $$\sum_{n=0}^\infty\frac{41673840n+4777111}{5780^n}W_n\left(\frac{1444}{1445}\right) =\frac{147758475}{\sqrt{95}\,\pi}$$ where $W_n(x)=\sum_{k=0}^n\binom nk\binom{n+k}k\binom{2k}k\binom{2(n-k)}{n-k}x^k$. To do so, we transform our series to series of the type $$\sum_{n=0}^\infty\frac{an+b}{m^n}\sum_{k=0}^n\binom nk^4$$ studied by Cooper in 2012. In addition, we pose $17$ new series for $1/\pi$ motivated by congruences; for example, we conjecture that $$\sum_{k=0}^\infty\frac{4290k+367}{3136^k}\binom{2k}kT_k(14,1)T_k(17,16)=\frac{5390}{\pi},$$ where $T_k(b,c)$ is the coefficient of $x^k$ in the expansion of $(x^2+bx+c)^k$.


title: Tur\'an-type results for intersection graphs of boxes

abstract: In this short note, we prove the following analog of the K\H{o}v\'ari-S\'os-Tur\'an theorem for intersection graphs of boxes. If $G$ is the intersection graph of $n$ axis-parallel boxes in $\mathbb{R}^{d}$ such that $G$ contains no copy of $K_{t,t}$, then $G$ has at most $ctn(\log n)^{2d+3}$ edges, where $c=c(d)>0$ only depends on $d$. Our proof is based on exploring connections between boxicity, separation dimension and poset dimension. Using this approach, we also show that a construction of Basit et al. of $K_{2,2}$-free incidence graphs of points and rectangles in the plane can be used to disprove a conjecture of Alon et al. We show that there exist graphs of separation dimension 4 having superlinear number of edges.


title: Routing Networks with Co-training for Continual Learning

abstract: The core challenge with continual learning is catastrophic forgetting, the phenomenon that when neural networks are trained on a sequence of tasks they rapidly forget previously learned tasks. It has been observed that catastrophic forgetting is most severe when tasks are dissimilar to each other. We propose the use of sparse routing networks for continual learning. For each input, these network architectures activate a different path through a network of experts. Routing networks have been shown to learn to route similar tasks to overlapping sets of experts and dissimilar tasks to disjoint sets of experts. In the continual learning context this behaviour is desirable as it minimizes interference between dissimilar tasks while allowing positive transfer between related tasks. In practice, we find it is necessary to develop a new training method for routing networks, which we call co-training which avoids poorly initialized experts when new tasks are presented. When combined with a small episodic memory replay buffer, sparse routing networks with co-training outperform densely connected networks on the MNIST-Permutations and MNIST-Rotations benchmarks.


title: Finite-Sample Guarantees for Wasserstein Distributionally Robust Optimization: Breaking the Curse of Dimensionality

abstract: Wasserstein distributionally robust optimization (DRO) aims to find robust and generalizable solutions by hedging against data perturbations in Wasserstein distance. Despite its recent empirical success in operations research and machine learning, existing performance guarantees for generic loss functions are either overly conservative due to the curse of dimensionality, or plausible only in large sample asymptotics. In this paper, we develop a non-asymptotic framework for analyzing the out-of-sample performance for Wasserstein robust learning and the generalization bound for its related Lipschitz and gradient regularization problems. To the best of our knowledge, this gives the first finite-sample guarantee for generic Wasserstein DRO problems without suffering from the curse of dimensionality. Our results highlight that Wasserstein DRO, with a properly chosen radius, balances between the empirical mean of the loss and the variation of the loss, measured by the Lipschitz norm or the gradient norm of the loss. Our analysis is based on two novel methodological developments that are of independent interest: 1) a new concentration inequality controlling the decay rate of large deviation probabilities by the variation of the loss and, 2) a localized Rademacher complexity theory based on the variation of the loss.


title: On the Identification of Fair Auditors to Evaluate Recommender Systems based on a Novel Non-Comparative Fairness Notion

abstract: Decision-support systems are information systems that offer support to people's decisions in various applications such as judiciary, real-estate and banking sectors. Lately, these support systems have been found to be discriminatory in the context of many practical deployments. In an attempt to evaluate and mitigate these biases, algorithmic fairness literature has been nurtured using notions of comparative justice, which relies primarily on comparing two/more individuals or groups within the society that is supported by such systems. However, such a fairness notion is not very useful in the identification of fair auditors who are hired to evaluate latent biases within decision-support systems. As a solution, we introduce a paradigm shift in algorithmic fairness via proposing a new fairness notion based on the principle of non-comparative justice. Assuming that the auditor makes fairness evaluations based on some (potentially unknown) desired properties of the decision-support system, the proposed fairness notion compares the system's outcome with that of the auditor's desired outcome. We show that the proposed fairness notion also provides guarantees in terms of comparative fairness notions by proving that any system can be deemed fair from the perspective of comparative fairness (e.g. individual fairness and statistical parity) if it is non-comparatively fair with respect to an auditor who has been deemed fair with respect to the same fairness notions. We also show that the converse holds true in the context of individual fairness. A brief discussion is also presented regarding how our fairness notion can be used to identify fair and reliable auditors, and how we can use them to quantify biases in decision-support systems.


title: Driftfusion: An open source code for simulating ordered semiconductor devices with mixed ionic-electronic conducting materials in one-dimension

abstract: The recent emergence of lead-halide perovskites as active layer materials for thin film semiconductor devices including solar cells, light emitting diodes, and memristors has motivated the development of several new drift-diffusion models that include the effects of both mobile electronic and ionic charge carriers. The aim of this work is to provide a comprehensive guide to Driftfusion, a versatile simulation tool built for simulating one-dimensional ordered semiconductor devices with mixed ionic-electronic conducting layers. Driftfusion enables users to model devices with multiple, distinct, material layers and up to four charge carrier species: electrons and holes by default plus up to two ionic species. The time-dependent carrier continuity equations are fully-coupled to Poisson's equation enabling transient optoelectronic device measurement protocols to be simulated. In addition to material and device-wide properties, users have direct access to adapt the physical models for carrier transport, generation and recombination. Furthermore, a graded-interface approach circumvents the requirement for boundary conditions at material interfaces and enables interface-specific properties, such as high rates of interfacial recombination, to be introduced.


title: Moduli of stable maps with fields

abstract: We construct a moduli space of stable maps with fields associated to a triple $(X,E,s)$ of a projective variety (or a DM stack with projective moduli space) a vector bundle and a section. We show the class coincides up to a sign with the virtual fundamental class on the moduli space of stable maps to $Z=Z(s)\subset X$. We show that this gives a generalization of the Quantum Lefschetz hyperplane principle. This generalizes similar comparison results of Chang--Li, Kim--Oh and Chang--Li and presents a different approach from Chen--Janda--Webb.


title: Greybody factor and quasinormal modes of Regular Black Holes

abstract: In this work, we investigate the quasinormal frequencies of a class of regular black hole solutions which generalize Bardeen and Hayward spacetimes. In particular, we analyze scalar, vector and gravitational perturbations of the black hole both with the semianalytic WKB method and their time-domain profiles. We analyze in detail the behavior of the spectrum depending on the parameter $p/q$ of the black hole, the quantum number of angular momentum and the $s$ number. In addition, we compare our results with the classical solution valid for $p = q = 1$.


title: Quantum Thermodynamics and Quantum Coherence Engines

abstract: Advantages of quantum effects in several technologies, such as computation and communication, have already been well appreciated, and some devices, such as quantum computers and communication links, exhibiting superiority to their classical counterparts have been demonstrated. The close relationship between information and energy motivates us to explore if similar quantum benefits can be found in energy technologies. Investigation of performance limits for a broader class of information-energy machines is the subject of the rapidly emerging field of quantum thermodynamics. Extension of classical thermodynamical laws to the quantum realm is far from trivial. This short review presents some of the recent efforts in this fundamental direction and focuses on quantum heat engines and their efficiency bounds when harnessing energy from non-thermal resources, specifically those containing quantum coherence and correlations.


title: Lifespan estimates for local solutions to the semilinear wave equation in Einstein-de Sitter spacetime

abstract: In this paper, we prove some blow-up results for the semilinear wave equation in generalized Einstein-de Sitter spacetime by using an iteration argument and we derive upper bound estimates for the lifespan. In particular, we will focus on the critical cases which require the employment of a slicing procedure in the iterative mechanism. Furthermore, in order to deal with the main critical case, we will introduce a non-autonomous and parameter-dependent Cauchy problem for a linear ODE of second-order, whose explicit solution will be determined by applying the theory of special functions.


title: On good approximations and Bowen-Series expansion

abstract: We consider the continued fraction expansion of real numbers under the action of a non-uniform lattice in PSL(2,R) and prove metric relations between the convergents and a natural geometric notion of good approximations.


title: Privacy-Preserving Machine Learning in Untrusted Clouds Made Simple

abstract: We present a practical framework to deploy privacy-preserving machine learning (PPML) applications in untrusted clouds based on a trusted execution environment (TEE). Specifically, we shield unmodified PyTorch ML applications by running them in Intel SGX enclaves with encrypted model parameters and encrypted input data to protect the confidentiality and integrity of these secrets at rest and during runtime. We use the open-source Graphene library OS with transparent file encryption and SGX-based remote attestation to minimize porting effort and seamlessly provide file protection and attestation. Our approach is completely transparent to the machine learning application: the developer and the end-user do not need to modify the ML application in any way.


title: Structural, electronic, and magnetic properties of nearly-ideal $J_{\rm eff}$ $=$ 1/2 iridium halides

abstract: Heavy transition metal magnets with $J_{\rm eff}$ $=$ 1/2 electronic ground states have attracted recent interest due to their penchant for hosting new classes of quantum spin liquids and superconductors. Unfortunately, model systems with ideal $J_{\rm eff}$ $=$ 1/2 states are scarce due to the importance of non-cubic local distortions in most candidate materials. In this work, we identify a family of iridium halide systems [i.e. K$_2$IrCl$_6$, K$_2$IrBr$_6$, (NH$_4$)$_2$IrCl$_6$, and Na$_2$IrCl$_6 \cdotp $ 6(H$_2$O)] with Ir$^{4+}$ electronic ground states in extremely close proximity to the ideal $J_{\rm eff}$ $=$ 1/2 limit, despite a variation in the low-temperature global crystal structures. We also find ordered magnetic ground states for the three anhydrous systems, with single crystal neutron diffraction on K$_2$IrBr$_6$ revealing Type-I antiferromagnetism. This spin configuration is consistent with expectations for significant Kitaev exchange in a face-centered-cubic magnet.


title: Learning Anatomical Segmentations for Tractography from Diffusion MRI

abstract: Deep learning approaches for diffusion MRI have so far focused primarily on voxel-based segmentation of lesions or white-matter fiber tracts. A drawback of representing tracts as volumetric labels, rather than sets of streamlines, is that it precludes point-wise analyses of microstructural or geometric features along a tract. Traditional tractography pipelines, which do allow such analyses, can benefit from detailed whole-brain segmentations to guide tract reconstruction. Here, we introduce fast, deep learning-based segmentation of 170 anatomical regions directly on diffusion-weighted MR images, removing the dependency of conventional segmentation methods on T 1-weighted images and slow pre-processing pipelines. Working natively in diffusion space avoids non-linear distortions and registration errors across modalities, as well as interpolation artifacts. We demonstrate consistent segmentation results between 0 .70 and 0 .87 Dice depending on the tissue type. We investigate various combinations of diffusion-derived inputs and show generalization across different numbers of gradient directions. Finally, integrating our approach to provide anatomical priors for tractography pipelines, such as TRACULA, removes hours of pre-processing time and permits processing even in the absence of high-quality T 1-weighted scans, without degrading the quality of the resulting tract estimates.


title: On the strong continuity of generalised semigroups

abstract: It is well known that weakly continuous semigroups defined over $\mathbb{R}_{+}$ are automatically strongly continuous. We extend this result to more generally defined semigroups, including multiparameter semigroups.


title: Sharp isoperimetric inequalities for infinite plane graphs with bounded vertex and face degrees

abstract: We give sharp bounds for isoperimetric constants of infinite plane graphs(tessellations) with bounded vertex and face degrees. For example if $G$ is a plane graph satisfying the inequalities $p_1 \leq \mbox{deg}\, v \leq p_2$ for $v \in V(G)$ and $q_1 \leq \mbox{deg}\, f \leq q_2$ for $f \in F(G)$, where $p_1, p_2, q_1$, and $q_2$ are natural numbers such that $1/p_i + 1/q_i \leq 1/2$, $i=1,2$, then we show that \[ \Phi (p_1, q_1) \leq \inf_S \frac{|\partial S|}{|V(S)|} \leq \Phi (p_2, q_2), \] where the infimum is taken over all finite nonempty subgraphs $S \subset G$, $\partial S$ is the set of edges connecting $S$ to $G \setminus S$, and $\Phi(p,q)$ is defined by \[ \Phi (p, q) = (p-2) \sqrt{1 - \frac{4}{(p-2)(q-2)}}. \] For $p_1=3$ this gives an affirmative answer for a conjecture by Lawrencenko, Plummer, and Zha from 2002, and for general $p_i$ and $q_i$ our result fully resolves a question in the book by Lyons and Peres from 2016, where they extended the conjecture of Lawrencenko et al. to the above form. We also prove a discrete analogue of Weil's isoperimetric theorem, extending a result of Angel, Benjamini, and Horesh from 2018, and give a positive answer for a problem asked by Angel et al. in the same paper.


title: Automated Model Selection for Time-Series Anomaly Detection

abstract: Time-series anomaly detection is a popular topic in both academia and industrial fields. Many companies need to monitor thousands of temporal signals for their applications and services and require instant feedback and alerts for potential incidents in time. The task is challenging because of the complex characteristics of time-series, which are messy, stochastic, and often without proper labels. This prohibits training supervised models because of lack of labels and a single model hardly fits different time series. In this paper, we propose a solution to address these issues. We present an automated model selection framework to automatically find the most suitable detection model with proper parameters for the incoming data. The model selection layer is extensible as it can be updated without too much effort when a new detector is available to the service. Finally, we incorporate a customized tuning algorithm to flexibly filter anomalies to meet customers' criteria. Experiments on real-world datasets show the effectiveness of our solution.


title: Assessment of unsteady flow predictions using hybrid deep learning based reduced order models

abstract: In this paper, we present two deep learning-based hybrid data-driven reduced order models for the prediction of unsteady fluid flows. The first model projects the high-fidelity time series data from a finite element Navier-Stokes solver to a low-dimensional subspace via proper orthogonal decomposition (POD). The time-dependent coefficients in the POD subspace are propagated by the recurrent net (closed-loop encoder-decoder updates) and mapped to a high-dimensional state via the mean flow field and POD basis vectors. This model is referred as POD-RNN. The second model, referred to as convolution recurrent autoencoder network (CRAN), employs convolutional neural networks (CNN) as layers of linear kernels with nonlinear activations, to extract low-dimensional features from flow field snapshots. The flattened features are advanced using a recurrent (closed-loop manner) net and up-sampled (transpose convoluted) gradually to high-dimensional snapshots. Two benchmark problems of the flow past a cylinder and flow past a side-by-side cylinder are selected as the test problems to assess the efficacy of these models. For the problem of flow past a single cylinder, the performance of both the models is satisfactory, with CRAN being a bit overkill. However, it completely outperforms the POD-RNN model for a more complicated problem of flow past side-by-side cylinders. Owing to the scalability of CRAN, we briefly introduce an observer-corrector method for the calculation of integrated pressure force coefficients on the fluid-solid boundary on a reference grid. This reference grid, typically a structured and uniform grid, is used to interpolate scattered high-dimensional field data as snapshot images. These input images are convenient in training CRAN. This motivates us to further explore the application of CRAN models for the prediction of fluid flows.


title: String Backgrounds of the Yang-Baxter Deformed $AdS_4\times\mathbb{CP}^3$ Superstring

abstract: We build string backgrounds for Yang-Baxter deformations of the $AdS_4\times\mathbb{CP}^3$ superstring generated by $r$-matrices satisfying the classical Yang-Baxter equation. We obtain the metric and the NS-NS two-form of the gravity dual corresponding to noncommutative and dipole deformations of ABJM theory, as well as a deformed background with Schr\"odinger symmetry. The first two backgrounds may also be found by TsT transformations while for the last background we get a new family of non-relativistic ABJM theories with Schr\"odinger symmetry.


title: Data Augmentation for Electrocardiogram Classification with Deep Neural Network

abstract: Electrocardiogram (ECG) is the most crucial monitoring modality to diagnose cardiovascular events. Precise and automatic detection of abnormal ECG patterns is beneficial to both physicians and patients. In the automatic detection of abnormal ECG patterns, deep neural networks (DNNs) have shown significant achievements. However, DNNs require large amount of labeled data, which are often expensive to obtain. On the other hand, recent research have shown by randomly combining data augmentations can improve image classification accuracy. Thus, in this work we explore data augmentation suitable for ECG data and propose ECG Augment. We show by introducing ECG Augment, we can improve classification of atrial fibrillation with single lead ECG data, without changing an architecture of DNN.


title: Performance Analysis of FEM Solvers on Practical Electromagnetic Problems

abstract: The paper presents a comparative analysis of different commercial and academic software. The comparison aims to examine how the integrated adaptive grid refinement methodologies can deal with challenging, electromagnetic-field related problems. For this comparison, two benchmark problems were examined in the paper. The first example is a solution of an L-shape domain like test problem, which has a singularity at a certain point in the geometry. The second problem is an induction heated aluminum rod, which accurate solution needs to solve a non-linear, coupled physical fields. The accurate solution of this problem requires applying adaptive mesh generation strategies or applying a very fine mesh in the electromagnetic domain, which can significantly increase the computational complexity. The results show that the fully-hp adaptive meshing strategies, which are integrated into Agros-suite, can significantly reduce the task's computational complexity compared to the automatic h-adaptivity, which is part of the examined, popular commercial solvers.


title: A Conservative High-Order Method Utilizing Dynamic Transfinite Mortar Elements for Flow Simulation on Curved Sliding Meshes

abstract: We present a high-order method for flow simulation on unstructured curved nonconforming sliding meshes. This method utilizes dynamic transfinite mortar elements to exchange flow information between the two sides of a sliding interface. The method is arbitrarily high-order accurate in space, provably conservative, and satisfies outflow condition. Moreover, it retains the accuracy of a time marching scheme, and thus allows substantial reduction of rotational speed effects when equipped with a high-order temporal scheme. The method's capability of simultaneously handling multiple rotational objects is also explored. Details on the implementation are provided as well.


title: Is Third-Party Provided Travel Time Helpful to Estimate Freeway Performance Measures?

abstract: Transportation agencies monitor freeway performance using various measures such as VMT (Vehicle Miles Traveled), VHD (Vehicle Hours of Delay), and VHT (Vehicle Hours Traveled). Public transportation agencies typically rely on point detector data to estimate these freeway performance measures. Point detectors, such as inductive loops cannot capture the travel time for a corridor, which can lead to inaccurate performance measure estimation. This research develops a hybrid method, which estimates of freeway performance measures using a mix of probe data from third-parties and data from traditional point detectors. Using a simulated I-210 model, the overall framework using multiple data sources is evaluated, and compared with the traditional point detector-based estimation method. In the traditional method, point speeds are estimated with the flow and occupancy values using the g-factors. Data from 5% of the total vehicle are used to generate the third-party vendor provided travel time data. The analysis is conducted for multiple scenarios, including peak and off-peak periods. Findings suggest that fusing data from both third-party vendors and point detectors can help estimate performance measures better, compared to the traditional method, in scenarios that have noticeable traffic demand on freeways.


title: A Lightweight CNN Model for Detecting Respiratory Diseases from Lung Auscultation Sounds using EMD-CWT-based Hybrid Scalogram

abstract: Listening to lung sounds through auscultation is vital in examining the respiratory system for abnormalities. Automated analysis of lung auscultation sounds can be beneficial to the health systems in low-resource settings where there is a lack of skilled physicians. In this work, we propose a lightweight convolutional neural network (CNN) architecture to classify respiratory diseases using hybrid scalogram-based features of lung sounds. The hybrid scalogram features utilize the empirical mode decomposition (EMD) and continuous wavelet transform (CWT). The proposed scheme's performance is studied using a patient independent train-validation set from the publicly available ICBHI 2017 lung sound dataset. Employing the proposed framework, weighted accuracy scores of 99.20% for ternary chronic classification and 99.05% for six-class pathological classification are achieved, which outperform well-known and much larger VGG16 in terms of accuracy by 0.52% and 1.77% respectively. The proposed CNN model also outperforms other contemporary lightweight models while being computationally comparable.


title: Learning Slab Classes to Alleviate Memory Holes in Memcached

abstract: We consider the problem of memory holes in slab allocators, where an item entered into memory occupies more memory than it actually requires due to a difference between the nearest larger slab class size and the size of the entered item. We solve this problem by using a greedy algorithm that analyses the pattern of the sizes of items previously entered into the memory and accordingly re-configuring the default slab classes to better suit the learned traffic pattern to minimize memory holes. Using this approach for a consistent data pattern, in our findings, has yielded significant reductions in memory wastage. We consider Memcached as it is one of the most widely used implementations of slab allocators today, and has native support to reconfigure its default slab classes.


title: Walk Extraction Strategies for Node Embeddings with RDF2Vec in Knowledge Graphs

abstract: As KGs are symbolic constructs, specialized techniques have to be applied in order to make them compatible with data mining techniques. RDF2Vec is an unsupervised technique that can create task-agnostic numerical representations of the nodes in a KG by extending successful language modelling techniques. The original work proposed the Weisfeiler-Lehman (WL) kernel to improve the quality of the representations. However, in this work, we show both formally and empirically that the WL kernel does little to improve walk embeddings in the context of a single KG. As an alternative to the WL kernel, we propose five different strategies to extract information complementary to basic random walks. We compare these walks on several benchmark datasets to show that the \emph{n-gram} strategy performs best on average on node classification tasks and that tuning the walk strategy can result in improved predictive performances.


title: Interval hulls of $N$-matrices and almost $P$-matrices

abstract: We establish a characterization of almost $P$-matrices via a sign non-reversal property. In this we are inspired by the analogous results for $N$-matrices. Next, the interval hull of two $m \times n$ matrices $A=(a_{ij})$ and $B = (b_{ij})$, denoted by $\mathbb{I}(A,B)$, is the collection of all matrices $C \in \mathbb{R}^{m \times n}$ such that each $c_{ij}$ is a convex combination of $a_{ij}$ and $b_{ij}$. Using the sign non-reversal property, we identify a finite subset of $\mathbb{I}(A,B)$ that determines if all matrices in $\mathbb{I}(A,B)$ are $N$-matrices/almost $P$-matrices. This provides a test for an entire class of matrices simultaneously to be $N$-matrices/almost $P$-matrices. We also establish analogous results for semipositive and minimally semipositive matrices. These characterizations may be considered similar in spirit to that of $P$-matrices by Bialas-Garloff [Linear Algebra Appl. 1984] and Rohn-Rex [SIMAX 1996], and of positive definite matrices by Rohn [SIMAX 1994].


title: Palatini Higgs and Coleman-Weinberg inflation with non-minimal coupling

abstract: We present the impact of non-minimal coupling $\xi \phi^2 R$ on the inflationary parameters by taking into account the models of single-field inflation with the inflaton that has a non-zero vacuum expectation value ($v$) after the period of inflation in Palatini gravity. We discuss the well-known symmetry breaking type potentials, namely the Higgs potential and Coleman-Weinberg potential. We show the inflationary predictions of these potentials, for both $\phi>v$ and $\phi<v$ inflation, the regions in the $v-\xi$ plane for which the values of $n_s$ and $r$ are in agreement with the recent measurements. We also show the linear inflation behavior as a solution of Coleman-Weinberg potential for $\xi v^2=1$ limit. Finally, we take into account the inflationary predictions of Coleman-Weinberg potential for preferred $\xi$ values as a function of $v$ in Palatini formalism.


title: Photon-jet correlations in p-p and Pb-Pb collisions using JETSCAPE framework

abstract: It is now well established that jet modification is a multistage effect; hence a single model alone cannot describe all facets of jet modification. The JETSCAPE framework is a multistage framework that uses several modules to simulate different stages of jet propagation through the QGP medium. These simulations require a set of parameters to ensure a smooth transition between stages. We fine tune these parameters to successfully describe a variety of observables, such as the nuclear modification factors of leading hadrons and jets, jet shape, and jet fragmentation function. Photons can be produced in the hard scattering or as radiation from quarks inside jets. In this work, we study photon-jet transverse momentum imbalance and azimuthal correlation for both $p-p$ and $Pb-Pb$ collision systems. All the photons produced in each event, including the photons from hard scattering, radiation from the parton shower, and radiation from hadronization are considered with an isolation cut to directly compare with experimental data. The simulations are conducted using the same set of tuned parameters as used for the jet analysis. No new parameters are introduced or tuned. We demonstrate a significantly improved agreement with photons from $Pb-Pb$ collisions compared to prior efforts. This work provides an independent, parameter free verification of the multistage evolution framework.


title: Fairness principles for insurance contracts in the presence of default risk

abstract: We use the theory of cooperative games for the design of fair insurance contracts. An insurance contract needs to specify the premium to be paid and a possible participation in the benefit (or surplus) of the company. It results from the analysis that when a contract is exposed to the default risk of the insurance company, ex-ante equilibrium considerations require a certain participation in the benefit of the company to be specified in the contracts. The fair benefit participation of agents appears as an outcome of a game involving the residual risks induced by the default possibility and using fuzzy coalitions.


title: Study of Short-Term Personalized Glucose Predictive Models on Type-1 Diabetic Children

abstract: Research in diabetes, especially when it comes to building data-driven models to forecast future glucose values, is hindered by the sensitive nature of the data. Because researchers do not share the same data between studies, progress is hard to assess. This paper aims at comparing the most promising algorithms in the field, namely Feedforward Neural Networks (FFNN), Long Short-Term Memory (LSTM) Recurrent Neural Networks, Extreme Learning Machines (ELM), Support Vector Regression (SVR) and Gaussian Processes (GP). They are personalized and trained on a population of 10 virtual children from the Type 1 Diabetes Metabolic Simulator software to predict future glucose values at a prediction horizon of 30 minutes. The performances of the models are evaluated using the Root Mean Squared Error (RMSE) and the Continuous Glucose-Error Grid Analysis (CG-EGA). While most of the models end up having low RMSE, the GP model with a Dot-Product kernel (GP-DP), a novel usage in the context of glucose prediction, has the lowest. Despite having good RMSE values, we show that the models do not necessarily exhibit a good clinical acceptability, measured by the CG-EGA. Only the LSTM, SVR and GP-DP models have overall acceptable results, each of them performing best in one of the glycemia regions.


title: Model Fusion to Enhance the Clinical Acceptability of Long-Term Glucose Predictions

abstract: This paper presents the Derivatives Combination Predictor (DCP), a novel model fusion algorithm for making long-term glucose predictions for diabetic people. First, using the history of glucose predictions made by several models, the future glucose variation at a given horizon is predicted. Then, by accumulating the past predicted variations starting from a known glucose value, the fused glucose prediction is computed. A new loss function is introduced to make the DCP model learn to react faster to changes in glucose variations. The algorithm has been tested on 10 \textit{in-silico} type-1 diabetic children from the T1DMS software. Three initial predictors have been used: a Gaussian process regressor, a feed-forward neural network and an extreme learning machine model. The DCP and two other fusion algorithms have been evaluated at a prediction horizon of 120 minutes with the root-mean-squared error of the prediction, the root-mean-squared error of the predicted variation, and the continuous glucose-error grid analysis. By making a successful trade-off between prediction accuracy and predicted-variation accuracy, the DCP, alongside with its specifically designed loss function, improves the clinical acceptability of the predictions, and therefore the safety of the model for diabetic people.


title: Design and Fabrication of Novel Digital Transcranial Electrical Stimulator for Medical and Psychiatry Applications

abstract: In this article, we design a novel Transcranial Electrical Stimulator for medical applications, which is very cheap and can produce the desired signals very accurately. Our fabricated stimulator generates all current signals related to Transcranial Electrical Stimulation (TES) methods, i.e. Transcranial Direct Current Stimulation (tDCS), Transcranial Pulsed Current Stimulation (tPCS), Cranial Electrotherapy Stimulation (CES), and Microcurrent Electrical Therapy (MET). The proposed device has been constructed of an advanced digital controller which makes it tunable. One of the major advantages of the device is its ability to generate Burst pulses.


title: Cox-nnet v2.0: improved neural-network based survival prediction extended to large-scale EMR dataset

abstract: Cox-nnet is a neural-network based prognosis prediction method, originally applied to genomics data. Here we propose the version 2 of Cox-nnet, with significant improvement on efficiency and interpretability, making it suitable to predict prognosis based on large-scale electronic medical records (EMR) datasets. We also add permutation-based feature importance scores and the direction of feature coefficients. Applying on an EMR dataset of OPTN kidney transplantation, Cox-nnet v2.0 reduces the training time of Cox-nnet up to 32 folds (n=10,000) and achieves better prediction accuracy than Cox-PH (p<0.05). Availability and implementation: Cox-nnet v2.0 is freely available to the public at https://github.com/lanagarmire/Cox-nnet-v2.0


title: On SkipGram Word Embedding Models with Negative Sampling: Unified Framework and Impact of Noise Distributions

abstract: SkipGram word embedding models with negative sampling, or SGN in short, is an elegant family of word embedding models. In this paper, we formulate a framework for word embedding, referred to as Word-Context Classification (WCC), that generalizes SGN to a wide family of models. The framework, utilizing some "noise examples", is justified through a theoretical analysis. The impact of noise distribution on the learning of the WCC embedding models is studied experimentally, suggesting that the best noise distribution is in fact the data distribution, in terms of both the embedding performance and the speed of convergence during training. Along our way, we discover several novel embedding models that outperform the existing WCC models.


title: Solving the k-sparse Eigenvalue Problem with Reinforcement Learning

abstract: We examine the possibility of using a reinforcement learning (RL) algorithm to solve large-scale eigenvalue problems in which the desired the eigenvector can be approximated by a sparse vector with at most $k$ nonzero elements, where $k$ is relatively small compare to the dimension of the matrix to be partially diagonalized. This type of problem arises in applications in which the desired eigenvector exhibits localization properties and in large-scale eigenvalue computations in which the amount of computational resource is limited. When the positions of these nonzero elements can be determined, we can obtain the $k$-sparse approximation to the original problem by computing eigenvalues of a $k\times k$ submatrix extracted from $k$ rows and columns of the original matrix. We review a previously developed greedy algorithm for incrementally probing the positions of the nonzero elements in a $k$-sparse approximate eigenvector and show that the greedy algorithm can be improved by using an RL method to refine the selection of $k$ rows and columns of the original matrix. We describe how to represent states, actions, rewards and policies in an RL algorithm designed to solve the $k$-sparse eigenvalue problem and demonstrate the effectiveness of the RL algorithm on two examples originating from quantum many-body physics.


title: Differential Brauer Monoids

abstract: The differential Brauer monoid of a differential commutative ring R s defined. Its elements are the isomorphism classes of differential Azumaya R algebras with operation from tensor product subject to the relation that two such algebras are equivalent if matrix algebras over them are differentially isomorphic. The Bauer monoid, which is a group, is the same thing without the differential requirement. The differential Brauer monoid is then determined from the Brauer monoids of R and its ring of constants and the submonoid whose underlying Azumaya algebras are matrix rings.


title: Phasic Policy Gradient

abstract: We introduce Phasic Policy Gradient (PPG), a reinforcement learning framework which modifies traditional on-policy actor-critic methods by separating policy and value function training into distinct phases. In prior methods, one must choose between using a shared network or separate networks to represent the policy and value function. Using separate networks avoids interference between objectives, while using a shared network allows useful features to be shared. PPG is able to achieve the best of both worlds by splitting optimization into two phases, one that advances training and one that distills features. PPG also enables the value function to be more aggressively optimized with a higher level of sample reuse. Compared to PPO, we find that PPG significantly improves sample efficiency on the challenging Procgen Benchmark.


title: Mitiq: A software package for error mitigation on noisy quantum computers

abstract: We introduce Mitiq, a Python package for error mitigation on noisy quantum computers. Error mitigation techniques can reduce the impact of noise on near-term quantum computers with minimal overhead in quantum resources by relying on a mixture of quantum sampling and classical post-processing techniques. Mitiq is an extensible toolkit of different error mitigation methods, including zero-noise extrapolation, probabilistic error cancellation, and Clifford data regression. The library is designed to be compatible with generic backends and interfaces with different quantum software frameworks. We describe Mitiq using code snippets to demonstrate usage and discuss features and contribution guidelines. We present several examples demonstrating error mitigation on IBM and Rigetti superconducting quantum processors as well as on noisy simulators.


title: Measuring the primordial gravitational-wave background in the presence of astrophysical foregrounds

abstract: Primordial gravitational waves are expected to create a stochastic background encoding information about the early Universe that may not be accessible by other means. However, the primordial background is obscured by an astrophysical foreground consisting of gravitational waves from compact binaries. We demonstrate a Bayesian method for estimating the primordial background in the presence of an astrophysical foreground. Since the background and foreground signal parameters are estimated simultaneously, there is no subtraction step, and therefore we avoid astrophysical contamination of the primordial measurement, sometimes referred to as "residuals". Additionally, since we include the non-Gaussianity of the astrophysical foreground in our model, this method represents the statistically optimal approach to the simultaneous detection of a multi-component stochastic background.


title: Pathways to Sustainable Planetary Science

abstract: Climate change is a major impending threat to the future of humanity. According to the International Panel on Climate Change (IPCC), our emissions are estimated to have caused 0.8 deg C-1.2 deg C of anthropogenic global warming (AGW) above pre-industrial levels. AGW is likely to reach 1.5 degrees C between 2030 and 2052 if it continues to increase at the current rate. As the climate change is driven by the release of carbon dioxide and other greenhouse gases (GHG) into the atmosphere, there is a broad consensus that the mitigation of climate change requires transition to low GHG emission energy sources, technologies and practices. Implementing such changes systematically from individual to community-wide scales together with the resulting cultural changes and leadership towards environmental consciousness and responsibility are crucial to mitigate the looming damage of AGW. Given planetary scientists' wide recognition of the realities of climate change, and the need for us to maintain credibility by leading by example, it is appropriate to make own professional behavior more environmentally responsible. While scientists are few in numbers, and planetary scientists far fewer, high volumes of academic travel to conferences, panels, colloquia, and research collaboration visits together with extensive use of large, energetically demanding infrastructures make the "carbon footprint" of scientists much higher than that of an average citizen. This White Paper focuses on how modifying our activities, particularly associated with academic travel, can affect the carbon footprint of the planetary science community, and it makes recommendations on how the community and the funding agencies could best participate in the cultural change required to mitigate the damage that AGW will cause.


title: Cephalogram Synthesis and Landmark Detection in Dental Cone-Beam CT Systems

abstract: Due to the lack of standardized 3D cephalometric analytic methodology, 2D cephalograms synthesized from 3D cone-beam computed tomography (CBCT) volumes are widely used for cephalometric analysis in dental CBCT systems. However, compared with conventional X-ray film based cephalograms, such synthetic cephalograms lack image contrast and resolution. In addition, the radiation dose during the scan for 3D reconstruction causes potential health risks. In this work, we propose a sigmoid-based intensity transform that uses the nonlinear optical property of X-ray films to increase image contrast of synthetic cephalograms. To improve image resolution, super resolution deep learning techniques are investigated. For low dose purpose, the pixel-to-pixel generative adversarial network (pix2pixGAN) is proposed for 2D cephalogram synthesis directly from two CBCT projections. For landmark detection in the synthetic cephalograms, an efficient automatic landmark detection method using the combination of LeNet-5 and ResNet50 is proposed. Our experiments demonstrate the efficacy of pix2pixGAN in 2D cephalogram synthesis, achieving an average peak signal-to-noise ratio (PSNR) value of 33.8 with reference to the cephalograms synthesized from 3D CBCT volumes. Pix2pixGAN also achieves the best performance in super resolution, achieving an average PSNR value of 32.5 without the introduction of checkerboard or jagging artifacts. Our proposed automatic landmark detection method achieves 86.7% successful detection rate in the 2 mm clinical acceptable range on the ISBI Test1 data, which is comparable to the state-of-the-art methods. The method trained on conventional cephalograms can be directly applied to landmark detection in the synthetic cephalograms, achieving 93.0% and 80.7% successful detection rate in 4 mm precision range for synthetic cephalograms from 3D volumes and 2D projections respectively.


title: Sympathetic cooling of a radio-frequency LC circuit to its ground state in an optoelectromechanical system

abstract: We present a complete theory for laser cooling of a macroscopic radio-frequency LC electrical circuit by means of an optoelectromechanical system, consisting of an optical cavity dispersively coupled to a nanomechanical oscillator, which is in turn capacitively coupled to the LC circuit of interest. We determine the optimal parameter regime where the LC resonator can be cooled down to its quantum ground state, which requires a large optomechanical cooperativity, and a larger electromechanical cooperativity. Moreover, comparable optomechanical and electromechanical coupling rates are preferable for reaching the quantum ground state.


title: Learning-Based Distributionally Robust Model Predictive Control of Markovian Switching Systems with Guaranteed Stability and Recursive Feasibility

abstract: We present a data-driven model predictive control scheme for chance-constrained Markovian switching systems with unknown switching probabilities. Using samples of the underlying Markov chain, ambiguity sets of transition probabilities are estimated which include the true conditional probability distributions with high probability. These sets are updated online and used to formulate a time-varying, risk-averse optimal control problem. We prove recursive feasibility of the resulting MPC scheme and show that the original chance constraints remain satisfied at every time step. Furthermore, we show that under sufficient decrease of the confidence levels, the resulting MPC scheme renders the closed-loop system mean-square stable with respect to the true-but-unknown distributions, while remaining less conservative than a fully robust approach.


title: Superconductivity-enhanced spin pumping: Role of Andreev resonances

abstract: We describe a simple hybrid superconductor$|$ferromagnetic-insulator structure manifesting spin-resolved Andreev bound states in which dynamic magnetization is employed to probe spin related physics. We show that, at low bias and below $T_c$, the transfer of spin angular momentum pumped by an externally driven ferromagnetic insulator is greatly affected by the formation of spin-resolved Andreev bound states. Our results indicate that these bound states capture the essential physics of condensate-facilitated spin flow. For finite thicknesses of the superconducting layer, comparable to the coherence length, resonant Andreev bound states render highly transmitting subgap spin transport channels. We point out that the resonant enhancement of the subgap transport channels establishes a prototype Fabry-P\'erot resonator for spin pumping.


title: A novel highly efficient Lagrangian model for massively multidomain simulations: parallel context

abstract: A new method for the simulation of evolving multi-domains problems has been introduced in a previous work (RealIMotion), Florez et al. (2020). In this article further developments of the model will be presented. The main focus here is a robust parallel implementation using a distributed-memory approach with the Message Passing Interface (MPI) library OpenMPI. The original 2D sequential methodology consists in a modified front-tracking approach where the main originality is that not only interfaces between domains are discretized but their interiors are also meshed. The interfaces are tracked based on the topological degree of each node on the mesh and the remeshing and topological changes of the domains are driven by selective local operations performed over an element patch. The accuracy and the performance of the sequential method has proven very promising in Florez et al. (2020). In this article a parallel implementation will be discussed and tested in context of motion by curvature flow for polycrystals, i.e. by considering Grain Growth (GG) mechanism. Results of the performance of the model are given and comparisons with other approaches in the literature are discussed.


title: (In)Existence of Equilibria for 2-Players, 2-Values Games with Concave Valuations

abstract: We consider 2-players, 2-values minimization games where the players' costs take on two values, $a,b$, $a<b$. The players play mixed strategies and their costs are evaluated by unimodal valuations. This broad class of valuations includes all concave, one-parameter functions $\mathsf{F}: [0,1]\rightarrow \mathbb{R}$ with a unique maximum point. Our main result is an impossibility result stating that: If the maximum is obtained in $(0,1)$ and $\mathsf{F}\left(\frac{1}{2}\right)\ne b$, then there exists a 2-players, 2-values game without $\mathsf{F}$-equilibrium. The counterexample game used for the impossibility result belongs to a new class of very sparse 2-players, 2-values bimatrix games which we call normal games. In an attempt to investigate the remaining case $\mathsf{F}\left(\frac{1}{2}\right) = b$, we show that: - Every normal, $n$-strategies game has an ${\mathsf{F}}$-equilibrium when ${\mathsf{F}}\left( \frac{1}{2} \right) = b$. We present a linear time algorithm for computing such an equilibrium. - For 2-players, 2-values games with 3 strategies we have that if $\mathsf{F}\left(\frac{1}{2}\right) \le b$, then every 2-players, 2-values, 3-strategies game has an $\mathsf{F}$-equilibrium; if $\mathsf{F}\left(\frac{1}{2}\right) > b$, then there exists a normal 2-players, 2-values, 3-strategies game without $\mathsf{F}$-equilibrium. To the best of our knowledge, this work is the first to provide an (almost complete) answer on whether there is, for a given concave function $\mathsf{F}$, a counterexample game without $\mathsf{F}$-equilibrium.


title: CuratorNet: Visually-aware Recommendation of Art Images

abstract: Although there are several visually-aware recommendation models in domains like fashion or even movies, the art domain lacks thesame level of research attention, despite the recent growth of the online artwork market. To reduce this gap, in this article we introduceCuratorNet, a neural network architecture for visually-aware recommendation of art images. CuratorNet is designed at the core withthe goal of maximizing generalization: the network has a fixed set of parameters that only need to be trained once, and thereafter themodel is able to generalize to new users or items never seen before, without further training. This is achieved by leveraging visualcontent: items are mapped to item vectors through visual embeddings, and users are mapped to user vectors by aggregating the visualcontent of items they have consumed. Besides the model architecture, we also introduce novel triplet sampling strategies to build atraining set for rank learning in the art domain, resulting in more effective learning than naive random sampling. With an evaluationover a real-world dataset of physical paintings, we show that CuratorNet achieves the best performance among several baselines,including the state-of-the-art model VBPR. CuratorNet is motivated and evaluated in the art domain, but its architecture and trainingscheme could be adapted to recommend images in other areas


title: Multimessenger parameter estimation of GW170817: from jet structure to the Hubble constant

abstract: The electromagnetic radiation that followed the neutron star merger event GW170817 revealed that gamma-ray burst afterglows from jets misaligned with our line of sight exhibit a light curve with slowly rising flux. The slope of the rising light curve depends sensitively on the angle of the observer with respect to the jet axis, which is likely to be perpendicular to the merger plane of the neutron star binary. Therefore, the afterglow emission can be used to constrain the inclination of the merging system. Here, we calculate the gamma-ray burst afterglow emission based on the realistic jet structure derived from general-relativistic magnetohydrodynamical simulations of a black hole torus system for the central engine of the gamma-ray burst. Combined with gravitational wave parameter estimation, we fit the multi-epoch afterglow emission of GW170817. We show that with such a jet model, the observing angle can be tightly constrained by multi messenger observations. The best fit observing angle of GW170817 is $\theta_{\rm v} = 0.38\pm 0.02$ rad. With such a constraint, we can break the degeneracy between inclination angle and luminosity distance in gravitational wave parameter estimation, and substantially increase the precision with which the Hubble constant is constrained by the standard siren method. Our estimation of the distance is $D_{\rm L}=43.4\pm 1\ \rm Mpc$ and the Hubble constant constraint is $69.5\pm 4\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$. As a result, multimessenger observations of short-duration gamma-ray bursts, combined with a good theoretical understanding of the jet structure, can be powerful probes of cosmological parameters.


title: Identification and estimation of Structural VARMA models using higher order dynamics

abstract: We use information from higher order moments to achieve identification of non-Gaussian structural vector autoregressive moving average (SVARMA) models, possibly non-fundamental or non-causal, through a frequency domain criterion based on a new representation of the higher order spectral density arrays of vector linear processes. This allows to identify the location of the roots of the determinantal lag matrix polynomials based on higher order cumulants dynamics and to identify the rotation of the model errors leading to the structural shocks up to sign and permutation. We describe sufficient conditions for global and local parameter identification that rely on simple rank assumptions on the linear dynamics and on finite order serial and component independence conditions for the structural innovations. We generalize previous univariate analysis to develop asymptotically normal and efficient estimates exploiting second and non-Gaussian higher order dynamics given a particular structural shocks ordering without assumptions on causality or invertibility. Bootstrap approximations to finite sample distributions and the properties of numerical methods are explored with real and simulated data.


title: Interpretation of structure in the di-$J/\psi$ spectrum

abstract: Structure in the di-$J/\psi$ mass spectrum observed by the LHCb experiment around 6.9 and 7.2 GeV is interpreted in terms of $J^{PC}=0^{++}$ and $2^{++}$ resonances between a $cc$ diquark and a $\bar c \bar c$ antidiquark, using a recently confirmed string-junction picture to calculate tetraquark masses. The main peak around 6.9 GeV is likely dominated by the $0^{++}(2S)$ state, a radial excitation of the $cc$-$\bar c \bar c$ tetraquark, which we predict at $6.871\pm 0.025$ GeV. The dip around 6.75 GeV is ascribed to the opening of the \Swave\ di-$\chi_{c0}$ channel, while the dip around 7.2 GeV could be correlated with the opening of the di-$\eta_c(2S)$ or $\Xi_{cc} \bar \Xi_{cc}$ channel. The low-mass part of the di-$J/\psi$ structure appears to require a broad resonance consistent with a predicted $2^{++}(1S)$ state with invariant mass around $M_{\rm inv} = 6400$ MeV. Implications for $bb \bar b \bar b$ tetraquarks are discussed.


title: Semi-discrete optimal transport methods for the semi-geostrophic equations

abstract: We give a new and constructive proof of the existence of global-in-time weak solutions of the 3-dimensional incompressible semi-geostrophic equations (SG) in geostrophic coordinates, for arbitrary initial measures with compact support. This new proof, based on semi-discrete optimal transport techniques, works by characterising discrete solutions of SG in geostrophic coordinates in terms of trajectories satisfying an ordinary differential equation. It is advantageous in its simplicity and its explicit relation to Eulerian coordinates through the use of Laguerre tessellations. Using our method, we obtain improved time-regularity for a large class of discrete initial measures, and we compute explicitly two discrete solutions. The method naturally gives rise to an efficient numerical method, which we illustrate by presenting simulations of a 2-dimensional semi-geostrophic flow in geostrophic coordinates generated using a numerical solver for the semi-discrete optimal transport problem coupled with an ordinary differential equation solver.


title: Explicit Tamagawa numbers for certain algebraic tori over number fields

abstract: Given a number field extension $K/k$ with an intermediate field $K^+$ fixed by a central element of the corresponding Galois group of prime order $p$, we build an algebraic torus over $k$ whose rational points are elements of $K^\times$ sent to $k^\times$ via the norm map $N_{K/K^+}$. The goal is to compute the Tamagawa number of that torus explicitly via Ono's formula that expresses it as a ratio of cohomological invariants. A fairly complete and detailed description of the cohomology of the character lattice of such a torus is given when $K/k$ is Galois. Partial results including the numerator are given when the extension is not Galois, or more generally when the torus is defined by an \'etale algebra. We also present tools developed in SAGE for this purpose, allowing us to build and compute the cohomology and explore the local-global principles for such an algebraic torus. Particular attention is given to the case when $[K:K^+]=2$ and $K$ is a CM-field. This case corresponds to tori in $\mathrm{GSp}_{2n}$, and most examples will be in that setting. This is motivated by the application to abelian varieties over finite fields and the Hasse principle for bilinear forms.


title: Smooth Converse Lyapunov-Barrier Theorems for Asymptotic Stability with Safety Constraints and Reach-Avoid-Stay Specifications

abstract: Stability and safety are two important aspects in safety-critical control of dynamical systems. It has been a well established fact in control theory that stability properties can be characterized by Lyapunov functions. Reachability properties can also be naturally captured by Lyapunov functions for finite-time stability. Motivated by safety-critical control applications, such as in autonomous systems and robotics, there has been a recent surge of interests in characterizing safety properties using barrier functions. Lyapunov and barrier functions conditions, however, are sometimes viewed as competing objectives. In this paper, we provide a unified theoretical treatment of Lyapunov and barrier functions in terms of converse theorems for stability properties with safety guarantees and reach-avoid-stay type specifications. We show that if a system (modeled as a dynamical system with measurable perturbations) possesses a stability with safety property, then there exists a smooth Lyapunov function to certify such a property. This Lyapunov function is shown to be defined on the entire set of initial conditions from which solutions satisfy this property. A similar but slightly weaker statement is made for reach-avoid-stay specifications. We show by a simple example that the latter statement cannot be strengthened without additional assumptions. We further extend the results for systems with control inputs and prove existence of converse Lyapunov-barrier functions for reach-and-avoid specifications. While the converse Lyapunov-barrier theorems are not constructive, as with classical converse Lyapunov theorems, we believe that the unified necessary and sufficient conditions with a single Lyapunov-barrier function are of theoretical interest and can hopefully shed some light on computational approaches.


title: not-so-BigGAN: Generating High-Fidelity Images on Small Compute with Wavelet-based Super-Resolution

abstract: State-of-the-art models for high-resolution image generation, such as BigGAN and VQVAE-2, require an incredible amount of compute resources and/or time (512 TPU-v3 cores) to train, putting them out of reach for the larger research community. On the other hand, GAN-based image super-resolution models, such as ESRGAN, can not only upscale images to high dimensions, but also are efficient to train. In this paper, we present not-so-big-GAN (nsb-GAN), a simple yet cost-effective two-step training framework for deep generative models (DGMs) of high-dimensional natural images. First, we generate images in low-frequency bands by training a sampler in the wavelet domain. Then, we super-resolve these images from the wavelet domain back to the pixel-space with our novel wavelet super-resolution decoder network. Wavelet-based down-sampling method preserves more structural information than pixel-based methods, leading to significantly better generative quality of the low-resolution sampler (e.g., 64x64). Since the sampler and decoder can be trained in parallel and operate on much lower dimensional spaces than end-to-end models, the training cost is substantially reduced. On ImageNet 512x512, our model achieves a Fr\'echet Inception Distance (FID) of 10.59 -- beating the baseline BigGAN model -- at half the compute (256 TPU-v3 cores).


title: Newtonian-like gravity with variable $G$

abstract: We propose a Lagrangian formulation for a varying $G$ Newtonian-like theory inspired by the Brans-Dicke gravity. Rather than imposing an {\it ad hoc} dependence for the gravitational coupling, as previously done in the literature, in our proposal the running of $G$ emerges naturally from the internal dynamical structure of the theory. We explore the features of the resulting gravitational field for static and spherically symmetric mass distributions as well as within the cosmological framework.


title: Vacuum decay in quadratic gravity

abstract: Metastable states decay at zero temperature through quantum tunneling at an exponentially small rate, which depends on the Coleman-de Luccia instanton, also known as bounce. In some theories, the bounce may not exist or its on-shell action may be ill-defined or infinite, thus hindering the vacuum decay process. In this paper, we test this possibility in modified theories of gravity interacting with a real scalar field. We consider an Einstein-Hilbert term with a non-minimally coupled scalar field and a quadratic Ricci scalar contribution. To tackle the problem we use a new analytic method, with which we prove that the scalar field on the bounce has a universal behaviour at large Euclidean radii, almost independently of the potential. Our main result is that the quadratic Ricci scalar prevents the decay, regardless of the other terms in the action.


title: Coagulation equations for aerosol dynamics

abstract: Binary coagulation is an important process in aerosol dynamics by which two particles merge to form a larger one. The distribution of particle sizes over time may be described by the so-called Smoluchowski's coagulation equation. This integrodifferential equation exhibits complex non-local behaviour that strongly depends on the coagulation rate considered. We first discuss well-posedness results for the Smoluchowski's equation for a large class of coagulation kernels as well as the existence and nonexistence of stationary solutions in the presence of a source of small particles. The existence result uses Schauder fixed point theorem, and the nonexistence result relies on a flux formulation of the problem and on power law estimates for the decay of stationary solutions with a constant flux. We then consider a more general setting. We consider that particles may be constituted by different chemicals, which leads to multi-component equations describing the distribution of particle compositions. We obtain explicit solutions in the simplest case where the coagulation kernel is constant by using a generating function. Using an approximation of the solution we observe that the mass localizes along a straight line in the size space for large times and large sizes.


title: Ties between Parametrically Polymorphic Type Systems and Finite Control Automata

abstract: We present a correspondence and bisimulation between variants of parametrically polymorphic type systems and variants of finite control automata, such as FSA, PDA, tree automata and Turing machine. Within this correspondence we show that two recent celebrated results on automatic generation of fluent API are optimal in certain senses, present new results on the studied type systems, formulate open problems, and present potential software engineering applications, other than fluent API generation, which may benefit from judicious use of type theory.


title: Surveying the side-chain network approach to protein structure and dynamics: The SARS-CoV-2 spike protein as an illustrative case

abstract: Network theory-based approaches provide valuable insights into the variations in global structural connectivity between differing dynamical states of proteins. Our objective is to review network-based analyses to elucidate such variations, especially in the context of subtle conformational changes. We present technical details of the construction and analyses of protein structure networks, encompassing both the non-covalent connectivity and dynamics. We examine the selection of optimal criteria for connectivity based on the physical concept of percolation. We highlight the advantages of using side-chain based network metrics in contrast to backbone measurements. As an illustrative example, we apply the described network approach to investigate the global conformational change between the closed and partially open states of the SARS-CoV-2 spike protein. This conformational change in the spike protein is crucial for coronavirus entry and fusion into human cells. Our analysis reveals global structural reorientations between the two states of the spike protein despite small changes between the two states at the backbone level. We also observe some differences at strategic locations in the structures, correlating with their functions, asserting the advantages of the side-chain network analysis. Finally we present a view of allostery as a subtle synergistic-global change between the ligand and the receptor, the incorporation of which would enhance the drug design strategies.


title: Gravitational-wave Merger Forecasting: Scenarios for the early detection and localization of compact-binary mergers with ground based observatories

abstract: We present the prospects for the early (pre-merger) detection and localization of compact-binary coalescences using gravitational waves over the next 10 years. Early warning can enable the direct observation of the prompt and early electromagnetic emission of a neutron star merger. We examine the capabilities of the ground based detectors at their "Design" sensitivity (2021-2022), the planned "A+" upgrade (2024-2026), and the envisioned "Voyager" concept (late 2020's). We find that for a fiducial rate of binary neutron star mergers of $1000 ~\mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, the Design, A+, and Voyager era networks can provide 18, 54, and 195s of warning for one source per year of observing, respectively, with a sky localization area $<$100 deg$^2$ at a $90\%$ credible level. At the same rate, the A+ and Voyager era networks will be able to provide 9 and 43s of warning, respectively, for a source with $<$10 deg$^2$ localization area. We compare the idealized search sensitivity to that achieved by the PyCBC Live search tuned for pre-merger detection. The gravitational-wave community will be prepared to produce pre-merger alerts. Our results motivate the operation of observatories with wide fields-of-view, automation, and the capability for fast slewing to observe simultaneously with the gravitational-wave network.