title: One-shot dynamical resource theory

abstract: A fundamental problem in resource theory is to study the manipulation of the resource. Focusing on a general dynamical resource theory of quantum channels, here we consider tasks of one-shot resource distillation and dilution with a single copy of the resource. For any target of unitary channel or pure state preparation channel, we establish a universal strategy to determine upper and lower bounds on rates that convert between any given resource and the target. We show that the rates are related to resource measures based on the channel robustness and the channel hypothesis testing entropy, with regularization factors of the target resource measures. The strategy becomes optimal with converged bounds when the channel robustness is finite and measures of the target resource collapse to the same value. The single-shot result also applies to asymptotic parallel manipulation of channels to obtain asymptotic resource conversion rates. We give several examples of dynamical resources, including the purity, classical capacity, quantum capacity, non-uniformity, coherence, and entanglement of quantum channels. Our results are applicable to general dynamical resource theories with potential applications in quantum communication, fault-tolerant quantum computing, and quantum thermodynamics.


title: Batch Group Normalization

abstract: Deep Convolutional Neural Networks (DCNNs) are hard and time-consuming to train. Normalization is one of the effective solutions. Among previous normalization methods, Batch Normalization (BN) performs well at medium and large batch sizes and is with good generalizability to multiple vision tasks, while its performance degrades significantly at small batch sizes. In this paper, we find that BN saturates at extreme large batch sizes, i.e., 128 images per worker, i.e., GPU, as well and propose that the degradation/saturation of BN at small/extreme large batch sizes is caused by noisy/confused statistic calculation. Hence without adding new trainable parameters, using multiple-layer or multi-iteration information, or introducing extra computation, Batch Group Normalization (BGN) is proposed to solve the noisy/confused statistic calculation of BN at small/extreme large batch sizes with introducing the channel, height and width dimension to compensate. The group technique in Group Normalization (GN) is used and a hyper-parameter G is used to control the number of feature instances used for statistic calculation, hence to offer neither noisy nor confused statistic for different batch sizes. We empirically demonstrate that BGN consistently outperforms BN, Instance Normalization (IN), Layer Normalization (LN), GN, and Positional Normalization (PN), across a wide spectrum of vision tasks, including image classification, Neural Architecture Search (NAS), adversarial learning, Few Shot Learning (FSL) and Unsupervised Domain Adaptation (UDA), indicating its good performance, robust stability to batch size and wide generalizability. For example, for training ResNet-50 on ImageNet with a batch size of 2, BN achieves Top1 accuracy of 66.512% while BGN achieves 76.096% with notable improvement.


title: Six-dimensional time-space crystalline structures

abstract: Time crystalline structures are characterized by regularity that single-particle or many-body systems manifest in the time domain, closely resembling the spatial regularity of ordinary space crystals. Here we show that time and space crystalline structures can be combined together and even six-dimensional time-space lattices can be realized. As an example, we demonstrate that such time-space crystalline structures can reveal the six-dimensional quantum Hall effect quantified by the third Chern number.

categories: cond-mat.quant-gas cond-mat.mes-hall physics.atom-ph quant-ph


title: Well-posedness and direct internal stability of coupled non-degenrate Kirchhoff system via heat conduction

