title: Exciting out-of-equilibrium states in macromolecules through light pumping

abstract: In the present paper we address the problem of the energy downconversion of the light absorbed by a protein into its internal vibrational modes. We consider the case in which the light receptors are fluorophores either naturally co-expressed with the protein or artificially covalently bound to some of its amino acids. In a recent work [Phys. Rev. X 8, 031061 (2018)], it has been experimentally found that by shining a laser light on the fluorophores attached to a protein the energy fed to it can be channeled into the normal mode of lowest frequency of vibration thus making the subunits of the protein coherently oscillate. Even if the phonon condensation phenomenon has been theoretically explained, the first step - the energy transfer from electronic excitation into phonon excitation - has been left open. The present work is aimed at filling this gap.


title: Fair clustering via equitable group representations

abstract: What does it mean for a clustering to be fair? One popular approach seeks to ensure that each cluster contains groups in (roughly) the same proportion in which they exist in the population. The normative principle at play is balance: any cluster might act as a representative of the data, and thus should reflect its diversity. But clustering also captures a different form of representativeness. A core principle in most clustering problems is that a cluster center should be representative of the cluster it represents, by being "close" to the points associated with it. This is so that we can effectively replace the points by their cluster centers without significant loss in fidelity, and indeed is a common "use case" for clustering. For such a clustering to be fair, the centers should "represent" different groups equally well. We call such a clustering a group-representative clustering. In this paper, we study the structure and computation of group-representative clusterings. We show that this notion naturally parallels the development of fairness notions in classification, with direct analogs of ideas like demographic parity and equal opportunity. We demonstrate how these notions are distinct from and cannot be captured by balance-based notions of fairness. We present approximation algorithms for group representative $k$-median clustering and couple this with an empirical evaluation on various real-world data sets.


title: Pulse-timing symmetry breaking in an excitable optical system with delay

abstract: Excitable systems with delayed feedback are important in areas from biology to neuroscience and optics. They sustain multistable pulsing regimes with different number of equidistant pulses in the feedback loop. Experimentally and theoretically, we report on the pulse-timing symmetry breaking of these regimes in an optical system. A bifurcation analysis unveils that this originates in a resonance phenomenon and that symmetry-broken states are stable in large regions of the parameter space. These results have impact in photonics for e.g. optical computing and versatile sources of optical pulses.


title: Disentangling User Interest and Conformity for Recommendation with Causal Embedding

abstract: Recommendation models are usually trained on observational interaction data. However, observational interaction data could result from users' conformity towards popular items, which entangles users' real interest. Existing methods tracks this problem as eliminating popularity bias, e.g., by re-weighting training samples or leveraging a small fraction of unbiased data. However, the variety of user conformity is ignored by these approaches, and different causes of an interaction are bundled together as unified representations, hence robustness and interpretability are not guaranteed when underlying causes are changing. In this paper, we present DICE, a general framework that learns representations where interest and conformity are structurally disentangled, and various backbone recommendation models could be smoothly integrated. We assign users and items with separate embeddings for interest and conformity, and make each embedding capture only one cause by training with cause-specific data which is obtained according to the colliding effect of causal inference. Our proposed methodology outperforms state-of-the-art baselines with remarkable improvements on two real-world datasets on top of various backbone models. We further demonstrate that the learned embeddings successfully capture the desired causes, and show that DICE guarantees the robustness and interpretability of recommendation.


title: Great Chiral Fluorescence from Optical Duality Silver Nanostructures Enabled by 3D Laser Printing

abstract: Featured by prominent flexibility and fidelity in producing sophisticated stereoscopic structures transdimensionally, three-dimensional (3D) laser printing technique has vastly extended the toolkit for delivering diverse functional devices. Yet chiral nanoemitters heavily resorting to artificial structures that manifest efficient emission and tightly confined light-mater interactions simultaneously remains alluring but dauntingly challenging for this technique at this moment. In this work, we assert the chiral photoluminescence is implemented from silver nanostructures of optical duality in one go via a twofold three-dimensional laser printing scheme. Such laser printing protocol allows the highly desired duality by simultaneously producing uniformly distributed fluorescent silver nanoclusters and aggregated plasmonic silver nanoparticles to tightly confine chiral interactions at the nanoscale. A helical emitter of 550 nm-helix-diameter as fabricated has seen a record-high luminescence anisotropic factor with the absolute value up to 0.58, which is two orders of magnitude greater than fluorescent chiral silver clusters. This method holds great promise for future versatile applications in chiroptical nanodevices.


title: Correlators in the supereigenvalue model in the Ramond sector

abstract: We investigate the supereigenvalue model in the Ramond sector. We prove that its partition function can be obtained by acting on elementary functions with exponents of the given operators. The Virasoro constraints for this supereigenvalue model are presented. The remarkable property of these bosonic constraint operators is that they obey the Witt algebra and null 3-algebra. The compact expression of correlators can be derived from these Virasoro constraints.


title: Gradient boosting machine with partially randomized decision trees

abstract: The gradient boosting machine is a powerful ensemble-based machine learning method for solving regression problems. However, one of the difficulties of its using is a possible discontinuity of the regression function, which arises when regions of training data are not densely covered by training points. In order to overcome this difficulty and to reduce the computational complexity of the gradient boosting machine, we propose to apply the partially randomized trees which can be regarded as a special case of the extremely randomized trees applied to the gradient boosting. The gradient boosting machine with the partially randomized trees is illustrated by means of many numerical examples using synthetic and real data.


title: Chaotic fluctuations in graphs with amplification

abstract: We consider a model for chaotic diffusion with amplification on graphs associated with piecewise-linear maps of the interval. We investigate the possibility of having power-law tails in the invariant measure by approximate solution of the Perron-Frobenius equation and discuss the connection with the generalized Lyapunov exponents $L(q)$. We then consider the case of open maps where trajectories escape and demonstrate that stationary power-law distributions occur when $L(q)=r$, with $r$ being the escape rate. The proposed system is a toy model for coupled active chaotic cavities or lasing networks and allows to elucidate in a simple mathematical framework the conditions for observing L\'evy statistical regimes and chaotic intermittency in such systems.


title: Imperfections and Corrections

abstract: After a review of linear imperfections and their causes, we discuss how to model them, the diagnostic equipment needed to monitor them, and the correction algorithms to fix the problem they cause. We first address linear systems - beam lines or linear accelerators. In a later part we cover circular systems, such as storage rings.


title: Electroluminescence and electron avalanching in two-phase detectors

abstract: Electroluminescence and electron avalanching are the physical effects used in two-phase argon and xenon detectors for dark mater search and neutrino detection, to amplify the primary ionization signal directly in cryogenic noble-gas media. We review the concepts of such light and charge signal amplification, including a combination thereof, both in the gas and in the liquid phase. Puzzling aspects of the physics of electroluminescence and electron avalanching in two-phase detectors are explained and detection techniques based on these effects are described.


title: Solenoidal extensions in domains with obstacles: explicit bounds and applications to Navier-Stokes equations

abstract: We introduce a new method for constructing solenoidal extensions of fairly general boundary data in (2d or 3d) cubes that contain an obstacle. This method allows us to provide explicit bounds for the Dirichlet norm of the extensions. It runs as follows: by inverting the trace operator, we first determine suitable extensions, not necessarily solenoidal, of the data; then we analyze the Bogovskii problem with the resulting divergence to obtain a solenoidal extension; finally, by solving a variational problem involving the infinity-Laplacian and using ad hoc cutoff functions, we find explicit bounds in terms of the geometric parameters of the obstacle. The natural applications of our results lie in the analysis of inflow-outflow problems, in which an explicit bound on the inflow velocity is needed to estimate the threshold for uniqueness in the stationary Navier-Stokes equations and, in case of symmetry, the stability of the obstacle immersed in the fluid flow.


title: Relativistic Corrections to the Growth of Structure in Modified Gravity

abstract: We present a method to introduce relativistic corrections including linear dark energy perturbations in Horndeski theory into Newtonian simulations based on the N-body gauge approach. We assume that standard matter species (cold dark matter, baryons, photons and neutrinos) are only gravitationally-coupled with the scalar field and we then use the fact that one can include modified gravity effects as an effective dark energy fluid in the total energy-momentum tensor. In order to compute the scalar field perturbations, as well as the cosmological background and metric perturbations, we use the Einstein-Boltzmann code \hiclass. As an example, we study the impact of relativistic corrections on the matter power spectrum in k-essence, a subclass of Horndeski theory, including the effects of massless and massive neutrinos. For massive neutrinos with $\sum m_{\nu} = 0.1$ eV, the corrections due to relativistic species (photons, neutrinos and dark energy) can introduce a maximum deviation of approximately $7\%$ to the power spectrum at $k \sim 10^{-3} \ \textrm{Mpc}^{-1}$ at $z=0$, for a scalar field with sound speed $c_{s}^{2}\sim 0.013$ during matter domination epoch. Our formalism makes it possible to test beyond $\Lambda$CDM models probed by upcoming large-scale structure surveys on very large scales.


title: Representing Pure Nash Equilibria in Argumentation

abstract: In this paper we describe an argumentation-based representation of normal form games, and demonstrate how argumentation can be used to compute pure strategy Nash equilibria. Our approach builds on Modgil's Extended Argumentation Frameworks. We demonstrate its correctness, prove several theoretical properties it satisfies, and outline how it can be used to explain why certain strategies are Nash equilibria to a non-expert human user.


title: Boosting Active Learning for Speech Recognition with Noisy Pseudo-labeled Samples

