title: Giant optical activity and Kerr effect in type-I and type-II Weyl semimetals

abstract: We explore optical activity in thin films and bulk of type-I and type II Weyl semimetals (WSM), and demonstrate the existence of a giant Kerr effect in both. In time-reversal symmetry broken WSM thin films, the polarization rotation is caused by the optical Hall conductivity including the anomalous Hall term. The Kerr angle is found to be $\propto Q/\omega$, with $Q$ and $\omega$ being the Weyl node separation and the optical frequency, respectively. In contrast, the optical activity in the bulk WSM is dominated by axion electrodynamics, which persists even in the Pauli blocked regime of no optical transitions. In bulk WSM, $Q$ acts analogous to the magnetization in magnetic materials, leading to large `polar Kerr effect' (linear in $Q$) when light is incident on WSM surface without Fermi arc states, and the `Voigt effect' (quadratic in $Q$), when light is incident on surface with Fermi arc states.


title: A spectral decomposition of the attractor of piecewise contracting maps of the interval

abstract: We study the asymptotic dynamics of piecewise contracting maps defined on a compact interval. For maps that are not necessarily injective, but have a finite number of local extrema and discontinuity points, we prove the existence of a decomposition of the support of the asymptotic dynamics into a finite number of minimal components. Each component is either a periodic orbit or a minimal Cantor set and such that the $\omega$-limit set of (almost) every point in the interval is exactly one of these components. Moreover, we show that each component is the $\omega$-limit set, or the closure of the orbit, of a one-sided limit of the map at a discontinuity point or at a local extremum.


title: EHAAS: Energy Harvesters As A Sensor for Place Recognition on Wearables

abstract: A wearable based long-term lifelogging system is desirable for the purpose of reviewing and improving users' lifestyle habits. Energy harvesting (EH) is a promising means for realizing sustainable lifelogging. However, present EH technologies suffer from instability of the generated electricity caused by changes of environment, e.g., the output of a solar cell varies based on its material, light intensity, and light wavelength. In this paper, we leverage this instability of EH technologies for other purposes, in addition to its use as an energy source. Specifically, we propose to determine the variation of generated electricity as a sensor for recognizing "places" where the user visits, which is important information in the lifelogging system. First, we investigate the amount of generated electricity of selected energy harvesting elements in various environments. Second, we design a system called EHAAS (Energy Harvesters As A Sensor) where energy harvesting elements are used as a sensor. With EHAAS, we propose a place recognition method based on machine-learning and implement a prototype wearable system. Our prototype evaluation confirms that EHAAS achieves a place recognition accuracy of 88.5% F-value for nine different indoor and outdoor places. This result is better than the results of existing sensors (3-axis accelerometer and brightness). We also clarify that only two types of solar cells are required for recognizing a place with 86.2% accuracy.


title: The VANDELS survey: the role of ISM and galaxy physical properties on the escape of Ly{\alpha} emission in z $\sim$ 3.5 star-forming galaxies

abstract: We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.


title: Dynamics of front propagation in a compliant channel

abstract: Front-propagating systems provide some of the most fundamental physical examples of interfacial instability and pattern formation. However, their nonlinear dynamics are rarely addressed. Here, we present an experimental study of air displacing a viscous fluid within a collapsed, compliant channel - a model system for pulmonary airway reopening. We show that compliance induces fingering instabilities absent in the rigid channel and we present the first experimental observations of the counter-intuitive 'pushing' behaviour previously predicted numerically, for which a reduction in air pressure results in faster flow. We find that pushing modes are unstable and moreover, that the dynamics of the air-fluid front involves a host of transient finger shapes over a significant range of experimental parameters.


title: KLEVER: An experiment to measure \boldmath{${\rm BR}(K_L \to \pi^0 \nu \overline{\nu})$} at the CERN SPS

abstract: The KLEVER experiment aims to measure ${\rm BR}(K_L \to \pi^0 \nu \overline{\nu})$, supplementing the ongoing NA62 measurement of ${\rm BR}(K^+ \to \pi^+ \nu \overline{\nu})$, to provide new input on CKM unitarity and potentially new physics. KLEVER is undergoing continuous development, with particular efforts focused on the design of the target and the beamline. As described here, adaptations are required relative to the K12 beamline in its current format, and a series of simulations has been performed to ensure that an adequate particle flux can be achieved while simultaneously suppressing problematic backgrounds.


title: Electric Field Tunable Correlated States and Magnetic Phase Transitions in Twisted Bilayer-Bilayer Graphene

abstract: The recent discovery of correlated insulator states and superconductivity in magic-angle twisted bilayer graphene has paved the way to the experimental investigation of electronic correlations in tunable flat band systems realized in twisted van der Waals heterostructures. This novel twist angle degree of freedom and control should be generalizable to other 2D systems, which may exhibit similar correlated physics behavior while at the same time enabling new techniques to tune and control the strength of electron-electron interactions. Here, we report on a new highly tunable correlated system based on small-angle twisted bilayer-bilayer graphene (TBBG), consisting of two rotated sheets of Bernal-stacked bilayer graphene. We find that TBBG exhibits a rich phase diagram, with tunable correlated insulators states that are highly sensitive to both twist angle and to the application of an electric displacement field, the latter reflecting the inherent polarizability of Bernal-stacked bilayer graphene. We find correlated insulator states that can be switched on and off by the displacement field at all integer electron fillings of the moir\'{e} unit cell. The response of these correlated states to magnetic fields points towards evidence of electrically switchable magnetism. Moreover, the strong dependence of the resistance at low temperature and near the correlated insulator states indicates possible proximity to a superconducting phase. Furthermore, in the regime of lower twist angles, TBBG shows multiple sets of flat bands near charge neutrality, resulting in numerous correlated states corresponding to half-filling of each of these flat bands. Our results pave the way to the exploration of novel twist-angle and electric-field controlled correlated phases of matter in novel multi-flat band twisted superlattices.


title: A Graph-structured Dataset for Wikipedia Research