abstract: In the paper under study, we consider the following coupled non-degenerate Kirchhoff system \begin{equation}\label{P} \left \{ \begin{aligned} &\displaystyle y_{tt}-\upvarphi\Big(\int_\Omega | \nabla y |^2\,dx\Big)\Delta y +\upalpha \Delta \uptheta=0, &\mbox{ in }&\; \Omega \times (0, +\infty)\\ &\displaystyle \uptheta_t-\Delta \uptheta-\upbeta \Delta y_t =0, &\mbox{ in }&\; \Omega \times (0, +\infty)\\ &\displaystyle y=\uptheta=0,\; &\mbox{ on }&\;\partial\Omega\times(0, +\infty)\\ %&\displaystyle y=0,\; &\mbox{ on }&\;\partial\Omega\times(0, +\infty)\\ %&\displaystyle \partial_\nu y=0, &\mbox{ on }&\;\Gamma_1\times(0, +\infty)\\ &\displaystyle y(\cdot, 0)=y_0, \; y_t(\cdot, 0)=y_1,\;\uptheta(\cdot, 0)=\uptheta_0, \; \; &\mbox{ in }&\; \Omega\\ \end{aligned} \right. \end{equation} where $\Omega$ is a bounded open subset of $\mathbb{R}^n$, $\upalpha$ and $\upbeta$ be two nonzero real numbers with the same sign and $\upvarphi$ is given by $\upvarphi(s)= \mathfrak{m}_0+\mathfrak{m}_1s$ with some positive constants $\mathfrak{m}_0$ and $\mathfrak{m}_1$. So we prove existence of solution and establish its exponential decay. The method used is based on multiplier technique and some integral inequalities due to Haraux and Komornik\cite{H1,KOM}.


title: Unsupervised embedding of trajectories captures the latent structure of scientific migration

abstract: Human migration and mobility drives major societal phenomena including epidemics, economies, innovation, and the diffusion of ideas. Although human mobility and migration have been heavily constrained by geographic distance throughout the history, advances and globalization are making other factors such as language and culture increasingly more important. Advances in neural embedding models, originally designed for natural language, provide an opportunity to tame this complexity and open new avenues for the study of migration. Here, we demonstrate the ability of the model word2vec to encode nuanced relationships between discrete locations from migration trajectories, producing an accurate, dense, continuous, and meaningful vector-space representation. The resulting representation provides a functional distance between locations, as well as a digital double that can be distributed, re-used, and itself interrogated to understand the many dimensions of migration. We show that the unique power of word2vec to encode migration patterns stems from its mathematical equivalence with the gravity model of mobility. Focusing on the case of scientific migration, we apply word2vec to a database of three million migration trajectories of scientists derived from the affiliations listed on their publication records. Using techniques that leverage its semantic structure, we demonstrate that embeddings can learn the rich structure that underpins scientific migration, such as cultural, linguistic, and prestige relationships at multiple levels of granularity. Our results provide a theoretical foundation and methodological framework for using neural embeddings to represent and understand migration both within and beyond science.


title: Two populations of LIGO-Virgo black holes

abstract: We analyse the LIGO-Virgo data, including the recently released GWTC-2 dataset, to test a hypothesis that the data contains more than one population of black holes. We perform a maximum likelihood analysis including a population of astrophysical black holes with a truncated power-law mass function whose merger rate follows from star formation rate, and a population of primordial black holes for which we consider log-normal and critical collapse mass functions. We find that primordial black holes alone are strongly disfavoured by the data, while the best fit is obtained for the template combining astrophysical and primordial merger rates. Alternatively, the data may hint towards two different astrophysical black hole populations. We also update the constraints on primordial black hole abundance from LIGO-Virgo observations finding that in the $2-400M_\odot$ mass range they must comprise less than 0.2% of dark matter.


title: A Time-Varying Fine Structure Constant from Naturally Ultralight Dark Matter

abstract: We present a class of models in which the coupling of the photon to an ultralight scalar field that has a time-dependent vacuum expectation value causes the fine structure constant to oscillate in time. The scalar field is assumed to constitute all or part of the observed dark matter. Its mass is protected against radiative corrections by a discrete $\mathbb{Z}_N$ exchange symmetry that relates the Standard Model to several copies to itself. The abundance of dark matter is set by the misalignment mechanism. We show that the oscillations in the fine structure constant are large enough to be observed in current and near-future experiments.


title: Neural Dynamic Policies for End-to-End Sensorimotor Learning

abstract: The current dominant paradigm in sensorimotor control, whether imitation or reinforcement learning, is to train policies directly in raw action spaces such as torque, joint angle, or end-effector position. This forces the agent to make decisions individually at each timestep in training, and hence, limits the scalability to continuous, high-dimensional, and long-horizon tasks. In contrast, research in classical robotics has, for a long time, exploited dynamical systems as a policy representation to learn robot behaviors via demonstrations. These techniques, however, lack the flexibility and generalizability provided by deep learning or reinforcement learning and have remained under-explored in such settings. In this work, we begin to close this gap and embed the structure of a dynamical system into deep neural network-based policies by reparameterizing action spaces via second-order differential equations. We propose Neural Dynamic Policies (NDPs) that make predictions in trajectory distribution space as opposed to prior policy learning methods where actions represent the raw control space. The embedded structure allows end-to-end policy learning for both reinforcement and imitation learning setups. We show that NDPs outperform the prior state-of-the-art in terms of either efficiency or performance across several robotic control tasks for both imitation and reinforcement learning setups. Project video and code are available at https://shikharbahl.github.io/neural-dynamic-policies/


title: Holding the Cosmos in Your Hand: Developing 3D Modeling and Printing Pipelines for Communications and Research

abstract: Three-dimensional (3D) visualization has opened up a Universe of possible scientific data representations. 3D printing has the potential to make seemingly abstract and esoteric data sets accessible, particularly through the lens of translating data into forms that can be explored in the tactile modality for people who are blind or visually impaired. This article will briefly review 3D modeling in astrophysics, astronomy, and planetary science, before discussing 3D printed astrophysical and planetary geophysical data sets and their current and potential applications with non-expert audiences. The article will also explore the prospective pipeline and benefits of other 3D data outputs in accessible scientific research and communications, including extended reality and data sonification.


title: Nutation in antiferromagnetic resonance

abstract: The effect of inertial spin dynamics is compared between ferromagnetic, antiferromagnetic and ferrimagnetic systems. The linear response to an oscillating external magnetic field is calculated within the framework of the inertial Landau--Lifshitz--Gilbert equation using analytical theory and computer simulations. Precession and nutation resonance peaks are identified, and it is demonstrated that the precession frequencies are reduced by the spin inertia, while the lifetime of the excitations is enhanced. The interplay between precession and nutation is found to be the most prominent in antiferromagnets, where the timescale of the exchange-driven sublattice dynamics is comparable to inertial relaxation times. Consequently, antiferromagnetic resonance techniques should be better suited for the search for intrinsical inertial spin dynamics on ultrafast timescales than ferromagnetic resonance.


title: A Unified Model for Gate Level Propagation Analysis

abstract: Classic hardware verification techniques (e.g., X-propagation and fault-propagation) and more recent hardware security verification techniques based on information flow tracking (IFT) aim to understand how information passes, affects, and otherwise modifies a circuit. These techniques all have separate usage scenarios, but when dissected into their core functionality, they relate in a fundamental manner. In this paper, we develop a common framework for gate level propagation analysis. We use our model to generate synthesizable propagation logic to use in standard EDA tools. To justify our model, we prove that Precise Hardware IFT is equivalent to gate level X-propagation and imprecise fault propagation. We also show that the difference between Precise Hardware IFT and fault propagation is not significant for 74X-series and '85 ISCAS benchmarks with more than 313 gates and the difference between imprecise hardware IFT and Precise Hardware IFT is almost always significant regardless of size.


title: Weight Update Skipping: Reducing Training Time for Artificial Neural Networks

abstract: Artificial Neural Networks (ANNs) are known as state-of-the-art techniques in Machine Learning (ML) and have achieved outstanding results in data-intensive applications, such as recognition, classification, and segmentation. These networks mostly use deep layers of convolution or fully connected layers with many filters in each layer, demanding a large amount of data and tunable hyperparameters to achieve competitive accuracy. As a result, storage, communication, and computational costs of training (in particular training time) become limiting factors to scale them up. In this paper, we propose a new training methodology for ANNs that exploits the observation of improvement of accuracy shows temporal variations which allow us to skip updating weights when the variation is minuscule. During such time windows, we keep updating bias which ensures the network still trains and avoids overfitting; however, we selectively skip updating weights (and their time-consuming computations). Such a training approach virtually achieves the same accuracy with considerably less computational cost, thus lower training time. We propose two methods for updating weights and evaluate them by analyzing four state-of-the-art models, AlexNet, VGG-11, VGG-16, ResNet-18 on CIFAR datasets. On average, our two proposed methods called WUS and WUS+LR reduced the training time (compared to the baseline) by 54%, and 50%, respectively on CIFAR-10; and 43% and 35% on CIFAR-100, respectively.


title: A focus on the Riemann's Hypothesis

abstract: Riemann's hypothesis, formulated in 1859, concerns the location of the zeros of Riemann's Zeta function. The history of the Riemann hypothesis is well known. In 1859, the German mathematician B. Riemann presented a paper to the Berlin Academy of Mathematic. In that paper, he proposed that this function, called Riemann-zeta function takes values 0 on the complex plane when s=0.5+it. This hypothesis has great significance for the world of mathematics and physics. This solutions would lead to innumerable completions of theorems that rely upon its truth. Over a billion zeros of the function have been calculated by computers and shown that all are on this line s = 0.5+it. In this paper, we initially show that Riemann's (Z\^eta) function and the analytical extension of this function called (Aleph)) are distinct. After extending this function in the complex plane except the point s=1, we will show the existence and then the uniqueness of real part zeros equal to 1/2.


title: Impact of weather factors on migration intention using machine learning algorithms

abstract: A growing attention in the empirical literature has been paid to the incidence of climate shocks and change in migration decisions. Previous literature leads to different results and uses a multitude of traditional empirical approaches. This paper proposes a tree-based Machine Learning (ML) approach to analyze the role of the weather shocks towards an individual's intention to migrate in the six agriculture-dependent-economy countries such as Burkina Faso, Ivory Coast, Mali, Mauritania, Niger, and Senegal. We perform several tree-based algorithms (e.g., XGB, Random Forest) using the train-validation-test workflow to build robust and noise-resistant approaches. Then we determine the important features showing in which direction they are influencing the migration intention. This ML-based estimation accounts for features such as weather shocks captured by the Standardized Precipitation-Evapotranspiration Index (SPEI) for different timescales and various socioeconomic features/covariates. We find that (i) weather features improve the prediction performance although socioeconomic characteristics have more influence on migration intentions, (ii) country-specific model is necessary, and (iii) international move is influenced more by the longer timescales of SPEIs while general move (which includes internal move) by that of shorter timescales.


title: Simulating hydrodynamics on noisy intermediate-scale quantum devices with random circuits

abstract: In a recent milestone experiment, Google's processor Sycamore heralded the era of "quantum supremacy" by sampling from the output of (pseudo-)random circuits. We show that such random circuits provide tailor-made building blocks for simulating quantum many-body systems on noisy intermediate-scale quantum (NISQ) devices. Specifically, we propose an algorithm consisting of a random circuit followed by a trotterized Hamiltonian time evolution to study hydrodynamics and to extract transport coefficients in the linear response regime. We numerically demonstrate the algorithm by simulating the buildup of spatiotemporal correlation functions in one- and two-dimensional quantum spin systems, where we particularly scrutinize the inevitable impact of errors present in any realistic implementation. Importantly, we find that the hydrodynamic scaling of the correlations is highly robust with respect to the size of the Trotter step, which opens the door to reach nontrivial time scales with a small number of gates. While errors within the random circuit are shown to be irrelevant, we furthermore unveil that meaningful results can be obtained for noisy time evolutions with error rates achievable on near-term hardware. Our work emphasizes the practical relevance of random circuits on NISQ devices beyond the abstract sampling task.


title: Fast Optical Transients from Stellar-Mass Black Hole Tidal Disruption Events in Young Star Clusters

abstract: Observational evidence suggests that the majority of stars may have been born in stellar clusters or associations. Within these dense environments, dynamical interactions lead to high rates of close stellar encounters. A variety of recent observational and theoretical indications suggest stellar-mass black holes may be present and play an active dynamical role in stellar clusters of all masses. In this study, we explore the tidal disruption of main sequence stars by stellar-mass black holes in young star clusters. We compute a suite of over 3000 independent $N$-body simulations that cover a range in cluster mass, metallicity, and half-mass radii. We find stellar-mass black hole tidal disruption events (TDEs) occur at an overall rate of up to roughly $200\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ in young stellar clusters in the local universe. These TDEs are expected to have several characteristic features, namely fast rise times of order a day, peak X-ray luminosities of at least $10^{44}\,\rm{erg\,s}^{-1}$, and bright optical luminosities (roughly $10^{41}-10^{44}\,\rm{erg\,s}^{-1}$) associated with reprocessing by a disk wind. In particular, we show these events share many features in common with the emerging class of Fast Blue Optical Transients.


title: Retarded Green's Function from Rotating AdS Black Holes: Emergent CFT$_2$ and Viscosity

abstract: Using the AdS/CFT correspondence we consider the retarded Green's function in the background of rotating near-extremal AdS$_4$ black holes. Following the canonical AdS/CFT dictionary into the asymptotic boundary we get a CFT$_3$ result. We also take a new route and zoom in on the near-horizon region, blow up this region and show that it yields a CFT$_2$ result. We argue that the decoupling of the near-horizon region is akin to the decoupling of the near-throat region of a D3-brane, which led to the original formulation of the AdS/CFT correspondence, thus implying that the Kerr/CFT correspondence follows as a decoupling of the standard AdS/CFT correspondence applied to rotating black holes. As a byproduct, we compute the shear viscosity to entropy density ratio for the strongly coupled boundary CFT$_3$, and find that it violates the $1 / (4 \pi)$ bound.


title: A medium-resolution spectrum of the exoplanet HIP 65426 b

abstract: Medium-resolution integral-field spectrographs (IFS) coupled with adaptive-optics such as Keck/OSIRIS, VLT/MUSE, or SINFONI are appearing as a new avenue for enhancing the detection and characterization capabilities of young, gas giant exoplanets at large heliocentric distances (>5 au). We analyzed K-band VLT/SINFONI medium-resolution (R_lambda~5577) observations of the young giant exoplanet HIP 65426 b. Our dedicated IFS data analysis toolkit (TExTRIS) optimized the cube building, star registration, and allowed for the extraction of the planet spectrum. A Bayesian inference with the nested sampling algorithm coupled with the self-consistent forward atmospheric models BT-SETTL15 and Exo-REM using the ForMoSA tool yields Teff=1560 +/- 100K, log(g)<4.40dex, [M/H]=0.05 +/- 0.22dex, and an upper limit on the C/O ratio (<0.55). The object is also re-detected with the so-called "molecular mapping" technique. The technique yields consistent atmospheric parameters, but the loss of the planet pseudo-continuum in the process degrades or modifies the constraints on these parameters. The solar to sub-solar C/O ratio suggests an enrichment by solids at formation if the planet was formed beyond the water snowline (>20 au) by core-accretion. However, a formation by gravitational instability can not be ruled out. The metallicity is compatible with the bulk enrichment of massive Jovian planets from the Bern planet population models. Finally, we measure a radial velocity of 26 +/- 15km/s compatible with our revised measurement on the star. This is the fourth imaged exoplanet for which a radial velocity can be evaluated, illustrating the potential of such observations for assessing the coevolution of imaged systems belonging to star forming regions, such as HIP 65426.


title: Left-right symmetry and electric dipole moments. A global analysis

abstract: We perform a global fit using results of searches for electric dipole moments (EDM) of diamagnetic systems within the context of the minimal left-right symmetric model. In this way, we disentangle the new "left-right" electroweak and $\bar{\theta}$ contributions that cannot be separated using a single EDM system. Although the fit is done for a specific model, the approach can be applied to any particle physics model. Finally, we revisit the constraint on the $D$ coefficient in $\beta$-decay and find that current EDM bounds do not preclude observation of this T-violating effect in a possible next generation $\beta$-decay experiment.


title: Evolution of Galaxy Star Formation and Metallicity: Impact on Double Compact Objects Mergers

abstract: We study the impact of different galaxy statistics and empirical metallicity scaling relations on the merging rates and on the properties of compact objects binaries. First, we analyze the similarities and differences of using the star formation rate functions or the stellar mass functions as galaxy statistics for the computation of the cosmic star formation rate density. Then we investigate the effects of adopting the Fundamental Metallicity Relation or a classic Mass Metallicity Relation to assign metallicity to galaxies with given properties. We find that when the Fundamental Metallicity Relation is exploited, the bulk of the star formation occurs at relatively high metallicities even at high redshift; the opposite holds when the Mass Metallicity Relation is employed, since in this case the metallicity at which most of the star formation takes place strongly decreases with redshift. We discuss the various reasons and possible biases originating this discrepancy. Finally, we show the impact that these different astrophysical prescriptions have on the merging rates and on the properties of compact objects binaries; specifically, we present results for the redshift dependent merging rates and for the chirp mass and time delay distributions of the merging binaries.


title: Capacity of a lossy photon channel with direct detection

abstract: We calculate numerically the capacity of a lossy photon channel assuming photon number resolving detection at the output. We consider scenarios of input Fock and coherent states ensembles and show that the latter always exhibits worse performance than the former. We obtain capacity of a discrete-time Poisson channel as a limiting behavior of the Fock states ensemble capacity. We show also that in the regime of a moderate number of photons and low losses the Fock states ensemble with direct detection is beneficial with respect to capacity limits achievable with quadrature detection.


title: SLITronomy: towards a fully wavelet-based strong lensing inversion technique

abstract: Strong gravitational lensing provides a wealth of astrophysical information on the baryonic and dark matter content of galaxies. It also serves as a valuable cosmological probe by allowing us to measure the Hubble constant independently of other methods. These applications all require the difficult task of inverting the lens equation and simultaneously reconstructing the mass profile of the lens along with the original light profile of the unlensed source. As there is no reason for either the lens or the source to be simple, we need methods that both invert the lens equation with a large number of degrees of freedom and also enforce a well-controlled regularisation that avoids the appearance of spurious structures. This can be beautifully accomplished by representing signals in wavelet space. Building on the Sparse Lens Inversion Technique (SLIT), in this work we present an improved sparsity-based method that describes lensed sources using wavelets and optimises over the parameters given an analytical lens mass profile. We apply our technique on simulated HST and E-ELT data, as well as on real HST images of lenses from the Sloan Lens ACS (SLACS) sample, assuming a lens model. We show that wavelets allow us to reconstruct lensed sources containing detailed substructures when using both present-day data and high-resolution images from future thirty-meter-class telescopes. Wavelets moreover provide a much more tractable solution in terms of quality and computation time compared to using a source model that combines smooth analytical profiles and shapelets. Requiring very little human interaction, our pixel-based technique fits into the effort to devise automated modelling schemes. It can be incorporated in the standard workflow of sampling analytical lens model parameters. The method, which we call SLITronomy, is freely available as a new plug-in to the modelling software Lenstronomy.


title: Quantum Quasi-Monte Carlo algorithm for out-of-equilibrium Green functions at long times

abstract: We extend the recently developed Quantum Quasi-Monte Carlo (QQMC) approach to obtain the full frequency dependence of Green functions in a single calculation. QQMC is a general approach for calculating high-order perturbative expansions in power of the electron-electron interaction strength. In contrast to conventional Markov chain Monte Carlo sampling, QQMC uses low-discrepancy sequences for a more uniform sampling of the multi-dimensional integrals involved and can potentially outperform Monte Carlo by several orders of magnitudes. A core concept of QQMC is the a priori construction of a "model function" that approximates the integrand and is used to optimize the sampling distribution. In this paper, we show that the model function concept extends to a kernel approach for the computation of Green functions. We illustrate the approach on the Anderson impurity model and show that the scaling of the error with the number of integrand evaluations $N$ is $\sim 1/N^{0.86}$ in the best cases, and comparable to Monte Carlo scaling $\sim 1/N^{0.5}$ in the worst cases. We find a systematic improvement over Monte Carlo sampling by at least two orders of magnitude while using a basic form of model function. Finally, we compare QQMC results with calculations performed with the Fork Tensor Product State (FTPS) method, a recently developed tensor network approach for solving impurity problems. Applying a simple Pad\'e approximant for the series resummation, we find that QQMC matches the FTPS results beyond the perturbative regime.


title: Group testing for overlapping communities

abstract: In this paper, we propose algorithms that leverage a known community structure to make group testing more efficient. We consider a population organized in connected communities: each individual participates in one or more communities, and the infection probability of each individual depends on the communities (s)he participates in. Use cases include students who participate in several classes, and workers who share common spaces. Group testing reduces the number of tests needed to identify the infected individuals by pooling diagnostic samples and testing them together. We show that making testing algorithms aware of the community structure, can significantly reduce the number of tests needed both for adaptive and non-adaptive group testing.


title: Adaptive Explicit Kernel Minkowski Weighted K-means

abstract: The K-means algorithm is among the most commonly used data clustering methods. However, the regular K-means can only be applied in the input space and it is applicable when clusters are linearly separable. The kernel K-means, which extends K-means into the kernel space, is able to capture nonlinear structures and identify arbitrarily shaped clusters. However, kernel methods often operate on the kernel matrix of the data, which scale poorly with the size of the matrix or suffer from the high clustering cost due to the repetitive calculations of kernel values. Another issue is that algorithms access the data only through evaluations of $K(x_i, x_j)$, which limits many processes that can be done on data through the clustering task. This paper proposes a method to combine the advantages of the linear and nonlinear approaches by using driven corresponding approximate finite-dimensional feature maps based on spectral analysis. Applying approximate finite-dimensional feature maps were only discussed in the Support Vector Machines (SVM) problems before. We suggest using this method in kernel K-means era as alleviates storing huge kernel matrix in memory, further calculating cluster centers more efficiently and access the data explicitly in feature space. These explicit feature maps enable us to access the data in the feature space explicitly and take advantage of K-means extensions in that space. We demonstrate our Explicit Kernel Minkowski Weighted K-mean (Explicit KMWK-mean) method is able to be more adopted and find best-fitting values in new space by applying additional Minkowski exponent and feature weights parameter. Moreover, it can reduce the impact of concentration on nearest neighbour search by suggesting investigate among other norms instead of Euclidean norm, includes Minkowski norms and fractional norms (as an extension of the Minkowski norms with p<1).


title: Policy Maker's Credibility with Predetermined Instruments for Forward-Looking Targets

abstract: The aim of the present paper is to provide criteria for a central bank of how to choose among different monetary-policy rules when caring about a number of policy targets such as the output gap and expected inflation. Special attention is given to the question if policy instruments are predetermined or only forward looking. Using the new-Keynesian Phillips curve with a cost-push-shock policy-transmission mechanism, the forward-looking case implies an extreme lack of robustness and of credibility of stabilization policy. The backward-looking case is such that the simple-rule parameters can be the solution of Ramsey optimal policy under limited commitment. As a consequence, we suggest to model explicitly the rational behavior of the policy maker with Ramsey optimal policy, rather than to use simple rules with an ambiguous assumption leading to policy advice that is neither robust nor credible.


title: Learning summary features of time series for likelihood free inference

abstract: There has been an increasing interest from the scientific community in using likelihood-free inference (LFI) to determine which parameters of a given simulator model could best describe a set of experimental data. Despite exciting recent results and a wide range of possible applications, an important bottleneck of LFI when applied to time series data is the necessity of defining a set of summary features, often hand-tailored based on domain knowledge. In this work, we present a data-driven strategy for automatically learning summary features from univariate time series and apply it to signals generated from autoregressive-moving-average (ARMA) models and the Van der Pol Oscillator. Our results indicate that learning summary features from data can compete and even outperform LFI methods based on hand-crafted values such as autocorrelation coefficients even in the linear case.


title: Persistent Laplacians: properties, algorithms and implications

abstract: We present a thorough study of the theoretical properties and devise efficient algorithms for the \emph{persistent Laplacian}, an extension of the standard combinatorial Laplacian to the setting of pairs (or, in more generality, sequences) of simplicial complexes $K \hookrightarrow L$, which was independently introduced by Lieutier et al. and by Wang et al. In particular, in analogy with the non-persistent case, we first prove that the nullity of the $q$-th persistent Laplacian $\Delta_q^{K,L}$ equals the $q$-th persistent Betti number of the inclusion $(K \hookrightarrow L)$. We then present an initial algorithm for finding a matrix representation of $\Delta_q^{K,L}$, which itself helps interpret the persistent Laplacian. We exhibit a novel relationship between the persistent Laplacian and the notion of Schur complement of a matrix which has several important implications. In the graph case, it both uncovers a link with the notion of effective resistance and leads to a persistent version of the Cheeger inequality. This relationship also yields an additional, very simple algorithm for finding (a matrix representation of) the $q$-th persistent Laplacian which in turn leads to a novel and fundamentally different algorithm for computing the $q$-th persistent Betti number for a pair $(K,L)$ which can be significantly more efficient than standard algorithms. Finally, we study persistent Laplacians for simplicial filtrations and present novel stability results for their eigenvalues. Our work brings methods from spectral graph theory, circuit theory, and persistent homology together with a topological view of the combinatorial Laplacian on simplicial complexes.


title: ACN-Sim: An Open-Source Simulator for Data-Driven Electric Vehicle Charging Research

abstract: ACN-Sim is a data-driven, open-source simulation environment designed to accelerate research in the field of smart electric vehicle (EV) charging. It fills the need in this community for a widely available, realistic simulation environment in which researchers can evaluate algorithms and test assumptions. ACN-Sim provides a modular, extensible architecture, which models the complexity of real charging systems, including battery charging behavior and unbalanced three-phase infrastructure. It also integrates with a broader ecosystem of research tools. These include ACN-Data, an open dataset of EV charging sessions, which provides realistic simulation scenarios and ACN-Live, a framework for field-testing charging algorithms. It also integrates with grid simulators like MATPOWER, PandaPower and OpenDSS, and OpenAI Gym for training reinforcement learning agents.


title: Effects of hybridization and spin-orbit coupling to induce odd frequency pairing in two-band superconductors

abstract: The effects of spin independent hybridization potential and spin orbit coupling on two band superconductor with equal time s-wave inter band pairing order parameter is investigated theoretically. To study symmetry classes in two band superconductors the Gorkov equations are solved analytically. By defining spin singlet and spin triplet s wave order parameter due to two band degree of freedom the symmetry classes of Cooper pair are studied. For spin singlet case it is shown that spin independent hybridization generates Cooper pair belongs to even frequency spin singlet even momentum even band parity (ESEE) symmetry class for both intraband and interband pairing correlations. For spin triplet order parameter, intraband pairing correlation generates odd frequency spin triplet even momentum even band parity (OTEE) symmetry class whereas, interband pairing correlation generates even frequency spin triplet even momentum odd band parity ETEO) class. For the spin singlet, spin orbit coupling generates pairing correlation that belongs to odd frequency spin singlet odd momentum even band parity (OSOE) symmetry class and even frequency spin singlet even momentum even band parity (ESEE) for intraband and interband pairing correlation respectively. In the spin triplet case for itraband and interband correlation, spin orbit coupling generates even-frequency spin triplet odd momentum even band parity (ETOE) and even frequency spin triplet even momentum odd band parity (ETEO) respectively.


