title: Description of growth and oscillation of solutions of complex LDE's

abstract: It is known that, equally well in the unit disc as in the whole complex plane, the growth of the analytic coefficients $A_0,\dotsc,A_{k-2}$ of \begin{equation*} f^{(k)} + A_{k-2} f^{(k-2)} + \dotsb + A_1 f'+ A_0 f = 0, \quad k\geq 2, \end{equation*} determines, under certain growth restrictions, not only the growth but also the oscillation of its non-trivial solutions, and vice versa. A uniform treatment of this principle is given in the disc $D(0,R)$, $0<R\leq \infty$, by using several measures for growth that are more flexible than those in the existing literature, and therefore permit more detailed analysis. In particular, results obtained are not restricted to cases where solutions are of finite (iterated) order of growth in the classical sense. The new findings are based on an accurate integrated estimate for logarithmic derivatives of meromorphic functions, which preserves generality in terms of three free parameters.


title: Magnetized Gas Collimation of Interstellar Outflow Scaled by Laser-produced Plasma

abstract: Young stellar objects/planetary nebula outflow anisotropies usually involve wind-wind interactions and magnetic collimation, but detailed structures of wind and magnetic fields inside collimation region remain undetermined. We numerically investigated its laboratory counterpart, based on poloidal field collimation in magnetocentrifugal launching model. Our analog consist of fast wind: Aluminum plasma generated by laser and magnetized ambient: molecular Helium and B=5-60 Tesla embedded field. Elevating magnetic field strength or decreasing gas density can alter expansion morphology, from sphere to prolonged cavity and ultimately to collimated jet. Outflow patterns can be quantitatively predicted based on the knowledge of its surroundings through a set of external Mach numbers. We conclude that such mixed gas and magnetic field dynamics are consistent with astronomical observations of protostars and planetary nebulae in certain evolution stages, and here provide a scalable framework allowing fitting of flow-field structures in astronomical unresolved regions by assuming their possible geometries on a repeatable laboratory platform.


title: Dimensional Regularization and Dispersive Two-Loop Calculations

abstract: The two-loop contributions are now often required by the precision experiments, yet are hard to express analytically while keeping precision. One way to approach this challenging task is via the dispersive approach, allowing to replace sub-loop diagram by effective propagator. This paper builds on our previous work, where we developed a general approach based on representation of many-point Passarino-Veltman functions in two-point function basis. In this work, we have extracted the UV-divergent poles of the Passarino-Veltman functions analytically and presented them as the dimensionally-regularized and multiply-subtracted dispersive sub-loop insertions, including self-energy, triangle, box and pentagon type.


title: Reservoir Computing on Spin-Torque Oscillator Array

abstract: We numerically study reservoir computing on a spin-torque oscillator (STO) array, describing the magnetization dynamics of the STO array by a nonlinear oscillator model. The STOs exhibit synchronized oscillation due to coupling by magnetic dipolar fields. We show that reservoir computing can be performed using the synchronized oscillation state. The performance can be improved by increasing the number of STOs. The performance becomes highest at the boundary between the synchronized and disordered states. Using an STO array, we can achieve higher performance than that of an echo-state network with similar number of units. This result indicates that STO arrays are promising for hardware implementation of reservoir computing.


title: On the density or measure of sets and their sumsets in the integers or the circle

abstract: Let $\mathrm{d}(A)$ be the asymptotic density (if it exists) of a sequence of integers $A$. For any real numbers $0\leq\alpha\leq\beta\leq 1$, we solve the question of the existence of a sequence $A$ of positive integers such that $\mathrm{d}(A)=\alpha$ and $\mathrm{d}(A+A)=\beta$. More generally we study the set of $k$-tuples $(\mathrm{d}(iA))_{1\leq i\leq k}$ for $A\subset \mathbb{N}$. This leads us to introduce subsets defined by diophantine constraints inside a random set of integers known as the set of ``pseudo $s$th powers''. We consider similar problems for subsets of the circle $\mathbb{R}/\mathbb{Z}$, that is, we partially determine the set of $k$-tuples $(\mu(iA))_{1\leq i\leq k}$ for $A\subset \mathbb{R}/\mathbb{Z}$.


title: The Poisson bracket invariant for open covers consisting of topological disks on surfaces

abstract: L. Buhovsky, A. Logunov and S. Tanny proved the (strong) Poisson bracket conjecture by Leonid Polterovich in dimension $2$. In this note, instead of open cover consisting of displaceable sets in their work, we consider open cover constituted of topological discs and give a necessary and sufficient condition that Poisson bracket invariants of these covers are positive.


title: Privacy-Preserving P2P Energy Market on the Blockchain

abstract: Quartierstrom creates a peer-to-peer marketplace for locally generated solar power. The marketplace is implemented as a smart contract on a permissioned blockchain governed by all prosumers. Two privacy-by-design concepts are presented which guarantee that the users individual load profile is not leaked to any third party despite using a blockchain. The first approach leverages UTXO based coin mixing protocols in combination with an account-based on-chain smart contract. The second approach relies on an off-chain smart contract running in trusted execution environments.


title: As simple as possible but not simpler in Multiple Criteria Decision Aiding: the ROR-SMAA level dependent Choquet integral approach

abstract: The level dependent Choquet integral has been proposed to handle decision making problems in which the importance and the interaction of criteria may depend on the level of the alternatives' evaluations. This integral is based on a level dependent capacity, which is a family of single capacities associated to each level of evaluation for the considered criteria. We present two possible formulations of the level dependent capacity where importance and interaction of criteria are constant inside each one of the subintervals in which the interval of evaluations for considered criteria is split or vary with continuity inside the whole interval of evaluations. Since, in general, there is not only one but many level dependent capacities compatible with the preference information provided by the Decision Maker, we propose to take into account all of them by using the Robust Ordinal Regression (ROR) and the Stochastic Multicriteria Acceptability Analysis (SMAA). On one hand, ROR defines a necessary preference relation (if an alternative a is at least as good as an alternative b for all compatible level dependent capacities), and a possible preference relation (if a is at least as good as b for at least one compatible level dependent capacity). On the other hand, considering a random sampling of compatible level dependent capacities, SMAA gives the probability that each alternative reaches a certain ranking position as well as the probability that an alternative is preferred to another. A real-world decision problem on rankings of universities is provided to illustrate the proposed methodology.


title: Small perturbations for a Duffing-like evolution equation involving non-commuting operators

abstract: We consider an abstract evolution equation with linear damping, a nonlinear term of Duffing type, and a small forcing term. The abstract problem is inspired by some models for damped oscillations of a beam subject to external loads or magnetic fields, and shaken by a transversal force. The main feature is that very natural choices of the boundary conditions lead to equations whose linear part involves two operators that do not commute. We extend to this setting the results that are known in the commutative case, namely that for asymptotically small forcing terms all solutions are eventually close to the three equilibrium points of the unforced equation, two stable and one unstable.


title: Suppressing the thermal SZ-induced variance in CMB-cluster lensing estimators

abstract: Accurate galaxy cluster mass measurements from the gravitational lensing of the cosmic microwave background temperature maps depend on mitigating potential biases from the cluster's own thermal Sunyaev-Zel'dovich (SZ) effect signal. Quadratic lensing estimators use a pair of maps to extract the lensing signal: a large scale gradient map and a small scale lensing map. The SZ bias can be eliminated by using an SZ-free map in the pair, with the gradient map being favored for signal-to-noise reasons. However, while this approach eliminates the bias, the SZ power in small scale lensing map adds extra variance that can become significant for high mass clusters and low noise surveys. In this work, we propose projecting out an SZ template to reduce the SZ variance. Any residual SZ signal after template fitting is uncorrelated with the SZ-free gradient map, and thus does not bias the mass measurements. For massive clusters above $4\times 10^{14}$ $M_{\odot}$ observed by the upcoming CMB-S4 or Simons Observatory experiments, we find that the template fitting approach would increase the cluster lensing signal-to-noise by a factor of 1.4.


title: Traces of reciprocal singular moduli

abstract: We show that the generating series of traces of reciprocal singular moduli is a mixed mock modular form of weight $3/2$ whose shadow is given by a linear combination of products of unary and binary theta functions. To prove these results, we extend the Kudla-Millson theta lift of Bruinier and Funke to meromorphic modular functions.


title: Realization of $\mathcal{PT}$-symmetric and $\mathcal{PT}$-symmetry broken states in static optical lattice potentials

abstract: R. Labouvie et al. [Phys. Rev. Lett. 116, 235302 (2016)] have described an experiment with a weakly interacting Bose-Einstein condensate trapped in a one-dimensional optical lattice with localized loss created by a focused electron beam. We show that by setting suitable initial currents between neighboring sites it is possible to create $\mathcal{PT}$-symmetric quasi-stationary and $\mathcal{PT}$-symmetry broken decaying states in an embedded two-mode subsystem. This subsystem exhibits gain provided by the coupling to the reservoir sites and localized loss due to the electron beam, and shows the same dynamics as a non-Hermitian two-mode system with symmetric real and antisymmetric imaginary time-independent potentials, except for a proportionality factor in the chemical potential. We also show that there are three other equivalent scenarios, and that the presence of a localized loss term significantly reduces the size of the condensate required for the realization.


title: TTbar deformation and the light-cone gauge

abstract: The homogeneous inviscid Burgers equation which determines the spectrum of a TTbar deformed model has a natural interpretation as the condition of the gauge invariance of the target space-time energy and momentum of a (non-critical) string theory quantised in a generalised uniform light-cone gauge which depends on the deformation parameter. As a simple application of the light-cone gauge interpretation we derive the TTbar deformed Lagrangian for a system of any number of scalars, fermions and chiral bosons with an arbitrary potential. We find that the TTbar deformation is driven by the canonical Noether stress-energy tensor but not the covariant one.


title: Spin dynamics of hot excitons in diluted magnetic semiconductors with spin-orbit interaction

abstract: We explore the impact of a Rashba-type spin-orbit interaction in the conduction band on the spin dynamics of hot excitons in diluted magnetic semiconductor quantum wells. In materials with strong spin-orbit coupling, we identify parameter regimes where spin-orbit effects greatly accelerate the spin decay and even change the dynamics qualitatively in the form of damped oscillations. Furthermore, we show that the application of a small external magnetic field can be used to either mitigate the influence of spin-orbit coupling or entirely remove its effects for fields above a material-dependent threshold.


