title: Task-Oriented Optimal Sequencing of Visualization Charts

abstract: A chart sequence is used to describe a series of visualization charts generated in the exploratory analysis by data analysts. It provides information details in each chart as well as a logical relationship among charts. While existing research targets on generating chart sequences that match human's perceptions, little attention has been paid to formulate task-oriented connections between charts in a chart design space. We present a novel chart sequencing method based on reinforcement learning to capture the connections between charts in the context of three major analysis tasks, including correlation analysis, anomaly detection, and cluster analysis. The proposed method formulates a chart sequencing procedure as an optimization problem, which seeks an optimal policy to sequencing charts for the specific analysis task. In our method, a novel reward function is introduced, which takes both the analysis task and the factor of human cognition into consideration. We conducted one case study and two user studies to evaluate the effectiveness of our method under the application scenarios of visualization demonstration, sequencing charts for reasoning analysis results, and making a chart design choice. The study results showed the power of our method.


title: Computational approaches to non-convex, sparsity-inducing multi-penalty regularization

abstract: In this work we consider numerical efficiency and convergence rates for solvers of non-convex multi-penalty formulations when reconstructing sparse signals from noisy linear measurements. We extend an existing approach, based on reduction to an augmented single-penalty formulation, to the non-convex setting and discuss its computational intractability in large-scale applications. To circumvent this limitation, we propose an alternative single-penalty reduction based on infimal convolution that shares the benefits of the augmented approach but is computationally less dependent on the problem size. We provide linear convergence rates for both approaches, and their dependence on design parameters. Numerical experiments substantiate our theoretical findings.


title: Fusion of Sensors Data in Automotive Radar Systems: A Spectral Estimation Approach

abstract: To accurately estimate locations and velocities of surrounding targets (cars) is crucial for advanced driver assistance systems based on radar sensors. In this paper we derive methods for fusing data from multiple radar sensors in order to improve the accuracy and robustness of such estimates. First we pose the target estimation problem as a multivariate multidimensional spectral estimation problem. The problem is multivariate since each radar sensor gives rise to a measurement channel. Then we investigate how the use of the cross-spectra affects target estimates. We see that the use of the magnitude of the cross-spectrum significantly improves the accuracy of the target estimates, whereas an attempt to compensate the phase lag of the cross-spectrum only gives marginal improvement. This paper may be viewed as a first step towards applying high-resolution methods that builds on multidimensional multivariate spectral estimation for sensor fusion.


title: Embedding-based system for the Text part of CALL v3 shared task

abstract: This paper presents a scoring system that has shown the top result on the text subset of CALL v3 shared task. The presented system is based on text embeddings, namely NNLM~\cite{nnlm} and BERT~\cite{Bert}. The distinguishing feature of the given approach is that it does not rely on the reference grammar file for scoring. The model is compared against approaches that use the grammar file and proves the possibility to achieve similar and even higher results without a predefined set of correct answers. The paper describes the model itself and the data preparation process that played a crucial role in the model training.


title: High-temperature droplet epitaxy of symmetric GaAs/AlGaAs quantum dots

abstract: We introduce a high-temperature droplet epitaxy procedure, based on the control of the arsenization dynamics of nanoscale droplets of liquid Ga on GaAs(111)A surfaces. The use of high temperatures for the self-assembly of droplet epitaxy quantum dots solves major issues related to material defects, introduced during the droplet epitaxy fabrication process, which limited its use for single and entangled photon sources for quantum photonics applications. We identify the region in the parameter space which allows quantum dots to self-assemble with the desired emission wavelength and highly symmetric shape while maintaining a high optical quality. The role of the growth parameters during the droplet arsenization is discussed and modelled.


title: Mesh Variational Autoencoders with Edge Contraction Pooling

abstract: 3D shape analysis is an important research topic in computer vision and graphics. While existing methods have generalized image-based deep learning to meshes using graph-based convolutions, the lack of an effective pooling operation restricts the learning capability of their networks. In this paper, we propose a novel pooling operation for mesh datasets with the same connectivity but different geometry, by building a mesh hierarchy using mesh simplification. For this purpose, we develop a modified mesh simplification method to avoid generating highly irregularly sized triangles. Our pooling operation effectively encodes the correspondence between coarser and finer meshes in the hierarchy. We then present a variational auto-encoder structure with the edge contraction pooling and graph-based convolutions, to explore probability latent spaces of 3D surfaces. Our network requires far fewer parameters than the original mesh VAE and thus can handle denser models thanks to our new pooling operation and convolutional kernels. Our evaluation also shows that our method has better generalization ability and is more reliable in various applications, including shape generation, shape interpolation and shape embedding.


title: The Core-Cusp Problem Revisited: ULDM vs. CDM

abstract: The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that ULDM may resolve the core-cusp discrepancy associated with pure NFW halos without recourse to baryonic physics. However, it has been demonstrated that the density of ULDM halos can exceed those of comparable NFW configurations at some radii and halo masses, apparently exacerbating the problem rather than solving it. This situation arises because, although solitonic cores are flat at their centres, they obey an inverse mass-radius scaling relationship. Meanwhile, the mass of the inner soliton increases with the total halo mass, and therefore the inner core becomes more peaked at large halo masses. We describe a parameterisation of the radial density profiles of ULDM halos that allows for environmental variability of the core-halo mass relation in order to investigate this issue in more detail. For halos up to 10^12 solar masses we find feasible ULDM profiles for which the central density is lower than their NFW counterparts at astrophysically accessible radii. However, comparisons to observed profiles do not strongly favour either option; both give reasonable fits to subsets of the data for some parameter choices. Consequently, we find that robust tests of the core-cusp problem in ULDM will require more comprehensive observational data and simulations that include baryonic feedback.


title: Bringing Schrodinger's Cat to Life with Non-Equilibrium Respiration

abstract: In this study, we have proposed a method based on non-equilibrium effects to generate the superposition of macroscopically distinguishable quantum states, known as Schrodinger cat states, by using a Mach-Zehnder interferometry type experiment. Interaction of the input number state with a Kerr medium in the presence of a couple of heat baths in different temperatures in interaction picture and without imposing Markov assumption is considered. We have shown that the study of dynamics of the cat states under non-equilibrium condition open a way for the robustness of quantum features against the destructive role of the environment even at high temperature limit. It is verified that mutual influence of the environments, far from equilibrium, on the open system, makes it possible to revive quantum beats for longer time intervals. Moreover, we have probed how the traits of the environment, like its temperature and the Ohmic, super-Ohmic or sub-Ohmic functionality of the spectral density, may affect the pattern of the oscillation between alive or dead states of the cat.


title: Quantum Calculus-based Volterra LMS for Nonlinear Channel Estimation

abstract: A novel adaptive filtering method called $q$-Volterra least mean square ($q$-VLMS) is presented in this paper. The $q$-VLMS is a nonlinear extension of conventional LMS and it is based on Jackson's derivative also known as $q$-calculus. In Volterra LMS, due to large variance of input signal the convergence speed is very low. With proper manipulation we successfully improved the convergence performance of the Volterra LMS. The proposed algorithm is analyzed for the step-size bounds and results of analysis are verified through computer simulations for nonlinear channel estimation problem.


title: Free-Lunch Saliency via Attention in Atari Agents

abstract: We propose a new approach to visualize saliency maps for deep neural network models and apply it to deep reinforcement learning agents trained on Atari environments. Our method adds an attention module that we call FLS (Free Lunch Saliency) to the feature extractor from an established baseline (Mnih et al., 2015). This addition results in a trainable model that can produce saliency maps, i.e., visualizations of the importance of different parts of the input for the agent's current decision making. We show experimentally that a network with an FLS module exhibits performance similar to the baseline (i.e., it is "free", with no performance cost) and can be used as a drop-in replacement for reinforcement learning agents. We also design another feature extractor that scores slightly lower but provides higher-fidelity visualizations. In addition to attained scores, we report saliency metrics evaluated on the Atari-HEAD dataset of human gameplay.


title: UVIT Observations of UV-Bright Stars in four Galactic Globular Clusters

abstract: We have performed photometric analysis of four Galactic globular clusters (GGCs): NGC 4147, NGC 4590, NGC 5053 and NGC 7492 using far-UV and near-UV filters of the Ultraviolet Imaging Telescope (UVIT) on-board AstroSat. With the help of color-magnitude diagrams (CMDs), we have identified 150 blue horizontal branch stars (BHBs), and 40 blue straggler stars (BSS) in the four GGCs. We study the temperature and radial distribution of BHBs and BSS for the four GGCs.


title: Radial Drift and Concurrent Ablation of Boulder-Sized Objects

abstract: Context. The composition of a protoplanetary disk at a given location does not only depend on temperature and pressure but also on the time dependent transport of matter, such as radial drift of solid bodies, which could release water and other volatile species before disintegration or accretion onto a larger body with potentially considerable implications for the composition of planets. Aims. We perform a parameter study focused on the water depletion of different sized bodies able to cross the water snowline by gas induced radial drift. Methods. Either the analytical Hertz-Knudsen-Langmuir sublimation formula assuming equilibrium temperature within the body or a more involved, numerical model for the internal thermal evolution is coupled with an alpha-disk model. Different properties of the disk and the embedded body are explored. Results. Bodies with radii up to 100 meter drift faster towards the central star than the water snowline, hence, cross it. The region that can be reached before complete disintegration - and is therefore polluted with H2O ice - extends to 10 percent closer to the star than the snowline location. The extent of this polluted region could be multiple times larger in the presence of a dust mantle, which is, however, unlikely to form due to frequent collisions with smaller-than centimeter sized objects. Conclusions. Given a significant abundance of meter sized boulders in protoplanetary disks, the transport of water by radial drift of these bodies towards regions closer to the star than the snowline is not negligible and this flux of volatiles can be estimated for a given distribution of solid body sizes and compositions. A simple expression for surface sublimation is applicable for a homogeneous body consisting of only dust and water ice without a dust mantle.


