title: Rank $N$ Vafa-Witten invariants, modularity and blow-up

abstract: We derive explicit expressions for the generating functions of refined Vafa-Witten invariants $\Omega(\gamma,y)$ of $\mathbb{P}^2$ of arbitrary rank $N$ and for their non-holomorphic modular completions. In the course of derivation we also provide: i) a generalization of the recently found generating functions of $\Omega(\gamma,y)$ and their completions for Hirzebruch and del Pezzo surfaces in the canonical chamber of the moduli space to a generic chamber; ii) a version of the blow-up formula expressed directly in terms of these generating functions and its reformulation in a manifestly modular form.


title: Hydrodynamical instability with noise in the Keplerian accretion discs: Modified Landau equation

abstract: Origin of hydrodynamical instability and turbulence in the Keplerian accretion disc as well as similar laboratory shear flows, e.g. plane Couette flow, is a long standing puzzle. These flows are linearly stable. Here we explore the evolution of perturbation in such flows in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could be result of thermal fluctuations (however small be), Brownian ratchet, grain-fluid interactions and feedback from outflows in astrophysical discs etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is modified Landau equation. We show that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter timescale in the presence of force, leading to a fast transition to turbulence. Interestingly, emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike original Landau equation based solution.


title: Decompositions of Ehrhart $h^*$-polynomials for rational polytopes

abstract: The Ehrhart quasipolynomial of a rational polytope $P$ encodes the number of integer lattice points in dilates of $P$, and the $h^*$-polynomial of $P$ is the numerator of the accompanying generating function. We provide two decomposition formulas for the $h^*$-polynomial of a rational polytope. The first decomposition generalizes a theorem of Betke and McMullen for lattice polytopes. We use our rational Betke--McMullen formula to provide a novel proof of Stanley's Monotonicity Theorem for the $h^*$-polynomial of a rational polytope. The second decomposition generalizes a result of Stapledon, which we use to provide rational extensions of the Stanley and Hibi inequalities satisfied by the coefficients of the $h^*$-polynomial for lattice polytopes. Lastly, we apply our results to rational polytopes containing the origin whose duals are lattice polytopes.


title: Gravitational Faraday and Spin-Hall Effects of Light

abstract: The gravitational Faraday and its dual spin-Hall effects of light arise in space-times of non-zero angular momentum. These effects were studied in stationary, asymptotically flat space-times. Here we study these effects in arbitrary, non-stationary, asymptotically flat space-times. These effects arise due to interaction between light polarisation and space-time angular momentum. As a result of such interaction, the phase velocity of left- and right-handed circularly polarised light becomes different, that results in the gravitational Faraday effect. This difference implies different dynamics of these components, that begin to propagate along different paths\textemdash the gravitational spin-Hall effect of light. Due to this effect, the gravitational field splits a multicomponent beam of unpolarized light and produces polarized gravitational rainbow. The component separation is an accumulative effect observed in long range asymptotics. To study this effect, we construct uniform eikonal expansion and derive dynamical equation describing this effect. To analyse the dynamical equation, we present it in the local space and time decomposition form. The spatial part of the equation presented in the related optical metric is analogous to the dynamical equation of a charged particle moving in magnetic field under influence of the Coriolis force.


title: Charging the Walking U(N)$\times$U(N) Higgs Theory as a Complex CFT

abstract: We apply a semi-classical method to compute the conformal field theory (CFT) data for the U(N)xU(N) non-abelian Higgs theory in four minus epsilon dimensions at its complex fixed point. The theory features more than one coupling and walking dynamics. Given our charge configuration, we identify a family of corresponding operators and compute their scaling dimensions which remarkably agree with available results from conventional perturbation theory validating the use of the state-operator correspondence for a complex CFT.


title: Overcoming Statistical Shortcuts for Open-ended Visual Counting

abstract: Machine learning models tend to over-rely on statistical shortcuts. These spurious correlations between parts of the input and the output labels does not hold in real-world settings. We target this issue on the recent open-ended visual counting task which is well suited to study statistical shortcuts. We aim to develop models that learn a proper mechanism of counting regardless of the output label. First, we propose the Modifying Count Distribution (MCD) protocol, which penalizes models that over-rely on statistical shortcuts. It is based on pairs of training and testing sets that do not follow the same count label distribution such as the odd-even sets. Intuitively, models that have learned a proper mechanism of counting on odd numbers should perform well on even numbers. Secondly, we introduce the Spatial Counting Network (SCN), which is dedicated to visual analysis and counting based on natural language questions. Our model selects relevant image regions, scores them with fusion and self-attention mechanisms, and provides a final counting score. We apply our protocol on the recent dataset, TallyQA, and show superior performances compared to state-of-the-art models. We also demonstrate the ability of our model to select the correct instances to count in the image. Code and datasets are available: https://github.com/cdancette/spatial-counting-network


title: Highly compact neutron stars and screening mechanisms. I. Equilibrium and stability

abstract: Modified theories of gravity that offer viable models for dark energy often rely on mechanisms that screen their effects in high density environments. From this perspective, it would appear that, once solar system constraints are satisfied, these theories would predict a trivial phenomenology for (much denser) neutron stars. In this work we explore the fact that in scalar-tensor theories the scalar degree of freedom does not couple to the mass density alone, but to the trace of the energy-momentum tensor - which can increase and eventually change sign as density and pressure build up in the core of neutron stars -, and investigate whether there could be a partial unscreening of the scalar field inside the most compact stars found in Nature. For this purpose, we construct neutron star solutions with realistic equations of state in theories with screening mechanisms and study their stability under radial perturbations. In particular, we find that stable solutions with a unscreened core can exist in chameleon models, while for the environmentally-dependent dilaton model a wealth of new, scalarized equilibrium solutions are found, some of which can be stable.


title: Phases and Quantum Phase Transitions in an Anisotropic Ferromagnetic Kitaev-Heisenberg-$\ \Gamma$ Magnet

abstract: We study the spin-$1/2$ ferromagnetic Heisenberg-Kitaev-$\Gamma$ model in the anisotropic (Toric code) limit to reveal the nature of the quantum phase transition between the gapped $Z_2$ quantum spin liquid and a spin ordered phase (driven by Heisenberg interactions) as well as a trivial paramagnet (driven by pseudo-dipolar interactions, $\Gamma$). The transitions are obtained by a simultaneous condensation of the Ising electric and magnetic charges-- the fractionalized excitations of the $Z_2$ quantum spin liquid. Both these transitions can be continuous and are examples of deconfined quantum critical points. Crucial to our calculations are the symmetry implementations on the soft electric and magnetic modes that become critical. In particular, we find strong constraints on the structure of the critical theory arising from time reversal and lattice translation symmetries with the latter acting as an anyon permutation symmetry that endows the critical theory with a manifestly self-dual structure. We find that the transition between the quantum spin liquid and the spin-ordered phase belongs to a self-dual modified Abelian Higgs field theory while that between the spin liquid and the trivial paramagnet belongs to a self-dual $Z_2$ gauge theory. We also study the effect of an external Zeeman field to show an interesting similarity between the polarised paramagnet obtained due to the Zeeman field and the trivial paramagnet driven the pseudo-dipolar interactions. Interestingly, both the spin liquid and the spin ordered phases have easily identifiable counterparts in the isotropic limit and the present calculations may shed insights into the corresponding transitions in the material relevant isotropic limit.


title: Gauges in Three-Dimensional Gravity and Holographic Fluids

abstract: Solutions to Einstein's vacuum equations in three dimensions are locally maximally symmetric. They are distinguished by their global properties and their investigation often requires a choice of gauge. Although analyses of this sort have been performed abundantly, several relevant questions remain. These questions include the interplay between the standard Bondi gauge and the Eddington--Finkelstein type of gauge used in the fluid/gravity holographic reconstruction of these spacetimes, as well as the Fefferman--Graham gauge, when available i.e. in anti de Sitter. The goal of the present work is to set up a thorough dictionary for the available descriptions with emphasis on the relativistic or Carrollian holographic fluids, which portray the bulk from the boundary in anti-de Sitter or flat instances. A complete presentation of residual diffeomorphisms with a preliminary study of their algebra accompanies the situations addressed here.


title: Fefferman-Graham and Bondi Gauges in the Fluid/Gravity Correspondence

abstract: In three-dimensional gravity, we discuss the relation between the Fefferman--Graham gauge, the Bondi gauge and the Eddington--Finkelstein type of gauge, often referred to as the derivative expansion, involved in the fluid/gravity correspondence. Starting with a negative cosmological constant, for each gauge, we derive the solution space and the residual gauge diffeomorphisms. We construct explicitly the diffeomorphisms that relate the various gauges, and establish the precise matching of their boundary data. We show that Bondi and Fefferman--Graham gauges are equivalent, while the fluid/gravity derivative expansion, originating from a partial gauge fixing, exhibits an extra unspecified function that encodes the boundary fluid velocity. The Bondi gauge turns out to describe a subspace of the derivative expansion's solution space, featuring a fluid in a specific hydrodynamic frame. We pursue our analysis with the Ricci-flat limit of the Bondi gauge and of the fluid/gravity derivative expansion. The relations between them persist in this limit, which is well-defined and non-trivial. Moreover, the flat limit of the derivative expansion maps to the ultra-relativistic limit on the boundary. This procedure allows to unravel the holographic properties of the Bondi gauge for vanishing cosmological constant, in terms of its boundary Carrollian dual fluid.


title: Quantum time dilation in atomic spectra

