title: An Evaluation of the Human-Interpretability of Explanation

abstract: Recent years have seen a boom in interest in machine learning systems that can provide a human-understandable rationale for their predictions or decisions. However, exactly what kinds of explanation are truly human-interpretable remains poorly understood. This work advances our understanding of what makes explanations interpretable under three specific tasks that users may perform with machine learning systems: simulation of the response, verification of a suggested response, and determining whether the correctness of a suggested response changes under a change to the inputs. Through carefully controlled human-subject experiments, we identify regularizers that can be used to optimize for the interpretability of machine learning systems. Our results show that the type of complexity matters: cognitive chunks (newly defined concepts) affect performance more than variable repetitions, and these trends are consistent across tasks and domains. This suggests that there may exist some common design principles for explanation systems.


title: Properties of minimal charts and their applications V: charts of type $(3,2,2)$

abstract: Let $\Gamma$ be a chart, and we denote by $\Gamma_m$ the union of all the edges of label $m$. A chart $\Gamma$ is of type $(3,2,2)$ if there exists a label $m$ such that $w(\Gamma)=7$, $w(\Gamma_m\cap\Gamma_{m+1})=3$, $w(\Gamma_{m+1}\cap\Gamma_{m+2})=2$, and $w(\Gamma_{m+2}\cap\Gamma_{m+3})=2$ where $w(G)$ is the number of white vertices in $G$. In this paper, we prove that there is no minimal chart of type $(3,2,2)$.


title: Analysis of hidden-bottom bb\bar{b}\bar{b} states

abstract: Motivated by the searching for $bb\bar{b}\bar{b}$ states at LHC recently, we calculate the ground-state energies of $bb\bar{b}\bar{b}$ states with quantum numbers $IJ^P=00^+,01^+,02^+$ in a nonrelativistic chiral quark model using the Gaussian expansion method. In our calculations, two structures, meson-meson and diquark-antidiquark, and their coupling, along with all possible color configurations are considered. It is expected that the studies shall be helpful for the experimental searching of fully-heavy exotic tetraquark states.


title: Descriptor system techniques and software tools

abstract: The role of the descriptor system representation as basis for reliable numerical computations for system analysis and synthesis, and in particular, for the manipulation of rational matrices, is discussed and available robust numerical software tools are described.


title: Algebraic surfaces with infinitely many twistor lines

abstract: We prove that a reduced and irreducible algebraic surface in $\mathbb{CP}^{3}$ containing infinitely many twistor lines cannot have odd degree. Then, exploiting the theory of quaternionic slice regularity and the normalization map of a surface, we give constructive existence results for even degrees.


title: Mathematics of nerve signals

abstract: Mathematical models describing the signals propagating in nerve fibres are described. The emphasis is on the mathematical structures of governing equations while the extremely rich physiological aspects are here not analysed. Based on models of single waves, a joint coupled model is presented which is able to describe the action potential and the accompanying mechanical effects togehter with temperature changes within one system of partial differential equations. The whole signal is an ensemble which includes primary and secondary components. The primary components of a signal are the action potential itself and longitudinal mechanical waves in axoplasm and surrounding biomembrane. These components are characterized by corresponding velocities. The secondary components of a signal are derived from primary components and include transverse displacement of a biomembrane and the temperature change. These secondary components have no independent velocities in the presented model.


title: A note on the Dirac operator with Kirchoff-type vertex conditions on noncompact metric graphs

abstract: In this note we present some properties of the Dirac operator on noncompact metric graphs with Kirchoff-type vertex conditions. In particular, we discuss the specific features of the spectrum of the operator and, finally, we give some further details on the associated quadratic form (and on the form domain).


title: Representations and the reduction theorem for ultragraph Leavitt path algebras

abstract: In this paper we study representations of ultragraph Leavitt path algebras via branching systems and, using partial skew ring theory, prove the reduction theorem for these algebras. We apply the reduction theorem to show that ultragraph Leavitt path algebras are semiprime and to completely describe faithfulness of the representations arising from branching systems, in terms of the dynamics of the branching systems. Furthermore, we study permutative representations and provide a sufficient criteria for a permutative representation of an ultragraph Leavitt path algebra to be equivalent to a representation arising from a branching system. We apply this criteria to describe a class of ultragraphs for which every representation (satisfying a mild condition) is permutative and has a restriction that is equivalent to a representation arising from a branching system.


title: Query Inseparability for ALC Ontologies

abstract: We investigate the problem whether two ALC ontologies are indistinguishable (or inseparable) by means of queries in a given signature, which is fundamental for ontology engineering tasks such as ontology versioning, modularisation, update, and forgetting. We consider both knowledge base (KB) and TBox inseparability. For KBs, we give model-theoretic criteria in terms of (finite partial) homomorphisms and products and prove that this problem is undecidable for conjunctive queries (CQs), but 2ExpTime-complete for unions of CQs (UCQs). The same results hold if (U)CQs are replaced by rooted (U)CQs, where every variable is connected to an answer variable. We also show that inseparability by CQs is still undecidable if one KB is given in the lightweight DL EL and if no restrictions are imposed on the signature of the CQs. We also consider the problem whether two ALC TBoxes give the same answers to any query over any ABox in a given signature and show that, for CQs, this problem is undecidable, too. We then develop model-theoretic criteria for Horn-ALC TBoxes and show using tree automata that, in contrast, inseparability becomes decidable and 2ExpTime-complete, even ExpTime-complete when restricted to (unions of) rooted CQs.


title: On "Quantum interference with slits" and its "revisited"

abstract: Marcella in 2002 published a "quantum-mechanical" treatment of the famous single- and double-slit interference experiment in classical wave optics by a simple assumption that quantum mechanical wave function is a constant anywhere within a slit. Rothman and Boughn in 2011 commented that Marcella introduced no quantum physics into the problem other than\ a symbol substitution $p=\hbar k$, and the used essentially the classical wave optics. In present comment, we point out that though Marcella made a fundamental mistake, the problem is nevertheless remediable to give the satisfactory quantum mechanical results which are in qualitatively agreement with experimental ones with material particles.


title: Network Parameter Learning Using Nonlinear Transforms, Local Representation Goals and Local Propagation Constraints

abstract: In this paper, we introduce a novel concept for learning of the parameters in a neural network. Our idea is grounded on modeling a learning problem that addresses a trade-off between (i) satisfying local objectives at each node and (ii) achieving desired data propagation through the network under (iii) local propagation constraints. We consider two types of nonlinear transforms which describe the network representations. One of the nonlinear transforms serves as activation function. The other one enables a locally adjusted, deviation corrective components to be included in the update of the network weights in order to enable attaining target specific representations at the last network node. Our learning principle not only provides insight into the understanding and the interpretation of the learning dynamics, but it offers theoretical guarantees over decoupled and parallel parameter estimation strategy that enables learning in synchronous and asynchronous mode. Numerical experiments validate the potential of our approach on image recognition task. The preliminary results show advantages in comparison to the state-of-the-art methods, w.r.t. the learning time and the network size while having competitive recognition accuracy.


title: Dual-band Dielectric Light-harvesting Nanoantennae Made by Nature

abstract: Mechanisms to use nanoparticles to separate sunlight into photovoltaic useful range and thermally useful range to increase the efficiency of solar cells and to dissipate heat radiatively are discussed based upon lessons we learnt from photosynthesis. We show that the dual-band maxima in the absorption spectrum of bacterial light harvestors not only are due to the bacteriochlorophylls involved but also come from the geometry of the light harvestor. Being able to manipulate these two bands arbitrarily enables us to fabricate the nanoparticles required. Such mechanisms are also useful for the design of remote power charging and light sensors.

categories: physics.optics cond-mat.mes-hall physics.bio-ph q-bio.BM quant-ph


title: Modeling the Mitral Valve

abstract: This work is concerned with modeling and simulation of the mitral valve, one of the four valves in the human heart. The valve is composed of leaflets, the free edges of which are supported by a system of chordae, which themselves are anchored to the papillary muscles inside the left ventricle. First, we examine valve anatomy and present the results of original dissections. These display the gross anatomy and information on fiber structure of the mitral valve. Next, we build a model valve following a design-based methodology, meaning that we derive the model geometry and the forces that are needed to support a given load, and construct the model accordingly. We incorporate information from the dissections to specify the fiber topology of this model. We assume the valve achieves mechanical equilibrium while supporting a static pressure load. The solution to the resulting differential equations determines the pressurized configuration of the valve model. To complete the model we then specify a constitutive law based on a stress-strain relation consistent with experimental data that achieves the necessary forces computed in previous steps. Finally, using the immersed boundary method, we simulate the model valve in fluid in a computer test chamber. The model opens easily and closes without leak when driven by physiological pressures over multiple beats. Further, its closure is robust to driving pressures that lack atrial systole or are much lower or higher than normal.


title: In Search of the Thermal Eccentricity Distribution

abstract: About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.


title: On dual codes in the Doob schemes