abstract: Wikipedia is a rich and invaluable source of information. Its central place on the Web makes it a particularly interesting object of study for scientists. Researchers from different domains used various complex datasets related to Wikipedia to study language, social behavior, knowledge organization, and network theory. While being a scientific treasure, the large size of the dataset hinders pre-processing and may be a challenging obstacle for potential new studies. This issue is particularly acute in scientific domains where researchers may not be technically and data processing savvy. On one hand, the size of Wikipedia dumps is large. It makes the parsing and extraction of relevant information cumbersome. On the other hand, the API is straightforward to use but restricted to a relatively small number of requests. The middle ground is at the mesoscopic scale when researchers need a subset of Wikipedia ranging from thousands to hundreds of thousands of pages but there exists no efficient solution at this scale. In this work, we propose an efficient data structure to make requests and access subnetworks of Wikipedia pages and categories. We provide convenient tools for accessing and filtering viewership statistics or "pagecounts" of Wikipedia web pages. The dataset organization leverages principles of graph databases that allows rapid and intuitive access to subgraphs of Wikipedia articles and categories. The dataset and deployment guidelines are available on the LTS2 website \url{https://lts2.epfl.ch/Datasets/Wikipedia/}.


title: A Girsanov approach to slow parameterizing manifolds in the presence of noise

abstract: We consider a three-dimensional slow-fast system arising in fluid dynamics with quadratic nonlinearity and additive noise. The associated deterministic system of this stochastic differential equation (SDE) exhibits a periodic orbit and a slow manifold. We show that in presence of noise, the deterministic slow manifold can be viewed as an approximate parameterization of the fast variable of the SDE in terms of the slow variables, for certain parameter regimes. We exploit this fact to obtain a two dimensional reduced model from the original stochastic system, which results into a Hopf normal form with additive noise. Both, the original as well as the reduced system admit ergodic invariant measures describing their respective long-time behaviour. It is then shown that for a suitable Wasserstein metric on a subset of the space of probability measures on the phase space, the discrepancy between the marginals along the radial component of each invariant measure is controlled by a parameterization defect measuring the quality of the parameterization. An important technical tool to arrive at this result is Girsanov's theorem that allows us to derive such error estimates in presence of an oscillatory instability. This approach is finally extended to parameter regimes for which the variable to parameterize is no longer evolving on a faster timescale than that of the resolved variables. There also, error estimates involving the Wasserstein metric are derived but this time for reduced systems obtained from stochastic parameterizing manifolds involving path-dependent coefficients to cope with such challenging regimes.


title: LMI Properties and Applications in Systems, Stability, and Control Theory

abstract: Linear matrix inequalities (LMIs) commonly appear in systems, stability, and control applications. Many analysis and synthesis problems in these areas can be solved as feasibility or optimization problems subject to LMI constraints. Although most well-known LMI properties and manipulation tricks, such as the Schur complement and the congruence transformation, can be found in standard references, many useful LMI properties are scattered throughout the literature. The purpose of this document is to collect and organize properties, tricks, and applications related to LMIs from a number of references together in a single document. In this sense, the document can be thought of as an "LMI encyclopedia" or "LMI cookbook." Proofs of the properties presented in this document are not included when they can be found in the cited references in the interest of brevity. Illustrative examples are included whenever necessary to fully explain a certain property. Multiple equivalent forms of LMIs are often presented to give the reader a choice of which form may be best suited for a particular problem at hand. The equivalency of some of the LMIs in this document may be straightforward to more experienced readers, but the authors believe that some readers may benefit from the presentation of multiple equivalent LMIs.


title: Contextual Bandits with Random Projection

abstract: Contextual bandits with linear payoffs, which are also known as linear bandits, provide a powerful alternative for solving practical problems of sequential decisions, e.g., online advertisements. In the era of big data, contextual data usually tend to be high-dimensional, which leads to new challenges for traditional linear bandits mostly designed for the setting of low-dimensional contextual data. Due to the curse of dimensionality, there are two challenges in most of the current bandit algorithms: the first is high time-complexity; and the second is extreme large upper regret bounds with high-dimensional data. In this paper, in order to attack the above two challenges effectively, we develop an algorithm of Contextual Bandits via RAndom Projection (\texttt{CBRAP}) in the setting of linear payoffs, which works especially for high-dimensional contextual data. The proposed \texttt{CBRAP} algorithm is time-efficient and flexible, because it enables players to choose an arm in a low-dimensional space, and relaxes the sparsity assumption of constant number of non-zero components in previous work. Besides, we provide a linear upper regret bound for the proposed algorithm, which is associated with reduced dimensions.


title: Comprehensive tunneling spectroscopy of quasi-freestanding MoS$_2$ on graphene on Ir(111)

abstract: We apply scanning tunneling spectroscopy to determine the bandgaps of mono-, bi- and trilayer MoS$_2$ grown on a graphene single crystal on Ir(111). Besides the typical scanning tunneling spectroscopy at constant height, we employ two additional spectroscopic methods giving extra sensitivity and qualitative insight into the $k$-vector of the tunneling electrons. Employing this comprehensive set of spectroscopic methods in tandem, we deduce a bandgap of $2.53\pm0.08$ eV for the monolayer. This is close to the predicted values for freestanding MoS$_2$ and larger than is measured for MoS$_2$ on other substrates. Through precise analysis of the `comprehensive' tunneling spectroscopy we also identify critical point energies in the mono- and bilayer MoS$_2$ band structures. These compare well with their calculated freestanding equivalents, evidencing the graphene/Ir(111) substrate as an excellent environment upon which to study the many feted electronic phenomena of monolayer MoS$_2$ and similar materials. Additionally, this investigation serves to expand the fledgling field of the comprehensive tunneling spectroscopy technique itself.


title: Resonantly induced friction in driven nanomechanical systems

abstract: We propose a new mechanism of friction in resonantly driven vibrational systems. The form of the friction force follows from the time- and spatial-symmetry arguments. We consider a microscopic mechanism of this resonant force in nanomechanical systems. The friction can be negative, leading to an instability of forced vibrations of a nanoresonator and the onset of self-sustained oscillations in the rotating frame.


title: Almost Tight Lower Bounds for Hard Cutting Problems in Embedded Graphs

abstract: We prove essentially tight lower bounds, conditionally to the Exponential Time Hypothesis, for two fundamental but seemingly very different cutting problems on surface-embedded graphs: the Shortest Cut Graph problem and the Multiway Cut problem. A cut graph of a graph $G$ embedded on a surface $S$ is a subgraph of $G$ whose removal from $S$ leaves a disk. We consider the problem of deciding whether an unweighted graph embedded on a surface of genus $g$ has a cut graph of length at most a given value. We prove a time lower bound for this problem of $n^{\Omega(g/\log g)}$ conditionally to ETH. In other words, the first $n^{O(g)}$-time algorithm by Erickson and Har-Peled [SoCG 2002, Discr.\ Comput.\ Geom.\ 2004] is essentially optimal. We also prove that the problem is W[1]-hard when parameterized by the genus, answering a 17-year old question of these authors. A multiway cut of an undirected graph $G$ with $t$ distinguished vertices, called terminals, is a set of edges whose removal disconnects all pairs of terminals. We consider the problem of deciding whether an unweighted graph $G$ has a multiway cut of weight at most a given value. We prove a time lower bound for this problem of $n^{\Omega(\sqrt{gt + g^2+t}/\log(g+t))}$, conditionally to ETH, for any choice of the genus $g\ge0$ of the graph and the number of terminals $t\ge4$. In other words, the algorithm by the second author [Algorithmica 2017] (for the more general multicut problem) is essentially optimal; this extends the lower bound by the third author [ICALP 2012] (for the planar case). Reductions to planar problems usually involve a grid-like structure. The main novel idea for our results is to understand what structures instead of grids are needed if we want to exploit optimally a certain value $g$ of the genus.


title: Tracking the ultrafast motion of an antiferromagnetic order parameter

abstract: The unique functionalities of antiferromagnets offer promising routes to advance information technology. Their compensated magnetic order leads to spin resonances in the THz-regime, which suggest the possibility to coherently control antiferromagnetic (AFM) devices orders of magnitude faster than traditional electronics. However, the required time resolution, complex sublattice interations and the relative inaccessibility of the AFM order parameter pose serious challenges to studying AFM spin dynamics. Here, we reveal the temporal evolution of an AFM order parameter directly in the time domain. We modulate the AFM order in hexagonal YMnO$_\mathrm{3}$ by coherent magnon excitation and track the ensuing motion of the AFM order parameter using time-resolved optical second-harmonic generation (SHG). The dynamic symmetry reduction by the moving order parameter allows us to separate electron dynamics from spin dynamics. As transient symmetry reductions are common to coherent excitations, we have a general tool for tracking the ultrafast motion of an AFM order parameter.


title: Iterated Extended Kalman Smoother-based Variable Splitting for $L_1$-Regularized State Estimation

abstract: In this paper, we propose a new framework for solving state estimation problems with an additional sparsity-promoting $L_1$-regularizer term. We first formulate such problems as minimization of the sum of linear or nonlinear quadratic error terms and an extra regularizer, and then present novel algorithms which solve the linear and nonlinear cases. The methods are based on a combination of the iterated extended Kalman smoother and variable splitting techniques such as alternating direction method of multipliers (ADMM). We present a general algorithmic framework for variable splitting methods, where the iterative steps involving minimization of the nonlinear quadratic terms can be computed efficiently by iterated smoothing. Due to the use of state estimation algorithms, the proposed framework has a low per-iteration time complexity, which makes it suitable for solving a large-scale or high-dimensional state estimation problem. We also provide convergence results for the proposed algorithms. The experiments show the promising performance and speed-ups provided by the methods.


title: Single-step Options for Adversary Driving

abstract: In this paper, we use reinforcement learning for safety driving in adversary settings. In our work, the knowledge in state-of-art planning methods is reused by single-step options whose action suggestions are compared in parallel with primitive actions. We show two advantages by doing so. First, training this reinforcement learning agent is easier and faster than training the primitive-action agent. Second, our new agent outperforms the primitive-action reinforcement learning agent, human testers as well as the state-of-art planning methods that our agent queries as skill options.


title: Neutrino Echoes from Multimessenger Transient Sources

abstract: The detection of the high-energy neutrino event, IceCube-170922A, demonstrated that multimessenger particle astrophysics triggered by neutrino alerts is feasible. We consider time delay signatures caused by secret neutrino interactions with the cosmic neutrino background and dark matter and suggest that these can be used as a novel probe of neutrino interactions beyond the Standard Model (BSM). The tests with BSM-induced neutrino echoes are distinct from existing constraints from the spectral modification and will be enabled by multimessenger observations of bright neutrino transients with future experiments such as IceCube-Gen2, KM3Net, and Hyper-Kamiokande. The constraints are complementary to those from accelerator and laboratory experiments and powerful for testing various particle models that explain tensions prevailing in the cosmological data.


title: Joint User Association and Resource Allocation in Heterogeneous Cellular Networks: Comparison of Two Modeling Approaches

abstract: The performance of different combinations of user association (UA) and resource allocation (RA) in heterogeneous cellular networks has been extensively studied using a classic modeling approach based on system snapshots. There have been also many studies focusing on the dynamics of the system using queueing models. These modeling approaches are rarely compared with each other though they each bring different insights to the design problem. In this paper, we consider a queueing model-based approach to study the interplay of UA and RA, and compare the results to those obtained using snapshot models. Specifically, we formulate three different joint UA and RA optimization problems corresponding to the following three performance metrics: the maximum achievable arrival rate, the average system delay, and the maximum per-user delay. These problems are non-convex integer programs. We have therefore developed numerical techniques to compute either their exact solutions or tight lower bounds. We obtain results for different combinations of RA and UA schemes, and compare the trends with those obtained via the snapshot approach. The trends on RA are very similar, which we take as a cross-validation of the two modeling approaches for this kind of problem. The trends on user association are somewhat different which indicates a lack of robustness of the results and the need for a careful validation of UA models.


title: Heterogeneous Prestressed Precast Beams Multiperiod Production Planning Problem: Modeling and Solution Methods

abstract: A prestressed precast beam is a type of beam that is stretched with traction elements. A common task in a factory of prestressed precast beams involves fulfilling, within a time horizon, the demand ordered by clients. A typical order includes beams of different lengths and types, with distinct beams potentially requiring different curing periods. We refer to the problem of planning such production as Heterogeneous Prestressed Precast Beams Multiperiod Production Planning (HPPBMPP). We formally define the HPPBMPP, argue its NP-hardness, and introduce four novel integer programming models for its solution and a size reduction heuristic (SRH). We propose six priority rules to produce feasible solutions. We perform computational tests on a set of synthetic instances that are based on data from a real-world scenario and discuss a case study. Our experiments suggest that the models can optimally solve small instances, while the SRH can produce high-quality solutions for most instances.


title: Praesepe (NGC 2632) and its tidal tails

abstract: Within a 400 pc sphere around the Sun, we search for Praesepe's tidal tails in the Gaia DR2 dataset. We used a modified convergent-point method to search for stars with space velocities close to the space velocity of the Praesepe cluster. We find a clear indication for the existence of Praesepe's tidal tails, both extending up to 165~pc from the centre of the cluster. A total of 1393 stars populate the cluster and its tails, giving a total mass of 794 M_Sun. We determined a tidal radius of 10.77 pc for the cluster and a tidal mass of 483 M_Sun. The corresponding half-mass radius is 4.8 pc. We also found clear indication for mass segregation in the cluster. The tidal tails outside 2 tidal radii are populated by 389 stars. The total contamination of our sample by field stars lies between 50 to 100 stars or 3.6 to 7.2 per cent. We used an astrometrically and photometrically clean sub-sample of Gaia DR2 which makes our Praesepe sample incomplete beyond M_G ~ 12.0 mag, which corresponds to about 0.25 M_Sun. A comparison with an N-body model of the cluster and its tails shows remarkably good coincidence. Three new white dwarfs are found in the tails.


title: Autocovariance Varieties of Moving Average Random Fields

abstract: We study the autocovariance functions of moving average random fields over the integer lattice $\mathbb{Z}^d$ from an algebraic perspective. These autocovariances are parametrized polynomially by the moving average coefficients, hence tracing out algebraic varieties. We derive dimension and degree of these varieties and we use their algebraic properties to obtain statistical consequences such as identifiability of model parameters. We connect the problem of parameter estimation to the algebraic invariants known as euclidean distance degree and maximum likelihood degree. Throughout, we illustrate the results with concrete examples. In our computations we use tools from commutative algebra and numerical algebraic geometry.


title: The impact of non-ideal surfaces on the solid-water interaction: a time-resolved adsorption study

abstract: The initial interaction of water with semiconductors determines the electronic structure of the solid-liquid interface. The exact nature of this interaction is, however, often unknown. Here, we study gallium phosphide-based surfaces exposed to H2O by means of in situ reflection anisotropy spectroscopy. We show that the introduction of typical imperfections in the form of surface steps or trace contaminants not only changes the dynamics of the interaction, but also its qualitative nature. This emphasises the challenges for the comparability of experiments with (idealised) electronic structure models for electrochemistry.


title: Universal Area Law in Turbulence

abstract: We re-visit the Area Law in Turbulence discovered many years ago \cite{M93} and verified recently in numerical experiments\cite{S19}. We derive this law in a simpler way, at the same time outlining the limits of its applicability. Using the PDF for velocity circulation as a functional of the loop in coordinate space, we obtain explicit formulas for vorticity correlations in presence of velocity circulation. These functions are related to the shape of the scaling function of the PDF as well as the shape of the minimal surface inside the loop. The background of velocity circulation does not eliminate turbulence but makes observable quantities in inertial range \textbf{calculable}. The scaling dimension of velocity circulation as a function of large area remains unknown. Numerical experiments \cite{S19} suggest transition for log-log derivative of circulation moments $\left<\Gamma^p\right>$ by the loop area from Kolmogorov index $\frac{2p}{3}$ at $p <4$ down to approximately $0.58 p$ for $4 \leq p \leq 10$ within available Reynolds numbers. We argue that Area Law applies to these moments only in the limit $p\rightarrow \infty$ when they are dominated by the tails of the PDF. So, these numerical experiments suggest that the scaling index in Area law is less then $\frac{2}{3}$.


title: Maximum Nullity and Forcing Number on Graphs with Maximum Degree at most Three

abstract: A dynamic coloring of the vertices of a graph $G$ starts with an initial subset $F$ of colored vertices, with all remaining vertices being non-colored. At each time step, a colored vertex with exactly one non-colored neighbor forces this non-colored neighbor to be colored. The initial set $F$ is called a forcing set of $G$ if, by iteratively applying the forcing process, every vertex in $G$ becomes colored. The forcing number of a graph $G$, denoted by $F(G)$, is the cardinality of a minimum forcing set of $G$. The maximum nullity of $G$, denoted by $M(G)$, is defined to be the largest possible nullity over all real symmetric matrices $A$ whose $a_{ij} \neq 0$ for $i \neq j$, whenever two vertices $u_{i}$ and $u_{j}$ of $G$ are adjacent. In this paper, we characterize all graphs $G$ of order $n$, maximum degree at most three, and $F(G)=3$. Also we classify these graphs with their maximum nullity.


title: Scaling methods for accelerating kinetic Monte Carlo simulations of chemical reaction networks

abstract: Various kinetic Monte Carlo algorithms become inefficient when some of the population sizes in a system are large, which gives rise to a large number of reaction events per unit time. Here, we present a new acceleration algorithm based on adaptive and heterogeneous scaling of reaction rates and stoichiometric coefficients. The algorithm is conceptually related to the commonly used idea of accelerating a stochastic simulation by considering a sub-volume $\lambda \Omega$ ($0<\lambda<1$) within a system of interest, which reduces the number of reaction events per unit time occurring in a simulation by a factor $1/\lambda$ at the cost of greater error in unbiased estimates of first moments and biased overestimates of second moments. Our new approach offers two unique benefits. First, scaling is adaptive and heterogeneous, which eliminates the pitfall of overaggressive scaling. Second, there is no need for an \emph{a priori} classification of populations as discrete or continuous (as in a hybrid method), which is problematic when discreteness of a chemical species changes during a simulation. The method requires specification of only a single algorithmic parameter, $N_c$, a global critical population size above which populations are effectively scaled down to increase simulation efficiency. The method, which we term partial scaling, is implemented in the open-source BioNetGen software package. We demonstrate that partial scaling can significantly accelerate simulations without significant loss of accuracy for several published models of biological systems. These models characterize activation of the mitogen-activated protein kinase ERK, prion protein aggregation, and T-cell receptor signaling.


title: Plug and play methods for magnetic resonance imaging (long version)

abstract: Magnetic Resonance Imaging (MRI) is a non-invasive diagnostic tool that provides excellent soft-tissue contrast without the use of ionizing radiation. Compared to other clinical imaging modalities (e.g., CT or ultrasound), however, the data acquisition process for MRI is inherently slow, which motivates undersampling and thus drives the need for accurate, efficient reconstruction methods from undersampled datasets. In this article, we describe the use of "plug-and-play" (PnP) algorithms for MRI image recovery. We first describe the linearly approximated inverse problem encountered in MRI. Then we review several PnP methods, where the unifying commonality is to iteratively call a denoising subroutine as one step of a larger optimization-inspired algorithm. Next, we describe how the result of the PnP method can be interpreted as a solution to an equilibrium equation, allowing convergence analysis from the equilibrium perspective. Finally, we present illustrative examples of PnP methods applied to MRI image recovery.


title: An optimization model for line planning and timetabling in automated urban metro subway networks

abstract: In this paper we present a Mixed Integer Nonlinear Programming model that we developed as part of a pilot study requested by the R&D company Metrolab in order to design tools for finding solutions for line planning and timetable situations in automated urban metro subway networks. Our model incorporates important factors in public transportation systems from both, a cost-oriented and a passenger-oriented perspective, as time-dependent demands, interchange stations, short-turns and technical features of the trains in use. The incoming flows of passengers are modeled by means of piecewise linear demand functions which are parameterized in terms of arrival rates and bulk arrivals. Decisions about frequencies, train capacities, short-turning and timetables for a given planning horizon are jointly integrated to be optimized in our model. Finally, a novel Math-Heuristic approach is proposed to solve the problem. The results of extensive computational experiments are reported to show its applicability and effectiveness to handle real-world subway networks


title: Totally Asynchronous Distributed Quadratic Programming with Independent Stepsizes and Regularizations

abstract: Quadratic programs arise in robotics, communications, smart grids, and many other applications. As these problems grow in size, finding solutions becomes much more computationally demanding, and new algorithms are needed to efficiently solve them. Targeting large-scale problems, we develop a multi-agent quadratic programming framework in which each agent updates only a small number of the total decision variables in a problem. Agents communicate their updated values to each other, though we do not impose any restrictions on the timing with which they do so, nor on the delays in these transmissions. Furthermore, we allow weak parametric coupling among agents, in the sense that they are free to independently choose their stepsizes, subject to mild restrictions. We show that these stepsize restrictions depend upon a problem's condition number. We further provide the means for agents to independently regularize the problem they solve, thereby improving condition numbers and, as we will show, convergence properties, while preserving agents' independence in selecting parameters. Simulation results are provided to demonstrate the success of this framework on a practical quadratic program.


title: The importance of better models in stochastic optimization

abstract: Standard stochastic optimization methods are brittle, sensitive to stepsize choices and other algorithmic parameters, and they exhibit instability outside of well-behaved families of objectives. To address these challenges, we investigate models for stochastic minimization and learning problems that exhibit better robustness to problem families and algorithmic parameters. With appropriately accurate models---which we call the aProx family---stochastic methods can be made stable, provably convergent and asymptotically optimal; even modeling that the objective is nonnegative is sufficient for this stability. We extend these results beyond convexity to weakly convex objectives, which include compositions of convex losses with smooth functions common in modern machine learning applications. We highlight the importance of robustness and accurate modeling with a careful experimental evaluation of convergence time and algorithm sensitivity.


title: Spectral enclosures for non-self-adjoint discrete Schr\"odinger operators

abstract: We study location of eigenvalues of one-dimensional discrete Schr\"odinger operators with complex $\ell^{p}$-potentials for $1\leq p\leq \infty$. In the case of $\ell^{1}$-potentials, the derived bound is shown to be optimal. For $p>1$, two different spectral bounds are obtained. The method relies on the Birman-Schwinger principle and various techniques for estimations of the norm of the Birman-Schwinger operator.


title: Column2Vec: Structural Understanding via Distributed Representations of Database Schemas

abstract: We present Column2Vec, a distributed representation of database columns based on column metadata. Our distributed representation has several applications. Using known names for groups of columns (i.e., a table name), we train a model to generate an appropriate name for columns in an unnamed table. We demonstrate the viability of our approach using schema information collected from open source applications on GitHub.


title: Exploring the Kondo effect of an extended impurity with chains of Co adatoms in a magnetic field

abstract: Motivated by recent STM experiments, we explore the magnetic field induced Kondo effect that takes place at symmetry protected level crossings in finite Co adatom chains. We argue that the effective two-level system realized at a level crossing acts as an extended impurity coupled to the conduction electrons of the substrate by a distribution of Kondo couplings at the sites of the chain. Using auxiliary-field quantum Monte Carlo simulations, which quantitatively reproduce the field dependence of the zero-bias signal, we show that a proper Kondo resonance is present at the sites where the effective Kondo coupling dominates. Our modeling and numerical simulations provide a theoretical basis for the interpretation of the STM spectrum in terms of level crossings of the Co adatom chains.


title: The radiative transport equation with heterogeneous cross-sections

abstract: We consider the classical integral equation reformulation of the radiative transport equation (RTE) in a heterogeneous medium, assuming isotropic scattering. We prove an estimate for the norm of the integral operator in this formulation which is explicit in the (variable) coefficients of the problem (also known as the cross-sections). This result uses only elementary properties of the transport operator and some classical functional analysis. As a corollary, we obtain a bound on the convergence rate of source iteration (a classical stationary iterative method for solving the RTE). We also obtain an estimate for the solution of the RTE which is explicit in its dependence on the cross-sections. The latter can be used to estimate the solution in certain Bochner norms when the cross-sections are random fields. Finally we use our results to give an elementary proof that the generalised eigenvalue problem arising in nuclear reactor safety has only real and positive eigenvalues.


title: Fast learning synapses with molecular spin valves via selective magnetic potentiation

abstract: We studied LSMO/Alq3/AlOx/Co molecular spin valves in view of their use as synapses in neuromorphic computing. In neuromorphic computing, the learning ability is embodied in specific changes of the synaptic weight. In this perspective, the relevant parameter is the conductance of the molecular spin valve, which plays the role of the synaptic weight. In this work we demonstrated that the conductance can be changes by the repeated application of voltage pulses. We studied the parameter space of the pulses in order to determine the most effective voltage and duration of the pulses. The conductance could also be modified by aligning the magnetizations of the ferromagnetic electrodes parallel or anti parallel to each other. This phenomenon, known as magnetoresistance, affects high conductance devices while leaving low conductance devices unaffected. We studied how this weight update rule affected the speed of reward-based learning in an actor-critic framework, compared to a linear update rule. This nonlinear update performed significantly better (50 learning trials; Epochs to reach a performance goal of 0.975 was 896+/-301 in the nonlinear case and 1076+/-484 in the nonlinear case; Welch t-test: p<0.05). The linear update resulted in more learning trails with very long convergence times, which was largely absent in the nonlinear update.


title: Solovay reduction and continuity

abstract: The objective of this study is a better understanding of the relationships between reduction and continuity. Solovay reduction is a variation of Turing reduction based on the distance of two real numbers. We characterize Solovay reduction by the existence of a certain real function that is computable (in the sense of computable analysis) and Lipschitz continuous. We ask whether there exists a reducibility concept that corresponds to H\"older continuity. The answer is affirmative. We introduce quasi Solovay reduction and characterize this new reduction via H\"older continuity. In addition, we separate it from Solovay reduction and Turing reduction and investigate the relationships between complete sets and partial randomness.


title: A unified approach to symmetry for semilinear equations associated to the Laplacian in $\mathbb{R}^N$

abstract: We show radial symmetry of positive solutions to the H\'{e}non equation $-\Delta u = |x|^{-\ell} u^q $ in $\mathbb{R}^N \setminus \{ 0\} $, where $\ell \geq 0$, $q>0$ and satisfy further technical conditions. A new ingredient is a maximum principle for open subsets of a half space. It allows to apply the Moving Plane Method once a slow decay of the solution at infinity has been established, that is $\lim _{|x|\to \infty } |x|^{\gamma } u(x) =L $, for some numbers $\gamma \in (0, N-2)$ and $L >0$. Moreover, some examples of non-radial solutions are given for $q> \frac{N+1}{N-3}$ and $N\geq 4$. We also establish radial symmetry for related and more general problems in $\mathbb{R}^N $ and $\mathbb{R}^N \setminus \{ 0\} $.


title: High sensitivity VLBI with SKA

abstract: The Square Kilometre Array (SKA), with the aim of achieving a collecting area of one square kilometre, will be the world's largest radio telescope. A scientific collaboration between 12 countries (with more to join), it will consist of one Observatory with 2 telescopes located in South Africa and Australia. The telescope deployment is planned in two phases, but even in its first stage (SKA1) it will already enable transformational science in a broad range of scientific objectives. The inclusion of SKA1 in the Global VLBI networks (SKA-VLBI) will provide access to very high angular resolution to SKA science programmes in anticipation of the science to be realized with the full telescope deployment (SKA2). This contribution provides an overview of the SKA Observatory VLBI capability, the key operational concepts and outlines the need to update the science use cases.


title: Pushing Purcell-enhancement beyond its limits

abstract: Purcell-enhanced emission from a coupled emitter-cavity system is a fundamental manifestation of cavity quantum electrodynamics. Starting from a theoretical description we derive a scheme for photon emission from an emitter coupled to a birefringent cavity that exceeds hitherto anticipated limitations. Based on a recent study and experimental investigation of the intra-cavity coupling of orthogonal polarisation modes in birefringent cavities, we now decouple the emitter and the photon prior to emission from the cavity mode. Effectively, this is "hiding" the emitter from the photon in the cavity to suppress re-excitation, increasing the overall emission through the cavity mirrors. In doing so we show that tailored cavity birefringence can offer significant advantages and that these are practically achievable within the bounds of present-day technology. It is found that birefringence can mitigate the tradeoff between stronger emitter-cavity coupling and efficient photon extraction. This allows for longer cavities to be constructed without a loss of performance -- a significant result for applications where dielectric mirrors interfere with any trapping fields confining the emitter. We then generalise our model to consider a variety of equivalent schemes. For instance, detuning a pair of ground states in a three-level emitter coupled to a cavity in a Lambda-system is shown to provide the same enhancement, and it can be combined with a birefringent cavity to further increase performance. Additionally, it is found that when directly connecting multiple ground states of the emitter to form a chain of coupled states, the extraction efficiency approaches its fundamental upper limit. The principles proposed in this work can be applied in multiple ways to any emitter-cavity system, paving the way to surpassing the traditional limits of such systems with technologies that exist today.


title: Molecular beam epitaxy growth of superconducting Sr2RuO4 films

abstract: We report growth of superconducting Sr2RuO4 films by oxide molecular beam epitaxy (MBE). Careful tuning of the Ru flux with an electron beam evaporator enables us to optimize growth conditions including the Ru/Sr flux ratio and also to investigate stoichiometry effects on the structural and transport properties. The highest onset transition temperature of about 1.1 K is observed for films grown in a slightly Ru-rich flux condition in order to suppress Ru deficiency. The realization of superconducting Sr2RuO4 films via oxide MBE opens up a new route to study the unconventional superconductivity of this material.


title: Anomalous Enhancement of Upper Critical Field in Sr2RuO4 Thin Films

abstract: We report large enhancement of upper critical field Hc2 observed in superconducting Sr2RuO4 thin films. Through dimensional crossover approaching two dimensions, Hc2 except the in-plane field direction is dramatically enhanced compared to bulks, following a definite relation distinct from bulk one between Hc2 and the transition temperature. The anomalous enhancement of Hc2 is highly suggestive of important changes of the superconducting properties, possibly accompanied with rotation of the triplet d-vector. Our findings will become a crucial step to further explore exotic properties by employing Sr2RuO4 thin films.


title: Superfield continuous spin equations of motion

abstract: We propose a description of %manifestly supersymmetric continuous spin representations in $4D,\mathcal{N}=1$ Minkowski superspace at the level of equations of motions. The usual continuous spin wave function is promoted to a chiral or a complex linear superfield which includes both the single-valued (span integer helicities) and the double-valued (span half-integer helicities) representations thus making their connection under supersymmetry manifest. The set of proposed superspace constraints for both superfield generate the expected Wigner's conditions for both representations.


title: On Representability of Algebraic Functions by Radicals

abstract: This preprint is dedicated to a self contained simple proof of the classical criteria for representability of algebraic functions of several complex variables by radicals. It also contains a criteria for representability of algebroidal functions by composition of single-valued analytic functions and radicals, and a result related to the 13-th Hilbert problem. This preprint is an extended version of the author's 1971 paper. It is written as a part of the comments to a new edition (in preparation) of the classical book `Integration in finite terms' by J.F.~Ritt.


title: On limiting trace inequalities for vectorial differential operators

abstract: We establish that trace inequalities $$\|D^{k-1}u\|_{L^{\frac{n-s}{n-1}}(\mathbb{R}^{n},d\mu)} \leq c \|\mu\|_{L^{1,n-s}(\mathbb{R}^{n})}^{\frac{n-1}{n-s}}\|\mathbb{A}[D]u\|_{L^{1}(\mathbb{R}^{n},d\mathscr{L}^{n})}$$ hold for vector fields $u\in C^{\infty}(\mathbb{R}^{n};\mathbb{R}^{N})$ if and only if the $k$-th order homogeneous linear differential operator $\mathbb{A}[D]$ on $\mathbb{R}^{n}$ is elliptic and cancelling, provided that $s<1$, and give partial results for $s=1$, where stronger conditions on $\mathbb{A}[D]$ are necessary. Here, $\|\mu\|_{L^{1,\lambda}}$ denotes the $(1,\lambda)$-Morrey norm of the measure $\mu$, so that such traces can be taken, for example, with respect to the Hausdorff measure $\mathscr{H}^{n-s}$ restricted to fractals of codimension $0<s<1$. The above class of inequalities give a systematic generalisation of Adams' trace inequalities to the limit case $p=1$ and can be used to prove trace embeddings for functions of bounded $\mathbb{A}$-variation, thereby comprising Sobolev functions and functions of bounded variation or deformation. We moreover establish a multiplicative version of the above inequality, which implies ($\mathbb{A}$-)strict continuity of the associated trace operators on $\text{BV}^{\mathbb{A}}$.


title: Aggressive Local Search for Constrained Optimal Control Problems with Many Local Minima

abstract: This paper is concerned with numerically finding a global solution of constrained optimal control problems with many local minima. The focus is on the optimal decentralized control (ODC) problem, whose feasible set is recently shown to have an exponential number of connected components and consequently an exponential number of local minima. The rich literature of numerical algorithms for nonlinear optimization suggests that if a local search algorithm is initialized in an arbitrary connected component of the feasible set, it would search only within that component and find a stationary point there. This is based on the fact that numerical algorithms are designed to generate a sequence of points (via searching for descent directions and adjusting the step size), whose corresponding continuous path is trapped in a single connected component. In contrast with this perception rooted in convex optimization, we numerically illustrate that local search methods for non-convex constrained optimization can obliviously jump between different connected components to converge to a global minimum, via an aggressive step size adjustment using backtracking and the Armijio rule. To support the observations, we prove that from almost every arbitrary point in any connected component of the feasible set, it is possible to generate a sequence of points using local search to jump to different components and converge to a global solution. However, due to the NP-hardness of the problem, such fine-tuning of the parameters of a local search algorithm may need prior knowledge or be time consuming. This paper offers the first result on escaping non-global local solutions of constrained optimal control problems with complicated feasible sets.


title: Modelling and optimization of the excitonic diffraction grating

abstract: Periodical spatial modulation of the excitonic resonance in a quantum well could lead to the formation of a new highly directional and resonant coherent optical response -- resonant diffraction. Such excitonic diffraction gratings were demonstrated in epitaxially grown quantum wells patterned by the low-dose ion beam irradiation before or after the growth. In this paper we present a theoretical model of the resonant diffraction formation based on the step-by-step approximation of the Maxwell equation solution. The resulting theory allows us to reliably describe experimental data, as well as to predict the way to increase the diffraction efficiency.


title: An Efficient Schmidt-EKF for 3D Visual-Inertial SLAM

abstract: It holds great implications for practical applications to enable centimeter-accuracy positioning for mobile and wearable sensor systems. In this paper, we propose a novel, high-precision, efficient visual-inertial (VI)-SLAM algorithm, termed Schmidt-EKF VI-SLAM (SEVIS), which optimally fuses IMU measurements and monocular images in a tightly-coupled manner to provide 3D motion tracking with bounded error. In particular, we adapt the Schmidt Kalman filter formulation to selectively include informative features in the state vector while treating them as nuisance parameters (or Schmidt states) once they become matured. This change in modeling allows for significant computational savings by no longer needing to constantly update the Schmidt states (or their covariance), while still allowing the EKF to correctly account for their cross-correlations with the active states. As a result, we achieve linear computational complexity in terms of map size, instead of quadratic as in the standard SLAM systems. In order to fully exploit the map information to bound navigation drifts, we advocate efficient keyframe-aided 2D-to-2D feature matching to find reliable correspondences between current 2D visual measurements and 3D map features. The proposed SEVIS is extensively validated in both simulations and experiments.


title: Z_2-genus of graphs and minimum rank of partial symmetric matrices

abstract: The \emph{genus} $\mathrm{g}(G)$ of a graph $G$ is the minimum $g$ such that $G$ has an embedding on the orientable surface $M_g$ of genus $g$. A drawing of a graph on a surface is \emph{independently even} if every pair of nonadjacent edges in the drawing crosses an even number of times. The \emph{$\mathbb{Z}_2$-genus} of a graph $G$, denoted by $\mathrm{g}_0(G)$, is the minimum $g$ such that $G$ has an independently even drawing on $M_g$. By a result of Battle, Harary, Kodama and Youngs from 1962, the graph genus is additive over 2-connected blocks. In 2013, Schaefer and \v{S}tefankovi\v{c} proved that the $\mathbb{Z}_2$-genus of a graph is additive over 2-connected blocks as well, and asked whether this result can be extended to so-called 2-amalgamations, as an analogue of results by Decker, Glover, Huneke, and Stahl for the genus. We give the following partial answer. If $G=G_1\cup G_2$, $G_1$ and $G_2$ intersect in two vertices $u$ and $v$, and $G-u-v$ has $k$ connected components (among which we count the edge $uv$ if present), then $|\mathrm{g}_0(G)-(\mathrm{g}_0(G_1)+\mathrm{g}_0(G_2))|\le k+1$. For complete bipartite graphs $K_{m,n}$, with $n\ge m\ge 3$, we prove that $\frac{\mathrm{g}_0(K_{m,n})}{\mathrm{g}(K_{m,n})}=1-O(\frac{1}{n})$. Similar results are proved also for the Euler $\mathbb{Z}_2$-genus. We express the $\mathbb{Z}_2$-genus of a graph using the minimum rank of partial symmetric matrices over $\mathbb{Z}_2$; a problem that might be of independent interest.


title: $\pi$-fluxes, semi-metals and flat bands in artificial materials

abstract: The possibility of engineering experimentally viable systems that realize gauge fluxes within plaquettes of hopping have been subject of search for decades due to vast amounts of theoretical study. This is of particular interest for topological band insulators, where it is known that such fluxes bind protected mid-gap states. These modes can hybridize in extended flux lattices, giving rise to semi-metallic bands that are highly tunable. We demonstrate that within artificial materials, local $\pi$-fluxes can be naturally realized. Consequently, we provide concrete set-ups to access this physics and analyze similar self-organized band structures and physical properties. Our work therefore does not only pinpoint simple systems exhibiting flat bands, but also opens up a route to study flux lattice models and associated effective theories in a novel scene.


title: Partial regularity for symmetric quasiconvex functionals on BD

abstract: We establish the first partial regularity results for (strongly) symmetric quasiconvex functionals of linear growth on BD, the space of functions of bounded deformation. By Rindler's foundational work (Lower semicontinuity for integral functionals in the space of functions of bounded deformation via rigidity and Young measures, Arch. Ration. Mech. Anal. 202 (2011), no. 1, 63-113), symmetric quasiconvexity is the foremost notion as to sequential weak*-lower semicontinuity of functionals on BD. The overarching main difficulty here is Ornstein's Non-Inequality, implying that the BD-case is genuinely different from the study of variational integrals on BV. In particular, this paper extends the recent work of Kristensen and the author (Partial regularity for BV-Minimizers, Arch. Ration. Mech. Anal. 232 (2019), Issue 3, 1429-1473) from the BV- to the BD-situation. Alongside, we establish partial regularity results for strongly quasiconvex functionals of superlinear growth by reduction to the full gradient case, which might be of independent interest.


title: TATi-Thermodynamic Analytics ToolkIt: TensorFlow-based software for posterior sampling in machine learning applications

abstract: With the advent of GPU-assisted hardware and maturing high-efficiency software platforms such as TensorFlow and PyTorch, Bayesian posterior sampling for neural networks becomes plausible. In this article we discuss Bayesian parametrization in machine learning based on Markov Chain Monte Carlo methods, specifically discretized stochastic differential equations such as Langevin dynamics and extended system methods in which an ensemble of walkers is employed to enhance sampling. We provide a glimpse of the potential of the sampling-intensive approach by studying (and visualizing) the loss landscape of a neural network applied to the MNIST data set. Moreover, we investigate how the sampling efficiency itself can be significantly enhanced through an ensemble quasi-Newton preconditioning method. This article accompanies the release of a new TensorFlow software package, the Thermodynamic Analytics ToolkIt, which is used in the computational experiments.


title: An Optimal Task Allocation Strategy for Heterogeneous Multi-Robot Systems

abstract: For a team of heterogeneous robots executing multiple tasks, we propose a novel algorithm to optimally allocate tasks to robots while accounting for their different capabilities. Motivated by the need that robot teams have in many real-world applications of remaining operational for long periods of time, we allow each robot to choose tasks taking into account the energy consumed by executing them, besides the global specifications on the task allocation. The tasks are encoded as constraints in an energy minimization problem solved at each point in time by each robot. The prioritization of a task over others -- effectively signifying the allocation of the task to that particular robot -- occurs via the introduction of slack variables in the task constraints. Moreover, the suitabilities of certain robots towards certain tasks are also taken into account to generate a task allocation algorithm for a team of robots with heterogeneous capabilities. The efficacy of the developed approach is demonstrated both in simulation and on a team of real robots.


title: Photometric Mesh Optimization for Video-Aligned 3D Object Reconstruction

abstract: In this paper, we address the problem of 3D object mesh reconstruction from RGB videos. Our approach combines the best of multi-view geometric and data-driven methods for 3D reconstruction by optimizing object meshes for multi-view photometric consistency while constraining mesh deformations with a shape prior. We pose this as a piecewise image alignment problem for each mesh face projection. Our approach allows us to update shape parameters from the photometric error without any depth or mask information. Moreover, we show how to avoid a degeneracy of zero photometric gradients via rasterizing from a virtual viewpoint. We demonstrate 3D object mesh reconstruction results from both synthetic and real-world videos with our photometric mesh optimization, which is unachievable with either na\"ive mesh generation networks or traditional pipelines of surface reconstruction without heavy manual post-processing.


title: Online Gaussian Process State-Space Model: Learning and Planning for Partially Observable Dynamical Systems

abstract: This paper proposes an online learning method of Gaussian process state-space model (GP-SSM). GP-SSM is a probabilistic representation learning scheme that represents unknown state transition and/or measurement models as Gaussian processes (GPs). While the majority of prior literature on learning of GP-SSM are focused on processing a given set of time series data, data may arrive and accumulate sequentially over time in most dynamical systems. Storing all such sequential data and updating the model over entire data incur large amount of computational resources in space and time. To overcome this difficulty, we propose a practical method, termed \textit{onlineGPSSM}, that incorporates stochastic variational inference (VI) and online VI with novel formulation. The proposed method mitigates the computational complexity without catastrophic forgetting and also support adaptation to changes in a system and/or a real environments. Furthermore, we present application of onlineGPSSM into the reinforcement learning (RL) of partially observable dynamical systems by integrating onlineGPSSM with Bayesian filtering and trajectory optimization algorithms. Numerical examples are presented to demonstrate applicability of the proposed method.


title: Cross-study Reliability of the Open Card Sorting Method

abstract: Information architecture forms the foundation of users' navigation experience. Open card sorting is a widely-used method to create information architectures based on users' groupings of the content. However, little is known about the method's cross-study reliability: Does it produce consistent content groupings for similar profile participants involved in different card sort studies? This paper presents an empirical evaluation of the method's cross-study reliability. Six card sorts involving 140 participants were conducted: three open sorts for a travel website, and three for an eshop. Results showed that participants provided highly similar card sorting data for the same content. A rather high agreement of the produced navigation schemes was also found. These findings provide support for the cross-study reliability of the open card sorting method.


title: Entropies and their Asymptotic Theory in the discrete case

abstract: We present some new nonparametric estimators of entropies and we establish almost sure consistency and central limit Theorems for some of the most important entropies in the discrete case. Our theorical results are validated by simulations.


title: Bachet's game with lottery moves

abstract: Bachet's game is a variant of the game of Nim. There are $n$ objects in one pile. Two players take turns to remove any positive number of objects not exceeding some fixed number $m$. The player who takes the last object loses. We consider a variant of Bachet's game in which each move is a lottery over set $\{1,2,\ldots, m\}$. The outcome of a lottery is the number of objects that player takes from the pile. We show that under some nondegenericity assumptions on the set of available lotteries the probability that the first player wins in subgame perfect Nash equilibrium converges to $1/2$ as $n$ tends to infinity.


title: Numerical Algorithmic Science and Engineering within Computer Science: Rationale, Foundations and Organization

abstract: A re-calibration is proposed for "numerical analysis" as it arises specifically within the broader, embracing field of modern computer science (CS). This would facilitate research into theoretical and practicable models of real-number computation at the foundations of CS, and it would also advance the instructional objectives of the CS field. Our approach is premised on the key observation that the great "watershed" in numerical computation is much more between finite- and infinite-dimensional numerical problems than it is between discrete and continuous numerical problems. A revitalized discipline for numerical computation within modern CS can more accurately be defined as "numerical algorithmic science & engineering (NAS&E), or more compactly, as "numerical algorithmics," its focus being the algorithmic solution of numerical problems that are either discrete, or continuous over a space of finite dimension, or a combination of the two. It is the counterpart within modern CS of the numerical analysis discipline, whose primary focus is the algorithmic solution of continuous, infinite-dimensional numerical problems and their finite-dimensional approximates, and whose specialists today have largely been repatriated to departments of mathematics. Our detailed overview of NAS&E from the viewpoints of rationale, foundations, and organization is preceded by a recounting of the role played by numerical analysts in the evolution of academic departments of computer science, in order to provide background for NAS&E and place the newly-emerging discipline within its larger historical context.


title: Modelling Diffusion through Statistical Network Analysis: A Simulation Study

abstract: The study of international relations by definition deals with interdependencies among countries. One form of interdependence between countries is the diffusion of country-level features, such as policies, political regimes, or conflict. In these studies, the outcome variable tends to be categorical, and the primary concern is the clustering of the outcome variable among connected countries. Statistically, such clustering is studied with spatial econometric models. This paper instead proposes the use of a statistical network approach to model diffusion with a binary outcome variable. Using statistical network instead of spatial econometric models allows for a more natural specification of the diffusion process, assuming autocorrelation in the outcomes rather than the corresponding latent variable, and it simplifies the inclusion of temporal dynamics, higher level interdependencies and interactions between network ties and country-level features. In our simulations, the performance of the Stochastic Actor-Oriented Model (SAOM) estimator is evaluated. Our simulation results show that spatial parameters and coefficients on additional covariates in a static binary spatial autoregressive model are accurately recovered when using SAOM, albeit on a different scale. To demonstrate the use of this model, the paper applies the model to the international diffusion of same-sex marriage.


title: Face Detection in Repeated Settings

abstract: Face detection is an important first step before face verification and recognition. In unconstrained settings it is still an open challenge because of the variation in pose, lighting, scale, background and location. However, for the purposes of verification we can have a control on background and location. Images are primarily captured in places such as the entrance to a sensitive building, in front of a door or some location where the background does not change. We present a correlation based face detection algorithm to detect faces in such settings, where we control the location, and leave lighting, pose, and scale uncontrolled. In these scenarios the results indicate that our algorithm is easy and fast to train, outperforms Viola and Jones face detection accuracy and is faster to test.


title: Evolution of Dark Energy Perturbations for Slotheon Field and Power Spectrum

abstract: Within the framework of modified gravity model namely Slotheon model, inspired by the theory of extra dimensions, we explore the behaviour of Dark Energy and the perturbations thereof. The Dark Energy and matter perturbations equations are then derived and solved numerically by defining certain dimensionless variables and properly chosen initial conditions. The results are compared with those for standard quintessence model and $\Lambda$CDM model. The matter power spectrum is obtained and also compared with that for $\Lambda$CDM model. It appears that Dark Energy in Slotheon model is more akin to that for $\Lambda$CDM model than the standard quintessence model.


title: Numerical consideration of quasilinear electron cloud dynamics in plasma

abstract: The dynamics of a hot electron cloud in the solar corona-like plasma based on the numerical solution of kinetic equations of weak turbulence theory is considered. Different finite difference schemes are examined to fit the exact analytical solutions of quasilinear equations in hydrodynamic limit (gas-dynamic solution). It is shown that the scheme suggested demonstrates correct asymptotic behavior and can be employed to solve initial value problems for an arbitrary initial electron distribution function.

categories: physics.plasm-ph astro-ph.HE astro-ph.SR physics.space-ph


title: Learning Hierarchical Representations of Electronic Health Records for Clinical Outcome Prediction

abstract: Clinical outcome prediction based on the Electronic Health Record (EHR) plays a crucial role in improving the quality of healthcare. Conventional deep sequential models fail to capture the rich temporal patterns encoded in the longand irregular clinical event sequences. We make the observation that clinical events at a long time scale exhibit strongtemporal patterns, while events within a short time period tend to be disordered co-occurrence. We thus propose differentiated mechanisms to model clinical events at different time scales. Our model learns hierarchical representationsof event sequences, to adaptively distinguish between short-range and long-range events, and accurately capture coretemporal dependencies. Experimental results on real clinical data show that our model greatly improves over previous state-of-the-art models, achieving AUC scores of 0.94 and 0.90 for predicting death and ICU admission respectively, Our model also successfully identifies important events for different clinical outcome prediction tasks


title: Gravitational Waves, Extreme Astrophysics, and Fundamental Physics with Precision Pulsar Timing

abstract: Precision pulsar timing at the level of tens to hundreds of nanoseconds allows detection of nanohertz gravitational waves (GWs) from supermassive binary black holes (SMBBHs) at the cores of merging galaxies and, potentially, from exotic sources such as cosmic strings. The same timing data used for GW astronomy also yield precision masses of neutron stars orbiting other compact objects, constraints on the equation of state of nuclear matter, and precision tests of General Relativity, the Strong Equivalence Principle, and alternative theories of gravity. Timing can also lead to stringent constraints on the photon mass and on changes in fundamental constants and could reveal low mass objects (rogue planets, dark matter clumps) that traverse pulsar lines of sight. Data sets also allow modeling of the density, magnetic field, and turbulence in the interstellar plasma. Roughly 100 millisecond pulsars (MSPs) are currently being timed at $\sim$GHz frequencies using the largest radio telescopes in the world. These data sets currently represent ~1000 pulsar-years and will increase dramatically over the next decade. These topics are presented as a program of key science with flowdown technical requirements for achieving the science.


title: Velocity Dependent Dark Matter Interactions in Single-Electron Resolution Semiconductor Detectors with Directional Sensitivity

abstract: We investigate the velocity and recoil momentum dependence of dark matter interactions with ordinary matter. In particular we focus on the single-electron resolution semiconductor detectors, which allow experimental assessment of sub-GeV dark matter masses. We find that, within a specific mass range depending on the detector material, the dark matter interactions result in a signal characterized by daily modulation. Furthermore, we find that the detailed structure of this modulation is sensitive to the velocity and momentum dependence of dark matter interactions. We identify the optimal mass range for the prevalence of these effects.


title: The Compact Linear Collider (CLIC) - Project Implementation Plan

abstract: The Compact Linear Collider (CLIC) is a TeV-scale high-luminosity linear $e^+e^-$ collider under development by international collaborations hosted by CERN. This document provides an overview of the design, technology, and implementation aspects of the CLIC accelerator. For an optimal exploitation of its physics potential, CLIC is foreseen to be built and operated in stages, at centre-of-mass energies of 380 GeV, 1.5 TeV and 3 TeV, for a site length ranging between 11 km and 50 km. CLIC uses a Two-Beam acceleration scheme, in which normal-conducting high-gradient 12 GHz accelerating structures are powered via a high-current Drive Beam. For the first stage, an alternative with X-band klystron powering is also considered. CLIC accelerator optimisation, technical developments, and system tests have resulted insignificant progress in recent years. Moreover, this has led to an increased energy efficiency and reduced power consumption of around 170 MW for the 380 GeV stage, together with a reduced cost estimate of approximately 6 billion CHF. The construction of the first CLIC energy stage could start as early as 2026 and first beams would be available by 2035, marking the beginning of a physics programme spanning 25-30 years and providing excellent sensitivity to Beyond Standard Model physics, through direct searches and via a broad set of precision measurements of Standard Model processes, particularly in the Higgs and top-quark sectors.


title: Approximate Information Tests on Statistical Submanifolds

abstract: Parametric inference posits a statistical model that is a specified family of probability distributions. Restricted inference, e.g., restricted likelihood ratio testing, attempts to exploit the structure of a statistical submodel that is a subset of the specified family. We consider the problem of testing a simple hypothesis against alternatives from such a submodel. In the case of an unknown submodel, it is not clear how to realize the benefits of restricted inference. To do so, we first construct information tests that are locally asymptotically equivalent to likelihood ratio tests. Information tests are conceptually appealing but (in general) computationally intractable. However, unlike restricted likelihood ratio tests, restricted information tests can be approximated even when the statistical submodel is unknown. We construct approximate information tests using manifold learning procedures to extract information from samples of an unknown (or intractable) submodel, thereby providing a roadmap for computational solutions to a class of previously impenetrable problems in statistical inference. Examples illustrate the efficacy of the proposed methodology.


title: A Cosmic UV/X-ray Background Model Update

abstract: We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model (Faucher-Giguere et al. 2009), the new model provides a better match to a large number of up-to-date empirical constraints, including: 1) new galaxy and AGN luminosity functions; 2) stellar spectra including binary stars; 3) obscured and unobscured AGN; 4) a measurement of the non-ionizing UV background; 5) measurements of the intergalactic HI and HeII photoionization rates at z~0-6; 6) the local X-ray background; and 7) improved measurements of the intergalactic opacity. In this model, AGN dominate the HI ionizing background at z<~3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z~2, the slow evolution of the HI ionization rate inferred from the high-redshift HI Lya forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of ~1% suffices to ionize the intergalactic medium at z=3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the HeII Lya forest in hydrodynamic simulations.


