title: Appearance and disappearance of ferromagnetism in ultra-thin LaMnO$_3$ on SrTiO$_3$ substrate: a viewpoint from first-principles

abstract: The intrinsic magnetic state (ferromagnetic or antiferromagnetic) of ultra-thin LaMnO$_3$ films on the mostly used SrTiO$_3$ substrate is a long-existing question under debate. Either strain effect or non-stoichiometry was argued to be responsible for the experimental ferromagnetism. In a recent experiment [Science \textbf{349}, 716 (2015)], one more mechanism, namely the self-doping due to polar discontinuity, was argued to be the driving force of ferromagnetism beyond the critical thickness. Here systematic first-principles calculations have been performed to check these mechanisms in ultra-thin LaMnO$_3$ films as well as superlattices. Starting from the very precise descriptions of both LaMnO$_3$ and SrTiO$_3$, it is found that the compressive strain is the dominant force for the appearance of ferromagnetism, while the open surface with oxygen vacancies leads to the suppression of ferromagnetism. Within LaMnO$_3$ layers, the charge reconstructions involve many competitive factors and certainly go beyond the intuitive polar catastrophe model established for LaAlO$_3$/SrTiO$_3$ heterostructures. Our study not only explains the long-term puzzle regarding the magnetism of ultra-thin LaMnO$_3$ films, but also shed light on how to overcome the notorious magnetic dead layer in ultra-thin manganites.


title: Convergence of a $B$-$E$ based finite element method for MHD models on Lipschitz domains

abstract: We discuss a class of magnetic-electric fields based finite element schemes for stationary magnetohydrodynamics (MHD) systems with two types of boundary conditions. We establish a key $L^{3}$ estimate for divergence-free finite element functions for a new type of boundary conditions. With this estimate and a similar one in [Hu&Xu,2018], we rigorously prove the convergence of Picard iterations and the finite element schemes with weak regularity assumptions. These results demonstrate the convergence of the finite element methods for singular solutions.


title: Unidirectional Random Growth with Resetting

abstract: We review and classify stochastic processes without detailed balance condition. We obtain stationary distributions and investigate their stability in terms of generalized entropic divergences beyond the Kullback-Leibler formula. A simple stochastic model with local growth rates and direct resetting to the ground state is investigated and applied to various networks, scientific citations and Facebook popularity, hadronic yields in high energy particle reactions, income and wealth distributions, biodiversity and settlement size distribution.


title: Theory of time-averaged neutral dynamics with environmental stochasticity

abstract: Competition is the main driver of population dynamics, which shapes the genetic composition of populations and the assembly of ecological communities. Neutral models assume that all the individuals are equivalent and that the dynamics is governed by demographic (shot) noise, with a steady state species abundance distribution (SAD) that reflects a mutation-extinction equilibrium. Recently, many empirical and theoretical studies emphasized the importance of environmental variations that affect coherently the relative fitness of entire populations. Here we consider two generic time-averaged neutral models, in both the relative fitness of each species fluctuates independently in time but its mean is zero. The first (model A) describes a system with local competition and linear fitness-dependence of the birth-death rates, while in the second (model B) the competition is global and the fitness dependence is nonlinear. Due to this nonlinearity, model B admits a noise-induced stabilization mechanism that facilitates the invasion of new mutants. A self-consistent mean-field approach is used to reduce the multi-species problem to two-species dynamics, and the large-$N$ asymptotics of the emerging set of Fokker-Planck equations is presented and solved. Our analytic expressions are shown to fit the SADs obtained from extensive Monte-Carlo simulations and from numerical solutions of the corresponding master equations.


title: Radiative Dirac neutrino mass, DAMPE dark matter and leptogenesis

abstract: We explain the electron-positron excess reported by the DAMPE collaboration recently in a radiative Dirac seesaw model where a dark $U(1)_X$ gauge symmetry can (i) forbid the tree-level Yukawa couplings of three right-handed neutrinos to the standard model lepton and Higgs doublets, (ii) predict the existence of three dark fermions for the gauge anomaly cancellation, (iii) mediate a testable scattering of the lightest dark fermion off the nucleons. Our model can also accommodate a successful leptogenesis to generate the cosmic baryon asymmetry.


title: Domain Wall Motion Driven by Laplace Pressure in CoFeB-MgO Nanodots with Perpendicular Anisotropy

abstract: We have studied the magnetization reversal of CoFeB-MgO nanodots with perpendicular anisotropy for size ranging from w=400 nm to 1 {\mu}m. Contrary to previous experiments, the switching field distribution is shifted toward lower magnetic fields as the size of the elements is reduced with a mean switching field varying as 1/w. We show that this mechanism can be explained by the nucleation of a pinned magnetic domain wall (DW) at the edges of the nanodots where damages are introduced by the patterning process. As the surface tension (Laplace pressure) applied on the DW increases when reducing the size of the nanodots, we demonstrate that the depinning field to reverse the entire elements varies as 1/w. These results suggest that the presence of DWs has to be considered in the switching process of nanoscale elements and open a path toward scalable spintronic devices.


title: Quantifying layer similarity in multiplex networks: a systematic study

abstract: Computing layer similarities is an important way of characterizing multiplex networks because various static properties and dynamic processes depend on the relationships between layers. We provide a taxonomy and experimental evaluation of approaches to compare layers in multiplex networks. Our taxonomy includes, systematizes and extends existing approaches, and is complemented by a set of practical guidelines on how to apply them.


title: Element Distinctness Revisited

abstract: The element distinctness problem is the problem of determining whether the elements of a list are distinct, that is, if $x=(x_1,...,x_N)$ is a list with $N$ elements, we ask whether the elements of $x$ are distinct or not. The solution in a classical computer requires $N$ queries because it uses sorting to check whether there are equal elements. In the quantum case, it is possible to solve the problem in $O(N^{2/3})$ queries. There is an extension which asks whether there are $k$ colliding elements, known as element $k$-distinctness problem. This work obtains optimal values of two critical parameters of Ambainis' seminal quantum algorithm [SIAM J.~Comput., 37, 210-239, 2007]. The first critical parameter is the number of repetitions of the algorithm's main block, which inverts the phase of the marked elements and calls a subroutine. The second parameter is the number of quantum walk steps interlaced by oracle queries. We show that, when the optimal values of the parameters are used, the algorithm's success probability is $1-O(N^{1/(k+1)})$, quickly approaching 1. The specification of the exact running time and success probability is important in practical applications of this algorithm.


title: Enclosure of the Numerical Range and Resolvent Estimates of Non-selfadjoint Operator Functions

abstract: In this paper we discuss the relationship between the numerical range of an extensive class of unbounded operator functions and the joint numerical range of the operator coefficients. Furthermore, we derive methods on how to find estimates of the joint numerical range. Those estimates are used to obtain explicitly computable enclosures of the numerical range of the operator function and resolvent estimates. The enclosure and upper estimate of the norm of the resolvent are optimal given the estimate of the joint numerical range.


title: Fitting and selecting scattering data

abstract: The main purpose of scattering experiments is to unveil the underlying structure of the colliding particles and their interaction. Typically one measures scattering observables (cross sections and polarizations) at discrete angles and energies and mutually consistent data may validate or falsify proposed theories or models. However, the accumulation of data from different laboratories while potentially improves the statistical significance it may sometimes generate mutually inconsistent data as a side-effect. Thus, some decision has to be made on what are the maximal amount of data which are mutually compatible. We show elastic $\pi N $ and $NN$ scattering as prominent examples where this selection is called for. We discuss how it can be done in a self-consistent manner invoking a principle of maximal consensus of the database and with the help of a sufficiently flexible model involving a minimal number of theoretical assumptions. In the NN case this has become possible with a combination of long distance field theoretical constraints at the hadronic level such as pion exchanges and electromagnetic effects and a coarse graining of the unknown interaction over the shortest de Broglie wavelength being probed in the scattering process.


title: Radially-Distorted Conjugate Translations

abstract: This paper introduces the first minimal solvers that jointly solve for affine-rectification and radial lens distortion from coplanar repeated patterns. Even with imagery from moderately distorted lenses, plane rectification using the pinhole camera model is inaccurate or invalid. The proposed solvers incorporate lens distortion into the camera model and extend accurate rectification to wide-angle imagery, which is now common from consumer cameras. The solvers are derived from constraints induced by the conjugate translations of an imaged scene plane, which are integrated with the division model for radial lens distortion. The hidden-variable trick with ideal saturation is used to reformulate the constraints so that the solvers generated by the Grobner-basis method are stable, small and fast. Rectification and lens distortion are recovered from either one conjugately translated affine-covariant feature or two independently translated similarity-covariant features. The proposed solvers are used in a \RANSAC-based estimator, which gives accurate rectifications after few iterations. The proposed solvers are evaluated against the state-of-the-art and demonstrate significantly better rectifications on noisy measurements. Qualitative results on diverse imagery demonstrate high-accuracy undistortions and rectifications. The source code is publicly available at https://github.com/prittjam/repeats.


title: Nature of optical excitations in the frustrated kagome compound Herbertsmithite

abstract: Optical conductivity measurements are combined with density functional theory calculations in order to understand the electrodynamic response of the frustrated Mott insulators Herbertsmithite $\mathrm{ZnCu_{3}(OH)_{6}Cl_{2}}$ and the closely-related kagome-lattice compound $\mathrm{Y_{3}Cu_{9}(OH)_{19}Cl_{8}}$. We identify these materials as charge-transfer rather than Mott-Hubbard insulators, similar to the high-$T_c$ cuprate parent compounds. The band edge is at 3.3 and 3.6 eV, respectively, establishing the insulating nature of these compounds. Inside the gap, we observe dipole-forbidden local electronic transitions between the Cu $3d$ orbitals in the range 1--2 eV. With the help of \textit{ab initio} calculations we demonstrate that the electrodynamic response in these systems is directly related to the role of on-site Coulomb repulsion: while charge-transfer processes have their origin on transitions between the ligand band and the Cu $3d$ upper Hubbard band, \textit{local} $d$-$d$ excitations remain rather unaffected by correlations.