abstract: The cost of annotating transcriptions for large speech corpora becomes a bottleneck to maximally enjoy the potential capacity of deep neural network-based automatic speech recognition models. In this paper, we present a new training pipeline boosting the conventional active learning approach targeting label-efficient learning to resolve the mentioned problem. Existing active learning methods only focus on selecting a set of informative samples under a labeling budget. One step further, we suggest that the training efficiency can be further improved by utilizing the unlabeled samples, exceeding the labeling budget, by introducing sophisticatedly configured unsupervised loss complementing supervised loss effectively. We propose new unsupervised loss based on consistency regularization, and we configure appropriate augmentation techniques for utterances to adopt consistency regularization in the automatic speech recognition task. From the qualitative and quantitative experiments on the real-world dataset and under real-usage scenarios, we show that the proposed training pipeline can boost the efficacy of active learning approaches, thus successfully reducing a sustainable amount of human labeling cost.


title: Learning Stabilizing Controllers for Unstable Linear Quadratic Regulators from a Single Trajectory

abstract: The principal task to control dynamical systems is to ensure their stability. When the system is unknown, robust approaches are promising since they aim to stabilize a large set of plausible systems simultaneously. We study linear controllers under quadratic costs model also known as linear quadratic regulators (LQR). We present two different semi-definite programs (SDP) which results in a controller that stabilizes all systems within an ellipsoid uncertainty set. We further show that the feasibility conditions of the proposed SDPs are \emph{equivalent}. Using the derived robust controller syntheses, we propose an efficient data dependent algorithm -- \textsc{eXploration} -- that with high probability quickly identifies a stabilizing controller. Our approach can be used to initialize existing algorithms that require a stabilizing controller as an input while adding constant to the regret. We further propose different heuristics which empirically reduce the number of steps taken by \textsc{eXploration} and reduce the suffered cost while searching for a stabilizing controller.


title: On the X-ray scattering pre-peak of linear mono-ols and the related micro-structure from computer simulations

abstract: The X-ray scattering intensities I(k) of linear alkanols OH(CH2)n-1CH3, obtained from experiments (methanol to 1-undecanol) and computer simulations (methanol to 1-nonanol) of different force field models, are comparatively studied, particularly in order to explain the origin and the properties of the scattering pre-peak in the k-vector range 0.3A^{-1}-1A^{-1}. The experimental I(k) show two apparent features: the pre-peak position kP decreases with increasing n, and more intriguingly, the amplitude AP goes through a maximum at 1-butanol (n=4). The first feature is well reproduced by all force field models, while the second shows a strong model dependence. The simulations reveal various shapes of clusters of the hydroxyl head-group, from n>2. kP is directly related to the size of the \emph{meta-objects} corresponding to such clusters surrounded by their alkyl tails. The explanation of the Ap turnover at n=4 is more involved, in terms of cancellations of atom-atom S(k) contributions related to domain ordering. The flexibility of the alkyl tails tend to reduce the cross contributions, thus revealing the crucial importance of this parameter in the models. Force fields with all-atom representation are less successful in reproducing the pre-peak features for smaller alkanols n<6, possibly because they blur the charge ordering process since all atoms bear partial charges. The analysis clearly shows that it is not possible to obtain a model free explanation of the features of I(k)


title: Separation of Geometrical and Topological Entanglement in Confined Polymers Driven Out-Of-Equilibrium

abstract: We use Brownian dynamics simulations and advanced topological profiling methods to characterize the out-of-equilibrium evolution of self-entanglement in linear polymers confined into nano-channels and under periodic compression. We distinguish two main forms of entanglement, geometrical and topological. The latter is measured by the number of (essential) crossings of the physical knot detected after a suitable bridging of the chain termini. The former is instead measured as the average number of times a linear chain appears to cross itself when viewed under all projections, and is irrespective of the physical knotted state. The key discovery of our work is that these two forms of entanglement are uncoupled and evolve with distinct dynamics. While geometrical entanglement is typically in phase with the compression-elongation cycles and it is primarily sensitive to its force f, the topological measure is mildly sensitive to cyclic modulation but strongly depends on both compression force f and duration k. The findings could assist the interpretation of experiments using fluorescence molecular tracers to track physical knots in polymers. Furthermore, we identify optimal regions in the experimentally-controllable parameter space in which to obtain more/less topological and geometrical entanglement; this may help designing polymers with targeted topology.


title: Design of a Near-Ideal Fault-Tolerant Routing Algorithm for Network-on-Chip-Based Multicores

abstract: With relentless CMOS technology downsizing Networks-on-Chips (NoCs) are inescapably experiencing escalating susceptibility to wearout and reduced reliability. While faults in processors and memories may be masked via redundancy, or mitigated via techniques such as task migration, NoCs are especially vulnerable to hardware faults as a single link breakdown may cause inter-tile communication to halt indefinitely, rendering the whole multicore chip inoperable. As such, NoCs impose the risk of becoming the pivotal point of failure in chip multicores that utilize them. Aiming towards seamless NoC operation in the presence of faulty links we propose Hermes, a near-ideal fault-tolerant routing algorithm that meets the objectives of exhibiting high levels of robustness, operating in a distributed mode, guaranteeing freedom from deadlocks, and evening-out traffic, among many. Hermes is a limited-overhead deadlock-free hybrid routing algorithm, utilizing load-balancing routing on fault-free paths to sustain high-throughput, while providing pre-reconfigured escape path selection in the vicinity of faults. Under such online mechanisms, Hermes's performance degrades gracefully with increasing faulty link counts, a crucially desirable response lacking in prior-art. Additionally, Hermes identifies non-communicating network partitions in scenarios where faulty links are topologically densely distributed such that packets being routed to physically isolated regions cause no network stagnation due to indefinite chained blockages starting at sub-network boundaries. An extensive experimental evaluation, including utilizing traffic workloads gathered from full-system chip multi-processor simulations, shows that Hermes improves network throughput by up to $3\times$ when compared against the state-of-the-art. Further, hardware synthesis results prove Hermes's efficacy.


title: Hybridizing the 1/5-th Success Rule with Q-Learning for Controlling the Mutation Rate of an Evolutionary Algorithm

abstract: It is well known that evolutionary algorithms (EAs) achieve peak performance only when their parameters are suitably tuned to the given problem. Even more, it is known that the best parameter values can change during the optimization process. Parameter control mechanisms are techniques developed to identify and to track these values. Recently, a series of rigorous theoretical works confirmed the superiority of several parameter control techniques over EAs with best possible static parameters. Among these results are examples for controlling the mutation rate of the $(1+\lambda)$~EA when optimizing the OneMax problem. However, it was shown in [Rodionova et al., GECCO'19] that the quality of these techniques strongly depends on the offspring population size $\lambda$. We introduce in this work a new hybrid parameter control technique, which combines the well-known one-fifth success rule with Q-learning. We demonstrate that our HQL mechanism achieves equal or superior performance to all techniques tested in [Rodionova et al., GECCO'19] and this -- in contrast to previous parameter control methods -- simultaneously for all offspring population sizes $\lambda$. We also show that the promising performance of HQL is not restricted to OneMax, but extends to several other benchmark problems.


title: Optimal Bounds in Normal Approximation for Many Interacting Worlds

abstract: In this paper, we use Stein's method to obtain optimal bounds, both in Kolmogorov and in Wasserstein distance, in the normal approximation for the empirical distribution of the ground state of a many-interacting-worlds harmonic oscillator proposed by Hall, Deckert, and Wiseman [Phys. Rev. X. (2014)]. Our bounds on the Wasserstein distance solve a conjecture of McKeague and Levin [Ann. Appl. Probab. (2016)].


title: On additive functions with additional derivation properties

abstract: The purpose of this paper is to introduce the notion of a generalized derivation which derivates a prescribed family of smooth vector-valued functions of several variables. The basic calculus rules are established and then a result derived which shows that if a function $f$ satisfies an addition theorem whose determining operation is derivable with respect to an additive function $d$, then the function $f$ is itself derivable with respect to $d$. As an application of this approach, new proof of a generalization of a recent result of Maksa is obtained. We also extend the result of Nishiyama and Horinouchi and formulate two open problems.


title: Learning Optimal Power Flow: Worst-Case Guarantees for Neural Networks

abstract: This paper introduces for the first time a framework to obtain provable worst-case guarantees for neural network performance, using learning for optimal power flow (OPF) problems as a guiding example. Neural networks have the potential to substantially reduce the computing time of OPF solutions. However, the lack of guarantees for their worst-case performance remains a major barrier for their adoption in practice. This work aims to remove this barrier. We formulate mixed-integer linear programs to obtain worst-case guarantees for neural network predictions related to (i) maximum constraint violations, (ii) maximum distances between predicted and optimal decision variables, and (iii) maximum sub-optimality. We demonstrate our methods on a range of PGLib-OPF networks up to 300 buses. We show that the worst-case guarantees can be up to one order of magnitude larger than the empirical lower bounds calculated with conventional methods. More importantly, we show that the worst-case predictions appear at the boundaries of the training input domain, and we demonstrate how we can systematically reduce the worst-case guarantees by training on a larger input domain than the domain they are evaluated on.


title: X-ray detection of ultrashort spin current pulses in synthetic antiferromagnets

abstract: We explore the ultrafast generation of spin currents in magnetic multilayer samples by applying fs laser pulses to one layer and measuring the magnetic response in the other layer by element-resolved x-ray spectroscopy. In Ni(5~nm)/Ru(2~nm)/Fe(4~nm), the Ni and Fe magnetization directions couple antiferromagnetically due to the RKKY interaction, but may be oriented parallel through an applied magnetic field. After exciting the top Ni layer with a fs laser pulse, we find that also the Fe layer underneath demagnetizes, with a $4.1 \pm 1.9$\% amplitude difference between parallel and antiparallel orientation of the Ni and Fe magnetizations. We attribute this difference to the influence of a spin current generated by the fs laser pulse that transfers angular momentum from the Ni into the Fe layer. Our results confirm that superdiffusive spin transport plays a role in determining the sub-ps demagnetization dynamics of synthetic antiferromagnetic layers, but also evidence large depolarization effects due to hot electron dynamics, which are independent of the relative alignment of the magnetization in Ni and Fe.