title: Stability, reheating and leptogenesis

abstract: In a minimal model of leptogenesis, the observed baryon asymmetry is realized after high-scale reheating into the lightest sterile neutrino. We consider constraints on this scenario from the stability of the Higgs vacuum during pre-heating. Depending on the reheat temperature, the lightest sterile neutrino may be in or out of thermal equilibrium at production. Demanding stability of the Higgs vacuum during pre-heating, we find strong constraints which primarily impact the parameter space of thermal leptogenesis.


title: The evolution of the cosmic molecular gas density

abstract: One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.


title: Constraining integrable AdS/CFT with factorized scattering

abstract: We consider (warped) AdS string backgrounds which allow for a GKP spinning string/null cusp solution. Integrability implies that the worldsheet S-matrix should factorize, which in turn constrains the form of the warp factor as a function of the coordinates of the internal space. This constraint is argued to rule out integrability for all supersymmetric AdS7 and AdS6 backgrounds as well as AdS5 Gaiotto-Maldacena backgrounds and a few highly supersymmetric AdS4 and AdS3 backgrounds.


title: Observational constraints on the survival of pristine stars

abstract: There is a longstanding discussion about whether low mass stars can form from pristine gas in the early Universe. A particular point of interest is whether we can find surviving pristine stars from the first generation in our local neighbourhood. We present here a simple analytical estimate that puts tighter constraints on the existence of such stars. In the conventional picture, should these stars have formed in significant numbers and have preserved their pristine chemical composition until today, we should have found them already. With the presented method most current predictions for survivor counts larger than zero can be ruled out.