title: Cascaded Channel Estimation for Large Intelligent Metasurface Assisted Massive MIMO

abstract: In this letter, we consider the problem of channel estimation for large intelligent metasurface (LIM) assisted massive multiple-input multiple-output (MIMO) systems. The main challenge of this problem is that the LIM integrated with a large number of low-cost metamaterial antennas can only passively reflect the incident signal by a certain phase shift, and does not have any signal processing capability. To deal with this, we introduce a general framework for the estimation of the transmitter-LIM and LIM-receiver cascaded channel, and propose a two-stage algorithm that includes a sparse matrix factorization stage and a matrix completion stage. Simulation results illustrate that the proposed method can achieve accurate channel estimation for LIM-assisted massive MIMO systems.


title: Simulation of Pedestrian Evacuation Using Conjugated Forces Model to Explore the Optimal Bypassing Strategy

abstract: The conjugated forces model (CFM) capable of reproducing bypassing behaviour is proposed and adopted to simulate pedestrian evacuation. Primarily, the concept of collision and evading and surpassing behaviour is particularly defined to describe the proposed CFM. A pedestrian whose behaviour can be described by the CFM will change his desired direction when confronted with a collision, and choose to actively bypass others. Secondly, pedestrian movement in a passageway is simulated to testify the capacity of the CFM to reproduce microscopic bypassing behaviour, and the bypassing parameter is then introduced to describe different movement strategies. Finally, Evacuation form a one-exit room is simulated to find the optimal bypassing strategy. It is concluded that the optimal bypassing strategy varies with different exit width and initial pedestrian number.


title: A survey of molecular cores in M17 SWex

abstract: A survey of molecular cores covering the infrared dark cloud known as the M17 southwest extension (M17 SWex) has been carried out with the 45 m Nobeyama Radio Telescope. Based on the N2H+ (J= 1-0) data obtained, we have identified 46 individual cores whose masses are in the range 43 to 3026 Mo. We examined the relationship between the physical parameters of the cores and those of young stellar objects (YSOs) associated with the cores found in the literature. The comparison of the virial mass and the core mass indicates that most of the cores can be gravitationally stable if we assume a large external pressure. Among the 46 cores, we found four massive cores with YSOs. They have large mass of >~ 1000Mo and line width of >~ 2.5 km s^-1 which are similar to those of clumps forming high mass stars. However, previous studies have shown that there is no active massive star formation in this region. Recent measurements of near-infrared polarization infer that the magnetic field around M17 SWex is likely to be strong enough to support the cores against self-gravity. We therefore suggest that the magnetic field may prevent the cores from collapsing, causing the low-level of massive star formation in M17 SWex.


title: Band alignment and scattering considerations for enhancing the thermoelectric power factor of complex materials: The case of Co-based half-Heuslers

abstract: Half-Heuslers, an emerging thermoelectric material group, has complex bandstructures with multiple bands that can be aligned through band engineering approaches, giving us an opportunity to improve their power factor. In this work, going beyond the constant relaxation time approximation, we perform an investigation of the benefits of band alignment in improving the thermoelectric power factor under different density of states dependent scattering scenarios. As a test case we consider the Co-based p-type half-Heuslers TiCoSb, NbCoSn and ZrCoSb. First, using simplified effective mass models combined with Boltzmann transport, we investigate the conditions of band alignment that are beneficial to the thermoelectric power factor under three different carrier scattering scenarios: i) the usual constant relaxation time approximation, ii) intra-band scattering restricted to the current valley with the scattering rates proportional to the density of states as dictated by Fermi's Golden Rule, and iii) both intra- and inter-band scattering across all available valleys, with the rates determined by the total density of states at the relevant energies. We demonstrate that the band-alignment outcome differs significantly depending on the scattering details. Next, using the density functional theory calculated bandstructures of the half-Heuslers we study their power factor behavior under strain induced band alignment. We show that strain can improve the power factor of half-Heuslers, but the outcome heavily depends on the curvatures of the bands involved, the specifics of the carrier scattering mechanisms, and the initial band separation. Importantly, we also demonstrate that band alignment is not always beneficial to the power factor.


title: A Riesz basis criterion for Schr\"odinger operators with boundary conditions dependent on the eigenvalue parameter

abstract: We establish a criterion for a set of eigenfunctions of the one-dimensional Schr\"{o}dinger operator with distributional potentials and boundary conditions containing the eigenvalue parameter to be a Riesz basis for $\mathscr{L}_2(0,\pi)$.


title: Cluster-GCN: An Efficient Algorithm for Training Deep and Large Graph Convolutional Networks

abstract: Graph convolutional network (GCN) has been successfully applied to many graph-based applications; however, training a large-scale GCN remains challenging. Current SGD-based algorithms suffer from either a high computational cost that exponentially grows with number of GCN layers, or a large space requirement for keeping the entire graph and the embedding of each node in memory. In this paper, we propose Cluster-GCN, a novel GCN algorithm that is suitable for SGD-based training by exploiting the graph clustering structure. Cluster-GCN works as the following: at each step, it samples a block of nodes that associate with a dense subgraph identified by a graph clustering algorithm, and restricts the neighborhood search within this subgraph. This simple but effective strategy leads to significantly improved memory and computational efficiency while being able to achieve comparable test accuracy with previous algorithms. To test the scalability of our algorithm, we create a new Amazon2M data with 2 million nodes and 61 million edges which is more than 5 times larger than the previous largest publicly available dataset (Reddit). For training a 3-layer GCN on this data, Cluster-GCN is faster than the previous state-of-the-art VR-GCN (1523 seconds vs 1961 seconds) and using much less memory (2.2GB vs 11.2GB). Furthermore, for training 4 layer GCN on this data, our algorithm can finish in around 36 minutes while all the existing GCN training algorithms fail to train due to the out-of-memory issue. Furthermore, Cluster-GCN allows us to train much deeper GCN without much time and memory overhead, which leads to improved prediction accuracy---using a 5-layer Cluster-GCN, we achieve state-of-the-art test F1 score 99.36 on the PPI dataset, while the previous best result was 98.71 by [16]. Our codes are publicly available at https://github.com/google-research/google-research/tree/master/cluster_gcn.


title: A novel Algorithm for Optimal Placement of Multiple Inertial Sensors to Improve the Sensing Accuracy

abstract: This paper proposes a novel algorithm to determine the optimal placement of redundant inertial sensors such as accelerometers and gyroscopes (gyros) for increasing the sensing accuracy. In this paper, we have proposed a novel iterative algorithm to find the optimal sensor configuration. The proposed algorithm utilizes the majorization-minimization (MM) algorithm and the duality principle to find the optimal configuration. Unlike the state-of-the-art which are mainly geometrical in nature and restricted to certain noise statistics, the proposed algorithm gives the exact positions of the sensors, and moreover, the proposed algorithm is independent of the nature of the noise at different sensors. The proposed alogrithm has been implemented and tested via numerical simulation in the MATLAB. The simulation results show that the algorithm converges to the optimal configurations and show the effectiveness of the proposed algorithm.


title: Squeezed vacuum states from a whispering gallery mode resonator

abstract: Squeezed vacuum states enable optical measurements below the quantum limit and hence are a valuable resource for applications in quantum metrology and also quantum communication. However, most available sources require high pump powers in the milliwatt range and large setups, which hinders real world applications. Furthermore, degenerate operation of such systems presents a challenge. Here, we use a compact crystalline whispering gallery mode resonator made of lithium niobate as a degenerate parametric oscillator. We demonstrate about 1.4 dB noise reduction below the shot noise level for only 300 $\mu\text{W}$ of pump power in degenerate single mode operation. Furthermore, we report a record pump threshold as low as 1.35 $\mu\text{W}$. Our results show that the whispering gallery based approach presents a promising platform for a compact and efficient source for nonclassical light.


title: Ready-to-use Fourier domain templates for compact binaries inspiraling along moderately eccentric orbits

abstract: We derive analytic expressions that provide Fourier domain gravitational wave (GW) response function for compact binaries inspiraling along moderately eccentric orbits. These expressions include amplitude corrections to the two GW polarization states that are accurate to the first post-Newtonian (PN) order. Additionally, our fully 3PN accurate GW phase evolution incorporates eccentricity effects up to sixth order at each PN order. Further, we develop a prescription to incorporate analytically the effects of 3PN accurate periastron advance in the GW phase evolution. This is how we provide a ready-to-use and efficient inspiral template family for compact binaries in moderately eccentric orbits. Preliminary GW data analysis explorations suggest that our template family should be required to construct analytic inspiral-merger-ringdown templates to model moderately eccentric compact binary coalescence.


title: Simultaneous Conjugacy Classes as Combinatorial Invariants of Finite Groups

abstract: Let $G$ be a finite group. We consider the problem of counting simultaneous conjugacy classes of $n$-tuples and simultaneous conjugacy classes of commuting $n$-tuples in $G$. Let $\alpha_{G,n}$ denote the number of simultaneous conjugacy classes of $n$-tuples, and $\beta_{G,n}$ the number of simultaneous conjugacy classes of commuting $n$-tuples in $G$. The generating functions $A_G(t) = \sum_{n\geq 0} \alpha_{G,n}t^n,$ and $B_G(t) = \sum_{n\geq 0} \beta_{G,n}t^n$ are rational functions of $t$. We show that $A_G(t)$ determines and is completely determined by the class equation of $G$. We show that $\alpha_{G,n}$ grows exponentially with growth factor equal to the cardinality of $G$, whereas $\beta_{G,n}$ grows exponentially with growth factor equal to the maximum cardinality of an abelian subgroup of $G$. The functions $A_G(t)$ and $B_G(t)$ may be regarded as combinatorial invariants of the finite group $G$. We study dependencies amongst these invariants and the notion of isoclinism for finite groups. We prove that the normalized functions $A_G(t/|G|)$ and $B_G(t/|G|)$ are invariants of isoclinism families.


title: Synthesis of complex organic molecules in soft x-ray irradiated ices