title: A rare and hidden attractor with noise in a biophysical Hodgkin-Huxley-type of model

abstract: In the present paper the result of numerical simulations of the model based on the Hodgkin-Huxley formalism with bistability between hidden and rare attractors in the presence of noise are studied.


title: Realizations of simple affine vertex algebras and their modules: the cases $\widehat{sl(2)}$ and $\widehat{osp(1,2)}$

abstract: We study embeddings of the simple admissible affine vertex algebras $V_k(sl(2))$ and $V_k(osp(1,2))$, $k \notin {\Bbb Z}_{\ge 0}$, into the tensor product of rational Virasoro and $N=1$ Neveu-Schwarz vertex algebra with lattice vertex algebras. We prove that the admissible affine vertex algebra $V_k(sl(2))$ can be embedded into vertex algebra $L^{Vir} (c_{p,p'}, 0) \otimes \Pi(0)$ where $L^{Vir} (c_{p,p'}, 0) $ is suitable minimal Virasoro vertex algebra and $\Pi(0)$ is a vertex algebra of lattice type. By using these realizations we construct a family of weight, logarithmic and Whittaker $\widehat{sl(2)}$ and $\widehat{osp(1,2)}$--modules. As an application, we construct all irreducible degenerate Whittaker modules for $V_k(sl(2))$.


title: Multivariate Time Series Classification with WEASEL+MUSE

abstract: Multivariate time series (MTS) arise when multiple interconnected sensors record data over time. Dealing with this high-dimensional data is challenging for every classifier for at least two aspects: First, an MTS is not only characterized by individual feature values, but also by the interplay of features in different dimensions. Second, this typically adds large amounts of irrelevant data and noise. We present our novel MTS classifier WEASEL+MUSE which addresses both challenges. WEASEL+MUSE builds a multivariate feature vector, first using a sliding-window approach applied to each dimension of the MTS, then extracts discrete features per window and dimension. The feature vector is subsequently fed through feature selection, removing non-discriminative features, and analysed by a machine learning classifier. The novelty of WEASEL+MUSE lies in its specific way of extracting and filtering multivariate features from MTS by encoding context information into each feature. Still the resulting feature set is small, yet very discriminative and useful for MTS classification. Based on a popular benchmark of 20 MTS datasets, we found that WEASEL+MUSE is among the most accurate classifiers, when compared to the state of the art. The outstanding robustness of WEASEL+MUSE is further confirmed based on motion gesture recognition data, where it out-of-the-box achieved similar accuracies as domain-specific methods.


title: First determination of $D^{*+}$-meson fragmentation functions and their uncertainties at next-to-next-to-leading order

abstract: We present, for the first time, a set of next-to-next-to-leading order (NNLO) fragmentation functions (FFs) describing the production of charmed-meson $D^{*+}$ from partons. Exploiting the universality and scaling violations of FFs, we extract the NLO and NNLO FFs through a global fit to all relevant data sets from single-inclusive $e^+e^-$ annihilation. The uncertainties for the resulting FFs as well as the corresponding observables are estimated using the Hessian approach. We evaluate the quality of the {\tt SKM18} FFs determined in this analysis by comparing with the recent results in literature and show how they describe the available data for single-inclusive $D^{*+}$-meson production in electron-positron annihilation. As a practical application, we apply the extracted FFs to make our theoretical predictions for the scaled-energy distributions of $D^{*+}$-mesons inclusively produced in top quark decays. We explore the implications of {\tt SKM18} for LHC phenomenology and show that our findings of this study can be introduced as a channel to indirect search for top-quark properties.


title: The signal of the stochastic gravitational wave background and the angular correlation of its energy density

abstract: The gravitational wave radiation emitted by all, resolved and unresolved, astrophysical sources in the observable universe generates a stochastic background. This background has a directional dependence inherited from the inhomogeneities of the matter distribution. This article proposes a new and independent derivation of the angular dependence of its energy density by focusing on the total gravitational wave signal produced by an ensemble of incoherent sources. This approach clarifies the origin of the angular correlation and the relation between the gravitational wave signal that can be measured by interferometers and the energy density of the stochastic background.


title: Quantum phase transition with dissipative frustration

abstract: We study the quantum phase transition of the one-dimensional phase model in the presence of dissipative frustration, provided by an interaction of the system with the environment through two non-commuting operators. Such a model can be realized in Josephson junction chains with shunt resistances and resistances between the chain and the ground. Using a self-consistent harmonic approximation, we determine the phase diagram at zero temperature which exhibits a quantum phase transition between an ordered phase, corresponding to the superconducting state, and a disordered phase, corresponding to the insulating state with localized superconducting charge. Interestingly, we find that the critical line separating the two phases has a non monotonic behavior as a function of the dissipative coupling strength. This result is a consequence of the frustration between (i) one dissipative coupling that quenches the quantum phase fluctuations favoring the ordered phase and (ii) one that quenches the quantum momentum (charge) fluctuations leading to a vanishing phase coherence. Moreover, within the self-consistent harmonic approximation, we analyze the dissipation induced crossover between a first and second order phase transition, showing that quantum frustration increases the range in which the phase transition is second order. The non monotonic behavior is reflected also in the purity of the system that quantifies the degree of correlation between the system and the environment, and in the logarithmic negativity as entanglement measure that encodes the internal quantum correlations in the chain.


title: Noninvasive Quantum Measurement of Arbitrary Operator Order by Engineered Non-Markovian Detectors

abstract: The development of solid-state quantum technologies requires the understanding of quantum measurements in interacting, non-isolated quantum systems. In general, a permanent coupling of detectors to a quantum system leads to memory effects that have to be taken into account in interpreting the measurement results. We analyze a generic setup of two detectors coupled to a quantum system and derive a compact formula in the weak-measurement limit that interpolates between an instantaneous (text-book type) and almost continuous - detector dynamics-dependent - measurement. A quantum memory effect that we term system-mediated detector-detector interaction is crucial to observe non-commuting observables simultaneously. Finally, we propose a mesoscopic double-dot detector setup in which the memory effect is tunable and that can be used to explore the transition to non-Markovian quantum measurements experimentally.


title: Local phase transitions in driven colloidal suspensions

abstract: Using dynamical density functional theory and Brownian dynamics simulations we investigate the influence of a driven tracer particle on the density distribution of a colloidal suspension at a thermodynamic statepoint close to the liquid side of the binodal. In bulk systems we find that a localized region of the colloid-poor phase, a 'cavitation bubble', forms behind the moving tracer. The extent of the cavitation bubble is investigated as a function of both the size and velocity of the tracer. The addition of a confining boundary enables us to investigate the interaction between the local phase instability at the substrate and that at the particle surface. When both the substrate and tracer interact repulsively with the colloids we observe the formation of a colloid-poor bridge between the substrate and the tracer. When a shear flow is applied parallel to the substrate the bridge becomes distorted and, at sufficiently high shear-rates, disconnects from the substrate to form a cavitation bubble.


title: Conformal Field Theory, Vertex Operator Algebras and Operator Algebras

abstract: We present recent progress in theory of local conformal nets which is an operator algebraic approach to study chiral conformal field theory. We emphasize representation theoretic aspects and relations to theory of vertex operator algebras which gives a different and algebraic formulation of chiral conformal field theory.

categories: math-ph math.CT math.MP math.OA math.QA math.RT


title: Static field limit of excitation probabilities in laser-atom interactions

abstract: We consider the interaction of atomic hydrogen, in its ground state, with an electromagnetic pulse whose duration is fixed in terms of the number of optical cycles. We study the probability of excitation of the atom in the static field limit i.e. for field frequencies going to zero. Despite the fact that the well known Born-Fock adiabatic theorem is valid only for a system whose energy spectrum is discrete, we show that it is still possible to use this theorem to derive, in the low frequency limit, an analytical formula which gives the probability of transition to any excited state of the atom as a function of the field intensity, the carrier envelope phase and the number of optical cycles within the pulse. The results for the probability of excitation to low-lying excited states, obtained with this formula, agree with those we get by solving the time-dependent Schroedinger equation. The domain of validity is discussed in detail.


title: A Stable SPH Discretization of the Elliptic Operator with Heterogeneous Coefficients

abstract: Smoothed particle hydrodynamics (SPH) has been extensively used to model high and low Reynolds number flows, free surface flows and collapse of dams, study pore-scale flow and dispersion, elasticity, and thermal problems. In different applications, it is required to have a stable and accurate discretization of the elliptic operator with homogeneous and heterogeneous coefficients. In this paper, the stability and approximation analysis of different SPH discretization schemes (traditional and new) of the diagonal elliptic operator for homogeneous and heterogeneous media are presented. The optimum and new discretization scheme is also proposed. This scheme enhances the Laplace approximation (Brookshaw's scheme (1985) and Schwaiger's scheme (2008)) used in the SPH community for thermal, viscous, and pressure projection problems with an isotropic elliptic operator. The numerical results are illustrated by numerical examples, where the comparison between different versions of the meshless discretization methods are presented.


title: Iron superhydrides FeH5 and FeH6: stability, electronic properties and superconductivity

abstract: Recently a big number of works devoted to search for new hydrides with record high-temperature superconductivity and at the same time the successful synthesis of potential high-TC superconducting FeH5 was reported. We present a systematic search for stable compounds in the Fe-H system using variable-composition version of the evolutionary algorithm USPEX. All known (FeH, FeH3, FeH5) and several new Fe3H5, Fe3H13 and FeH6 iron hydrides were found to be stable, resulting in a very complex phase diagram with rich structural relationships between phases. We calculate electronic properties of two potentially high-TC FeH5 and FeH6 phases in the pressure range from 150 to 300 GPa. Indeed, hydrogen-rich FeH5 and FeH6 phases were found to be superconducting within Bardeen-Cooper-Schrieffer theory, with TC values of up to 46 K.


title: Generalizing Gale's theorem on backward induction and domination of strategies