title: Imprints of Mass Accretion History on the Shape of the Intracluster Medium and the $T_X-M$ Relation

abstract: We use a statistical sample of galaxy clusters from a large cosmological $N$-body$+$hydrodynamics simulation to examine the relation between morphology, or shape, of the X-ray emitting intracluster medium (ICM) and the mass accretion history of the galaxy clusters. We find that the mass accretion rate (MAR) of a cluster is correlated with the ellipticity of the ICM. The correlation is largely driven by material accreted in the last $\sim 4.5$~Gyr, indicating a characteristic time-scale for relaxation of cluster gas. Furthermore, we find that the ellipticity of the outer regions ($R\sim R_{\rm 500c}$) of the ICM is correlated with the overall MAR of clusters, while ellipticity of the inner regions ($\lesssim 0.5 R_{\rm 500c}$) is sensitive to recent major mergers with mass ratios of $\geq 1:3$. Finally, we examine the impact of variations in cluster mass accretion history on the X-ray observable-mass scaling relations. We show that there is a {\it continuous\/} anti-correlation between the residuals in the $T_x-M$ relation and cluster MARs, within which merging and relaxed clusters occupy extremes of the distribution rather than form two peaks in a bi-modal distribution, as was often assumed previously. Our results indicate the systematic uncertainties in the X-ray observable-mass relations can be mitigated by using the information encoded in the apparent ICM ellipticity.