abstract: We study the chemical evolution of H2O:CO:NH3 ice mixtures irradiated with soft X-rays, in the range 250-1250 eV. We identify many nitrogen-bearing molecules such as e.g., OCN-, NH4+ , HNCO, CH3CN, HCONH2, and NH2COCONH2. Several infrared features are compatible with glycine or its isomers. During the irradiation, we detected through mass spectroscopy many species desorbing the ice. Such findings support either the infrared identifications and reveal less abundant species with not clear infrared features. Among them, m/z = 57 has been ascribed to methyl isocyanate (CH3NCO), a molecule of prebiotic relevance, recently detected in protostellar environments. During the warm up after the irradiation, several infrared features including 2168 cm-1 band of OCN-, 1690 cm-1 band of formamide, and the 1590 cm-1 band associated to three different species, HCOO-, CH3NH2 and NH3+CH2COO survive up to room temperature. Interestingly, many high masses have been also detected. Possible candidates are methyl-formate, (m/z = 60, HCOOCH3), ethanediamide (m/z = 88, NH2COCONH2), and N-acetyl-L-aspartic acid (m=z = 175). This latter species is compatible with the presence of the m/z = 43, 70 and 80 fragments. Photo-desorption of organics is relevant for the detection of such species in the gas-phase of cold environments, where organic synthesis in ice mantles should dominate. We estimate the gas-phase enrichment of some selected species in the light of a protoplanetary disc model around young solar-type stars.


title: System size and Multiplicity dependence of Chemical freeze-out parameters at the Large Hadron Collider Energies

abstract: The collision system and multiplicity dependence of chemical freeze-out temperature ($T_{\rm ch}$) and strangeness saturation factor ($\gamma_{s}$) are obtained by studying the particle ratios at the Large Hadron Collider (LHC) energies. Here, we consider the new results in pp at 13 TeV, p+Pb at $\sqrt{s_{\rm NN}}$ = 5.02 TeV, Xe+Xe at $\sqrt{s_{\rm NN}}$ = 5.44 TeV and Pb+Pb at $\sqrt{s_{\rm NN}}$ = 5.02 TeV along with the earlier results in pp at $\sqrt{s}$ = 7 TeV and Pb+Pb at $\sqrt{s_{\rm NN}}$ = 2.76 TeV. A statistical thermal model is used to extract the chemical freeze-out parameters in different multiplicity classes. To understand the particle production from small to large collision systems two ensembles namely, canonical and grand canonical have been considered in this study. A clear observation of multiplicity dependence of $T_{\rm ch}$ and $\gamma_{s}$ is observed. The values obtained in high-multiplicity pp collisions are found to be similar to the peripheral Pb+Pb collisions. A final state midrapidity charged particle multiplicity density of around 20-30 appears to be a threshold below which, the chemical freeze-out temperature is lower than the kinetic freeze-out temperature.


title: A novel Multiplicative Polynomial Kernel for Volterra series identification

abstract: Volterra series are especially useful for nonlinear system identification, also thanks to their capability to approximate a broad range of input-output maps. However, their identification from a finite set of data is hard, due to the curse of dimensionality. Recent approaches have shown how regularized kernel-based methods can be useful for this task. In this paper, we propose a new regularization network for Volterra models identification. It relies on a new kernel given by the product of basic building blocks. Each block contains some unknown parameters that can be estimated from data using marginal likelihood optimization. In comparison with other algorithms proposed in the literature, numerical experiments show that our approach allows to better select the monomials that really influence the system output, much increasing the prediction capability of the model.


title: Guiding Inferences in Connection Tableau by Recurrent Neural Networks

abstract: We present a dataset and experiments on applying recurrent neural networks (RNNs) for guiding clause selection in the connection tableau proof calculus. The RNN encodes a sequence of literals from the current branch of the partial proof tree to a hidden vector state; using it, the system selects a clause for extending the proof tree. The training data and learning setup are described, and the results are discussed and compared with state of the art using gradient boosted trees. Additionally, we perform a conjecturing experiment in which the RNN does not just select an existing clause, but completely constructs the next tableau goal.


title: Topology analysis for anomalous Hall effect in the non-collinear antiferromagnetic states of Mn$_3\mathit{A}$N ($A$ = Ni, Cu, Zn, Ga, Ge, Pd, In, Sn, Ir, Pt)

abstract: We investigate topological features of electronic structures which produce large anomalous Hall effect in the non-collinear antiferromagnetic metallic states of anti-perovskite manganese nitrides by first-principles calculations. We first predict the stable magnetic structures of these compounds to be non-collinear antiferromagnetic structures characterized by either $T_{1g}$ or $T_{2g}$ irreducible representation by evaluating the total energy for all of the magnetic structures classified according to the symmetry and multipole moments. The topology analysis is next performed for the Wannier tight-binding models obtained from the first-principles band structures. Our results reveal the small Berry curvature induced through the coupling between occupied and unoccupied states with the spin-orbit coupling, which is widely spread around the Fermi surface in the Brillouin zone, dominantly contributes after the $k$-space integration to the anomalous Hall conductivity, while the local divergent Berry curvature around Weyl points has a rather small contribution to the anomalous Hall conductivity.


title: An effective scanning method of the NMSSM parameter space

abstract: The next-to-minimal supersymmetric standard model (NMSSM) naturally provides a 125 GeV Higgs boson without the need for large loop corrections from multi-TeV stop quarks. Furthermore, the NMSSM provides an electroweak scale dark matter candidate consistent with all experimental data, like relic density and non-observation of direct dark matter signals with the present experimental sensitivity. However, more free parameters are introduced in the NMSSM, which are strongly correlated. A simple parameter scan without knowing the correlation matrix is not efficient and can miss significant regions of the parameter space. We introduce a new technique to sample the NMSSM parameter space, which takes into account the correlations. For this we project the 7D NMSSM parameter space onto the 3D Higgs boson mass parameter space. The reduced dimensionality allows for a non-random sampling and therefore a complete coverage of the allowed NMSSM parameters. In addition, the parameter correlations and possible deviations of the signal strengths of the observed 125 Higgs boson from the SM values are easily predicted.


title: Diagnosing the performance of human mobility models at small spatial scales using volunteered geographic information

abstract: Accurate modelling of local population movement patterns is a core contemporary concern for urban policymakers, affecting both the short term deployment of public transport resources and the longer term planning of transport infrastructure. Yet, while macro-level population movement models (such as the gravity and radiation models) are well developed, micro-level alternatives are in much shorter supply, with most macro-models known to perform badly in smaller geographic confines. In this paper we take a first step to remedying this deficit, by leveraging two novel datasets to analyse where and why macro-level models of human mobility break down at small scales. In particular, we use an anonymised aggregate dataset from a major mobility app and combine this with freely available data from OpenStreetMap concerning land-use composition of different areas around the county of Oxfordshire in the United Kingdom. We show where different models fail, and make the case for a new modelling strategy which moves beyond rough heuristics such as distance and population size towards a detailed, granular understanding of the opportunities presented in different areas of the city.


title: Multivariate Alexander quandles, III. Sublinks

abstract: If $L$ is a classical link then the multivariate Alexander quandle, $Q_A(L)$, is a substructure of the multivariate Alexander module, $M_A(L)$. In the first paper of this series we showed that if two links $L$ and $L'$ have $Q_A(L) \cong Q_A(L')$, then after an appropriate re-indexing of the components of $L$ and $L'$, there will be a module isomorphism $M_A(L) \cong M_A(L')$ of a particular type, which we call a"Crowell equivalence." In the present paper we show that $Q_A(L)$ (up to quandle isomorphism) is a strictly stronger link invariant than $M_A(L)$ (up to re-indexing and Crowell equivalence). This result follows from the fact that $Q_A(L)$ determines the $Q_A$ quandles of all the sublinks of $L$, up to quandle isomorphisms.


title: Utilizing the redundant constraints for the uplift payment elimination

abstract: A power market with non-convexities may not have an equilibrium price for power that provides economic stability of the centralized dispatch outcome. In this case, the market players are entitled to receive the uplift payments that compensate the economic profit lost when following the centralized dispatch. We consider a special class of the (possibly non-linear) redundant constraints that are redundant not only on the feasible set of the centralized dispatch optimization problem (and, therefore, do not change the centralized dispatch outcome) but also on the larger set obtained when the power balance constraint is relaxed. We show that the Lagrangian relaxation of these redundant constraints may reduce the uplift payments without changing the duality gap. For any given market price (or a pricing algorithm that sets the producer revenue as a function of its output volume) in a uninode multi-period power market with fixed load, we explicitly construct a family of the redundant constraints that do not change the maximum profit of the producer and result in zero uplift payment. We show that the introduction and subsequent Lagrangian relaxation of just one redundant constraint in the centralized dispatch problem suffice to eliminate the uplift payments for all the producers. In the case of the convex hull pricing method, the introduction of these redundant constraints affects neither the duality gap nor the market price for power. The results can be straightforwardly generalized to a power market with the price-sensitive load.


title: Spin Detection in Robotic Table Tennis

abstract: In table tennis, the rotation (spin) of the ball plays a crucial role. A table tennis match will feature a variety of strokes. Each generates different amounts and types of spin. To develop a robot that can compete with a human player, the robot needs to detect spin, so it can plan an appropriate return stroke. In this paper we compare three methods to estimate spin. The first two approaches use a high-speed camera that captures the ball in flight at a frame rate of 380 Hz. This camera allows the movement of the circular brand logo printed on the ball to be seen. The first approach uses background difference to determine the position of the logo. In a second alternative, we train a CNN to predict the orientation of the logo. The third method evaluates the trajectory of the ball and derives the rotation from the effect of the Magnus force. This method gives the highest accuracy and is used for a demonstration. Our robot successfully copes with different spin types in a real table tennis rally against a human opponent.


title: Light bending and gravitational lensing in Brans-Dicke theory