abstract: In 1953 Gale noticed that for every n-person game in extensive form with perfect information modeled by a rooted treesome special Nash equilibrium in pure strategies can be found by an algorithm of successive elimination of leaves, which is now called backward induction. He also noticed the same procedure, performed for the normal form of this game, turns into successive elimination of dominated strategies of the players that results in a single strategy profile (x_1,..., x_n), which is called a domination equilibrium (DE) and appears to be a Nash-equilibrium (NE) too. In other words, the game in normal form obtained from a positional game with perfect information is dominance-solvable (DS) and also Nash-solvable (NS). Yet, an arbitrary game in normal form may be not DS. We strengthen Gale's results as follows. Consider several successive eliminations of dominated strategies that begins with X = X_1 x ... x X_n and ends in X' = X'_1 x ... x X'_n. We will call X' a D-box of X. Our main (but obvious) lemma claims that for any i =1,..., n} and for any strategy x_i in X_i its projection to a D-box X' is dominated by a strategy x'_i in X'_i. It follows that any DE is an NE and, hence, DS implies NS. It is enough to apply the lemma in case when X' consists of a single strategy profile. The same lemma implies that the domination procedure is well-defined. A D-box X' is called terminal if it is domination-free, that is, it contains no pair of strategies such that one of them is dominated by the other. Any two terminal D-boxes X' and X" of X are equal. More precisely, there exist $n$ permutations \pi = (\pi_1, ..., \pi_n)$, with \pi_i : X_i to X_i for i in I, that transform X' into X", that is, \pi(X') = X" and the payoffs are respected. We also recall some published results on dominance-solvable game forms.


title: A limit field for orthogonal range searches in two-dimensional random point search trees

abstract: We consider the cost of general orthogonal range queries in random quadtrees. The cost of a given query is encoded into a (random) function of four variables which characterize the coordinates of two opposite corners of the query rectangle. We prove that, when suitably shifted and rescaled, the random cost function converges uniformly in probability towards a random field that is characterized as the unique solution to a distributional fixed-point equation. We also state similar results for $2$-d trees. Our results imply for instance that the worst case query satisfies the same asymptotic estimates as a typical query, and thereby resolve an old question of Chanzy, Devroye and Zamora-Cura [\emph{Acta Inf.}, 37:355--383, 2000]


title: Descent Representations of Generalized Coinvariant Algebras

abstract: The coinvariant algebra $R_n$ is a well-studied $\mathfrak{S}_n$-module that is a graded version of the regular representation of $\mathfrak{S}_n$. Using a straightening algorithm on monomials and the Garsia-Stanton basis, Adin, Brenti, and Roichman gave a description of the Frobenius image of $R_n$, graded by partitions, in terms of descents of standard Young tableaux. Motivated by the Delta Conjecture of Macdonald polynomials, Haglund, Rhoades, and Shimozono gave an extension of the coinvariant algebra $R_{n,k}$ and an extension of the Garsia-Stanton basis. Chan and Rhoades further extend these results from $\mathfrak{S}_n$ to the complex reflection group $G(r,1,n)$ by defining a $G(r,1,n)$ module $S_{n,k}$ that generalizes the coinvariant algebra for $G(r,1,n)$. We extend the results of Adin, Brenti, and Roichman to $R_{n,k}$ and $S_{n,k}$.


title: 3D Object Imaging through Scattering Media

abstract: Human ability to visualize an image is usually hindered by optical scattering. Recent extensive studies have promoted imaging technique through turbid materials to a reality where color image can be restored behind scattering media in real time. The big challenge now is to recover a 3D object in a large field of view with depth resolving ability. Here, we reveal a new physical relationship between speckles generated from objects at different planes. With a single given point spread function, 3D imaging through scattering media is achieved even beyond the depth of field (DOF). Experimental testing of standard scattering media shows that the original DOF can be extended up to 5 times and the physical mechanism is depicted. This extended 3D imaging is expected to have important applications in science, technology, bio-medical, security and defense.


title: Raga Identification using Repetitive Note Patterns from prescriptive notations of Carnatic Music

abstract: Carnatic music, a form of Indian Art Music, has relied on an oral tradition for transferring knowledge across several generations. Over the last two hundred years, the use of prescriptive notations has been adopted for learning, sight-playing and sight-singing. Prescriptive notations offer generic guidelines for a raga rendition and do not include information about the ornamentations or the gamakas, which are considered to be critical for characterizing a raga. In this paper, we show that prescriptive notations contain raga attributes and can reliably identify a raga of Carnatic music from its octave-folded prescriptive notations. We restrict the notations to 7 notes and suppress the finer note position information. A dictionary based approach captures the statistics of repetitive note patterns within a raga notation. The proposed stochastic models of repetitive note patterns (or SMRNP in short) obtained from raga notations of known compositions, outperforms the state of the art melody based raga identification technique on an equivalent melodic data corresponding to the same compositions. This in turn shows that for Carnatic music, the note transitions and movements have a greater role in defining the raga structure than the exact note positions.


title: Nanoladder cantilevers made from diamond and silicon

abstract: We present a "nanoladder" geometry that minimizes the mechanical dissipation of ultrasensitive cantilevers. A nanoladder cantilever consists of a lithographically patterned scaffold of rails and rungs with feature size $\sim$ 100 nm. Compared to a rectangular beam of the same dimensions, the mass and spring constant of a nanoladder are each reduced by roughly two orders of magnitude. We demonstrate a low force noise of $158 (+62)(-42)\,$zN and $190 (+42)(-33)\,$zN in a one-Hz bandwidth for devices made from silicon and diamond, respectively, measured at temperatures between 100--150 mK. As opposed to bottom-up mechanical resonators like nanowires or nanotubes, nanoladder cantilevers can be batch-fabricated using standard lithography, which is a critical factor for applications in scanning force microscopy.


title: Why So Many Published Sensitivity Analyses Are False. A Systematic Review of Sensitivity Analysis Practices

abstract: Sensitivity analysis (SA) has much to offer for a very large class of applications, such as model selection, calibration, optimization, quality assurance and many others. Sensitivity analysis offers crucial contextual information regarding a prediction by answering the question "Which uncertain input factors are responsible for the uncertainty in the prediction?" SA is distinct from uncertainty analysis (UA), which instead addresses the question "How uncertain is the prediction?" As we discuss in the present paper much confusion exists in the use of these terms. A proper uncertainty analysis of the output of a mathematical model needs to map what the model does when the input factors are left free to vary over their range of existence. A fortiori, this is true of a sensitivity analysis. Despite this, most UA and SA still explore the input space; moving along mono-dimensional corridors which leave the space of variation of the input factors mostly unscathed. We use results from a bibliometric analysis to show that many published SA fail the elementary requirement to properly explore the space of the input factors. The results, while discipline-dependent, point to a worrying lack of standards and of recognized good practices. The misuse of sensitivity analysis in mathematical modelling is at least as serious as the misuse of the p-test in statistical modelling. Mature methods have existed for about two decades to produce a defensible sensitivity analysis. We end by offering a rough guide for proper use of the methods.


title: Saturation model of DIS : an update

abstract: We present the results of new fits to the recently extracted data on $F_2$ at low $x$ with the GBW saturation model and its modification to cover high values of $Q^2$. We find that the model stands the test of time and gives a good description of the data with slightly modified parameters. All the essential elements of the model, especially the saturation scale, are retained.


title: A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

abstract: An explicit moving boundary method for the numerical solution of time-dependent hyperbolic conservation laws on grids produced by the intersection of complex geometries with a regular Cartesian grid is presented. As it employs directional operator splitting, implementation of the scheme is rather straightforward. Extending the method for static walls from Klein et al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme calculates fluxes needed for a conservative update of the near-wall cut-cells as linear combinations of standard fluxes from a one-dimensional extended stencil. Here the standard fluxes are those obtained without regard to the small sub-cell problem, and the linear combination weights involve detailed information regarding the cut-cell geometry. This linear combination of standard fluxes stabilizes the updates such that the time-step yielding marginal stability for arbitrarily small cut-cells is of the same order as that for regular cells. Moreover, it renders the approach compatible with a wide range of existing numerical flux-approximation methods. The scheme is extended here to time dependent rigid boundaries by reformulating the linear combination weights of the stabilizing flux stencil to account for the time dependence of cut-cell volume and interface area fractions. The two-dimensional tests discussed include advection in a channel oriented at an oblique angle to the Cartesian computational mesh, cylinders with circular and triangular cross-section passing through a stationary shock wave, a piston moving through an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil profile.


title: Finite-Range Corrections to the Thermodynamics of the One-Dimensional Bose Gas

abstract: The Lieb-Liniger equation of state accurately describes the zero-temperature universal properties of a dilute one-dimensional Bose gas in terms of the s-wave scattering length. For weakly-interacting bosons we derive non-universal corrections to this equation of state taking into account finite-range effects of the inter-atomic potential. Within the finite-temperature formalism of functional integration we find a beyond-mean-field equation of state which depends on scattering length and effective range of the interaction potential. Our analytical results, which are obtained performing dimensional regularization of divergent zero-point quantum fluctuations, show that for the one-dimensional Bose gas thermodynamic quantities like pressure and sound velocity are modified by changing the ratio between the effective range and the scattering length.


title: Infinite loop spaces and positive scalar curvature in the presence of a fundamental group

abstract: This is a continuation of our previous work with Botvinnik on the nontriviality of the secondary index invariant on spaces of metrics of positive scalar curvature, in which we take the fundamental group of the manifolds into account. We show that the secondary index invariant associated to the vanishing of the Rosenberg index can be highly nontrivial, for positive scalar curvature Spin manifolds with torsionfree fundamental groups which satisfy the Baum--Connes conjecture. For example, we produce a compact Spin 6-manifold such that its space of positive scalar curvature metrics has each rational homotopy group infinite dimensional. At a more technical level, we introduce the notion of "stable metrics" and prove a basic existence theorem for them, which generalises the Gromov--Lawson surgery technique, and we also give a method for rounding corners of manifold with positive scalar curvature metrics.


title: The massless integer superspin multiplets revisited