title: Entanglement negativity as a universal non-Markovianity witness

abstract: In order to engineer an open quantum system and its evolution, it is essential to identify and control the memory effects. These are formally attributed to the non-Markovianity of dynamics that manifests itself by the evolution being indivisible in time, a property which can be witnessed by a non-monotonic behavior of contractive functions or correlation measures. We show that by monitoring directly the entanglement behavior of a system in a tripartite setting it is possible to witness all invertible non-Markovian dynamics, as well as all (also non-invertible) qubit evolutions. This is achieved by using negativity, a computable measure of entanglement, which in the usual bipartite setting is not a universal non-Markovianity witness. We emphasize further the importance of multipartite states by showing that non-Markovianity cannot be faithfully witnessed by any contractive function of single qubits. We support our statements by an explicit example of eternally non-Markovian qubit dynamics, for which negativity can witness non-Markovianity at arbitrary time scales.


title: Massive Higher Spins: Effective Theory and Consistency

abstract: We construct the effective field theory for a single massive higher-spin particle in flat spacetime. Positivity bounds of the S-matrix force the cutoff of the theory to be well below the naive strong-coupling scale, forbid any potential and make therefore higher-derivative operators important even at low energy. As interesting application, we discuss in detail the massive spin-3 theory and show that an extended Galileon-like symmetry of the longitudinal modes, even with spin, emerges at high energy.