abstract: As an important candidate theory of gravity, Brans-Dicke theory has been widely studied. In this paper, we investigate light bending and gravitational lensing by compact objects in Brans-Dicke theory in weak gravitational field. Firstly, we present a general formalism for calculating higher-order corrections to light bending angle and lensing observables for a static, spherically symmetric and flat spacetime, in which the metric is given in the isotropic coordinates. Secondly, we apply the general formalism to Brans-Dicke theory and get the corresponding light bending angle and lensing observables. Our results show that, although the sums over the low-order correction terms in magnifications of the primary and secondary images do not dependent on the theories of gravity, the sums over correction terms with order higher than three do. Moreover, we show that the total magnification has a non-vanishing first-order correction, rather than a vanishing contribution concluded in the literature. We find that the corrections to lensing observables of BD theory close to those of GR when the parameter $\omega$ tends to $+\infty$ from $-\frac32$, while opposition occurs when $\omega$ tends to $-2$ from $-\infty$.


title: Photometric analysis and evolutionary stages of the contact binary V2790 Ori

abstract: A photometric analysis and evolutionary stages of the contact binary V2790~Ori is presented. The $BVR_\mathrm{C}$ observations were carried out at the Thai National Observatory. The photometric light curves were fitted to provide fundamental parameters, required to use in the analysis of evolutionary stages of the binary. The results show that V2790~Ori is a W-type contact system with a mass ratio of $q$~=~2.932. The orbital period increase is found at a rate of $\mathrm{d}P\slash \mathrm{d}t~=~$1.03$\times10^{-7}~$d~yr$^{-1}$. It implies that a rate of mass transfer from the secondary component to the primary one is $\mathrm{d}m_2\slash \mathrm{d}t~=~$6.31$\times10^{-8}~M_{\odot}$~yr$^{-1}$. Furthermore, we find that from the detached phase to the contact phase, mass of the evolved secondary component has been lost 1.188$\pm$0.110~$\mathrm{M}_{\odot}$, i.e., mass lost by the system of 0.789$\pm$0.073~$\mathrm{M}_{\odot}$ and mass transfer to the primary of 0.399$\pm$0.037~$\mathrm{M}_{\odot}$. Since the time of the first overflow, the angular momentum loss is found to be 72.2~$\%$ of $J_\mathrm{FOF}$, causing the orbit and Roche surface to shrink until the present time.


title: Segregation--Assisted Spinodal and Transient Spinodal Phase Separation at Grain Boundaries

abstract: Segregation to grain boundaries affects their cohesion, corrosion and embrittlement and plays a critical role in heterogeneous nucleation. In order to quantitatively study segregation and phase separation at grain boundaries, we derive a density-based phase-field model. In this model, we describe the grain boundary free energy based on available bulk thermodynamic data while an atomic grain boundary density is obtained using atomistic simulations. To benchmark the performance of the model, we study Mn grain boundary segregation in the Fe--Mn system. 3D simulation results are compared against atom probe tomography measurements. We show that a continuous increase in the alloy composition results in a discontinuous jump in the Mn grain boundary segregation. This jump corresponds to an interfacial spinodal phase separation. For alloy compositions above the interfacial spinodal, we found a transient spinodal phase separation phenomenon which opens opportunities for knowledge-based microstructure design through the chemical manipulation of grain boundaries. The proposed density-based model provides a powerful tool to study thermodynamics and kinetics of segregation and phase separation at grain boundaries.


title: $^{42}$Ca and $^{50}$Ca with the (Many- and Few-body) Unified Method

abstract: A new method unifying many and few-body aspects of nuclear structure has recently been introduced \cite{hov18}. This method combines the many-body description of a core and the few-body structure of this core surrounded by two valence nucleons. For this reason this method is expected to work specially well when applied to nuclei close to the driplines, where the few-body halo structure with one or more nucleons outside the core is established. In this work we apply the new method to nuclei close to the valley of stability, with $^{42}$Ca and $^{50}$Ca as illustrations. We compare the results from uncorrelated mean-field calculations with the ones obtained with the unified method allowing arbitrary correlations in the valence space. We find that the unified method provides results rather similar, although distinguishable, to the Hartree-Fock calculations. The correlations are much less pronounced than at the driplines, which initially were targets for the unified method. The halo structure is not artificially maintained, but the correlations are here demonstrated to be applicable to well-bound nuclei. Excited states built on valence degrees of freedom are calculated for the same nuclei.


title: Extension of the King-Hele orbit contraction method for accurate, semi-analytical propagation of non-circular orbits

abstract: Numerical integration of orbit trajectories for a large number of initial conditions and for long time spans is computationally expensive. Semi-analytical methods were developed to reduce the computational burden. An elegant and widely used method of semi-analytically integrating trajectories of objects subject to atmospheric drag was proposed by King-Hele (KH). However, the analytical KH contraction method relies on the assumption that the atmosphere density decays strictly exponentially with altitude. If the actual density profile does not satisfy the assumption of a fixed scale height, as is the case for Earth's atmosphere, the KH method introduces potentially large errors for non-circular orbit configurations. In this work, the KH method is extended to account for such errors by using a newly introduced atmosphere model derivative. By superimposing exponentially decaying partial atmospheres, the superimposed KH method can be applied accurately while considering more complex density profiles. The KH method is further refined by deriving higher order terms during the series expansion. A variable boundary condition to choose the appropriate eccentricity regime, based on the series truncation errors, is introduced. The accuracy of the extended analytical contraction method is shown to be comparable to numerical Gauss-Legendre quadrature. Propagation using the proposed method compares well against non-averaged integration of the dynamics, while the computational load remains very low.


title: Fusion hierarchies, $T$-systems and $Y$-systems for the dilute $A_2^{(2)}$ loop models