title: Spin-nematic order induced superconductivity

abstract: We explore a spin-fermion model with fermion-spin-quadrupolar interaction. In a nematic phase, this interaction reduces to a four-fermion interaction that is a basis of superconductivity. When the coupling constant is positive the superconductivity is p-wave with spin-parallel paired fermions. When it is negative the superconductivity is p-wave and fermions are spin-antiparallel paired. For a system with zero chemical potential, even a very small coupling can bind fermions into bound state that leads to the superconductivity. When the chemical potential is non-zero the system possesses quantum critical transition from normal spin-nematic phase to phase where superconductivity coexists with spin-nematicity. The value of the quantum critical fermion-spin-nematicity coupling constant depends on the chemical potential.


title: Dynamical Gaseous Rings in Global Simulations of Protoplanetary Disk Formation

abstract: Global numerical simulations of protoplanetary disk formation and evolution were conducted in thin-disk limit, where the model included magnetically layered disk structure, a self-consistent treatment for the infall from cloud core as well as the smallest possible inner computational boundary. We compared the evolution of a layered disk with a fully magnetically active disk. We also studied how the evolution depends on the parameters of the layered disk model - the MRI triggering temperature and active layer thickness - as well as the mass of the prestellar cloud core. With the canonical values of parameters a dead zone formed within the inner $\approx$ 15 au region of the magnetically layered disk. The dead zone was not a uniform structure and long-lived, axisymmetric, gaseous rings ubiquitously formed within this region due to the action of viscous torques. The rings showed a remarkable contrast in the disk environment as compared to a fully magnetically active disk and were characterized by high surface density and low effective viscosity. Multiple gaseous rings could form simultaneously in the dead zone region which were highly dynamical and showed complex, time-dependent behavior such as inward migration, vortices, gravitational instability and large-scale spiral waves. An increase in MRI triggering temperature had only marginal effects, while changes in active layer thickness as well as the initial cloud core mass had significant effects on the structure and evolution of the inner disk. Dust with large fragmentation barrier could be trapped in the rings, which may play a key role in planet formation.


title: Time window to constrain the corner value of the global seismic-moment distribution

abstract: It is well accepted that, at the global scale, the Gutenberg-Richter (GR) law describing the distribution of earthquake magnitude or seismic moment has to be modified at the tail to properly account for the most extreme events. It is debated, though, how much additional time of earthquake recording will be necessary to properly constrain this tail. Using the global CMT catalog, we study how three modifications of the GR law that incorporate a corner-value parameter are compatible with the size of the largest observed earthquake in a given time window. Current data lead to a rather large range of parameter values (e.g., corner magnitude from 8.6 to 10.2 for the so-called tapered GR distribution). Updating this estimation in the future will strongly depend on the maximum magnitude observed, but, under reasonable assumptions, the range will be substantially reduced by the end of this century, contrary to claims in previous literature.


title: Machine learning and Kolmogorov analysis to reveal gravitational lenses

abstract: We present an automated approach to detect and extract information from the astronomical datasets on the shapes of such objects as galaxies, star clusters and, especially, elongated ones such as the gravitational lenses. First, the Kolmogorov stochasticity parameter is used to retrieve the sub-regions that worth further attention. Then we turn to image processing and machine learning Principal Component Analysis algorithm to retrieve the sought objects and reveal the information on their morphologies. We show the capability of our automated method to identify distinct objects, including of and to classify them based on the input parameters. A catalog of possible lensing objects is retrieved as an output of the software, then their inspection is performed for the candidates that survive the filters applied.


title: Ripser: efficient computation of Vietoris-Rips persistence barcodes

abstract: We present an algorithm for the computation of Vietoris-Rips persistence barcodes and describe its implementation in the software Ripser. The method relies on implicit representations of the coboundary operator and the filtration order of the simplices, avoiding the explicit construction and storage of the filtration coboundary matrix. Moreover, it makes use of apparent pairs, a simple but powerful method for constructing a discrete gradient field from a total order on the simplices of a simplicial complex, which is also of independent interest. Our implementation shows substantial improvements over previous software both in time and memory usage.


title: Towards model based control of the Vertical Gradient Freeze crystal growth process

abstract: In this contribution tracking control designs using output feedback are presented for a two-phase Stefan problem arising in the modeling of the Vertical Gradient Freeze process. The two-phase Stefan problem, consisting of two coupled free boundary problems, is a vital part of many crystal growth processes due to the temporally varying extent of the solid and liquid domains during growth. After discussing the special needs of the process, collocated as well as flatness-based state feedback designs are carried out. To render the setup complete, an observer design is performed, using a flatness-based approximation of the original distributed parameter system. The quality of the provided approximations as well as the performance of the open and closed loop control setups is analysed in several simulations.


title: Modularity and Projection of Bipartite Networks

abstract: This paper investigates community detection by modularity maximisation on bipartite networks. In particular we are interested in how the operation of projection, using one node set of the bipartite network to infer connections between nodes in the other set, interacts with community detection. We first define a notion of modularity appropriate for a projected bipartite network and outline an algorithm for maximising it in order to partition the network. Using both real and synthetic networks we compare the communities found by five different algorithms, where each algorithm maximises a different modularity function and sees different aspects of the bipartite structure. Based on these results we suggest a simple heuristic for finding communities in bipartite networks.


title: Statistical analysis of the azimuthal asymmetry in the $J/\psi$ leptoproduction in unpolarized $ep$ collisions

abstract: In this paper, we study the azimuthal asymmetry in the $J/\psi$ leptoproduction in unpolarized $ep$ collisions. There are two independent azimuthal asymmetry modulations, namely $\mathrm{cos}(\psi)$ and $\mathrm{cos}(2\psi)$, where $\psi$ is the azimuthal angle of the lepton scattering plane with respect to the hadron-interacting plane. We calculate the two modulations as functions of four kinematic variables, and find that they provide a very good laboratory to distinguish several models describing the heavy quarkonium production, including the color-singlet (CS) model, the nonrelativistic QCD (NRQCD) associated with the $^1S_0^{[8]}$ dominance picture, and the NRQCD in which the values of all the three color-octet (CO) long-distance matrix elements are of the same order. In order to make definite conclusions, we restrict our calculation in a specific kinematic region, where the CS and CO mechanisms can be distinguished by scrutinizing the values of the $\mathrm{cos}(\psi)$ modulation, while the $^1S_0^{[8]}$ dominance picture can be tested by measuring the values of the $\mathrm{cos}(2\psi)$ modulation. Calculating their values and carrying out a meticulous statistical analysis, we find that at an integrated luminosity $\mathcal{L}=1000\mathrm{pb}^{-1}$, the statistical uncertainties of the two quantities are small enough to tell the three models apart. When this experiment is implemented at the future $ep$ colliders such as the EIC, crucial information for the $J/\psi$ production mechanism might be discovered.


title: Elastic Properties of Symmetric Liquid-Liquid Interfaces

abstract: The mean ($\kappa$) and Gaussian ($\bar{\kappa}$) bending rigidities of liquid-liquid interfaces, of importance for shape fluctuations and topology of interfaces, respectively, are not yet established: even their signs are debated. Using the Scheutjens Fleer variant of the self-consistent field theory, we implemented a model for a symmetric L/L interface and obtained high precision (mean field) results in the grand canonical $(\mu, V, T)$-ensemble. We report positive values for both moduli when the system is close to critical where the rigidities show the same scaling behavior as the interfacial tension $\gamma$. At strong segregation, when the interfacial width becomes of the order of the segment size, $\bar{\kappa}$ turns negative. The length scale $\lambda \equiv \sqrt{\kappa/\gamma}$ is of order the segment size for all strengths of interaction; yet the $1/\sqrt{N}$ chain length correction reduces $\lambda$ significantly when the chain length $N$ is small.


title: V.I. Arnold's "pointwise" KAM Theorem

abstract: We review V.I. Arnold's 1963 celebrated paper \cite{ARV63} {\sl Proof of A.N. Kolmogorov's theorem on the conservation of conditionally periodic motions with a small variation in the Hamiltonian}, and prove that, optimizing Arnold's scheme, one can get "sharp" asymptotic quantitative conditions (as $\varepsilon\to 0$, $\varepsilon$ being the strength of the perturbation). All constants involved are explicitly computed.


title: Cross-Router Covert Channels

abstract: Many organizations protect secure networked devices from non-secure networked devices by assigning each class of devices to a different logical network. These two logical networks, commonly called the host network and the guest network, use the same router hardware, which is designed to isolate the two networks in software. In this work we show that logical network isolation based on host and guest networks can be overcome by the use of cross-router covert channels. Using specially-crafted network traffic, these channels make it possible to leak data between the host network and the guest network, and vice versa, through the use of the router as a shared medium. We performed a survey of routers representing multiple vendors and price points, and discovered that all of the routers we surveyed are vulnerable to at least one class of covert channel. Our attack can succeed even if the attacker has very limited permissions on the infected device, and even an iframe hosting malicious JavaScript code can be used for this purpose. We provide several metrics for the effectiveness of such channels, based on their pervasiveness, rate and covertness, and discuss possible ways of identifying and preventing these leakages.


title: Testing convexity of functions over finite domains