title: ASASSN-18tb: A Most Unusual Type Ia Supernova Observed by TESS and SALT

abstract: We present photometric and spectroscopic observations of the unusual Type Ia supernova ASASSN-18tb, including a series of SALT spectra obtained over the course of nearly six months and the first observations of a supernova by the Transiting Exoplanet Survey Satellite (TESS). We confirm a previous observation by Kollmeier et al. (2019) showing that ASASSN-18tb is the first relatively normal Type Ia supernova to exhibit clear broad ($\sim1000$ km s$^{-1}$) H$\alpha$ emission in its nebular phase spectra. We find that this event is best explained as a sub-Chandrasekhar mass explosion with $M_{Ni} \approx 0.3\; \rm{M}_\odot$. Despite the strong H$\alpha$ signature at late times, we find that the early rise of the supernova shows no evidence for deviations from a single-component power-law and is best fit with a moderately shallow power-law of index $1.69\pm0.04$. We find that the H$\alpha$ luminosity remains approximately constant after its initial detection at phase +37 d, and that the H$\alpha$ velocity evolution does not trace that of the Fe~III$~\lambda4660$ emission. These suggest that the H$\alpha$ emission arises from circumstellar medium (CSM) rather than swept up material from a non-degenerate companion. However, ASASSN-18tb is strikingly different from other known CSM-interacting Type Ia supernovae in a number of significant ways. Those objects typically show an H$\alpha$ luminosity two orders of magnitude higher than what is seen in ASASSN-18tb, pushing them away from the empirical light-curve relations that define "normal" Type Ia supernovae. Conversely, ASASSN-18tb exhibits a fairly typical light curve and luminosity for an underluminous or transitional SN Ia, with $M_R \approx -18.1$ mag. Moreover, ASASSN-18tb is the only SN Ia showing H$\alpha$ from CSM interaction to be discovered in an early-type galaxy.