abstract: The fusion hierarchy, $T$-system and $Y$-system of functional equations are the key to integrability for 2d lattice models. We derive these equations for the generic dilute $A_2^{(2)}$ loop models. The fused transfer matrices are associated with nodes of the infinite dominant integral weight lattice of $s\ell(3)$. For generic values of the crossing parameter $\lambda$, the $T$- and $Y$-systems do not truncate. For the case $\frac{\lambda}{\pi}=\frac{(2p'-p)}{4p'}$ rational so that $x=\mathrm{e}^{\mathrm{i}\lambda}$ is a root of unity, we find explicit closure relations and derive closed finite $T$- and $Y$-systems. The TBA diagrams of the $Y$-systems and associated Thermodynamic Bethe Ansatz (TBA) integral equations are not of simple Dynkin type. They involve $p'+2$ nodes if $p$ is even and $2p'+2$ nodes if $p$ is odd and are related to the TBA diagrams of $A_2^{(1)}$ models at roots of unity by a ${\Bbb Z}_2$ folding which originates from the addition of crossing symmetry. In an appropriate regime, the known central charges are $c=1-\frac{6(p-p')^2}{pp'}$. Prototypical examples of the $A_2^{(2)}$ loop models, at roots of unity, include critical dense polymers ${\cal DLM}(1,2)$ with central charge $c=-2$, $\lambda=\frac{3\pi}{8}$ and loop fugacity $\beta=0$ and critical site percolation on the triangular lattice ${\cal DLM}(2,3)$ with $c=0$, $\lambda=\frac{\pi}{3}$ and $\beta=1$. Solving the TBA equations for the conformal data will determine whether these models lie in the same universality classes as their $A_1^{(1)}$ counterparts. More specifically, it will confirm the extent to which bond and site percolation lie in the same universality class as logarithmic conformal field theories.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: A large data theory for nonlinear wave on the Schwarzschild background

abstract: We study both of the scattering and Cauchy problems for the semilinear wave equation with null quadratic form on the Schwarzschild background. Prescribing the scattering data that are given by the short pulse data on the future null infinity and are trivial on the future event horizon, we construct a class of globally smooth solutions backwards up to any finite time and show that the wave travels in such a way that almost all of the (large) energy is focusing in an outgoing null strip, while little radiates out of this strip. In reverse, considering a class of Cauchy data with large energy norms, there exists a unique and global solution in the future development. And most of the wave packet is confined in an incoming null strip and reflected to the future event horizon, whereas little is transmitted to the future null infinity.


title: H.E.S.S. and Suzaku observations of the Vela X pulsar wind nebula

abstract: Pulsar Wind Nebulae (PWNe) represent the most prominent population of Galactic very-high-energy gamma-ray sources and are thought to be an efficient source of leptonic cosmic rays. Vela X is a nearby middle-aged PWN, which shows bright X-ray and TeV gamma-ray emission toward an elongated structure called the cocoon. Since TeV emission is likely inverse-Compton emission of electrons while X-ray emission is synchrotron radiation of the same electrons, we aim to derive the properties of the relativistic particles and of magnetic fields with minimal modelling. We use data from the Suzaku XIS to derive the spectra from three compact regions in Vela X covering distances from 0.3 pc to 4 pc from the pulsar along the cocoon. We obtain gamma-ray spectra of the same regions from H.E.S.S. observations and fit a radiative model to the multi-wavelength spectra. The TeV electron spectra and magnetic field strengths are consistent within the uncertainties for the three regions, with energy densities of the order $10^{-12}\rm\,erg\,cm^{-3}$. The data indicate the presence of a cutoff in the electron spectrum at energies of $\sim$100 TeV and a magnetic field strength of $\sim$$6\,\rm\mu G$. Constraints on the presence of turbulent magnetic fields are weak. The pressure of TeV electrons and magnetic fields in the cocoon is dynamically negligible, requiring the presence of another dominant pressure component. Sub-TeV electrons cannot account completely for the missing pressure, that may be provided either by relativistic ions or from mixing of the ejecta with the pulsar wind. The electron spectra are consistent with expectations from transport scenarios dominated either by advection via the reverse shock or by diffusion. Constraints on turbulent magnetic fields and the shape of the electron cutoff can be improved by spectral measurements in the energy range $\gtrsim 10\rm\,keV$. (abridged)


title: Stratified sampling and bootstrapping for approximate Bayesian computation

abstract: Approximate Bayesian computation (ABC) is computationally intensive for complex model simulators. To exploit expensive simulations, data-resampling via bootstrapping can be employed to obtain many artificial datasets at little cost. However, when using this approach within ABC, the posterior variance is inflated, thus resulting in biased posterior inference. Here we use stratified Monte Carlo to considerably reduce the bias induced by data resampling. We also show empirically that it is possible to obtain reliable inference using a larger than usual ABC threshold. Finally, we show that with stratified Monte Carlo we obtain a less variable ABC likelihood. Ultimately we show how our approach improves the computational efficiency of the ABC samplers. We construct several ABC samplers employing our methodology, such as rejection and importance ABC samplers, and ABC-MCMC samplers. We consider simulation studies for static (Gaussian, g-and-k distribution, Ising model, astronomical model) and dynamic models (Lotka-Volterra). We compare against state-of-art sequential Monte Carlo ABC samplers, synthetic likelihoods, and likelihood-free Bayesian optimization. For a computationally expensive Lotka-Volterra case study, we found that our strategy leads to a more than 10-fold computational saving, compared to a sampler that does not use our novel approach.


title: Families of two-dimensional Coulomb gases on an ellipse: correlation functions and universality

abstract: We investigate a one-parameter family of Coulomb gases in two dimensions, which are confined to an ellipse, due to a hard wall constraint, and are subject to an additional external potential. At inverse temperature $\beta=2$ we can use the technique of planar orthogonal polynomials, borrowed from random matrix theory, to explicitly determine all $k$-point correlation functions for a fixed number of particles $N$. These are given by the determinant of the kernel of the corresponding orthogonal polynomials, which in our case are the Gegenbauer polynomials, or a subset of the asymmetric Jacobi polynomials, depending on the choice of external potential, as shown in a companion paper recently published by three of the authors. In the rotationally invariant case, when the ellipse becomes the unit disc, our findings agree with that of the ensemble of truncated unitary random matrices. The thermodynamical large-$N$ limit is investigated in the local scaling regime in the bulk and at the edge of the spectrum at weak and strong non-Hermiticity. We find new universality classes in these limits and recover the sine- and Bessel-kernel in the Hermitian limit. The limiting global correlation functions of particles in the interior of the ellipse are more difficult to obtain but found in the special cases corresponding to the Chebyshev polynomials.


title: Strong quantum entanglement via a controllable four wave mixing mechanism in an optomechanical system

abstract: We propose an approach to generate strong quantum entanglement by the controllable four wave mixing mechanism in a single cavity, weak coupling optomechanical system. The optomechanical system is driven by a strong two tone pump field and a weak signal field, simultaneously. The two tone pump field consists of a lower and an upper sideband, which couple with the optical cavity and mechanical resonator, and generate the beam splitter and two mode squeezing interactions under the rotating wave approximation. This interaction mechanism modifies the effective susceptibility of the optomechanical cavity and optomechanically induces a four wave mixing process. Strong quantum entanglement can be generated between the signal and four wave mixing fields with an entangled degree over 16 dB in realistic optomechanical systems. The generation scheme of the quantum entanglement is quite robust against thermal mechanical noise, and entanglement above 3 dB can persist at room temperature in the weak coupling regime.


title: Microscopic three-cluster study of light exotic nuclei

abstract: I develop a microscopic three-cluster model for exotic light nuclei. I use the hyperspherical formalism, associated with the Generator Coordinate Method. This model is well adapted to halo nuclei, since the long-range part of the radial wave functions is accurately reproduced. The core wave functions are described in the shell model, including excited states. This technique provides large bases, expressed in terms of projected Slater determinants. Matrix elements involve seven-dimension integrals, and therefore require long calculation times. I apply the model to 11Li, 14Be, 15B, and 17N described by two neutrons surrounding a 9Li, 12Be, 13B and 15N core, respectively. The 17Ne (as 15O+p+p) and 15Ne (as 13O+p+p) mirror nuclei are briefly discussed. I present the spectra and some spectroscopic properties, such as r.m.s. radii or E2 transition probabilities. I also analyze the importance of core excitations.


title: Sharp Bounds for Oscillatory Integral Operators with Homogeneous Polynomial Phases

abstract: We obtain sharp $L^p$ bounds for oscillatory integral operators with generic homogeneous polynomial phases in several variables. The phases considered in this paper satisfy the rank one condition which is an important notion introduced by Greenleaf, Pramanik and Tang. Under certain additional assumptions, we can establish sharp damping estimates with critical exponents to prove endpoint $L^p$ estimates.


title: Burst particle creation in gravitational collapse to a horizonless compact object

abstract: In the previous paper [Harada, Cardoso, and Miyata, Phys.\ Rev.\ D {\bf 99} (2019), 044039], it is shown that a hollow transmissive shell collapsing to an ultracompact object of radius very close to its horizon radius generally emits transient Hawking radiation followed by a couple of bursts separated each other by a long time interval. In the current paper, we expand the previous work in two independent directions: changing boundary conditions and specifying the equations of state of the matter. First, we introduce a perfectly reflective surface collapsing to an ultracompact object and find that this model also emits transient Hawking radiation that is followed only by a single burst. Second, we introduce two different collapse dynamics to an ultracompact object and specify the corresponding matter equations of state. We find that transient Hawking radiation is quite commonly seen in early times, while the subsequent bursts strongly depend on the boundary condition and the equation of state or the braking behavior of the surface.


title: Search for opposite sign muon-tau pair and a b-jet at LHC in the context of flavor anomalies

abstract: Extant anomalies in several semileptonic $B$-meson decays argue for physics beyond the Standard Model. Measurements of both neutral-current decays (such as $R_K$, $R_{K^*}$ and $B_s\rightarrow \phi\mu\mu$) as well as charged-current ones- $R(D)$ and $R(D^*)$ -provide strong hints for the violation of lepton flavor universality. Recent studies have shown that a class of effective field theory (EFT) models may explain such anomalies in terms of only a few parameters which can be determined phenomenologically. In this literature, we examine such resolutions in the context of the requisite $(\bar{s}$ $b$)($\bar\tau\ \tau)$ operator, and look for its signals at the 13 TeV LHC, with a final state of one $b$-jet, and an oppositely charged $\mu$-$\tau$ pair, with the muon coming from the decay of one of the $\tau$ leptons. We obtain discovery and exclusion limits on the model parameters as a function of luminosity at the 13 TeV LHC.


title: Limit theorems for Jacobi ensembles with large parameters

abstract: Consider Jacobi random matrix ensembles with the distributions $$c_{k_1,k_2,k_3}\prod_{1\leq i< j \leq N}\left(x_j-x_i\right)^{k_3}\prod_{i=1}^N \left(1-x_i\right)^{\frac{k_1+k_2}{2}-\frac{1}{2}}\left(1+x_i\right)^{\frac{k_2}{2}-\frac{1}{2}} dx$$ of the eigenvalues on the alcoves $$A:=\{x\in\mathbb R^N| \> -1\leq x_1\le ...\le x_N\leq 1\}.$$ For $(k_1,k_2,k_3)=\kappa\cdot (a,b,1)$ with $a,b>0$ fixed, we derive a central limit theorem for the distributions above for $\kappa\to\infty$. The drift and the inverse of the limit covariance matrix are expressed in terms of the zeros of classical Jacobi polynomials. We also rewrite the CLT in trigonometric form and determine the eigenvalues and eigenvectors of the limit covariance matrices. These results are related to corresponding limits for $\beta$-Hermite and $\beta$-Laguerre ensembles for $\beta\to\infty$ by Dumitriu and Edelman and by Voit.


title: Are all the frames equally important?

abstract: In this work, we address the problem of measuring and predicting temporal video saliency - a metric which defines the importance of a video frame for human attention. Unlike the conventional spatial saliency which defines the location of the salient regions within a frame (as it is done for still images), temporal saliency considers importance of a frame as a whole and may not exist apart from context. The proposed interface is an interactive cursor-based algorithm for collecting experimental data about temporal saliency. We collect the first human responses and perform their analysis. As a result, we show that qualitatively, the produced scores have very explicit meaning of the semantic changes in a frame, while quantitatively being highly correlated between all the observers. Apart from that, we show that the proposed tool can simultaneously collect fixations similar to the ones produced by eye-tracker in a more affordable way. Further, this approach may be used for creation of first temporal saliency datasets which will allow training computational predictive algorithms. The proposed interface does not rely on any special equipment, which allows to run it remotely and cover a wide audience.


title: Spontaneous emission of color centers at 4eV in hexagonal boron nitride under hydrostatic pressure

abstract: The light emission properties of color centers emitting in 3.3-4 eV region are investigated for hydrostatic pressures ranging up to 5GPa at liquid helium temperature. The light emission energy decreases with pressure less sensitively than the bandgap. This behavior at variance from the shift of the bandgap is typical of deep traps. Interestingly, hydrostatic pressure reveals the existence of levels that vary differently under pressure (smaller increase of the emission wavelength compared to the rest of the levels in this energy region or even decrease of it) with pressure. This discovery enriches the physics of the color centers operating in the UV in hBN.


title: Robust Online Algorithms for Dynamic Problems

abstract: Online algorithms that allow a small amount of migration or recourse have been intensively studied in the last years. They are essential in the design of competitive algorithms for dynamic problems, where objects can also depart from the instance. In this work, we give a general framework to obtain so called robust online algorithms for these dynamic problems: these online algorithm achieve an asymptotic competitive ratio of $\gamma+\epsilon$ with migration $O(1/\epsilon)$, where $\gamma$ is the best known offline asymptotic approximation ratio. In order to use our framework, one only needs to construct a suitable online algorithm for the static online case, where items never depart. To show the usefulness of our approach, we improve upon the best known robust algorithms for the dynamic versions of generalizations of Strip Packing and Bin Packing, including the first robust algorithm for general $d$-dimensional Bin Packing and Vector Packing.


title: ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

abstract: ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are useful in a wide range of application areas. Applications powered by ExaHyPE can be run on a student's laptop, but are also able to exploit thousands of processor cores on state-of-the-art supercomputers. The engine is able to dynamically increase the accuracy of the simulation using adaptive mesh refinement where required. Due to the robustness and shock capturing abilities of ExaHyPE's numerical methods, users of the engine can simulate linear and non-linear hyperbolic PDEs with very high accuracy. Users can tailor the engine to their particular PDE by specifying evolved quantities, fluxes, and source terms. A complete simulation code for a new hyperbolic PDE can often be realised within a few hours - a task that, traditionally, can take weeks, months, often years for researchers starting from scratch. In this paper, we showcase ExaHyPE's workflow and capabilities through real-world scenarios from our two main application areas: seismology and astrophysics.


title: Sharp estimate on the inner distance in planar domains

abstract: We show that the inner distance inside a bounded planar domain is at most the one-dimensional Hausdorff measure of the boundary of the domain. We prove this sharp result by establishing an improved Painlev\'e length estimate for connected sets and by using the metric removability of totally disconnected sets, proven by Kalmykov, Kovalev, and Rajala. We also give a totally disconnected example showing that for general sets the Painlev\'e length bound $\kappa(E) \le\pi \mathcal{H}^1(E)$ is sharp.


title: Topological Classification of Excitations in Quadratic Bosonic Systems

abstract: We investigate the topological classification of excitations in quadratic bosonic systems with an excitation band gap. Time-reversal, charge-conjugation, and parity symmetries in bosonic systems are introduced to realize a ten-fold symmetry classification. We find a specific decomposition of the quadratic bosonic Hamiltonian and use it to prove that each quadratic bosonic system is homotopic to a direct sum of two single-particle subsystems. The topological classification table is thus derived via inheriting from that of Atland-Zirnbauer classes and unique topological phases of bosons are predicted. Finally, concrete topological models are proposed to demonstrate the peculiarity of bosonic excitations.


title: Stochastic thermodynamics with arbitrary interventions

abstract: We extend the theory of stochastic thermodynamics in three directions: (i) instead of a continuously monitored system we consider measurements only at an arbitrary set of discrete times, (ii) we allow for imperfect measurements and incomplete information in the description, and (iii) we treat arbitrary manipulations (e.g. feedback control operations) which are allowed to depend on the entire measurement record. For this purpose we define for a driven system in contact with a single heat bath the four key thermodynamic quantities: internal energy, heat, work and entropy along a single `trajectory' for a causal model. The first law at the trajectory level and the second law on average is verified. We highlight the special case of Bayesian or `bare' measurements (incomplete information, but no average disturbance) which allows us to compare our theory with the literature and to derive a general inequality for the estimated free energy difference in Jarzynksi-type experiments. An analysis of a recent Maxwell demon experiment using real-time feedback control is also given. As a mathematical tool, we prove a classical version of Stinespring's dilation theorem, which might be of independent interest.


title: Deep Transfer Learning Methods for Colon Cancer Classification in Confocal Laser Microscopy Images

abstract: Purpose: The gold standard for colorectal cancer metastases detection in the peritoneum is histological evaluation of a removed tissue sample. For feedback during interventions, real-time in-vivo imaging with confocal laser microscopy has been proposed for differentiation of benign and malignant tissue by manual expert evaluation. Automatic image classification could improve the surgical workflow further by providing immediate feedback. Methods: We analyze the feasibility of classifying tissue from confocal laser microscopy in the colon and peritoneum. For this purpose, we adopt both classical and state-of-the-art convolutional neural networks to directly learn from the images. As the available dataset is small, we investigate several transfer learning strategies including partial freezing variants and full fine-tuning. We address the distinction of different tissue types, as well as benign and malignant tissue. Results: We present a thorough analysis of transfer learning strategies for colorectal cancer with confocal laser microscopy. In the peritoneum, metastases are classified with an AUC of 97.1 and in the colon, the primarius is classified with an AUC of 73.1. In general, transfer learning substantially improves performance over training from scratch. We find that the optimal transfer learning strategy differs for models and classification tasks. Conclusions: We demonstrate that convolutional neural networks and transfer learning can be used to identify cancer tissue with confocal laser microscopy. We show that there is no generally optimal transfer learning strategy and model as well as task-specific engineering is required. Given the high performance for the peritoneum, even with a small dataset, application for intraoperative decision support could be feasible.


title: Solving Poisson equations by the MN-curve approach

abstract: In this paper we apply the newly born choice theory of the shape parameters contained in the smooth radial basis functions to solve Poisson equations. Some people complain that Luh's choice theory, based on harmonic analysis, is mathematically complicated and applies only to function interpolations. Here we aim at presenting an easily accessible approach to solving differential equations with the choice theory which proves to be successful, not only by its easy accessibility, but also by its striking accuracy and efficiency.


title: Dissipation and Semigroup on $H^k_n$: Non-cutoff Linearized Boltzmann Operator with Soft Potential

abstract: In this paper, we find that the linearized collision operator $L$ of the non-cutoff Boltzmann equation with soft potential generates a strongly continuous semigroup on $H^k_n$, with $k,n\in\mathbb{R}$. In the theory of Boltzmann equation without angular cutoff, the weighted Sobolev space plays a fundamental role. The proof is based on pseudo-differential calculus and in general, for a specific class of Weyl quantization, the $L^2$ dissipation implies $H^k_n$ dissipation. This kind of estimate is also known as the G{\aa}rding's inequality.


title: Investigating an angular correlation between nearby starburst galaxies and UHECRs with the Telescope Array experiment

abstract: The arrival directions of cosmic rays detected by the Pierre Auger Observatory (Auger) with energies above 39~EeV were recently reported to correlate with the positions of 23 nearby starburst galaxies (SBGs): in their best-fit model, 9.7\% of the cosmic-ray flux originates from these objects and undergoes angular diffusion on a $12.9^\circ$~scale. On the other hand, some of the SBGs on their list, including the brightest one (M82), are at northern declinations outside the Auger field of view. Data from detectors in the northern hemisphere would be needed to look for cosmic-ray excesses near these objects. In this work, we tested the Auger best-fit model against data collected by the Telescope Array (TA) in a 9-year period, without trying to re-optimize the model parameters for our dataset in order not to introduce statistical penalties. The resulting test statistic (double log-likelihood ratio) was $-1.00$, corresponding to $1.1\sigma$ significance among isotropically generated random datasets, and to $-1.4\sigma$ significance among ones generated assuming the Auger best-fit model. In other words, our data is still insufficient to conclusively rule out either hypothesis. The ongoing fourfold expansion of TA will collect northern hemisphere data with much more statistics, improving our ability to discriminate between different flux models.


title: Transient hydrophobic exposure in the molecular dynamics of Abeta peptide at low water concentration

abstract: Abeta is a disordered peptide central to Alzheimer's Disease. Aggregation of Abeta has been widely explored, but its molecular crowding less so. The synaptic cleft where Abeta locates only holds 60-70 water molecules along its width. We subjected Abeta40 to 100 different simulations with variable water cell size. We show that even for this disordered aggregation-prone peptide, many properties are not cell-size dependent, i.e. a small cell is easily justified. The radius of gyration, intra-peptide, and peptide-water hydrogen bonds are well-sampled by short (50 ns) time scales at any cell size. Abeta is mainly disordered with 0-30% alpha helix but undergoes consistent alpha-beta transitions up to 14% strand in 5-10% of the simulations regardless of cell size. The similar prevalence in long and short simulations indicate small diffusion barriers for structural transitions in contrast to folded globular proteins, which we suggest is a defining hallmark of intrinsically disordered proteins. Importantly, the hydrophobic surface increases significantly in small cells (confidence level 95%, two-tailed t-test), as does the variation in exposure and backbone conformations (>40% and >27% increased standard deviations). Whereas hydrophilic exposure dominates hydrophobic exposure in large cells, this tendency breaks down at low water concentration. We interpret these findings as a concentration-dependent hydrophobic effect, with the small water layer unable to keep the protein unexposed, an effect mainly caused by the layered water-water interactions, not by the peptide dynamics. The exposure correlates with radius of gyration (R2 0.35-0.50) and could be important in crowded environments, e.g. the synaptic cleft.


title: A Linearly Convergent Proximal Gradient Algorithm for Decentralized Optimization

abstract: Decentralized optimization is a powerful paradigm that finds applications in engineering and learning design. This work studies decentralized composite optimization problems with non-smooth regularization terms. Most existing gradient-based proximal decentralized methods are known to converge to the optimal solution with sublinear rates, and it remains unclear whether this family of methods can achieve global linear convergence. To tackle this problem, this work assumes the non-smooth regularization term is common across all networked agents, which is the case for many machine learning problems. Under this condition, we design a proximal gradient decentralized algorithm whose fixed point coincides with the desired minimizer. We then provide a concise proof that establishes its linear convergence. In the absence of the non-smooth term, our analysis technique covers the well known EXTRA algorithm and provides useful bounds on the convergence rate and step-size.


title: A Comparative Review of Microservices and Monolithic Architectures

abstract: Microservices' architecture is getting attention in the academic community and the industry, and mostly is compared with monolithic architecture. Plenty of the results of these research papers contradict each other regarding the performance of these architectures. Therefore, these two architectures are compared in this paper, and some specific configurations of microservices' applications are evaluated as well in the term of service discovery. Monolithic architecture in concurrency testing showed better performance in throughput by 6% when compared to microservices architecture. The load testing scenario did not present significant difference between the two architectures. Furthermore, a third test comparing microservices applications built with different service discovery technologies such as Consul and Eureka showed that applications with Consul presented better results in terms of throughput.


title: Method for Assessing the Fidelity of Optical Diffraction Tomography Reconstruction Methods

abstract: We use a spatial light modulator in a diffraction tomographic system to assess the accuracy of different refractive index reconstruction algorithms. Optical phase conjugation principles through complex media, allows us to quantify the error for different refractive index reconstruction algorithms without access to the ground truth. To our knowledge, this is the first assessment technique that uses structured illumination experimentally to test the accuracy of different reconstruction schemes.


title: PymePix: A python library for SPIDR readout of Timepix3

abstract: PymePix is a new Python 3 library that provides control and acquisition for the Timepix3-SPIDR hardware. The rich set of data-structures and intuitive routines reduces time and coding effort to quickly configure, acquire, and visualize data from Timepix3. The highly extensible high-performance data-pipeline allows for alteration of the Timepix3 datastream into a form that is convinient for the user. This library is intended to be easily inserted into a standard scientific software stack as well as to allow for more direct interaction of Timepix3 with interactive flavors of Python. Included with the library are two example programs using PymePix: pymepix-acq is a command line control and acquisition program that can capture UDP packets and decode them into pixels and triggers. The second is PymePix-Viewer, an online control and data-acquisition program for general use, but with features geared toward mass-spectroscopy and ion imaging.


title: Generating Hypergraphs, Decomposability and Classification of Two-Step Nilpotent Lie Algebras

abstract: In 1973, Gauger proposed a generator-relation method and a duality theory for two-step nilpotent Lie algebras. Based upon these, he classified two-step nilpotent Lie algebras of dimension $8$. In 1999, Galitski and Timashev continued this approach to classify two-step nilpotent Lie algebras of dimension $9$. Their results were partially improved by Ren and Zhu in 2011, Yan and Deng in 2013. Some decomposable two-step nilpotent Lie algebras were excluded in the case of dimension $8$. In this paper, we define generating hypergraph for a two-step nilpotent Lie algebra. The two-step nilpotent Lie algebra is decomposable if and only if its generating hypergraph is not connected under certain bases. Using this result, we identify some decomposable two-step nilpotent Lie algebras in dimension $9$. We give a direct proof that the five two-step nilpotent Lie algebras for dimension $8$, classified by Ren and Zhu in 2011, are all indecomposable. We also introduce a conventional nomenclature for two-step nilpotent Lie algebras of dimension $n = 8, 9$, classified by Ren, Zhu, Yan and Deng, etc.


title: The extremal number of longer subdivisions

abstract: For a multigraph $F$, the $k$-subdivision of $F$ is the graph obtained by replacing the edges of $F$ with pairwise internally vertex-disjoint paths of length $k+1$. Conlon and Lee conjectured that if $k$ is even, then the $(k-1)$-subdivision of any multigraph has extremal number $O(n^{1+\frac{1}{k}})$, and moreover, that for any simple graph $F$ there exists $\varepsilon>0$ such that the $(k-1)$-subdivision of $F$ has extremal number $O(n^{1+\frac{1}{k}-\varepsilon})$. In this paper, we prove both conjectures.


title: Bounding Selmer groups for the Rankin--Selberg convolution of Coleman families

abstract: Let $f$ and $g$ be two cuspidal modular forms and let $\mathcal{F}$ be a Coleman family passing through $f$, defined over an open affinoid subdomain $V$ of weight space $\mathcal{W}$. Using ideas of Pottharst, under certain hypotheses on $f$ and $g$ we construct a coherent sheaf over $V \times \mathcal{W}$ which interpolates the Bloch-Kato Selmer group of the Rankin-Selberg convolution of two modular forms in the critical range (i.e. the range where the $p$-adic $L$-function $L_p$ interpolates critical values of the global $L$-function). We show that the support of this sheaf is contained in the vanishing locus of $L_p$.


title: Pele's hairs and exotic multiply twinned graphite closed-shell microcrystals in meteoritic dust of Chelyabinsk superbolide

abstract: When a space body enters Earth's atmosphere, its surface is exposed to high pressure and temperatures. The airflow tears off small droplets from the meteoroid forming a cloud of meteorite dust. Can new materials be synthesized in these unique conditions (high temperature, pressure, gaseous atmosphere, catalysts)? As a rule, meteoritic dust dissipates in the atmosphere without a trace or is mixed with terrestrial soil. The Chelyabinsk superbolide, the biggest in the 21st century, which exploded on February 15, 2013 above snowy fields of the Southern Urals, was an exception. Two new types of materials were found during an in-depth study: thread-like structures that were not previously associated with meteorite falls and unique carbon crystals with a size of several micrometers that were not observed before. The nature of formation of thread-like structures is fully similar to the nature of formation of Pele's hair during eruptions of terrestrial volcanos. Multiple twin growth mechanism of formation of closed shell graphite microcrystals was proposed based on DFT and classical/ab initio MD simulations. It was found that among several possible embryo carbon nanoclusters the C60 fullerene and polyhexacyclooctadecane -C18H12- may be the main suspects responsible for the formation of closed shell quasi-spherical and hexagonal rod graphite microcrystals.


title: Risk-Sensitive Credit Portfolio Optimization under Partial Information and Contagion Risk

abstract: This paper investigates the finite horizon risk-sensitive portfolio optimization in a regime-switching credit market with physical and information-induced default contagion. It is assumed that the underlying regime-switching process has countable states and is unobservable. The stochastic control problem is formulated under partial observations of asset prices and sequential default events. By establishing a martingale representation theorem based on incomplete and phasing out filtration, we connect the control problem to a quadratic BSDE with jumps, in which the driver term is non-standard and carries the conditional filter as an infinite-dimensional parameter. By proposing some truncation techniques and proving a uniform a priori estimates, we obtain the existence of a solution to the BSDE using the convergence of solutions associated to some truncated BSDEs. The verification theorem can be concluded with the aid of our BSDE results, which in turn yields the uniqueness of the solution to the BSDE.


title: Well-Posedness of Sparse Frequency Estimation

abstract: The problem of estimating the frequencies of an exponential sum has been studied extensively over the last years. It can be understood as a sparse estimation problem, as it strives to identify the sparse representation of a signal using exponentials. In this paper, we are interested in its intrinsic stability properties. We derive a bound very similar to the restricted isometry property. Conditional well-posedness follows: Any exponential sum with samples close to the unknown ground truth has close frequencies as well, provided that it satisfies our model assumptions. The most important assumption is that the frequencies are well-separated. Furthermore, we show that the presented bound is sharp and gives rise to improved estimates of condition numbers of certain Vandermonde matrices.


title: Deep Reinforcement Learning Based Parameter Control in Differential Evolution

abstract: Adaptive Operator Selection (AOS) is an approach that controls discrete parameters of an Evolutionary Algorithm (EA) during the run. In this paper, we propose an AOS method based on Double Deep Q-Learning (DDQN), a Deep Reinforcement Learning method, to control the mutation strategies of Differential Evolution (DE). The application of DDQN to DE requires two phases. First, a neural network is trained offline by collecting data about the DE state and the benefit (reward) of applying each mutation strategy during multiple runs of DE tackling benchmark functions. We define the DE state as the combination of 99 different features and we analyze three alternative reward functions. Second, when DDQN is applied as a parameter controller within DE to a different test set of benchmark functions, DDQN uses the trained neural network to predict which mutation strategy should be applied to each parent at each generation according to the DE state. Benchmark functions for training and testing are taken from the CEC2005 benchmark with dimensions 10 and 30. We compare the results of the proposed DE-DDQN algorithm to several baseline DE algorithms using no online selection, random selection and other AOS methods, and also to the two winners of the CEC2005 competition. The results show that DE-DDQN outperforms the non-adaptive methods for all functions in the test set; while its results are comparable with the last two algorithms.


title: Numerical calculation of 5-loop QED contributions to the electron anomalous magnetic moment

abstract: Preliminary numerical results of calculating the 5-loop QED contributions to the electron (and muon) anomalous magnetic moment are presented. The results include the total contribution of the 5-loop Feynman diagrams without lepton loops and the contributions of nine gauge-invariant classes that form that set. A discrepancy with known results is revealed. The contributions of the gauge-invariant classes are presented for the first time. The method of the computation is briefly described (with the corresponding references). The calculation is based on the following elements: 1) a subtraction procedure for removing infrared (IR) and ultraviolet (UV) divergences point by point in Feynman parametric representation before integration; 2) a nonadaptive Monte Carlo integration method that is founded on probability density functions (PDF) that are constructed for each Feynman diagram individually using its combinatorial structure; 3) a GPU-based numerical integration with the help of the supercomputer "Govorun" (JINR, Dubna, Russia).