abstract: We propose a new off-shell formulation for the massless ${\cal N}=1$ supersymmetric multiplet of integer superspin $s$ in four dimensions, where $s =2,3,\dots$ (the $s=1$ case corresponds to the gravitino multiplet). Its gauge freedom matches that of the superconformal superspin-$s$ multiplet described in arXiv:1701.00682. The gauge-invariant action involves two compensating multiplets in addition to the superconformal superspin-$s$ multiplet. Upon imposing a partial gauge fixing, this action reduces to the one describing the so-called longitudinal formulation for the massless superspin-$s$ multiplet. Our new model is shown to possess a dual realisation obtained by applying a superfield Legendre transformation. We present a non-conformal higher spin supercurrent multiplet associated with the new integer superspin theory. This fermionic supercurrent is shown to occur in the Fayet-Sohnius model for a massive ${\cal N}=2$ hypermultiplet. We also give a new off-shell realisation for the massless gravitino multiplet.


title: A framework for on-line calibration of LINAC devices

abstract: General description of an on-line procedure of calibration for IGRT (Image Guided Radiotherapy) is given. The algorithm allows to improve targeting cancer by estimating its position in space and suggests appropriate correction of the position of the patient. The description is given in the Geometric Algebra language which significantly simplifies calculations and clarifies presentation.


title: Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence

abstract: Helicity is the scalar product between velocity and vorticity and, just like energy, its integral is an in-viscid invariant of the three-dimensional incompressible Navier-Stokes equations. However, space-and time-discretization methods typically corrupt this property, leading to violation of the inviscid conservation principles. This work investigates the discrete helicity conservation properties of spectral and finite-differencing methods, in relation to the form employed for the convective term. Effects due to Runge-Kutta time-advancement schemes are also taken into consideration in the analysis. The theoretical results are proved against inviscid numerical simulations, while a scale-dependent analysis of energy, helicity and their non-linear transfers is performed to further characterize the discretization errors of the different forms in forced helical turbulence simulations.


title: A radio spectral index map and catalogue at 147-1400 MHz covering 80 per cent of the sky

abstract: The radio spectral index is a powerful probe for classifying cosmic radio sources and understanding the origin of the radio emission. Combining data at 147 MHz and 1.4 GHz from the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS), we produced a large-area radio spectral index map of ~80 per cent of the sky (Dec > -40 deg), as well as a radio spectral index catalogue containing 1,396,515 sources, of which 503,647 are not upper or lower limits. Almost every TGSS source has a detected counterpart, while this is true only for 36 per cent of NVSS sources. We released both the map and the catalogue to the astronomical community. The catalogue is analysed to discover systematic behaviours in the cosmic radio population. We find a differential spectral behaviour between faint and bright sources as well as between compact and extended sources. These trends are explained in terms of radio galaxy evolution. We also confirm earlier reports of an excess of steep-spectrum sources along the galactic plane. This corresponds to 86 compact and steep-spectrum source in excess compared to expectations. The properties of this excess are consistent with normal non-recycled pulsars, which may have been missed by pulsation searches due to larger than average scattering along the line of sight.


title: Direct Segmented Sonification of Characteristic Features of the Data Domain

abstract: Sonification and audification create auditory displays of datasets. Audification translates data points into digital audio samples and the auditory display's duration is determined by the playback rate. Like audification, auditory graphs maintain the temporal relationships of data while using parameter mappings (typically data-to-frequency) to represent the ordinate values. Such direct approaches have the advantage of presenting the data stream `as is' without the imposed interpretations or accentuation of particular features found in indirect approaches. However, datasets can often be subdivided into short non-overlapping variable length segments that each encapsulate a discrete unit of domain-specific significant information and current direct approaches cannot represent these. We present Direct Segmented Sonification (DSSon) for highlighting the segments' data distributions as individual sonic events. Using domain knowledge to segment data, DSSon presents segments as discrete auditory gestalts while retaining the overall temporal regime and relationships of the dataset. The method's structural decoupling from the sound stream's formation means playback speed is independent of the individual sonic event durations, thereby offering highly flexible time compression/stretching to allow zooming into or out of the data. Demonstrated by three models applied to biomechanical data, DSSon displays high directness, letting the data `speak' for themselves.


title: On the normalized $p$-parabolic equation in arbitrary domains

abstract: The boundary regularity for the normalized $p$-parabolic equation $u_t =\frac{1}{p}|Du|^{2-p}\Delta_pu$ is studied. Perron's method is used to construct solutions in arbitrary domains. We classify the regular boundary points in terms of barrier functions, and prove an Exterior Sphere result. A fundamental solution is identified. A Petrovsky criterion is established, and we examine the convergence of solutions as $p \to \infty$.


title: Electron correlation effects and two-photon absorption in diamond shaped graphene quantum dots

abstract: In quasi-1D $\pi$-conjugated polymers such as \emph{trans}-polyacetylene and polyenes, electron correlation effects determine the "reversed" excited state ordering in which the lowest two-photon $2A_{g}$ state lies below the lowest one-photon $1B_{u}$ state. In this work, we present conclusive theoretical evidence of reversed excited state ordering in fairly 2D $\pi$-conjugated systems, namely, diamond-shaped graphene quantum dots (DQDs). Our electron correlated calculations show that DQDs begin to exhibit reversed excited ordering with increasing size, in disagreement with independent-particle picture. This signals the onset of strong correlation effects which renders them nonluminescent. Further, we calculate and analyze the two-photon absorption (TPA) spectra as well as photoinduced absorption (PA) spectra of these systems and find excellent agreement with the available experimental results. Our investigations demonstrate that unlike a strictly 1D system like\emph{trans}-polyacetylene, the non-linear and excited state absorptions in DQDs are highly intricate, with several even parity states responsible for strong absorptions. Our results could play an important role in the design of graphene-based non-linear optical devices.


title: Studies of mesic atoms and nuclei

abstract: $K^-$ mesons offer a unique setting where mesic atoms have been studied both experimentally and theoretically, thereby placing constraints on the possible existence and properties of meson-nuclear quasibound states. Here we review progress in this field made recently by the Jerusalem--Prague Collaboration using near-threshold $K^-N$ scattering amplitudes generated in several meson--baryon coupled channels models inspired by a chiral EFT approach. Our own procedure of handling subthreshold kinematics self consistently is used to transform these free-space energy dependent amplitudes to in-medium density dependent amplitudes from which $K^-$ optical potentials are derived. To fit the world data of kaonic atoms, these single-nucleon optical potentials are augmented by multi-nucleon terms. It is found that only two of the studied models reproduce also the single-nucleon absorption fractions available from old bubble chamber experiments. These two models are then checked for possible $K^-$ nuclear quasibound states, despite realizing that $K^-$ optical potentials are not constrained by kaonic atom data at densities exceeding half nuclear-matter density. We find that when such states exist, their widths are invariably above 100 MeV, forbiddingly large to allow observation. Multi-nucleon absorption is found to be substantial in this respect. This suggests that observable strongly bound $K^-$ mesons are limited to the very light systems, such as $K^-pp$.


title: Quantum Gates to other Universes

abstract: We present a microscopic model of a bridge connecting two large Anti-de-Sitter Universes. The Universes admit a holographic description as three-dimensional ${\cal N}=4$ supersymmetric gauge theories based on large linear quivers, and the bridge is a small rank-$n$ gauge group that acts as a messenger. On the gravity side, the bridge is a piece of a highly-curved AdS$_5\times$S$_5$ throat carrying $n$ units of five-form flux. We derive a universal expression for the mixing of the two massless gravitons: $M^2 \simeq 3n^2 (\kappa_4^2 + \kappa_4^{\prime\,2})/16\pi^2$, where $M$ is the mass splitting of the gravitons, $\kappa_4^2, \kappa_4^{\prime\,2}$ are the effective gravitational couplings of the AdS$_4$ Universes, and $n$ is the quantized charge of the gate. This agrees with earlier results based on double-trace deformations, with the important difference that the effective coupling is here quantized. We argue that the apparent non-localities of holographic double-trace models are resolved by integrating-in the (scarce) degrees of freedom of the gate.


title: Demonstration of a 2x2 programmable phase plate for electrons

abstract: First results on the experimental realisation of a 2x2 programmable phase plate for electrons are presented. The design consists of an array of electrostatic einzel lenses that influence the phase of electron waves passing through 4 separately controllable aperture holes. This functionality is demonstrated in a conventional transmission electron microscope operating at 300~kV and results are in very close agreement with theoretical predictions. The dynamic creation of a set of electron probes with different phase symmetry is demonstrated, thereby bringing adaptive optics in TEM one step closer to reality. The limitations of the current design and how to overcome these in the future are discussed. Simulations show how further evolved versions of the current proof of concept might open new and exciting application prospects for beam shaping and aberration correction.


title: Two Practical Methods for Coronal Intensity Determination

abstract: Determining the relative brightness of the solar corona is one of the most critical stages in solar eclipse studies. For this purpose, images taken with different exposures and polarization angles in white-light observations are used. The composite image of each polarization angle is produced by combining the images of different exposures. With the help of the intensity calibration function of these images, the relative intensity of the corona can be calculated. The total brightness of the solar corona is calculated using Stokes parameters obtained from intensity values of three polarization angles. In this study, two methods are presented: the first is used to obtain the intensity calibration function of the photographic material using calibration images, and the second is used to calculate the combined intensity values of images taken with different polarization angles.


title: Spin and Center of mass comparison between the PN approach and the asymptotic formulation

abstract: In this work we analyse the similarities and differences between the equations of motion for the center of mass and intrinsic angular momentum for isolated sources of gravitational radiation obtained by two different formulations. One approach is based on the asymptotic formulation of the GR whereas the other relies on Post- Newtonian methods. Several conclusions are obtained which could be useful for further developments in both approaches.


title: Model-independent analysis of the DAMPE excess