title: Critical Viscosity of a Fluctuating Superconductor

abstract: We consider a fluctuating superconductor in the vicinity of the transition temperature, $T_c$. The fluctuation shear viscosity is calculated. In two dimensions, the leading correction to viscosity is negative and scales as $\delta \eta(T) \propto \ln(T-T_c)$. Critical hydrodynamics of the fluctuating superconductor involves two fluids -- a fluid of fluctuating pairs and a quasiparticle fluid of single-electron excitations. The pair viscosity (Aslamazov-Larkin) term is shown to be zero. The (density of states) correction to viscosity of single-electron excitations is negative, which is due to fluctuating pairing that results in a reduction of electron density. Scattering of electrons off of the fluctuations gives rise to an enhanced quasiparticle scattering and another (Maki-Thomson) negative correction to viscosity. Our results suggest that fluctuating superconductors provide a promising platform to investigate low-viscosity electronic media and may potentially host fermionic/electronic turbulence. Some experimental probes of two-fluid critical hydrodynamics are proposed such as time-of-flight measurement of turbulent energy cascades in critical cold atom superfluids and magnetic dynamos in three-dimensional fluctuating superconductors.


title: Enhanced Multi-Qubit Phase Estimation in Noisy Environments by Local Encoding

abstract: The first generation of multi-qubit quantum technologies will consist of noisy, intermediate-scale devices for which active error correction remains out of reach. To exploit such devices, it is thus imperative to use passive error protection that meets a careful trade-off between noise protection and resource overhead. Here, we experimentally demonstrate that single-qubit encoding can significantly enhance the robustness of entanglement and coherence of four-qubit graph states against local noise with a preferred direction. In particular, we explicitly show that local encoding provides a significant practical advantage for phase estimation in noisy environments. This demonstrates the efficacy of local unitary encoding under realistic conditions, with potential applications in multi-qubit quantum technologies for metrology, multi-partite secrecy and error correction.