title: Less Memory, Faster Speed: Refining Self-Attention Module for Image Reconstruction

abstract: Self-attention (SA) mechanisms can capture effectively global dependencies in deep neural networks, and have been applied to natural language processing and image processing successfully. However, SA modules for image reconstruction have high time and space complexity, which restrict their applications to higher-resolution images. In this paper, we refine the SA module in self-attention generative adversarial networks (SAGAN) via adapting a non-local operation, revising the connectivity among the units in SA module and re-implementing its computational pattern, such that its time and space complexity is reduced from $\text{O}(n^2)$ to $\text{O}(n)$, but it is still equivalent to the original SA module. Further, we explore the principles behind the module and discover that our module is a special kind of channel attention mechanisms. Experimental results based on two benchmark datasets of image reconstruction, verify that under the same computational environment, two models can achieve comparable effectiveness for image reconstruction, but the proposed one runs faster and takes up less memory space.


title: Operator norm and numerical radius analogues of Cohen's inequality

abstract: Let $D$ be an invertible multiplication operator on $L^2(X, \mu)$, and let $A$ be a bounded operator on $L^2(X, \mu)$. In this note we prove that $\|A\|^2 \le \|D A\| \, \|D^{-1} A\|$, where $\|\cdot\|$ denotes the operator norm. If, in addition, the operators $A$ and $D$ are positive, we also have $w(A)^2 \le w(D A) \, w(D^{-1} A)$, where $w$ denotes the numerical radius.