abstract: The Dark Matter Particle Explorer (DAMPE) recently released measurements of the electron spectrum with a hint of a narrow peak at about 1.4 TeV. We investigate dark matter (DM) models that could produce such a signal by annihilation in a nearby subhalo whilst simultaneously satisfying constraints from DM searches. In our model-independent approach, we consider all renormalizable interactions via a spin 0 or 1 mediator between spin 0 or 1/2 DM particles and the Standard Model leptons. We find that of the 20 combinations, 10 are ruled out by velocity or helicity suppression of the annihilation cross section to fermions. The remaining 10 models, though, evade constraints from the relic density, collider and direct detection searches, and include models of spin 0 and 1/2 DM coupling to a spin 0 or 1 mediator. We delineate the regions of mediator mass and couplings that could explain the DAMPE excess. In all cases the mediator is required to be heaver than about 2 TeV by LEP limits.


title: A tool for visualizing the execution of programs and stack traces especially suited for novice programmers

abstract: Software engineering education and training have obstacles caused by a lack of basic knowledge about a process of program execution. The article is devoted to the development of special tools that help to visualize the process. We analyze existing tools and propose a new approach to stack and heap visualization. The solution is able to overcome major drawbacks of existing tools and suites well for analysis of programs written in Java and C/C++.


title: Elastohydrodynamic lift at a soft wall

abstract: We study experimentally the motion of non-deformable microbeads in a linear shear flow close to a wall bearing a thin and soft polymer layer. Combining microfluidics and 3D optical tracking, we demonstrate that the steady-state bead/surface distance increases with the flow strength. Moreover, such lift is shown to result from flow-induced deformations of the layer, in quantitative agreement with theoretical predictions from elastohydrodynamics. This study thus provides the first experimental evidence of "soft lubrication" at play at small scale, in a system relevant {\it e.g.} to the physics of blood microcirculation.


title: 3DContextNet: K-d Tree Guided Hierarchical Learning of Point Clouds Using Local and Global Contextual Cues

abstract: Classification and segmentation of 3D point clouds are important tasks in computer vision. Because of the irregular nature of point clouds, most of the existing methods convert point clouds into regular 3D voxel grids before they are used as input for ConvNets. Unfortunately, voxel representations are highly insensitive to the geometrical nature of 3D data. More recent methods encode point clouds to higher dimensional features to cover the global 3D space. However, these models are not able to sufficiently capture the local structures of point clouds. Therefore, in this paper, we propose a method that exploits both local and global contextual cues imposed by the k-d tree. The method is designed to learn representation vectors progressively along the tree structure. Experiments on challenging benchmarks show that the proposed model provides discriminative point set features. For the task of 3D scene semantic segmentation, our method significantly outperforms the state-of-the-art on the Stanford Large-Scale 3D Indoor Spaces Dataset (S3DIS).


title: On the proof of the $C^0$-inextendibility of the Schwarzschild spacetime

abstract: This article presents a streamlined version of the author's original proof of the $C^0$-inextendibility of the maximal analytic Schwarzschild spacetime. Firstly, we deviate from the original proof by using the result, recently established in collaboration with Galloway and Ling, that given a $C^0$-extension of a globally hyperbolic spacetime, one can find a timelike geodesic that leaves this spacetime. This result much simplifies the proof of the inextendibility through the exterior region of the Schwarzschild spacetime. Secondly, we give a more flexible and shorter argument for the inextendibility through the interior region. Furthermore, we present a small new structural result for the boundary of a globally hyperbolic spacetime within a $C^0$-extension which serves as a new and simpler starting point for the proof.


title: Asymptotic structure of spacetime and the Newman-Penrose formalism: a brief review

abstract: A brief review about the Newman-Penrose formalism and the asymptotic structure of the spacetime is given. The goal of this review is to describe the latest developments in these topics and make a summary of the most important articles published by Newman and collaborators. Additionally, we discuss some aspects of this approach, and we compute the spin coefficients and the Weyl scalars for a general stationary axisymmetric spacetimes in a tetrad basis different from that defined by the principal null geodesic directions.


title: Short-distance breakdown of the Higgs mechanism and the robustness of the BCS theory for charged superconductors

abstract: Through the Higgs mechanism, the long-range Coulomb interaction eliminates the low-energy Goldstone phase mode in superconductors and transfers spectral weight all the way up to the plasma frequency. Here we show that the Higgs mechanism breaks down for length scales shorter than the superconducting coherence length while it stays intact, even at high energies, in the long-wavelength limit. This effect is a consequence of the composite nature of the Higgs field of superconductivity and the broken Lorentz invariance in a solid. Most importantly, the breakdown of the Higgs mechanism inside the superconducting coherence volume is crucial to ensure the stability of the BCS mean-field theory in the weak-coupling limit. We also show that changes in the gap equation due to plasmon-induced fluctuations can lead to significant corrections to the mean-field theory and reveal that changes in the density-fluctuation spectrum of a superconductor are not limited to the vicinity of the gap.


title: Learning to Learn from Weak Supervision by Full Supervision

abstract: In this paper, we propose a method for training neural networks when we have a large set of data with weak labels and a small amount of data with true labels. In our proposed model, we train two neural networks: a target network, the learner and a confidence network, the meta-learner. The target network is optimized to perform a given task and is trained using a large set of unlabeled data that are weakly annotated. We propose to control the magnitude of the gradient updates to the target network using the scores provided by the second confidence network, which is trained on a small amount of supervised data. Thus we avoid that the weight updates computed from noisy labels harm the quality of the target network model.


title: Sources of gravitational waves in asymptotically flat Einstein-Maxwell spacetime

abstract: In this work, the dynamic of isolated systems in general relativity is described when gravitational radiation and electromagnetic fields are present. In this construction, the asymptotic fields received at null infinity together with the regularized null cone cuts equation, and the center of mass of an asymptotically flat Einstein-Maxwell spacetime are used. A set of equations are derived in the low speed regime, linking their time evolution to the emitted gravitational radiation and to the Maxwell fields received at infinity. These equations should be useful when describing the dynamic of compact sources, such as the final moments of binary coalescence and the evolution of the final black hole. Additionally, we compare our equations with those coming from a similar approach given by Newman, finding some differences in the motion of the center of mass and spin of the gravitational system.


title: Computing the nucleon charge and axial radii directly at $Q^2=0$ in lattice QCD

abstract: We describe a procedure for extracting momentum derivatives of nucleon matrix elements on the lattice directly at $Q^2=0$. This is based on the Rome method for computing momentum derivatives of quark propagators. We apply this procedure to extract the nucleon isovector magnetic moment and charge radius as well as the isovector induced pseudoscalar form factor at $Q^2=0$ and the axial radius. For comparison, we also determine these quantities with the traditional approach of computing the corresponding form factors, i.e. $G^v_E(Q^2)$ and $G_M^v(Q^2)$ for the case of the vector current and $G_P^v(Q^2)$ and $G_A^v(Q^2)$ for the axial current, at multiple $Q^2$ values followed by $z$-expansion fits. We perform our calculations at the physical pion mass using a 2HEX-smeared Wilson-clover action. To control the effects of excited-state contamination, the calculations were done at three source-sink separations and the summation method was used. The derivative method produces results consistent with those from the traditional approach but with larger statistical uncertainties especially for the isovector charge and axial radii.


title: MR image reconstruction using deep density priors

abstract: Algorithms for Magnetic Resonance (MR) image reconstruction from undersampled measurements exploit prior information to compensate for missing k-space data. Deep learning (DL) provides a powerful framework for extracting such information from existing image datasets, through learning, and then using it for reconstruction. Leveraging this, recent methods employed DL to learn mappings from undersampled to fully sampled images using paired datasets, including undersampled and corresponding fully sampled images, integrating prior knowledge implicitly. In this article, we propose an alternative approach that learns the probability distribution of fully sampled MR images using unsupervised DL, specifically Variational Autoencoders (VAE), and use this as an explicit prior term in reconstruction, completely decoupling the encoding operation from the prior. The resulting reconstruction algorithm enjoys a powerful image prior to compensate for missing k-space data without requiring paired datasets for training nor being prone to associated sensitivities, such as deviations in undersampling patterns used in training and test time or coil settings. We evaluated the proposed method with T1 weighted images from a publicly available dataset, multi-coil complex images acquired from healthy volunteers (N=8) and images with white matter lesions. The proposed algorithm, using the VAE prior, produced visually high quality reconstructions and achieved low RMSE values, outperforming most of the alternative methods on the same dataset. On multi-coil complex data, the algorithm yielded accurate magnitude and phase reconstruction results. In the experiments on images with white matter lesions, the method faithfully reconstructed the lesions. Keywords: Reconstruction, MRI, prior probability, machine learning, deep learning, unsupervised learning, density estimation


title: Domain imaging across the magneto-structural phase transition in Fe$_{1+y}$Te

abstract: The investigation of the magnetic phase transitions in the parent compounds of Fe-based superconductors is regarded essential for an understanding of the pairing mechanism in the related superconducting compounds. Even though the chemical and electronic properties of these materials are often strongly inhomogeneous on a nanometer length scale, studies of the magnetic phase transitions using spatially resolved experimental techniques are still scarce. Here, we present a real space spin-resolved scanning tunneling microscopy investigation of the surface of Fe$_{1+y}$Te single crystals with different excess Fe content, $y$, which are continuously driven through the magnetic phase transition. For Fe$_{1.08}$Te, the transition into the low-temperature monoclinic commensurate antiferromagnetic phase is accompanied by the sudden emergence of ordering into four rotational domains with different orientations of the monoclinic lattice and of the antiferromagnetic order, showing how structural and magnetic order are intertwined. In the low-temperature phase of Fe$_{1.12}$Te one type of the domain boundaries disappears, and the transition into the paramagnetic phase gets rather broad, which is assigned to the formation of a mixture of orthorhombic and monoclinic phases.


title: Heating and cooling of coronal loops with turbulent suppression of parallel heat conduction