title: Semi-supervised time series classification method for quantum computing

abstract: In this paper we develop methods to solve two problems related to time series (TS) analysis using quantum computing: reconstruction and classification. We formulate the task of reconstructing a given TS from a training set of data as an unconstrained binary optimization (QUBO) problem, which can be solved by both quantum annealers and gate-model quantum processors. We accomplish this by discretizing the TS and converting the reconstruction to a set cover problem, allowing us to perform a one-versus-all method of reconstruction. Using the solution to the reconstruction problem, we show how to extend this method to perform semi-supervised classification of TS data. We present results indicating our method is competitive with current semi- and unsupervised classification techniques, but using less data than classical techniques.


title: Confinement-induced demixing and crystallization

abstract: We simulate a strongly size-disperse hard-sphere fluid confined between two parallel, hard walls. We find that confinement induces crystallization into n-layered hexagonal lattices and a novel honeycomb-shaped structure, facilitated by fractionation. The onset of freezing prevents the formation of a stable glass phase and occurs at much smaller packing fraction than in bulk. Varying the wall separation triggers solid-to-solid transitions and a systematic change of the size-distribution of crystalline particles, which we rationalize using a semi-quantitative theory. We show that the crystallization can be exploited in a wedge geometry to demix particles of different sizes.


title: Towards Reliable Real-time Opera Tracking: Combining Alignment with Audio Event Detectors to Increase Robustness

abstract: Recent advances in real-time music score following have made it possible for machines to automatically track highly complex polyphonic music, including full orchestra performances. In this paper, we attempt to take this to an even higher level, namely, live tracking of full operas. We first apply a state-of-the-art audio alignment method based on online Dynamic Time-Warping (OLTW) to full-length recordings of a Mozart opera and, analyzing the tracker's most severe errors, identify three common sources of problems specific to the opera scenario. To address these, we propose a combination of a DTW-based music tracker with specialized audio event detectors (for applause, silence/noise, and speech) that condition the DTW algorithm in a top-down fashion, and show, step by step, how these detectors add robustness to the score follower. However, there remain a number of open problems which we identify as targets for ongoing and future research.


title: Low-cost Retina-like Robotic Lidars Based on Incommensurable Scanning

abstract: High performance lidars are essential in autonomous robots such as self-driving cars, automated ground vehicles and intelligent machines. Traditional mechanical scanning lidars offer superior performance in autonomous vehicles, but the potential mass application is limited by the inherent manufacturing difficulty. We propose a robotic lidar sensor based on incommensurable scanning that allows straightforward mass production and adoption in autonomous robots. Some unique features are additionally permitted by this incommensurable scanning. Similar to the fovea in human retina, this lidar features a peaked central angular density, enabling in applications that prefers eye-like attention. The incommensurable scanning method of this lidar could also provide a much higher resolution than conventional lidars which is beneficial in robotic applications such as sensor calibration. Examples making use of these advantageous features are demonstrated.


title: Wave Propagation of Visual Stimuli in Focus of Attention

abstract: Fast reactions to changes in the surrounding visual environment require efficient attention mechanisms to reallocate computational resources to most relevant locations in the visual field. While current computational models keep improving their predictive ability thanks to the increasing availability of data, they still struggle approximating the effectiveness and efficiency exhibited by foveated animals. In this paper, we present a biologically-plausible computational model of focus of attention that exhibits spatiotemporal locality and that is very well-suited for parallel and distributed implementations. Attention emerges as a wave propagation process originated by visual stimuli corresponding to details and motion information. The resulting field obeys the principle of "inhibition of return" so as not to get stuck in potential holes. An accurate experimentation of the model shows that it achieves top level performance in scanpath prediction tasks. This can easily be understood at the light of a theoretical result that we establish in the paper, where we prove that as the velocity of wave propagation goes to infinity, the proposed model reduces to recently proposed state of the art gravitational models of focus of attention.


title: The interplay between randomness and structure during learning in RNNs

abstract: Recurrent neural networks (RNNs) trained on low-dimensional tasks have been widely used to model functional biological networks. However, the solutions found by learning and the effect of initial connectivity are not well understood. Here, we examine RNNs trained using gradient descent on different tasks inspired by the neuroscience literature. We find that the changes in recurrent connectivity can be described by low-rank matrices, despite the unconstrained nature of the learning algorithm. To identify the origin of the low-rank structure, we turn to an analytically tractable setting: training a linear RNN on a simplified task. We show how the low-dimensional task structure leads to low-rank changes to connectivity. This low-rank structure allows us to explain and quantify the phenomenon of accelerated learning in the presence of random initial connectivity. Altogether, our study opens a new perspective to understanding trained RNNs in terms of both the learning process and the resulting network structure.


title: Fast Scrambling under an RG-flow

abstract: The notion of information scrambling is tied to the long-time behaviour of a system and therefore is related to its infra-red dynamics. In fast scramblers, information spreads as a logarithmic function of the number of degrees of freedom. Ordinarily, under an RG-flow, scrambling will become faster in the IR since many degrees of freedom are integrated out, as a consequence of the c-theorem. In this article, we consider a class of Holographic quantum field theories (QFT), which are strongly coupled large $N$ gauge theories with large number of adjoint and fundamental matter, in which scrambling slows down in the IR. This happens since more degrees of freedom are inserted in the IR, compared to the UV, in a precise sense. For generic large $N$ gauge theories, we also explore general, perturbative flow features of the corresponding Lyapunov exponent, based on the Callan-Symanzik equation.


title: Numerical solution for Fokker-Planck equation using a two-level scheme

abstract: A numerical solution to the Fokker-Planck equation using a two-level scheme is presented. The Fokker-Planck (FP) equation is of parabolic type equation govern the time evolution of probability density function of the stochastic processes. The FP equation also preserves the positivity and conservative of the total probability. A Chang-Cooper discretization scheme is used to ensure the positiveness and conservation of the total probability with second-order accuracy. We investigate a two-level scheme with factor-three-coarsening strategy and have a significant reduction in computations and CPU time. Numerical experiments are performed to validate the efficiency and second-order accuracy of the proposed two-level algorithm with backward time-difference schemes.


title: Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds

abstract: Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas. However, exact local energy-momentum conservation law for the electromagnetic gyrokinetic system has not been found despite continuous effort. Without such a local conservation law, energy-momentum can be instantaneously transported across spacetime, which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory. Standard Noether's procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds. To overcome this difficulty, we developed a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact local conservation laws, in particular the energy-momentum conservation law, for the electromagnetic gyrokinetic system. A weak Euler-Lagrange equation is established to replace the standard Euler-Lagrange equation for the particles. It is discovered that an induced weak Euler-Lagrange current enters the local conservation laws. And it is the new physics captured by the high-order field theory on heterogeneous manifolds.


title: Coupled-mode theory for microresonators with quadratic nonlinearity

abstract: We use Maxwell's equations to derive several models describing the interaction of the multi-mode fundamental field and its second harmonic in a ring microresonator with quadratic nonlinearity and quasi-phase-matching. We demonstrate how multi-mode three-wave mixing sums entering nonlinear polarisation response can be calculated via Fourier transforms of products of the field envelopes. Quasi-phase-matching gratings with arbitrary profiles are incorporated seamlessly into our models. We also introduce several levels of approximations allowing to account for dispersion of nonlinear coefficients and demonstrate how coupled-mode equations can be reduced to the envelope Lugiato-Lefever-like equations with self-steepening terms. An estimate for the $\chi^{(2)}$ induced cascaded Kerr nonlinearity, in the regime of imperfect phase-matching, puts it above the intrinsic Kerr effect by several orders of magnitude.


title: Bayesian analysis of mixture autoregressive models covering the complete parameter space

abstract: Mixture autoregressive (MAR) models provide a flexible way to model time series with predictive distributions which depend on the recent history of the process and are able to accommodate asymmetry and multimodality. Bayesian inference for such models offers the additional advantage of incorporating the uncertainty in the estimated models into the predictions. We introduce a new way of sampling from the posterior distribution of the parameters of MAR models which allows for covering the complete parameter space of the models, unlike previous approaches. We also propose a relabelling algorithm to deal a posteriori with label switching. We apply our new method to simulated and real datasets, discuss the accuracy and performance of our new method, as well as its advantages over previous studies. The idea of density forecasting using MCMC output is also introduced.


title: Robust and scalable h-adaptive aggregated unfitted finite elements for interface elliptic problems

abstract: This work introduces a novel, fully robust and highly-scalable, $h$-adaptive aggregated unfitted finite element method for large-scale interface elliptic problems. The new method is based on a recent distributed-memory implementation of the aggregated finite element method atop a highly-scalable Cartesian forest-of-trees mesh engine. It follows the classical approach of weakly coupling nonmatching discretisations at the interface to model internal discontinuities at the interface. We propose a natural extension of a single-domain parallel cell aggregation scheme to problems with a finite number of interfaces; it straightforwardly leads to aggregated finite element spaces that have the structure of a Cartesian product. We demonstrate, through standard numerical analysis and exhaustive numerical experimentation on several complex Poisson and linear elasticity benchmarks, that the new technique enjoys the following properties: well-posedness, robustness with respect to cut location and material contrast, optimal ($h$-adaptive) approximation properties, high scalability and easy implementation in large-scale finite element codes. As a result, the method offers great potential as a useful finite element solver for large-scale interface problems modelled by partial differential equations.


title: On the Principle of Accountability: Challenges for Smart Homes & Cybersecurity