title: Modeling Voters in Multi-Winner Approval Voting

abstract: In many real world situations, collective decisions are made using voting and, in scenarios such as committee or board elections, employing voting rules that return multiple winners. In multi-winner approval voting (AV), an agent submits a ballot consisting of approvals for as many candidates as they wish, and winners are chosen by tallying up the votes and choosing the top-$k$ candidates receiving the most approvals. In many scenarios, an agent may manipulate the ballot they submit in order to achieve a better outcome by voting in a way that does not reflect their true preferences. In complex and uncertain situations, agents may use heuristics instead of incurring the additional effort required to compute the manipulation which most favors them. In this paper, we examine voting behavior in single-winner and multi-winner approval voting scenarios with varying degrees of uncertainty using behavioral data obtained from Mechanical Turk. We find that people generally manipulate their vote to obtain a better outcome, but often do not identify the optimal manipulation. There are a number of predictive models of agent behavior in the COMSOC and psychology literature that are based on cognitively plausible heuristic strategies. We show that the existing approaches do not adequately model real-world data. We propose a novel model that takes into account the size of the winning set and human cognitive constraints, and demonstrate that this model is more effective at capturing real-world behaviors in multi-winner approval voting scenarios.


title: Quantum-Assisted Optical Interferometers: Instrument Requirements