title: Thermodynamics of Lorentzian Taub-NUT spacetimes

abstract: The thermodynamics of the Taub-NUT solution has been predominantly studied in the Euclidean sector, upon imposing the condition for the absence of Misner strings. Such thermodynamics is quite exceptional: the periodicity of the Euclidean time is restricted and thence the NUT charge cannot be independently varied, the entropy is not equal to a quarter of the area, and the thermodynamic volume can be negative. In this paper we revisit this paradigm and study the thermodynamics of the Lorentzian Taub-NUT solution, maintaining (as recently shown relatively harmless) Misner strings. We argue that in order to formulate a full cohomogeneity first law where the NUT parameter can be independently varied, it is natural to introduce a new charge together with its conjugate quantity. We consider two scenarios: one in which the entropy is given by the Lorentzian version of the Noether charge, the other in which the entropy is given by the standard Bekenstein--Hawking area law. In both cases consistent thermodynamics with positive thermodynamic volume can be formulated.


title: New Insights into the Physical Conditions and Internal Structure of a Candidate Proto-Globular Cluster

abstract: We present $\sim$0.1" resolution ($\sim$10 pc) ALMA observations of a molecular cloud identified in the merging Antennae galaxies with the potential to form a globular cluster, nicknamed the ``Firecracker.' Since star formation has not yet begun at an appreciable level, this cloud provides an example of what the birth environment of a globular cluster may have looked like before stars form and disrupt the natal conditions. Using emission from $^{12}$CO(2-1), $^{12}$CO(3-2), $^{13}$CO(2-1), HCN(4-3), and HCO$^+$(4-3), we are able to resolve the cloud's structure and find that it has a characteristic radius of 22 pc and a mass of 1--9$\times10^6 M_\odot$. We also constrain the abundance ratios of $^{12}$CO/$^{13}$CO and H$_2$/\twelveCO. Based on the calculated mass, we determine that the commonly used CO-to-H$_2$ conversion factor varies spatially, with average values in the range $X_{CO}=(0.12-1.1)\times10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We demonstrate that if the cloud is bound (as is circumstantially suggested by its bright, compact morphology), an external pressure of $P/k > 10^8$ K cm$^{-3}$ is required. This would be consistent with theoretical expectations that globular cluster formation requires high pressure environments. The position-velocity diagram of the cloud and its surrounding material suggests that this pressure may be produced by the collision of filaments. The radial profile of the column density can be fit with both a Gaussian and Bonnor-Ebert profile. The relative line strengths of HCN and HCO$^+$ in this region suggest that these molecular lines can be used as tracers for the evolutionary stage of a cluster.


title: Astro2020 Science White Paper: The Local Relics of of Supermassive Black Hole Seeds

abstract: We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.


title: Gradient based sample selection for online continual learning

abstract: A continual learning agent learns online with a non-stationary and never-ending stream of data. The key to such learning process is to overcome the catastrophic forgetting of previously seen data, which is a well known problem of neural networks. To prevent forgetting, a replay buffer is usually employed to store the previous data for the purpose of rehearsal. Previous works often depend on task boundary and i.i.d. assumptions to properly select samples for the replay buffer. In this work, we formulate sample selection as a constraint reduction problem based on the constrained optimization view of continual learning. The goal is to select a fixed subset of constraints that best approximate the feasible region defined by the original constraints. We show that it is equivalent to maximizing the diversity of samples in the replay buffer with parameters gradient as the feature. We further develop a greedy alternative that is cheap and efficient. The advantage of the proposed method is demonstrated by comparing to other alternatives under the continual learning setting. Further comparisons are made against state of the art methods that rely on task boundaries which show comparable or even better results for our method.


title: Vortex lattice in a uniform Bose-Einstein condensate in a box trap

abstract: We study numerically the vortex-lattice formation in a rapidly rotating uniform quasi-two-dimensional Bose-Einstein condensate (BEC) in a box trap. We consider two types of boxes: square and circle. In a square-shaped 2D box trap, when the number of generated vortices is the square of an integer, the vortices are found to be arranged in a perfect square lattice, although deviations near the center are found when the number of generated vortices is arbitrary. In case of a circular box trap, the generated vortices in the rapidly rotating BEC lie on concentric closed orbits. Near the center, these orbits have the shape of polygons, whereas near the periphery the orbits are circles. The circular box trap is equivalent to the rotating cylindrical bucket used in early experiment(s) with liquid He II. The number of generated vortices in both cases is in qualitative agreement with Feynman's universal estimate. The numerical simulation for this study is performed by a solution of the underlying mean-field Gross-Pitaevskii (GP) equation in the rotating frame, where the wave function for the generated vortex lattice is a stationary state. Consequently, the imaginary-time propagation method can be used for a solution of the GP equation, known to lead to an accurate numerical solution. We also demonstrated the dynamical stability of the vortex lattices in real-time propagation upon a small change of the angular frequency of rotation, using the converged imaginary-time wave function as the initial state.


title: Against "Reality" in Physics

abstract: The concept of "reality" is often raised in the context of philosophical foundations of physics or interpretations of quantum mechanics. When this term is so raised, it is a warning to me that I am about to be led down a rabbit hole. Such diversions usually lead nowhere unless you consider endless discussions of Schrodinger's cat, wave function collapse, quantum non-locality, and parallel universes to be useful. A prime example is the famous Einstein, Podolsky, and Rosen paper wherein they concluded that the quantum wave function cannot provide a complete description of physical reality. In this essay I suggest that, in physics discourse, the term "reality" should be avoided at all costs.


title: Performance Measurement for Deep Bayesian Neural Network

abstract: Deep Bayesian neural network has aroused a great attention in recent years since it combines the benefits of deep neural network and probability theory. Because of this, the network can make predictions and quantify the uncertainty of the predictions at the same time, which is important in many life-threatening areas. However, most of the recent researches are mainly focusing on making the Bayesian neural network easier to train, and proposing methods to estimate the uncertainty. I notice there are very few works that properly discuss the ways to measure the performance of the Bayesian neural network. Although accuracy and average uncertainty are commonly used for now, they are too general to provide any insight information about the model. In this paper, we would like to introduce more specific criteria and propose several metrics to measure the model performance from different perspectives, which include model calibration measurement, data rejection ability and uncertainty divergence for samples from the same and different distributions.


title: Discovery of polarizations from ground state absorption lines: tracer of sub-Gauss magnetic field on 89Her

abstract: We report the identification of the polarization of ground state absorption lines from post-AGB 89~Hercules. Two ground state neutral iron lines are found to have counterintuitive high-amplitude polarizations and an unchanged polarization direction through the orbital period, as opposed to the pattern of polarizations of absorption lines from excited states, which are synchronized with the orbital phase owing to optical pumping. This can be explained with magnetic realignment on the ground state. The 3D mean magnetic field is thereby unveiled from the degree and direction of the polarizations of the two iron lines. The field strength is also constrained to be $\lesssim 100\,{\rm mG}$. Our result has thus improved the accuracy by orders of magnitude compared to the previous $10\,{\rm G}$ upper limit set by non-detection of the Zeeman effect.


title: Parabolic Minkowski convolutions of viscosity solutions to fully nonlinear equations

abstract: This paper is concerned with the Minkowski convolution of viscosity solutions of fully nonlinear parabolic equations. We adopt this convolution to compare viscosity solutions of initial-boundary value problems in different domains. As a consequence, we can for instance obtain parabolic power concavity of solutions to a general class of parabolic equations. Our results apply to the Pucci operator, the normalized $q$-Laplacians with $1<q\leq\infty$, the Finsler Laplacian and more general quasilinear operators.