title: The Cartesian Grid Active Flux Method with Adaptive Mesh Refinement

abstract: We present the first implementation of the Active Flux method on adaptively refined Cartesian grids. The Active Flux method is a third order accurate finite volume method for hyperbolic conservation laws, which is based on the use of point values as well as cell average values of the conserved quantities. The resulting method has a compact stencil in space and time and good stability properties. The method is implemented as a new solver in ForestClaw, a software for parallel adaptive mesh refinement of patch-based solvers. On each Cartesian grid patch the single grid Active Flux method can be applied. The exchange of data between grid patches is organised via ghost cells. The local stencil in space and time and the availability of the point values that are used for the reconstruction, leads to an efficient implementation. The resulting method is third order accurate, conservative and allows the use of subcycling in time.


title: Probability of generation by random permutations of given cycle type

abstract: Suppose $\pi$ and $\pi'$ are two random elements of $S_n$ with constrained cycle types such that $\pi$ has $x n^{1/2}$ fixed points and $yn/2$ two-cycles, and likewise $\pi'$ has $x' n^{1/2}$ fixed points and $y'n/2$ two-cycles. We show that the events that $G = \langle \pi, \pi' \rangle$ is transitive and $G \geq A_n$ both have probability approximately \[(1 - yy')^{1/2} \exp\left(- \frac{xx' + \frac12 x^2 y' + \frac12 {x'}^2 y}{1 - yy'}\right),\] provided $(x, x')$ is not close to $(0, \infty)$ or $(\infty, 0)$. This formula is derived from some preliminary results in a recent paper (arXiv:1904.12180) of the authors. As an application, we show that two uniformly random elements of uniformly random conjugacy classes of $S_n$ generate the group with probability about 51%.


title: Constraining super-light sterile neutrinos at Borexino and KamLAND

abstract: The presence of a super-light sterile neutrino can lead to a dip in the survival probability of solar neutrinos, and explain the suppression of the upturn in the low energy solar neutrino data. In this work, we systematically study the survival probabilities in the 3+1 framework by taking into account of the non-adiabatic transitions and the coherence effect. We obtain an analytic equation that can predict the position of the dip. We also place constraints on the parameter space of sterile neutrinos by using the latest Borexino and KamLAND data. We find that the low and high energy neutrino data at Borexino are sensitive to different regions in the sterile neutrino parameter space. In the case with only $\theta_{01}$ being nonzero, the $\rm{{}^{8}B}$ data sets the strongest bounds at $\Delta m_{01}^{2} \approx (1.1\sim2.2)\Delta m_{21}^{2}$, while the low energy neutrino data is more sensitive to other mass-squared regions. The lowest bounds on $\Delta m_{01}^{2}$ from the $\rm{pp}$ data can reach $10^{-12} \ \rm{eV^{2}}$ because of the coherence effect. Also, due to the presence of non-adiabatic transitions, the bounds in the range of $10^{-9} \ \textrm{eV}^{2} \lesssim \Delta m_{01}^{2} \lesssim 10^{-5} \ \textrm{eV}^{2}$ become weaker as $\Delta m_{01}^{2}$ or $\sin^{2}2\theta_{01}$ decreases. We also find that in the case with only $\theta_{02}$ or $\theta_{03}$ being nonzero, the low energy solar neutrino data set similar but weaker bounds as compared to the case with only $\theta_{01}$ being nonzero. However, the bounds from the high energy solar data and the KamLAND data are largely affected by the sterile mixing angles.


title: An automatic pipeline for atlas-based fetal and neonatal brain segmentation and analysis

abstract: The automatic segmentation of perinatal brain structures in magnetic resonance imaging (MRI) is of utmost importance for the study of brain growth and related complications. While different methods exist for adult and pediatric MRI data, there is a lack for automatic tools for the analysis of perinatal imaging. In this work, a new pipeline for fetal and neonatal segmentation has been developed. We also report the creation of two new fetal atlases, and their use within the pipeline for atlas-based segmentation, based on novel registration methods. The pipeline is also able to extract cortical and pial surfaces and compute features, such as curvature, thickness, sulcal depth, and local gyrification index. Results show that the introduction of the new templates together with our segmentation strategy leads to accurate results when compared to expert annotations, as well as better performances when compared to a reference pipeline (developing Human Connectome Project (dHCP)), for both early and late-onset fetal brains.


title: Probabilistic multivariate early warning signals

abstract: A broad range of natural and social systems from human microbiome to financial markets can go through critical transitions, where the system suddenly collapses to another stable configuration. Critical transitions can be unexpected, with potentially catastrophic consequences. Anticipating them early and accurately can facilitate controlled system manipulation and mitigation of undesired outcomes. Obtaining reliable predictions have been difficult, however, as often only a small fraction of the relevant variables can be monitored, and even minor perturbations can induce drastic changes in fragile states of a complex system. Data-driven indicators have been proposed as an alternative to prediction and signal an increasing risk of forthcoming transitions. Autocorrelation and variance are examples of generic indicators that tend to increase at the vicinity of an approaching tipping point across a range of systems. An important shortcoming in these and other widely studied indicators is that they deal with simplified one-dimensional representations of complex systems. Here, we demonstrate that a probabilistic data aggregation strategy can provide new ways to improve early warning detection by more efficiently utilizing the available information from multivariate time series. In particular, we consider a probabilistic variant of a vector autoregression model as a novel early warning indicator and argue that it has theoretical advantages related to model regularization, treatment of uncertainties, and parameter interpretation. We evaluate the performance against alternatives in a simulation benchmark and show improved sensitivity in EWS detection in a common ecological model encompassing multiple interacting species.


title: Rotating Love: The dynamical tides of spinning Newtonian stars

abstract: We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilising a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order.


title: A natural mechanism for a SM-like Higgs boson in the 2HDM without decoupling

abstract: The properties of the Higgs boson discovered at the Large Hadron Collider are very well described by the Standard Model (SM). Thus, any theory that invokes an extended Higgs sector must explain why the neutral scalar observed at the LHC so closely resembles the SM Higgs boson. In this talk, I review the Higgs alignment limit, in which one neutral scalar state of the Higgs sector is SM-like. An approximate Higgs alignment can be achieved "naturally" either via decoupling or via an approximate symmetry. Using the two-Higgs doublet model as a prototype for an extended Higgs sector, I examine the symmetries of the scalar potential and their soft breakings that may be responsible for the SM-like properties of the observed Higgs boson, and I demonstrate how to extend such (softly-broken) symmetries to the Yukawa sector of the model.


title: Is climate change time reversible?

abstract: This paper proposes strategies to detect time reversibility in stationary stochastic processes by using the properties of mixed causal and noncausal models. It shows that they can also be used for non-stationary processes when the trend component is computed with the Hodrick-Prescott filter rendering a time-reversible closed-form solution. This paper also links the concept of an environmental tipping point to the statistical property of time irreversibility and assesses fourteen climate indicators. We find evidence of time irreversibility in $GHG$ emissions, global temperature, global sea levels, sea ice area, and some natural oscillation indices. While not conclusive, our findings urge the implementation of correction policies to avoid the worst consequences of climate change and not miss the opportunity window, which might still be available, despite closing quickly.


title: Low-Resolution Imaging FMCW Lidar

abstract: We demonstrate the imaging capability of a frequency modulated continuous wave lidar based on a fiber bundle. The lidar constructs velocity and range images for hard targets at a rate of 60 Hz. The sensing range is up to 30 m with 20 mW of output power. The instrument employs custom electronics with seven parallel heterodyne receivers. An example of image recovery is presented on 6-pixel "pictures" of a spinning disk and a drone hovering in the air. In experiments, we also tested the laser tuning linearity correction with a phase-locked loop. We see the practicality of such a low-resolution system as a boost in scanning rate of conventional lidars or for direct target imaging with a further upgrade of pixel count.


title: On the $f$-vectors of $r$-multichain subdivisions

abstract: For a poset $P$ and an integer $r\geq 1$, let $P_r$ be a collection of all $r$-multichains in $P$. Corresponding to each strictly increasing map $\i:[r]\rightarrow [2r]$, there is an order $\preceq_{\i}$ on $P_r$. Let $\D(G_{\i}(P_r))$ be the clique complex of the graph $G_{\i}$ associated to $P_r$ and $\i$. In a recent paper \cite{NW}, it is shown that $\D(G_{\i}(P_r))$ is a subdivision of $P$ for a class of strictly increasing maps. In this paper, we show that all these subdivisions have the same $f$-vector. We give an explicit description of the transformation matrices from the $f$- and $h$-vectors of $\Delta$ to the $f$- and $h$-vectors of these subdivisions when $P$ is a poset of faces of $\D$. We study two important subdivisions Cheeger-M\"{u}ller-Schrader's subdivision and the $r$-colored barycentric subdivision which fall in our class of $r$-multichain subdivisions.


title: Chemical transformer compression for accelerating both training and inference of molecular modeling

abstract: Transformer models have been developed in molecular science with excellent performance in applications including quantitative structure-activity relationship (QSAR) and virtual screening (VS). Compared with other types of models, however, they are large, which results in a high hardware requirement to abridge time for both training and inference processes. In this work, cross-layer parameter sharing (CLPS), and knowledge distillation (KD) are used to reduce the sizes of transformers in molecular science. Both methods not only have competitive QSAR predictive performance as compared to the original BERT model, but also are more parameter efficient. Furthermore, by integrating CLPS and KD into a two-state chemical network, we introduce a new deep lite chemical transformer model, DeLiCaTe. DeLiCaTe captures general-domains as well as task-specific knowledge, which lead to a 4x faster rate of both training and inference due to a 10- and 3-times reduction of the number of parameters and layers, respectively. Meanwhile, it achieves comparable performance in QSAR and VS modeling. Moreover, we anticipate that the model compression strategy provides a pathway to the creation of effective generative transformer models for organic drug and material design.


title: A least-squares Galerkin approach to gradient recovery for Hamilton-Jacobi-Bellman equation with Cordes coefficients

abstract: We propose a conforming finite element method to approximate the strong solution of the second order Hamilton-Jacobi-Bellman equation with Dirichlet boundary and coefficients satisfying Cordes condition. We show the convergence of the continuum semismooth Newton method for the fully nonlinear Hamilton-Jacobi-Bellman equation. Applying this linearization for the equation yields a recursive sequence of linear elliptic boundary value problems in nondivergence form. We deal numerically with such BVPs via the least-squares gradient recovery of Lakkis & Mousavi [2021, arxiv:1909.00491]. We provide an optimal-rate apriori and aposteriori error bounds for the approximation. The aposteriori error are used to drive an adaptive refinement procedure. We close with computer experiments on uniform and adaptive meshes to reconcile the theoretical findings.


title: GraphSPME: Markov Precision Matrix Estimation and Asymptotic Stein-Type Shrinkage

abstract: GraphSPME is an open source Python, R and C++ header-only package implement-ing non-parametric sparse precision matrix estimation along with asymptotic Stein-type shrinkage estimation of the covariance matrix. The user defines a potential neighbourhood structure and provides data that potentially are p >> n. This paper introduces a novel approach for finding the optimal order (that data allows to estimate) of a potential Markov property. The algorithm is implemented in the package, alleviating the problem of users making Markov assumptions and implementing corresponding complex higher-order neighbourhood structures. Estimation is made accurate and stable by simultaneously utilising both Markov properties and Stein-type shrinkage. Asymptotic results on Stein-type shrinkage ensure that non-singular well conditioned matrices are obtained in an automatic manner. Final symmetry conversion creates symmetric positive definite estimates. Furthermore, the estimation routine is made efficient and scalable to very high-dimensional problems (~10^7) by utilising the sparse nature of the precision matrix under Markov assumptions. Implementation wise, the sparsity is exploited by employing the sparsity possibilities made available by the Eigen C++ linear-algebra library. The package and examples are available at https://github.com/equinor/GraphSPME


title: Continuum limit of parton distribution functions from the pseudo-distribution approach on the lattice

abstract: Precise quantification of the structure of nucleons is one of the crucial aims of hadronic physics for the coming years. The expected progress related to ongoing and planned experiments should be accompanied by calculations of partonic distributions from lattice QCD. While key insights from the lattice are expected to come for distributions that are difficult to access experimentally, it is important that lattice QCD can reproduce the well-known unpolarized parton distribution functions (PDFs) with full control over systematic uncertainties. One of the novel methods for accessing the partonic $x$-dependence is the pseudo-distribution approach, which employs matrix elements of a spatially-extended nonlocal Wilson-line operator of length $z$. In this paper, we address the issue of discretization effects, related to the necessarily nonzero value of the lattice spacing $a$, which start at first order in $a$ as a result of the nonlocal operator. We use twisted mass fermions simulated at three values of the lattice spacing, at a pion mass of 370 MeV, and extract the continuum limit of isovector unpolarized PDFs. We also test, for the first time in the pseudo-distribution approach, the effects of the recently derived two-loop matching. Finally, we address the issue of the reliability of the extraction with respect to the maximal value of $z$.


title: Systematical study of $\Omega_c$-like molecular states from interactions $\Xi_c^{(',*)}\bar{K}^{(*)}$ and $\Xi^{(*)}D^{(*)}$

abstract: In this work, the $\Omega_c$-like molecular states are systematically investigated in a quasipotential Bethe-Salpeter equation approach. The relevant interactions $\Xi_c^{(*,')}\bar{K}^{(*)}$, $\Xi^{(*)}D^{(*)}$, and $\Omega^{(*)}_c(\pi/\eta/\rho/\omega)$ are described by light meson exchanges with the help of the effective Lagrangians with SU(3), chiral, and heavy quark symmetries. The obtained potential kernels of considered interactions are inserted into the quasipotential Bethe-Salpeter equation, and coupled-channel calculations are performed to find possible molecular states and its couplings to the channels considered. The results suggest that an isoscalar state can be produced from the $\Xi^*_c\bar{K}$ interaction with spin parity $3/2^-$, which can be related to state $\Omega_c(3120)$. And its isoscalar partner is predicted with a dominant decay in the $\Omega_c^*\pi$ channel. The isoscalar and isovector states with $1/2^-$ can be produced from the $\Xi'_c\bar{K}$ interaction with a threshold close to the mass of the $\Omega_c(3050)$ and $\Omega_c(3065)$. Their couplings to the $\Xi_c \bar{K}$ channel are very weak, and the isovector one has strong coupling to $\Omega_c\pi$. High precision measurement is helpful to confirm or search such molecular states. Experimental search of states with higher masses generated from interactions $\Xi^{(*,')}_c \bar{K}^*$ and $\Xi^* D^{(*)}$ are also suggested by the current results.


title: Violations of Quark-Hadron Duality in Low-Energy Determinations of $\alpha_s$

abstract: Using the spectral functions measured in $\tau$ decays, we investigate the actual numerical impact of duality violations on the extraction of the strong coupling. These effects are tiny in the standard $\alpha_s(m_\tau^2)$ determinations from integrated distributions of the hadronic spectrum with pinched weights, or from the total $\tau$ hadronic width. The pinched-weight factors suppress very efficiently the violations of duality, making their numerical effects negligible in comparison with the larger perturbative uncertainties. However, combined fits of $\alpha_s$ and duality-violation parameters, performed with non-protected weights, are subject to large systematic errors associated with the assumed modelling of duality-violation effects. These uncertainties have not been taken into account in the published analyses, based on specific models of quark-hadron duality.


title: Characterization of the Gray-Wyner Rate Region for Multivariate Gaussian Sources: Optimality of Gaussian Auxiliary RV

abstract: Examined in this paper, is the Gray and Wyner achievable lossy rate region for a tuple of correlated multivariate Gaussian random variables (RVs) $X_1 : \Omega \rightarrow {\mathbb R}^{p_1}$ and $X_2 : \Omega \rightarrow {\mathbb R}^{p_2}$ with respect to square-error distortions at the two decoders. It is shown that among all joint distributions induced by a triple of RVs $(X_1,X_2, W)$, such that $W : \Omega \rightarrow {\mathbb W} $ is the auxiliary RV taking continuous, countable, or finite values, the Gray and Wyner achievable rate region is characterized by jointly Gaussian RVs $(X_1,X_2, W)$ such that $W $ is an $n$-dimensional Gaussian RV. It then follows that the achievable rate region is parametrized by the three conditional covariances $Q_{X_1,X_2|W}, Q_{X_1|W}, Q_{X_2|W}$ of the jointly Gaussian RVs. Furthermore, if the RV $W$ makes $X_1$ and $X_2$ conditionally independent, then the corresponding subset of the achievable rate region, is simpler, and parametrized by only the two conditional covariances $Q_{X_1|W}, Q_{X_2|W}$. The paper also includes the characterization of the Pangloss plane of the Gray-Wyner rate region along with the characterizations of the corresponding rate distortion functions, their test-channel distributions, and structural properties of the realizations which induce these distributions.


title: Fundamental Laws of Binary Classification

abstract: Finding discriminant functions of minimum risk binary classification systems is a novel geometric locus problem -- which requires solving a system of fundamental locus equations of binary classification -- subject to deep-seated statistical laws. We show that a discriminant function of a minimum risk binary classification system is the solution of a locus equation that represents the geometric locus of the decision boundary of the system, wherein the discriminant function is connected to the decision boundary by an exclusive principal eigen-coordinate system -- at which point the discriminant function is represented by a geometric locus of a novel principal eigenaxis -- structured as a dual locus of likelihood components and principal eigenaxis components. We demonstrate that a minimum risk binary classification system acts to jointly minimize its eigenenergy and risk by locating a point of equilibrium, at which point critical minimum eigenenergies exhibited by the system are symmetrically concentrated in such a manner that the novel principal eigenaxis of the system exhibits symmetrical dimensions and densities, so that counteracting and opposing forces and influences of the system are symmetrically balanced with each other -- about the geometric center of the locus of the novel principal eigenaxis -- whereon the statistical fulcrum of the system is located. Thereby, a minimum risk binary classification system satisfies a state of statistical equilibrium -- so that the total allowed eigenenergy and the expected risk exhibited by the system are jointly minimized within the decision space of the system -- at which point the system exhibits the minimum probability of classification error.


title: Metaverse: Security and Privacy Issues

abstract: The metaverse promises a host of bright opportunities for business, economics, and society. Though, a number of critical aspects are still to be considered and the analysis of their impact is almost non-existent. In this paper, we provide several contributions. We start by analysing the foundations of the metaverse, later we focus on the novel privacy and security issues introduced by this new paradigm, and finally we broaden the scope of the contribution highlighting some of the far-reaching yet logical implications of the metaverse on a number of domains, not all of them in tech. Throughout the paper, we also discuss possible research directions. We believe that the provided holistic view on the foundations, technology, and issues related to the metaverse-with a focus on security and privacy-, other than being an interesting contribution on its own, could also pave the way for a few multidisciplinary research avenues.


title: A new galaxy spectral energy distribution model consistent with the evolution of dust

abstract: The spectral energy distribution (SED) of galaxies provides fundamental information on the related physical processes. However, the SED is significantly affected by dust in its interstellar medium. Dust is mainly produced by asymptotic giant branch stars and Type II supernovae. In addition, the dust mass increases through the metal accretion, and the grain size changes by the collisions between the grains. The contribution of each process and the extinction depend on the size distribution. Therefore, the SED model should treat the evolution of the dust mass and size distribution. In spite of the importance of dust evolution, many previous SED models have not considered the evolution of the total mass and size distribution in a physically consistent manner. In this work, we constructed a new radiative transfer SED model, based on our dust evolution model consistent with the chemical evolution. To reduce the computational cost, we adopted the mega-grain and the one-dimensional plane parallel galaxy approximation. As a fiducial case, we calculated Milky Way-like galaxy SEDs at various ages under the closed-box model. We found that a galaxy at the age of 100~Myr does not produce small grains such as polycyclic aromatic hydrocarbons. After 1~Gyr, we observed a drastic increase of infrared emission and attenuation caused by a rapid increase of dust mass. This phenomenon can be treated appropriately for the first time by our new model. This model can be used for the SED fitting to a galaxy at any stage of evolution.


title: Qualitative Differences Between Evolutionary Strategies and Reinforcement Learning Methods for Control of Autonomous Agents

abstract: In this paper we analyze the qualitative differences between evolutionary strategies and reinforcement learning algorithms by focusing on two popular state-of-the-art algorithms: the OpenAI-ES evolutionary strategy and the Proximal Policy Optimization (PPO) reinforcement learning algorithm -- the most similar methods of the two families. We analyze how the methods differ with respect to: (i) general efficacy, (ii) ability to cope with sparse rewards, (iii) propensity/capacity to discover minimal solutions, (iv) dependency on reward shaping, and (v) ability to cope with variations of the environmental conditions. The analysis of the performance and of the behavioral strategies displayed by the agents trained with the two methods on benchmark problems enable us to demonstrate qualitative differences which were not identified in previous studies, to identify the relative weakness of the two methods, and to propose ways to ameliorate some of those weakness. We show that the characteristics of the reward function has a strong impact which vary qualitatively not only for the OpenAI-ES and the PPO but also for alternative reinforcement learning algorithms, thus demonstrating the importance of optimizing the characteristic of the reward function to the algorithm used.


title: A $(3+\varepsilon)$-Approximate Correlation Clustering Algorithm in Dynamic Streams

abstract: Grouping together similar elements in datasets is a common task in data mining and machine learning. In this paper, we study streaming and parallel algorithms for correlation clustering, where each pair of elements is labeled either similar or dissimilar. The task is to partition the elements and the objective is to minimize disagreements, that is, the number of dissimilar elements grouped together and similar elements that get separated. Our main contribution is a semi-streaming algorithm that achieves a $(3 + \varepsilon)$-approximation to the minimum number of disagreements using a single pass over the stream. In addition, the algorithm also works for dynamic streams. Our approach builds on the analysis of the PIVOT algorithm by Ailon, Charikar, and Newman [JACM'08] that obtains a $3$-approximation in the centralized setting. Our design allows us to sparsify the input graph by ignoring a large portion of the nodes and edges without a large extra cost as compared to the analysis of PIVOT. This sparsification makes our technique applicable in several models of massive graph processing, such as semi-streaming and Massively Parallel Computing (MPC), where sparse graphs can typically be handled much more efficiently. Our work improves on the approximation ratio of the recent single-pass $5$-approximation algorithm and on the number of passes of the recent $O(1/\varepsilon)$-pass $(3 + \varepsilon)$-approximation algorithm [Behnezhad, Charikar, Ma, Tan FOCS'22, SODA'23]. Our algorithm is also more robust and can be applied in dynamic streams. Furthermore, it is the first single pass $(3 + \varepsilon)$-approximation algorithm that uses polynomial post-processing time.


title: Random walks and rank one isometries on CAT(0) spaces

abstract: Let $G$ be a discrete group, $\mu$ a measure on $G$ and $X$ a proper CAT(0) space. We show that if $G$ acts non-elementarily with a rank one element on $X$, then the pushforward $\{Z_n o \}_n$ to $X$ of the random walk generated by $\mu$ converges almost surely to a rank one point of the boundary. We also show that in this context, there is a unique stationary measure on the visual boundary $\partial_\infty X$ of $X$, and that the drift of the random walk is almost surely positive.


title: The young stars in the Galactic Center

abstract: We present a large ${\sim 30" \times 30"}$ spectroscopic survey of the Galactic Center using the SINFONI IFU at the VLT. Combining observations of the last two decades we compile spectra of over $2800$ stars. Using the Bracket-$\gamma$ absorption lines we identify $195$ young stars, extending the list of known young stars by $79$. In order to explore the angular momentum distribution of the young stars, we introduce an isotropic cluster prior. This prior reproduces an isotropic cluster in a mathematically exact way, which we test through numerical simulations. We calculate the posterior angular momentum space as function of projected separation from Sgr~A*. We find that the observed young star distribution is substantially different from an isotropic cluster. We identify the previously reported feature of the clockwise disk and find that its angular momentum changes as function of separation from the black hole, and thus confirm a warp of the clockwise disk ($p \sim 99.2\%$). At large separations, we discover three prominent overdensities of angular momentum. One overdensity has been reported previously, the counter-clockwise disk. The other two are new. Determining the likely members of these structures, we find that as many as $75\%$ of stars can be associated with one of these features. Stars belonging to the warped clockwise-disk show a top heavy K-band luminosity function, while stars belonging to the larger separation features do not. Our observations are in good agreement with the predictions of simulations of in-situ star formation, and argue for common formation of these structures.


title: Exact Exponents for Concentration and Isoperimetry in Product Polish Spaces

abstract: In this paper, we derive variational formulas for the asymptotic exponents (i.e., convergence rates) of the concentration and isoperimetric functions in the product Polish probability space under certain mild assumptions. These formulas are expressed in terms of relative entropies (which are from information theory) and optimal transport cost functionals (which are from optimal transport theory). Hence, our results verify an intimate connection among information theory, optimal transport, and concentration of measure or isoperimetric inequalities. In the concentration regime, the corresponding variational formula is in fact a dimension-free bound in the sense that this bound is valid for any dimension. A cardinality bound for the alphabet of the auxiliary random variable in the expression of the asymptotic isoperimetric exponent is provided, which makes the expression computable by a finite-dimensional program for the finite alphabet case. We lastly apply our results to obtain an isoperimetric inequality in the classic isoperimetric setting, which is asymptotically sharp under certain conditions. The proofs in this paper are based on information-theoretic and optimal transport techniques.


title: Modeling the impact of civilian firearm ownership in the evolution of violent crimes

abstract: We propose a simple mathematical model to describe the evolution of violent crimes. For such purpose, we built a model based on ordinary differential equations that take into account the number of violent crimes and the number of legal and illegal guns. The dynamics is governed by probabilities, modeling for example the police action, the risk perception regarding crimes that leads to increase of ownership of legal guns, and so on. Our analytical and numerical results show that, in addition to the rise of criminality due to the presence of illegal guns, the increase of legal guns leads to a fast increase of violent crimes, suggesting that the access of firearms by civilians is not a good option regarding the control of crimes.


title: Cell-Free MmWave Massive MIMO Systems with Low-Capacity Fronthaul Links and Low-Resolution ADC/DACs

abstract: In this paper, we consider the uplink channel estimation phase and downlink data transmission phase of cell-free millimeter wave (mmWave) massive multiple-input multiple-output (MIMO) systems with low-capacity fronthaul links and low-resolution analog-to-digital converters/digital-to-analog converters (ADC/DACs). In cell-free massive MIMO, a control unit dictates the baseband processing at a geographical scale, while the base stations communicate with the control unit through fronthaul links. Unlike most of previous works in cell-free massive MIMO with finite-capacity fronthaul links, we consider the general case where the fronthaul capacity and ADC/DAC resolution are not necessarily the same. In particular, the fronthaul compression and ADC/DAC quantization occur independently where each one is modeled based on the information theoretic argument and additive quantization noise model (AQNM). Then, we address the codebook design problem that aims to minimize the channel estimation error for the independent and identically distributed (i.i.d.) and colored compression noise cases. Also, we propose an alternating optimization (AO) method to tackle the max-min fairness problem. In essence, the AO method alternates between two subproblems that correspond to the power allocation and codebook design problems. The AO method proposed for the zero-forcing (ZF) precoder is guaranteed to converge, whereas the one for the maximum ratio transmission (MRT) precoder has no such guarantee. Finally, the performance of the proposed schemes is evaluated by the simulation results in terms of both energy and spectral efficiency. The numerical results show that the proposed scheme for the ZF precoder yields spectral and energy efficiency 28% and 15% higher than that of the best baseline.


title: On the operators of Hardy-Littlewood-P\'olya type

abstract: In this paper we introduce and study several new Hardy-Littlewood-P\'olya-type operators. In particular, we study a Hardy-Littlewood-P\'olya-type operator induced by a positive Borel measure on $[0,1)$. We establish some sufficient and necessary conditions for the boundedness (compactness) of these operators. We also determine the exact values of the norms of the Hardy-Littlewood-P\'olya-type operators for certain special cases.


title: An easy technique for focus characterization and optimization of XUV and soft x-ray pulses

abstract: For many applications of extreme ultraviolet (XUV) and x-ray pulses a small focus size is crucial to reach the required intensity or spatial resolution. In this article, we present a simple way to characterize an XUV focus with a resolution of 1.85 micrometer. Furthermore, this technique is applied for the measurement and optimization of the focus of an ellipsoidal mirror for photon energies ranging from 18 to 150 eV generated by high-order harmonics. We envisage a broad range of applications of this approach with sub-micrometer resolution from high-harmonic sources via synchrotrons to free-electron lasers.


title: Volumetric-mapping-based inverse design of 3D architected materials and mobility control by topology reconstruction

abstract: The recent development of modular origami structures has ushered in a new era for active metamaterials with multiple degrees of freedom (multi-DOF). Notably, no systematic inverse design approach for volumetric modular origami structures has been reported. Moreover, very few topologies of modular origami have been studied for the design of active metamaterials with multi-DOF. Herein, we develop an inverse design method and reconfigurable algorithm for constructing 3D active architected structures - we synthesize modular origami structures that can be volumetrically mapped to a target 3D shape. We can control the reconfigurability by reconstructing the topology of the architected structures. Our inverse design based on volumetric mapping with mobility control by topology reconstruction can be used to construct architected metamaterials with any 3D complex shape that are also transformable with multi-DOF. Our work opens a new path toward 3D reconfigurable structures based on volumetric inverse design. This work is significant for the design of 3D active metamaterials and 3D morphing devices for automotive, aerospace, and biomedical engineering applications.


title: Graphene/fluorinated graphene systems for a wide spectrum of electronics application

abstract: Heterostructures prepared from graphene and fluorographene (FG) using the technology of 2D printing on solid and flexible substrates were fabricated and studied. Excellent stability of printed graphene layers and, to a lesser degree, composite graphene: PEDOT: PSS layers were shown. Extraordinary properties of FG as an insulating layer for graphene-based heterostructures at fluorination degree above 30% were demonstrated. It is shown that the leakage current in thin (20-40 nm) films is normally smaller than 10^-8 A/cm2, the breakdown field being greater than 108 V/cm. In hybrid structures with printed FG layers in which graphene was transferred onto, or capsulated with, an FG layer, an increase in charge-carrier mobility and material conductivity amounting to 5-6 times was observed. The spectrum of future applications of FG layers can be further extended due to the possibility of obtaining, from weakly fluorinated graphene (< 20%), functional layers exhibiting a negative differential resistance behavior and, at fluorination degrees of 20-23%, field-effect-transistor channels with current modulation reaching several orders. Composite or bilayer films based on fluorographene and V2O5 or polyvinyl alcohol exhibit a stable resistive switching behavior. On the whole, graphene/FG heterostructures enjoy huge potential for their use in a wide spectrum of application, including flexible electronics.


title: Assessing the Limits of the Distributional Hypothesis in Semantic Spaces: Trait-based Relational Knowledge and the Impact of Co-occurrences

abstract: The increase in performance in NLP due to the prevalence of distributional models and deep learning has brought with it a reciprocal decrease in interpretability. This has spurred a focus on what neural networks learn about natural language with less of a focus on how. Some work has focused on the data used to develop data-driven models, but typically this line of work aims to highlight issues with the data, e.g. highlighting and offsetting harmful biases. This work contributes to the relatively untrodden path of what is required in data for models to capture meaningful representations of natural language. This entails evaluating how well English and Spanish semantic spaces capture a particular type of relational knowledge, namely the traits associated with concepts (e.g. bananas-yellow), and exploring the role of co-occurrences in this context.


title: Spin coherence times of point defects in two-dimensional materials from first principles

abstract: The spin coherence times of 69 triplet defect centers in 45 different 2D host materials are calculated using the cluster correlation expansion (CCE) method with parameters of the spin Hamiltonian obtained from density functional theory (DFT). Several of the triplets are found to exhibit extraordinarily large spin coherence times making them interesting for quantum information processing. The dependence of the spin coherence time on various factors, including the hyperfine coupling strength, the dipole-dipole coupling, and the nuclear g-factors, are systematically investigated. The analysis shows that the spin coherence time is insensitive to the atomistic details of the defect center and rather is dictated by the nuclear spin properties of the host material. Symbolic regression is then used to derive a simple expression for spin coherence time, which is validated on a test set of 55 doublet defects unseen by the regression model. The simple expression permits order-of-magnitude estimates of the spin coherence time without expensive first principles calculations.


title: Optimal flat functions in Carleman-Roumieu ultraholomorphic classes in sectors

abstract: We construct optimal flat functions in Carleman-Roumieu ultraholomorphic classes associated to general strongly nonquasianalytic weight sequences, and defined on sectors of suitably restricted opening. A general procedure is presented in order to obtain linear continuous extension operators, right inverses of the Borel map, for the case of regular weight sequences in the sense of Dyn'kin. Finally, we discuss some examples (including the well-known $q$-Gevrey case) where such optimal flat functions can be obtained in a more explicit way.


title: Evidence of enhanced Zn-diffusion observed during the growth of Inverted Metamorphic Solar Cells

abstract: Zinc-diffusion can induce multiple failures in the electrical performance of a multijunction solar cell. In this work, we show an important Zn-diffusion from the AlGaInP back-surface-field layer to the emitter of the GaInP top cell of an inverted multijunction solar cell. Through the analysis of different doping profiles, we provide strong evidence that the diffusion mechanism is (1) triggered by the growth of the tunnel junction cathode and (2) involves point defects. We analyze the implications of Zn-diffusion on the bandgap, the rear-passivation and the minority carrier quality of the GaInP solar subcell by relating the electrical performance of different samples to its corresponding doping profile.


title: On the Phases of a Semi-Sectorial Matrix

abstract: In this paper, we extend the definition of phases of sectorial matrices to those of semi-sectorial matrices, which are possibly singular. Properties of the phases are also extended, including those of the Moore-Penrose generalized inverse, compressions and Schur complements, matrix sums and products. In particular, a majorization relation is established between the phases of the nonzero eigenvalues of $AB$ and the phases of the compressions of $A$ and $B$, which leads to a generalized matrix small phase theorem. For the matrices which are not necessarily semi-sectorial, we define their (largest and smallest) essential phases via diagonal similarity transformation. An explicit expression for the essential phases of a Laplacian matrix of a directed graph is obtained.


title: Multivector Contractions Revisited, Part I

abstract: We reorganize, simplify and expand the theory of contractions or interior products of multivectors, and related topics like Hodge star duality. Many results are generalized and new ones are given, like: geometric characterizations of blade contractions and regressive products, higher-order graded Leibniz rules, determinant formulas, improved complex star operators, etc. Different contractions and conventions found in the literature are discussed and compared, in special those of Clifford Geometric Algebra. Applications of the theory are developed in a follow-up paper.


title: Reduction of detection limit and quantification uncertainty due to interferent by neural classification with abstention

abstract: Many measurements in the physical sciences can be cast as counting experiments, where the number of occurrences of a physical phenomenon informs the prevalence of the phenomenon's source. Often, detection of the physical phenomenon (termed signal) is difficult to distinguish from naturally occurring phenomena (termed background). In this case, the discrimination of signal events from background can be performed using classifiers, and they may range from simple, threshold-based classifiers to sophisticated neural networks. These classifiers are often trained and validated to obtain optimal accuracy, however we show that the optimal accuracy classifier does not generally coincide with a classifier that provides the lowest detection limit, nor the lowest quantification uncertainty. We present a derivation of the detection limit and quantification uncertainty in the classifier-based counting experiment case. We also present a novel abstention mechanism to minimize the detection limit or quantification uncertainty \emph{a posteriori}. We illustrate the method on two data sets from the physical sciences, discriminating Ar-37 and Ar-39 radioactive decay from non-radioactive events in a gas proportional counter, and discriminating neutrons from photons in an inorganic scintillator and report results therefrom.

categories: physics.data-an cs.LG hep-ex physics.ins-det


title: AnySeq/GPU: A Novel Approach for Faster Sequence Alignment on GPUs

abstract: In recent years, the rapidly increasing number of reads produced by next-generation sequencing (NGS) technologies has driven the demand for efficient implementations of sequence alignments in bioinformatics. However, current state-of-the-art approaches are not able to leverage the massively parallel processing capabilities of modern GPUs with close-to-peak performance. We present AnySeq/GPU-a sequence alignment library that augments the AnySeq1 library with a novel approach for accelerating dynamic programming (DP) alignment on GPUs by minimizing memory accesses using warp shuffles and half-precision arithmetic. Our implementation is based on the AnyDSL compiler framework which allows for convenient zero-cost abstractions through guaranteed partial evaluation. We show that our approach achieves over 80% of the peak performance on both NVIDIA and AMD GPUs thereby outperforming the GPU-based alignment libraries AnySeq1, GASAL2, ADEPT, and NVBIO by a factor of at least 3.6 while achieving a median speedup of 19.2x over these tools across different alignment scenarios and sequence lengths when running on the same hardware. This leads to throughputs of up to 1.7 TCUPS (tera cell updates per second) on an NVIDIA GV100, up to 3.3 TCUPS with half-precision arithmetic on a single NVIDIA A100, and up to 3.8 TCUPS on an AMD MI100.


title: Noise-Tolerant Learning for Audio-Visual Action Recognition

abstract: Recently, video recognition is emerging with the help of multi-modal learning, which focuses on integrating distinct modalities to improve the performance or robustness of the model. Although various multi-modal learning methods have been proposed and offer remarkable recognition results, almost all of these methods rely on high-quality manual annotations and assume that modalities among multi-modal data provide semantically relevant information. Unfortunately, the widely used video datasets are usually coarse-annotated or collected from the Internet. Thus, it inevitably contains a portion of noisy labels and noisy correspondence. To address this challenge, we use the audio-visual action recognition task as a proxy and propose a noise-tolerant learning framework to find anti-interference model parameters against both noisy labels and noisy correspondence. Specifically, our method consists of two phases that aim to rectify noise by the inherent correlation between modalities. First, a noise-tolerant contrastive training phase is performed to make the model immune to the possible noisy-labeled data. To alleviate the influence of noisy correspondence, we propose a cross-modal noise estimation component to adjust the consistency between different modalities. As the noisy correspondence existed at the instance level, we further propose a category-level contrastive loss to reduce its interference. Second, in the hybrid-supervised training phase, we calculate the distance metric among features to obtain corrected labels, which are used as complementary supervision to guide the training. Extensive experiments on a wide range of noisy levels demonstrate that our method significantly improves the robustness of the action recognition model and surpasses the baselines by a clear margin.


title: Modeling the Optical to Ultraviolet Polarimetric Variability from Thomson Scattering in Colliding Wind Binaries

abstract: Massive star binaries are critical laboratories for measuring masses and stellar wind mass-loss rates. A major challenge is inferring viewing inclination and extracting information about the colliding wind interaction (CWI) region. Polarimetric variability from electron scattering in the highly ionized winds provides important diagnostic information about system geometry. We combine for the first time the well-known generalized treatment of \citet{brown_polarisation_1978} for variable polarization from binaries with the semi-analytic solution for the geometry and surface density CWI shock interface between the winds based on Canto et al 1996. Our calculations include some simplifications in the form of inverse square-law wind densities and the assumption of axisymmetry, but in so doing arrive at several robust conclusions. One is that when the winds are nearly equal (e.g., O\,+\,O binaries), the polarization has a relatively mild decline with binary separation. Another is that despite Thomson scattering being a gray opacity, the continuum polarization can show chromatic effects at ultraviolet wavelengths but will be mostly constant at longer wavelengths. Finally, when one wind dominates the other, as for example in WR+OB binaries, the polarization is expected to be larger at wavelengths where the OB component is more luminous, and generally smaller at wavelengths where the WR component is more luminous. This behavior arises because from the perspective of the WR star, the distortion of the scattering envelope from spherical is a minor perturbation situated far from the WR star. By contrast, the polarization contribution from the OB star is dominated by the geometry of the CWI shock.


title: Scalable Vehicle Re-Identification via Self-Supervision

abstract: As Computer Vision technologies become more mature for intelligent transportation applications, it is time to ask how efficient and scalable they are for large-scale and real-time deployment. Among these technologies is Vehicle Re-Identification which is one of the key elements in city-scale vehicle analytics systems. Many state-of-the-art solutions for vehicle re-id mostly focus on improving the accuracy on existing re-id benchmarks and often ignore computational complexity. To balance the demands of accuracy and computational efficiency, in this work we propose a simple yet effective hybrid solution empowered by self-supervised training which only uses a single network during inference time and is free of intricate and computation-demanding add-on modules often seen in state-of-the-art approaches. Through extensive experiments, we show our approach, termed Self-Supervised and Boosted VEhicle Re-Identification (SSBVER), is on par with state-of-the-art alternatives in terms of accuracy without introducing any additional overhead during deployment. Additionally we show that our approach, generalizes to different backbone architectures which facilitates various resource constraints and consistently results in a significant accuracy boost.


title: Rethinking Reinforcement Learning based Logic Synthesis

abstract: Recently, reinforcement learning has been used to address logic synthesis by formulating the operator sequence optimization problem as a Markov decision process. However, through extensive experiments, we find out that the learned policy makes decisions independent from the circuit features (i.e., states) and yields an operator sequence that is permutation invariant to some extent in terms of operators. Based on these findings, we develop a new RL-based method that can automatically recognize critical operators and generate common operator sequences generalizable to unseen circuits. Our algorithm is verified on both the EPFL benchmark, a private dataset and a circuit at industrial scale. Experimental results demonstrate that it achieves a good balance among delay, area and runtime, and is practical for industrial usage.


title: Approximate dynamic programming for profit estimation of connected hydro reservoirs

abstract: In this paper, we study the operational problem of connected hydro power reservoirs which involves sequential decision-making in an uncertain and dynamic environment. The problem is traditionally formulated as a stochastic dynamic program accounting for the uncertainty of electricity prices and reservoir inflows. This formulation suffers from the curse of dimensionality, as the state space explodes with the number of reservoirs and the history of prices and inflows. To avoid computing the expectation of future value functions, the proposed model takes advantage of the so-called post-decision state. To further tackle the dimensionality issue, we propose an approximate dynamic programming approach that estimates the future value of water using a linear approximation architecture. When the time series of prices and inflows follow autoregressive processes, our approximation provides an upper bound on the future value function. We use an offline training algorithm based on the historical data of prices and inflows and run both in-sample and out-of-sample simulations. Two realistic test systems of cascade and network connected reservoirs serve to demonstrate the computational tractability of our approach. In particular, we provide numerical evidence of convergence and quality of solutions. For our test systems, our results show that profit estimation is improved by 20% when including inflows in the linear approximation.


title: The logic of planetary combination in Vettius Valens

abstract: The Anthologies of the second-century astrologer Vettius Valens (120-c.175 CE) is the most extensive surviving practical astrological text from the period. Despite this, the theoretical underpinnings of the Anthologies have been understudied; in general, the work has been overshadowed by Ptolemy's contemporaneous Tetrabiblos. While the Tetrabiblos explicitly aims to present a systematic account of astrology, Valens' work is often characterised as a miscellaneous collection, of interest to historians only for the evidence it preserves about the practical methods used in casting horoscopes. In this article, we argue that the Anthologies is also an invaluable resource for engagement with the conceptual basis of astrology. As a case study, we take a section of Anthologies Book 1 which lists the possible astrological effects of planets, both alone and in 'combinations' of two and three. We demonstrate that analysing Valens' descriptions quantitatively with textual analysis reveals a consistent internal logic of planetary combination. By classifying descriptive terms as positive or negative, we show that the resulting 'sentiment' of planetary combinations is well-correlated with their component parts. Furthermore, we find that the sentiment of three-planet combinations is more strongly correlated with the average sentiment of their three possible component pairs than with the average sentiment of individual planets, suggesting an iterative combinatorial logic. Recognition of this feature of astrological practice has been neglected compared to the mathematical methods for calculating horoscopes. We argue that this analysis not only provides evidence that the astrological lore detailed in Valens is more consistent than is often assumed, but is also indicative of a wider methodological technique in practical astrology: combinatorial reasoning from existing astrological lore.


title: Analysis of Distributed Ledger Technologies for Industrial Manufacturing

abstract: In recent years, industrial manufacturing has undergone massive technological changes that embrace digitalization and automation towards the vision of intelligent manufacturing plants. With the aim of maximizing efficiency and profitability in production, an important goal is to enable flexible manufacturing, both, for the customer (desiring more individualized products) and for the manufacturer (to adjust to market demands). Manufacturing-as-a-service can support this through manufacturing plants that are used by different tenants who utilize the machines in the plant, which are offered by different providers. To enable such pay-per-use business models, Distributed Ledger Technology (DLT) is a viable option to establish decentralized trust and traceability. Thus, in this paper, we study potential DLT technologies for an efficient and intelligent integration of DLT-based solutions in manufacturing environments. We propose a general framework to adapt DLT in manufacturing, then we introduce the use case of shared manufacturing, which we utilize to study the communication and computation efficiency of selected DLTs in resource-constrained wireless IoT networks.


title: CGR-CUSUM: A Continuous time Generalized Rapid Response Cumulative Sum chart

abstract: Rapidly detecting problems in the quality of care is of utmost importance for the well-being of patients. Without proper inspection schemes, such problems can go undetected for years. Cumulative sum (CUSUM) charts have proven to be useful for quality control, yet available methodology for survival outcomes is limited. The few available continuous time inspection charts usually require the researcher to specify an expected increase in the failure rate in advance, thereby requiring prior knowledge about the problem at hand. Misspecifying parameters can lead to false positive alerts and large detection delays. To solve this problem, we take a more general approach to derive the new Continuous time Generalized Rapid response CUSUM (CGR-CUSUM) chart. We find an expression for the approximate average run length (average time to detection) and illustrate the possible gain in detection speed by using the CGR-CUSUM over other commonly used monitoring schemes on a real life data set from the Dutch Arthroplasty Register as well as in simulation studies. Besides the inspection of medical procedures, the CGR-CUSUM can also be used for other real time inspection schemes such as industrial production lines and quality control of services.


title: Wall-induced translation of a rotating particle in viscoelastic fluid

abstract: Shear-thinning and viscoelasticity are two non-Newtonian fluid properties widely existing in biological fluids. In this study, we found that the translation motion of a rotating particle near a wall speed up firstly, and then slows down with enhancement of fluid viscoelasticity, which is different from the behavior reported in shear thinning fluid (Chen et al. J. Fluid Mech. 2021, 927). Our research is carried out by numerical simulation of Navier-Stokes equations combined with Oldroyd-B constitutive model. This work is expected to be helpful to understand the movement of a rotating sphere near a wall in complex fluids comprehensively.


title: Bidirectional Optimisation for Load Shaping within Coupled Microgrids

abstract: We address the problem of load shaping within a network of coupled microgrids (MGs) in a bilevel optimisation framework. To this end, we consider the charging/discharging rates of residential energy storage devices within each MG on the lower level and the power exchange among neighbouring MGs on the upper level as optimisation variables. We improve a previously developed model such that the maximal amount of exchanged power does not depend on the power demand, thus, increasing the flexibility within the network, and adapt the corresponding bidirectional optimisation scheme accordingly. For efficiency, standard distributed optimisation routines are used for the optimisation on the lower level; the power exchange problem on the upper level is replaced by parallelisable small-scale quadratic programmings. We prove global convergence of the optimisation scheme and illustrate the potential of the approach in a numerical case study based on real-world data.


title: Three-dimensional buoyant hydraulic fracture growth: constant release from a point source

abstract: Hydraulic fractures propagating at depth are subjected to buoyant forces caused by the density contrast between fluid and solid. This paper is concerned with the analysis of the transition from an initially radial towards an elongated buoyant growth -- a critical topic for understanding the extent of vertical hydraulic fractures in the upper Earth crust. Using fully coupled numerical simulations and scaling arguments, we show that a single dimensionless number governs buoyant hydraulic fracture growth: the dimensionless viscosity of a radial hydraulic fracture at the time when buoyancy becomes of order one. It quantifies if the transition to buoyancy occurs when the growth of the radial hydraulic fracture is either still in the regime dominated by viscous flow dissipation or is already in the regime where fracture energy dissipation dominates. A family of fracture shapes emerge at late time from finger-like (toughness regime) to inverted elongated cudgel-like (viscous regime). 3D toughness dominated buoyant fractures exhibit a finger-like shape with a constant volume toughness dominated head and a viscous tail having a constant uniform horizontal breadth: there is no further horizontal growth past the onset of buoyancy. However, if the transition to buoyancy occurs while in the viscosity dominated regime, both vertical and horizontal growths continue to match scaling arguments. As soon as the fracture toughness is not strictly zero, horizontal growth stops when the dimensionless horizontal toughness becomes of order one. The horizontal breadth follows the predicted scaling.


title: Photoelastic Stress Response of Complex 3D-Printed Particle Shapes

abstract: While stress visualization within 3-dimensional particles would greatly advance our understanding of the behaviors of complex particles, traditional photoelastic methods suffer from a lack of available technology for producing suitable complex particles. Recently, 3D-printing has created new possibilities for enhancing the scope of stress analysis within physically representative granules. Here, we investigate and evaluate opportunities offered by 3D-printing a single particle with a complex external shape with photoelastic properties. We report the results of X-ray computed tomography and 3D-printing, combined with traditional photoelastic analysis, to visualize strain for particles ranging from simple 2D discs to complex 3D printed coffee beans, including with internal voids. We find that the relative orientation of the print layers and the loading force affects the optical response of the discs, but without a significant difference in their mechanical properties. Furthermore, we present semi-quantitative measurements of stresses within 3D-printed complex particles. The paper outlines the potential limitations and areas of future interest for stress visualization of 3-dimensional particles.


title: Model Agnostic Local Explanations of Reject

abstract: The application of machine learning based decision making systems in safety critical areas requires reliable high certainty predictions. Reject options are a common way of ensuring a sufficiently high certainty of predictions made by the system. While being able to reject uncertain samples is important, it is also of importance to be able to explain why a particular sample was rejected. However, explaining general reject options is still an open problem. We propose a model agnostic method for locally explaining arbitrary reject options by means of interpretable models and counterfactual explanations.


title: Opinion polarization in human communities can emerge as a natural consequence of beliefs being interrelated

abstract: The emergence of opinion polarization within human communities -- the phenomenon that individuals within a society tend to develop conflicting attitudes related to the greatest diversity of topics -- has been a focus of interest for decades, both from theoretical and modelling points of view. Regarding modelling attempts, an entire scientific field -- opinion dynamics -- has emerged in order to study this and related phenomena. Within this framework, agents' opinions are usually represented by a scalar value which undergoes modification due to interaction with other agents. Under certain conditions, these models are able to reproduce polarization -- a state increasingly familiar to our everyday experience. In the present paper, an alternative explanation is suggested along with its corresponding model. More specifically, we demonstrate that by incorporating the following two well-known human characteristics into the representation of agents: (1) in the human brain beliefs are interconnected, and (2) people strive to maintain a coherent belief system; polarization immediately occurs under exposure to news and information. Furthermore, the model accounts for the proliferation of fake news, and shows how opinion polarization is related to various cognitive biases.


title: Gaussian Processes and Bayesian Optimization for High Precision Experiments

abstract: High-precision measurements require optimal setups and analysis tools to achieve continuous improvements. Systematic corrections need to be modeled with high accuracy and known uncertainty to reconstruct underlying physical phenomena. To this end, we present Gaussian processes for modeling experiments and usage with Bayesian optimization, on the example of an electron energy detector, achieving optimal performance. We demonstrate the method's strengths and outline stochastic variational Gaussian processes for physics applications with large data sets, enabling new solutions for current problems.


title: Attacking and Defending Deep Reinforcement Learning Policies

abstract: Recent studies have shown that deep reinforcement learning (DRL) policies are vulnerable to adversarial attacks, which raise concerns about applications of DRL to safety-critical systems. In this work, we adopt a principled way and study the robustness of DRL policies to adversarial attacks from the perspective of robust optimization. Within the framework of robust optimization, optimal adversarial attacks are given by minimizing the expected return of the policy, and correspondingly a good defense mechanism should be realized by improving the worst-case performance of the policy. Considering that attackers generally have no access to the training environment, we propose a greedy attack algorithm, which tries to minimize the expected return of the policy without interacting with the environment, and a defense algorithm, which performs adversarial training in a max-min form. Experiments on Atari game environments show that our attack algorithm is more effective and leads to worse return of the policy than existing attack algorithms, and our defense algorithm yields policies more robust than existing defense methods to a range of adversarial attacks (including our proposed attack algorithm).


title: KnowGraph-PM: a Knowledge Graph based Pricing Model for Semiconductors Supply Chains

abstract: Semiconductor supply chains are described by significant demand fluctuation that increases as one moves up the supply chain, the so-called bullwhip effect. To counteract, semiconductor manufacturers aim to optimize capacity utilization, to deliver with shorter lead times and exploit this to generate revenue. Additionally, in a competitive market, firms seek to maintain customer relationships while applying revenue management strategies such as dynamic pricing. Price change potentially generates conflicts with customers. In this paper, we present KnowGraph-PM, a knowledge graph-based dynamic pricing model. The semantic model uses the potential of faster delivery and shorter lead times to define premium prices, thus entail increased profits based on the customer profile. The knowledge graph enables the integration of customer-related information, e.g., customer class and location to customer order data. The pricing algorithm is realized as a SPARQL query that relies on customer profile and order behavior to determine the corresponding price premium. We evaluate the approach by calculating the revenue generated after applying the pricing algorithm. Based on competency questions that translate to SPARQL queries, we validate the created knowledge graph. We demonstrate that semantic data integration enables customer-tailored revenue management.


title: Thermodynamics as Combinatorics: A Toy Theory

abstract: We discuss a simple toy model which allows, in a natural way, for deriving central facts from thermodynamics such as its fundamental laws, including Carnot's version of the second principle. Our viewpoint represents thermodynamic systems as binary strings, and it links their temperature to their Hamming weight. From this, we can reproduce the possibility of negative temperatures, the notion of equilibrium as the co\"incidence of two notions of temperature - statistical versus structural -, as well as the zeroth law of thermodynamics (transitivity of the thermal-equilibrium relation), which we find to be redundant, as other authors, yet at the same time not to be universally valid.


title: Testing Born-Infeld $f(T)$ teleparallel gravity through Sgr A$^\star$ observations

abstract: We use observational data from the S2 star orbiting around the Galactic Center to constrain a black hole solution of extended teleparallel gravity models. Subsequently, we construct the shadow images of Sgr A$^{\star}$ black hole. In particular, we constrain the parameter $\alpha=1/\lambda$ which appears in the Born-Infeld $f(T)$ model. In the strong gravity regime we find that the shadow radius increases with the increase of the parameter $\alpha$. Specifically, from the S2 star observations, we find within 1$\sigma$ that the parameter must lie between $0 \leq \alpha/M^2 \leq 6 \times 10^{-4}$. Consequently, we used the best fit parameters to model the shadow images of Sgr A$^{\star}$ black hole and then using the Gauss-Bonnet theorem we analysed the deflection angle for leading order expansions of the parameter $\alpha$. It is found that within the parameter range, these observables are very close to the Schwarzschild case. Furthermore, using the best fit parameters for the Born-Infeld $f(T)$ model we show the angular diameter is consistent with recent observations for the Sgr A$^{\star}$ black hole angular diameter $(51.8 \pm 2.3) \mu$arcsec and difficult to be distinguished from the GR. For the deflection angle of light, in leading order terms, we find that the deflection angle expressed in the ADM mass coincides with the GR, but the ADM mass in the Born-Infeld $f(T)$ gravity increases with the increase of $\alpha$ and the overall deflection angle is expected to me greater in $f(T)$ gravity. As a consequence of this fact, we have shown that the electromagnetic intensity observed in shadow images is smaller compared to GR.


title: Analysis of The Vehicle Routing Problem Solved via Hybrid Quantum Algorithms in Presence of Noisy Channels

abstract: The vehicle routing problem (VRP) is an NP-hard optimization problem that has been an interest of research for decades in science and industry. The objective is to plan routes of vehicles to deliver goods to a fixed number of customers with optimal efficiency. Classical tools and methods provide good approximations to reach the optimal global solution. Quantum computing and quantum machine learning provide a new approach to solving combinatorial optimization of problems faster due to inherent speedups of quantum effects. Many solutions of VRP are offered across different quantum computing platforms using hybrid algorithms such as quantum approximate optimization algorithm and quadratic unconstrained binary optimization. In this work, we build a basic VRP solver for 3 and 4 cities using the variational quantum eigensolver on a fixed ansatz. The work is further extended to evaluate the robustness of the solution in several examples of noisy quantum channels. We find that the performance of the quantum algorithm depends heavily on what noise model is used. In general, noise is detrimental, but not equally so among different noise sources.


title: A bootstrap approach for validating the number of groups identified by latent class growth models

abstract: The use of longitudinal finite mixture models such as group-based trajectory modeling has seen a sharp increase during the last decades in the medical literature. However, these methods have been criticized especially because of the data-driven modelling process which involves statistical decision-making. In this paper, we propose an approach that uses bootstrap to sample observations with replacement from the original data to validate the number of groups identified and to quantify the uncertainty in the number of groups. The method allows investigating the statistical validity and the uncertainty of the groups identified in the original data by checking if the same solution is also found across the bootstrap samples. In a simulation study, we examined whether the bootstrap-estimated variability in the number of groups reflected the replication-wise variability. We also compared the replication-wise variability to the Bayesian posterior probability. We evaluated the ability of three commonly used adequacy criteria (average posterior probability, odds of correct classification and relative entropy) to identify uncertainty in the number of groups. Finally, we illustrated the proposed approach using data from the Quebec Integrated Chronic Disease Surveillance System to identify longitudinal medication patterns between 2015 and 2018 in older adults with diabetes.


title: Aphorisms on Epidemiological Modelling

abstract: Epidemiological modelling is critiqued towards a scientific practice of negativity in the context of Scotland's Centre of Expertise on Animal Disease Outbreaks (EPIC). The paratactical approach to the melancholy science is invoked to problematise One Health, the intra-pandemic modelling culture, and to delineate an inkling of the negative in EPIC's work.


title: Taming Continuous Posteriors for Latent Variational Dialogue Policies

abstract: Utilizing amortized variational inference for latent-action reinforcement learning (RL) has been shown to be an effective approach in Task-oriented Dialogue (ToD) systems for optimizing dialogue success. Until now, categorical posteriors have been argued to be one of the main drivers of performance. In this work we revisit Gaussian variational posteriors for latent-action RL and show that they can yield even better performance than categoricals. We achieve this by simplifying the training procedure and propose ways to regularize the latent dialogue policy to retain good response coherence. Using continuous latent representations our model achieves state of the art dialogue success rate on the MultiWOZ benchmark, and also compares well to categorical latent methods in response coherence.


title: A Precis of Language Models are not Models of Language

abstract: Natural Language Processing is one of the leading application areas in the current resurgence of Artificial Intelligence, spearheaded by Artificial Neural Networks. We show that despite their many successes at performing linguistic tasks, Large Neural Language Models are ill-suited as comprehensive models of natural language. The wider implication is that, in spite of the often overbearing optimism about AI, modern neural models do not represent a revolution in our understanding of cognition.


title: Relating Information and Proof

abstract: In mathematics information is a number that measures uncertainty (entropy) based on a probabilistic distribution, often of an obscure origin. In real life language information is a datum, a statement, more precisely, a formula. But such a formula should be justified by a proof. I try to formalize this perception of information. The measure of informativeness of a proof is based on the set of proofs related to the formulas under consideration. This set of possible proofs (`a knowledge base') defines a probabilistic measure, and entropic weight is defined using this measure. The paper is mainly conceptual, it is not clear where and how this approach can be applied.


title: A review of ontologies for smart and continuous commissioning

abstract: Smart and continuous commissioning (SCCx) of buildings can result in a significant reduction in the gap between design and operational performance. Ontologies play an important role in SCCx as they facilitate data readability and reasoning by machines. A better understanding of ontologies is required in order to develop and incorporate them in SCCx. This paper critically reviews the state-of-the-art research on building data ontologies since 2014 within the SCCx domain through sorting them based on building data types, general approaches, and applications. The data types of two main domains of building information modeling and building management system have been considered in the majority of existing ontologies. Three main applications are evident from a critical analysis of existing ontologies: (1) key performance indicator calculation, (2) building performance improvement, and (3) fault detection and diagnosis. The key gaps found in the literature review are a holistic ontology for SCCx and insight on how such approaches should be evaluated. Based on these findings, this study provides recommendations for future necessary research including: identification of SCCx-related data types, assessment of ontology performance, and creation of open-source approaches.


title: MATLAB implementation of hp finite elements on rectangles

abstract: A simple MATLAB implementation of hierarchical shape functions on 2D rectangles is explained and available for download. Global shape functions are ordered for a given polynomial degree according to the indices of the nodes, edges, or elements to which they belong. For a uniform p-refinement, the hierarchical structure enables an effective assembly of mass and stiffness matrices. A solution of a boundary value problem is approximated for various levels of uniform h and p refinements.


title: EFX Allocations: Simplifications and Improvements

abstract: The existence of EFX allocations is a fundamental open problem in discrete fair division. Given a set of agents and indivisible goods, the goal is to determine the existence of an allocation where no agent envies another following the removal of any single good from the other agent's bundle. Since the general problem has been illusive, progress is made on two fronts: $(i)$ proving existence when the number of agents is small, $(ii)$ proving existence of relaxations of EFX. In this paper, we improve results on both fronts (and simplify in one of the cases). We prove the existence of EFX allocations with three agents, restricting only one agent to have an MMS-feasible valuation function (a strict generalization of nice-cancelable valuation functions introduced by Berger et al. which subsumes additive, budget-additive and unit demand valuation functions). The other agents may have any monotone valuation functions. Our proof technique is significantly simpler and shorter than the proof by Chaudhury et al. on existence of EFX allocations when there are three agents with additive valuation functions and therefore more accessible. Secondly, we consider relaxations of EFX allocations, namely, approximate-EFX allocations and EFX allocations with few unallocated goods (charity). Chaudhury et al. showed the existence of $(1-\epsilon)$-EFX allocation with $O((n/\epsilon)^{\frac{4}{5}})$ charity by establishing a connection to a problem in extremal combinatorics. We improve their result and prove the existence of $(1-\epsilon)$-EFX allocations with $\tilde{O}((n/ \epsilon)^{\frac{1}{2}})$ charity. In fact, some of our techniques can be used to prove improved upper-bounds on a problem in zero-sum combinatorics introduced by Alon and Krivelevich.


title: Distribution Estimation for Probabilistic Loops

abstract: We present an algorithmic approach to estimate the value distributions of random variables of probabilistic loops whose statistical moments are (partially) known. Based on these moments, we apply two statistical methods, Maximum Entropy and Gram-Charlier series, to estimate the distributions of the loop's random variables. We measure the accuracy of our distribution estimation by comparing the resulting distributions using exact and estimated moments of the probabilistic loop, and performing statistical tests. We evaluate our method on several probabilistic loops with polynomial updates over random variables drawing from common probability distributions, including examples implementing financial and biological models. For this, we leverage symbolic approaches to compute exact higher-order moments of loops as well as use sampling-based techniques to estimate moments from loop executions. Our experimental results provide practical evidence of the accuracy of our method for estimating distributions of probabilistic loop outputs.


title: ecpc: An R-package for generic co-data models for high-dimensional prediction

abstract: High-dimensional prediction considers data with more variables than samples. Generic research goals are to find the best predictor or to select variables. Results may be improved by exploiting prior information in the form of co-data, providing complementary data not on the samples, but on the variables. We consider adaptive ridge penalised generalised linear and Cox models, in which the variable specific ridge penalties are adapted to the co-data to give a priori more weight to more important variables. The R-package ecpc originally accommodated various and possibly multiple co-data sources, including categorical co-data, i.e. groups of variables, and continuous co-data. Continuous co-data, however, was handled by adaptive discretisation, potentially inefficiently modelling and losing information. Here, we present an extension to the method and software for generic co-data models, particularly for continuous co-data. At the basis lies a classical linear regression model, regressing prior variance weights on the co-data. Co-data variables are then estimated with empirical Bayes moment estimation. After placing the estimation procedure in the classical regression framework, extension to generalised additive and shape constrained co-data models is straightforward. Besides, we show how ridge penalties may be transformed to elastic net penalties with the R-package squeezy. In simulation studies we first compare various co-data models for continuous co-data from the extension to the original method. Secondly, we compare variable selection performance to other variable selection methods. Moreover, we demonstrate use of the package in several examples throughout the paper.


title: Robust structured light in atmospheric turbulence

abstract: Structured light is routinely used in free space optical communication channels, both classical and quantum, where information is encoded in the spatial structure of the mode for increased bandwidth. Unlike polarisation, the spatial structure of light is perturbed through such channels by atmospheric turbulence, and consequently, much attention has focused on whether one mode type is more robust than another, but with seemingly inconclusive and contradictory results. Both real-world and experimentally simulated turbulence conditions have revealed that free-space structured light modes are perturbed in some manner by turbulence, resulting in both amplitude and phase distortions. Here, we present complex forms of structured light which are invariant under propagation through the atmosphere: the true eigenmodes of atmospheric turbulence. We provide a theoretical procedure for obtaining these eigenmodes and confirm their invariance both numerically and experimentally. Although we have demonstrated the approach on atmospheric turbulence, its generality allows it to be extended to other channels too, such as underwater and in optical fibre.


title: Black strings from dark matter

abstract: In this paper, we obtain two different static black string solutions by considering as sources axisymmetric dark matter distributions in 3+1 dimensions. These solutions tend asymptotically to the usual static and uncharged black string vacuum solution predicted by General Relativity (GR). We show that both the solutions present an event horizon each, like the vacuum solution, which is larger than the horizon of the latter. Then, we obtain the Hawking temperature associated with the black string solutions. Differently from what occurs with the static black string in the vacuum, we find that there exists a linear density of mass (or tension) remnant associated with a vanishing Hawking temperature for the obtained solutions. Thus, we analyze how the presence of dark matter affects the occurrence of the remnants. Further, we calculate other thermodynamic quantities, namely entropy, heat capacity, and free energy per length unit, showing that thermal phase transitions can occur in the presence of dark matter. We also analyze the weak (and null) energy conditions and conclude that the dark matter does not behave like an exotic fluid. Finally, we obtain the corresponding stationary solutions, determining their tensions as functions of both the mass and angular momentum of the black strings.


title: A boundary maximum principle for stationary pairs of varifolds with fixed contact angle

abstract: In this note, we establish a boundary maximum principle for a class of stationary pairs of varifolds satisfying a fixed contact angle condition in any compact Riemannian manifold with smooth boundary.


title: Localization of n-exangulated categories

abstract: Nakaoka-Ogawa-Sakai considered the localization of an extriangulated category. This construction unified the Serre quotient of abelian categories and the Verdier quotient of triangulated categories. Recently, Herschend-Liu-Nakaoka defined $n$-exangulated categories as a higher dimensional analogue of extriangulated categories. Let $\mathcal C$ be an $n$-exangulated category and $\mathcal{F}$ be a multiplicative system satisfying mild assumption. In this article, we give a necessary and sufficient condition for the localization of $\mathcal C$ be an $n$-exangulated category. This way gives a new class of $n$-exangulated categories which are neither $n$-exact nor $(n+2)$-angulated in general. Moreover, our result also generalizes work by Nakaoka-Ogawa-Sakai.


title: Dual pairs of operators, harmonic analysis of singular non-atomic measures and Krein-Feller diffusion

abstract: We show that a Krein-Feller operator is naturally associated to a fixed measure $\mu$, assumed positive, $\sigma$-finite, and non-atomic. Dual pairs of operators are introduced, carried by the two Hilbert spaces, $L^{2}\left(\mu\right)$ and $L^{2}\left(\lambda\right)$, where $\lambda$ denotes Lebesgue measure. An associated operator pair consists of two specific densely defined (unbounded) operators, each one contained in the adjoint of the other. This then yields a rigorous analysis of the corresponding $\mu$-Krein-Feller operator as a closable quadratic form. As an application, for a given measure $\mu$, including the case of fractal measures, we compute the associated diffusion, semigroup, Dirichlet forms, and $\mu$-generalized heat equation.


title: A Fast Attention Network for Joint Intent Detection and Slot Filling on Edge Devices

abstract: Intent detection and slot filling are two main tasks in natural language understanding and play an essential role in task-oriented dialogue systems. The joint learning of both tasks can improve inference accuracy and is popular in recent works. However, most joint models ignore the inference latency and cannot meet the need to deploy dialogue systems at the edge. In this paper, we propose a Fast Attention Network (FAN) for joint intent detection and slot filling tasks, guaranteeing both accuracy and latency. Specifically, we introduce a clean and parameter-refined attention module to enhance the information exchange between intent and slot, improving semantic accuracy by more than 2%. FAN can be implemented on different encoders and delivers more accurate models at every speed level. Our experiments on the Jetson Nano platform show that FAN inferences fifteen utterances per second with a small accuracy drop, showing its effectiveness and efficiency on edge devices.


title: Pareto-Optimal Fairness-Utility Amortizations in Rankings with a DBN Exposure Model

abstract: In recent years, it has become clear that rankings delivered in many areas need not only be useful to the users but also respect fairness of exposure for the item producers. We consider the problem of finding ranking policies that achieve a Pareto-optimal tradeoff between these two aspects. Several methods were proposed to solve it; for instance a popular one is to use linear programming with a Birkhoff-von Neumann decomposition. These methods, however, are based on a classical Position Based exposure Model (PBM), which assumes independence between the items (hence the exposure only depends on the rank). In many applications, this assumption is unrealistic and the community increasingly moves towards considering other models that include dependences, such as the Dynamic Bayesian Network (DBN) exposure model. For such models, computing (exact) optimal fair ranking policies remains an open question. We answer this question by leveraging a new geometrical method based on the so-called expohedron proposed recently for the PBM (Kletti et al., WSDM'22). We lay out the structure of a new geometrical object (the DBN-expohedron), and propose for it a Carath\'eodory decomposition algorithm of complexity $O(n^3)$, where $n$ is the number of documents to rank. Such an algorithm enables expressing any feasible expected exposure vector as a distribution over at most $n$ rankings; furthermore we show that we can compute the whole set of Pareto-optimal expected exposure vectors with the same complexity $O(n^3)$. Our work constitutes the first exact algorithm able to efficiently find a Pareto-optimal distribution of rankings. It is applicable to a broad range of fairness notions, including classical notions of meritocratic and demographic fairness. We empirically evaluate our method on the TREC2020 and MSLR datasets and compare it to several baselines in terms of Pareto-optimality and speed.


title: Twisting the near-field radiative heat switch in hyperbolic antiferromagnets

abstract: We study the twisted control of the near-field radiative heat transfer between two hyperbolic antiferromagnetic insulators under external magnetic fields. We show that the near-field heat flux can be affected by both the twist angle $\theta$ and the magnitude of the applied magnetic field with different broken symmetries. Irrespective of twist angle, the external magnetic field causes the radiative heat flux to change nonmonotonically, and the minimum heat flux can be found with the magnetic fields of about 1.5 T. Such nonmonotonic behavior is due to the fact that the magnetic field can radically change the nature of the magnon polaritons with time reversal symmetry breaking. The field not only affects the topological structure of surface magnon polaritons, but also induces the volume magnon polaritons that progressively dominate the heat transfer as the field increases. We further propose a twist-induced thermal switch device with inversion symmetry breaking, which can severely regulate radiative heat flux through different magnetic fields. Our findings account for a characteristic modulation of radiative heat transfer with implications for applications in dynamic thermal management.


title: Generalizing to Evolving Domains with Latent Structure-Aware Sequential Autoencoder

abstract: Domain generalization aims to improve the generalization capability of machine learning systems to out-of-distribution (OOD) data. Existing domain generalization techniques embark upon stationary and discrete environments to tackle the generalization issue caused by OOD data. However, many real-world tasks in non-stationary environments (e.g. self-driven car system, sensor measures) involve more complex and continuously evolving domain drift, which raises new challenges for the problem of domain generalization. In this paper, we formulate the aforementioned setting as the problem of evolving domain generalization. Specifically, we propose to introduce a probabilistic framework called Latent Structure-aware Sequential Autoencoder (LSSAE) to tackle the problem of evolving domain generalization via exploring the underlying continuous structure in the latent space of deep neural networks, where we aim to identify two major factors namely covariate shift and concept shift accounting for distribution shift in non-stationary environments. Experimental results on both synthetic and real-world datasets show that LSSAE can lead to superior performances based on the evolving domain generalization setting.


title: The relation between a generalized Fibonacci sequence and the length of Cunningham chains

abstract: Let $p$ be a prime number. A chain $\{p,2p+1,4p+3,\cdots,(p+1)2^{l(p)-1}-1\}$ is called the Cunningham chain generated by $p$ if all elements are prime number and $(p+1)2^{l(p)}-1$ is composite. Then $l(p)$ is called the length of the Cunningham chain. It is conjectured by Bateman and Horn in 1962 that the number of prime $p\leq N$ such that $l(p)\geq k$ is asymptotically equal to $B_k N/(\log N)^k$ with a real $B_k>0$ for all natural number $k$. This suggests that $l(p)=\Omega(\log p/\log\log p)$. However, so far no good estimation is known. It has not even been proven whether $\limsup_{p\to\infty} l(p)$ is infinite or not. All we know is that $l(p)=5$ if $p=2$ and $l(p)<p$ for odd $p$ by Fermat's little theorem. Let $\alpha\geq3$ be an integer. In this article, a generalized Fibonacci sequence $\mathcal{F}_\alpha=\{F_n\}_{n=0}^\infty$ is defined as $F_0=0,F_1=1, F_{n+2}=\alpha F_{n+1}+F_n (n\geq0)$, and ${}_{\mathcal{F}_\alpha}\sigma(n)=\sum_{d\mid n, 0<d\in\mathcal{F}_\alpha}d$ is called a divisor function on $\mathcal{F}_\alpha$. Then we obtain an interesting relation between the iteration of ${}_{\mathcal{F}_\alpha}\sigma$ and the length of Cunningham chains. For two primes $p$ and $q$, the fact $p=2q+1$ or $2q-1$ is equivalent to ${}_{\mathcal{F}_\alpha}\sigma({}_{\mathcal{F}_\alpha}\sigma(F_p))={}_{\mathcal{F}_\alpha}\sigma(F_q)$ for some $\alpha$. By this relation, we get $l(p)\ll\log p$ under a certain condition. It seems that this sufficient condition is plausible by numerical test. Furthermore, the condition, written in terms of prime numbers, can be replaced by the condition written in terms of natural numbers. This implies that the problem of upper estimation of $l(p)$ is reduced to that on natural numbers.


title: Designing an Optimal LSST Deep Drilling Program for Cosmology with Type Ia Supernovae

abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time is forecast to collect a large sample of Type Ia supernovae (SNe Ia) that could be instrumental in unveiling the nature of Dark Energy. The feat, however, requires measuring the two components of the Hubble diagram - distance modulus and redshift - with a high degree of accuracy. Distance is estimated from SNe Ia parameters extracted from light curve fits, where the average quality of light curves is primarily driven by survey parameters such as the cadence and the number of visits per band. An optimal observing strategy is thus critical for measuring cosmological parameters with high accuracy. We present in this paper a three-stage analysis aiming at quantifying the impact of the Deep Drilling (DD) strategy parameters on three critical aspects of the survey: the redshift completeness (originating from the Malmquist cosmological bias), the number of well-measured SNe Ia, and the cosmological measurements. Analyzing the current LSST survey simulations, we demonstrate that the current DD survey plans are characterized by a low completeness ($z~\sim$ 0.55-0.65), and irregular and low cadences (few days) that dramatically decrease the size of the well-measured SNe Ia sample. We then propose a modus operandi that provides the number of visits (per band) required to reach higher redshifts. The results of this approach are used to design a set of optimized DD surveys for SNe Ia cosmology. We show that most accurate cosmological measurements are achieved with Deep Rolling surveys characterized by a high cadence (one day), a rolling strategy (each field observed at least two seasons), and two sets of fields: ultra-deep ($z \gtrsim 0.8$) and deep ($z \gtrsim 0.6$) fields. We also demonstrate that a deterministic scheduler including a gap recovery mechanism is critical to achieve a high quality DD survey required for SNe Ia cosmology.


title: Photosphere Recession and Luminosity of Homologous Explosions Revisited

abstract: By assuming the photosphere located at the outmost edge of the ejecta, Arnett et al. (1980, 1982, 1989) presented the light curves of homologous explosions in supernovae analytically and numerically to include recombination effects. Actually as homologous expansion proceeds, the photosphere recedes deeper into the ejecta. In this situation, the photosphere radius increases at early times and decreases later on which can be described by a simple method proposed by Liu et al. (2018). To study how the photosphere recession effect the luminosity evolution, we impose a boundary condition on the photosphere to determine the spatial and time distribution of the temperature of the ejecta which is clarified to be reasonable. We find that the photosphere recession reduce the luminosity compared with the previous result without the recession, which can be tested with observations of Type-IIP supernovae.


title: Long-term memory and synapse-like dynamics in two-dimensional nanofluidic channels

abstract: Fine-tuned ion transport across nanoscale pores is key to many biological processes such as neurotransmission. Recent advances have enabled the confinement of water and ions to two dimensions, unveiling transport properties unreachable at larger scales and triggering hopes to reproduce the ionic machinery of biological systems. Here we report experiments demonstrating the emergence of memory in the transport of aqueous electrolytes across (sub)nanoscale channels. We unveiled two types of nanofluidic memristors, depending on channel material and confinement, with memory from minutes to hours. We explained how large timescales could emerge from interfacial processes like ionic self-assembly or surface adsorption. Such behavior allowed us to implement Hebbian learning with nanofluidic systems. This result lays the ground for biomimetic computations on aqueous electrolytic chips.