abstract: It has been recently suggested that optical interferometers may not require a phase-stable optical link between the stations if instead sources of quantum-mechanically entangled pairs could be provided to them, enabling extra-long baselines and benefiting numerous topics in astrophysics and cosmology. We developed a new variation of this idea, proposing that photons from two different sources could be interfered at two decoupled stations, requiring only a slow classical connection between them. We show that this approach could allow high-precision measurements of the relative astrometry of the two sources, with a simple estimate giving angular resolution of $10 \ \mu$as in a few hours' observation of two bright stars. We also give requirements on the instrument for these observations, in particular on its temporal and spectral resolution. Finally, we discuss possible technologies for the instrument implementation and first proof-of-principle experiments.


title: Cross-Modal Generalization: Learning in Low Resource Modalities via Meta-Alignment

abstract: The natural world is abundant with concepts expressed via visual, acoustic, tactile, and linguistic modalities. Much of the existing progress in multimodal learning, however, focuses primarily on problems where the same set of modalities are present at train and test time, which makes learning in low-resource modalities particularly difficult. In this work, we propose algorithms for cross-modal generalization: a learning paradigm to train a model that can (1) quickly perform new tasks in a target modality (i.e. meta-learning) and (2) doing so while being trained on a different source modality. We study a key research question: how can we ensure generalization across modalities despite using separate encoders for different source and target modalities? Our solution is based on meta-alignment, a novel method to align representation spaces using strongly and weakly paired cross-modal data while ensuring quick generalization to new tasks across different modalities. We study this problem on 3 classification tasks: text to image, image to audio, and text to speech. Our results demonstrate strong performance even when the new target modality has only a few (1-10) labeled samples and in the presence of noisy labels, a scenario particularly prevalent in low-resource modalities.


title: Dynamical environments of MU69 and similar objects

abstract: We explore properties of the long-term dynamics of particles (moonlets, fragments, debris etc.) around KBO 2014 MU69 (Arrokoth), as well as around similar contact-binary objects potentially present in the Kuiper belt. The chaotic diffusion of particles inside the Hill sphere of MU69 (or, generally, a similar object) is studied by means of construction of appropriate stability diagrams and by application of analytical approaches generally based on the Kepler map theory. The formation and evolution of the particle clouds, due to the chaotic diffusion inside the Hill sphere, are studied and the cloud lifetimes are estimated.


title: Classification of affine normal $SL_2$-varieties with a dense orbit

abstract: In this note we give a classification of all affine normal $SL_2$-varieties containing an open dense orbit over an algebraically closed field of characteristic zero. Such a classification was first obtained by Popov. Here we provide an alternative approach.


title: MCMC Confidence Intervals and Biases

abstract: The recent paper "Simple confidence intervals for MCMC without CLTs" by J.S. Rosenthal, showed the derivation of a simple MCMC confidence interval using only Chebyshev's inequality, not CLT. That result required certain assumptions about how the estimator bias and variance grow with the number of iterations $n$. In particular, the bias is $o(1/\sqrt{n})$. This assumption seemed mild. It is generally believed that the estimator bias will be $O(1/n)$ and hence $o(1/\sqrt{n})$. However, questions were raised by researchers about how to verify this assumption. Indeed, we show that this assumption might not always hold. In this paper, we seek to simplify and weaken the assumptions in the previously mentioned paper, to make MCMC confidence intervals without CLTs more widely applicable.


title: Examining the end of inflation with primordial black holes mass distribution and gravitational waves

abstract: We explicitly construct a double-field inflationary model, which satisfies the latest Planck constraints at the cosmic microwave background (CMB) scales and produces the whole dark matter energy density as primordial black holes (PBHs), in the mass range $10^{-17}~M_{\odot}\lesssim M_{{}_{\rm PBH}}\lesssim 10^{-13}~M_{\odot}$. The PBHs can be produced after the end of slow-roll inflation from the bubbles of true vacuum that nucleate during the course of inflation. Obtaining PBHs in this mass range enforces the scale of inflation to be extremely low, $10^{-7} \lesssim H \lesssim 10^{-3} ~{\rm GeV}$, which makes the efforts to observe gravitational waves at the CMB scales futile, although it is high enough to allow for a successful big bang nucleosynthesis. We will show that the shape of the mass distribution of the PBHs is dependent on how inflation ends and the Universe settles from the metastable direction to the true one. End of inflation can also be probed by examining the gravitational waves spectrum. In particular, we show that if exit from the rolling metastable direction to the true vacuum of the potential happens through a first-order phase transition after the end of slow-roll inflation, it leaves behind a stochastic gravitational wave background (SGWB), which is potentially observable by the Laser Interferometer Space Antenna. Examining the mass distribution of PBHs and possible SGWB from the end of inflation, we may be able to gain invaluable information about the end of inflation.


title: Encoding the latent posterior of Bayesian Neural Networks for uncertainty quantification

abstract: Bayesian neural networks (BNNs) have been long considered an ideal, yet unscalable solution for improving the robustness and the predictive uncertainty of deep neural networks. While they could capture more accurately the posterior distribution of the network parameters, most BNN approaches are either limited to small networks or rely on constraining assumptions such as parameter independence. These drawbacks have enabled prominence of simple, but computationally heavy approaches such as Deep Ensembles, whose training and testing costs increase linearly with the number of networks. In this work we aim for efficient deep BNNs amenable to complex computer vision architectures, e.g. ResNet50 DeepLabV3+, and tasks, e.g. semantic segmentation, with fewer assumptions on the parameters. We achieve this by leveraging variational autoencoders (VAEs) to learn the interaction and the latent distribution of the parameters at each network layer. Our approach, Latent-Posterior BNN (LP-BNN), is compatible with the recent BatchEnsemble method, leading to highly efficient ({in terms of computation and} memory during both training and testing) ensembles. LP-BNN s attain competitive results across multiple metrics in several challenging benchmarks for image classification, semantic segmentation and out-of-distribution detection.


title: On-Device Sentence Similarity for SMS Dataset

abstract: Determining the sentence similarity between Short Message Service (SMS) texts/sentences plays a significant role in mobile device industry. Gauging the similarity between SMS data is thus necessary for various applications like enhanced searching and navigation, clubbing together SMS of similar type when given a custom label or tag is provided by user irrespective of their sender etc. The problem faced with SMS data is its incomplete structure and grammatical inconsistencies. In this paper, we propose a unique pipeline for evaluating the text similarity between SMS texts. We use Part of Speech (POS) model for keyword extraction by taking advantage of the partial structure embedded in SMS texts and similarity comparisons are carried out using statistical methods. The proposed pipeline deals with major semantic variations across SMS data as well as makes it effective for its application on-device (mobile phone). To showcase the capabilities of our work, our pipeline has been designed with an inclination towards one of the possible applications of SMS text similarity discussed in one of the following sections but nonetheless guarantees scalability for other applications as well.


title: Signatures for Berezinsky-Kosterlitz-Thouless critical behaviour in the planar antiferromagnet BaNi$_2$V$_2$O$_8$