title: Schr\"odinger cats and steady states in subharmonic generation with Kerr nonlinearities

abstract: We discuss general properties of the equilibrium state of parametric down-conversion in superconducting quantum circuits with detunings and Kerr anharmonicities, in the strongly nonlinear regime. By comparing moments of the steady state and those of a Schr\"odinger cat, we show that true Schr\"odinger cats cannot survive in the steady state if there is any single-photon loss. A delta-function 'cat-like' steady-state distribution can be formed, but this only exists in the limit of an extremely large nonlinearity. The steady state is a mixed state, which is more complex than a mixture or linear combination of delta-functions, and whose purity is reduced by driving. We expect this general behaviour to occur in other driven, dissipative quantum subharmonic non-equilibrium open systems.


title: Integrated conditional moment test and beyond: when the number of covariates is divergent

abstract: The classic integrated conditional moment test is a promising method for testing regression model misspecification. However, it severely suffers from the curse of dimensionality. To extend it to handle the testing problem for parametric multi-index models with diverging number of covariates, we investigate three issues in inference in this paper. First, we study the consistency and asymptotically linear representation of the least squares estimator of the parameter matrix at faster rates of divergence than those in the literature for nonlinear models. Second, we propose, via sufficient dimension reduction techniques, an adaptive-to-model version of the integrated conditional moment test. We study the asymptotic properties of the new test under both the null and alternative hypothesis to examine its ability of significance level maintenance and its sensitivity to the global and local alternatives that are distinct from the null at the fastest possible rate in hypothesis testing. Third, we derive the consistency of the bootstrap approximation for the new test in the diverging dimension setting. The numerical studies show that the new test can very much enhance the performance of the original ICM test in high-dimensional scenarios. We also apply the test to a real data set for illustrations.