abstract: This chapter introduces the Accountability Principle and its role in data protection governance. We focus on what accountability means in the context of cybersecurity management in smart homes, considering the EU General Data Protection Law requirements to secure personal data. This discussion sits against the backdrop of two key new developments in data protection law. Firstly, the law is moving into the home, due to narrowing of the so called household exemption. Concurrently, household occupants may now have legal responsibilities to comply with the GDPR, as they find themselves jointly responsible for compliance, as they are possibly held to determine the means and purposes of data collection with IoT device vendors. As a complex socio-technical space, we consider the interactions between accountability requirements and the competencies of this new class of domestic data controllers (DDCs). Specifically, we consider the value and limitations of edge-based security analytics to manage smart home cybersecurity risks, reviewing a range of prototypes and studies of their use. We also reflect on interpersonal power dynamics in the domestic setting e.g. device control; existing social practices around privacy and security management in smart homes; and usability issues that may hamper DDCs ability to rely on such solutions. We conclude by reflecting on 1) the need for collective security management in homes and 2) the increasingly complex divisions of responsibility in smart homes between device users, account holders, IoT device/software/firmware vendors, and third parties.


title: V-Dream: Immersive Exploration of Generative Design Solution Space

abstract: Generative Design workflows have introduced alternative paradigms in the domain of computational design, allowing designers to generate large pools of valid solutions by defining a set of goals and constraints. However, analyzing and narrowing down the generated solution space, which usually consists of various high-dimensional properties, has been a major challenge in current generative workflows. By taking advantage of the interactive unbounded spatial exploration, and the visual immersion offered in virtual reality platforms, we propose V-Dream, a virtual reality generative analysis framework for exploring large-scale solution spaces. V-Dream proposes a hybrid search workflow in which a spatial stochastic search approach is combined with a recommender system allowing users to pick desired candidates and eliminate the undesired ones iteratively. In each cycle, V-Dream reorganizes the remaining options in clusters based on the defined features. Moreover, our framework allows users to inspect design solutions and evaluate their performance metrics in various hierarchical levels, assisting them in narrowing down the solution space through iterative cycles of search/select/re-clustering of the solutions in an immersive fashion. Finally, we present a prototype of our proposed framework, illustrating how users can navigate and narrow down desired solutions from a pool of over 16000 monitor stands generated by Autodesk's Dreamcatcher software.


title: Riesz means on locally symmetric spaces

abstract: We prove that for a certain class of $n$ dimensional rank one locally symmetric spaces, if $f \in L^p$, $1\leq p \leq 2$, then the Riesz means of order $z$ of $f$ converge to $f$ almost everywhere, for $\operatorname{Re}z> (n-1)(1/p-1/2).$


title: The precise representative for the gradient of the Riesz potential of a finite measure

abstract: Given a finite nonnegative Borel measure $m$ in $\mathbb{R}^{d}$, we identify the Lebesgue set $\mathcal{L}(V_{s}) \subset \mathbb{R}^{d}$ of the vector-valued function $$V_{s}(x) = \int_{\mathbb{R}^{d}}\frac{x - y}{|x - y|^{s + 1}} \mathrm{d}m(y), $$ for any order $0 < s < d$. We prove that $a \in \mathcal{L}(V_{s})$ if and only if the integral above has a principal value at $a$ and $$\lim_{r \to 0}{\frac{m(B_{r}(a))}{r^{s}}} = 0.$$ In that case, the precise representative of $V_{s}$ at $a$ coincides with the principal value of the integral. We also study the existence of Lebesgue points for the Cauchy integral of the intrinsic probability measure associated with planar Cantor sets, which leads to challenging new questions.


title: Gapped dilatons in scale invariant superfluids

abstract: We study a paradigmatic model in field theory where a global $U(1)$ and scale symmetries are jointly and spontaneously broken. At zero density the model has a non-compact flat direction, which at finite density needs to be slightly lifted. The resulting low-energy spectrum is composed by a standard gapless $U(1)$ Nambu-Goldstone mode and a light dilaton whose gap is determined by the chemical potential and corrected by the couplings. Even though $U(1)$ and scale symmetries commute, there is a mixing between the $U(1)$ Nambu-Goldstone and the dilaton that is crucial to recover the expected dynamics of a conformal fluid and leads to a phonon propagating at the speed of sound. The results rely solely on an accurate study of the Ward-Takahashi identities and are checked against standard fluctuation computations. We extend our results to a boosted superfluid, and comment the relevance of our findings to condensed matter applications.

categories: hep-th cond-mat.quant-gas cond-mat.str-el hep-ph


title: Electronic and magnetic properties of graphene quantum dots with two charged vacancies

abstract: Electronic and magnetic properties of a system of two charged vacancies in hexagonal shaped graphene quantum dots are investigated using a mean-field Hubbard model as a function of the Coulomb potential strength $\beta$ of the charge impurities and the distance R between them. For $\beta=0$, the magnetic properties of the vacancies are dictated by Lieb's rules where the opposite (same) sub-lattice vacancies are coupled antiferromagnetically (ferromagnetically) and exhibit Fermi oscillations. Here, we demonstrate the emergence of a non-magnetic regime within the subcritical region: as the Coulomb potential strength is increased to $\beta \sim 0.1$ , before reaching the frustrated atomic collapse regime, the magnetization is strongly suppressed and the ground state total spin is given by $S_{z}=0$ both for opposite and same sublattice vacancy configurations. When long-range electron-electronz in-teractions are included within extended mean-field Hubbard model, the critical value for the frustrated collapse increases from $\beta_{cf} \sim 0.28$ to $\beta_{cf} \sim 0.36$ for $R < 27 \ \r{A} $ .


title: On the diversity of the jet production efficiency in Swift/BAT AGNs

abstract: We address the very large diversity of the jet production efficiency in active galactic nuclei (AGNs) by using data on low redshift AGNs selected from the Swift/BAT catalog and having black hole (BH) masses larger than $10^{8.5}\,M_{\odot}$. Most of these AGNs accrete at intermediate rates and have bolometric luminosities dominated by mid-IR radiation. Our sample contains $14\%$ radio-loud (RL), $6\%$ radio-intermediate (RI), and $80\%$ radio-quiet (RQ) AGNs. All RL objects are found to have extended radio structures and most of them have classical FR II morphology. Converting their radio loudness to the jet production efficiency, we find that the median of this efficiency is on the order of $(\epsilon_d/0.1)\%$, where $\epsilon_d=L_{\rm bol}/\dot{M}c^2$ is the radiation efficiency of the accretion disk. Without knowing the contribution of jets to the radio emission in the RQ AGNs, we are only able to estimate their efficiencies using upper limits. Their median is found to be $0.002(\epsilon_d/0.1)\%$. Our results suggest that some threshold conditions must be satisfied to allow production of strong, relativistic jets in RL AGNs. We discuss several possible scenarios and argue that the production of collimated, relativistic jets must involve the Blandford-Znajek mechanism and can be activated only in those AGNs whose lifetime is longer than the time required to enter the magnetically arrested disk (MAD). Presuming that MAD is required to collimate relativistic jets, we expect that the weak nonrelativistic jets observed in some RQ AGNs are produced by accretion disks rather than by rotating BHs.


title: The Brownian disk viewed from a boundary point

abstract: We provide a new construction of Brownian disks in terms of forests of continuous random trees equipped with nonnegative labels corresponding to distances from a distinguished point uniformly distributed on the boundary of the disk. This construction shows in particular that distances from the distinguished point evolve along the boundary as a five-dimensional Bessel bridge. As an important ingredient of our proofs, we show that the uniform measure on the boundary, as defined in the earlier work of Bettinelli and Miermont, is the limit of the suitably normalized volume measure on a small tubular neighborhood of the boundary. Our construction also yields a simple proof of the equivalence between the two definitions of the Brownian half-plane.


title: Earth and Mars -- distinct inner Solar System products

abstract: Composition of terrestrial planets records planetary accretion, core-mantle and crust-mantle differentiation, and surface processes. Here we compare the compositional models of Earth and Mars to reveal their characteristics and formation processes. Earth and Mars are equally enriched in refractory elements (1.9 $\times$ CI), although Earth is more volatile-depleted and less oxidized than Mars. Their chemical compositions were established by nebular fractionation, with negligible contributions from post-accretionary losses of moderately volatile elements. The degree of planetary volatile element depletion might correlate with the abundances of chondrules in the accreted materials, planetary size, and their accretion timescale, which provides insights into composition and origin of Mercury, Venus, the Moon-forming giant impactor, and the proto-Earth. During its formation before and after the nebular disk's lifetime, the Earth likely accreted more chondrules and less matrix-like materials than Mars and chondritic asteroids, establishing its marked volatile depletion. A giant impact of an oxidized, differentiated Mars-like (i.e., composition and mass) body into a volatile-depleted, reduced proto-Earth produced a Moon-forming debris ring with mostly a proto-Earth's mantle composition. Chalcophile and some siderophile elements in the silicate Earth added by the Mars-like impactor were extracted into the core by a sulfide melt. In contrast, the composition of Mars indicates its rapid accretion of lesser amounts of chondrules under nearly uniform oxidizing conditions. Mars' rapid cooling and early loss of its dynamo likely led to the absence of plate tectonics and surface water, and the present-day low surface heat flux. These similarities and differences between the Earth and Mars made the former habitable and the other inhospitable to uninhabitable.


title: Understanding Recurrent Neural Networks Using Nonequilibrium Response Theory

abstract: Recurrent neural networks (RNNs) are brain-inspired models widely used in machine learning for analyzing sequential data. The present work is a contribution towards a deeper understanding of how RNNs process input signals using the response theory from nonequilibrium statistical mechanics. For a class of continuous-time stochastic RNNs (SRNNs) driven by an input signal, we derive a Volterra type series representation for their output. This representation is interpretable and disentangles the input signal from the SRNN architecture. The kernels of the series are certain recursively defined correlation functions with respect to the unperturbed dynamics that completely determine the output. Exploiting connections of this representation and its implications to rough paths theory, we identify a universal feature -- the response feature, which turns out to be the signature of tensor product of the input signal and a natural support basis. In particular, we show that SRNNs, with only the weights in the readout layer optimized and the weights in the hidden layer kept fixed and not optimized, can be viewed as kernel machines operating on a reproducing kernel Hilbert space associated with the response feature.