abstract: We investigate the critical properties of the spin-$1$ honeycomb antiferromagnet BaNi$_2$V$_2$O$_8$, both below and above the ordering temperature $T_N$ using neutron diffraction and muon spin rotation measurements. Our results characterize BaNi$_2$V$_2$O$_8$ as a two-dimensional (2D) antiferromagnet across the entire temperature range, displaying a series of crossovers from 2D Ising-like to 2D XY and then to 2D Heisenberg behavior with increasing temperature. In particular, the extracted critical exponent of the order parameter reveals a narrow temperature regime close to $T_N$, in which the system behaves as a 2D XY antiferromagnet. Above $T_N$, evidence for Berezinsky-Kosterlitz-Thouless behavior driven by vortex excitations is obtained from the scaling of the correlation length. Our experimental results are in accord with classical and quantum Monte Carlo simulations performed for microscopic magnetic model Hamiltonians for BaNi$_2$V$_2$O$_8$.


title: MPG: A Multi-ingredient Pizza Image Generator with Conditional StyleGANs

abstract: Multilabel conditional image generation is a challenging problem in computer vision. In this work we propose Multi-ingredient Pizza Generator (MPG), a conditional Generative Neural Network (GAN) framework for synthesizing multilabel images. We design MPG based on a state-of-the-art GAN structure called StyleGAN2, in which we develop a new conditioning technique by enforcing intermediate feature maps to learn scalewise label information. Because of the complex nature of the multilabel image generation problem, we also regularize synthetic image by predicting the corresponding ingredients as well as encourage the discriminator to distinguish between matched image and mismatched image. To verify the efficacy of MPG, we test it on Pizza10, which is a carefully annotated multi-ingredient pizza image dataset. MPG can successfully generate photo-realist pizza images with desired ingredients. The framework can be easily extend to other multilabel image generation scenarios.


title: Interferometric 4D-STEM for Lattice Distortion and Interlayer Spacing Measurements in Bilayer and Trilayer Two-dimensional Materials

abstract: Van der Waals materials composed of stacks of individual atomic layers have attracted considerable attention due to their exotic electronic properties that can be altered by, for example, manipulating the twist angle of bilayer materials or the stacking sequence of trilayer materials. To fully understand and control the unique properties of these few-layer materials, a technique that can provide information about their local in-plane structural deformations, twist direction, and out-of-plane structure is needed. In principle, interference in overlap regions of Bragg disks originating from separate layers of a material encodes three-dimensional information about the relative positions of atoms in the corresponding layers. Here, we describe an interferometric four-dimensional scanning transmission electron microscopy technique that utilizes this phenomenon to extract precise structural information from few-layer materials with nm-scale resolution. We demonstrate how this technique enables measurement of local pm-scale in-plane lattice distortions as well as twist direction and average interlayer spacings in bilayer and trilayer graphene, and therefore provides a means to better understand the interplay between electronic properties and precise structural arrangements of few-layer 2D materials.


title: Gluing Noncommutative Twistor Spaces

abstract: We describe a general procedure, based on Gerstenhaber-Schack complexes, for extending to quantized twistor spaces the Donaldson-Friedman gluing of twistor spaces via deformation theory of singular spaces. We consider in particular various possible quantizations of twistor spaces that leave the underlying spacetime manifold classical, including the geometric quantization of twistor spaces originally constructed by the second author, as well as some variants based on noncommutative geometry. We discuss specific aspects of the gluing construction for these different quantization procedures.


title: High order steady-state diffusion approximations

abstract: We derive and analyze new diffusion approximations of stationary distributions of Markov chains that are based on second- and higher-order terms in the expansion of the Markov chain generator. Our approximations achieve a higher degree of accuracy compared to diffusion approximations widely used for the past fifty years, while retaining a similar computational complexity. To support our approximations, we present a combination of theoretical and numerical results across three different models. Our approximations are derived recursively through Stein/Poisson equations, and the theoretical results are proved using Stein's method.


title: A Survey on Deep Learning for Human Mobility

abstract: The study of human mobility is crucial due to its impact on several aspects of our society, such as disease spreading, urban planning, well-being, pollution, and more. The proliferation of digital mobility data, such as phone records, GPS traces, and social media posts, combined with the predictive power of artificial intelligence, triggered the application of deep learning to human mobility. Existing surveys focus on single tasks, data sources, mechanistic or traditional machine learning approaches, while a comprehensive description of deep learning solutions is missing. This survey provides a taxonomy of mobility tasks, a discussion on the challenges related to each task and how deep learning may overcome the limitations of traditional models, a description of the most relevant solutions to the mobility tasks described above and the relevant challenges for the future. Our survey is a guide to the leading deep learning solutions to next-location prediction, crowd flow prediction, trajectory generation, and flow generation. At the same time, it helps deep learning scientists and practitioners understand the fundamental concepts and the open challenges of the study of human mobility.


title: Numerical solution of a time-fractional nonlinear Rayleigh-Stokes problem

abstract: We study a semilinear fractional-in-time Rayleigh-Stokes problem for a generalized second-grade fluid with a Lipschitz continuous nonlinear source term and initial data $u_0\in\dot{H}^\nu(\Omega)$, $\nu\in[0,2]$. We discuss stability of solutions and provide regularity results. Two spatially semidiscrete schemes are analyzed based on standard Galerkin and lumped mass finite element methods, respectively. Further, a fully discrete scheme is obtained by applying a convolution quadrature in time generated by the backward Euler method, and optimal error estimates are derived for smooth and nonsmooth initial data. Finally, numerical examples are provided to illustrate the theoretical results.


title: $Z_1$-oscillation in the retarding force of metals for slow ions: Comparative study of theoretical modelings

abstract: Theoretical calculations, applying elements of the orbital-based density functional method for self-consistent screening in an electron gas to the stopping power of metals for slow projectiles are reviewed, by focusing on differences in modelings. New, two-channel-based results on the average retarding force are presented. All theoretical results are compared with pioneering experimental data, obtained for carbon target, on $Z_1$-oscillation in stopping at random collisional condition.


title: Automatic image-based respiratory signal extraction in real-time CMR

abstract: Purpose: To develop a fully automatic method for extraction and directionality determination of respiratory signal in free-breathing, real-time (RT) cardiac MRI. Methods: The respiratory signal is extracted by a principal component analysis method from RT cine images. Then, a two-step procedure is used to determine the directionality (sign) of the respiratory signal. First, the signal polarity of all slices is made consistent with a reference slice. Second, a global sign correction is performed by maximizing the correlation of the respiratory signal with the zeroth-moment center curve. The proposed method is evaluated in multi-slice RT cine from eleven volunteers and two patients. The motion in a manually selected region-of-interest is used as reference. Results: The extracted respiratory signal using the proposed method exhibits high, positive correlation with the reference in all cases and is more robust compared to a recently proposed method. To demonstrate its clinical utility, the method is used to identify heartbeats from inspiratory and expiratory phases and, in turn, to account for respiration-induced changes in the cardiac output. Conclusions: The proposed method enables fully automatic extraction and directionality determinations of respiratory signal from RT cardiac cine images, allowing accurate cardiac function quantification.


title: A Small Satellite Version of a Broad-band Soft X-ray Polarimeter

abstract: We describe a new implementation of a broad-band soft X-ray polarimeter, substantially based on a previous design. This implementation, the Pioneer Soft X-ray Polarimeter (PiSoX) is a SmallSat, designed for NASA's call for Astrophysics Pioneers, small missions that could be CubeSats, balloon experiments, or SmallSats. As in the REDSoX Polarimeter, the grating arrangement is designed optimally for the purpose of polarimetry with broad-band focussing optics by matching the dispersion of the spectrometer channels to laterally graded multilayers (LGMLs). The system can achieve polarization modulation factors over 90%. For PiSoX, the optics are lightweight Si mirrors in a one-bounce parabolic configuration. High efficiency, blazed gratings from opposite sectors are oriented to disperse to a LGML forming a channel covering the wavelength range from 35 to 75 Angstroms (165 - 350 eV). Upon satellite rotation, the intensities of the dispersed spectra, after reflection and polarizing by the LGMLs, give the three Stokes parameters needed to determine a source's linear polarization fraction and orientation. The design can be extended to higher energies as LGMLs are developed further. We describe examples of the potential scientific return from instruments based on this design.


title: Simultaneous averaging to zero by unitary mixing operators

abstract: We show that if every element a vector subspace of a C*-algebra can be averaged to zero by means of unitary mixing operators, then all the elements of the subspace can be simultaneously averaged to zero by a net of unitary mixing operators. Moreover, such subspaces admit a simple description in terms of commutators and kernels of states on the C*-algebra. We apply this result to center-valued expectations in C*-algebras with the Dixmier property.


title: Stochastic processes and host-parasite coevolution: linking coevolutionary dynamics and DNA polymorphism data

abstract: Between-species coevolution, and in particular antagonistic host-parasite coevolution, is a major process shaping within-species diversity. In this paper we investigate the role of various stochastic processes affecting the outcome of the deterministic coevolutionary models. Specifically, we assess 1) the impact of genetic drift and mutation on the maintenance of polymorphism at the interacting loci, and 2) the change in neutral allele frequencies across the genome of both coevolving species due to co-demographic population size changes. We find that genetic drift decreases the likelihood to observe classic balancing selection signatures, and that for most realistic values of the coevolutionary parameters, balancing selection signatures cannot be seen at the host loci. Further, we reveal that contrary to classic expectations, fast changes in parasite population size due to eco-evo feedbacks can be tracked by the allelic site-frequency spectrum measured at several time points. Changes in host population size are, however, less pronounced and thus not observable. Finally, we also review several understudied stochastic processes occurring in host-parasite coevolution which are of importance to predict maintenance of polymorphism at the underlying loci and the genome-wide nucleotide diversity of host and parasite populations.


title: A software decoder implementation for H.266/VVC video coding standard