abstract: Using the "enthalpy-based thermal evolution of loops" (EBTEL) model, we investigate the hydrodynamics of the plasma in a flaring coronal loop in which heat conduction is limited by turbulent scattering of the electrons that transport the thermal heat flux. The EBTEL equations are solved analytically in each of the two (conduction-dominated and radiation-dominated) cooling phases. Comparison of the results with typical observed cooling times in solar flares shows that the turbulent mean free-path $\lambda_T$ lies in a range corresponding to a regime in which classical (collision-dominated) conduction plays at most a limited role. We also consider the magnitude and duration of the heat input that is necessary to account for the enhanced values of temperature and density at the beginning of the cooling phase and for the observed cooling times. We find through numerical modeling that in order to produce a peak temperature $\simeq 1.5 \times 10^7$~K and a 200~s cooling time consistent with observations, the flare heating profile must extend over a significant period of time; in particular, its lingering role must be taken into consideration in any description of the cooling phase. Comparison with observationally-inferred values of post-flare loop temperatures, densities, and cooling times thus leads to useful constraints on both the magnitude and duration of the magnetic energy release in the loop, as well as on the value of the turbulent mean free-path $\lambda_T$.

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


title: Non-Left-Orderable Surgeries on 1-Bridge Braids

abstract: Boyer, Gordon, and Watson have conjectured that an irreducible rational homology 3-sphere is an L-space if and only if its fundamental group is not left-orderable. Since Dehn surgeries on knots in $S^3$ can produce large families of L-spaces, it is natural to examine the conjecture on these 3-manifolds. Greene, Lewallen, and Vafaee have proved that all 1-bridge braids are L-space knots. In this paper, we consider three families of 1-bridge braids. First we calculate the knot groups and peripheral subgroups. We then verify the conjecture on the three cases by applying the criterion developed by Christianson, Goluboff, Hamann, and Varadaraj, when they verified the same conjecture for certain twisted torus knots and generalized the criteria of Clay and Watson and of Ichihara and Temma.


title: Semi-Distributed Demand Response Solutions for Smart Homes

abstract: The Internet of Things (IoT) paradigm brings an opportunity for advanced Demand Response (DR) solutions. It enables visibility and control on the various appliances that may consume, store or generate energy within a home. It has been shown that a centralized control on the appliances of a set of households leads to efficient DR mechanisms; unfortunately, such solutions raise privacy and scalability issues. In this chapter we propose an approach that deals with these issues. Specifically, we introduce a scalable two-levels control system where a centralized controller allocates power to each house on one side and, each household implements a DR local solution on the other side. A limited feedback to the centralized controller allows to enhance the performance with little impact on privacy. The solution is proposed for the general framework of capacity markets.


title: Roadmap of left-right models based on GUTs

abstract: We perform a detailed study of the grand unified theories $SO(10)$ and $E(6)$ with left-right intermediate gauge symmetries of the form $SU(N)_L\otimes SU(N)_R \otimes \mathcal{G}$. Proton decay lifetime constrains the unification scale to be $ \gtrsim 10^{16}$ GeV and, as discussed in this paper, unwanted cosmological relics can be evaded if the intermediate symmetry scale is $ \gtrsim 10^{12}$ GeV. With these conditions, we study the renormalization group evolution of the gauge couplings and do a comparative analysis of all possible left-right models where unification can occur. Both the D-parity conserved and broken scenarios, as well as the supersymmetric (SUSY) and Non-supersymmetric (Non-SUSY) versions, are considered. In addition to the fermion and scalar representations at each stage of the symmetry breaking, contributing to the $\beta$-functions, we list the intermediate left-right groups which successfully meet these requirements. We make use of the dimension-5 kinetic mixing effective operators for achieving unification and large intermediate scale. A significant result in the supersymmetric case is that to achieve successful unification for some breaking patterns, the scale of SUSY breaking needs to be at least a few TeV. In some of these cases, the intermediate scale can be as low as $\sim 10^{12}$ GeV, for SUSY scale to be $\sim 30$ TeV. This has important consequences in the collider searches for SUSY particles and phenomenology of the lightest neutralino as dark matter.


title: Two-sided Facility Location

abstract: Recent years have witnessed the rise of many successful e-commerce marketplace platforms like the Amazon marketplace, AirBnB, Uber/Lyft, and Upwork, where a central platform mediates economic transactions between buyers and sellers. Motivated by these platforms, we formulate a set of facility location problems that we term Two-sided Facility location. In our model, agents arrive at nodes in an underlying metric space, where the metric distance between any buyer and seller captures the quality of the corresponding match. The platform posts prices and wages at the nodes, and opens a set of facilities to route the agents to. The agents at any facility are assumed to be matched. The platform ensures high match quality by imposing a distance constraint between a node and the facilities it is routed to. It ensures high service availability by ensuring flow to the facility is at least a pre-specified lower bound. Subject to these constraints, the goal of the platform is to maximize the social surplus (or gains from trade) subject to weak budget balance, i.e., profit being non-negative. We present an approximation algorithm for this problem that yields a $(1 + \epsilon)$ approximation to surplus for any constant $\epsilon > 0$, while relaxing the match quality (i.e., maximum distance of any match) by a constant factor. We use an LP rounding framework that easily extends to other objectives such as maximizing volume of trade or profit. We justify our models by considering a dynamic marketplace setting where agents arrive according to a stochastic process and have finite patience (or deadlines) for being matched. We perform queueing analysis to show that for policies that route agents to facilities and match them, ensuring a low abandonment probability of agents reduces to ensuring sufficient flow arrives at each facility.


title: Light yield and radiation hardness studies of scintillator strips with a filler

abstract: Detectors based on polystyrene scintillator strips with WLS fiber readout are widely used to register charged particles in many high-energy physics experiments. The fibers are placed into grooves or holes along the strip. The detection efficiency of these devices can be significantly increased by improving the optical contact between the scintillator and the fiber by adding an optical filler into the groove/hole. This work is devoted to the study of the light yield of a 5-m-long scintillator strip with a 1.2-mm-diameter Kuraray Y11(200)~MC WLS fiber inserted into the strip's co-extruded hole filled with synthetic silicon resin SKTN-MED(E). The light yield was studied using cosmic muons and a $^{60}$Co radioactive source. Radiation hardness study of viscous fillers and short strip samples were performed on the IBR-2 pulsed research reactor of fast neutrons at JINR.


title: Who wins the Miss Contest for Imputation Methods? Our Vote for Miss BooPF

abstract: Missing data is an expected issue when large amounts of data is collected, and several imputation techniques have been proposed to tackle this problem. Beneath classical approaches such as MICE, the application of Machine Learning techniques is tempting. Here, the recently proposed missForest imputation method has shown high imputation accuracy under the Missing (Completely) at Random scheme with various missing rates. In its core, it is based on a random forest for classification and regression, respectively. In this paper we study whether this approach can even be enhanced by other methods such as the stochastic gradient tree boosting method, the C5.0 algorithm or modified random forest procedures. In particular, other resampling strategies within the random forest protocol are suggested. In an extensive simulation study, we analyze their performances for continuous, categorical as well as mixed-type data. Therein, MissBooPF, a combination of the stochastic gradient tree boosting method together with the parametrically bootstrapped random forest method, appeared to be promising. Finally, an empirical analysis focusing on credit information and Facebook data is conducted.


title: Causal horizons in a bouncing universe

abstract: As our understanding of the past in a bouncing universe is limited, it becomes difficult to propose a cosmological model which can give some understanding of the causal structure of the bouncing universe. In this article we address the issue related to the particle horizon problem in the bouncing universe models. It is shown that in many models the particle horizon does not exist, and consequently the horizon problem is trivially solved. In some cases a bouncing universe can have a particle horizon and we specify the conditions for its existence. In the absence of a particle horizon the Hubble surface specifies the causal structure of a bouncing universe. We specify the complex relationship between the Hubble surface and the particle horizon when the particle horizon exists. The article also address the issue related to the event horizon in a bouncing universe. A toy example of a bouncing universe is first presented where we specify the conditions which dictate the presence of a particle horizon. Next we specify the causal structures of three widely used bouncing models. The first case is related to quintom matter bounce model, the second one is loop quantum cosmology based bounce model and lastly $f(R)$ gravity induced bounce model. We present a brief discussion on the horizon problem in bouncing cosmologies. We point out that the causal structure of the various bounce models fit our general theoretical predictions.


title: CONFLLVM: A Compiler for Enforcing Data Confidentiality in Low-Level Code

abstract: We present an instrumenting compiler for enforcing data confidentiality in low-level applications (e.g. those written in C) in the presence of an active adversary. In our approach, the programmer marks secret data by writing lightweight annotations on top-level definitions in the source code. The compiler then uses a static flow analysis coupled with efficient runtime instrumentation, a custom memory layout, and custom control-flow integrity checks to prevent data leaks even in the presence of low-level attacks. We have implemented our scheme as part of the LLVM compiler. We evaluate it on the SPEC micro-benchmarks for performance, and on larger, real-world applications (including OpenLDAP, which is around 300KLoC) for programmer overhead required to restructure the application when protecting the sensitive data such as passwords. We find that performance overheads introduced by our instrumentation are moderate (average 12% on SPEC), and the programmer effort to port OpenLDAP is only about 160 LoC.


title: Predicting the hypervelocity star population in Gaia

abstract: Hypervelocity stars (HVSs) are amongst the fastest objects in our Milky Way. These stars are predicted to come from the Galactic center (GC) and travel along unbound orbits across the Galaxy. In the coming years, the ESA satellite Gaia will provide the most complete and accurate catalogue of the Milky Way, with full astrometric parameters for more than $1$ billion stars. In this paper, we present the expected sample size and properties (mass, magnitude, spatial, velocity distributions) of HVSs in the Gaia stellar catalogue. We build three Gaia mock catalogues of HVSs anchored to current observations, exploring different ejection mechanisms and GC stellar population properties. In all cases, we predict hundreds to thousands of HVSs with precise proper motion measurements within a few tens of kpc from us. For stars with a relative error in total proper motion below $10 \%$, the mass range extends to ~$10 M_{\odot}$ but peaks at ~$1$ $M_\odot$. The majority of Gaia HVSs will therefore probe a different mass and distance range compared to the current non-Gaia sample. In addition, a subset of a few hundreds to a few thousands of HVSs with $M$ ~ $3$ $M_\odot$ will be bright enough to have a precise measurement of the three-dimensional velocity from Gaia alone. Finally, we show that Gaia will provide more precise proper motion measurements for the current sample of HVS candidates. This will help identifying their birthplace narrowing down their ejection location, and confirming or rejecting their nature as HVSs. Overall, our forecasts are extremely encouraging in terms of quantity and quality of HVS data that can be exploited to constrain both the Milky Way potential and the GC properties.