title: Regular-Triangle Trimer and Charge Order Preserving the Anderson Condition in the Pyrochlore Structure of CsW$_2$O$_6$

abstract: Since the discovery of the Verwey transition in magnetite, transition metal compounds with pyrochlore structures have been intensively studied as a platform for realizing remarkable electronic phase transitions. We report the discovery of a unique phase transition that preserves the cubic symmetry of the beta-pyrochlore oxide CsW$_2$O$_6$, where each of W 5d electrons are confined in regular-triangle W3 trimers. This trimer formation is an unprecedented self-organization of d electrons, which can be resolved into a charge order satisfying the Anderson condition in a nontrivial way, orbital order caused by the distortion of WO6 octahedra, and the formation of a spin-singlet pair in a regular-triangle trimer. Electronic instability due to the unusual three-dimensional nesting of Fermi surfaces and the localized nature of the 5d electrons characteristic of the pyrochlore oxides were found to play important roles in this unique charge-orbital-spin coupled phenomenon.


title: Fully $S$-idempotent modules

abstract: Let R be a commutative ring with identity and S be a multiplicatively closed subset of R. The aim of this paper is to introduce the notion of fully S-idempotent modules as a generalization of fully idempotent modules and investigate some properties of this class of modules.


title: Synthesis and study of highly dense and smooth TiN thin films

abstract: This study aims towards a systematic reciprocity of the tunable synthesis parameters - partial pressure of N$_2$ gas, ion energy (\Ei) and Ti interface in TiN thin film samples deposited using ion beam sputtering at ambient temperature (300\,K). At the optimum partial pressure of N$_2$ gas, samples were prepared with or without Ti interface at \Ei~=~1.0 or 0.5\,keV. They were characterized using x-ray reflectivity (XRR) to deduce thickness, roughness and density. The roughness of TiN thin films was found to be below 1\,nm, when deposited at the lower \Ei~of 0.5\,keV and when interfaced with a layer of Ti. Under these conditions, the density of TiN sample reaches to 5.80($\pm$0.03)\,g~cm$^{-3}$, a value highest hitherto for any TiN sample. X-ray diffraction and electrical resistivity measurements were performed. It was found that the cumulative effect of the reduction in \Ei~from 1.0 to 0.5\,keV and the addition of Ti interface favors (111) oriented growth leading to dense and smooth TiN films and a substantial reduction in the electrical resistivity. The reduction in \Ei~has been attributed to the surface kinetics mechanism (simulated using SRIM) where the available energy of the sputtered species (\Esp) leaving the target at \Ei~= 0.5\,keV is the optimum value favoring the growth of defects free homogeneously distributed films. The electronic structure of samples was probed using N K-edge absorption spectroscopy and the information about the crystal field and spin-orbit splitting confirmed TiN phase formation. In essence, through this work, we demonstrate the role of \Esp~and Ti interface in achieving highly dense and smooth TiN thin films with low resistivity without the need of a high temperature or substrate biasing during the thin film deposition process.


title: A Survey of Syntactic-Semantic Parsing Based on Constituent and Dependency Structures

abstract: Syntactic and semantic parsing has been investigated for decades, which is one primary topic in the natural language processing community. This article aims for a brief survey on this topic. The parsing community includes many tasks, which are difficult to be covered fully. Here we focus on two of the most popular formalizations of parsing: constituent parsing and dependency parsing. Constituent parsing is majorly targeted to syntactic analysis, and dependency parsing can handle both syntactic and semantic analysis. This article briefly reviews the representative models of constituent parsing and dependency parsing, and also dependency graph parsing with rich semantics. Besides, we also review the closely-related topics such as cross-domain, cross-lingual and joint parsing models, parser application as well as corpus development of parsing in the article.


title: Symbolic dynamics in a binary asteroid system

abstract: We highlight the existence of a topological horseshoe arising from a a--priori stable model of the binary asteroid dynamics. The inspection is numerical and uses correctly aligned windows, as described in a recent paper by A. Gierzkiewicz and P. Zgliczy\'nski, combined with a recent analysis of an associated secular problem.


title: Empirical Relations for the Sensitivities of Solar-like Oscillations to Magnetic Perturbations

abstract: Oscillation mode frequencies of stars are typically treated as static for a given stellar model. However, in reality they can be perturbed by time varying sources such as magnetic fields and flows. We calculate the sensitivities of radial p-mode oscillations of a set of models for masses between 0.7-3.0 M$_{\odot}$ from the main sequence to the early asymptotic giant branch. We fit these mode sensitivities with polynomials in fundamental stellar parameters for six stages of stellar evolution. We find that the best-fitting relations differ from those proposed in the literature and change between stages of stellar evolution. Together with a measure of the strength of the perturbation, e.g., of the level of magnetic activity, the presented relations can be used for assessing whether a star's observed oscillation frequencies are likely to be close to the unperturbed ground state or whether they should be adjusted.


title: Room Temperature Magnetocaloric effect in CrTe1-xSbx Alloys

abstract: We investigate the magnetocaloric effect (MCE), relative cooling power (RCP) and crystalline structure in Sb substituted CrTe_{1-x}Sb_{x} (0 \leq x \leq 0.2) alloy. The Rietveld refinement of the XRD pattern of CrTe1-xSbx showed the emerging of pure hexagonal NiAs structure with P63/mmc (194) space group with increasing Sb substitution. We detect a slight increase in the basal plane a-lattice parameter, with a much larger reduction in the c-axis. Magnetic isotherms were measured in the temperature range of 50-400K. The results revealed an increase in the maximum entropy change |S_{M}(T,H)| with Sb-substitutions in the temperature range (~285-325K). Moreover, The RCP values increased by about 33% with 20% Sb substitutions. These findings suggest that CrTe_{1-x}Sb_{x} alloys can be used in room temperature magnetic cooling at fraction of the coast of pure Gd element the porotype magnetic material for magnetic refrigeration.


title: On the Time Trend of COVID-19: A Panel Data Study

abstract: In this paper, we study the trending behaviour of COVID-19 data at country level, and draw attention to some existing econometric tools which are potentially helpful to understand the trend better in future studies. In our empirical study, we find that European countries overall flatten the curves more effectively compared to the other regions, while Asia & Oceania also achieve some success, but the situations are not as optimistic elsewhere. Africa and America are still facing serious challenges in terms of managing the spread of the virus, and reducing the death rate, although in Africa the virus spreads slower and has a lower death rate than the other regions. By comparing the performances of different countries, our results incidentally agree with Gu et al. (2020), though different approaches and models are considered. For example, both works agree that countries such as USA, UK and Italy perform relatively poorly; on the other hand, Australia, China, Japan, Korea, and Singapore perform relatively better.


title: Mechanism of Instrumental Game Theory in The Legal Process via Stochastic Options Pricing Induction

abstract: Economic theory has provided an estimable intuition in understanding the perplexing ideologies in law, in the areas of economic law, tort law, contract law, procedural law and many others. Most legal systems require the parties involved in a legal dispute to exchange information through a process called discovery. The purpose is to reduce the relative optimisms developed by asymmetric information between the parties. Like a head or tail phenomenon in stochastic processes, uncertainty in the adjudication affects the decisions of the parties in a legal negotiation. This paper therefore applies the principles of aleatory analysis to determine how negotiations fail in the legal process, introduce the axiological concept of optimal transaction cost and formulates a numerical methodology based on backwards induction and stochastic options pricing economics in estimating the reasonable and fair bargain in order to induce settlements thereby increasing efficiency and reducing social costs.


title: Influence of Incremental Constraints on Energy Consumption and Static Scheduling Time for Moldable Tasks with Deadline

abstract: Static scheduling of independent, moldable tasks on parallel machines with frequency scaling comprises decisions on core allocation, assignment, frequency scaling and ordering, to meet a deadline and minimize energy consumption. Constraining some of these decisions reduces the solution space, i.e. may increase energy consumption, but may also reduce scheduling time or give the chance to tackle larger task sets. We investigate the influence of different constraints that lead from an unrestricted scheduler via two intermediate steps to the crown scheduler, by presenting integer linear programs for all four schedulers. We compare scheduling time and energy consumption for a benchmark suite of synthetic task sets of different sizes. Our results indicate that the final step towards the crown scheduler -- the execution order constraint -- is responsible for faster scheduling when task sets are small, and lower energy consumption when we deal with large task sets.


title: Dataset for Automatic Summarization of Russian News

abstract: Automatic text summarization has been studied in a variety of domains and languages. However, this does not hold for the Russian language. To overcome this issue, we present Gazeta, the first dataset for summarization of Russian news. We describe the properties of this dataset and benchmark several extractive and abstractive models. We demonstrate that the dataset is a valid task for methods of text summarization for Russian. Additionally, we prove the pretrained mBART model to be useful for Russian text summarization.


title: On interplay of surface tension and inertial stabilization mechanisms in the interface dynamics with interfacial mass flux