abstract: Quantum time dilation occurs when a clock moves in a superposition of relativistic momentum wave packets. We utilize the lifetime of an excited hydrogen-like atom as a clock to demonstrate how quantum time dilation manifests in a spontaneous emission process. The resulting emission rate differs when compared to the emission rate of an atom prepared in a mixture of momentum wave packets at order $v^2/c^2$. This effect is accompanied by a quantum correction to the Doppler shift due to the coherence between momentum wave packets. This quantum Doppler shift affects the spectral line shape at order $v/c$. However, its effect on the decay rate is suppressed when compared to the effect of quantum time dilation. We argue that spectroscopic experiments offer a technologically feasible platform to explore the effects of quantum time dilation.


title: Socially Fair k-Means Clustering

abstract: We show that the popular k-means clustering algorithm (Lloyd's heuristic), used for a variety of scientific data, can result in outcomes that are unfavorable to subgroups of data (e.g., demographic groups). Such biased clusterings can have deleterious implications for human-centric applications such as resource allocation. We present a fair k-means objective and algorithm to choose cluster centers that provide equitable costs for different groups. The algorithm, Fair-Lloyd, is a modification of Lloyd's heuristic for k-means, inheriting its simplicity, efficiency, and stability. In comparison with standard Lloyd's, we find that on benchmark datasets, Fair-Lloyd exhibits unbiased performance by ensuring that all groups have equal costs in the output k-clustering, while incurring a negligible increase in running time, thus making it a viable fair option wherever k-means is currently used.


title: NLO QCD+EW predictions for $tHj$ and $tZj$ production at the LHC

abstract: In this work we calculate the cross sections for the hadroproduction of a single top quark or antiquark in association with a Higgs ($tHj$) or $Z$ boson ($tZj$) at NLO QCD+EW accuracy. In the case of $tZj$ production we consider both the case of the $Z$ boson undecayed and the complete final state $t \ell^+ \ell^- j$, including off-shell and non-resonant effects. We perform our calculation in the five-flavour-scheme (5FS), without selecting any specific production channel ($s$-, $t$- or $tW$ associated). Moreover, we provide a more realistic estimate of the theory uncertainty by carefully including the differences between the four-flavour-scheme (4FS) and 5FS predictions. The difficulties underlying this procedure in the presence of EW corrections are discussed in detail. We find that NLO EW corrections are in general within the NLO QCD theory uncertainties only if the flavour scheme uncertainty (4FS vs.~5FS) is taken into account. For the case of $t \ell^+ \ell^- j$ production we also investigate differences between NLO QCD+EW predictions and NLO QCD predictions matched with a parton shower simulation including multiple photon emissions.


title: Positivity among P-partition generating functions

abstract: We seek simple conditions on a pair of labeled posets that determine when the difference of their $(P,\omega)$-partition enumerators is $F$-positive, i.e., positive in Gessel's fundamental basis. This is a quasisymmetric analogue of the extensively studied problem of finding conditions on a pair of skew shapes that determine when the difference of their skew Schur functions is Schur-positive. We determine necessary conditions and separate sufficient conditions for $F$-positivity, and show that a broad operation for combining posets preserves positivity properties. We conclude with classes of posets for which we have conditions that are both necessary and sufficient.


title: Flexible Mixture Priors for Large Time-varying Parameter Models

abstract: Time-varying parameter (TVP) models often assume that the TVPs evolve according to a random walk. This assumption, however, might be questionable since it implies that coefficients change smoothly and in an unbounded manner. In this paper, we relax this assumption by proposing a flexible law of motion for the TVPs in large-scale vector autoregressions (VARs). Instead of imposing a restrictive random walk evolution of the latent states, we carefully design hierarchical mixture priors on the coefficients in the state equation. These priors effectively allow for discriminating between periods where coefficients evolve according to a random walk and times where the TVPs are better characterized by a stationary stochastic process. Moreover, this approach is capable of introducing dynamic sparsity by pushing small parameter changes towards zero if necessary. The merits of the model are illustrated by means of two applications. Using synthetic data we show that our approach yields precise parameter estimates. When applied to US data, the model reveals interesting patterns of low-frequency dynamics in coefficients and forecasts well relative to a wide range of competing models.


title: Probing phase transitions of holographic entanglement entropy with fixed area states

abstract: Recent results suggest that new corrections to holographic entanglement entropy should arise near phase transitions of the associated Ryu-Takayanagi (RT) surface. We study such corrections by decomposing the bulk state into fixed-area states and conjecturing that a certain `diagonal approximation' will hold. In terms of the bulk Newton constant $G$, this yields a correction of order $O(G^{-1/2})$ near such transitions, which is in particular larger than generic corrections from the entanglement of bulk quantum fields. However, the correction becomes exponentially suppressed away from the transition. The net effect is to make the entanglement a smooth function of all parameters, turning the RT `phase transition' into a crossover already at this level of analysis. We illustrate this effect with explicit calculations (again assuming our diagonal approximation) for boundary regions given by a pair of disconnected intervals on the boundary of the AdS$_3$ vacuum and for a single interval on the boundary of the BTZ black hole. In a natural large-volume limit where our diagonal approximation clearly holds, this second example verifies that our results agree with general predictions made by Murthy and Srednicki in the context of chaotic many-body systems. As a further check on our conjectured diagonal approximation, we show that it also reproduces the $O(G^{-1/2})$ correction found Penington et al for an analogous quantum RT transition. Our explicit computations also illustrate the cutoff-dependence of fluctuations in RT-areas.


title: LiHoF$_4$: Cuboidal Demagnetizing Factor in an Ising Ferromagnet

abstract: The demagnetizing factor has an important effect on the physics of ferromagnets. For cuboidal samples it depends on susceptibility and the historic problem of determining this function continues to generate theoretical and experimental challenges. To test a recent theory, we measure the magnetic susceptibility of the Ising dipolar ferromagnet LiHoF$_4$, using samples of varying aspect ratio, and we reconsider the demagnetizing transformation necessary to obtain the intrinsic material susceptibility. Our experimental results confirm that the microscopic details of the material significantly affect the transformation, as predicted. In particular, we find that the uniaxial Ising spins require a demagnetizing transformation that differs from the one needed for Heisenberg spins and that use of the wrong demagnetizing transformation would result in unacceptably large errors in the measured physical properties of the system. Our results further shed light on the origin of the mysterious `flat' susceptibility of ordered ferromagnets by demonstrating that the intrinsic susceptibility of the ordered ferromagnetic phase is infinite, regardless of sample shape.


title: Faster Secure Data Mining via Distributed Homomorphic Encryption

abstract: Due to the rising privacy demand in data mining, Homomorphic Encryption (HE) is receiving more and more attention recently for its capability to do computations over the encrypted field. By using the HE technique, it is possible to securely outsource model learning to the not fully trustful but powerful public cloud computing environments. However, HE-based training scales badly because of the high computation complexity. It is still an open problem whether it is possible to apply HE to large-scale problems. In this paper, we propose a novel general distributed HE-based data mining framework towards one step of solving the scaling problem. The main idea of our approach is to use the slightly more communication overhead in exchange of shallower computational circuit in HE, so as to reduce the overall complexity. We verify the efficiency and effectiveness of our new framework by testing over various data mining algorithms and benchmark data-sets. For example, we successfully train a logistic regression model to recognize the digit 3 and 8 within around 5 minutes, while a centralized counterpart needs almost 2 hours.


title: Housing Market Prediction Problem using Different Machine Learning Algorithms: A Case Study

abstract: Developing an accurate prediction model for housing prices is always needed for socio-economic development and well-being of citizens. In this paper, a diverse set of machine learning algorithms such as XGBoost, CatBoost, Random Forest, Lasso, Voting Regressor, and others, are being employed to predict the housing prices using public available datasets. The housing datasets of 62,723 records from January 2015 to November 2019 are obtained from Florida Volusia County Property Appraiser website. The records are publicly available and include the real estate or economic database, maps, and other associated information. The database is usually updated weekly according to the State of Florida regulations. Then, the housing price prediction models using machine learning techniques are developed and their regression model performances are compared. Finally, an improved housing price prediction model for assisting the housing market is proposed. Particularly, a house seller or buyer, or a real estate broker can get insight in making better-informed decisions considering the housing price prediction. The empirical results illustrate that based on prediction model performance, Coefficient of Determination (R2), Mean Square Error (MSE), Mean Absolute Error (MAE), and computational time, the XGBoost algorithm performs superior to the other models to predict the housing price.


title: Extensively Matching for Few-shot Learning Event Detection

abstract: Current event detection models under super-vised learning settings fail to transfer to newevent types. Few-shot learning has not beenexplored in event detection even though it al-lows a model to perform well with high gener-alization on new event types. In this work, weformulate event detection as a few-shot learn-ing problem to enable to extend event detec-tion to new event types. We propose two novelloss factors that matching examples in the sup-port set to provide more training signals to themodel. Moreover, these training signals can beapplied in many metric-based few-shot learn-ing models. Our extensive experiments on theACE-2005 dataset (under a few-shot learningsetting) show that the proposed method can im-prove the performance of few-shot learning


title: Supermassive black holes in cosmological simulations I: M_BH-M_star relation and black hole mass function

abstract: The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.


title: Nearly Optimal Robust Method for Convex Compositional Problems with Heavy-Tailed Noise

abstract: In this paper, we propose robust stochastic algorithms for solving convex compositional problems of the form $f(\E_\xi g(\cdot; \xi)) + r(\cdot)$ by establishing {\bf sub-Gaussian confidence bounds} under weak assumptions about the tails of noise distribution, i.e., {\bf heavy-tailed noise} with bounded second-order moments. One can achieve this goal by using an existing boosting strategy that boosts a low probability convergence result into a high probability result. However, piecing together existing results for solving compositional problems suffers from several drawbacks: (i) the boosting technique requires strong convexity of the objective; (ii) it requires a separate algorithm to handle non-smooth $r$; (iii) it also suffers from an additional polylogarithmic factor of the condition number. To address these issues, we directly develop a single-trial stochastic algorithm for minimizing optimal strongly convex compositional objectives, which has a nearly optimal high probability convergence result matching the lower bound of stochastic strongly convex optimization up to a logarithmic factor. To the best of our knowledge, this is the first work that establishes nearly optimal sub-Gaussian confidence bounds for compositional problems under heavy-tailed assumptions.


title: Towards Recurrent Autoregressive Flow Models

abstract: Stochastic processes generated by non-stationary distributions are difficult to represent with conventional models such as Gaussian processes. This work presents Recurrent Autoregressive Flows as a method toward general stochastic process modeling with normalizing flows. The proposed method defines a conditional distribution for each variable in a sequential process by conditioning the parameters of a normalizing flow with recurrent neural connections. Complex conditional relationships are learned through the recurrent network parameters. In this work, we present an initial design for a recurrent flow cell and a method to train the model to match observed empirical distributions. We demonstrate the effectiveness of this class of models through a series of experiments in which models are trained on three complex stochastic processes. We highlight the shortcomings of our current formulation and suggest some potential solutions.


title: Spectral Inclusion and Pollution for a Class of Dissipative Perturbations

abstract: Spectral inclusion and spectral pollution results are proved for sequences of linear operators of the form $T_0 + i \gamma s_n$ on a Hilbert space, where $s_n$ is strongly convergent to the identity operator and $\gamma > 0$. We work in both an abstract setting and a more concrete Sturm-Liouville framework. The results provide rigorous justification for a method of computing eigenvalues in spectral gaps.


title: The pressure-enhanced superconducting phase of Sr$_x$-Bi$_2$Se$_3$ probed by hard point contact spectroscopy

abstract: The superconducting systems emerging from topological insulators upon metal ion intercalation or application of high pressure are ideal for investigation of possible topological superconductivity. In this context, Sr-intercalated Bi$_2$Se$_3$ is specially interesting because it displays pressure induced re-entrant superconductivity where the high pressure phase shows almost two times higher $T_c$ than the ambient superconducting phase ( $T_C\sim$ 2.9 K). Interestingly, unlike the ambient phase, the pressure-induced superconducting phase shows strong indication of unconventional superconductivity. However, since the pressure-induced phase remains inaccessible to spectroscopic techniques, the detailed study of the phase remained an unattained goal. Here we show that the high-pressure phase can be realized under a mesoscopic point contact, where transport spectroscopy can be used to probe the spectroscopic properties of the pressure-induced phase. We find that the point contact junctions on the high-pressure phase show unusual response to magnetic field supporting the possibility of unconventional superconductivity.


title: Mass enhancement in 3d and s-p perovskites from symmetry breaking

abstract: In some d-electron oxides the measured effective mass m(exptl) has long been known to be significantly larger than the model effective mass m(model) deduced from mean-field band theory, i.e., m(exptl) = beta * m(model), where beta > 1 is the "mass enhancement", or "mass renormalization" factor. Previous applications of density functional theory (DFT), based on the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement, a fact that was taken as evidence of strong electronic correlation being the exclusive enabling physics. The current paper reports that known modalities of energy-lowering symmetry-broken spin and structural effects included in mean-field DFT show mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO3, SrTiO3, BaTiO3, and LaMnO3 as well as p-electron perovskites CsPbI3 and SrBiO3, including metals (SrVO3) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. The paper analyzes the contributions of different symmetry-broken modalities to mass enhancement, finds common effects in the range of d- as well as p-electron perovskites. Thus, symmetry breaking in mean-field theory could describe effects that were previously attributed exclusively to complex correlated treatments. Indeed, broken symmetries are often experimentally observed, e.g., octahedral tilting, Jahn-Teller distortions, or bond disproportionation, and often provide intuitive explanations in terms of degeneracy removal due to lowering of structural symmetry, or as shifting of band energies in explainable directions.


title: Epitaxy, exfoliation, and strain-induced magnetism in rippled Heusler membranes

abstract: Single-crystalline membranes of functional materials enable the tuning of properties via extreme strain states; however, conventional routes for producing membranes require the use of sacrificial layers and chemical etchants, which can both damage the membrane and limit the ability to make them ultrathin. Here we demonstrate the epitaxial growth of the cubic Heusler compound GdPtSb on graphene-terminated Al$_2$O$_3$ substrates. Despite the presence of the graphene interlayer, the Heusler films have epitaxial registry to the underlying sapphire, as revealed by x-ray diffraction, reflection high energy electron diffraction, and transmission electron microscopy. The weak Van der Waals interactions of graphene enable mechanical exfoliation to yield free-standing GdPtSb membranes, which form ripples when transferred to a flexible polymer handle. Whereas unstrained GdPtSb is antiferromagnetic, measurements on rippled membranes show a spontaneous magnetic moment at room temperature, with a saturation magnetization of 5.2 bohr magneton per Gd. First-principles calculations show that the coupling to homogeneous strain is too small to induce ferromagnetism, suggesting a dominant role for strain gradients. Our membranes provide a novel platform for tuning the magnetic properties of intermetallic compounds via strain (piezomagnetixm and magnetostriction) and strain gradients (flexomagnetism).


title: Ab Initio Modeling of Phonon-Assisted Relaxation of Electrons and Excitons in Semiconductor Nanocrystals for Multiexciton Generation

abstract: Electron-phonon and exciton-phonon interactions in nanoclusters are formulated and computed under the framework of GW-BSE (Bethe-Salpeter equation) approach. The phonon effect is modeled with the two-particle representation for the first time. The nonradiative relaxation rates of electrons and excitons are calculated. It is uncovered that both single-phonon relaxation and multiple-phonon relaxation are significant in nanocrystals, and correspond to two types of physical processes that have totally different spectral lineshapes. Furthermore, the multiple-phonon relaxation always occurs and its rates are comparable to the corresponding single-phonon relaxation rates for both electrons and excitons in the system studied (Si46). The inelastic scattering rates of electrons and excitons are also calculated based on many-body Green function theory. For the electronic states in Si46, the inelastic scattering decay is predicted to be a primary decay mechanism for multiexciton relaxation, and nonradiative relaxation rates are larger than inelastic scattering rates for most excitonic states in Si46.


title: Capturing Label Characteristics in VAEs

abstract: We present a principled approach to incorporating labels in VAEs that captures the rich characteristic information associated with those labels. While prior work has typically conflated these by learning latent variables that directly correspond to label values, we argue this is contrary to the intended effect of supervision in VAEs-capturing rich label characteristics with the latents. For example, we may want to capture the characteristics of a face that make it look young, rather than just the age of the person. To this end, we develop the CCVAE, a novel VAE model and concomitant variational objective which captures label characteristics explicitly in the latent space, eschewing direct correspondences between label values and latents. Through judicious structuring of mappings between such characteristic latents and labels, we show that the CCVAE can effectively learn meaningful representations of the characteristics of interest across a variety of supervision schemes. In particular, we show that the CCVAE allows for more effective and more general interventions to be performed, such as smooth traversals within the characteristics for a given label, diverse conditional generation, and transferring characteristics across datapoints.


title: Is Network the Bottleneck of Distributed Training?

abstract: Recently there has been a surge of research on improving the communication efficiency of distributed training. However, little work has been done to systematically understand whether the network is the bottleneck and to what extent. In this paper, we take a first-principles approach to measure and analyze the network performance of distributed training. As expected, our measurement confirms that communication is the component that blocks distributed training from linear scale-out. However, contrary to the common belief, we find that the network is running at low utilization and that if the network can be fully utilized, distributed training can achieve a scaling factor of close to one. Moreover, while many recent proposals on gradient compression advocate over 100x compression ratio, we show that under full network utilization, there is no need for gradient compression in 100 Gbps network. On the other hand, a lower speed network like 10 Gbps requires only 2x--5x gradients compression ratio to achieve almost linear scale-out. Compared to application-level techniques like gradient compression, network-level optimizations do not require changes to applications and do not hurt the performance of trained models. As such, we advocate that the real challenge of distributed training is for the network community to develop high-performance network transport to fully utilize the network capacity and achieve linear scale-out.


title: A Streaming Machine Learning Framework for Online Aggression Detection on Twitter

abstract: The rise of online aggression on social media is evolving into a major point of concern. Several machine and deep learning approaches have been proposed recently for detecting various types of aggressive behavior. However, social media are fast paced, generating an increasing amount of content, while aggressive behavior evolves over time. In this work, we introduce the first, practical, real-time framework for detecting aggression on Twitter via embracing the streaming machine learning paradigm. Our method adapts its ML classifiers in an incremental fashion as it receives new annotated examples and is able to achieve the same (or even higher) performance as batch-based ML models, with over 90% accuracy, precision, and recall. At the same time, our experimental analysis on real Twitter data reveals how our framework can easily scale to accommodate the entire Twitter Firehose (of 778 million tweets per day) with only 3 commodity machines. Finally, we show that our framework is general enough to detect other related behaviors such as sarcasm, racism, and sexism in real time.


title: On the origin of the quantum group symmetry in 3d quantum gravity

abstract: It is well-known that quantum groups are relevant to describe the quantum regime of 3d gravity. They encode a deformation of the gauge symmetries parametrized by the value of the cosmological constant. They appear as a form of regularization either through the quantization of the Chern-Simons formulation or the state sum approach of Turaev-Viro. Such deformations are perplexing from a continuum and classical picture since the action is defined in terms of undeformed gauge invariance. We present here a novel way to derive from first principle and from the classical action such quantum group deformation. The argument relies on two main steps. First we perform a canonical transformation, which deformed the gauge invariance and the boundary symmetries, and makes them depend on the cosmological constant. Second we implement a discretization procedure relying on a truncation of the degrees of freedom from the continuum.


title: Kinematics of an isolated galaxy, CIG 993

abstract: We present scanning Fabry-Perot interferometric observations of CIG 993, an apparently isolated luminous infrared galaxy also exhibiting luminous blue compact and Wolf-Rayet galaxy features, as well as high star formation rate. Our high resolution observations of the H$\alpha$ emission line allowed us to derive the radial velocity field, as well as the velocity dispersion $\sigma$, and the residual velocity fields of the galaxy. This galaxy exhibits several kinematical components. On one hand, the velocity gradients detected on the velocity field can be associated with a rotating disk -contrary to previous results with less spectral resolution-. However, the velocity field, the $\sigma$ and residual velocity field show significant deviations from circular motions in the central part of the galaxy that matches a region with high number of Wolf-Rayet and O stars, coincident with the blue luminous component of the galaxy. We find narrow and broad velocity components for the ionized gas in the central part of the galaxy. The broad component is evidence of a central outflow related with the on-going burst of stellar formation. The morpho-kinematical analysis of the galaxy indicates we are only seeing the brightest parts of the galaxy which correspond to the bulge, a central bar and the beginning of the disk. We believe CIG 993 is a disk galaxy harboring important star-forming processes most likely caused by a relatively recent interaction. This could imply that small encounters could change the global characteristics of a galaxy without disturbing the main rotation disk motion, nor the morphology of the galaxy.


title: Right-truncated Archimedean and related copulas

abstract: The copulas of random vectors with standard uniform univariate margins truncated from the right are considered and a general formula for such right-truncated conditional copulas is derived. This formula is analytical for copulas that can be inverted analytically as functions of each single argument. This is the case, for example, for Archimedean and related copulas. The resulting right-truncated Archimedean copulas are not only analytically tractable but can also be characterized as tilted Archimedean copulas. This finding allows one, for example, to more easily derive analytical properties such as the coefficients of tail dependence or sampling procedures of right-truncated Archimedean copulas. As another result, one can easily obtain a limiting Clayton copula for a general vector of truncation points converging to zero; this is an important property for (re)insurance and a fact already known in the special case of equal truncation points, but harder to prove without aforementioned characterization. Furthermore, right-truncated Archimax copulas with logistic stable tail dependence functions are characterized as tilted outer power Archimedean copulas and an analytical form of right-truncated nested Archimedean copulas is also derived.


title: Simple and Principled Uncertainty Estimation with Deterministic Deep Learning via Distance Awareness

abstract: Bayesian neural networks (BNN) and deep ensembles are principled approaches to estimate the predictive uncertainty of a deep learning model. However their practicality in real-time, industrial-scale applications are limited due to their heavy memory and inference cost. This motivates us to study principled approaches to high-quality uncertainty estimation that require only a single deep neural network (DNN). By formalizing the uncertainty quantification as a minimax learning problem, we first identify input distance awareness, i.e., the model's ability to quantify the distance of a testing example from the training data in the input space, as a necessary condition for a DNN to achieve high-quality (i.e., minimax optimal) uncertainty estimation. We then propose Spectral-normalized Neural Gaussian Process (SNGP), a simple method that improves the distance-awareness ability of modern DNNs, by adding a weight normalization step during training and replacing the output layer with a Gaussian process. On a suite of vision and language understanding tasks and on modern architectures (Wide-ResNet and BERT), SNGP is competitive with deep ensembles in prediction, calibration and out-of-domain detection, and outperforms the other single-model approaches.


title: Nash SIR: An Economic-Epidemiological Model of Strategic Behavior During a Viral Epidemic

abstract: This paper develops a Nash-equilibrium extension of the classic SIR model of infectious-disease epidemiology ("Nash SIR"), endogenizing people's decisions whether to engage in economic activity during a viral epidemic and allowing for complementarity in social-economic activity. An equilibrium epidemic is one in which Nash equilibrium behavior during the epidemic generates the epidemic. There may be multiple equilibrium epidemics, in which case the epidemic trajectory can be shaped through the coordination of expectations, in addition to other sorts of interventions such as stay-at-home orders and accelerated vaccine development. An algorithm is provided to compute all equilibrium epidemics.


title: Exergames for telerehabilitation

abstract: Recent advancements in technology have improved the connectivity between humans enhancing the transfer of information. Leveraging these technological marvels in the healthcare industry has led to the development of telehealth allowing patients and clinicians to receive and administer treatment remotely. Telerehabilitation is a subset of telehealth that facilitates remote rehabilitation treatment for patients. Providing rehabilitative services to the aging baby boomer population requires tech-savvy solutions to augment the therapists and clinicians for effective remote monitoring and tele-medicine. Hence, this thesis develops easy-to-use exergames for low-cost mechatronic devices targeting rehabilitation of post-stroke patients. Specifically, it demonstrates wearable inertial sensors for exergames consisting of an animated virtual coach for providing patients with instructions for performing range of motion exercises. Next, a gaming environment is developed for task-specific rehabilitation such as eating. Finally, exergames are developed for rehabilitation of pincer grasping. In addition to gamified interfaces providing an engaging rehabilitation experience to the user, the data acquired from the mechatronic devices facilitate data-driven telerehabilitation.


title: Pendant Drop Tensiometry: A Machine Learning Approach

abstract: Modern pendant drop tensiometry relies on numerical solution of the Young-Laplace equation and allow to determine the surface tension from a single picture of a pendant drop with high precision. Most of these techniques solve the Young-Laplace equation many times over to find the material parameters that provide a fit to a supplied image of a real droplet. Here we introduce a machine learning approach to solve this problem in a computationally more efficient way. We train a deep neural network to determine the surface tension of a given droplet shape using a large training set of numerically generated droplet shapes. We show that the deep learning approach is superior to the current state of the art shape fitting approach in speed and precision, in particular if shapes in the training set reflect the sensitivity of the droplet shape with respect to surface tension. In order to derive such an optimized training set we clarify the role of the Worthington number as quality indicator in conventional shape fitting and in the machine learning approach. Our approach demonstrates the capabilities of deep neural networks in the material parameter determination from rheological deformation experiments in general.

categories: cond-mat.soft cs.LG physics.chem-ph physics.comp-ph


title: Interface learning of multiphysics and multiscale systems

abstract: Complex natural or engineered systems comprise multiple characteristic scales, multiple spatiotemporal domains, and even multiple physical closure laws. To address such challenges, we introduce an interface learning paradigm and put forth a data-driven closure approach based on memory embedding to provide physically correct boundary conditions at the interface. To enable the interface learning for hyperbolic systems by considering the domain of influence and wave structures into account, we put forth the concept of upwind learning towards a physics-informed domain decomposition. The promise of the proposed approach is shown for a set of canonical illustrative problems. We highlight that high-performance computing environments can benefit from this methodology to reduce communication costs among processing units in emerging machine learning ready heterogeneous platforms toward exascale era.


title: From biophysical to integrate-and-fire modelling

abstract: This paper proposes a methodology to extract a low-dimensional integrate-and-fire model from an arbitrarily detailed single-compartment biophysical model. The method aims at relating the modulation of maximal conductance parameters in the biophysical model to the modulation of parameters in the proposed integrate-and-fire model. The approach is illustrated on two well-documented examples of cellular neuromodulation: the transition between Type I and Type II excitability and the transition between spiking and bursting.


title: Constraint-Based Regularization of Neural Networks

abstract: We propose a method for efficiently incorporating constraints into a stochastic gradient Langevin framework for the training of deep neural networks. Constraints allow direct control of the parameter space of the model. Appropriately designed, they reduce the vanishing/exploding gradient problem, control weight magnitudes and stabilize deep neural networks and thus improve the robustness of training algorithms and the generalization capabilities of the trained neural network. We present examples of constrained training methods motivated by orthogonality preservation for weight matrices and explicit weight normalizations. We describe the methods in the overdamped formulation of Langevin dynamics and the underdamped form, in which momenta help to improve sampling efficiency. The methods are explored in test examples in image classification and natural language processing.


title: Mimicking the TYC strategy: Weak Allee effects, and a nonhyperbolic extinction boundary

abstract: The Trojan Y Chromosome strategy (TYC) is a genetic biocontrol strategy designed to alter the sex ratio of a target invasive population by reducing the number of females over time. Recently an alternative strategy is introduced, that mimics the TYC strategy by harvesting females whilst stocking males . We consider the FHMS strategy, with a weak Allee effect, and show that the extinction boundary need note be hyperbolic. To the best of our knowledge, this is the first example of a non-hyperbolic extinction boundary in mating models, structured by sex. Next, we consider the spatially explicit model and show that the weak Allee effect is both sufficient and necessary for Turing patterns to occur. We discuss the applicability of our results to large scale biocontrol, as well as compare and contrast our results to the case with a strong Allee effect.


title: Tuning Magnetic Order in the van der Waals Metal Fe5GeTe2 by Cobalt Substitution

abstract: Fe5-xGeTe2 is a van der Waals material with one of the highest reported bulk Curie temperatures, $T_C$ ~ 310K. In this study, theoretical calculations and experiments are utilized to demonstrate that the magnetic ground state is highly sensitive to local atomic arrangements and the interlayer stacking. Cobalt substitution is found to be an effective way to manipulate the magnetic properties while also increasing the ordering temperature. In particular, cobalt substitution up to 30% enhances $T_C$ and changes the magnetic anisotropy, while approximately 50% cobalt substitution yields an antiferromagnetic state. Single crystal x-ray diffraction evidences a structural change upon increasing the cobalt concentration, with a rhombohedral cell observed in the parent material and a primitive cell observed for ~46% cobalt content relative to iron. First principles calculations demonstrate that it is a combination of high cobalt content and the concomitant change to primitive layer stacking that produces antiferromagnetic order. These results illustrate the sensitivity of magnetism in Fe5-xGeTe2 to composition and structure, and emphasize the important role of structural order/disorder and layer stacking in cleavable magnetic materials.


title: Heterogeneous partition of cellular blood-borne nanoparticles through microvascular bifurcations

abstract: Blood flowing through microvascular bifurcations has been an active research topic for many decades, while the partitioning pattern of nanoscale solutes in the blood remains relatively unexplored. Here, we demonstrate a multiscale computational framework for direct numerical simulation of the nanoparticle (NP) partitioning through physiologically-relevant vascular bifurcations in the presence of red blood cells (RBCs). The computational framework is established by embedding a newly-developed particulate suspension inflow/outflow boundary condition into a multiscale blood flow solver. The computational framework is verified by recovering a tubular blood flow without a bifurcation and validated against the experimental measurement of an intravital bifurcation flow. The classic Zweifach-Fung (ZF) effect is shown to be well captured by the method. Moreover, we observe that NPs exhibit a ZF-like heterogeneous partition in response to the heterogeneous partition of the RBC phase. The NP partitioning prioritizes the high-flow-rate daughter branch except for extreme (large or small) suspension flow partition ratios under which the complete phase separation tends to occur. By analyzing the flow field and the particle trajectories, we show that the ZF-like heterogeneity in NP partition can be explained by the RBC-entrainment effect caused by the deviation of the flow separatrix preceded by the tank-treading of RBCs near the bifurcation junction. The recovery of homogeneity in the NP partition under extreme flow partition ratios is due to the plasma skimming of NPs in the cell-free layer. These findings, based on the multiscale computational framework, provide biophysical insights to the heterogeneous distribution of NPs in microvascular beds that are observed pathophysiologically.


title: The Unusual Eruption of the Extragalactic Classical Nova M31N 2017-09a

abstract: M31N 2017-09a is a classical nova and was observed for some 160 days following its initial eruption, during which time it underwent a number of bright secondary outbursts. The light-curve is characterized by continual variation with excursions of at least 0.5 magnitudes on a daily time-scale. The lower envelope of the eruption suggests that a single power-law can describe the decline rate. The eruption is relatively long with $t_2 = 111$, and $t_3 = 153$ days.


title: Markovian RNN: An Adaptive Time Series Prediction Network with HMM-based Switching for Nonstationary Environments

abstract: We investigate nonlinear regression for nonstationary sequential data. In most real-life applications such as business domains including finance, retail, energy and economy, timeseries data exhibits nonstationarity due to the temporally varying dynamics of the underlying system. We introduce a novel recurrent neural network (RNN) architecture, which adaptively switches between internal regimes in a Markovian way to model the nonstationary nature of the given data. Our model, Markovian RNN employs a hidden Markov model (HMM) for regime transitions, where each regime controls hidden state transitions of the recurrent cell independently. We jointly optimize the whole network in an end-to-end fashion. We demonstrate the significant performance gains compared to vanilla RNN and conventional methods such as Markov Switching ARIMA through an extensive set of experiments with synthetic and real-life datasets. We also interpret the inferred parameters and regime belief values to analyze the underlying dynamics of the given sequences.


title: Simulations of M87 and Sgr A* Imaging with the Millimetron Space Observatory on near-Earth orbits

abstract: High resolution imaging of supermassive black holes shadows is a direct way to verify the theory of general relativity at extreme gravity conditions. Very Long Baseline Interferometry (VLBI) observations at millimeter/sub-millimeter wavelengths can provide such angular resolution for supermassive black holes, located in Sgr A* and M87. Recent VLBI observations of M87 with the Event Horizon Telescope (EHT) has shown such capabilities. The maximum obtainable spatial resolution of EHT is limited by Earth diameter and atmospheric phase variations. In order to improve the image resolution longer baselines are required. Radioastron space mission has successfully demonstrated the capabilities of Space-Earth VLBI with baselines much larger than Earth diameter. Millimetron is a next space mission of the Russian Space Agency that will operate at millimeter wavelengths. Nominal orbit of the observatory will be located around Lagrangian L2 point of the Sun-Earth system. In order to optimize the VLBI mode, we consider a possible second stage of the mission that could use near-Earth high elliptical orbit (HEO). In this contribution a set of near-Earth orbits is used for the synthetic space-ground VLBI observations of Sgr A* and M87 in joint Millimetron and EHT configuration. General-relativistic magnetohydrodynamic models (GRMHD) for black hole environment of Sgr A* and M87 are used for static and dynamic imaging simulations at 230 GHz. A comparison preformed between ground and space-ground baselines demonstrates that joint observations with Millimetron and EHT significantly improve the image resolution and allow the EHT+Millimetron to obtain snapshot images of Sgr A* probing dynamics at fast timescales.


title: Learning-Based Real-Time Event Identification Using Rich Real PMU Data

abstract: A large-scale deployment of phasor measurement units (PMUs) that reveal the inherent physical laws of power systems from a data perspective enables an enhanced awareness of power system operation. However, the high-granularity and non-stationary nature of PMU time series and imperfect data quality could bring great technical challenges to real-time system event identification. To address these issues, this paper proposes a two-stage learning-based framework. At the first stage, a Markov transition field (MTF) algorithm is exploited to extract the latent data features by encoding temporal dependency and transition statistics of PMU data in graphs. Then, a spatial pyramid pooling (SPP)-aided convolutional neural network (CNN) is established to efficiently and accurately identify operation events. The proposed method fully builds on and is also tested on a large real dataset from several tens of PMU sources (and the corresponding event logs), located across the U.S., with a time span of two consecutive years. The numerical results validate that our method has high identification accuracy while showing good robustness against poor data quality.


title: Neuromechanical Mechanisms of Gait Adaptation in C. elegans: Relative Roles of Neural and Mechanical Coupling

abstract: Understanding principles of neurolocomotion requires the synthesis of neural activity, sensory feedback, and biomechanics. The nematode \textit{C. elegans} is an ideal model organism for studying locomotion in an integrated neuromechanical setting because its neural circuit has a well-characterized modular structure and its undulatory forward swimming gait adapts to the surrounding fluid with a shorter wavelength in higher viscosity environments. This adaptive behavior emerges from the neural modules interacting through a combination of mechanical forces, neuronal coupling, and sensory feedback mechanisms. However, the relative contributions of these coupling modes to gait adaptation are not understood. The model consists of repeated neuromechanical modules that are coupled through the mechanics of the body, short-range proprioception, and gap-junctions. The model captures the experimentally observed gait adaptation over a wide range of mechanical parameters, provided that the muscle response to input from the nervous system is faster than the body response to changes in internal and external forces. The modularity of the model allows the use of the theory of weakly coupled oscillators to identify the relative roles of body mechanics, gap-junctional coupling, and proprioceptive coupling in coordinating the undulatory gait. The analysis shows that the wavelength of body undulations is set by the relative strengths of these three coupling forms. In a low-viscosity fluid environment, the competition between gap-junctions and proprioception produces a long wavelength undulation, which is only achieved in the model with sufficiently strong gap-junctional coupling.The experimentally observed decrease in wavelength in response to increasing fluid viscosity is the result of an increase in the relative strength of mechanical coupling, which promotes a short wavelength.


title: A block coordinate descent optimizer for classification problems exploiting convexity

abstract: Second-order optimizers hold intriguing potential for deep learning, but suffer from increased cost and sensitivity to the non-convexity of the loss surface as compared to gradient-based approaches. We introduce a coordinate descent method to train deep neural networks for classification tasks that exploits global convexity of the cross-entropy loss in the weights of the linear layer. Our hybrid Newton/Gradient Descent (NGD) method is consistent with the interpretation of hidden layers as providing an adaptive basis and the linear layer as providing an optimal fit of the basis to data. By alternating between a second-order method to find globally optimal parameters for the linear layer and gradient descent to train the hidden layers, we ensure an optimal fit of the adaptive basis to data throughout training. The size of the Hessian in the second-order step scales only with the number weights in the linear layer and not the depth and width of the hidden layers; furthermore, the approach is applicable to arbitrary hidden layer architecture. Previous work applying this adaptive basis perspective to regression problems demonstrated significant improvements in accuracy at reduced training cost, and this work can be viewed as an extension of this approach to classification problems. We first prove that the resulting Hessian matrix is symmetric semi-definite, and that the Newton step realizes a global minimizer. By studying classification of manufactured two-dimensional point cloud data, we demonstrate both an improvement in validation error and a striking qualitative difference in the basis functions encoded in the hidden layer when trained using NGD. Application to image classification benchmarks for both dense and convolutional architectures reveals improved training accuracy, suggesting possible gains of second-order methods over gradient descent.


title: Improvements in Computation and Usage of Joint CDFs for the N-Dimensional Order Statistic

abstract: Order statistics provide an intuition for combining multiple lists of scores over a common index set. This intuition is particularly valuable when the lists to be combined cannot be directly compared in a sensible way. We describe here the advantages of a new method for using joint CDFs of such order statistics to combine score lists. We also present, with proof, a new algorithm for computing such joint CDF values, with runtime linear in the size of the combined list.


title: Sustainable Recreational Fishing Using a Novel Electrical Muscle Stimulation (EMS) Lure and Ensemble Network Algorithm to Maximize Catch and Release Survivability

abstract: With 200-700 million anglers in the world, sportfishing is nearly five times more common than commercial trawling. Worldwide, hundreds of thousands of jobs are linked to the sportfishing industry, which generates billions of dollars for water-side communities and fisheries conservatories alike. However, the sheer popularity of recreational fishing poses threats to aquatic biodiversity that are hard to regulate. For example, as much as 25% of overfished populations can be traced to anglers. This alarming statistic is explained by the average catch and release mortality rate of 43%, which primarily results from hook-related injuries and careless out-of-water handling. The provisional-patented design proposed in this paper addresses both these problems separately First, a novel, electrical muscle stimulation based fishing lure is proposed as a harmless and low cost alternative to sharp hooks. Early prototypes show a constant electrical current of 90 mA applied through a 200g European perch's jaw can support a reeling tension of 2N - safely within the necessary ranges. Second, a fish-eye camera bob is designed to wirelessly relay underwater footage to a smartphone app, where an ensemble convolutional neural network automatically classifies the fish's species, estimates its length, and cross references with local and state fishing regulations (ie. minimum size, maximum bag limit, and catch season). This capability reduces overfishing by helping anglers avoid accidentally violating guidelines and eliminates the need to reel the fish in and expose it to negligent handling. IN conjunction, this cheap, lightweight, yet high-tech invention is a paradigm shift in preserving a world favorite pastime; while at the same time making recreational fishing more sustainable.


title: Using Weighted P-Values in Fisher's Method

abstract: Fisher's method prescribes a way to combine p-values from multiple experiments into a single p-value. However, the original method can only determine a combined p-value analytically if all constituent p-values are weighted equally. Here we present, with proof, a method to combine p-values with arbitrary weights.


title: Crossed morphisms, (integration of) post-Lie algebras and the post-Lie Magnus expansion

abstract: In the first part of this letter it will be shown that the post-Lie Magnus expansion can be interpreted as a crossed morphism between two (local) Lie group. The second part will be devoted to present two combinatorial methods, both based on special tubings on planar trees, to compute the coefficients of this remarkable formal series.


title: The (black hole mass)-(color) relations for early- and late-type galaxies: red and blue sequences

abstract: [Abridged] Tight correlations between supermassive black hole (SMBH) mass ($M_{\rm BH}$) and the properties of the host galaxy have useful implications for our understanding of the growth of SMBHs and evolution of galaxies. Here, we present newly observed correlations between $M_{\rm BH}$ and the host galaxy total UV$-$ [3.6] color ($\mathcal{C_{\rm UV,tot}}$, Pearson's r = $0.6-0.7$) for a sample of 67 galaxies (20 early-type galaxies and 47 late-type galaxies) with directly measured $M_{\rm BH}$ in the GALEX/S$^{4}$G survey. The colors are carefully measured in a homogeneous manner using the galaxies' FUV, NUV and 3.6 $\micron$ magnitudes and their multi-component structural decompositions in the literature. We find that more massive SMBHs are hosted by (early- and late-type) galaxies with redder colors, but the $M_{\rm BH}- \mathcal{C_{\rm UV,tot}}$ relations for the two morphological types have slopes that differ at $\sim 2 \sigma$ level. Early-type galaxies define a red sequence in the $M_{\rm BH}- \mathcal{C_{\rm UV,tot}}$ diagrams, while late-type galaxies trace a blue sequence. Within the assumption that the specific star formation rate of a galaxy (sSFR) is well traced by $L_{\rm UV}/L_{\rm 3.6}$, it follows that the SMBH masses for late-type galaxies exhibit a steeper dependence on sSFR than those for early-type galaxies. The $M_{\rm BH}- \mathcal{C_{\rm UV,tot}}$ and $M_{\rm BH}-L_{\rm 3.6,tot}$ relations for the sample galaxies reveal a comparable level of vertical scatter in the log $M_{\rm BH}$ direction, roughly $5\%-27\%$ more than the vertical scatter of the $M_{\rm BH}-\sigma$ relation. Our $M_{\rm BH}- \mathcal{C_{\rm UV,tot}}$ relations suggest different channels of SMBH growth for early- and late-type galaxies, consistent with their distinct formation and evolution scenarios.


title: Smoothed Analysis of Online and Differentially Private Learning

abstract: Practical and pervasive needs for robustness and privacy in algorithms have inspired the design of online adversarial and differentially private learning algorithms. The primary quantity that characterizes learnability in these settings is the Littlestone dimension of the class of hypotheses [Ben-David et al., 2009, Alon et al., 2019]. This characterization is often interpreted as an impossibility result because classes such as linear thresholds and neural networks have infinite Littlestone dimension. In this paper, we apply the framework of smoothed analysis [Spielman and Teng, 2004], in which adversarially chosen inputs are perturbed slightly by nature. We show that fundamentally stronger regret and error guarantees are possible with smoothed adversaries than with worst-case adversaries. In particular, we obtain regret and privacy error bounds that depend only on the VC dimension and the bracketing number of a hypothesis class, and on the magnitudes of the perturbations.


title: Segment Distribution around the Center of Gravity of Branched Polymers

abstract: Mathematical expressions for mass distributions around the center of gravity are derived for branched polymers with the help of the Isihara formula. We introduce the Gaussian approximation for the end-to-end vector, $\vec{r}_{G\nu_{i}}$, from the center of gravity to the $i$th mass point on the $\nu$th arm. Then, for star polymers, the result is \begin{equation} \varphi_{star}(s)=\frac{1}{N}\sum_{\nu=1}^{f}\sum_{i=1}^{N_{\nu}}\left(\frac{d}{2\pi\left\langle r_{G\nu_{i}}^{2}\right\rangle}\right)^{d/2}\exp\left(-\frac{d}{2\left\langle r_{G\nu_{i}}^{2}\right\rangle}s^{2}\right)\notag \end{equation} for a sufficiently large $N$, where $f$ denotes the number of arms. It is found that the resultant $\varphi_{star}(s)$ is, unfortunately, not Gaussian. For dendrimers \begin{equation} \varphi_{dend}(s)=\sum_{h=1}^{g}\omega_{h}\left(\frac{d}{2pi\left\langle r_{G_{h}}^{2}\right\rangle}\right)^{d/2}\exp\left(-\frac{d}{2\left\langle r_{G_{h}}^{2}\right\rangle}s^{2}\right)\notag \end{equation} where $\omega_{h}$ denotes the weight fraction of masses in the $h$th generation on a dendrimer constructed from $g$ generations, so that $\sum_{h=1}^{g}\omega_{h}=1$. To be specific, $\omega_{1}=1/N$ and $\omega_{h}=(f-1)^{h-2}/N$ for $h\ge 2$. These distributions can be described by the same grand sum of each Gaussian function for the end-to-end distance from the center of gravity to each mass point. Note that for a large $f$ and $g$, the statistical weight of younger generations becomes dominant. As a consequence, the mass distribution of unperturbed dendrimers approaches the Gaussian form in the limit of a large $f$ and $g$. It is shown that the radii of gyration of dendrimers increase logarithmically with $N$, which leading to the exponent, $\nu_{0}=0$. An example of randomly branched polymers is also discussed.


title: A two-dimensional multi-class traffic flow model

abstract: The aim of this work is to introduce a two-dimensional macroscopic traffic model for multiple populations of vehicles. Starting from the paper [20], where a two-dimensional model for a single class of vehicles is proposed, we extend the dynamics to a multi-class model leading to a coupled system of conservation laws in two space dimensions. Besides the study of the Riemann problems we also present a Lax-Friedrichs type discretization scheme recovering the theoretical results by means of numerical tests. We calibrate the multi-class model with real data and compare the fitted model to the real trajectories. Finally, we test the ability of the model to simulate the overtaking of vehicles.


title: Interpreting the Latent Space of GANs via Correlation Analysis for Controllable Concept Manipulation

abstract: Generative adversarial nets (GANs) have been successfully applied in many fields like image generation, inpainting, super-resolution and drug discovery, etc., by now, the inner process of GANs is far from been understood. To get deeper insight of the intrinsic mechanism of GANs, in this paper, a method for interpreting the latent space of GANs by analyzing the correlation between latent variables and the corresponding semantic contents in generated images is proposed. Unlike previous methods that focus on dissecting models via feature visualization, the emphasis of this work is put on the variables in latent space, i.e. how the latent variables affect the quantitative analysis of generated results. Given a pretrained GAN model with weights fixed, the latent variables are intervened to analyze their effect on the semantic content in generated images. A set of controlling latent variables can be derived for specific content generation, and the controllable semantic content manipulation be achieved. The proposed method is testified on the datasets Fashion-MNIST and UT Zappos50K, experiment results show its effectiveness.


title: Fast Simulation of Magnetic Field Gradients for Optimization of Pulse Sequences

abstract: We study how to simulate, efficiently, pulse field gradients (PFG) used in nuclear magnetic resonance (NMR). An efficient simulation requires discretization in time and space. We study both discretizations and provide a guideline to choose best discretization values depending on the precision required experimentally. We provide a theoretical study and simulation showing the minimum number of divisions we need in space for simulating, with high precision, a sequence composed of several unitary evolution and PFG. We show that the fast simulation of PFG allow us to optimize sequences composed of PFG, radio-frequency pulses and free evolution, to implement non-unitary evolution (quantum channels). As an evidence of the success of our work, we performed two types of experiments. First, we implement two quantum channels and compare the results with their theoretical predictions. In the second experiment, we used the fast simulation of PFG to optimize and implement a sequence to prepare pseudo pure state with better signal to noise ratio than any known procedure till now.


title: Angular Independent Photonic Pigments via the Controlled Micellization of Amphiphilic Bottlebrush Block Copolymers

abstract: Photonic materials with angular independent structural colour are highly desirable because they offer the broad viewing angles required for application as colorants in paints, cosmetics, textiles or displays. However, they are challenging to fabricate as they require isotropic nanoscale architectures with only short-range correlation. In this article, porous microparticles with such a structure are produced in a single, scalable step from an amphiphilic bottlebrush block copolymer. This is achieved by exploiting a novel controlled micellization self-assembly mechanism within emulsified toluene-in-water droplets. By restricting water permeation through the droplet interface, the size of the pores can be precisely addressed, resulting in structurally coloured pigments. Furthermore, the reflected colour can be tuned to reflect across the full visible spectrum using only a single polymer (Mn = 290 kDa) by altering the initial emulsification conditions. Such photonic pigments have several key advantages over their crystalline analogues, as they provide isotropic structural coloration that suppresses iridescence and improves colour purity without the need for either refractive index matching or the inclusion of a broadband absorber.


title: M2Net: Multi-modal Multi-channel Network for Overall Survival Time Prediction of Brain Tumor Patients

abstract: Early and accurate prediction of overall survival (OS) time can help to obtain better treatment planning for brain tumor patients. Although many OS time prediction methods have been developed and obtain promising results, there are still several issues. First, conventional prediction methods rely on radiomic features at the local lesion area of a magnetic resonance (MR) volume, which may not represent the full image or model complex tumor patterns. Second, different types of scanners (i.e., multi-modal data) are sensitive to different brain regions, which makes it challenging to effectively exploit the complementary information across multiple modalities and also preserve the modality-specific properties. Third, existing methods focus on prediction models, ignoring complex data-to-label relationships. To address the above issues, we propose an end-to-end OS time prediction model; namely, Multi-modal Multi-channel Network (M2Net). Specifically, we first project the 3D MR volume onto 2D images in different directions, which reduces computational costs, while preserving important information and enabling pre-trained models to be transferred from other tasks. Then, we use a modality-specific network to extract implicit and high-level features from different MR scans. A multi-modal shared network is built to fuse these features using a bilinear pooling model, exploiting their correlations to provide complementary information. Finally, we integrate the outputs from each modality-specific network and the multi-modal shared network to generate the final prediction result. Experimental results demonstrate the superiority of our M2Net model over other methods.


title: Implementation of a Quantum Engine Fuelled by Information

abstract: Here, we show the implementation of a complete cycle of a quantum engine fuelled by information. This engine is a quantum version of the Szilard engine, where information is used to extract heat from the environment and fully convert it into work. In our experiment, this work is used to make a weight, initially in the ground state, reach its excited state. We measure the energy and the state of each component of the engine, after each step of the cycle, and compare them with the theoretical prediction to show that the cycle is implemented with high precision. We also perform experiments to show that the engine is well isolated from the environment after the heat extraction, and we measure the entropy of the weight to show the full conversion of heat into work. Thus, we successfully demonstrate that information can be used as a fuel for single-reservoir engines.


title: MoFlow: An Invertible Flow Model for Generating Molecular Graphs

abstract: Generating molecular graphs with desired chemical properties driven by deep graph generative models provides a very promising way to accelerate drug discovery process. Such graph generative models usually consist of two steps: learning latent representations and generation of molecular graphs. However, to generate novel and chemically-valid molecular graphs from latent representations is very challenging because of the chemical constraints and combinatorial complexity of molecular graphs. In this paper, we propose MoFlow, a flow-based graph generative model to learn invertible mappings between molecular graphs and their latent representations. To generate molecular graphs, our MoFlow first generates bonds (edges) through a Glow based model, then generates atoms (nodes) given bonds by a novel graph conditional flow, and finally assembles them into a chemically valid molecular graph with a posthoc validity correction. Our MoFlow has merits including exact and tractable likelihood training, efficient one-pass embedding and generation, chemical validity guarantees, 100\% reconstruction of training data, and good generalization ability. We validate our model by four tasks: molecular graph generation and reconstruction, visualization of the continuous latent space, property optimization, and constrained property optimization. Our MoFlow achieves state-of-the-art performance, which implies its potential efficiency and effectiveness to explore large chemical space for drug discovery.


title: An Online Method for A Class of Distributionally Robust Optimization with Non-Convex Objectives

abstract: In this paper, we propose a practical online method for solving a class of distributionally robust optimization (DRO) with non-convex objectives, which has important applications in machine learning for improving the robustness of neural networks. In the literature, most methods for solving DRO are based on stochastic primal-dual methods. However, primal-dual methods for DRO suffer from several drawbacks: (1) manipulating a high-dimensional dual variable corresponding to the size of data is time expensive; (2) they are not friendly to online learning where data is coming sequentially. To address these issues, we consider a class of DRO with an KL divergence regularization on the dual variables, transform the min-max problem into a compositional minimization problem, and propose practical duality-free online stochastic methods without requiring a large mini-batch size. We establish the state-of-the-art complexities of the proposed methods with and without a Polyak-\L ojasiewicz (PL) condition of the objective. Empirical studies on large-scale deep learning tasks (i) demonstrate that our method can speed up the training by more than 2 times than baseline methods and save days of training time on a large-scale dataset with $\sim$ 265K images, and (ii) verify the supreme performance of DRO over Empirical Risk Minimization (ERM) on imbalanced datasets. Of independent interest, the proposed method can be also used for solving a family of stochastic compositional problems with state-of-the-art complexities.


title: Is T Tauri North a 'Classical' T Tauri star?

abstract: We present high-resolution $H$- and $K$-band spectroscopic observations of the archetypal T Tauri star T Tau North. Synthetic spectral modeling is used to derive the $K$-band temperature, surface gravity, magnetic field strength, and rotational velocity for this star. The $K$-band spectroscopic temperature measured is $\rm T_{K\text{-band}} = 3976 \pm 90 K$, which is $\sim 1000 \, \rm K$ cooler than the temperature measured from optical observations. Our $K$-band temperature measurement for T Tau N is confirmed using equivalent width line ratio vs. temperature relations in the $H$-band, from which a $\rm T_{H\text{-band}} = 4085 \pm 155 K$ is derived. This optical vs. IR temperature difference is interpreted as cool or hot spots, or both, covering a significant part of the surface of T Tau N. The gravity derived for T Tau N $\log{g}=3.45 \pm 0.14$ is lower than the gravity of nearly every other star in a sample of 24 classical T Tauri stars in Taurus. Combining these temperature and gravity results with magnetic stellar evolutionary models, we find the age of T Tau N to be less than 1 Myr old. These results suggest that T Tau N is in an earlier evolutionary stage than most classical T Tauri stars in Taurus, arguing that it is a protostar ejected from the embedded southern binary system shortly after its formation.


title: Detecting Pulsar Polarization below 100 MHz with the Long Wavelength Array

abstract: Using the first station of the Long Wavelength Array (LWA1), we examine polarized pulsar emission between 25 and 88 MHz. Polarized light from pulsars undergoes Faraday rotation as it passes through the magnetized interstellar medium. Observations from low frequency telescopes are ideal for obtaining precise rotation measures (RMs) because the effect of Faraday rotation is proportional to the square of the observing wavelength. With these RMs, we obtained polarized pulse profiles to see how polarization changes in the 25-88 MHz range. The RMs were also used to derive values for the electron density weighted average Galactic magnetic field along the line of sight. We present rotation measures and polarization profiles of 15 pulsars acquired using data from LWA1. These results provide new insight into low-frequency polarization characteristics and pulsar emission heights, and complement measurements at higher frequencies.


title: Self-supervised Learning on Graphs: Deep Insights and New Direction

abstract: The success of deep learning notoriously requires larger amounts of costly annotated data. This has led to the development of self-supervised learning (SSL) that aims to alleviate this limitation by creating domain specific pretext tasks on unlabeled data. Simultaneously, there are increasing interests in generalizing deep learning to the graph domain in the form of graph neural networks (GNNs). GNNs can naturally utilize unlabeled nodes through the simple neighborhood aggregation that is unable to thoroughly make use of unlabeled nodes. Thus, we seek to harness SSL for GNNs to fully exploit the unlabeled data. Different from data instances in the image and text domains, nodes in graphs present unique structure information and they are inherently linked indicating not independent and identically distributed (or i.i.d.). Such complexity is a double-edged sword for SSL on graphs. On the one hand, it determines that it is challenging to adopt solutions from the image and text domains to graphs and dedicated efforts are desired. On the other hand, it provides rich information that enables us to build SSL from a variety of perspectives. Thus, in this paper, we first deepen our understandings on when, why, and which strategies of SSL work with GNNs by empirically studying numerous basic SSL pretext tasks on graphs. Inspired by deep insights from the empirical studies, we propose a new direction SelfTask to build advanced pretext tasks that are able to achieve state-of-the-art performance on various real-world datasets. The specific experimental settings to reproduce our results can be found in \url{https://github.com/ChandlerBang/SelfTask-GNN}.


title: A parallel hybrid implementation of the 2D acoustic wave equation

abstract: In this paper, we propose a hybrid parallel programming approach for a numerical solution of a two-dimensional acoustic wave equation using an implicit difference scheme for a single computer. The calculations are carried out in an implicit finite difference scheme. First, we transform the differential equation into an implicit finite-difference equation and then using the ADI method, we split the equation into two sub-equations. Using the cyclic reduction algorithm, we calculate an approximate solution. Finally, we change this algorithm to parallelize on GPU, GPU+OpenMP, and Hybrid (GPU+OpenMP+MPI) computing platforms. The special focus is on improving the performance of the parallel algorithms to calculate the acceleration based on the execution time. We show that the code that runs on the hybrid approach gives the expected results by comparing our results to those obtained by running the same simulation on a classical processor core, CUDA, and CUDA+OpenMP implementations.


title: Oscillating vector solitary waves in soft laminates

abstract: Vector solitary waves are nonlinear waves of coupled polarizations that propagate with constant velocity and shape. In mechanics, they hold the potential to control locomotion, mitigate shocks and transfer information, among other functionalities. Recently, such elastic waves were numerically observed in compressible rubber-like laminates. Here, we conduct numerical experiments to characterize the possible vector solitary waves in these laminates, and expose a new type of waves whose amplitude and velocity oscillate periodically without dispersing in time. This oscillation is a manifestation of a periodic transfer of energy between the two wave polarizations, which we consider as internal mode of the solitary wave. We find that the vector solitary waves propagate faster at higher amplitudes, and determine a lower bound for their velocity. We describe a procedure for identifying which initial strains generate such vector solitary waves. This procedure also enables an additional classification between tensile and compressive solitary waves, according to the way that the axial strain changes as the waves propagate.


title: Scanning Gate Microscopy of Localized States in a gate-defined Bilayer Graphene Channel

abstract: We use Scanning Gate Microscopy to demonstrate the presence of localized states arising from potential inhomogeneities in a 50nm-wide, gate-defined conducting channel in encapsulated bilayer graphene. When imaging the channel conductance under the influence of a local tip-induced potential, we observe ellipses of enhanced conductance as a function of the tip position. These ellipses allow us to infer the location of the localized states and to study their dependence on the displacement field. For large displacement fields, we observe that localized states tend to occur halfway into the channel. All our observations can be well explained within the framework of stochastic Coulomb blockade.


title: Using Sentiment Information for Preemptive Detection of Toxic Comments in Online Conversations

abstract: The challenge of automatic detection of toxic comments online has been the subject of a lot of research recently, but the focus has been mostly on detecting it in individual messages after they have been posted. Some authors have tried to predict if a conversation will derail into toxicity using the features of the first few messages. In this paper, we combine that approach with previous work on toxicity detection using sentiment information, and show how the sentiments expressed in the first messages of a conversation can help predict upcoming toxicity. Our results show that adding sentiment features does help improve the accuracy of toxicity prediction, and also allow us to make important observations on the general task of preemptive toxicity detection.


title: Special unitary particle pusher for extreme fields

abstract: An exact momentum update useful for particle codes was previously given. The expressions involved were unwieldy. By treating the momentum as a representation of SU(2)xSU(2) instead of SO(3,1), a more compact expression for the exact momentum update is obtained. An expansion in powers of the timestep can be formulated such that invariance properties are exactly preserved, and push rates comparable to the standard Boris pusher are obtained.


title: Are you wearing a mask? Improving mask detection from speech using augmentation by cycle-consistent GANs

abstract: The task of detecting whether a person wears a face mask from speech is useful in modelling speech in forensic investigations, communication between surgeons or people protecting themselves against infectious diseases such as COVID-19. In this paper, we propose a novel data augmentation approach for mask detection from speech. Our approach is based on (i) training Generative Adversarial Networks (GANs) with cycle-consistency loss to translate unpaired utterances between two classes (with mask and without mask), and on (ii) generating new training utterances using the cycle-consistent GANs, assigning opposite labels to each translated utterance. Original and translated utterances are converted into spectrograms which are provided as input to a set of ResNet neural networks with various depths. The networks are combined into an ensemble through a Support Vector Machines (SVM) classifier. With this system, we participated in the Mask Sub-Challenge (MSC) of the INTERSPEECH 2020 Computational Paralinguistics Challenge, surpassing the baseline proposed by the organizers by 2.8%. Our data augmentation technique provided a performance boost of 0.9% on the private test set. Furthermore, we show that our data augmentation approach yields better results than other baseline and state-of-the-art augmentation methods.


title: The Essential Role of Empirical Validation in Legislative Redistricting Simulation

abstract: As granular data about elections and voters become available, redistricting simulation methods are playing an increasingly important role when legislatures adopt redistricting plans and courts determine their legality. These simulation methods are designed to yield a representative sample of all redistricting plans that satisfy statutory guidelines and requirements such as contiguity, population parity, and compactness. A proposed redistricting plan can be considered gerrymandered if it constitutes an outlier relative to this sample according to partisan fairness metrics. Despite their growing use, an insufficient effort has been made to empirically validate the accuracy of the simulation methods. We apply a recently developed computational method that can efficiently enumerate all possible redistricting plans and yield an independent uniform sample from this population. We show that this algorithm scales to a state with a couple of hundred geographical units. Finally, we empirically examine how existing simulation methods perform on realistic validation data sets.


title: Down type iso-singlet quarks at the HL-LHC and FCC-hh

abstract: We study the discovery potential of down type iso-singlet quarks, $D$, predicted by the $E_6$ GUT model in the ${pp\rightarrow D\bar{D}\rightarrow ZZd\bar{d} \rightarrow \ell^+\ell^-\ell^+\ell^- d\bar{d}}$ channel at the HL-LHC and FCC-hh colliders. The analysis is performed using a high level analysis description language and its runtime interpreter. The study shows that, using solely this channel, HL-LHC can discover $D$ quarks up to a mass of 710 GeV whereas FCC-hh up to 2430 GeV with data collected in their complete run periods.


title: Determining The Magnetic Potential In The Fractional Magnetic Calder\'on Problem

abstract: We determine both the magnetic potential and the electric potential from the exterior partial measurements of the Dirichlet-to-Neumann map in the fractional linear magnetic Calder\'on problem by using an integral identity. We also determine both the magnetic potential and the nonlinearity in the fractional semilinear magnetic Calder\'on problem by using a first order linearization.


title: Application of Local Approaches to the Assessment of Fatigue Test results obtained for Welded Joints at Sub-Zero Temperatures

abstract: Several studies have found significant increase in the fatigue strength of welded joints of structural steels at sub-zero temperatures. This study addresses the research by investigating the applicability of local fatigue assessment methods to welded joints exposed to sub-zero temperatures. For this purpose, fatigue test results of two fillet weld details with weld toe and weld root failure are evaluated at a range of temperatures using a variety of structural hot-spot and notch stress approaches, then are compared to the nominal stress approach. Large differences in prediction accuracy are found for the analysed assessment methods and both failure locations.


title: Extremal overlap-free and extremal $\beta$-free binary words

abstract: An overlap-free (or $\beta$-free) word $w$ over a fixed alphabet $\Sigma$ is extremal if every word obtained from $w$ by inserting a single letter from $\Sigma$ at any position contains an overlap (or a factor of exponent at least $\beta$, respectively). We find all lengths which admit an extremal overlap-free binary word. For every extended real number $\beta$ such that $2^+\leq\beta\leq 8/3$, we show that there are arbitrarily long extremal $\beta$-free binary words.


title: Identification of the Majorana edge modes in tight-binding systems based on the Krylov method

abstract: Low dimensional structures in the non-trivial topological phase can host the in-gap Majorana bound states, identified experimentally as zero-bias peaks in differential conductance. Theoretical methods for studying Majorana modes are mostly based on the bulk-boundary correspondence or exact diagonalization of finite systems via, e.g., Bogoliubov-de Gennes formalism. In this paper, we develop an efficient method for identifying the Majorana in-gap (edge) states via looking for extreme eigenvalues of symmetric matrices. The presented approach is based on the Krylov method and allows for study the spatial profile of the modes as well as the spectrum of the system. The advantage of this method is the calculation cost, which shows linear dependence on the number of lattice sites. The latter problem may be solved for very large clusters of arbitrary shape/geometry. In order to demonstrate the efficiency of our approach, we study two- and three-dimensional clusters described by the Kitaev and Rashba models for which we determine the number of Majorana modes and calculate their spatial structures. Additionally, we discuss the impact of the system size on the physical properties of the topological phase of the magnetic nanoisland deposited on the superconducting surface. In this case, we have shown that the eigenvalues of the in-gap states depend on the length of the system edge.


title: Herschel 158$ \mu$m [CII] Observations of "CO-dark" Gas in the Perseus Giant Molecular Cloud

abstract: We present observations of velocity-resolved [CII] 158 $\mu$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integrated intensity of the [CII] emission remains relatively constant across each boundary, despite the observed range in optical extinction between 1 mag and 10 mag. This flat profile indicates a constant heating and cooling rate within both regions observed. The integrated intensity of [CII] emission is reasonably well correlated with the neutral hydrogen (HI) column density, as well as total gas column density. This, in addition to the 80$'$ (7 pc) extent of the [CII] emission from cloud center, suggests that the HI envelope plays a dominant role in explaining the [CII] emission emanating from Perseus. We compare the [CII] and $^{12}$CO integrated intensities with predictions from a 1-D, two-sided slab PDR model and show that a simple core $+$ envelope, equilibrium model without an additional "CO-dark" H$_2$ component can reproduce observations well. Additional observations are needed to disentangle how much of the [CII] emission is associated with the "CO-dark" H$_2$ gas, as well as constrain spatial variations of the dust-to-gas ratio across Perseus.


title: Mapping dark matter and finding filaments: calibration of lensing analysis techniques on simulated data

abstract: We quantify the performance of mass mapping techniques on mock imaging and gravitational lensing data of galaxy clusters. The optimum method depends upon the scientific goal. We assess measurements of clusters' radial density profiles, departures from sphericity, and their filamentary attachment to the cosmic web. We find that mass maps produced by direct (KS93) inversion of shear measurements are unbiased, and that their noise can be suppressed via filtering with MRLens. Forward-fitting techniques, such as Lenstool, suppress noise further, but at a cost of biased ellipticity in the cluster core and over-estimation of mass at large radii. Interestingly, current searches for filaments are noise-limited by the intrinsic shapes of weakly lensed galaxies, rather than by the projection of line-of-sight structures. Therefore, space-based or balloon-based imaging surveys that resolve a high density of lensed galaxies, could soon detect one or two filaments around most clusters.