abstract: Versatile Video Coding Standard (H.266/VVC) was completed by Joint Video Expert Team (JVET) of ITU-T and ISO/IEC, in July 2020. This new ITU recommendation/international standard is a successor to the well-known H.265/HEVC video coding standard with roughly doubled compression efficiency, but also at the cost of an increased computational complexity. The complexity of H.266/VVC decoder processing modules is studied in this paper. An optimized decoder implementation using SIMD instruction extensions and additional parallel processing including data and task level parallelism is presented, which can achieve real-time decoding of 4K 60fps VVC bitstreams on an x86 based CPU.


title: An epidemic model for COVID-19 transmission in Argentina: Exploration of the alternating quarantine and massive testing strategies

abstract: The COVID-19 pandemic has challenged authorities at different levels of government administration around the globe. When faced with diseases of this severity, it is useful for the authorities to have prediction tools to estimate in advance the impact on the health system and the human, material, and economic resources that will be necessary. In this paper, we construct an extended Susceptible-Exposed-Infected-Recovered model that incorporates the social structure of Mar del Plata, the $4^\circ$ most inhabited city in Argentina and head of the Municipality of General Pueyrred\'on. Moreover, we consider detailed partitions of infected individuals according to the illness severity, as well as data of local health resources, to bring these predictions closer to the local reality. Tuning the corresponding epidemic parameters for COVID-19, we study an alternating quarantine strategy, in which a part of the population can circulate without restrictions at any time, while the rest is equally divided into two groups and goes on successive periods of normal activity and lockdown, each one with a duration of $\tau$ days. Besides, we implement a random testing strategy over the population. We found that $\tau = 7$ is a good choice for the quarantine strategy since it matches with the weekly cycle as it reduces the infected population. Focusing on the health system, projecting from the situation as of September 30, we foresee a difficulty to avoid saturation of ICU, given the extremely low levels of mobility that would be required. In the worst case, our model estimates that four thousand deaths would occur, of which 30\% could be avoided with proper medical attention. Nonetheless, we found that aggressive testing would allow an increase in the percentage of people that can circulate without restrictions, being the equipment required to deal with the additional critical patients relatively low.


title: Light scattering by a metallic nanoparticle coated with a nematic liquid crystal

abstract: We study the optical properties of gold nanoparticles coated with a nematic liquid crystal whose director field is distributed around the nanoparticle according to the anchoring conditions at the surface of the nanoparticle. The distribution of the nematic liquid crystal is obtained by minimization of the corresponding Frank free-energy functional whilst the optical response is calculated by the discrete-dipole approximation. We find, in particular, that the anisotropy of the nematic liquid-crystal coating does not affect much the (isotropic) optical response of the nanoparticle. However, for strong anchoring of the nematic liquid-crystal molecules on the surface of nanoparticle, the inhomogeneity of the coating which is manifested by a ring-type singularity (disclination or Saturn ring), produces an enhancement of the extinction cross spectrum over the entire visible spectrum.


title: Chaotic dynamics in the presence of medical malpractice litigation: a topological proof via linked twist maps for two evolutionary game theoretic contexts

abstract: In the present work we reconsider the evolutionary game theoretic models by Antoci et al. (2016, 2018) describing the dynamic outcomes arising from the interactions between patients and physicians, whose behavior is subject to clinical and legal risks. In particular, Antoci et al. (2016) analyzed the case of positive defensive medicine, while Antoci et al. (2018) dealt with the case of negative defensive medicine. We show that, when the models admit a nonisochronous center, it is possible to prove the existence of chaotic dynamics for the Poincar\'e map associated with those systems via the method of Linked Twist Maps (LTMs). To such aim we exploit in both frameworks, using a similar rationale, the periodic dependence on time of a model parameter influencing the position of the center and describing some risk associated with certain medical interventions, whose seasonal variation is empirically grounded. We also provide the ecological interpretation of the same model, proposed by Harvie et al. (2007) and connected with intraspecific competition and environmental carrying capacity in predator-prey settings. In this case, it is sensible to assume a seasonal variation both for the carrying capacities and for the intrinsic growth rates of the two populations. Although such parameters influence the shape of the orbits, but do not affect the center position, we show that it is still possible to prove the existence of chaotic dynamics for the associated Poincar\'e map via the LTMs technique dealing with a different geometrical configuration for orbits in the phase plane.


title: Orientation Matters: 6-DoF Autonomous Camera Movement for Minimally Invasive Surgery

abstract: We propose a new method for six-degree-of-freedom (6-DoF) autonomous camera movement for minimally invasive surgery, which, unlike previous methods, takes into account both the position and orientation information from structures in the surgical scene. In addition to locating the camera for a good view of the manipulated object, our autonomous camera takes into account workspace constraints, including the horizon and safety constraints. We developed a simulation environment to test our method on the "wire chaser" surgical training task from validated training curricula in conventional laparoscopy and robot-assisted surgery. Furthermore, we propose, for the first time, the application of the proposed autonomous camera method in video-based surgical skill assessment, an area where videos are typically recorded using fixed cameras. In a study with N=30 human subjects, we show that video examination of the autonomous camera view as it tracks the ring motion over the wire leads to more accurate user error (ring touching the wire) detection than when using a fixed camera view, or camera movement with a fixed orientation. Our preliminary work suggests that there are potential benefits to autonomous camera positioning informed by scene orientation, and this can direct designers of automated endoscopes and surgical robotic systems, especially when using chip-on-tip cameras that can be wristed for 6-DoF motion.


title: Maximizing the Influence: Analytic Activation Probability Computation Approach

abstract: Influence maximization is the problem of finding a subset of the most influential individuals in a network. The impact of social networks on the dissemination of information and the development of viral marketing has made this problem as the subject of many studies. Influence maximization is $\NP$-hard and many greedy algorithms have been proposed to solve the problem. In this paper we propose a greedy algorithm that approximates the influence by using a novel analytic activation probability computation method. We propose a nonlinear system of equations to compute the activation probabilities which is more accurate than state-of-the-arts. Moreover, we propose a method to compute the activation probabilities efficiently, in order to reduce the running time of the algorithm. We examine our proposed methods on some real-world data sets. Our experiments demonstrate that the proposed algorithms outperform the previous works by computing more accurate activation probabilities. In addition, Our efficient method has much improvements on the running-time, so it can be easily scaled on large networks.


title: MinMax Mean-Field Team Approach for a Leader-Follower Network: A Saddle-Point Strategy

abstract: This paper investigates a soft-constrained MinMax control problem of a leader-follower network. The network consists of one leader and an arbitrary number of followers that wish to reach consensus with minimum energy consumption in the presence of external disturbances. The leader and followers are coupled in the dynamics and cost function. Two non-classical information structures are considered: mean-field sharing and intermittent mean-field sharing, where the mean-field refers to the aggregate state of the followers. In mean-field sharing, every follower observes its local state, the state of the leader and the mean field while in the intermittent mean-field sharing, the mean-field is only observed at some (possibly no) time instants. A social welfare cost function is defined, and it is shown that a unique saddle-point strategy exists which minimizes the worst-case value of the cost function under mean-field sharing information structure. The solution is obtained by two scalable Riccati equations, which depend on a prescribed attenuation parameter, serving as a robustness factor. For the intermittent mean-field sharing information structure, an approximate saddle-point strategy is proposed, and its converges to the saddle-point is analyzed. Two numerical examples are provided to demonstrate the efficacy of the obtained results.


title: The DeRisk database: Extreme Design Waves for Offshore Wind Turbines

abstract: The estimation of extreme loads from waves is an essential part of the design of an offshore wind turbine. Standard design codes suggest to either use simplified methods based on regular waves, or to perform fully nonlinear computations. The former might not provide an accurate representation of the extreme waves, while the latter is computationally too intensive for design iterations. We address these limitations by using the fully nonlinear solver OceanWave3D to establish the DeRisk database, a large dataset of extreme waves kinematics in a two-dimensional domain. From the database, which is open and freely available, a designer can extract fully-nonlinear wave kinematics for a wave condition and water depth of interest by identifying a suitable computation in the database and, if needed, by Froude-scaling the kinematics. The nonlinear solver is validated against the DeRisk model experiments at two different water depths, $33.0 [m]$ and $20.0 [m]$, and an excellent agreement is found for the analyzed cases. The experiments are used to calibrate OceanWave3D's numerical breaking filter constant, and the best agreement is found for $\beta=0.5$. We compare the experimental static force with predictions by the DeRisk database and the Rainey force model, and with state-of-the-art industrial practices. For milder storms, we find a good agreement in the predicted extreme force between the present methodology and the standard methodologies. At the deep location and for stronger storms, the largest loads are given by slamming loads due to breaking waves. In this condition, the database methodology is less accurate than the embedded stream function method and more accurate than the WiFi JIP methodology, providing generally nonconservative estimates. For strong storms at the shallower location, where wave breaking is less dominating, the database methodology is the most accurate overall.


title: Ranking-based Convolutional Neural Network Models for Peptide-MHC Binding Prediction

abstract: T-cell receptors can recognize foreign peptides bound to major histocompatibility complex (MHC) class-I proteins, and thus trigger the adaptive immune response. Therefore, identifying peptides that can bind to MHC class-I molecules plays a vital role in the design of peptide vaccines. Many computational methods, for example, the state-of-the-art allele-specific method MHCflurry, have been developed to predict the binding affinities between peptides and MHC molecules. In this manuscript, we develop two allele-specific Convolutional Neural Network (CNN)-based methods named ConvM and SpConvM to tackle the binding prediction problem. Specifically, we formulate the problem as to optimize the rankings of peptide-MHC bindings via ranking-based learning objectives. Such optimization is more robust and tolerant to the measurement inaccuracy of binding affinities, and therefore enables more accurate prioritization of binding peptides. In addition, we develop a new position encoding method in ConvM and SpConvM to better identify the most important amino acids for the binding events. Our experimental results demonstrate that our models significantly outperform the state-of-the-art methods including MHCflurry with an average percentage improvement of 6.70% on AUC and 17.10% on ROC5 across 128 alleles.


title: Revisiting the Figure of Merit of Concentrated Solar Power Receivers