abstract: This work focuses on the interfacial dynamics with interfacial mass flux in the presence of acceleration and surface tension. We employ the general matrix method to find the fundamental solutions for the linearized boundary value problem conserving mass, momentum and energy. We find that the dynamics can be stable or unstable depending on the values of the acceleration, the surface tension and the density ratio. In the stable regime, the flow has the non-perturbed fields in the bulk, is shear-free at the interface, and has the constant interface velocity. The dynamics is unstable only when it is accelerated, and when the acceleration value exceeds a threshold combining contributions of the inertial stabilization mechanism and the surface tension. The properties of this instability unambiguously differentiate it from other fluid instabilities. Particularly, its velocity field has potential and vortical components in the bulk and is shear free at the interface. Its dynamics describes the standing wave with the growing amplitude, and has the growing interface velocity. For strong accelerations, this fluid instability of the conservative dynamics has the fastest growth-rate and the largest stabilizing surface tension value when compared to the classical Landau and Rayleigh-Taylor dynamics. We find the values of the initial perturbation wavelength at which the fluid instability can be stabilized and at which it has the fastest growth. We identify theory benchmarks for experiments and simulations in high energy density plasmas and its outcomes for application problems in nature and technology.


title: Elucidating the mechanism of step-emulsification

abstract: Three-dimensional, time-dependent direct simulations of step emulsification micro-devices highlight two essential mechanisms for droplet formation: first, the onset of an adverse pressure gradient driving a back-flow of the continuous phase from the external reservoir to the micro-channel. Second, the striction of the flowing jet which leads to its subsequent rupture. It is also shown that such a rupture is delayed and eventually suppressed by increasing the flow speed of the dispersed phase within the channel, due to the stabilising effect of dynamic pressure. This suggests a new criterion for dripping-jetting transition, based on local values of the Capillary and Weber numbers.


title: A Workflow for Exploring Ligand Dissociation from a Macromolecule: Efficient Random Acceleration Molecular Dynamics Simulation and Interaction Fingerprints Analysis of Ligand Trajectories

abstract: The dissociation of ligands from proteins and other biomacromolecules occurs over a wide range of timescales. For most pharmaceutically relevant inhibitors, these timescales are far beyond those that are accessible by conventional molecular dynamics (MD) simulation. Consequently, to explore ligand egress mechanisms and compute dissociation rates, it is necessary to enhance the sampling of ligand unbinding. Random Acceleration MD (RAMD) is a simple method to enhance ligand egress from a macromolecular binding site, which enables the exploration of ligand egress routes without prior knowledge of the reaction coordinates. Furthermore, the tauRAMD procedure can be used to compute the relative residence times of ligands. When combined with a machine-learning analysis of protein-ligand interaction fingerprints (IFPs), molecular features that affect ligand unbinding kinetics can be identified. Here, we describe the implementation of RAMD in GROMACS 2020, which provides significantly improved computational performance, with scaling to large molecular systems. For the automated analysis of RAMD results, we developed MD-IFP, a set of tools for the generation of IFPs along unbinding trajectories and for their use in the exploration of ligand dynamics. We demonstrate that the analysis of ligand dissociation trajectories by mapping them onto the IFP space enables the characterization of ligand dissociation routes and metastable states. The combined implementation of RAMD and MD-IFP provides a computationally efficient and freely available workflow that can be applied to hundreds of compounds in a reasonable computational time and will facilitate the use of tauRAMD in drug design.


title: Real-time reconstruction of long-lived particles at LHCb using FPGAs

abstract: Finding tracks downstream of the magnet at the earliest LHCb trigger level is not part of the baseline plan of the upgrade trigger, on account of the significant CPU time required to execute the search. Many long-lived particles, such as $K^0_S$ and strange baryons, decay after the vertex track detector, so that their reconstruction efficiency is limited. We present a study of the performance of a future innovative real-time tracking system based on FPGAs, developed within a R\&D effort in the context of the LHCb Upgrade Ib (LHC Run~4), dedicated to the reconstruction of the particles downstream of the magnet in the forward tracking detector (Scintillating Fibre Tracker), that is capable of processing events at the full LHC collision rate of 30 MHz.


title: The EuroSys 2020 Online Conference: Experience and lessons learned