title: Testing general relativity with supermassive black holes using X-ray reflection spectroscopy

abstract: In this paper, we review our current efforts to test General Relativity in the strong field regime by studying the reflection spectrum of supermassive black holes. So far we have analyzed 11 sources with observations of NuSTAR, Suzaku, Swift, and XMM-Newton. Our results are consistent with general relativity, according to which the spacetime metric around astrophysical black holes should be well approximated by the Kerr solution. We discuss the systematic uncertainties in our model and we present a preliminary study on the impact of some of them on the measurement of the spacetime metric.


title: Matrix liberation process II: Relation to orbital free entropy

abstract: We investigate the concept of orbital free entropy from the viewpoint of matrix liberation process. We will show that many basic questions around the definition of orbital free entropy are reduced to the question of full large deviation principle for the matrix liberation process. We will also obtain a large deviation upper bound for a certain family of random matrices that is an essential ingredient to define the orbital free entropy. The resulting rate function is made up into a new approach to free mutual information.


title: Stability of the Couette flow under the 2D steady Navier-Stokes flow

abstract: In this note, we investigate the stability property of shear flows under the 2D stationary Navier-Stokes equations, and we obtain that the Couette flow $(y,0)$ is stable under the space of $\mathcal{D}^{1,q}(\mathbb{R}^2)$ for any $1<q<\infty$ and unstable in the space of $\mathcal{D}^{1,\infty}(\mathbb{R}^2)$. A key observation is the anisotropic cut-off function. We also consider the Poiseuille flow $(y^2,0)$, which is stable in $\mathcal{D}^{1,q}(\mathbb{R}^2)$ with $\frac43<q\leq4.$


title: Flow and air-entrainment around partially submerged vertical cylinders

abstract: In this study, a partially submerged vertical cylinder is moved at constant velocity through water, which is initially at rest. During the motion, the wake behind the cylinder induces free-surface deformation. Eleven cylinders, with diameters from $D=1.4$ to 16 cm, were tested at two different conditions: (i) constant immersed height $h$ and (ii) constant $h/D$. The range of translation velocities and diameters are in the regime of turbulent wake with experiments carried out for $4500<Re<240 \,000$ and $0.2<Fr<2.4$, where $Re$ and $Fr$ are the Reynolds and Froude numbers based on $D$. The focus here is on drag force measurements and relatively strong free-surface deformation up to air-entrainment. Specifically, two modes of air-entraiment have been uncovered: (i) in the cavity along the cylinder wall and (ii) in the wake of the cylinder. A scaling for the critical velocity for air-entrainment in the cavity has been observed in agreement with a simple model. Furthermore, for $Fr>1.2$, the drag force varies linearly with $Fr$.


title: High-order harmonic generation in gapped bilayer graphene

abstract: Microscopic nonlinear quantum theory of interaction of coherent electromagnetic radiation with gapped bilayer graphene is developed. The Liouville-von Neumann equation for the density matrix is solved numerically at the multiphoton excitation regime. The developed theory of interaction of charged carriers with strong driving wave field is valid near the Dirac points of the Brillouin zone. We consider the harmonic generation process in the nonadiabatic regime of interaction when the Keldysh parameter is of the order of unity. On the basis of numerical solutions, we examine the rates of odd and even high-harmonics at the particle-hole annihilation in the field of a strong pump wave of arbitrary polarization. Obtained results show that the gapped bilayer graphene can serve as an effective medium for generation of even and odd high harmonics in the THz and far infrared domains of frequencies.


title: Occurrence of Great Magnetic Storms on 6-8 March 1582

abstract: Although knowing the occurrence frequency of severe space weather events is important for a modern society, it is insufficiently known due to the lack of magnetic or sunspot observations, before the Carrington event in 1859 known as one of the largest events during the last two centuries. Here, we show that a severe magnetic storm occurred on 8 March 1582 based on auroral records in East Asia. The equatorward boundary of auroral visibility reached 28.8{\deg} magnetic latitude. The equatorward boundary of the auroral oval is estimated to be 33.0{\deg} invariant latitude (ILAT), which is comparable to the storms on 25/26 September 1909 (~31.6{\deg} ILAT, minimum Dst of -595 nT), 28/29 August 1859 (~36.5{\deg} ILAT), and 13/14 March 1989 (~40{\deg} ILAT, minimum Dst of -589 nT). Assuming that the equatorward boundary is a proxy for the scale of magnetic storms, we presume that the storm on March 1582 was severe. We also found that the storm on March 1582 lasted, at least, for three days by combining European records. The auroral oval stayed at mid-latitude for the first two days and moved to low-latitude (in East Asia) for the last day. It is plausible that the storm was caused by a series of ICMEs (interplanetary coronal mass ejections). We can reasonably speculate that a first ICME could have cleaned up interplanetary space to make the following ICMEs more geo-effective, as probably occurred in the Carrington and Halloween storms.


abstract: In this note, we study the relation between Fontaine-Laffaille modules and strongly divisible modules, without assuming the main theorem of Fontaine-Laffaille (but we need to assume the main results concerning strongly divisible modules). This in particular gives a new proof for the main theorem of Fontaine-Laffaille (for $p>2$).


title: FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs

abstract: Complex algebraic calculations can be performed by reconstructing analytic results from numerical evaluations over finite fields. We describe FiniteFlow, a framework for defining and executing numerical algorithms over finite fields and reconstructing multivariate rational functions. The framework employs computational graphs, known as dataflow graphs, to combine basic building blocks into complex algorithms. This allows to easily implement a wide range of methods over finite fields in high-level languages and computer algebra systems, without being concerned with the low-level details of the numerical implementation. This approach sidesteps the appearance of large intermediate expressions and can be massively parallelized. We present applications to the calculation of multi-loop scattering amplitudes, including the reduction via integration-by-parts identities to master integrals or special functions, the computation of differential equations for Feynman integrals, multi-loop integrand reduction, the decomposition of amplitudes into form factors, and the derivation of integrable symbols from a known alphabet. We also release a proof-of-concept C++ implementation of this framework, with a high-level interface in Mathematica.


title: Limit of torsion semi-stable Galois representations with unbounded weights

abstract: Let $K$ be a complete discrete valuation field of characteristic $(0, p)$ with perfect residue field, and let $T$ be an integral $\mathbb{Z}_p$-representation of $\mathrm{Gal}(\overline{K}/K)$. A theorem of T. Liu says that if $T/p^n T$ is torsion semi-stable (resp. crystalline) of uniformly bounded Hodge-Tate weights for all $n \geq 1$, then $T$ is also semi-stable (resp. crystalline). In this note, we show that we can relax the condition of "uniformly bounded Hodge-Tate weights" to an unbounded (log-)growth condition.


title: The TrieJax Architecture: Accelerating Graph Operations Through Relational Joins

abstract: Graph pattern matching (e.g., finding all cycles and cliques) has become an important component in many critical domains such as social networks, biology, and cyber-security. This development motivated research to develop faster algorithms that target graph pattern matching. In recent years, the database community has shown that mapping graph pattern matching problems to a new class of relational join algorithms provides an efficient framework for computing these problems. In this paper, we argue that this new class of relational join algorithms is highly amenable to specialized hardware acceleration thanks to two fundamental properties: improved locality and inherent concurrency. The improved locality is a result of the provably bound number of intermediate results these algorithms generate, which results in smaller working sets. In addition, their inherent concurrency can be leveraged for effective hardware acceleration and hiding memory latency. We demonstrate the hardware amenability of this new class of algorithms by introducing TrieJax, a hardware accelerator for graph pattern matching. The TrieJax design leverages the improved locality and high concurrency properties to dramatically accelerate graph pattern matching, and can be tightly integrated into existing manycore processors. We evaluate TrieJax on a set standard graph pattern matching queries and datasets. Our evaluation shows that TrieJax outperforms recently proposed hardware accelerators for graph and database processing that do not employ the new class of algorithms by 7-63x on average (up to 539x), while consuming 15-179x less energy (up to 1750x). systems that do incorporate modern relational join algorithms by 9-20x on average (up to 45x), while consuming 59-110x less energy (up to 372x).


title: An iterative scheme for feature based positioning using a weighted dissimilarity measure

abstract: We propose an iterative scheme for feature-based positioning using a new weighted dissimilarity measure with the goal of reducing the impact of large errors among the measured or modeled features. The weights are computed from the location-dependent standard deviations of the features and stored as part of the reference fingerprint map (RFM). Spatial filtering and kernel smoothing of the kinematically collected raw data allow efficiently estimating the standard deviations during RFM generation. In the positioning stage, the weights control the contribution of each feature to the dissimilarity measure, which in turn quantifies the difference between the set of online measured features and the fingerprints stored in the RFM. Features with little variability contribute more to the estimated position than features with high variability. Iterations are necessary because the variability depends on the location, and the location is initially unknown when estimating the position. Using real WiFi signal strength data from extended test measurements with ground truth in an office building, we show that the standard deviations of these features vary considerably within the region of interest and are neither simple functions of the signal strength nor of the distances from the corresponding access points. This is the motivation to include the empirical standard deviations in the RFM. We then analyze the deviations of the estimated positions with and without the location-dependent weighting. In the present example the maximum radial positioning error from ground truth are reduced by 40% comparing to kNN without the weighted dissimilarity measure.


title: Fluid pumping and active flexoelectricity can promote lumen nucleation in cell assemblies

abstract: We discuss the physical mechanisms that promote or suppress the nucleation of a fluid-filled lumen inside a cell assembly or a tissue. We discuss lumen formation in a continuum theory of tissue material properties in which the tissue is described as a two-fluid system to account for its permeation by the interstitial fluid, and we include fluid pumping as well as active electric effects. Considering a spherical geometry and a polarized tissue, our work shows that fluid pumping and tissue flexoelectricity play a crucial role in lumen formation. We furthermore explore the large variety of long-time states that are accessible for the cell aggregate and its lumen. Our work reveals a role of the coupling of mechanical, electrical and hydraulic phenomena in tissue lumen formation.