abstract: The Doob scheme $D(m,n'+n'')$ is a metric association scheme defined on $E_4^m \times F_4^{n'}\times Z_4^{n''}$, where $E_4=GR(4^2)$ or, alternatively, on $Z_4^{2m} \times Z_2^{2n'} \times Z_4^{n''}$. We prove the MacWilliams identities connecting the weight distributions of a linear or additive code and its dual. In particular, for each case, we determine the dual scheme, on the same set but with different metric, such that the weight distribution of an additive code $C$ in the Doob scheme $D(m,n'+n'')$ is related by the MacWilliams identities with the weight distribution of the dual code $C^\perp$ in the dual scheme. We note that in the case of a linear code $C$ in $E_4^m \times F_4^{n'}$, the weight distributions of $C$ and $C^\perp$ in the same scheme are also connected.


title: Multiband gravitational-wave event rates and stellar physics

abstract: Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected number of multiband detections first by using information from current LIGO/Virgo data, and then through population synthesis simulations of binary stars. We estimate that few to tens of LISA detections can be used to predict mergers detectable on the ground. Conversely, hundreds of events could potentially be extracted from the LISA data stream using prior information from ground detections. In general, the merger signal of binaries observable by LISA is strong enough to be unambiguously identified by both current and future ground-based detectors. Therefore third-generation detectors will not increase the number of multiband detections compared to LIGO/Virgo. We use population synthesis simulations of isolated binary stars to explore some of the stellar physics that could be constrained with multiband events, and we show that specific formation pathways might be overrepresented in multiband events compared to ground-only detections.


title: On $(2n/3-1)$-resilient $(n,2)$-functions

abstract: A $\{00,01,10,11\}$-valued function on the vertices of the $n$-cube is called a $t$-resilient $(n,2)$-function if it has the same number of $00$s, $01$s, $10$s and $11$s among the vertices of every subcube of dimension $t$. The Friedman and Fon-Der-Flaass bounds on the correlation immunity order say that such a function must satisfy $t\le 2n/3-1$; moreover, the $(2n/3-1)$-resilient $(n,2)$-functions correspond to the equitable partitions of the $n$-cube with the quotient matrix $[[0,r,r,r],[r,0,r,r],[r,r,0,r],[r,r,r,0]]$, $r=n/3$. We suggest constructions of such functions and corresponding partitions, show connections with Latin hypercubes and binary $1$-perfect codes, characterize the non-full-rank and the reducible functions from the considered class, and discuss the possibility to make a complete characterization of the class.


title: On multifold packings of radius-1 balls in Hamming graphs

abstract: A $\lambda$-fold $r$-packing (multiple radius-$r$ covering) in a Hamming metric space is a code $C$ such that the radius-$r$ balls centered in $C$ cover each vertex of the space by not more (not less, respectively) than $\lambda$ times. The well-known $r$-error-correcting codes correspond to the case $\lambda=1$, while in general multifold $r$-packing are related with list decodable codes. We (a) propose asymptotic bounds for the maximum size of a $q$-ary $2$-fold $1$-packing as $q$ grows; (b) prove that a $q$-ary distance-$2$ MDS code of length $n$ is an optimal $n$-fold $1$-packing if $q\ge 2n$; (c) derive an upper bound for the size of a binary $\lambda$-fold $1$-packing and a lower bound for the size of a binary multiple radius-$1$ covering (the last bound allows to update the small-parameters table); (d) classify all optimal binary $2$-fold $1$-packings up to length $9$, in particular, establish the maximum size $96$ of a binary $2$-fold $1$-packing of length $9$; (e) prove some properties of $1$-perfect unitrades, which are a special case of $2$-fold $1$-packings. Keywords: Hamming graph, multifold ball packings, two-fold ball packings, list decodable codes, multiple coverings, completely regular codes, linear programming bound


title: Constraining Gas Motions in the Intra-Cluster Medium

abstract: The detailed velocity structure of the diffuse X-ray emitting intra-cluster medium (ICM) remains one of the last missing key ingredients in understanding the microphysical properties of these hot baryons and constraining our models of the growth and evolution of structure on the largest scales in the Universe. Direct measurements of the gas velocities from the widths and shifts of X-ray emission lines were recently provided for the central region of the Perseus Cluster of galaxies by $Hitomi$, and upcoming high-resolution X-ray microcalorimeters onboard $XRISM$ and $Athena$ are expected to extend these studies to many more systems. In the mean time, several other direct and indirect methods have been proposed for estimating the velocity structure in the ICM, ranging from resonant scattering to X-ray surface brightness fluctuation analysis, the kinematic Sunyaev-Zeldovich effect, or using optical line emitting nebulae in the brightest cluster galaxies as tracers of the motions of the ambient plasma. Here, we review and compare the existing estimates of the velocities of the hot baryons, as well as the various overlapping physical processes that drive motions in the ICM, and discuss the implications of these measurements for constraining the viscosity and identifying the source of turbulence in clusters of galaxies.


title: Localization as an entanglement phase transition in boundary-driven Anderson models

abstract: The Anderson localization transition is one of the most well studied examples of a zero temperature quantum phase transition. On the other hand, many open questions remain about the phenomenology of disordered systems driven far out of equilibrium. Here we study the localization transition in the prototypical three-dimensional, noninteracting Anderson model when the system is driven at its boundaries to induce a current carrying non-equilibrium steady state. Recently we showed that the diffusive phase of this model exhibits extensive mutual information of its non-equilibrium steady-state density matrix. We show that that this extensive scaling persists in the entanglement and at the localization critical point, before crossing over to a short-range (area-law) scaling in the localized phase. We introduce an entanglement witness for fermionic states that we name the mutual coherence, which, for fermionic Gaussian states, is also a lower bound on the mutual information. Through a combination of analytical arguments and numerics, we determine the finite-size scaling of the mutual coherence across the transition. These results further develop the notion of entanglement phase transitions in open systems, with direct implications for driven many-body localized systems, as well as experimental studies of driven-disordered systems.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.mes-hall cond-mat.quant-gas quant-ph


title: Spectrum of light- and heavy-baryons

abstract: A symmetry-preserving truncation of the strong-interaction bound-state equations is used to calculate the spectrum of ground-state $J=1/2^+$, $3/2^+$ $(qq^\prime q^{\prime\prime})$-baryons, where $q, q^\prime, q^{\prime\prime} \in \{u,d,s,c,b\}$, their first positive-parity excitations and parity partners. Using two parameters, a description of the known spectrum of 39 such states is obtained, with a mean-absolute-relative-difference between calculation and experiment of 3.6(2.7)%. From this foundation, the framework is subsequently used to predict the masses of 90 states not yet seen empirically.


title: Towards Machine-assisted Meta-Studies: The Hubble Constant

abstract: We present an approach for automatic extraction of measured values from the astrophysical literature, using the Hubble constant for our pilot study. Our rules-based model -- a classical technique in natural language processing -- has successfully extracted 298 measurements of the Hubble constant, with uncertainties, from the 208,541 available arXiv astrophysics papers. We have also created an artificial neural network classifier to identify papers in arXiv which report novel measurements. From the analysis of our results we find that reporting measurements with uncertainties and the correct units is critical information when distinguishing novel measurements in free text. Our results correctly highlight the current tension for measurements of the Hubble constant and recover the $3.5\sigma$ discrepancy -- demonstrating that the tool presented in this paper is useful for meta-studies of astrophysical measurements from a large number of publications.


title: Supernovae feedback propagation: the role of turbulence

abstract: Modelling the propagation of supernova (SN) bubbles, in terms of energy, momentum and spatial extent, is critical for simulations of galaxy evolution which do not capture these scales. To date, small scale models of SN feedback predict that the evolution of above-mentioned quantities can be solely parameterised by average quantities of the surrounding gas, such as density. However, most of these studies neglect the turbulent motions of this medium. In this paper, we study the propagation and evolution of SNe in turbulent environments. We confirm that the time evolution of injected energy and momentum can be characterised by the average density. However, the details of the density structure of the interstellar medium play a crucial role in the spatial extent of the bubble, even at a given average density. We demonstrate that spherically symmetric models of SN bubbles do not model well their spatial extent, and therefore cannot not be used to design sub-grid models of SNe feedback at galactic and cosmological scales.


title: Antiferromagnetism and chiral d-wave superconductivity from an effective $t-J-D$ model for twisted bilayer graphene

abstract: Starting from the strong-coupling limit of an extended Hubbard model, we develop a spin-fermion theory to study the insulating phase and pairing symmetry of the superconducting phase in twisted bilayer graphene. Assuming that the insulating phase is an anti-ferromagnetic insulator, we show that fluctuations of the anti-ferromagnetic order in the conducting phase can mediate superconducting pairing. Using a self-consistent mean-field analysis, we find that the pairing wave function has a chiral d-wave symmetry. Consistent with this observation, we show explicitly the existence of chiral Majorana edge modes by diagonalizing our proposed Hamiltonian on a finite-sized system. These results establish twisted bilayer graphene as a promising platform to realize topological superconductivity.


title: The Three Causes of Low-Mass Assembly Bias

abstract: We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masses and for alternative formation time definitions. We show that splashback subhaloes are responsible for two thirds of the assembly bias signal, but do not account for the entire effect. After splashback subhaloes have been removed, we find that the remaining assembly bias signal is due to a relatively small fraction ($\lesssim 10\%$) of haloes in dense regions. We test a number of additional physical processes thought to contribute to assembly bias and demonstrate that the two key processes are the slowing of mass growth by large-scale tidal fields and by the high velocities of ambient matter in sheets and filaments. We also rule out several other proposed physical causes of halo assembly bias. Based on our results, we argue that there are three processes that contribute to assembly bias of low-mass halos: large-scale tidal fields, gravitational heating due to the collapse of large-scale structures, and splashback subhaloes located outside the virial radius.


title: Bridging the Semantic Gap with SQL Query Logs in Natural Language Interfaces to Databases

abstract: A critical challenge in constructing a natural language interface to database (NLIDB) is bridging the semantic gap between a natural language query (NLQ) and the underlying data. Two specific ways this challenge exhibits itself is through keyword mapping and join path inference. Keyword mapping is the task of mapping individual keywords in the original NLQ to database elements (such as relations, attributes or values). It is challenging due to the ambiguity in mapping the user's mental model and diction to the schema definition and contents of the underlying database. Join path inference is the process of selecting the relations and join conditions in the FROM clause of the final SQL query, and is difficult because NLIDB users lack the knowledge of the database schema or SQL and therefore cannot explicitly specify the intermediate tables and joins needed to construct a final SQL query. In this paper, we propose leveraging information from the SQL query log of a database to enhance the performance of existing NLIDBs with respect to these challenges. We present a system Templar that can be used to augment existing NLIDBs. Our extensive experimental evaluation demonstrates the effectiveness of our approach, leading up to 138% improvement in top-1 accuracy in existing NLIDBs by leveraging SQL query log information.


title: Categorical Semantics for Time Travel

abstract: We introduce a general categorical framework to reason about quantum theory and other process theories living in spacetimes where Closed Timelike Curves (CTCs) are available, allowing resources to travel back in time and provide computational speedups. Our framework is based on a weakening of the definition of traced symmetric monoidal categories, obtained by dropping the yanking axiom and the requirement that the trace be defined on all morphisms. We show that the two leading models for quantum theory with closed timelike curves---namely the P-CTC model of Lloyd et al. and the D-CTC model of Deutsch---are captured by our framework, and in doing so we provide the first compositional description of the D-CTC model. Our description of the D-CTC model results in a process theory which respects the constraints of relativistic causality: this is in direct contrast to the P-CTC model, where CTCs are implemented by a trace and allow post-selection to be performed deterministically.


title: Compressed Diffusion

abstract: Diffusion maps are a commonly used kernel-based method for manifold learning, which can reveal intrinsic structures in data and embed them in low dimensions. However, as with most kernel methods, its implementation requires a heavy computational load, reaching up to cubic complexity in the number of data points. This limits its usability in modern data analysis. Here, we present a new approach to computing the diffusion geometry, and related embeddings, from a compressed diffusion process between data regions rather than data points. Our construction is based on an adaptation of the previously proposed measure-based Gaussian correlation (MGC) kernel that robustly captures the local geometry around data points. We use this MGC kernel to efficiently compress diffusion relations from pointwise to data region resolution. Finally, a spectral embedding of the data regions provides coordinates that are used to interpolate and approximate the pointwise diffusion map embedding of data. We analyze theoretical connections between our construction and the original diffusion geometry of diffusion maps, and demonstrate the utility of our method in analyzing big datasets, where it outperforms competing approaches.


title: Impurity states and Localization in Bilayer Graphene: the Low Impurity Concentration Regime

abstract: We study the problem of non-magnetic impurities adsorbed on bilayer graphene in the diluted regime. We analyze the impurity spectral densities for various concentrations and gate fields. We also analyze the effect of the adsorbate on the local density of states (LDOS) of the different C atoms in the structure and present some evidence of strong localization for the electronic states with energies close to the Dirac point.


title: Predictability of missing links in complex networks

abstract: Predicting missing links in real networks is an important problem in network science to which considerable efforts have been devoted, giving as a result a vast plethora of link prediction methods in the literature. In this work, we take a different point of view on the problem and study the theoretical limitations to the predictability of missing links. In particular, we hypothesise that there is an irreducible uncertainty in link prediction on real networks as a consequence of the random nature of their formation process. By considering ensembles defined by well-known network models, we prove analytically that even the best possible link prediction method for an ensemble, given by the ranking of the ensemble connection probabilities, yields a limited precision. This result suggests a theoretical limitation to the predictability of links in real complex networks. Finally, we show that connection probabilities inferred by fitting network models to real networks allow to estimate an upper-bound to the predictability of missing links, and we further propose a method to approximate such bound from incomplete instances of real-world networks.


title: Scattered light shadows in warped protoplanetary discs

abstract: Three-dimensional hydrodynamic numerical simulations have demonstrated that the structure of a protoplanetary disc may be strongly affected by a planet orbiting in a plane that is misaligned to the disc. When the planet is able to open a gap, the disc is separated into an inner, precessing disc and an outer disc with a warp. In this work, we compute infrared scattered light images to investigate the observational consequences of such an arrangement. We find that an inner disc misaligned by a less than a degree to the outer disc is indeed able to cast a shadow at larger radii. In our simulations a planet of around 6 Jupiter masses inclined by around 2 degrees is enough to warp the disc and cast a shadow with a depth of more than 10% of the average flux at that radius. We also demonstrate that warp in the outer disc can cause a variation in the azimuthal brightness profile at large radii. Importantly, this latter effect is a function of the distance from the star and is most prominent in the outer disc. We apply our model to the TW Hya system, where a misaligned, precessing inner disc has been invoked to explain an recently observed shadow in the outer disc. Consideration of the observational constraints suggest that an inner disc precessing due to a misaligned planet is an unlikely explanation for the features found in TW Hya.


title: Signatures of Bimodality in Nebular Phase Type Ia Supernova Spectra

abstract: One observational prediction for Type Ia supernova (SNe Ia) explosions produced through white dwarf-white dwarf collisions is the presence of bimodal velocity distributions for the $^{56}$Ni decay products, although this signature can also be produced by an off-center ignition in a delayed detonation explosion. These bimodal velocity distributions can manifest as double-peaked or flat-topped spectral features in late-time spectroscopic observations for favorable viewing angles. We present nebular-phase spectroscopic observations of 17 SNe Ia obtained with the Large Binocular Telescope (LBT). Combining these observations with an extensive search of publicly available archival data, we collect a total sample of 48 SNe Ia and classify them based on whether they show compelling evidence for bimodal velocity profiles in three features associated with $^{56}$Ni decay products: the [Fe II] and [Fe III] feature at $\sim5300$ \r{A}, the [Co III] $\lambda5891$ feature, and the [Co III] and [Fe II] feature at $\sim6600$ \r{A}. We identify 9 bimodal SNe in our sample, and we find that these SNe have average peak $M_V$ about 0.3 mag fainter than those which do not. This is consistent with theoretical predictions for explosions created by nearly head-on collisions of white dwarfs due to viewing angle effects and $^{56}$Ni yields.


title: BLOCK: Bilinear Superdiagonal Fusion for Visual Question Answering and Visual Relationship Detection

abstract: Multimodal representation learning is gaining more and more interest within the deep learning community. While bilinear models provide an interesting framework to find subtle combination of modalities, their number of parameters grows quadratically with the input dimensions, making their practical implementation within classical deep learning pipelines challenging. In this paper, we introduce BLOCK, a new multimodal fusion based on the block-superdiagonal tensor decomposition. It leverages the notion of block-term ranks, which generalizes both concepts of rank and mode ranks for tensors, already used for multimodal fusion. It allows to define new ways for optimizing the tradeoff between the expressiveness and complexity of the fusion model, and is able to represent very fine interactions between modalities while maintaining powerful mono-modal representations. We demonstrate the practical interest of our fusion model by using BLOCK for two challenging tasks: Visual Question Answering (VQA) and Visual Relationship Detection (VRD), where we design end-to-end learnable architectures for representing relevant interactions between modalities. Through extensive experiments, we show that BLOCK compares favorably with respect to state-of-the-art multimodal fusion models for both VQA and VRD tasks. Our code is available at https://github.com/Cadene/block.bootstrap.pytorch.


title: Euler flows and singular geometric structures

abstract: Tichler proved that a manifold admitting a smooth closed one-form fibers over a circle. More generally a manifold admitting $k$ independent closed one-forms fibers over a torus $T^k$. In this article we explain a version of this construction for manifolds with boundary using the techniques of $b$-calculus. We explore new applications of this idea to Fluid Dynamics and more concretely in the study of stationary solutions of the Euler equations. In the study of Euler flows on manifolds, two dichotomic situations appear. For the first one, in which the Bernoulli function is not constant, we provide a new proof of Arnold's structure theorem and describe $b$-symplectic structures on some of the singular sets of the Bernoulli function. When the Bernoulli function is constant, a correspondence between contact structures with singularities and what we call $b$-Beltrami fields is established, thus mimicking the classical correspondence between Beltrami fields and contact structures. These results provide a new technique to analyze the geometry of steady fluid flows on non-compact manifolds with cylindrical ends.


title: Your Gameplay Says It All: Modelling Motivation in Tom Clancy's The Division

abstract: Is it possible to predict the motivation of players just by observing their gameplay data? Even if so, how should we measure motivation in the first place? To address the above questions, on the one end, we collect a large dataset of gameplay data from players of the popular game Tom Clancy's The Division. On the other end, we ask them to report their levels of competence, autonomy, relatedness and presence using the Ubisoft Perceived Experience Questionnaire. After processing the survey responses in an ordinal fashion we employ preference learning methods based on support vector machines to infer the mapping between gameplay and the reported four motivation factors. Our key findings suggest that gameplay features are strong predictors of player motivation as the best obtained models reach accuracies of near certainty, from 92% up to 94% on unseen players.


title: Remarks on the influence of quantum vacuum fluctuations over a charged test particle near a conducting wall

abstract: Quantum vacuum fluctuations of the electromagnetic field in empty space seem not to produce observable effects over the motion of a charged test particle. However, when a change in the background vacuum state is implemented, as for instance when a conducting boundary is introduced, dispersions of the particle velocity may occur. As a consequence, besides the existence of classical effects due to the interaction between particle and boundary, there will be a quantum contribution to the motion of the particle whose magnitude depends on how fast the transition between the different vacuum states occurs. Here this issue is revisited and a smooth transition with a controllable switching time between the vacuum states of the system is implemented. Dispersions of the particle velocity in both, zero and finite temperature regimes are examined. More than just generalizing previous results for specific configurations, new effects are unveiled. Particularly, it is shown that the well known vacuum dominance reported to occur arbitrarily near the wall is a consequence of assumed idealizations. The use of a controllable switching enable us to conclude that thermal effects can be as important as or even stronger than vacuum effects arbitrarily near the wall. Additionally, the residual effect predicted to occur in the late time regime was here shown to be linked to the duration of the transition. In this sense, such effect is understood to be a sort of particle energy exchanging due to the vacuum state transition. Furthermore, in certain arrangements a sort of cooling effect over the motion of the particle can occur, i.e., the kinetic energy of the particle is lessen by a certain amount due to subvacuum quantum fluctuations.


title: Coupling phenomena and collective effects in resonant meta-molecules supporting plasmonic and magnetic functionalities: a review

abstract: We review both the fundamental aspects and the applications of functional magneto-optic and opto-magnetic metamaterials displaying collective and coupling effects on the nanoscale, where the concepts of optics and magnetism merge to produce unconventional phenomena. The use of magnetic materials instead of the usual noble metals allows for an additional degree of freedom for the control of electromagnetic field properties, as well as it allows light to interact with the spins of the electrons and to actively manipulate the magnetic properties of such nanomaterials. In this context, we explore the concepts of near-field coupling of plasmon modes in magnetic meta-molecules, as well as the effect of excitation of surface lattice resonances in magneto-plasmonic crystals. Moreover, we discuss how these coupling effects can be exploited to artificially enhance optical magnetism in plasmonic meta-molecules and crystals. Finally, we highlight some of the present challenges and provide a perspective on future directions of the research towards photon-driven fast and efficient nanotechnologies bridging magnetism and optics beyond current limits.


title: Perils and Challenges of Social Media and Election Manipulation Analysis: The 2018 US Midterms

abstract: One of the hallmarks of a free and fair society is the ability to conduct a peaceful and seamless transfer of power from one leader to another. Democratically, this is measured in a citizen population's trust in the electoral system of choosing a representative government. In view of the well documented issues of the 2016 US Presidential election, we conducted an in-depth analysis of the 2018 US Midterm elections looking specifically for voter fraud or suppression. The Midterm election occurs in the middle of a 4 year presidential term. For the 2018 midterms, 35 senators and all the 435 seats in the House of Representatives were up for re-election, thus, every congressional district and practically every state had a federal election. In order to collect election related tweets, we analyzed Twitter during the month prior to, and the two weeks following, the November 6, 2018 election day. In a targeted analysis to detect statistical anomalies or election interference, we identified several biases that can lead to wrong conclusions. Specifically, we looked for divergence between actual voting outcomes and instances of the #ivoted hashtag on the election day. This analysis highlighted three states of concern: New York, California, and Texas. We repeated our analysis discarding malicious accounts, such as social bots. Upon further inspection and against a backdrop of collected general election-related tweets, we identified some confounding factors, such as population bias, or bot and political ideology inference, that can lead to false conclusions. We conclude by providing an in-depth discussion of the perils and challenges of using social media data to explore questions about election manipulation.


title: Euclid Preparation IV. Impact of undetected galaxies on weak-lensing shear measurements

abstract: In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.


title: Gaussian Conditional Random Fields for Classification

abstract: Gaussian conditional random fields (GCRF) are a well-known used structured model for continuous outputs that uses multiple unstructured predictors to form its features and at the same time exploits dependence structure among outputs, which is provided by a similarity measure. In this paper, a Gaussian conditional random fields model for structured binary classification (GCRFBC) is proposed. The model is applicable to classification problems with undirected graphs, intractable for standard classification CRFs. The model representation of GCRFBC is extended by latent variables which yield some appealing properties. Thanks to the GCRF latent structure, the model becomes tractable, efficient and open to improvements previously applied to GCRF regression models. In addition, the model allows for reduction of noise, that might appear if structures were defined directly between discrete outputs. Additionally, two different forms of the algorithm are presented: GCRFBCb (GCRGBC - Bayesian) and GCRFBCnb (GCRFBC - non Bayesian). The extended method of local variational approximation of sigmoid function is used for solving empirical Bayes in Bayesian GCRFBCb variant, whereas MAP value of latent variables is the basis for learning and inference in the GCRFBCnb variant. The inference in GCRFBCb is solved by Newton-Cotes formulas for one-dimensional integration. Both models are evaluated on synthetic data and real-world data. It was shown that both models achieve better prediction performance than unstructured predictors. Furthermore, computational and memory complexity is evaluated. Advantages and disadvantages of the proposed GCRFBCb and GCRFBCnb are discussed in detail.


title: Power Loading based on Portfolio Theory for Densified Millimeter-Wave Small-Cell Communications

abstract: We experimentally demonstrate a novel scheme of power loading based on portfolio theory for millimeter-wave small-cell densification. By exploiting the statistical characteristics of interference, this approach improves the average throughput by 91% and reduces the variance.


title: Topography of (exo)planets

abstract: Current technology is not able to map the topography of rocky exoplanets, simply because the objects are too faint and far away to resolve them. Nevertheless, indirect effect of topography should be soon observable thanks to photometry techniques, and the possibility of detecting specular reflections. In addition, topography may have a strong effect on Earth-like exoplanet climates because oceans and mountains affect the distribution of clouds \citep{Houze2012}. Also topography is critical for evaluating surface habitability \citep{Dohm2015}. We propose here a general statistical theory to describe and generate realistic synthetic topographies of rocky exoplanetary bodies. In the solar system, we have examined the best-known bodies: the Earth, Moon, Mars and Mercury. It turns out that despite their differences, they all can be described by multifractral statistics, although with different parameters. Assuming that this property is universal, we propose here a model to simulate 2D spherical random field that mimics a rocky planetary body in a stellar system. We also propose to apply this model to estimate the statistics of oceans and continents to help to better assess the habitability of distant worlds.


title: Weak antilocalization in quasi-two-dimensional electronic states of epitaxial LuSb thin films

abstract: Observation of large non-saturating magnetoresistance in rare-earth monopnictides has raised enormous interest in understanding the role of its electronic structure. Here, by a combination of molecular-beam epitaxy, low-temperature transport, angle-resolved photoemssion spectroscopy, and hybrid density functional theory we have unveiled the bandstructure of LuSb, where electron-hole compensation is identified as a mechanism responsible for large magnetoresistance in this topologically trivial compound. In contrast to bulk single crystal analogues, quasi-two-dimensional behavior is observed in our thin films for both electron and holelike carriers, indicative of dimensional confinement of the electronic states. Introduction of defects through growth parameter tuning results in the appearance of quantum interference effects at low temperatures, which has allowed us to identify the dominant inelastic scattering processes and elucidate the role of spin-orbit coupling. Our findings open up new possibilities of band structure engineering and control of transport properties in rare-earth monopnictides via epitaxial synthesis.


title: Coarse-graining the dynamics of immersed and driven fiber assemblies

abstract: An important class of fluid-structure problems involve the dynamics of ordered arrays of immersed, flexible fibers. While specialized numerical methods have been developed to study fluid-fiber systems, they become infeasible when there are many, rather than a few, fibers present, nor do these methods lend themselves to analytical calculation. Here, we introduce a coarse-grained continuum model, based on local-slender body theory, for elastic fibers immersed in a viscous Newtonian fluid. It takes the form of an anisotropic Brinkman equation whose skeletal drag is coupled to elastic forces. This model has two significant benefits: (1) the density effects of the fibers in a suspension become analytically manifest, and (2) it allows for the rapid simulation of dense suspensions of fibers in regimes inaccessible to standard methods. As a first validation, without fitting parameters, we achieve very reasonable agreement with 3D Immersed Boundary simulations of a bed of anchored fibers bent by a shear flow. Secondly, we characterize the effect of density on the relaxation time of fiber beds under oscillatory shear, and find close agreement to results from full numerical simulations. We then study buckling instabilities in beds of fibers, using our model both numerically and analytically to understand the role of fiber density and the structure of buckling transitions. We next apply our model to study the flow-induced bending of inclined fibers in a channel, as has been recently studied as a flow rectifier, examining the nature of the internal flows within the bed, and the emergence of inhomogeneous permeability. Finally, we extend the method to study a simple model of metachronal waves on beds of actuated fibers, as a model for ciliary beds. Our simulations reproduce qualitatively the pumping action of coordinated waves of compression through the bed.


title: On the monotonicity of copula-based conditional distributions

abstract: In this paper, we find necessary and sufficient conditions so that copula-based conditional distributions of a response variable with respect to covariates, are ordered with respect to the simple stochastic order introduced by Lehmann. These conditions do not depend on the marginal distributions of the random variables. As a result, we have conditions to ensure that the conditional mean and conditional quantiles are monotonic with respect to the covariates.


title: Shape Analysis, Lebesgue Integration and Absolute Continuity Connections

abstract: As shape analysis of the form presented in Srivastava and Klassen's textbook 'Functional and Shape Data Analysis' is intricately related to Lebesgue integration and absolute continuity, it is advantageous to have a good grasp of the latter two notions. Accordingly, in these notes we review basic concepts and results about Lebesgue integration and absolute continuity. In particular, we review fundamental results connecting them to each other and to the kind of shape analysis, or more generally, functional data analysis presented in the aforementioned textbook, in the process shedding light on important aspects of all three notions. Many well-known results, especially most results about Lebesgue integration and some results about absolute continuity, are presented without proofs. However, a good number of results about absolute continuity and most results about functional data and shape analysis are presented with proofs. Actually, most missing proofs can be found in Royden's 'Real Analysis' and Rudin's 'Principles of Mathematical Analysis' as it is on these classic textbooks and Srivastava and Klassen's textbook that a good portion of these notes are based. However, if the proof of a result does not appear in the aforementioned textbooks, nor in some other known publication, or if all by itself it could be of value to the reader, an effort has been made to present it accordingly.


title: Base Station Distance Adaptive LEACH

abstract: For some applications, we need to deploy a network of sensors in working field to sense the environment and send collected data to a base station for processing; these sensors depend on non-rechargeable batteries, so the routing protocols for such network of sensors need to be efficient. LEACH is one of these protocols which is a hierarchical routing protocol and helps in saving energy in wireless sensor networks. Enhanced LEACH depends on a mathematical model to calculate the estimated average energy in each round. consequently, utilizing the node remaining energy to ensure rotating cluster head role over all the nodes. It also depends on a mathematical model to calculate base station distance from work field whereas LEACH does not take into its account any consideration for remaining energy of node. In this paper, we enhance LEACH (work efficiency in homogeneous networks) to adapt with base-station distance, thus more energy saving for certain distances from base-station. The obtained simulation results show that enhanced LEACH saves energy better than LEACH and increase network stability and reliability when base-station is inside working field and consume the same energy as LEACH when base-station is outside work field.


title: Quantifiable & Comparable Evaluations of Cyber Defensive Capabilities: A Survey & Novel, Unified Approach

abstract: Metrics and frameworks to quantifiably assess security measures have arisen from needs of three distinct research communities - statistical measures from the intrusion detection and prevention literature, evaluation of cyber exercises, e.g.,red-team and capture-the-flag competitions, and economic analyses addressing cost-versus-security tradeoffs. In this paper we provide two primary contributions to the security evaluation literature - a representative survey, and a novel framework for evaluating security that is flexible, applicable to all three use cases, and readily interpretable. In our survey of the literature we identify the distinct themes from each community's evaluation procedures side by side and flesh out the drawbacks and benefits of each. The evaluation framework we propose includes comprehensively modeling the resource, labor, and attack costs in dollars incurred based on expected resource usage, accuracy metrics, and time. This framework provides a unified approach in that it incorporates the accuracy and performance metrics, which dominate intrusion detection evaluation, the time to detection and impact to data and resources of an attack, favored by educational competitions' metrics, and the monetary cost of many essential security components used in financial analysis. Moreover, it is flexible enough to accommodate each use case, easily interpretable and comparable, and comprehensive in terms of costs considered.Finally, we provide two examples of the framework applied to real-world use cases. Overall, we provide a survey and a grounded, flexible framework with multiple concrete examples for evaluating security which can address the needs of three currently distinct communities.


title: Comparison and Experimental Validation of Predictive Models for Soft, Fiber-Reinforced Actuators

abstract: Successful soft robot modeling approaches appearing in recent literature have been based on a variety of distinct theories, including traditional robotic theory, continuum mechanics, and machine learning. Though specific modeling techniques have been developed for and validated against already realized systems, their strengths and weaknesses have not been explicitly compared against each other. In this paper, we show how three distinct model structures ---a lumped-parameter model, a continuum mechanical model, and a neural network--- compare in capturing the gross trends and specific features of the force generation of soft robotic actuators. In particular, we study models for Fiber Reinforced Elastomeric Enclosures (FREEs), which are a popular choice of soft actuator and that are used in several soft articulated systems, including soft manipulators, exoskeletons, grippers, and locomoting soft robots. We generated benchmark data by testing eight FREE samples that spanned broad design and kinematic spaces and compared the models on their ability to predict the loading-deformation relationships of these samples. This comparison shows the predictive capabilities of each model on individual actuators and each model's generalizability across the design space. While the neural net achieved the highest peak performance, the first principles-based models generalized best across all actuator design parameters tested. The results highlight the essential roles of mathematical structure and experimental parameter determination in building high-performing, generalizable soft actuator models with varying effort invested in system identification.


title: Chaos-preserving reduction of the spin-flip model for VCSELs: failure of the adiabatic elimination of the spin-population difference

abstract: When studying the dynamics of Vertical-Cavity Surface-Emitting Lasers, and their polarization properties, the spin-flip model appears to be the simplest model qualitatively reproducing all dynamical features that have been observed experimentally. Nonetheless, because of the fast time-scale of the spin-relaxation processes, the specific role and the importance of the spin-population difference - which is one of the specific feature of the spin-flip model - has been continuously questioned. In fact, the debate regarding the possible adiabatic elimination of the spin-population remains fully open. In this paper, our goal is to bring new light into this issue by demonstrating that this variable is essential to preserve the most complex dynamical features predicted by the spin-flip model, such as polarization chaos, and, therefore, needs to be conserved. To do so, we first perform a detailed analysis, focusing on the chaotic dynamics, to determine the minimal embedding dimension for the spin-flip model. As the latter confirms that a reduction of the model could be envisaged, we then consider the adiabatic elimination of the spin-population difference to highlight and explain its failure to reproduce essential dynamical features obtained in the original model.


title: The Dynamics of Canalizing Boolean Networks

abstract: Boolean networks are a popular modeling framework in computational biology to capture the dynamics of molecular networks, such as gene regulatory networks. It has been observed that many published models of such networks are defined by regulatory rules driving the dynamics that have certain so-called canalizing properties. In this paper, we investigate the dynamics of a random Boolean network with such properties using analytical methods and simulations. From our simulations, we observe that Boolean networks with higher canalizing depth have generally fewer attractors, the attractors are smaller, and the basins are larger, with implications for the stability and robustness of the models. These properties are relevant to many biological applications. Moreover, our results show that, from the standpoint of the attractor structure, high canalizing depth, compared to relatively small positive canalizing depth, has a very modest impact on dynamics. Motivated by these observations, we conduct mathematical study of the attractor structure of a random Boolean network of canalizing depth one (i.e., the smallest positive depth). For every positive integer $\ell$, we give an explicit formula for the limit of the expected number of attractors of length $\ell$ in an $n$-state random Boolean network as $n$ goes to infinity.


title: Probabilistic Discriminative Learning with Layered Graphical Models

abstract: Probabilistic graphical models are traditionally known for their successes in generative modeling. In this work, we advocate layered graphical models (LGMs) for probabilistic discriminative learning. To this end, we design LGMs in close analogy to neural networks (NNs), that is, they have deep hierarchical structures and convolutional or local connections between layers. Equipped with tensorized truncated variational inference, our LGMs can be efficiently trained via backpropagation on mainstream deep learning frameworks such as PyTorch. To deal with continuous valued inputs, we use a simple yet effective soft-clamping strategy for efficient inference. Through extensive experiments on image classification over MNIST and FashionMNIST datasets, we demonstrate that LGMs are capable of achieving competitive results comparable to NNs of similar architectures, while preserving transparent probabilistic modeling.


title: Tagged-Particle Statistics in Single-File Motion with Random-Acceleration and Langevin Dynamics

abstract: In the simplest model of single-file diffusion, $N$ point particles wander on a segment of the $x$ axis of length $L$, with hard core interactions, which prevent passing, and with overdamped Brownian dynamics, $\lambda\dot{x}=\eta(t)$, where $\eta(t)$ has the form of Gaussian white noise with zero mean. In 1965 Harris showed that in the limit $N\to\infty$, $L\to\infty$ with constant $\rho=N/L$, the mean square displacement of a tagged particle grows subdiffusively, as $t^{1/2}$, for long times. Recently, it has been shown that the proportionality constants of the $t^{1/2}$ law for randomly-distributed initial positions of the particles and for equally-spaced initial positions are not the same, but have ratio $\sqrt{2}$. In this paper we consider point particles on the $x$ axis, which collide elastically, and which move according to (i) random-acceleration dynamics $\ddot{x}=\eta(t)$ and (ii) Langevin dynamics $\ddot{x}+\lambda\dot{x}=\eta(t)$. The mean square displacement and mean-square velocity of a tagged particle are analyzed for both types of dynamics and for random and equally-spaced initial positions and Gaussian-distributed initial velocities. We also study tagged particle statistics, for both types of dynamics, in the spreading of a compact cluster of particles, with all of the particles initially at the origin.


title: Full-sky weak lensing: a nonlinear post-Friedmann treatment

abstract: We present a full-sky derivation of weak lensing observables in the Post-Friedmann (PF) formalism. Weak lensing has the characteristic of mixing small scales and large scales since it is affected by inhomogeneities integrated along the photon trajectory. With the PF formalism, we develop a modelling of lensing observables which encompasses both leading order relativistic effects and effects that are due to the fully non-linear matter distribution at small scales. We derive the reduced shear, convergence and rotation up to order $1/c^4$ in the PF approximation, accounting for scalar, vector and tensor perturbations, as well as galaxies' peculiar velocities. We discuss the various contributions that break the Kaiser-Squires relation between the shear and the convergence at different orders. We pay particular attention to the impact of the frame-dragging vector potential on lensing observables and we discuss potential ways to measure this effect in future lensing surveys.


title: Predicting Toxicity from Gene Expression with Neural Networks

abstract: We train a neural network to predict chemical toxicity based on gene expression data. The input to the network is a full expression profile collected either in vitro from cultured cells or in vivo from live animals. The output is a set of fine grained predictions for the presence of a variety of pathological effects in treated animals. When trained on the Open TG-GATEs database it produces good results, outperforming classical models trained on the same data. This is a promising approach for efficiently screening chemicals for toxic effects, and for more accurately evaluating drug candidates based on preclinical data.


title: Image Reconstruction from Undersampled Confocal Microscopy Data using Multiresolution Based Maximum Entropy Regularization

abstract: We consider the problem of reconstructing 2D images from randomly under-sampled confocal microscopy samples. The well known and widely celebrated total variation regularization, which is the L1 norm of derivatives, turns out to be unsuitable for this problem; it is unable to handle both noise and under-sampling together. This issue is linked with the notion of phase transition phenomenon observed in compressive sensing research, which is essentially the break-down of total variation methods, when sampling density gets lower than certain threshold. The severity of this breakdown is determined by the so-called mutual incoherence between the derivative operators and measurement operator. In our problem, the mutual incoherence is low, and hence the total variation regularization gives serious artifacts in the presence of noise even when the sampling density is not very low. There has been very few attempts in developing regularization methods that perform better than total variation regularization for this problem. We develop a multi-resolution based regularization method that is adaptive to image structure. In our approach, the desired reconstruction is formulated as a series of coarse-to-fine multi-resolution reconstructions; for reconstruction at each level, the regularization is constructed to be adaptive to the image structure, where the information for adaption is obtained from the reconstruction obtained at coarser resolution level. This adaptation is achieved by using maximum entropy principle, where the required adaptive regularization is determined as the maximizer of entropy subject to the information extracted from the coarse reconstruction as constraints. We demonstrate the superiority of the proposed regularization method over existing ones using several reconstruction examples.


title: Contactless Cardiac Arrest Detection Using Smart Devices

abstract: Out-of-hospital cardiac arrest (OHCA) is a leading cause of death worldwide. Rapid diagnosis and initiation of cardiopulmonary resuscitation (CPR) is the cornerstone of therapy for victims of cardiac arrest. Yet a significant fraction of cardiac arrest victims have no chance of survival because they experience an unwitnessed event, often in the privacy of their own homes. An under-appreciated diagnostic element of cardiac arrest is the presence of agonal breathing, an audible biomarker and brainstem reflex that arises in the setting of severe hypoxia. Here, we demonstrate that a support vector machine (SVM) can classify agonal breathing instances in real-time within a bedroom environment. Using real-world labeled 9-1-1 audio of cardiac arrests, we train the SVM to accurately classify agonal breathing instances. We obtain an area under the curve (AUC) of 0.998 and an operating point with an overall sensitivity and specificity of 97.03% (95% CI: 96.62 -- 97.41%) and 98.20% (95% CI: 97.87 -- 98.49%). We achieve a false positive rate between 0% -- 0.10% over 82 hours (117,895 audio segments) of polysomnographic sleep lab data that includes snoring, hypopnea, central and obstructive sleep apnea events. We demonstrate the effectiveness of our contactless system in identifying real-world cardiac arrest-associated agonal breathing instances and successfully evaluate our classifier using commodity smart devices (Amazon Echo and Apple iPhone).


title: Superconducting quantum refrigerator: Breaking and rejoining Cooper pairs with magnetic field cycles

abstract: We propose a solid state refrigeration technique based on repeated adiabatic magnetization/demagnetization cycles of a superconductor which acts as the working substance. The gradual cooling down of a substrate (normal metal) in contact with the working substance is demonstrated for different initial temperatures of the substrate. Excess heat is given to a hot large-gap superconductor. The on-chip refrigerator works in a cyclic manner because of an effective thermal switching mechanism: Heat transport between N/N versus N/S junctions is asymmetric because of the appearance of the energy gap. This switch permits selective cooling of the metal. We find that this refrigeration technique can cool down a 0.3cm$^{3}$ block of Cu by almost two orders of magnitude starting from 200mK, and down to about 1mK starting from the base temperature of a dilution fridge (10mK). The corresponding cooling power for a 1cm$\times$1cm interface are 25 nW and 0.06 nW respectively, which scales with the area of the interface.


title: Completeness of infinitary heterogeneous logic

abstract: Given a regular cardinal $\kappa$ such that $\kappa^{<\kappa}=\kappa$ (e.g., if the Generalized Continuum Hypothesis holds), we develop a proof system for classical infinitary logic that includes heterogeneous quantification (i.e., infinite alternate sequences of quantifiers) within the language $\mathcal{L}_{\kappa^+, \kappa}$, where there are conjunctions and disjunctions of at most $\kappa$ any formulas and quantification (including the heterogeneous one) is applied to less than $\kappa$ many variables. This type of quantification is interpreted in $\mathcal{S}et$ using the usual second-order formulation in terms of strategies for games, and the axioms are based on a stronger variant of the axiom of determinacy for game semantics. Within this axiom system we prove the soundness and completeness theorem with respect to a class of set-valued structures that we call well-determined. Although this class is more restrictive than the class of determined structures in Takeuti's determinate logic, the completeness theorem works in our case for a wider variety of formulas of $\mathcal{L}_{\kappa^+, \kappa}$, and the category of well-determined models of heterogeneous theories is accessible. Our system is formulated within the sequent style of categorical logic and we do not need to impose any specific requirements on the proof trees, disregarding thus the eigenvariable conditions needed in Takeuti's system. We also investigate intuitionistic systems with heterogeneous quantifiers and prove analogously a completeness theorem with respect to well-determined structures in categories in general, in $\kappa$-Grothendieck toposes in particular, and, when $\kappa^{<\kappa}=\kappa$, also in Kripke models. Finally, we consider an extension of our system in which heterogeneous quantification with bounded quantifiers is expressible, and extend our completeness results to that case.


title: Unified description of pairing and quarteting correlations within the particle-hole-boson approach

abstract: We study the description of single-species and isovector pairing correlations in the framework of the projected-BCS (PBCS) and the Quartet Condensation Model (QCM) from a particle-hole perspective and we introduce the representation of the QCM quartet condensate state in terms of particle-hole excitations with respect to the Hartree-Fock state. We also present a new bosonic approximation for both PBCS and QCM. In each case, the starting point is the reformulation of the pair/quartet condensate state in terms of particle-hole excitations with respect to the Hartree-Fock state. The main simplification of our approach is the assumption that the pair operators corresponding to both particle and hole states obey bosonic commutation relations. This simplifies tremendously the computations and allows for an analytic derivation of the averaged Hamiltonian on the condenstate state as a function of the mixing amplitudes. We study both the pure bosonic approach and the renormalized version, and compare the particle-hole bosonic version to the naive prescription of applying the boson approximation directly to the original condensate state. We compare the fermionic and the renormalized particle-hole bosonic approach in the case of a picket fence model of doubly degenerate states and in a realistic shell model space with an effective interaction for the $N=Z$ nuclei above $^{100}$Sn.


title: Separable Resolution-of-the-Identity with All-Electron Gaussian Bases: Application to Cubic-scaling RPA

abstract: We explore a separable resolution-of-the-identity formalism built on quadratures over limited sets of real-space points designed for all-electron calculations. Our implementation preserves in particular the use of common atomic orbitals and their related auxiliary basis sets. The set up of the present density fitting scheme, i.e. the calculation of the system specific quadrature weights, scales cubically with respect to the system size. Extensive accuracy tests are presented for the Fock exchange and MP2 correlation energies. We finally demonstrate random phase approximation (RPA) correlation energy calculations with a scaling that is cubic in terms of operations, quadratic in memory, with a small crossover with respect to our standard RI-RPA implementation.


title: Bayesian active learning for optimization and uncertainty quantification in protein docking

abstract: Motivation: Ab initio protein docking represents a major challenge for optimizing a noisy and costly "black box"-like function in a high-dimensional space. Despite progress in this field, there is no docking method available for rigorous uncertainty quantification (UQ) of its solution quality (e.g. interface RMSD or iRMSD). Results: We introduce a novel algorithm, Bayesian Active Learning (BAL), for optimization and UQ of such black-box functions and flexible protein docking. BAL directly models the posterior distribution of the global optimum (or native structures for protein docking) with active sampling and posterior estimation iteratively feeding each other. Furthermore, we use complex normal modes to represent a homogeneous Euclidean conformation space suitable for high-dimension optimization and construct funnel-like energy models for encounter complexes. Over a protein docking benchmark set and a CAPRI set including homology docking, we establish that BAL significantly improve against both starting points by rigid docking and refinements by particle swarm optimization, providing for one third targets a top-3 near-native prediction. BAL also generates tight confidence intervals with half range around 25% of iRMSD and confidence level at 85%. Its estimated probability of a prediction being native or not achieves binary classification AUROC at 0.93 and AUPRC over 0.60 (compared to 0.14 by chance); and also found to help ranking predictions. To the best of our knowledge, this study represents the first uncertainty quantification solution for protein docking, with theoretical rigor and comprehensive assessment. Source codes are available at https://github.com/Shen-Lab/BAL.


title: Carleman estimates with sharp weights and boundary observability for wave operators with critically singular potentials

abstract: We establish a new family of Carleman inequalities for wave operators on cylindrical spacetime domains containing a potential that is critically singular, diverging as an inverse square on all the boundary of the domain. These estimates are sharp in the sense that they capture both the natural boundary conditions and the natural $H^1$-energy. The proof is based around three key ingredients: the choice of a novel Carleman weight with rather singular derivatives on the boundary, a generalization of the classical Morawetz inequality that allows for inverse-square singularities, and the systematic use of derivative operations adapted to the potential. As an application of these estimates, we prove a boundary observability property for the associated wave equations.


title: Conformal deformations preserving the Finslerian $R$-Einstein criterion

abstract: Given a Finslerian metric $F$ on a $C^4$-manifold, conformal deformations of $F$ preserving the $R$-Einstein criterion are presented. In particular, locally conformal invariance between two Finslerian $R$-Einstein metrics is characterized.


title: On some spectral properties of pseudo-differential operators on T

abstract: In this paper we use Riesz spectral Theory and Gershgorin Theory to obtain explicit information concerning the spectrum of pseudo-differential operators defined on the unit circle $\mathbb{T} := \mathbb{R}/ 2 \pi \mathbb{ Z}$. For symbols in the H\"ormander class $S^m_{1 , 0} (\mathbb{T} \times \mathbb{Z})$, we provide a sufficient and necessary condition to ensure that the corresponding pseudo-differential operator is a Riesz operator in $L^p (\mathbb{T})$, $1< p < \infty$, extending in this way compact operators characterisation and Ghoberg's lemma to $L^p (\mathbb{T})$. We provide an example of a non-compact Riesz pseudo-differential operator in $L^p (\mathbb{T})$, $1<p<2$. Also, for pseudo-differential operators with symbol satisfying some integrability condition, it is defined its associated matrix in terms of the Fourier coefficients of the symbol, and this matrix is used to give necessary and sufficient conditions for $L^2$-boundedness without assuming any regularity on the symbol, and to locate the spectrum of some operators.


title: Optimal mini-batch and step sizes for SAGA

abstract: Recently it has been shown that the step sizes of a family of variance reduced gradient methods called the JacSketch methods depend on the expected smoothness constant. In particular, if this expected smoothness constant could be calculated a priori, then one could safely set much larger step sizes which would result in a much faster convergence rate. We fill in this gap, and provide simple closed form expressions for the expected smoothness constant and careful numerical experiments verifying these bounds. Using these bounds, and since the SAGA algorithm is part of this JacSketch family, we suggest a new standard practice for setting the step sizes and mini-batch size for SAGA that are competitive with a numerical grid search. Furthermore, we can now show that the total complexity of the SAGA algorithm decreases linearly in the mini-batch size up to a pre-defined value: the optimal mini-batch size. This is a rare result in the stochastic variance reduced literature, only previously shown for the Katyusha algorithm. Finally we conjecture that this is the case for many other stochastic variance reduced methods and that our bounds and analysis of the expected smoothness constant is key to extending these results.


title: Novel topological beam-splitting in photonic crystals

abstract: We create a passive wave splitter, created purely by geometry, to engineer three-way beam splitting in electromagnetism in transverse electric polarisation. We do so by considering arrangements of Indium Phosphide dielectric pillars in air, in particular we place several inclusions within a cell that is then extended periodically upon a square lattice. Hexagonal lattice structures more commonly used in topological valleytronics but, as we discuss, three-way splitting is only possible using a square, or rectangular, lattice. To achieve splitting and transport around a sharp bend we use accidental, and not symmetry-induced, Dirac cones. Within each cell pillars are either arranged around a triangle or square; we demonstrate the mechanism of splitting and why it does not occur for one of the cases. The theory is developed and full scattering simulations demonstrate the effectiveness of the proposed designs.


title: muSR and Magnetometry Study of the Type-I Superconductor BeAu

abstract: We present muon spin rotation and relaxation (muSR) measurements as well as demagnetising field corrected magnetisation measurements on polycrystalline samples of the noncentrosymmetric superconductor BeAu. From muSR measurements in a transverse field, we determine that BeAu is a type-I superconductor with Hc = 256 Oe, amending the previous understanding of the compound as a type-II superconductor. To account for demagnetising effects in magnetisation measurements, we produce an ellipsoidal sample, for which a demagnetisation factor can be calculated. After correcting for demagnetising effects, our magnetisation results are in agreement with our muSR measurements. Using both types of measurements we construct a phase diagram from T = 30 mK to Tc = 3.25 K. We then study the effect of hydrostatic pressure and find that 450 MPa decreases Tc by 34 mK, comparable to the change seen in type-I elemental superconductors Sn, In and Ta, suggesting BeAu is far from a quantum critical point accessible by the application of pressure.


title: Globally regular solutions to dyonic anti-de Sitter $\mathfrak{su}(\infty)$ Einstein-Yang-Mills theory -- Existence and characterising charges

abstract: In this work, we find new static, spherically symmetric, dyonic, globally regular exact solutions to $\mathfrak{su}(\infty)$ Einstein-Yang-Mills theory with a negative cosmological constant $\Lambda$, in the regime that $|\Lambda|$ is very large. In this regime, we also prove that dyonic globally regular solutions may be uniquely characterised by a countably infinite set of effective global charges; and that dyon solutions may be distinguished from dyonic black hole solutions to the same field equations by their ADM masses. These solutions have potential modelling applications for certain exotic gravitational objects.


title: Trading beams for bandwidth: Imaging with randomized beamforming

abstract: We study the problem of actively imaging a range-limited far-field scene using an antenna array. We describe how the range limit imposes structure in the measurements across multiple wavelengths. This structure allows us to introduce a novel trade-off: the number of spatial array measurements (i.e., beams that have to be formed) can be reduced significantly lower than the number array elements if the scene is illuminated with a broadband source. To take advantage of this trade-off, we use a small number of "generic" linear combinations of the array outputs, instead of the phase offsets used in conventional beamforming. We provide theoretical justification for the proposed trade-off without making any strong structural assumptions on the target scene (such as sparsity) except that it is range limited. In proving our theoretical results, we take inspiration from the sketching literature. We also provide simulation results to establish the merit of the proposed signal acquisition strategy. Our proposed method results in a reduction in the number of required spatial measurements in an array imaging system and hence can directly impact their speed and cost of operation.


title: Phase diagrams of a 2D Ising spin-pseudospin model

abstract: We consider the competition of magnetic and charge ordering in high-Tc cuprates within the framework of the simplified static 2D spin-pseudospin model. This model is equivalent to the 2D dilute antiferromagnetic (AFM) Ising model with charged impurities. We present the mean-field results for the system under study and make a brief comparison with classical Monte Carlo (MC) calculations. Numerical simulations show that the cases of strong exchange and strong charge correlation differ qualitatively. For a strong exchange, the AFM phase is unstable with respect to the phase separation (PS) into the charge and spin subsystems, which behave like immiscible quantum liquids. An analytical expression was obtained for the PS temperature.


title: Fundamental Spin Interactions Underlying the Magnetic Anisotropy in the Kitaev Ferromagnet CrI$_3$

abstract: We lay the foundation for determining the microscopic spin interactions in two-dimensional (2D) ferromagnets by combining angle-dependent ferromagnetic resonance (FMR) experiments on high quality CrI$_3$ single crystals with theoretical modeling based on symmetries. We discover that the Kitaev interaction is the strongest in this material with $K \sim -5.2$ meV, 25 times larger than the Heisenberg exchange $J \sim -0.2$ meV, and responsible for opening the $\sim$5 meV gap at the Dirac points in the spin-wave dispersion. Furthermore, we find that the symmetric off-diagonal anisotropy $\Gamma \sim -67.5$ $\mu$eV, though small, is crucial for opening a $\sim$0.3 meV gap in the magnon spectrum at the zone center and stabilizing ferromagnetism in the 2D limit. The high resolution of the FMR data further reveals a $\mu$eV-scale quadrupolar contribution to the $S=3/2$ magnetism. Our identification of the underlying exchange anisotropies opens paths toward 2D ferromagnets with higher $T_\text{C}$ as well as magnetically frustrated quantum spin liquids based on Kitaev physics.


title: From genome to phenome: Predicting multiple cancer phenotypes based on somatic genomic alterations via the genomic impact transformer

abstract: Cancers are mainly caused by somatic genomic alterations (SGAs) that perturb cellular signaling systems and eventually activate oncogenic processes. Therefore, understanding the functional impact of SGAs is a fundamental task in cancer biology and precision oncology. Here, we present a deep neural network model with encoder-decoder architecture, referred to as genomic impact transformer (GIT), to infer the functional impact of SGAs on cellular signaling systems through modeling the statistical relationships between SGA events and differentially expressed genes (DEGs) in tumors. The model utilizes a multi-head self-attention mechanism to identify SGAs that likely cause DEGs, or in other words, differentiating potential driver SGAs from passenger ones in a tumor. GIT model learns a vector (gene embedding) as an abstract representation of functional impact for each SGA-affected gene. Given SGAs of a tumor, the model can instantiate the states of the hidden layer, providing an abstract representation (tumor embedding) reflecting characteristics of perturbed molecular/cellular processes in the tumor, which in turn can be used to predict multiple phenotypes. We apply the GIT model to 4,468 tumors profiled by The Cancer Genome Atlas (TCGA) project. The attention mechanism enables the model to better capture the statistical relationship between SGAs and DEGs than conventional methods, and distinguishes cancer drivers from passengers. The learned gene embeddings capture the functional similarity of SGAs perturbing common pathways. The tumor embeddings are shown to be useful for tumor status representation, and phenotype prediction including patient survival time and drug response of cancer cell lines.


title: Inaccuracy of Spatial Derivatives in Simulations of Supersonic Turbulence

abstract: We examine the accuracy of spatial derivatives computed from numerical simulations of supersonic turbulence. Two sets of simulations, carried out using a finite-volume code that evolves the hydrodynamic equations with an approximate Riemann solver and a finite-difference code that solves the Navier-Stokes equations, are tested against a number of criteria based on the continuity equation, including exact results at statistically steady state. We find that the spatial derivatives in the Navier-Stokes runs are accurate and satisfy all the criteria. In particular, they satisfy our exact results that the conditional mean velocity divergence, $\langle \nabla \cdot {\bs u}|s\rangle$, where $s$ is the logarithm of density, and the conditional mean of the advection of $s$, $\langle {\bs u} \cdot \nabla s|s\rangle$, vanish at steady state for all density values, $s$. On the other hand, the Riemann solver simulations fail all the tests that require accurate evaluation of spatial derivatives, resulting in apparent violation of the continuity equation, even if the solver enforces mass conservation. In particular, analysis of the Riemann simulations may lead to the incorrect conclusion that the $\pdv$ work tends to preferentially convert kinetic energy into thermal energy, inconsistent with the exact result that the energy exchange by $\pdv$ work is symmetric in barotropic supersonic turbulence at steady state. The inaccuracy of spatial derivatives is a general problem in the post-processing of simulations of supersonic turbulence with Riemann solvers. Solutions from such simulations must be used with caution in post-processing studies concerning the spatial gradients.


title: Minimax Testing of Identity to a Reference Ergodic Markov Chain

abstract: We exhibit an efficient procedure for testing, based on a single long state sequence, whether an unknown Markov chain is identical to or $\varepsilon$-far from a given reference chain. We obtain nearly matching (up to logarithmic factors) upper and lower sample complexity bounds for our notion of distance, which is based on total variation. Perhaps surprisingly, we discover that the sample complexity depends solely on the properties of the known reference chain and does not involve the unknown chain at all, which is not even assumed to be ergodic.


title: Fibrations, unique path lifting, and continuous monodromy

abstract: Given a path-connected space $X$ and $H\leq\pi_1(X,x_0)$, there is essentially only one construction of a map $p_H:(\widetilde{X}_H,\widetilde{x}_0)\rightarrow(X,x_0)$ with connected and locally path-connected domain that can possibly have the following two properties: $(p_{H})_{\#}\pi_1(\widetilde{X}_H,\widetilde{x}_0)=H$ and $p_H$ has the unique lifting property. $\widetilde{X}_H$ consists of equivalence classes of paths starting at $x_0$, appropriately topologized, and $p_H$ is the endpoint projection. For $p_H$ to have these two properties, $T_1$ fibers are necessary and unique path lifting is sufficient. However, $p_H$ always admits the standard lifts of paths. We show that $p_H$ has unique path lifting if it has continuous (standard) monodromies toward a $T_1$ fiber over $x_0$. Assuming, in addition, that $H$ is locally quasinormal (e.g., if $H$ is normal) we show that $X$ is homotopically path Hausdorff relative to $H$. We show that $p_H$ is a fibration if $X$ is locally path connected, $H$ is locally quasinormal, and all (standard) monodromies are continuous.


title: Star Cluster Catalogs for the LEGUS Dwarf Galaxies

abstract: We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.


title: Study of some holomorphic curves in $\C^3$ and their projection into the complex projectve space $\C P^2$

abstract: We study holomorphic curves $f:\C\longrightarrow \C^3$ avoiding four complex hyperplanes and a real subspace of real dimension four or five in $\C^3$. We show that the projection of $f$ into the complex projective space $\C P^2$ is not necessarily constant.


title: Characterizing Input Methods for Human-to-robot Demonstrations

abstract: Human demonstrations are important in a range of robotics applications, and are created with a variety of input methods. However, the design space for these input methods has not been extensively studied. In this paper, focusing on demonstrations of hand-scale object manipulation tasks to robot arms with two-finger grippers, we identify distinct usage paradigms in robotics that utilize human-to-robot demonstrations, extract abstract features that form a design space for input methods, and characterize existing input methods as well as a novel input method that we introduce, the instrumented tongs. We detail the design specifications for our method and present a user study that compares it against three common input methods: free-hand manipulation, kinesthetic guidance, and teleoperation. Study results show that instrumented tongs provide high quality demonstrations and a positive experience for the demonstrator while offering good correspondence to the target robot.


title: Cloud Atlas: High-Contrast Time-Resolved Observations of Planetary-Mass Companions

abstract: Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope.


title: Modelling of Mg II lines in solar prominences

abstract: Observations of the Mg II h and k lines in solar prominences with IRIS reveal a wide range of line shapes from simple non-reversed profiles to typical double-peaked reversed profiles with many other complex line shapes possible. The physical conditions responsible for this variety are not well understood. Our aim is to understand how physical conditions inside a prominence slab influence shapes and properties of emergent Mg II line profiles. We compute the spectrum of Mg II lines using a one-dimensional non-LTE radiative transfer code for two large grids of model atmospheres (isothermal isobaric, and with a transition region). The influence of the plasma parameters on the emergent spectrum is discussed in detail. Our results agree with previous studies. We present several dependencies between observables and prominence parameters which will help with interpretation of observations. A comparison with known limits of observed line parameters suggests that most observed prominences emitting in Mg II h and k lines are cold, low pressure, and optically thick structures. Our results indicate that there are good correlations between the Mg II k line intensities and the intensities of hydrogen lines, as well as the emission measure. One-dimensional non-LTE radiative transfer codes are well-suited to understand the main characteristics of the Mg II h and k line profiles in solar prominences, but more advanced codes will be necessary for detailed comparisons.


title: Learning Triggers for Heterogeneous Treatment Effects

abstract: The causal effect of a treatment can vary from person to person based on their individual characteristics and predispositions. Mining for patterns of individual-level effect differences, a problem known as heterogeneous treatment effect estimation, has many important applications, from precision medicine to recommender systems. In this paper we define and study a variant of this problem in which an individual-level threshold in treatment needs to be reached, in order to trigger an effect. One of the main contributions of our work is that we do not only estimate heterogeneous treatment effects with fixed treatments but can also prescribe individualized treatments. We propose a tree-based learning method to find the heterogeneity in the treatment effects. Our experimental results on multiple datasets show that our approach can learn the triggers better than existing approaches.


title: An open and parallel multiresolution framework using block-based adaptive grids

abstract: A numerical approach for solving evolutionary partial differential equations in two and three space dimensions on block-based adaptive grids is presented. The numerical discretization is based on high-order, central finite-differences and explicit time integration. Grid refinement and coarsening are triggered by multiresolution analysis, i.e. thresholding of wavelet coefficients, which allow controlling the precision of the adaptive approximation of the solution with respect to uniform grid computations. The implementation of the scheme is fully parallel using MPI with a hybrid data structure. Load balancing relies on space filling curves techniques. Validation tests for 2D advection equations allow to assess the precision and performance of the developed code. Computations of the compressible Navier-Stokes equations for a temporally developing 2D mixing layer illustrate the properties of the code for nonlinear multi-scale problems. The code is open source.