abstract: The 15th European Conference on Computer Systems (EuroSys'20) was organized as a virtual (online) conference on April 27-30, 2020. The main EuroSys'20 track took place April 28-30, 2020, preceded by five workshops (EdgeSys'20, EuroDW'20, EuroSec'20, PaPoC'20, SPMA'20) on April 27, 2020. The decision to hold a virtual (online) conference was taken in early April 2020, after consultations with the EuroSys community and internal discussions about potential options, eventually allowing about three weeks for the organization. This paper describes the choices we made to organize EuroSys'20 as a virtual (online) conference, the challenges we addressed, and the lessons learned.


title: Stability of a parametric harmonic oscillator with dichotomic noise

abstract: The harmonic oscillator is a powerful model that can appear as a limit case when examining a nonlinear system. A well known fact is, that without driving, the inclusion of a friction term makes the origin of the phase space -- which is a fixpoint of the system -- linearly stable. In this work we include a telegraph process as perturbation of the oscillator's frequency, for example to describe the motion of a particle with fluctuating charge gyrating in an external magnetic field. Increasing intensity of this colored noise is capable of changing the quality of the fixed point. To characterize the stability of the system, we use a stability measure, that describes the growth of the displacement of the system's phase space position and express it in a closed form. We expand the respective exponent for light friction and low noise intensity and compare both, the exact analytic solution and the expansion to numerical values. Our findings allow stability predictions for several physical systems.


title: HPRA: Hyperedge Prediction using Resource Allocation

abstract: Many real-world systems involve higher-order interactions and thus demand complex models such as hypergraphs. For instance, a research article could have multiple collaborating authors, and therefore the co-authorship network is best represented as a hypergraph. In this work, we focus on the problem of hyperedge prediction. This problem has immense applications in multiple domains, such as predicting new collaborations in social networks, discovering new chemical reactions in metabolic networks, etc. Despite having significant importance, the problem of hyperedge prediction hasn't received adequate attention, mainly because of its inherent complexity. In a graph with $n$ nodes the number of potential edges is $\mathcal{O}(n^{2})$, whereas in a hypergraph, the number of potential hyperedges is $\mathcal{O}(2^{n})$. To avoid searching through such a huge space, current methods restrain the original problem in the following two ways. One class of algorithms assume the hypergraphs to be $k$-uniform. However, many real-world systems are not confined only to have interactions involving $k$ components. Thus, these algorithms are not suitable for many real-world applications. The second class of algorithms requires a candidate set of hyperedges from which the potential hyperedges are chosen. In the absence of domain knowledge, the candidate set can have $\mathcal{O}(2^{n})$ possible hyperedges, which makes this problem intractable. We propose HPRA - Hyperedge Prediction using Resource Allocation, the first of its kind algorithm, which overcomes these issues and predicts hyperedges of any cardinality without using any candidate hyperedge set. HPRA is a similarity-based method working on the principles of the resource allocation process. In addition to recovering missing hyperedges, we demonstrate that HPRA can predict future hyperedges in a wide range of hypergraphs.


title: Distributed prediction of unsafe reconfiguration scenarios of modular robotic Programmable Matter

abstract: We present a distributed framework for predicting whether a planned reconfiguration step of a modular robot will mechanically overload the structure, causing it to break or lose stability under its own weight. The algorithm is executed by the modular robot itself and based on a distributed iterative solution of mechanical equilibrium equations derived from a simplified model of the robot. The model treats inter-modular connections as beams and assumes no-sliding contact between the modules and the ground. We also provide a procedure for simplified instability detection. The algorithm is verified in the Programmable Matter simulator VisibleSim, and in real-life experiments on the modular robotic system Blinky Blocks.


title: Constraining HeII Reionization Detection Uncertainties via Fast Radio Bursts

abstract: Context. The increased detection rate of Fast Radio Bursts (FRBs) makes it likely to get samples of sizes $\mathcal{O}(10^2)$ to $\mathcal{O}(10^3)$ in the near future. Because of their extragalactic origin can help us in understanding the epoch of helium reionization. Aims. We try to identify the epoch of Helium II (HeII) reionization, via the observations of early FRBs in range of $z=3$ to $4$. Methods. We build a model of FRB Dispersion Measure following the HeII reionization model, density fluctuation in large scale structure, host galaxy interstellar medium and local environment of FRB contribution. The model is fit to the ideal intergalactic medium (IGM) dispersion measure model to check the goodness of constraining the HeII reionization via FRB measurement statistics. Conclusion. We report our findings under two categories, accuracy in detection of HeII reionization via FRBs assuming no uncertainty in the redshift measurement and alternatively assuming a varied level of uncertainty in redshift measurement of the FRBs. We show that under the first case, a detection of $N\sim\mathcal{O} (10^2)$ FRBs give an uncertainty of $\sigma (z_{r, fit})\sim0.5$ from the fit model, and a detection of $N\sim\mathcal{O} (10^3)$ gives an uncertainty of $\sigma (z_{r, fit})\sim0.1$. While assuming a redshift uncertainty of level $5-20\%$, changes the $\sigma (z_{r, fit})\sim0.5$ to $0.6$ for $N\sim 100$ and $\sigma (z_{r, fit})\sim0.1$ to $0.15$ for $N \sim 1000$ case.


title: Characterising Ion-Irradiated FeCr: Hardness, Thermal Diffusivity and Lattice Strain

abstract: Ion-irradiated FeCr alloys are useful for understanding and predicting neutron-damage in the structural steels of future nuclear reactors. Previous studies have largely focused on the structure of irradiation-induced defects, probed by transmission electron microscopy (TEM), as well as changes in mechanical properties. Across these studies, a wide range of irradiation conditions has been employed on samples with different processing histories, which complicates the analysis of the relationship between defect structures and material properties. Furthermore, key properties, such as irradiation-induced changes in thermal transport and lattice strain, are little explored. Here we present a systematic study of Fe3Cr, Fe5Cr and Fe10Cr binary alloys implanted with 20 MeV Fe$^{3+}$ ions to nominal doses of 0.01 dpa and 0.1 dpa at room temperature. Nanoindentation, transient grating spectroscopy (TGS) and X-ray micro-beam Laue diffraction were used to study the changes in hardness, thermal diffusivity and strain in the material as a function of damage and Cr content. Our results suggest that Cr leads to an increased retention of irradiation-induced defects, causing substantial changes in hardness and lattice strain. However, thermal diffusivity varies little with increasing damage and instead degrades significantly with increasing Cr content in the material. We find significant lattice strains even in samples exposed to a nominal displacement damage of 0.01 dpa. The defect density predicted from the lattice strain measurements is significantly higher than that observed in previous TEM studies, suggesting that TEM may not fully capture the irradiation-induced defect population.


title: On the growth of linear recurrences in function fields

abstract: Let $ (G_n)_{n=0}^{\infty} $ be a non-degenerate linear recurrence sequence with power sum representation $ G_n = a_1(n) \alpha_1^n + \cdots + a_t(n) \alpha_t^n $. In this paper we will prove a function field analogue of the well known result that in the number field case, under some non-restrictive conditions, for $ n $ large enough the inequality $ \vert G_n\vert \geq \left( \max_{j=1,\ldots,t} \vert \alpha_j\vert \right)^{n(1-\varepsilon)} $ holds true.


title: Norm form equations with solutions taking values in a multi-recurrence

abstract: We are interested in solutions of a norm form equation that takes values in a given multi-recurrence. We show that among the solutions there are only finitely many values in each component which lie in the given multi-recurrence unless the recurrence is of precisely described exceptional shape. This gives a variant of the question on arithmetic progressions in the solution set of norm form equations.


title: Dominant tournament families

abstract: For a tournament $H$ with $h$ vertices, its typical density is $h!2^{-\binom{h}{2}}/aut(H)$, i.e. this is the expected density of $H$ in a random tournament. A family ${\mathcal F}$ of $h$-vertex tournaments is {\em dominant} if for all sufficiently large $n$, there exists an $n$-vertex tournament $G$ such that the density of each element of ${\mathcal F}$ in $G$ is larger than its typical density by a constant factor. Characterizing all dominant families is challenging already for small $h$. Here we characterize several large dominant families for every $h$. In particular, we prove the following for all $h$ sufficiently large: (i) For all tournaments $H^*$ with at least $5\log h$ vertices, the family of all $h$-vertex tournaments that contain $H^*$ as a subgraph is dominant. (ii) The family of all $h$-vertex tournaments whose minimum feedback arc set size is at most $\frac{1}{2}\binom{h}{2}-h^{3/2}\sqrt{\ln h}$ is dominant. For small $h$, we construct a dominant family of $6$ (i.e. $50\%$ of the) tournaments on $5$ vertices and dominant families of size larger than $40\%$ for $h=6,7,8,9$. For all $h$, we provide an explicit construction of a dominant family which is conjectured to obtain an absolute constant fraction of the tournaments on $h$ vertices. Some additional intriguing open problems are presented.


title: A Better Alternative to Error Feedback for Communication-Efficient Distributed Learning

abstract: Modern large-scale machine learning applications require stochastic optimization algorithms to be implemented on distributed compute systems. A key bottleneck of such systems is the communication overhead for exchanging information across the workers, such as stochastic gradients. Among the many techniques proposed to remedy this issue, one of the most successful is the framework of compressed communication with error feedback (EF). EF remains the only known technique that can deal with the error induced by contractive compressors which are not unbiased, such as Top-$K$. In this paper, we propose a new and theoretically and practically better alternative to EF for dealing with contractive compressors. In particular, we propose a construction which can transform any contractive compressor into an induced unbiased compressor. Following this transformation, existing methods able to work with unbiased compressors can be applied. We show that our approach leads to vast improvements over EF, including reduced memory requirements, better communication complexity guarantees and fewer assumptions. We further extend our results to federated learning with partial participation following an arbitrary distribution over the nodes, and demonstrate the benefits thereof. We perform several numerical experiments which validate our theoretical findings.


title: Differentiable Language Model Adversarial Attacks on Categorical Sequence Classifiers

abstract: An adversarial attack paradigm explores various scenarios for the vulnerability of deep learning models: minor changes of the input can force a model failure. Most of the state of the art frameworks focus on adversarial attacks for images and other structured model inputs, but not for categorical sequences models. Successful attacks on classifiers of categorical sequences are challenging because the model input is tokens from finite sets, so a classifier score is non-differentiable with respect to inputs, and gradient-based attacks are not applicable. Common approaches deal with this problem working at a token level, while the discrete optimization problem at hand requires a lot of resources to solve. We instead use a fine-tuning of a language model for adversarial attacks as a generator of adversarial examples. To optimize the model, we define a differentiable loss function that depends on a surrogate classifier score and on a deep learning model that evaluates approximate edit distance. So, we control both the adversability of a generated sequence and its similarity to the initial sequence. As a result, we obtain semantically better samples. Moreover, they are resistant to adversarial training and adversarial detectors. Our model works for diverse datasets on bank transactions, electronic health records, and NLP datasets.


title: Conserved quantities, continuation and compactly supported solutions of some shallow water models

abstract: A proof that strong solutions of the Dullin-Gottwald-Holm equation vanishing on an open set of the (1+1) space-time are identically zero is presented. In order to do it, we use a geometrical approach based on the conserved quantities of the equation to prove a unique continuation result for its solutions. We show that this idea can be applied to a large class of equations of the Camassa-Holm type.


title: Bragg diffraction of large organic molecules

abstract: We demonstrate Bragg diffraction of the antibiotic ciprofloxacin and the dye molecule phthalocyanine at a thick optical grating. The observed patterns show a single dominant diffraction order with the expected dependence on the incidence angle as well as oscillating population transfer between the undiffracted and diffracted beams. We achieve an equal-amplitude splitting of $14 \hbar k$ (photon momenta) and maximum momentum transfer of $18 \hbar k$. This paves the way for efficient, large-momentum beam splitters and mirrors for hot and complex molecules.


title: On the time-dependent electrolyte Seebeck effect

abstract: Single-ion Soret coefficients $\alpha_{i}$ characterize the tendency of ions in an electrolyte solution to move in a thermal gradient. When these coefficients differ between cations and anions, an electric field can be generated. For this so-called electrolyte Seebeck effect to occur, the different thermodiffusive fluxes need to be blocked by boundaries -- electrodes, for example. Local charge neutrality is then broken in the Debye-length vicinity of the electrodes. Confusingly, many authors point to these regions as the source of the thermoelectric field yet ignore them in derivations of the time-dependent Seebeck coefficient $S(t)$, giving a false impression that the electrolyte Seebeck effect is purely a bulk phenomenon. Without enforcing local electroneutrality, we derive $S(t)$ generated by a binary electrolyte with arbitrary ionic valencies subject to a time-dependent thermal gradient. Next, we experimentally measure $S(t)$ for five acids, bases, and salts near titanium electrodes. For the steady state we find $S\approx2~\mathrm{mV~K}^{-1}$ for many electrolytes, roughly one order of magnitude larger than predictions based on literature $\alpha_{i}$. We fit our expression for $S(t)$ to the experimental data, treating the $\alpha_{i}$ as fit parameters, and also find larger-than-literature values, accordingly.


title: Electron-positron annihilation to photons at $O(\alpha^3)$ revisited

abstract: We apply the modern multiloop methods to the calculation of the total cross sections of electron-positron annihilation to 2 and 3 photons exactly in $s/m^2$ with the accuracy $O(\alpha^3)$. Examining the asymptotics of our results, we find agreement with Ref. [Andreassi_et_al_1962] and discover mistakes in the results of Refs. [Eidelman&Kuraev_1978,Berends&Kleiss_1981]. This mistake is due to the terms, omitted in differential cross section in Refs. \[Eidelman&Kuraev_1978,Berends&Kleiss_1981], which are peaked in the kinematic region with all three photons being quasi-parallel to the collision axis. After restoring these terms, we find an agreement of the corrected result of Ref. [Berends&Kleiss_1981] with our result.


title: Tidal Love numbers of Proca stars

abstract: Proca stars (arXiv:1508.05395) are everywhere regular, asymptotically flat self-gravitating solitons composed of complex, massive vector bosons, forming a macroscopic, star-sized, Bose-Einstein condensate. They have been suggested as a classical, dynamical model of exotic compact objects, black hole mimickers and dark matter candidates. In spherical symmetry, they are qualitatively similar to their scalar cousins, the standard (scalar) boson stars in all respects studied so far. Here we study the tidal deformability and the quadrupolar tidal Love numbers, both electric-type and magnetic-type, of spherically symmetric Proca stars. The equations for the perturbations are derived and the numerical values of the Love numbers for some concrete background solutions are computed. We observe that both Love numbers are qualitatively similar to the ones of the scalar boson stars; in particular, the electric-type (magnetic-type) quadrupolar Love numbers are positive (negative) for Proca stars. Quantitatively, for the same compactness, the electric and magnetic Love numbers in the Proca case are closer (in magnitude) than in the scalar case, with the electric Love numbers having a slightly larger magnitude than the magnetic ones.


title: Plastic deformation of rough metallic surfaces

abstract: We present experimental and theoretical results for the surface topography of a plastically deformed metallic (aluminum) block. When a hard spherical body (here a steel-, silica glass- or silicon nitride ball) with a smooth surface is indented in a metal block with a nominally flat, but still rough, surface, a spherical-cup-like indentation result due to plastic flow. The surface roughness in the indented region is, however, not entirely flattened. The long wavelength (macroasperity) content of the roughness result from the roughness on the original (aluminum) surface, but now plastically deformed. The roughness at short length scale, in the plastically deformed macroasperity contact regions, result from the roughness on the hard ball, and from inhomogeneous plastic flow. We model the contact mechanics using the boundary element method, combined with a simple numerical procedure to take into account the plastic flow. The theory can semi-quantitatively describe the modification of the roughness by the plastic flow. Since the fluid leakage of metallic seals in most cases is determined by the long wavelength roughness, we predict that the leakage can be estimated based on the elastoplastic contact mechanics model employed here. The plastic deformations of surfaces of some glassy polymers is very different from what we observed for aluminum, which we attribute to strong work-hardening and to inhomogeneous plastic flow for the polymers. Thus the numerical procedure to account for the plastic flow proposed here cannot be applied to glassy polymers in general.


title: CPT-violating z=3 Horava-Lifshitz QED and generation of the Carroll-Field-Jackiw term

abstract: We study the $z=3$ Horava-Lifshitz QED with a CPT-breaking term, characterized by the axial vector $b_{\mu}$, and perform the one-loop calculations. Explicitly, we demonstrate that just as in the usual Lorentz-breaking QED, in our case, the Carroll-Field-Jackiw term is finite but ambiguous.


title: Parametrized test of parity-violating gravity using GWTC-1 events

abstract: Parity-violating (PV) gravity has recently attracted interest in several aspects. One of them is the axion-graviton coupling to test the axion-dark matter model. Moreover, by extending Chern-Simons (CS) gravity to include derivatives of a scalar field up to the second order, a more general class of PV gravity theory, which we call the CNCL model, has been proposed~[M. Crisostomi {\it et al.}, Phys. Rev. D, {\bf 97}, 044034 (2018)]. The model can be further extended by including even higher derivatives of the scalar field and/or higher curvature terms. In this paper, we discuss the effect of parity violation in the gravitational sector on the propagation of gravitational waves from binary coalescence by introducing a model-independent parametrization of modification. Our parametrization includes the CNCL model as well as CS gravity. The effect of parity violation on the gravitational waveform is maximum when the source binary orientation to our line of sight is edge-on, while the modified waveform reduces to the parity-symmetric one when the source is face-on. We perform a search for the signature of such modification by using the LIGO/Virgo O1/O2 catalog. We find that the catalog data is consistent with general relativity and obtain constraints on parity violation in gravity for various post-Newtonian order modifications for the first time. The obtained constraint on CS gravity is consistent with the results in previous works. On the other hand, the constraint on the CNCL model that we obtain is tighter than the previous results by roughly 7 orders of magnitude.


title: Existence of weak solutions for inhomogeneous generalized Navier-Stokes equations

abstract: We prove existence of weak solutions for the fully inhomogeneous, stationary generalized Navier-Stokes equations for shear-thinning fluids. Our proof is based on the theory of pseudomonotone operators and the Lipschitz truncation method, whose application is presented as a general result. Our approach requires a smallness and a regularity assumption on the data; we show that this is inevitable in the framework of pseudomonotone operators.


title: Time series copula models using d-vines and v-transforms

abstract: An approach to modelling volatile financial return series using stationary d-vine copula processes combined with Lebesgue-measure-preserving transformations known as v-transforms is proposed. By developing a method of stochastically inverting v-transforms, models are constructed that can describe both stochastic volatility in the magnitude of price movements and serial correlation in their directions. In combination with parametric marginal distributions it is shown that these models can rival and sometimes outperform well-known models in the extended GARCH family.


title: Three Lyman-alpha emitting filaments converging to a massive galaxy group at z=2.91: discussing the case for cold gas infall

abstract: We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^13Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The nebula morphological and kinematics properties and the prevalence of blueshifted components in the Lya spectra are consistent with a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also photoionization with ensuing recombination still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside an actively star-forming massive halo at high redshift.


title: Mapping quantum random walks onto a Markov chain by mapping a unitary transformation to a higher dimension of an irreducible matrix

abstract: Here, a new two-dimensional process, discrete in time and space, that yields the results of both a random walk and a quantum random walk, is introduced. This model describes the population distribution of four coin states |1>,-|1>, |0> -|0> in space without interference, instead of two coin states |1>, |0> .For the case of no boundary conditions, the model is similar to a Markov chain with a stochastic matrix, i.e., it conserves the population distribution of the four coin states, and by using a proper transformation, yield probability distributions of the two quantum states |1>, |0> in space, similar to a unitary operator. Numerical results for a quantum random walk on infinite and finite lines are introduced.


title: Deep Learning-based Single Image Face Depth Data Enhancement

abstract: Face recognition can benefit from the utilization of depth data captured using low-cost cameras, in particular for presentation attack detection purposes. Depth video output from these capture devices can however contain defects such as holes or general depth inaccuracies. This work proposes a deep learning face depth enhancement method in this context of facial biometrics, which adds a security aspect to the topic. U-Net-like architectures are utilized, and the networks are compared against hand-crafted enhancer types, as well as a similar depth enhancer network from related work trained for an adjacent application scenario. All tested enhancer types exclusively use depth data as input, which differs from methods that enhance depth based on additional input data such as visible light color images. Synthetic face depth ground truth images and degraded forms thereof are created with help of PRNet, to train multiple deep learning enhancer models with different network sizes and training configurations. Evaluations are carried out on the synthetic data, on Kinect v1 images from the KinectFaceDB, and on in-house RealSense D435 images. These evaluations include an assessment of the falsification for occluded face depth input, which is relevant to biometric security. The proposed deep learning enhancers yield noticeably better results than the tested preexisting enhancers, without overly falsifying depth data when non-face input is provided, and are shown to reduce the error of a simple landmark-based PAD method.


title: The absolute magnitudes $M_J$, the binary fraction, and the binary mass ratios of M7 to M9.5 dwarfs

abstract: We use the large homogeneous sample of late M dwarfs, M7 to M9.5, of Ahmed & Warren (2019) matched to $\textit{Gaia}$ DR2, to measure the relation between absolute magnitude and spectral type, and to infer the multiplicity fraction of the population, and the distribution of mass ratios in the binary systems. Binaries are identified photometrically as overluminous sources. In order to define a sample that is unbiased with respect to multiplicity we use distance limits that are a function of $G-J$ colour to define a volume-complete sample of 2706 systems. The $G-J$ colours are very precise, with random errors all less than 0.02. We measure absolute magnitudes $M_J$ that are on average 0.5 mag. brighter than previous determinations. We find evidence that the discrepancies arise from differences in spectral types in different samples. The measured binary fraction is $16.5\pm0.8\%$, of which $98\%$ are unresolved: both values are consistent with results of previous studies. The distribution of excess flux in the binaries, compared to the singles, is used to infer the mass ratio distribution $f(q)\propto q^\gamma$, where $q=M_s/M_p$. We infer a very steep distribution over this spectral range, with $\gamma>10$ ($99\%$ probability). This says that unresolved ultracool M dwarf binaries reside almost exclusively in equal mass systems, and implies that the spectral types of the unresolved binaries match to with 0.5 spectral subtypes. The intrinsic scatter in absolute magnitude $M_J$ for ultracool M dwarfs at fixed $G-J$ colour is measured to be 0.21 mag.


title: Managing spectral properties and Schmidt mode content of squeezed vacuum light using sum-frequency converter

abstract: Capabilities of quantum optical SFG-gate seeded by squeezed light are investigated in the frame of frequency Schmidt modes. Methods to manage and manipulate extensively the properties and mode content of squeezed light are developed. Possibilities to block and select any certain Schmidt mode of squeezed light with conservation of non-classical properties are demonstrated. The significant phase sensitivity of the gate is shown and the ways to manage the spectral distribution of the output light due to the phase effects and variable coupling between modes in the gate are demonstrated. The effect of swapping between modes in the gate is found. It allows to enhance squeezed light in a set of modes without loss of photon correlations which is important for further experiments and new applications.


title: Mixture of Conditional Gaussian Graphical Models for unlabelled heterogeneous populations in the presence of co-factors

abstract: Conditional correlation networks, within Gaussian Graphical Models (GGM), are widely used to describe the direct interactions between the components of a random vector. In the case of an unlabelled Heterogeneous population, Expectation Maximisation (EM) algorithms for Mixtures of GGM have been proposed to estimate both each sub-population's graph and the class labels. However, we argue that, with most real data, class affiliation cannot be described with a Mixture of Gaussian, which mostly groups data points according to their geometrical proximity. In particular, there often exists external co-features whose values affect the features' average value, scattering across the feature space data points belonging to the same sub-population. Additionally, if the co-features' effect on the features is Heterogeneous, then the estimation of this effect cannot be separated from the sub-population identification. In this article, we propose a Mixture of Conditional GGM (CGGM) that subtracts the heterogeneous effects of the co-features to regroup the data points into sub-population corresponding clusters. We develop a penalised EM algorithm to estimate graph-sparse model parameters. We demonstrate on synthetic and real data how this method fulfils its goal and succeeds in identifying the sub-populations where the Mixtures of GGM are disrupted by the effect of the co-features.


title: Vortices in nonequilibrium photon condensates

abstract: We present a theoretical study of vortices in arrays of photon condensates. Even when interactions are negligible, as is the case in current experiments, pumping and losses can lead to a finite vortex core size. While some properties of photon condensate vortices, such as their self-acceleration and the generation of vortex pairs by a moving vortex, resemble those in interacting polariton condensates far from equilibrium, in several aspects they differ from previously studied systems: the vortex core size is determined by the balance between pumping and tunneling, the core appears oblate in the direction of its motion and new vortex pairs can spontaneously nucleate in the core region.


title: A numerical algorithm to computationally solve the Hemker problem using Shishkin meshes

abstract: A numerical algorithm is presented to solve a benchmark problem proposed by Hemker. The algorithm incorporates asymptotic information into the design of appropriate piecewise-uniform Shishkin meshes. Moreover, different co-ordinate systems are utilized due to the different geometries and associated layer structures that are involved in this problem. Numerical results are presented to demonstrate the effectiveness of the proposed numerical algorithm.