abstract: The figure of merit (FOM) is a widely used metric to characterize the performance of concentrated solar power (CSP) receivers by comparing the amount of solar thermal energy retained by the receiver to the incident concentrated solar radiation. However, the FOM is a strong function of the concentration factor and receiver temperature, thus direct comparison of FOM values measured under disparate operating conditions is inappropriate. To remedy this problem, the present study proposes a new metric called the receiver effectiveness calculated by normalizing the actual FOM with its theoretical maximum. The receiver effectiveness can be employed for comparing receiver performances regardless of their operating conditions, and can be treated as more-like the second law efficiency of thermodynamics. In addition, a theoretical limit of the CSP plant efficiency is also examined by combining the maximum FOM and the Carnot efficiency for different concentration factors and receiver temperatures. The calculated maximum CSP plant efficiency clearly indicate that optimizing FOM does not always lead to a better CSP plant performance. Along with the FOM, the proposed receiver effectiveness and maximum CSP system efficiency should be considered as complementary metrics to evaluate the performance of the CSP system.


title: Long-range nontopological edge currents in charge-neutral graphene

abstract: Van der Waals heterostructures display a rich variety of unique electronic properties. To identify novel transport mechanisms, nonlocal measurements have been widely used, wherein a voltage is measured at contacts placed far away from the expected classical flow of charge carriers. This approach was employed in search of dissipationless spin and valley transport, topological charge-neutral currents, hydrodynamic flows and helical edge modes. Monolayer, bilayer, and few-layer graphene, transition-metal dichalcogenides, and moire superlattices were found to display pronounced nonlocal effects. However, the origin of these effects is hotly debated. Graphene, in particular, exhibits giant nonlocality at charge neutrality, a prominent behavior that attracted competing explanations. Utilizing superconducting quantum interference device on a tip (SQUID-on-tip) for nanoscale thermal and scanning gate imaging, we demonstrate that the commonly-occurring charge accumulation at graphene edges leads to giant nonlocality, producing narrow conductive channels that support long-range currents. Unexpectedly, while the edge conductance has little impact on the current flow in zero magnetic field, it leads to field-induced decoupling between edge and bulk transport at moderate fields. The resulting giant nonlocality both at charge neutrality and away from it produces exotic flow patterns in which charges can flow against the global electric field. We have visualized surprisingly intricate patterns of nonlocal currents, which are sensitive to edge disorder. The observed one-dimensional edge transport, being generic and nontopological, is expected to support nonlocal transport in many electronic systems, offering insight into numerous controversies in the literature and linking them to long-range guided electronic states at system edges.


title: Kolmogorov operator with the vector field in Nash class

abstract: We consider divergence-form parabolic equation with measurable uniformly elliptic matrix and the vector field in a large class containing, in particular, the vector fields in $L^p$, $p>d$, as well as some vector fields that are not even in $L_{\rm loc}^{2+\varepsilon}$, $\varepsilon>0$. We establish H\"{o}lder continuity of the bounded soutions, sharp two-sided Gaussian bound on the heat kernel, Harnack inequality.


title: Polynomial-time trace reconstruction in the low deletion rate regime

abstract: In the \emph{trace reconstruction problem}, an unknown source string $x \in \{0,1\}^n$ is transmitted through a probabilistic \emph{deletion channel} which independently deletes each bit with some fixed probability $\delta$ and concatenates the surviving bits, resulting in a \emph{trace} of $x$. The problem is to reconstruct $x$ given access to independent traces. Trace reconstruction of arbitrary (worst-case) strings is a challenging problem, with the current state of the art for poly$(n)$-time algorithms being the 2004 algorithm of Batu et al. \cite{BKKM04}. This algorithm can reconstruct an arbitrary source string $x \in \{0,1\}^n$ in poly$(n)$ time provided that the deletion rate $\delta$ satisfies $\delta \leq n^{-(1/2 + \varepsilon)}$ for some $\varepsilon > 0$. In this work we improve on the result of \cite{BKKM04} by giving a poly$(n)$-time algorithm for trace reconstruction for any deletion rate $\delta \leq n^{-(1/3 + \varepsilon)}$. Our algorithm works by alternating an alignment-based procedure, which we show effectively reconstructs portions of the source string that are not "highly repetitive", with a novel procedure that efficiently determines the length of highly repetitive subwords of the source string.


title: Experimental Evaluation of Algorithm-Assisted Human Decision-Making: Application to Pretrial Public Safety Assessment

abstract: Despite an increasing reliance on fully-automated algorithmic decision-making in our day-to-day lives, human beings still make highly consequential decisions. As frequently seen in business, healthcare, and public policy, recommendations produced by algorithms are provided to human decision-makers to guide their decisions. While there exists a fast-growing literature evaluating the bias and fairness of such algorithmic recommendations, an overlooked question is whether they help humans make better decisions. We develop a statistical methodology for experimentally evaluating the causal impacts of algorithmic recommendations on human decisions. We also show how to examine whether algorithmic recommendations improve the fairness of human decisions and derive the optimal decision rules under various settings. We apply the proposed methodology to preliminary data from the first-ever randomized controlled trial that evaluates the pretrial Public Safety Assessment (PSA) in the criminal justice system. A goal of the PSA is to help judges decide which arrested individuals should be released. On the basis of the preliminary data available, we find that providing the PSA to the judge has little overall impact on the judge's decisions and subsequent arrestee behavior. However, our analysis yields some potentially suggestive evidence that the PSA may help avoid unnecessarily harsh decisions for female arrestees regardless of their risk levels while it encourages the judge to make stricter decisions for male arrestees who are deemed to be risky. In terms of fairness, the PSA appears to increase the gender bias against males while having little effect on any existing racial differences in judges' decision. Finally, we find that the PSA's recommendations might be unnecessarily severe unless the cost of a new crime is sufficiently high.


title: Search for Coherent Elastic Scattering of Solar $^8$B Neutrinos in the XENON1T Dark Matter Experiment

abstract: We report on a search for nuclear recoil signals from solar $^8$B neutrinos elastically scattering off xenon nuclei in XENON1T data, lowering the energy threshold from 2.6 keV to 1.6 keV. We develop a variety of novel techniques to limit the resulting increase in backgrounds near the threshold. No significant $^8$B neutrino-like excess is found in an exposure of 0.6 t $\times$ y. For the first time, we use the non-detection of solar neutrinos to constrain the light yield from 1-2 keV nuclear recoils in liquid xenon, as well as non-standard neutrino-quark interactions. Finally, we improve upon world-leading constraints on dark matter-nucleus interactions for dark matter masses between 3 GeV/c$^2$ and 11 GeV/c$^2$ by as much as an order of magnitude.


title: Network Group Testing

abstract: We consider the problem of identifying infected individuals in a population of size N. We introduce a group testing approach that uses significantly fewer than N tests when infection prevalence is low. The most common approach to group testing, Dorfman testing, groups individuals randomly. However, as communicable diseases spread from individual to individual through underlying social networks, our approach utilizes network information to improve performance. Network grouping, which groups individuals by community, weakly dominates Dorfman testing in terms of the expected number of tests used. Network grouping's outperformance is determined by the strength of community structure in the network. When networks have strong community structure, network grouping achieves the lower bound for two-stage testing procedures. As an empirical example, we consider the scenario of a university testing its population for COVID-19. Using social network data from a Danish university, we demonstrate network grouping requires significantly fewer tests than Dorfman. In contrast to many proposed group testing approaches, network grouping is simple for practitioners to implement. In practice, individuals can be grouped by family unit, social group, or work group.


title: Deterministic Random Number Generator Attack against the Kirchhoff-Law-Johnson-Noise Secure Key Exchange Protocol

abstract: This paper demonstrates the vulnerability of the Kirchhoff-Law-Johnson-Noise (KLJN) secure key exchanger to compromised random number generator(s) even if these random numbers are used solely to generate the noises emulating the Johnson noise of Alice's and Bob's resistors. The attacks shown are deterministic in the sense that Eve's knowledge of Alice's and/or Bob's random numbers is basically deterministic. Moreover, no statistical evaluation is needed, except for rarely occurring events of negligible, random waiting time and verification time. We explore two situations. In the first case, Eve knows both Alice's and Bob's random noises. We show that, in this situation, Eve can quickly crack the secure key bit by using Ohm's Law. In the other situation, Eve knows only Bob's random noise. Then Eve first can learn Bob's resistance value by using Ohm's Law. Therefore, she will have the same knowledge as Bob, thus at the end of the bit exchange period, she will know Alice's bit.


title: Velocity of viscous fingers in miscible displacement

abstract: The paper investigates the linear growth of the mixing zone during polymer slug injection into a water reservoir. The velocities of the slug front and of the boundaries of the mixing zone are analyzed as key parameters. Using two different numerical methods (finite volumes and finite elements), the impact of the slug size, reservoir dimensions, Peclet number, and viscosity curve shape on the corresponding velocities is examined. Notwithstanding the realization of the solution by two computational schemes, the simulation results coincide with sufficient accuracy. The numerically obtained velocities are compared with theoretical estimates within the transverse flow equilibrium approximation and Koval model. Based on the comparison pattern, recommendations are presented on the use of specific analytical methods for estimating the growth rate of the mixing zone depending on the characteristics of the polymer.


title: The ThreeHundred: the structure and properties of cosmic filaments in the outskirts of galaxy clusters