abstract: We establish new upper and lower bounds on the number of queries required to test convexity of functions over various discrete domains. 1. We provide a simplified version of the non-adaptive convexity tester on the line. We re-prove the upper bound $O(\frac{\log(\epsilon n)}{\epsilon})$ in the usual uniform model, and prove an $O(\frac{\log n}{\epsilon})$ upper bound in the distribution-free setting. 2. We show a tight lower bound of $\Omega(\frac{\log(\epsilon n)}{\epsilon})$ queries for testing convexity of functions $f: [n] \rightarrow \mathbb{R}$ on the line. This lower bound applies to both adaptive and non-adaptive algorithms, and matches the upper bound from item 1, showing that adaptivity does not help in this setting. 3. Moving to higher dimensions, we consider the case of a stripe $[3] \times [n]$. We construct an \emph{adaptive} tester for convexity of functions $f\colon [3] \times [n] \to \mathbb R$ with query complexity $O(\log^2 n)$. We also show that any \emph{non-adaptive} tester must use $\Omega(\sqrt{n})$ queries in this setting. Thus, adaptivity yields an exponential improvement for this problem. 4. For functions $f\colon [n]^d \to \mathbb R$ over domains of dimension $d \geq 2$, we show a non-adaptive query lower bound $\Omega((\frac{n}{d})^{\frac{d}{2}})$.


title: Brauer-Manin obstruction for Erd\H{o}s-Straus surfaces

abstract: We study the failure of the integral Hasse principle and strong approximation for the Erd\H{o}s-Straus conjecture using the Brauer-Manin obstruction.


title: Modelling atmospheric escape and MgII near-ultraviolet absorption of the highly irradiated hot Jupiter WASP-12b

abstract: We present two-dimensional multi-fluid numerical modelling of the upper atmosphere of the hot Jupiter WASP-12b. The model includes hydrogen chemistry, and self-consistently describes the expansion of the planetary upper atmosphere and mass loss due to intensive stellar irradiation, assuming a weakly magnetized planet. We simulate the planetary upper atmosphere and its interaction with the stellar wind (SW) with and without the inclusion of tidal force and consider different XUV irradiation conditions and SW parameters. With the inclusion of tidal force, even for a fast SW, the escaping planetary material forms two streams, propagating towards and away from the star. The atmospheric escape and related mass loss rate reaching the value of 10^12 gs^-1 appear to be mostly controlled by the stellar gravitational pull. We computed the column density and dynamics of MgII ions considering three different sets of SW parameters and XUV fluxes. The simulations enable to compute the absorption at the position of the Mg h line and to reproduce the times of ingress and egress. In case of a slow SW and without accounting for tidal force, the high orbital velocity leads to the formation of a shock approximately in the direction of the planetary orbital motion. In this case, mass loss is proportional to the stellar XUV flux. At the same time, ignoring of tidal effects for WASP-12b is a strong simplification, so the scenario with a shock, altogether is an unrealistic one.


title: Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals

abstract: A recently developed finite element approach for fully numerical atomic structure calculations [S. Lehtola, Int. J. Quantum Chem. 119, e25945 (2019)] is extended to the description of atoms with spherically symmetric densities via fractionally occupied orbitals. Specialized versions of Hartree-Fock as well as local density and generalized gradient approximation density functionals are developed, allowing extremely rapid calculations at the basis set limit on the ground and low-lying excited states even for heavy atoms. The implementation of range-separation based on the Yukawa or complementary error function (erfc) kernels is also described, allowing complete basis set benchmarks of modern range-separated hybrid functionals with either integer or fractional occupation numbers. Finally, computation of atomic effective potentials at the local density or generalized gradient approximation levels for the superposition of atomic potentials (SAP) approach [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)] that has been shown to be a simple and efficient way to initialize electronic structure calculations is described. The present numerical approach is shown to afford beyond microhartree accuracy with a small number of numerical basis functions, and to reproduce literature results for the ground states of atoms and their cations for $1 \leq Z \leq 86 $. Our results indicate that the literature values deviate by up to 10 {\mu}Eh from the complete basis set limit. The numerical scheme for the erfc kernel is shown to work by comparison to results from large Gaussian basis set calculations from the literature. Spin-restricted ground states are reported for Hartree-Fock and Hartree-Fock-Slater calculations with fractional occupations for $1 \leq Z \leq 118$.


title: Escaping orbits are also rare in the almost periodic Fermi-Ulam ping-pong

abstract: We study the one-dimensional Fermi-Ulam ping-pong problem with a Bohr almost periodic forcing function and show that the set of initial condition leading to escaping orbits typically has Lebesgue measure zero.


title: A Constraint Model for the Tree Decomposition of a Graph

abstract: We present a constraint model for the problem of producing a tree decomposition of a graph. The inputs to the model are a simple graph G, the number of nodes in the desired tree decomposition and the maximum cardinality of each node in that decomposition. Via a sequence of decision problems, the model allows us to find the tree width of a graph whilst delivering a tree decomposition of that width, i.e. a witness.


title: Global performance and capabilities of the instruments on the 3.6-m Devasthal Optical Telescope

abstract: The recently commissioned 3.6-m Devasthal optical telescope has been used for various tests and science observations using three main instruments, namely, a charge-coupled device camera, a near-infrared camera, and an optical imager-cum-spectrograph. The published results from these instruments assert that the performance of the telescope at the Devasthal site is at par with the expectations. These back-end instruments open up vast opportunities for high-sensitivity observations of the celestial sky with the telescope. This paper provides a summary of the existing back-end instruments and attempts to highlight the importance of the Devasthal optical telescope in synergy with other telescopes operating at different wavelengths.


title: Environmental Disorder Regulation of Invasion and Genetic Loss

abstract: Many physical and natural systems, including the population of species, evolve in habitats with spatial stochastic variations of the individuals' motility. We study here the effect of those fluctuations on invasion and genetic loss. A Langevin equation for the \textit{position} and \textit{border} of the invasion front is obtained. A striking result is that small/large fluctuations of diffusivity suppress/intensify genetic loss. Our findings reveal the potential role of environmental fluctuations as a regulating factor for genetic loss and provide a simple explanation for the regional differences in the intensity of genetic drift observed during the final stages of human evolution and in tumor mutational landscapes.


title: Implication of Kinetic Alfven Waves to Magnetic Field Turbulence Spectra: Earth's Magnetosheath

abstract: In the present paper, we investigate the power-law behaviour of the magnetic field spectra in the Earths magnetosheath region using Cluster spacecraft data under solar minimum condition. The power spectral density of the magnetic field data and spectral slopes at various frequencies are analysed. Propagation angle and compressibility are used to test the nature of turbulent fluctuations. The magnetic field spectra have the spectral slopes between -1.5 to 0 down to spatial scales of 20 ion gyroradius and show clear evidence of a transition to steeper spectra for small scales with a second power-law, having slopes between -2.6 to -1.8. At low frequencies, f_sc<0.3f_ci(where f_ci is ion gyro-frequency), propagation angle approximately 90 degrees to the mean magnetic field, B_0, and compressibility shows a broad distribution, 0.1 < R > 0.9. On the other hand at f_sc>10f_ci, the propagation angle exhibits a broad range between 30-90 degree while 'R' has a small variation: 0.2 < R > 0.5. We conjecture that at high frequencies, the perpendicularly propagating Alfven waves could partly explain the statistical analysis of spectra. To support our prediction of kinetic Alfven wave-dominated spectral slope behaviour at high frequency, we also present a theoretical model and simulate the magnetic field turbulence spectra due to the nonlinear evolution of kinetic Alfven waves. The present study also shows the analogy between the observational and simulated spectra.


title: de Sitter duality and logarithmic decay of dark energy

abstract: We investigate infrared dynamics of four-dimensional Einstein gravity in de Sitter space. We set up a general framework to investigate dynamical scaling relations in quantum/classical gravitational theories. The conformal mode dependence of Einstein gravity is renormalized to the extent that general covariance is not manifest. We point out that the introduction of an inflaton is necessary as a counterterm. We observe and postulate a duality between quantum effects in Einstein gravity and classical evolutions in an inflation (or quintessence) model. The effective action of Einstein gravity can be constructed as an inflation model with manifest general covariance. We show that $g=G_N H^2/\pi$: the only dimensionless coupling of the Hubble parameter $H^2$ and the Newton's coupling $G_N$ in Einstein gravity is screened by the infrared fluctuations of the conformal mode. We evaluate the one-loop $\beta$ function of $g$ with respect to the cosmic time $\log Ht$ as $\beta(g)=-(1/2)g^2$, i.e., $g$ is asymptotically free toward the future. The exact $\beta$ function with the backreaction of $g$ reveals the existence of the ultraviolet fixed point. It indicates that the de Sitter expansion stared at the Planck scale with a minimal entropy $S=2$. We have identified the de Sitter entropy $1/g$ with the von Neumann entropy of the conformal zero mode. The former evolves according to the screening of $g$ and the Gibbons-Hawking formula. The latter is found to increase by diffusion in the stochastic process at the horizon in a consistent way. Our Universe is located very close to the fixed point $g=0$ with a large entropy. We discuss possible physical implications of our results such as logarithmic decay of dark energy.


title: Uniform bounds on harmonic Beltrami differentials and Weil-Petersson curvatures

abstract: In this article we show that for every finite area hyperbolic surface $X$ of type $(g,n)$ and any harmonic Beltrami differential $\mu$ on $X$, then the magnitude of $\mu$ at any point of small injectivity radius is uniform bounded from above by the ratio of the Weil-Petersson norm of $\mu$ over the square root of the systole of $X$ up to a uniform positive constant multiplication. We apply the uniform bound above to show that the Weil-Petersson Ricci curvature, restricted at any hyperbolic surface of short systole in the moduli space, is uniformly bounded from below by the negative reciprocal of the systole up to a uniform positive constant multiplication. As an application, we show that the average total Weil-Petersson scalar curvature over the moduli space is uniformly comparable to $-g$ as the genus $g$ goes to infinity.


title: The existence of a nontrivial weak solution to a double critical problem involving fractional Laplacian in ${\R}^n$ with a Hardy term