title: On the selection of models for the high-frequency quasiperiodic oscillations of black holes: The mass problem

abstract: High-frequency quasiperiodic oscillations of the X-ray flux comming from black-hole binaries have the potential to provide precise estimates of the mass and the spin of the central black hole, if an adequate model for them was available. There are several models in the literature but none of them is commonly accepted. One way to test the available models is to confront their predictions for the masses of black holes with results obtained through other methods. Here, we study the mass bounds that nine of the most commonly used models provide for the three microquasars with known masses GRS 1915+105, GRO 1655-40 and XTE 1550-564 which display high-frequency quasiperiodic oscillations in their X-ray spectra. We also propose a statistical method for the assessment of the average success of the models. The results allow us to discard five of them. Here, "the mass problem" designates the conflict between their predictions and the reference masses.


title: A note on power generalized extreme value distribution and its properties

abstract: Similar to the generalized extreme value (GEV) family, the generalized extreme value distributions under power normalization are introduced by Roudsari (1999) and Barakat et al. (2013). In this article, we study the asymptotic behavior of GEV laws under power normalization and derive expressions for the kth moments, entropy, ordering in dispersion, rare event estimation and application of real data set. We also show that, under some conditions, the Shannon entropy and variance of GEV families are ordered.


title: Locally optimal control of continuous variable entanglement

abstract: We consider a system of two bosonic modes each subject to the dynamics induced by a thermal Markovian environment and we identify instantaneous, local symplectic controls that minimise the loss of entanglement in the Gaussian regime. By minimising the decrease of the logarithmic negativity at every instant in time, it will be shown that a non-trivial, finite amount of local squeezing helps to counter the effect of decoherence during the evolution. We also determine optimal control routines in the more restrictive scenario where the control operations are applied on only one of the two modes. We find that applying an instantaneous control only at the beginning of the dynamics, i.e. preparing an appropriate initial state, is the optimal strategy for states with symmetric correlations and when the dynamics is the same on both modes. More generally, even in asymmetric cases, the delayed decay of entanglement resulting from the optimal preparation of the initial state with no further action turns out to be always very close to the optimised control where multiple operations are applied during the evolution. Our study extends directly to mono-symmetric systems of any number of modes, i.e. to systems that are invariant under any local permutation of the modes within any one partition, as they are locally equivalent to two-mode systems.


title: Tensor-network study of quantum phase transition on Sierpi\'nski fractal

abstract: The transverse-field Ising model on the Sierpi\'nski fractal, which is characterized by the fractal dimension $\log_2^{~} 3 \approx 1.585$, is studied by a tensor-network method, the Higher-Order Tensor Renormalization Group. We analyze the ground-state energy and the spontaneous magnetization in the thermodynamic limit. The system exhibits the second-order phase transition at the critical transverse field $h_{\rm c}^{~} = 1.865$. The critical exponents $\beta \approx 0.198$ and $\delta \approx 8.7$ are obtained. Complementary to the tensor-network method, we make use of the real-space renormalization group and improved mean-field approximations for comparison.


title: Flat bands in fractal-like geometry

abstract: We report the presence of multiple flat bands in a class of two-dimensional (2D) lattices formed by Sierpinski gasket (SPG) fractal geometries as the basic unit cells. Solving the tight-binding Hamiltonian for such lattices with different generations of a SPG network, we find multiple degenerate and non-degenerate completely flat bands, depending on the configuration of parameters of the Hamiltonian. Moreover, we find a generic formula to determine the number of such bands as a function of the generation index $\ell$ of the fractal geometry. We show that the flat bands and their neighboring dispersive bands have remarkable features, the most interesting one being the spin-1 conical-type spectrum at the band center without any staggered magnetic flux, in contrast to the Kagome lattice. We furthermore investigate the effect of the magnetic flux in these lattice settings and show that different combinations of fluxes through such fractal unit cells lead to richer spectrum with a single isolated flat band or gapless electron- or hole-like flat bands. Finally, we discuss a possible experimental setup to engineer such fractal flat band network using single-mode laser-induced photonic waveguides.


title: KIBS Innovative Entrepreneurship Networks on Social Media

abstract: The analysis of the use of social media for innovative entrepreneurship in the context has received little attention in the literature, especially in the context of Knowledge Intensive Business Services (KIBS). Therefore, this paper focuses on bridging this gap by applying text mining and sentiment analysis techniques to identify the innovative entrepreneurship reflected by these companies in their social media. Finally, we present and analyze the results of our quantitative analysis of 23.483 posts based on eleven Spanish and Italian consultancy KIBS Twitter Usernames and Keywords using data interpretation techniques such as clustering and topic modeling. This paper suggests that there is a significant gap between the perceived potential of social media and the entrepreneurial behaviors at the social context in business-to-business (B2B) companies.


title: On negatively curved bundles with hyperbolic fibers outside the Igusa stable range

abstract: We prove that the Teichm\"{u}ller space $\mathcal{T}^{<0}(M)$ of negatively curved metrics on a hyperbolic manifold $M$ has nontrivial $i$-th rational homotopy groups for some $i> \dim M$. Moreover, some elements of infinite order in $\pi_i B\mbox{Diff}(M)$ can be represented by bundles over $S^i$ with fiberwise negatively curved metrics.


title: Low-Complexity Statistically Robust Precoder/Detector Computation for Massive MIMO Systems

abstract: Massive MIMO is a variant of multiuser MIMO in which the number of antennas at the base station (BS) $M$ is very large and typically much larger than the number of served users (data streams) $K$. Recent research has illustrated the system-level advantages of such a system and in particular the beneficial effect of increasing the number of antennas $M$. These benefits, however, come at the cost of dramatic increase in hardware and computational complexity. This is partly due to the fact that the BS needs to compute suitable beamforming vectors in order to coherently transmit/receive data to/from each user, where the resulting complexity grows proportionally to the number of antennas $M$ and the number of served users $K$. Recently, different algorithms based on tools from random matrix theory in the asymptotic regime of $M,K \to \infty$ with $\frac{K}{M} \to \rho \in (0,1)$ have been proposed to reduce such complexity. The underlying assumption in all these techniques, however, is that the exact statistics (covariance matrix) of the channel vectors of the users is a priori known. This is far from being realistic, especially that in the high-dim regime of $M\to \infty$, estimation of the underlying covariance matrices is well known to be a very challenging problem. In this paper, we propose a novel technique for designing beamforming vectors in a massive MIMO system. Our method is based on the randomized Kaczmarz algorithm and does not require knowledge of the statistics of the users channel vectors. We analyze the performance of our proposed algorithm theoretically and compare its performance with that of other competitive techniques via numerical simulations. Our results indicate that our proposed technique has a comparable performance while it does not require the knowledge of the statistics of the users channel vectors.


title: A new look at sulphur chemistry in hot cores and corinos

abstract: Sulphur-bearing species are often used to probe the evolution of hot cores since their abundances are particularly sensitive to physical and chemical variations. However, the chemistry of sulphur is not well understood in these regions, notably because observations of several hot cores displayed a large variety of sulphur compositions, and because the reservoir of sulphur in dense clouds, in which hot cores form, is still poorly constrained. In order to help disentangled its complexity, we present a fresh comprehensive review of sulphur chemistry in hot cores along with a study of its sensibility to temperature and pre-collapse chemical composition. In parallel, we analyse the discrepencies that result from the use of two different types of models (static and dynamic) to highlight the sensitivity to the choice of model to be used in astrochemical studies. Our results show that the pre-collapse chemical composition is a critical parameter for sulphur chemistry in hot cores and could explain the different sulphur compositions observed. We also report that differences in abundances for a given species between the static and dynamic models can reach six orders of magnitude in the hot core, which reveals the key role of the choice of model in astrochemical studies.


title: FPS-SFT: A Multi-dimensional Sparse Fourier Transform Based on the Fourier Projection-slice Theorem

abstract: We propose a multi-dimensional (M-D) sparse Fourier transform inspired by the idea of the Fourier projection-slice theorem, called FPS-SFT. FPS-SFT extracts samples along lines (1-dimensional slices from an M-D data cube), which are parameterized by random slopes and offsets. The discrete Fourier transform (DFT) along those lines represents projections of M-D DFT of the M-D data onto those lines. The M-D sinusoids that are contained in the signal can be reconstructed from the DFT along lines with a low sample and computational complexity provided that the signal is sparse in the frequency domain and the lines are appropriately designed. The performance of FPS-SFT is demonstrated both theoretically and numerically. A sparse image reconstruction application is illustrated, which shows the capability of the FPS-SFT in solving practical problems.


title: The identity of information: how deterministic dependencies constrain information synergy and redundancy

abstract: Understanding how different information sources together transmit information is crucial in many domains. For example, understanding the neural code requires characterizing how different neurons contribute unique, redundant, or synergistic pieces of information about sensory or behavioral variables. Williams and Beer (2010) proposed a partial information decomposition (PID) which separates the mutual information that a set of sources contains about a set of targets into nonnegative terms interpretable as these pieces. Quantifying redundancy requires assigning an identity to different information pieces, to assess when information is common across sources. Harder et al. (2013) proposed an identity axiom stating that there cannot be redundancy between two independent sources about a copy of themselves. However, Bertschinger et al. (2012) showed that with a deterministically related sources-target copy this axiom is incompatible with ensuring PID nonnegativity. Here we study systematically the effect of deterministic target-sources dependencies. We introduce two synergy stochasticity axioms that generalize the identity axiom, and we derive general expressions separating stochastic and deterministic PID components. Our analysis identifies how negative terms can originate from deterministic dependencies and shows how different assumptions on information identity, implicit in the stochasticity and identity axioms, determine the PID structure. The implications for studying neural coding are discussed.