abstract: Galaxy cluster outskirts are described by complex velocity fields induced by diffuse material collapsing towards filaments, gas and galaxies falling into clusters, and gas shock processes triggered by substructures. A simple scenario that describes the large-scale tidal fields of the cosmic web is not able to fully account for this variety, nor for the differences between gas and collisionless dark matter. We have studied the filamentary structure in zoom-in resimulations centred on 324 clusters from The ThreeHundred project, focusing on differences between dark and baryonic matter. This paper describes the properties of filaments around clusters out to five $R_{200}$, based on the diffuse filament medium where haloes had been removed. For this, we stack the remaining particles of all simulated volumes to calculate the average profiles of dark matter and gas filaments. We find that filaments increase their thickness closer to nodes and detect signatures of gas turbulence at a distance of $\sim 2 h^{-1}\rm{Mpc}$ from the cluster. These are absent in dark matter. Both gas and dark matter collapse towards filament spines at a rate of $\sim 200 h^{-1} \rm{km ~ s^{-1}} $. We see that gas preferentially enters the cluster as part of filaments, and leaves the cluster centre outside filaments. We further see evidence for an accretion shock just outside the cluster. For dark matter, this preference is less obvious. We argue that this difference is related to the turbulent environment. This indicates that filaments act as highways to fuel the inner regions of clusters with gas and galaxies.


title: Some New Results Concerning Power Graphs and Enhanced Power Graphs of Groups

abstract: The directed power graph $\vec{\mathcal P}(\mathbf G)$ of a group $\mathbf G$ is the simple digraph with vertex set $G$ such that $x\rightarrow y$ if $y$ is a power of $x$. The power graph of $\mathbf G$, denoted by $\mathcal P(\mathbf G)$, is the underlying simple graph. The enhanced power graph $\mathcal P_e(\mathbf G)$ of $\mathbf G$ is the simple graph with vertex set $G$ in which two elements are adjacent if they generate a cyclic subgroup. In this paper, it is proven that, if two groups have isomorphic power graphs, then they have isomorphic enhanced power graphs, too. It is known that any finite nilpotent group of order divisible by at most two primes has perfect enhanced power graph. We investigated whether the same holds for all finite groups, and we have obtained a negative answer to that question. Further, we proved that, for any $n\geq 0$ and prime numbers $p$ and $q$, every group of order $p^nq$ and $p^2q^2$ has perfect enhanced power graph. We also give a complete characterization of symmetric and alternative groups with perfect enhanced graphs.


title: Supersonic turbulent flows over sinusoidal rough walls

abstract: Direct numerical simulations were performed to characterize fully developed supersonic turbulent channel flows over isothermal rough walls. The effect of roughness was incorporated using a level-set/volume-of-fluid immersed boundary method. Turbulence statistics of five channel flows are compared, including one reference case with both walls smooth and four cases with smooth top walls and bottom walls with two-dimensional (2D) and three-dimensional (3D) sinusoidal roughnesses. Results reveal a strong dependence of the turbulence on the roughness topography and the associated shock patterns. Specifically, the 2D geometries generate strong oblique shock waves that propagate across the channel and are reflected back to the rough-wall side. These strong shocks are absent in the smooth-wall channel and are significantly weaker in cases with 3D roughness geometries, replaced by weak shocklets. At the impingement locations of the shocks on the top wall in the 2D roughness cases, localized augmentations of turbulence shear production are observed. Such regions of augmented production also exist for the 3D cases, at a much weaker level. The oblique shock waves are thought to be responsible for a more significant entropy generation for cases with 2D surfaces than those with 3D ones, leading to a higher irreversible heat generation and consequently higher temperature values in 2D roughness cases. In the present supersonic channels, the effects of roughness extend beyond the near-wall layer due to the shocks. This suggests that outer layer similarity may not fully apply to a rough-wall supersonic turbulent flow.


title: A categorical K\"unneth formula for constructible Weil sheaves

abstract: We prove a K\"unneth-type equivalence of derived categories of lisse and constructible Weil sheaves on schemes in characteristic $p > 0$ for various coefficients, including finite discrete rings, algebraic field extensions $E \supset \mathbf Q_\ell$, $\ell \ne p$ and their rings of integers $O_E$. We also consider a variant for ind-construtible sheaves which applies to the cohomology of moduli stacks of shtukas over global function fields.


title: Physics Guided Machine Learning Methods for Hydrology

abstract: Streamflow prediction is one of the key challenges in the field of hydrology due to the complex interplay between multiple non-linear physical mechanisms behind streamflow generation. While physically-based models are rooted in rich understanding of the physical processes, a significant performance gap still remains which can be potentially addressed by leveraging the recent advances in machine learning. The goal of this work is to incorporate our understanding of physical processes and constraints in hydrology into machine learning algorithms, and thus bridge the performance gap while reducing the need for large amounts of data compared to traditional data-driven approaches. In particular, we propose an LSTM based deep learning architecture that is coupled with SWAT (Soil and Water Assessment Tool), an hydrology model that is in wide use today. The key idea of the approach is to model auxiliary intermediate processes that connect weather drivers to streamflow, rather than directly mapping runoff from weather variables which is what a deep learning architecture without physical insight will do. The efficacy of the approach is being analyzed on several small catchments located in the South Branch of the Root River Watershed in southeast Minnesota. Apart from observation data on runoff, the approach also leverages a 200-year synthetic dataset generated by SWAT to improve the performance while reducing convergence time. In the early phases of this study, simpler versions of the physics guided deep learning architectures are being used to achieve a system understanding of the coupling of physics and machine learning. As more complexity is introduced into the present implementation, the framework will be able to generalize to more sophisticated cases where spatial heterogeneity is present.


title: Appearance of objectivity for NV centers interacting with dynamically polarized nuclear environment

abstract: Quantum-to-classical transition still eludes a full understanding. Out of its multiple aspects, one has recently gained an increased attention - the appearance of objective world out of the quantum. One particularly idea is that objectivity appears thanks to specific quantum state structures formation during the evolution, known as Spectrum Broadcast Structures (SBS). Despite that quite some research was already performed on this strongest and most fundamental form of objectivity, its practical realization in a concrete physical medium has not been analyzed so far. In this work, we study the possibility to simulate objectivization process via SBS formation using widely studied Nitrogen-Vacancy centers in diamonds. Assuming achievable limits of dynamical polarization technique, we show that for high, but experimentally viable polarizations ($p>0.5$) of nuclear spins and for magnetic fields lower than $\approx \! 20$ Gauss the state of the NV center and its nearest polarized environment approaches reasonably well an SBS state.


title: On the Resource Allocation for Political Campaigns

abstract: In an election campaign, candidates must decide how to optimally allocate their efforts/resources optimally among the regions of a country. As a result, the outcome of the election will depend on the players' strategies and the voters' preferences. In this work, we present a zero-sum game where two candidates decide how to invest a fixed resource in a set of regions, while considering their sizes and biases. We explore the Majority System (MS) as well as the Electoral College (EC) voting systems. We prove equilibrium existence and uniqueness under MS in a deterministic model; in addition, their closed form expressions are provided when fixing the subset of regions and relaxing the non-negative investing constraint. For the stochastic case, we use Monte Carlo simulations to compute the players' payoffs, together with its gradient and hessian. For the EC, given the lack of Equilibrium in pure strategies, we propose an iterative algorithm to find Equilibrium in mixed strategies in a subset of the simplex lattice. We illustrate numerical instances under both election systems, and contrast players' equilibrium strategies. Finally, we show that polarization induces candidates to focus on larger regions with negative biases under MS, whereas candidates concentrate on swing states under EC.


title: Asymptotic behaviour of ancestral lineages in subcritical continuous-state branching populations

abstract: Consider the population model with infinite size associated to subcritical continuous-state branching processes (CSBP). Individuals reproduce independently according to the same subcritical offspring distribution. We study the long-term behaviour of the ancestral lineages as time goes to the past and show that the flow of ancestral lineages, properly renormalized, converges almost surely to the inverse of a drift-free subordinator whose Laplace exponent is explicit in terms of the branching mechanism. We provide an interpretation in terms of the genealogy of the population. In particular, we show that the inverse subordinator is partitioning the current population into ancestral families with distinct common ancestors. When Grey's condition is satisfied, the population comes from a discrete set of ancestors and the ancestral families are i.i.d and distributed according to the quasi-stationary distribution of the CSBP conditioned on non-extinction. When Grey's condition is not satisfied, the population comes from a continuum of ancestors which is described as the set of increase points $\mathscr{S}$ of the limiting inverse subordinator. The Hausdorff dimension of $\mathscr{S}$ is given. The proof is based on a general result for stochastically monotone processes of independent interest, which relates $\theta$-invariant measures and $\theta$-invariant functions for a process and its Siegmund dual.


title: Strong Near-Infrared Carbon Absorption in the Transitional Type Ia SN 2015bp

abstract: Unburned carbon is potentially a powerful probe of Type Ia supernova (SN) explosion mechanisms. We present comprehensive optical and near-infrared (NIR) data on the "transitional" Type Ia SN 2015bp. An early NIR spectrum ($t = -$9.9 days with respect to B-band maximum) displays a striking C I $\lambda1.0693\,\mu \rm{m}$ line at $11.9 \times 10^3$~km s$^{-1}$, distinct from the prominent Mg II $\lambda1.0927\,\mu \rm{m}$ feature, which weakens toward maximum light. SN 2015bp also displays a clear C II $\lambda6580$A notch early ($t = -10.9$ days) at $13.2 \times 10^3$~km s$^{-1}$, consistent with our NIR carbon detection. At $M_B = -$18.46, SN 2015bp is less luminous than a normal SN Ia and, along with iPTF13ebh, is the second member of the transitional subclass to display prominent early-time NIR carbon absorption. We find it unlikely that the C I feature is misidentified He I $\lambda1.0830\,\mu\rm{m}$ because this feature grows weaker toward maximum light, while the helium line produced in some double-detonation models grows stronger at these times. Intrigued by these strong NIR carbon detections, but lacking NIR data for other SNe Ia, we investigated the incidence of optical carbon in the sample of nine transitional SNe Ia with early-time data ($t \lesssim-$4 days). We find that four display C II $\lambda$6580A, while two others show tentative detections, in line with the SN Ia population as a whole. We conclude that at least $\sim$50% of transitional SNe Ia in our sample do not come from sub-Chandrasekhar mass explosions due to the clear presence of carbon in their NIR and optical spectra.