abstract: In this paper, we consider the existence of nontrivial weak solutions to a double critical problem involving fractional Laplacian with a Hardy term: \begin{equation} \label{eq0.1} (-\Delta)^{s}u-{\gamma} {\frac{u}{|x|^{2s}}}= {\frac{{|u|}^{ {2^{*}_{s}}(\beta)-2}u}{|x|^{\beta}}}+ \big [ I_{\mu}* F_{\alpha}(\cdot,u) \big](x)f_{\alpha}(x,u), \ \ u \in {\dot{H}}^s(\R^{n}) \end{equation} where $s \in(0,1)$, $0\leq \alpha,\beta<2s<n$, $\mu \in (0,n)$, $\gamma<\gamma_{H}$, $I_{\mu}(x)=|x|^{-\mu}$, $F_{\alpha}(x,u)=\frac{ {|u(x)|}^{ {2^{\#}_{\mu} }(\alpha)} }{ {|x|}^{ {\delta_{\mu} (\alpha)} } }$, $f_{\alpha}(x,u)=\frac{ {|u(x)|}^{{ 2^{\#}_{\mu} }(\alpha)-2}u(x) }{ {|x|}^{ {\delta_{\mu} (\alpha)} } }$, $2^{\#}_{\mu} (\alpha)=(1-\frac{\mu}{2n})\cdot 2^{*}_{s} (\alpha)$, $\delta_{\mu} (\alpha)=(1-\frac{\mu}{2n})\alpha$, ${2^{*}_{s}}(\alpha)=\frac{2(n-\alpha)}{n-2s}$ and $\gamma_{H}=4^s\frac{\Gamma^2(\frac{n+2s}{4})} {\Gamma^2(\frac{n-2s}{4})}$. We show that problem (\ref{eq0.1}) admits at least a weak solution under some conditions. To prove the main result, we develop some useful tools based on a weighted Morrey space. To be precise, we discover the embeddings \begin{equation} \label{eq0.2} {\dot{H}}^s(\R^{n}) \hookrightarrow {L}^{2^*_{s}(\alpha)}(\R^{n},|y|^{-\alpha}) \hookrightarrow L^{p,\frac{n-2s}{2}p+pr}(\R^{n},|y|^{-pr}) \end{equation} where $s \in (0,1)$, $0<\alpha<2s<n$, $p\in[1,2^*_{s}(\alpha))$, $r=\frac{\alpha}{ 2^*_{s}(\alpha) }$; We also establish an improved Sobolev inequality. By using mountain pass lemma along with an improved Sobolev inequality, we obtain a nontrivial weak solution to problem (\ref{eq0.1}) in a direct way. It is worth while to point out that the improved Sobolev inequality could be applied to simplify the proof of the main results in \cite{NGSS} and \cite{RFPP}.


title: Deriving phase field crystal theory from dynamical density functional theory: consequences of the approximations

abstract: Phase field crystal (PFC) theory, extensively used for modelling the structure of solids, can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Standard derivations neglect a term of form $\nabla\cdot[n\nabla L n]$, where $n$ is the scaled density profile and $L$ is a linear operator. We show that this term makes a significant contribution to the stability of the crystal, and dropping this term from the theory forces another approximation, that of replacing the logarithmic term from the ideal gas contribution to the free energy with its truncated Taylor expansion, to yield a polynomial in $n$. However, the consequences of doing this are the presence of an additional spinodal in the phase diagram, so the liquid is predicted first to freeze and then to melt again as the density is increased; and other periodic structures are erroneously predicted to be thermodynamic equilibria. A second approximation is to replace $L$ by a gradient expansion. This leads to the possibility of solutions failing to exist above a certain value of the average density. We illustrate these conclusions with a simple model two-dimensional fluid. The consequences of the PFC approximations are that the phase diagram is both qualitatively incorrect, in that it has a stripe phase, and quantitatively incorrect (by orders of magnitude) regarding the properties of the crystal. Thus, although PFC models are successful as phenomenological models of crystallisation, we find it impossible to derive the PFC model as an accurate approximation to DDFT, without introducing spurious artefacts. However, making a simple one-mode approximation for the logarithm of the density distribution is surprisingly accurate, which gives a tantalising hint that accurate PFC-type theories may instead be derived for the field $\log(\rho(x))$, rather than for the density profile itself.


title: Unraveling the hidden organisation of urban systems and their mobility flows

abstract: Increasing evidence suggests that cities are complex systems, with structural and dynamical features responsible for a broad spectrum of emerging phenomena. Here we use a unique data set of human flows and couple it with information on the underlying street network to study, simultaneously, the structural and functional organisation of 10 world megacities. We quantify the efficiency of flow exchange between areas of a city in terms of integration and segregation using well defined measures. Results reveal unexpected complex patterns that shed new light on urban organisation. Large cities tend to be more segregated and less integrated, while their overall topological organisation resembles that of small world networks. At the same time, the heterogeneity of flows distribution might act as a catalyst for further integrating a city. Our analysis unravels how human behaviour influences, and is influenced by, the urban environment, suggesting quantitative indicators to control integration and segregation of human flows that can be used, among others, for restriction policies to adopt during emergencies and, as an interesting byproduct, allows us to characterise functional (dis)similarities of different metropolitan areas, countries, and cultures.


title: A Note on the Possibility of Self-Reference in Mathematics

abstract: In this paper we propose an interpretation for self-referential propositions in a "meta-model" N* of ZF. This meta-model N* is considered as an informal model of arithmetic that mathematicians often use when working with number theory. Specifically, we assume that within this meta-model, the axiom system ZF is applied, interpretations for sentences can be offered, and natural language can be used. We show that under the proposed interpretation, some types of self-referential propositions that are considered legitimate in mathematics turn N* into an inconsistent model, and examine the connection of this result to a certain interpretation of godel's first incompleteness theorem. Some general problems which follow from the above discussion are then addressed.


title: Migrant mobility flows characterized with digital data

abstract: Monitoring migration flows is crucial to respond to humanitarian crisis and to design efficient policies. This information usually comes from surveys and border controls, but timely accessibility and methodological concerns reduce its usefulness. Here, we propose a method to detect migration flows worldwide using geolocated Twitter data. We focus on the migration crisis in Venezuela and show that the calculated flows are consistent with official statistics at country level. Our method is versatile and far-reaching, as it can be used to study different features of migration as preferred routes, settlement areas, mobility through several countries, spatial integration in cities, etc. It provides finer geographical and temporal resolutions, allowing the exploration of issues not contemplated in official records. It is our hope that these new sources of information can complement official ones, helping authorities and humanitarian organizations to better assess when and where to intervene on the ground.


title: Performance Study of Strongly Coupled Magnetic Resonance

abstract: Strongly Coupled Magnetic Resonance (SCMR) uses electromagnetic resonance in order to efficiently transfer power wirelessly over mid-range distances. Since the energy exchange capability of resonant objects higher than non-resonant objects, strongly coupled systems are able to achieve more efficient energy transfer than other wireless power transfer systems. The paper presents detailed experimental and simulated analysis of the performance of the SCMR system. A prototype of the SCMR system was implemented and experiments were conducted to analyze the performance of the system. Finally, the resonant frequency of the system was experimentally verified and the factors influencing the wireless power transfer were also studied


title: Probing the era of reheating for reconstructed inflationary potential in the RS II braneworld

abstract: We analyse the epoch of reheating after an inflationary phase in the Randal Sundrum(RS) Type-$\rm II$ braneworld, where we did not consider any particular model of inflation, but rather reconstructed the inflationary potential solving the flow equations using Monte Carlo (MC) approach. It is shown numerically that a potential conceived through the MC reconstruction technique can be represented by an effective potential as a function of the number of e-foldings($N$). Then, the epoch of reheating is studied for this reconstructed potential. The relation between the reheating temperature ($T_{\rm reh}$) and the 5-dimensional Planck mass($M_5$) is established. Moreover, it is argued that there is a stringent bound on the critical reheating temperature that also translates to a tight bound on $M_5$ .


title: Prediction of the spin triplet two-electron quantum dots in Si: towards controlled quantum simulations of magnetic systems

abstract: Ground state of two-electron quantum dots in single-valley materials like GaAs is always a spin singlet regardless of what the potential and interactions are. This statement cannot be generalized to the multi-valley materials like $n$-doped Si. Here we calculate numerically the spectrum of a two-electron Si quantum dot and show that the dot with the lateral size of several nm can have the spin triplet ground state which is impossible in the single-valley materials. Predicted singlet-triplet level crossing in two-electron Si quantum dots can potentially establish the platform for quantum simulation of magnetic many body systems based on quantum dots. We suggest several examples of such systems that open a way to controlled quantum simulations within the condensed matter setting.


title: A nano-carbon route to rare earth free permanent magnetism

abstract: High coercivity magnets are an important resource for renewable energy, electric vehicles and memory technologies. Most hard magnetic materials incorporate rare-earths such as neodymium and samarium, but the concerns about the environmental impact and supply stability of these materials is prompting research into alternatives. Here, we present a hybrid bilayer of cobalt and the nano-carbon molecule C60 which exhibits significantly enhanced coercivity with minimal reduction in magnetisation. We demonstrate how this anisotropy enhancing effect cannot be described by existing models of molecule-metal magnetic interfaces. We outline a new form of magnetic anisotropy, arising from asymmetric magneto-electric coupling in the metal-molecule interface. Because this phenomenon arises from pi-d hybrid orbitals, we propose calling this effect pi-anisotropy. While the critical temperature of this effect is currently limited by the rotational degree of freedom of the chosen molecule, C60, we describe how surface functionalisation would allow for the design of room-temperature, carbon based hard magnetic films.


title: Quantile-Frequency Analysis and Spectral Divergence Metrics for Diagnostic Checks of Time Series With Nonlinear Dynamics

abstract: Nonlinear dynamic volatility has been observed in many financial time series. The recently proposed quantile periodogram offers an alternative way to examine this phenomena in the frequency domain. The quantile periodogram is constructed from trigonometric quantile regression of time series data at different frequencies and quantile levels. It is a useful tool for quantile-frequency analysis (QFA) of nonlinear serial dependence. This paper introduces a number of spectral divergence metrics based on the quantile periodogram for diagnostic checks of financial time series models and model-based discriminant analysis. The parametric bootstrapping technique is employed to compute the $p$-values of the metrics. The usefulness of the proposed method is demonstrated empirically by a case study using the daily log returns of the S\&P 500 index over three periods of time together with their GARCH-type models. The results show that the QFA method is able to provide additional insights into the goodness of fit of these financial time series models that may have been missed by conventional tests. The results also show that the QFA method offers a more informative way of discriminant analysis for detecting regime changes in time series.


title: Reconfigurable sub-micron spin-wave majority gate with electrical transducers

abstract: Spin waves are excitations in ferromagnetic media that have been proposed as information carriers in hybrid spintronic devices with much lower operation power than conventional charge-based electronics. Their wave nature can be exploited in majority gates by using interference for computation. However, a scalable spin-wave majority gate that can be co-integrated alongside conventional electronics is still lacking. Here, we demonstrate a sub-micron inline spin-wave majority gate with fan-out. Time-resolved imaging of the magnetization dynamics by scanning transmission x-ray microscopy illustrates the device operation. All-electrical spin-wave spectroscopy further demonstrates majority gates with sub-micron dimensions, reconfigurable input and output ports, and frequency-division multiplexing. Challenges for hybrid spintronic computing systems based on spin-wave majority gates are discussed.


title: A Repairable System Supported by Two Spare Units and Serviced by Two Types of Repairers

abstract: We study a one-unit repairable system, supported by two identical spare units on cold standby, and serviced by two types of repairers. The model applies, for instance, to ANSI (American National Standard Institute) centrifugal pumps in a chemical plant. The failed unit undergoes repair either by an in-house repairer within a random or deterministic patience time, or else by a visiting expert repairer. The expert repairs one or all failed units before leaving, and does so faster but at a higher cost rate than the regular repairer. Four models arise depending on the number of repairs done by the expert and the nature of the patience time. We compare these models based on the limiting availability $A_{\infty}$, and the limiting profit per unit time $\omega$, using semi-Markov processes, when all distributions are exponential. As anticipated, to maximize $A_{\infty}$, the expert should repair all failed units. To maximize $\omega$, a suitably chosen deterministic patience time is better than a random patience time. Furthermore, given all cost parameters, we determine the optimum number of repairs the expert should complete, and the optimum patience time given to the regular repairer in order to maximize $\omega$.


title: Automated Corrosion Detection Using Crowd Sourced Training for Deep Learning

abstract: The automated detection of corrosion from images (i.e., photographs) or video (i.e., drone footage) presents significant advantages in terms of corrosion monitoring. Such advantages include access to remote locations, mitigation of risk to inspectors, cost savings and monitoring speed. The automated detection of corrosion requires deep learning to approach human level artificial intelligence (A.I.). The training of a deep learning model requires intensive image labelling, and in order to generate a large database of labelled images, crowd sourced labelling via a dedicated website was sought. The website (corrosiondetector.com) permits any user to label images, with such labelling then contributing to the training of a cloud based A.I. model - with such a cloud-based model then capable of assessing any fresh (or uploaded) image for the presence of corrosion. In other words, the website includes both the crowd sourced training process, but also the end use of the evolving model. Herein, the results and findings from the website (corrosiondetector.com) over the period of approximately one month, are reported.


title: Actions of monoidal categories and representations of Cartan type Lie algebras

abstract: Using crossed homomorphisms, we show that the category of weak representations (resp. admissible representations) of Lie-Rinehart algebras (resp. Leibniz pairs) is a left module category over the monoidal category of representations of Lie algebras. In particular, the corresponding bifunctor of monoidal categories is established to give new weak representations (resp. admissible representations) of Lie-Rinehart algebras (resp. Leibniz pairs). This generalizes and unifies various existing constructions of representations of many Lie algebras by using this new bifunctor. We construct some crossed homomorphisms in different situations and use our actions of monoidal categories to recover some known constructions of representations of various Lie algebras, also to obtain new representations for generalized Witt algebras and their Lie subalgebras. The cohomology theory of crossed homomorphisms between Lie algebras is introduced and used to study linear deformations of crossed homomorphisms.

categories: math.RT math-ph math.DG math.MP math.QA math.RA


title: About quasiperiodic tilings, possessing five-fold symmetry

abstract: A group-theoretical approach to the construction of quasiperiodic tilings of a Euclidean plane, possessing five-fold symmetry, is applied. Of the infinitely many of variants of quasiperiodic partitions of the plane, possessing the dihedral group of symmetry D5, special attention is paid to those that, can be obtained by one universal set, with consists five different tiles. Geometric characteristics of tiles from this set are determined. It is shown, by this set of tiles it is possible to carry out topologically different tilings of the plane, possessing rotating symmetries of both the fifth and tenth orders.


title: Tweets Can Tell: Activity Recognition using Hybrid Long Short-Term Memory Model

abstract: This paper presents techniques to detect the "offline" activity a person is engaged in when she is tweeting (such as dining, shopping or entertainment), in order to create a dynamic profile of the user, for uses such as better targeting of advertisements. To this end, we propose a hybrid LSTM model for rich contextual learning, along with studies on the effects of applying and combining multiple LSTM based methods with different contextual features. The hybrid model is shown to outperform a set of baselines and state-of-the-art methods. Finally, this paper presents an orthogonal validation with a real-case application. Our model generates an offline activity analysis for the followers of several well-known accounts, which is quite representative of the expected characteristics of these accounts.


title: Efficient Estimation by Fully Modified GLS with an Application to the Environmental Kuznets Curve

abstract: This paper develops the asymptotic theory of a Fully Modified Generalized Least Squares estimator for multivariate cointegrating polynomial regressions. Such regressions allow for deterministic trends, stochastic trends and integer powers of stochastic trends to enter the cointegrating relations. Our fully modified estimator incorporates: (1) the direct estimation of the inverse autocovariance matrix of the multidimensional errors, and (2) second order bias corrections. The resulting estimator has the intuitive interpretation of applying a weighted least squares objective function to filtered data series. Moreover, the required second order bias corrections are convenient byproducts of our approach and lead to standard asymptotic inference. We also study several multivariate KPSS-type of tests for the null of cointegration. A comprehensive simulation study shows good performance of the FM-GLS estimator and the related tests. As a practical illustration, we reinvestigate the Environmental Kuznets Curve (EKC) hypothesis for six early industrialized countries as in Wagner et al. (2020).


title: Security measurement of a medical communication scheme based on chaos and DNA coding

abstract: To encrypt sensitive information existing in a color DICOM images, a medical privacy protection scheme (called as MPPS) based on chaos and DNA coding was proposed by using two coupled chaotic systems to produce cryptographic primitives. Relying on some empirical analyses and experimental results, the designers of MPPS claimed that it can withstand a chosen-plaintext attack and some other classic attacking models. However, this statement is groundless. In this paper, we investigate the essential properties of MPPS and DNA coding, and we then propose an efficient chosen-plaintext attack to disclose its equivalent secret-key. The attack only needs $\lceil \log_{256}(3\cdot M\cdot N)\rceil+4$ pair of chosen plain-images and the corresponding cipher-images, where $M \times N$ and ``3" are the size of the RGB color image and the number of color channels, respectively. In addition, the other claimed superiorities are questioned from the perspective of modern cryptography. Both theoretical and experimental results are presented to support the efficiency of the proposed attack and the other reported security faults. The proposed cryptanalysis results will promote the proper application of DNA encoding to protect multimedia privacy data, especially that in a DICOM image.


title: Attractive conical surfaces create infinitely many bound states

abstract: In this paper we study spectral properties of a three-dimensional Schr\"odinger operator $-\Delta+V$ with a potential $V$ given, modulo rapidly decaying terms, by a function of the distance of $x \in \mathbb{R}^3$ to an infinite conical hypersurface with a smooth cross-section. As a main result we show that there are infinitely many discrete eigenvalues accumulating at the bottom of the essential spectrum which itself is identified as the ground-state energy of a certain one-dimensional operator. Most importantly, based on a result of Kirsch and Simon we are able to establish the asymptotic behavior of the eigenvalue counting function using an explicit spectral-geometric quantity associated with the cross-section. This shows a universal character of some previous results on conical layers and $\delta$-potentials created by conical surfaces.


title: Key Features of the Coupled Hand-operated Balanced Manipulator (HOBM) and Lightweight Robot (LWR)

abstract: The paper deals with coupled systems including hand-operated balanced manipulators and lightweight robots. The aim of such a cooperation is to displace heavy payloads with less powerful robots. In other term, in the coupled system for handling of heavy payloads by a HOBM an operator is replaced by a LWR. The advantages of the coupled HOBM and LWR are disclosed and the optimal design of the cooperative workspace is discussed. Behavior of the coupled system in a static mode when the velocities the HOBM are limited does not present any special problems. In this mode, the inertial forces are significantly lower than the gravitational one. The payload is completely balanced by the HOBM and the LWR assumes the prescribed displacements with low load. However, in a dynamic mode, the HOBM with massive links creates additional loads on the LWR, which can be significant. The present study considers a method for determination of inertia effects of the HOBM on the LWR. The given numerical simulations show the significant increasing of the input torques due to the inertia forces of the HOBM. Behavior of the HOBM with cable lift and the LWR is also examined.


title: Embedding-based Silhouette Community Detection

abstract: Mining complex data in the form of networks is of increasing interest in many scientific disciplines. Network communities correspond to densely connected subnetworks, and often represent key functional parts of real-world systems. In this work, we propose Silhouette Community Detection (SCD), an approach for detecting communities, based on clustering of network node embeddings, i.e. real valued representations of nodes derived from their neighborhoods. We investigate the performance of the proposed SCD approach on 234 synthetic networks, as well as on a real-life social network. Even though SCD is not based on any form of modularity optimization, it performs comparably or better than state-of-the-art community detection algorithms, such as the InfoMap and Louvain algorithms. Further, we demonstrate how SCD's outputs can be used along with domain ontologies in semantic subgroup discovery, yielding human-understandable explanations of communities detected in a real-life protein interaction network. Being embedding-based, SCD is widely applicable and can be tested out-of-the-box as part of many existing network learning and exploration pipelines.


title: Bononia, the Roman Bologna: Archaeoastronomy and Chronology

abstract: In an article written by Giulio Magli on the orientation of the Roman towns, Bononia, the Roman Bologna, is given as a specific example to support Magli's thesis on the existence of preferred solar alignments of the urban layout. Assuming that the Roman towns had been oriented to the sunrise on a given day of the year, Magli suggested possible preferred alignments according to Roman festivals, in particular the festival of Terminalia. Of Bononia, we know the date of foundation as Roman colony in 189 BC, given by Livy. We will show that, according to Roman chronology and Republican calendar, it is impossible that Bononia had been oriented to the sunrise on the day given by Livy. The discrepancy is huge. Moreover, the direction of the decumanus cannot match the dates of Terminalia for 189 BC. However, if we consider that the layout that we see today is that of a recolonization under Octavian, we can have a certain agreement between the direction of the decumanus and the sunrise on the day of Terminalia (within three days), and a perfect agreement with the day of the festival of Armilustrium. In the proposed discussion, we will stress in particular the problem of the discrepancy between the historical dates of Roman chronology and the Julian dates, those that we can find according to an astronomical analysis. This problem is general and must be properly considered in any archaeoastronomical analysis of Roman towns.


title: Next-Generation High-Resolution Vector-Borne Disease Risk Assessment

abstract: Vector-borne diseases cause more than 1 million deaths annually. Estimates of epidemic risk at high spatial resolutions can enable effective public health interventions. Our goal is to identify the risk of importation of such diseases into vulnerable cities at the granularity of neighborhoods. Conventional models cannot achieve such spatial resolution, especially in real-time. Besides, they lack real-time data on demographic heterogeneity, which is vital for accurate risk estimation. Social media, such as Twitter, promise data from which demographic and spatial information could be inferred in real-time. On the other hand, such data can be noisy and inaccurate. Our novel approach leverages Twitter data, using machine learning techniques at multiple spatial scales to overcome its limitations, to deliver results at the desired resolution. We validate our method against the Zika outbreak in Florida in 2016. Our main contribution lies in proposing a novel approach that uses machine learning on social media data to identify the risk of vector-borne disease importation at a sufficiently fine spatial resolution to permit effective intervention. This will lead to a new generation of epidemic risk assessment models, promising to transform public health by identifying specific locations for targeted intervention.


title: On The Inverse Relaxation Approach To Supercapacitors Characterization

abstract: A novel inverse relaxation technique for supercapacitor characterization is developed, modeled numerically, and experimentally tested on a number of commercial supercapacitors. It consists in shorting a supercapacitor for a short time $\tau$, then switching to the open circuit regime and measuring an initial rebound and long-time relaxation. The results obtained are: the ratio of "easy" and "hard" to access capacitance and the dependence $C(\tau)$, that determines what the capacitance the system responds at time-scale $\tau$; it can be viewed as an alternative to used by some manufacturers approach to characterize a supercapacitor by fixed capacitance and time-scale dependent internal resistance. Among the advantages of proposed technique is that it does not require a source of fixed current, what simplifies the setup and allows a high discharge current regime. The approach can be used as a replacement of low-frequency impedance measurements and the ones of IEC 62391 type, it can be effectively applied to characterization of supercapacitors and other relaxation type systems with porous internal structure. The technique can be completely automated by a microcontroller to measure, analyze, and output the results.


title: First Law of Entanglement Entropy in Flat-Space Holography

abstract: According to flat/Bondi-Metzner-Sachs invariant field theories (BMSFT) correspondence, asymptotically flat spacetimes in $(d+1)$-dimensions are dual to $d$-dimensional BMSFTs. In this duality, similar to the Ryu-Takayanagi proposal in the AdS/CFT correspondence, the entanglement entropy of subsystems in the field theory side is given by the area of some particular surfaces in the gravity side. In this paper we find the holographic counterpart of the first law of entanglement entropy (FLEE) in a two-dimensional BMSFT. We show that FLEE for the BMSFT perturbed states which are descried by three-dimensional flat-space cosmology, corresponds to the integral of a particular one-form on a closed curve. This curve consists of BMSFT interval and also null and spacelike geodesics in the bulk gravitational theory. Exterior derivative of this form is zero when it is calculated for the flat-space cosmology. However, for a generic perturbation of three-dimensional global Minkowski spacetime, the exterior derivative of one-form yields Einstein equation. This is the first step for constructing bulk geometry by using FLEE in the flat/BMSFT correspondence.


title: Etat de l'art des syst\`emes robotis\'es en vue d'une application pour la chirurgie otologique

abstract: This article deals with a patents state-of-the-art linked to our robotic system which holds an endoscope. In a first part, we analyze how the endoscope operates as well as we discover the environment where we wish to use it: the middle ear. In a second step, we establish a state-of-the-art on patents existing in the robotic and surgical fields. The approach is based on a specific search methodology using keywords and classifications. Thanks to all the patents found, we gather them according their country and date of publication, and we are able to analyze the results. This study brings to light technical solutions invented in recent years and allows us to find an innovative research axis.


title: The Elliptic Kashiwara-Vergne Lie algebra in low weights

abstract: In this paper, we study the elliptic Kashiwara-Vergne Lie Algebra $\mathfrak{krv}$, which is a certain Lie subalgebra of the Lie algebra of derivations of the free Lie algebra in two generators. It has a natural bigrading, such that the Lie bracket is of bidegree $(-1,-1)$. After recalling the graphical interpretation of this Lie algebra, we examine low degree elements of $\mathfrak{krv}$. More precisely, we find that $\mathfrak{krv}^{(2,j)}$ is one-dimensional for even $j$ and zero $j$ odd. We also compute $\operatorname{dim}(\mathfrak{krv})^{(3,m)} = \lfloor\frac{m-1}{2}\rfloor - \lfloor\frac{m-1}{3}\rfloor$. In particular, we show that in those degrees there are no odd elements and also confirm Enriquez' conjecture in those degrees.


title: Recoiling supermassive black holes in analytical and numerical galaxy potential

abstract: We follow trajectories of recoiling supermassive black holes (SMBHs) in analytical and numerical models of galaxy merger remnants with masses of $10^{11} \rm{M_{sun}}$ and $10^{12} \rm{M_{sun}}$. We construct various merger remnant galaxies in order to investigate how the central SMBH mass and the mass ratio of progenitor galaxies influence escape velocities of recoiling SMBHs. Our results show that static analytical models of major merger remnant galaxies overestimate the SMBHs escape velocities. During major mergers violent relaxation leads to the decrease of galaxy mass and lower potential at large remnant radii. This process is not depicted in static analytical potential but clearly seen in our numerical models. Thus, the evolving numerical model is a more realistic description of dynamical processes in galaxies with merging SMBHs. We find that SMBH escape velocities in numerical major merger remnant galaxies can be up to 25 per cent lower compared to those in analytical models. Consequently, SMBHs in numerical models generally reach greater galactocentric distances and spend more time on bound orbits outside of the galactic nuclei. Thus, numerical models predict a greater number of spatially-offset active galactic nuclei (AGNs).


title: Grasp Type Estimation for Myoelectric Prostheses using Point Cloud Feature Learning

abstract: Prosthetic hands can help people with limb difference to return to their life routines. Commercial prostheses, however have several limitations in providing an acceptable dexterity. We approach these limitations by augmenting the prosthetic hands with an off-the-shelf depth sensor to enable the prosthesis to see the object's depth, record a single view (2.5-D) snapshot, and estimate an appropriate grasp type; using a deep network architecture based on 3D point clouds called PointNet. The human can act as the supervisor throughout the procedure by accepting or refusing the suggested grasp type. We achieved the grasp classification accuracy of up to 88%. Contrary to the case of the RGB data, the depth data provides all the necessary object shape information, which is required for grasp recognition. The PointNet not only enables using 3-D data in practice, but it also prevents excessive computations. Augmentation of the prosthetic hands with such a semi-autonomous system can lead to better differentiation of grasp types, less burden on user, and better performance.


title: Analyse du besoin en assistance robotique dans la chirurgie de l'oreille

abstract: Otologic surgery has some specificities compared to others surgeries. The anatomic working space is small, with various anatomical structures to preserve, like ossicles or facial nerve. This requires the use of microscope or endoscope. The microscope let the surgeon use both hands, but allows only direct vision. The endoscope leaves only one hand to the surgeon to use his tools, but provides a "fish-eye" vision. The rise of endoscopy these past few years has led to the development of numerous devices for the surgeon: the Robotol, first otological robot designed to performed some movements and hold an endoscope, or the Endofix Exo. Both devices need the hand of the surgeon to be moved. No robotic device allows the endoscope to be directed autonomously while the surgeon keeps both hands free to work, just like when he is working with a microscope. The objective of our work is to define the specific needs of the otological assistance surgery.


title: New eigenvalue estimates involving Bessel functions

abstract: Given a compact Riemannian manifold (M n , g) with boundary $\partial$M , we give an estimate for the quotient $\partial$M f d$\mu$ g M f d$\mu$ g , where f is a smooth positive function defined on M that satisfies some inequality involving the scalar Laplacian. By the mean value lemma established in [37], we provide a differential inequality for f which, under some curvature assumptions, can be interpreted in terms of Bessel functions. As an application of our main result, a direct proof is given of the Faber-Krahn inequalities for Dirichlet and Robin Laplacian. Also, a new estimate is established for the eigenvalues of the Dirac operator that involves a positive root of Bessel function besides the scalar curvature. Independently, we extend the Robin Laplacian on functions to differential forms. We prove that this natural extension defines a self-adjoint and elliptic operator whose spectrum is discrete and consists of positive real eigenvalues. In particular, we characterize its first eigenvalue and provide a lower bound of it in terms of Bessel functions.


title: Universality of the Poincar\'e gravitational form factor constraints

abstract: Relativistic spin states are convention dependent. In this work we prove that the zero momentum-transfer limits of the leading two form factors in the decomposition of the energy-momentum tensor matrix elements are independent of this choice. In particular, we demonstrate that these constraints are insensitive to whether the corresponding states are massive or not, and that they arise purely due to the Poincar\'e covariance of the states.


title: "Sufficiently Advanced Technology" for Gravitational Wave Detection

abstract: The science fiction writer Arthur C. Clarke wrote, "Any sufficiently advanced technology is indistinguishable from magic." While not magical by any means, the technology used to detect gravitational waves starting in 2015 is surely sufficiently advanced to be remarkable by any ordinary standard. That technology was developed over a period of almost six decades; the people who were directly involved numbered in the thousands. In this article, I give an idiosyncratic account of the history, with a focus on the question of how people learned what measurement technology would be "sufficiently advanced" to succeed in detecting gravitational waves.


title: Tripartite Heterogeneous Graph Propagation for Large-scale Social Recommendation

abstract: Graph Neural Networks (GNNs) have been emerging as a promising method for relational representation including recommender systems. However, various challenging issues of social graphs hinder the practical usage of GNNs for social recommendation, such as their complex noisy connections and high heterogeneity. The oversmoothing of GNNs is an obstacle of GNN-based social recommendation as well. Here we propose a new graph embedding method Heterogeneous Graph Propagation (HGP) to tackle these issues. HGP uses a group-user-item tripartite graph as input to reduce the number of edges and the complexity of paths in a social graph. To solve the oversmoothing issue, HGP embeds nodes under a personalized PageRank based propagation scheme, separately for group-user graph and user-item graph. Node embeddings from each graph are integrated using an attention mechanism. We evaluate our HGP on a large-scale real-world dataset consisting of 1,645,279 nodes and 4,711,208 edges. The experimental results show that HGP outperforms several baselines in terms of AUC and F1-score metrics.


title: Modelling Regional Crime Risk using Directed Graph of Check-ins

abstract: The location-based social network, Foursquare, reflects the human activities of a city. The mobility dynamics inferred from Foursquare helps us understanding urban social events like crime In this paper, we propose a directed graph from the aggregated movement between regions using Foursquare data. We derive region risk factor from the movement direction, quantity and crime history in different periods of the day. Later, we propose a new set of features, DIrected graph Flow FEatuRes (DIFFER) which are associated with region risk factor. The reliable correlations between DIFFER and crime count are observed. We verify the effectiveness of the DIFFER in monthly crime count using Linear, XGBoost, and Random Forest regression in two cities, Chicago and New York City.


title: Linking Physicians to Medical Research Results via Knowledge Graph Embeddings and Twitter

abstract: Informing professionals about the latest research results in their field is a particularly important task in the field of health care, since any development in this field directly improves the health status of the patients. Meanwhile, social media is an infrastructure that allows public instant sharing of information, thus it has recently become popular in medical applications. In this study, we apply Multi Distance Knowledge Graph Embeddings (MDE) to link physicians and surgeons to the latest medical breakthroughs that are shared as the research results on Twitter. Our study shows that using this method physicians can be informed about the new findings in their field given that they have an account dedicated to their profession.


title: Multiplex graph matching matched filters

abstract: We consider the problem of detecting a noisy induced multiplex template network in a larger multiplex background network. Our approach, which extends the framework of Sussman et al. (2019) to the multiplex setting, leverages a multiplex analogue of the classical graph matching problem to use the template as a matched filter for efficiently searching the background for candidate template matches. The effectiveness of our approach is demonstrated both theoretically and empirically, with particular attention paid to the potential benefits of considering multiple channels.


title: Hyperlink Regression via Bregman Divergence

abstract: A collection of $U \: (\in \mathbb{N})$ data vectors is called a $U$-tuple, and the association strength among the vectors of a tuple is termed as the \emph{hyperlink weight}, that is assumed to be symmetric with respect to permutation of the entries in the index. We herein propose Bregman hyperlink regression (BHLR), which learns a user-specified symmetric similarity function such that it predicts the tuple's hyperlink weight from data vectors stored in the $U$-tuple. BHLR is a simple and general framework for hyper-relational learning, that minimizes Bregman-divergence (BD) between the hyperlink weights and estimated similarities defined for the corresponding tuples; BHLR encompasses various existing methods, such as logistic regression ($U=1$), Poisson regression ($U=1$), link prediction ($U=2$), and those for representation learning, such as graph embedding ($U=2$), matrix factorization ($U=2$), tensor factorization ($U \geq 2$), and their variants equipped with arbitrary BD. Nonlinear functions (e.g., neural networks), can be employed for the similarity functions. However, there are theoretical challenges such that some of different tuples of BHLR may share data vectors therein, unlike the i.i.d. setting of classical regression. We address these theoretical issues, and proved that BHLR equipped with arbitrary BD and $U \in \mathbb{N}$ is (P-1) statistically consistent, that is, it asymptotically recovers the underlying true conditional expectation of hyperlink weights given data vectors, and (P-2) computationally tractable, that is, it is efficiently computed by stochastic optimization algorithms using a novel generalized minibatch sampling procedure for hyper-relational data. Consequently, theoretical guarantees for BHLR including several existing methods, that have been examined experimentally, are provided in a unified manner.


title: An extension of the second order dynamical system that models Nesterov's convex gradient method

abstract: In this paper we deal with a general second order continuous dynamical system associated to a convex minimization problem with a Fr\`echet differentiable objective function. We show that inertial algorithms, such as Nesterov's algorithm, can be obtained via the natural explicit discretization from our dynamical system. Our dynamical system can be viewed as a perturbed version of the heavy ball method with vanishing damping, however the perturbation is made in the argument of the gradient of the objective function. This perturbation seems to have a smoothing effect for the energy error and eliminates the oscillations obtained for this error in the case of the heavy ball method with vanishing damping, as some numerical experiments show. We prove that the value of the objective function in a generated trajectory converges in order O(1/t^2) to the global minimum of the objective function. Moreover, we obtain that a trajectory generated by the dynamical system converges to a minimum point of the objective function.


title: Alternative Blockmodelling

abstract: Many approaches have been proposed to discover clusters within networks. Community finding field encompasses approaches which try to discover clusters where nodes are tightly related within them but loosely related with nodes of other clusters. However, a community network configuration is not the only possible latent structure in a graph. Core-periphery and hierarchical network configurations are valid structures to discover in a relational dataset. On the other hand, a network is not completely explained by only knowing the membership of each node. A high level view of the inter-cluster relationships is needed. Blockmodelling techniques deal with these two issues. Firstly, blockmodelling allows finding any network configuration besides to the well-known community structure. Secondly, blockmodelling is a summary representation of a network which regards not only membership of nodes but also relations between clusters. Finally, a unique summary representation of a network is unlikely. Networks might hide more than one blockmodel. Therefore, our proposed problem aims to discover a secondary blockmodel representation of a network that is of good quality and dissimilar with respect to a given blockmodel. Our methodology is presented through two approaches, (a) inclusion of cannot-link constraints and (b) dissimilarity between image matrices. Both approaches are based on non-negative matrix factorisation NMF which fits the blockmodelling representation. The evaluation of these two approaches regards quality and dissimilarity of the discovered alternative blockmodel as these are the requirements of the problem.


title: Beyond the mini-solar maximum of solar cycle 24: Declining solar magnetic fields and the response of the terrestrial magnetosphere

abstract: The present study examines the response of the terrestrial magnetosphere to the long-term steady declining trends observed in solar magnetic fields and solar wind micro-turbulence levels since mid-1990's that has been continuing beyond the mini-solar maximum of cycle 24. A detailed analysis of the response of the terrestrial magnetosphere has been carried out by studying the extent and shape of the Earth's magnetopause and bow shock over the past four solar cycles. We estimate sub-solar stand-off distance of the magnetopause and bow shock, and the shape of the magnetopause using numerical as well as empirical models. The computed magnetopause and bow shock stand-off distances have been found to be increasing steadily since around mid-1990's, consistent with the steady declining trend seen in solar magnetic fields and solar wind micro-turbulence levels. Similarly, we find an expansion in the shape of the magnetopause since 1996. The implications of the increasing trend seen in the magnetopause and bow shock stand-off distances are discussed and a forecast of the shape of the magnetopause in 2020, the minimum of cycle 24, has been made. Importantly, we also find two instances between 1968 and 1991 when the magnetopause stand-off distance dropped to values close to 6.6 earth radii, the geostationary orbit, for duration ranging from 9$-$11 hours and one event in 2005, post 1995 when the decline in photospheric fields began. Though there have been no such events since 2005, it represents a clear and present danger to our satellite systems.


title: Lorentz forces induce inhomogeneity and fluxes in active systems

abstract: We consider the nonequilibrium dynamics of a charged active Brownian particle in the presence of a space dependent magnetic field. It has recently been shown that the Lorentz force induces a particle flux perpendicular to density gradients, thus preventing a diffusive description of the dynamics. Whereas a passive system will eventually relax to an equilibrium state, unaffected by the magnetic field, an active system subject to a spatially varying Lorentz force settles into a nonequilibrium steady state characterized by an inhomogeneous density and divergence-free bulk fluxes. A macroscopic flux of charged active particles is induced by the gradient of the magnetic field only and does not require additional symmetric breaking such as density or potential gradients. This stands in marked contrast to similar phenomena in condensed matter such as the classical Hall effect. In a confined geometry we observe circulating fluxes, which can be reversed by inverting the direction of the magnetic field. Our theoretical approach, based on coarse-graining of the Fokker-Planck equation, yields analytical results for the density, fluxes, and polarization in the steady state, all of which are validated by direct computer simulation. We demonstrate that passive tracer particles can be used to measure the essential effects of the Lorentz force on the active particle bath, and we discuss under which conditions the effects of the flux could be observed experimentally.


title: Criteria for single photon sources with variable nonclassicality threshold

abstract: Single photon sources are necessary for optical quantum technology. The nonclassicality of emitted single photons manifests itself in diverse experiments and applications. The applications already require different nonclassical aspects of single photons, however, a suitable hierarchy of criteria is missing. We propose variable experimental tests, based on adjustable linear optical networks and single photon detectors, giving a hierarchy of such nonclassicality criteria. The hierarchy goes beyond Hanbury - Brown - Twiss test of photon antibunching and allows us to compare faithfully quality and performance of single-photon sources via their nonclassical properties.


title: A Simple and Effective Approach for Fine Tuning Pre-trained Word Embeddings for Improved Text Classification

abstract: This work presents a new and simple approach for fine-tuning pretrained word embeddings for text classification tasks. In this approach, the class in which a term appears, acts as an additional contextual variable during the fine tuning process, and contributes to the final word vector for that term. As a result, words that are used distinctively within a particular class, will bear vectors that are closer to each other in the embedding space and will be more discriminative towards that class. To validate this novel approach, it was applied to three Arabic and two English datasets that have been previously used for text classification tasks such as sentiment analysis and emotion detection. In the vast majority of cases, the results obtained using the proposed approach, improved considerably.


title: First-principles study of the charge ordered phase in $\kappa$-D$_3$(Cat-EDT-TTF/ST)$_2$: Stability of $\pi$-electron deuterium coupled ordering in hydrogen-bonded molecular conductors

abstract: We study the electronic and structural properties of the low-temperature ordered phase of hydrogen-bonded molecular conductors, $\kappa$-D$_3$(Cat-EDT-TTF)$_2$ and its selenium-substituted analog $\kappa$-D$_3$(Cat-EDT-ST)$_2$, by means of first-principles density functional theory~(DFT) calculations. In these compounds, the charge ordering in the $\pi$-electron system is coupled with the ordering of the displacements in the deuteriums forming the hydrogen-bond, equally shared by two oxygens in the high-temperature phase. While the structural optimization within the standard DFT method based on the generalized gradient approximation fails to reproduce the structural stability of the charge-ordered (CO) phase, we show that a hybrid functional of Heyd, Scuseria, and Ernzerhof can reproduce structural characters of the CO phase, owing to the more localized nature of the wave functions. Furthermore, using the ability of the hybrid functional to predict the electronic and structural properties, we find a stable noncentrosymmetric CO phase with another pattern of deuterium ordering.


title: Flat Chern Band From Twisted Bilayer MnBi$_2$Te$_4$

abstract: We construct a continuum model for the Moir\'e superlattice of twisted bilayer MnBi$_2$Te$_4$, and study the band structure of the bilayer in both ferromagnetic (FM) and antiferromagnetic (AFM) phases. We find the system exhibits highly tunable Chern bands with Chern number up to $3$. We show that a twist angle of $1^\circ$ turns the highest valence band into a flat band with Chern number $\pm1$ that is isolated from all other bands in both FM and AFM phases. This result provides a promising platform for realizing time-reversal breaking correlated topological phases, such as fractional Chern insulator and $p+ip$ topological superconductor. In addition, our calculation indicates that the twisted stacking facilitates the emergence of quantum anomalous Hall effect in MnBi$_2$Te$_4$.


title: Confident Head Circumference Measurement from Ultrasound with Real-time Feedback for Sonographers

abstract: Manual estimation of fetal Head Circumference (HC) from Ultrasound (US) is a key biometric for monitoring the healthy development of fetuses. Unfortunately, such measurements are subject to large inter-observer variability, resulting in low early-detection rates of fetal abnormalities. To address this issue, we propose a novel probabilistic Deep Learning approach for real-time automated estimation of fetal HC. This system feeds back statistics on measurement robustness to inform users how confident a deep neural network is in evaluating suitable views acquired during free-hand ultrasound examination. In real-time scenarios, this approach may be exploited to guide operators to scan planes that are as close as possible to the underlying distribution of training images, for the purpose of improving inter-operator consistency. We train on free-hand ultrasound data from over 2000 subjects (2848 training/540 test) and show that our method is able to predict HC measurements within 1.81$\pm$1.65mm deviation from the ground truth, with 50% of the test images fully contained within the predicted confidence margins, and an average of 1.82$\pm$1.78mm deviation from the margin for the remaining cases that are not fully contained.


title: Simplicial degree in complex networks. Applications of Topological Data Analysis to Network Science

abstract: Network Science provides a universal formalism for modelling and studying complex systems based on pairwise interactions between agents. However, many real networks in the social, biological or computer sciences involve interactions among more than two agents, having thus an inherent structure of a simplicial complex. We propose new notions of higher-order degrees of adjacency for simplices in a simplicial complex, allowing any dimensional comparison among them and their faces, which as far as we know were lacked in the literature. We introduce multi-parameter boundary and coboundary operators in an oriented simplicial complex and also a novel multi-combinatorial Laplacian is defined, which generalises the graph and combinatorial Laplacian. To illustrate the potential applications of these theoretical results, we perform a structural analysis of higher-order connectivity in simplicial-complex networks by studying the associated distributions with these simplicial degrees in 17 real-world datasets coming from different domains such as coauthor networks, cosponsoring Congress bills, contacts in schools, drug abuse warning networks, e-mail networks or publications and users in online forums. We find rich and diverse higher-order connectivity structures and observe that datasets of the same type reflect similar higher-order collaboration patterns. Furthermore, we show that if we use what we have called the maximal simplicial degree (which counts the distinct maximal communities in which our simplex and all its strict sub-communities are contained), then its degree distribution is, in general, surprisingly different from the classical node degree distribution.


title: Economic Power, Population, and the Size of Astronomical Community

abstract: The number of astronomers for a country registered to the IAU is known to have a correlation with the GDP. However, the robustness of this relationship can be doubted, because the fraction of astronomers joining the IAU differs from country to country. Here we revisit this correlation by using the recent data updated as of 2017, and then we find a similar correlation by using the total enumeration of astronomers and astrophysicists with PhD degrees and working in each country, instead of adopting the number of IAU members. We confirm the existence of two subgroup in the correlation. One group consists of European advanced countries having long history of modern astronomy, while the other group consists of countries having experienced recent rapid economic development. In order to find causation in the correlation, we obtain the long-term variations of the number of astronomers, population, and the GDP for a number of countries to find that the number of astronomers per citizen for recently developing countries has increased more rapidly as GDP per capita increased, than that for fully developed countries. We collect a demographic data of the Korean astronomical community. From these findings we estimate the proper size of the Korean astronomical community by considering the society's economic power and population. The current number of PhD astronomers working in Korea is approximately 310, but it should be 550 that is large enough to be comparable and competitive to the sizes of Spainish, Canadian, and Japanese astronomical communities. We discuss on the way how to overcome the vulnerability of the Korean astronomical community, based on the statistics of national R&D expenditure structure comparing with that of other major advanced countries.


title: Phase-Apodized-Pupil Lyot Coronagraphs for Arbitrary Telescope Pupils

abstract: The phase-apodized-pupil Lyot coronagraph (PAPLC) is a pairing of the apodized-pupil Lyot coronagraph (APLC) and the apodizing phase plate (APP) coronagraph. We describe a numerical optimization method to obtain globally-optimal solutions for the phase apodizers for arbitrary telescope pupils, based on the linear map between complex-amplitude transmission of the apodizer and the electric field in the post-coronagraphic focal plane. PAPLCs with annular focal-plane masks and point-symmetric dark zones perform analogous to their corresponding APLCs. However with a knife-edge focal-plane mask and one-sided dark zones, the PAPLC yields inner working angles as close as $1.4\lambda/D$ at contrasts of $10^{-10}$ and maximum post-coronagraphic throughput of >75% for telescope apertures with central obscurations of up to 30%. We present knife-edge PAPLC designs optimized for the VLT/SPHERE instrument and the LUVOIR-A aperture. These designs show that the knife-edge PAPLC retains its performance, even for realistic telescope pupils with struts, segments and non-circular outer edges.


title: A modelling methodology for social interaction experiments

abstract: Analysis of temporal network data arising from online interactive social experiments is not possible with standard statistical methods because the assumptions of these models, such as independence of observations, are not satisfied. In this paper, we outline a modelling methodology for such experiments where, as an example, we analyse data collected using the Virtual Interaction Application (VIAPPL) --- a software platform for conducting experiments that reveal how social norms and identities emerge through social interaction. We apply our model to show that ingroup favouritism and reciprocity are present in the experiments, and to quantify the strengthening of these behaviours over time. Our method enables us to identify participants whose behaviour is markedly different from the norm. We use the method to provide a visualisation of the data that highlights the level of ingroup favouritism, the strong reciprocal relationships, and the different behaviour of participants in the game. While our methodology was developed with VIAPPL in mind, its usage extends to any type of social interaction data.


title: Towards a Comprehensive Bibliography for SETI

abstract: In this work, we motivate, describe, and announce a living bibliography for academic papers and other works published in the Search for Extraterrestrial Intelligence (SETI). The bibliography makes use of bibliographic groups (bibgroups) in the NASA Astrophysics Data System (ADS), allowing it to be accessed and searched by any interested party, and is composed only of works which have a presence on the ADS. We establish criteria that describe the scope of our bibliography, which we define as any academic work which broadly: 1) advances knowledge within SETI, 2) deals with topics that are fundamentally related to or about SETI, or 3) is useful for the better understanding of SETI, and which has a presence on ADS. We discuss the future work needed to continue the development of the bibliography. The bibliography can be found by using the bibgroup field (bibgroup: SETI) in the ADS search engine.