title: ArCLight - a Compact Dielectric Large-Area Photon Detector

abstract: ArCLight is a novel device for detecting scintillation light over large areas with Photon Detection Efficiency (PDE) of the order of a few percent. Its robust technological design allows for efficient use in large-volume particle detectors, such as Liquid Argon Time Projection Chambers (LArTPCs) or liquid scintillator detectors. Due to its dielectric structure it can be placed inside volumes with high electric field. It could potentially replace vacuum PhotoMultiplier Tubes (PMTs) in applications where high PDE is not required. The photon detection efficiency for a 10x10cm2 detector prototype was measured to be in the range of 0.8% to 2.2% across the active area.


title: Sharp cohomological bound for uniformly quasiregularly elliptic manifolds

abstract: We show that if a compact, connected, and oriented $n$-manifold $M$ without boundary admits a non-constant non-injective uniformly quasiregular self-map, then the dimension of the real singular cohomology ring $H^*(M; \mathbb{R})$ of $M$ is bounded from above by $2^n$. This is a positive answer to a dynamical counterpart of the Bonk-Heinonen conjecture on the cohomology bound for quasiregularly elliptic manifolds. The proof is based on an intermediary result that, if $M$ is not a rational homology sphere, then each such uniformly quasiregular self-map on $M$ has a Julia set of positive Lebesgue measure.


title: An Empirical Study of Teaching Methodologies and Learning Outcomes for Online and in-class Networking Course Sections

abstract: To enhance student learning, we demonstrate an experimental study to analyze student learning outcomes in online and in-class sections of a core data communications course of the Undergraduate IT program in the Information Sciences and Technology (IST) Department at George Mason University (GMU). In this study, student performance is evaluated based on course assessments. This includes home and lab assignments, skill-based assessment, and traditional midterm exam across all 4 sections of the course. All sections have analogous content, assessment plan and teaching methodologies. Student demographics such as exam type and location preferences that may play an important role in their learning process are considered in our study. We had to collect vast amount of data from the learning management system (LMS), Blackboard (BB) Learn, in order to compare and study the results of several assessment outcomes for all students within their respective section and amongst students of other sections. We then tried to understand whether demographics have any influence on student performance by correlating individual student survey response to his/her performance in the class. The numerical results up to mid-semester reveal remarkable insights on student success in the online and face-to-face (F2F) sections.


title: Approximating Connected Safe Sets in Weighted Trees

abstract: For a graph $G$ and a non-negative integral weight function $w$ on the vertex set of $G$, a set $S$ of vertices of $G$ is $w$-safe if $w(C)\geq w(D)$ for every component $C$ of the subgraph of $G$ induced by $S$ and every component $D$ of the subgraph of $G$ induced by the complement of $S$ such that some vertex in $C$ is adjacent to some vertex of $D$. The minimum weight $w(S)$ of a $w$-safe set $S$ is the safe number $s(G,w)$ of the weighted graph $(G,w)$, and the minimum weight of a $w$-safe set that induces a connected subgraph of $G$ is its connected safe number $cs(G,w)$. Bapat et al. showed that computing $cs(G,w)$ is NP-hard even when $G$ is a star. For a given weighted tree $(T,w)$, they described an efficient $2$-approximation algorithm for $cs(T,w)$ as well as an efficient $4$-approximation algorithm for $s(T,w)$. Addressing a problem they posed, we present a PTAS for the connected safe number of a weighted tree. Our PTAS partly relies on an exact pseudopolynomial time algorithm, which also allows to derive an asymptotic FPTAS for restricted instances. Finally, we extend a bound due to Fujita et al. from trees to block graphs.


title: Mean-Square Performance Analysis of Noise-Robust Normalized Subband Adaptive Filter Algorithm

abstract: This paper studies the statistical models of the noise-robust normalized subband adaptive filter (NR-NSAF) algorithm in the mean and mean square deviation senses involving transient-state and steady-state behavior by resorting to the method of the vectorization operation and the Kronecker product. The analysis method does not require the Gaussian input signal. Moreover, the proposed analysis removes the paraunitary assumption imposed on the analysis filter banks as in the existing analyses of subband adaptive algorithms. Simulation results in various conditions demonstrate the effectiveness of our theoretical analysis. For a special form of the algorithm, the proposed steady-state expression is also better accurate than the previous analysis.


title: On density of subgraphs of halved cubes

abstract: Let $\mathcal S$ be a family of subsets of a set $X$ of cardinality $m$ and $\text{VC-dim}(\mathcal S)$ be the Vapnik-Chervonenkis dimension of $\mathcal S$. Haussler, Littlestone, and Warmuth (Inf. Comput., 1994) proved that if $G_1(\mathcal S)=(V,E)$ is the subgraph of the hypercube $Q_m$ induced by $\mathcal S$ (called the 1-inclusion graph of $\mathcal S$), then $\frac{|E|}{|V|}\le \text{VC-dim}({\mathcal S})$. Haussler (J. Combin. Th. A, 1995) presented an elegant proof of this inequality using the shifting operation. In this note, we adapt the shifting technique to prove that if $\mathcal S$ is an arbitrary set family and $G_{1,2}(\mathcal S)=(V,E)$ is the 1,2-inclusion graph of $\mathcal S$ (i.e., the subgraph of the square $Q^2_m$ of the hypercube $Q_m$ induced by $\mathcal S$), then $\frac{|E|}{|V|}\le \binom{d}{2}$, where $d:=\text{cVC-dim}^*(\mathcal S)$ is the clique-VC-dimension of $\mathcal S$ (which we introduce in this paper). The 1,2-inclusion graphs are exactly the subgraphs of halved cubes and comprise subgraphs of Johnson graphs as a subclass.


title: Real elliptic curves and cevian geometry

abstract: We study the elliptic curve $E_a: (ax+1)y^2+(ax+1)(x-1)y+x^2-x=0$, which we call the geometric normal form of an elliptic curve. We show that any elliptic curve whose $j$-invariant is real is isomorphic to a curve $E_a$ in geometric normal form, and show that for $a \notin \{0, -1, -9\}$, the points on $E_a$, minus a set of $6$ points, can be characterized in terms of the cevian geometry of a triangle.


title: Riemann hypothesis

abstract: This work is dedicated to the promotion of the results Hadamard, Landau E., Walvis A., Estarmann T and Paul R. Chernoff for pseudo zeta functions. The properties of zeta functions are studied, these properties can lead to new regularities of zeta functions.


title: Scalable synthesis of safety certificates from data with application to learning-based control

abstract: The control of complex systems faces a trade-off between high performance and safety guarantees, which in particular restricts the application of learning-based methods to safety-critical systems. A recently proposed framework to address this issue is the use of a safety controller, which guarantees to keep the system within a safe region of the state space. This paper introduces efficient techniques for the synthesis of a safe set and control law, which offer improved scalability properties by relying on approximations based on convex optimization problems. The first proposed method requires only an approximate linear system model and Lipschitz continuity of the unknown nonlinear dynamics. The second method extends the results by showing how a Gaussian process prior on the unknown system dynamics can be used in order to reduce conservatism of the resulting safe set. We demonstrate the results with numerical examples, including an autonomous convoy of vehicles.


title: mVMC - Open-source software for many-variable variational Monte Carlo method

abstract: mVMC (many-variable Variational Monte Carlo) is an open-source software based on the variational Monte Carlo method applicable for a wide range of Hamiltonians for interacting fermion systems. In mVMC, we introduce more than ten thousands variational parameters and simultaneously optimize them by using the stochastic reconfiguration (SR) method. In this paper, we explain basics and user interfaces of mVMC. By using mVMC, users can perform the calculation by preparing only one input file of about ten lines for widely studied quantum lattice models, and can also perform it for general Hamiltonians by preparing several additional input files. We show the benchmark results of mVMC for the Hubbard model, the Heisenberg model, and the Kondo-lattice model. These benchmark results demonstrate that mVMC provides ground-state and low-energy-excited-state wave functions for interacting fermion systems with high accuracy.


title: Leader-Based Optimal Coordination Control for the Consensus Problem of Multiagent Differential Games via Fuzzy Adaptive Dynamic Programming

abstract: In this paper, a new on-line scheme is presented to design the optimal coordination control for the consensus problem of multi-agent differential games by fuzzy adaptive dynamic programming (FADP), which brings together game theory, generalized fuzzy hyperbolic model (GFHM) and adaptive dynamic programming. In general, the optimal coordination control for multi-agent differential games is the solution of the coupled Hamilton-Jacobi (HJ) equations. Here, for the first time, GFHMs are used to approximate the solution (value functions) of the coupled HJ equations, based on policy iteration (PI) algorithm. Namely, for each agent, GFHM is used to capture the mapping between the local consensus error and local value function. Since our scheme uses the single-network rchitecture for each agent (which eliminates the action network model compared with dual-network architecture), it is a more reasonable architecture for multi-agent systems. Furthermore, the approximation solution is utilized to obtain the optimal coordination controls. Finally, we give the stability analysis for our scheme, and prove the weight estimation error and the local consensus error are uniformly ultimately bounded. Further, the control node trajectory is proven to be cooperative uniformly ultimately bounded.


title: High-Order method for Darcy flows in fractured porous media

abstract: We develop a novel Hybrid High-Order method for the simulation of Darcy flows in fractured porous media. The discretization hinges on a mixed formulation in the bulk region and on a primal formulation inside the fracture. Salient features of the method include a seamless treatment of nonconforming discretizations of the fracture, as well as the support of arbitrary approximation orders on fairly general meshes. For the version of the method corresponding to a polynomial degree k \v{e} 0, we prove convergence in h^{k+1} of the discretization error measured in an energy-like norm. In the error estimate, we explicitly track the dependence of the constants on the problem data, showing that the method is fully robust with respect to the heterogeneity of the permeability coefficients, and it exhibits only a mild dependence on the square root of the local anisotropy of the bulk permeability. The numerical validation on a comprehensive set of test cases confirms the theoretical results.