title: The gate-tunable strong and fragile topology of multilayer-graphene on a transition metal dichalcogenide

abstract: We analyze the phase diagram of multilayer-graphene sandwiched between identical transition metal dichalcogenides. Recently realized in all van-der-Wall heterostructures, these sandwiches induce sizable (1-15 meV) spin orbit coupling in the graphene, offering a way to engineer topological band-structures in a pristine and gate-tunable platform. We find a rich phase diagram that depends on the number of layers $N$ and the gate-tunable perpendicular electric field. For $N > 1$ and odd, the system is a strong 2D topological insulator with a gap equal to the strength of proximity-induced Ising spin-orbit coupling, which reverts to a trivial phase at moderate electric fields. For $N$-even, the low energy bands exhibit a recently proposed form of "fragile" crystalline topology, as well as electric-field tuned symmetry-protected phase transitions between distinct atomic insulators. Hence AB-stacked bilayer and ABC-stacked trilayer graphene are predicted to provide controllable experimental realizations of fragile and strong topology.


title: First stellar spectroscopy in Leo P

abstract: We present the first stellar spectroscopy in the low-luminosity (M_V ~-9.3 mag), dwarf galaxy Leo P. Its significantly low oxygen abundance (3% solar) and relative proximity (~1.6 Mpc) make it a unique galaxy to investigate the properties of massive stars with near-primordial compositions akin to those in the early Universe. From our VLT-MUSE spectroscopy we find the first direct evidence for an O-type star in the prominent HII region, providing an important test case to investigate the potential environmental dependence of the upper end of the initial mass function in the dwarf galaxy regime. We classify 14 further sources as massive stars (and 17 more as candidate massive stars), most likely B-type objects. From comparisons with published evolutionary models we argue that the absolute visual magnitudes of massive stars in very metal-poor systems such as Leo P and I Zw 18 may be fainter by ~0.5 mag compared to Galactic stars. We also present spectroscopy of two carbon stars identified previously as candidate asymptotic-giant-branch stars. Two of three further candidate asymptotic-giant-branch stars display CaII absorption, confirming them as cool, evolved stars; we also recover CaII absorption in the stacked data of the next brightest 16 stars in the upper red giant branch. These discoveries will provide targets for future observations to investigate the physical properties of these objects and to calibrate evolutionary models of luminous stars at such low metallicity. The MUSE data also reveal two 100pc-scale ring structures in Halpha emission, with the HII region located on the northern edge of the southern ring. Lastly, we report serendipitous observations of 20 galaxies, with redshifts ranging from z=0.39, to a close pair of star-forming galaxies at z=2.5.


title: Compensating for Interference in Sliding Window Detection Processes using a Bayesian Paradigm

abstract: Sliding window detectors are non-coherent decision processes, designed in an attempt to control the probability of false alarm, for application to radar target detection. In earlier low resolution radar systems it was possible to specify such detectors quite easily, due to the Gaussian nature of clutter returns, in an X-band maritime surveillance radar context. As radar resolution improved with corresponding developments in modern technology, it became difficult to construct sliding window detectors with the constant false alarm rate property. However, over the last eight years this situation has been rectified, due to improved understanding of the way in which such detectors should be constructed. This paper examines the Bayesian approach to the construction of such detectors. In particular, the design of sliding window detectors, with the constant false alarm rate property, with the capacity to manage interfering targets, will be outlined.


title: Transfer Free Energies and Partitioning of Small Molecules in Collapsed PNIPAM Polymers

abstract: A central quantity in the design of functional hydrogels used as nanocarrier systems, for instance for drug delivery or adaptive nanocatalysis, is the partition ratio, which quantifies the uptake of a molecular substance by the polymer matrix. By employing all-atom molecular dynamics simulations, we study the solvation and partitioning (with respect to bulk water) of small subnanometer-sized solutes in a dense matrix of collapsed Poly(N-isopropylacrylamide) (PNIPAM) polymers above the lower critical solution temperature (LCST) in aqueous solution. We examine the roles of the solute's polarity and its size on the solubility properties in the thermoresponsive polymer. We show that the transfer free energies of nonpolar solutes from bulk water into the polymer are favorable and scale in a good approximation with the solute's surface area. Even for small solute size variation, partitioning can vary over orders of magnitude. A polar nature of the solute, on the other hand, generally opposes the transfer, at least for alkyl solutes. Finally, we find a strong correlation between the transfer free energies in the gel and the adsorption free energies on a single extended polymer chain, which enables us to relate the partition ratios in the swollen and collapsed state of a PNIPAM gel.


title: PELICAN: deeP architecturE for the LIght Curve ANalysis

abstract: We developed a deeP architecturE for the LIght Curve ANalysis (PELICAN) for the characterization and the classification of light curves. It takes light curves as input, without any additional features. PELICAN can deal with the sparsity and the irregular sampling of light curves. It is designed to remove the problem of non-representativeness between the training and test databases coming from the limitations of the spectroscopic follow-up. We applied our methodology on different supernovae light curve databases. First, we evaluated PELICAN on the Supernova Photometric Classification Challenge for which we obtained the best performance ever achieved with a non-representative training database, by reaching an accuracy of 0.811. Then we tested PELICAN on simulated light curves of the LSST Deep Fields for which PELICAN is able to detect 87.4% of supernovae Ia with a precision higher than 98%, by considering a non-representative training database of 2k light curves. PELICAN can be trained on light curves of LSST Deep Fields to classify light curves of LSST main survey, that have a lower sampling rate and are more noisy. In this scenario, it reaches an accuracy of 96.5% with a training database of 2k light curves of the Deep Fields. It constitutes a pivotal result as type Ia supernovae candidates from the main survey might then be used to increase the statistics without additional spectroscopic follow-up. Finally we evaluated PELICAN on real data from the Sloan Digital Sky Survey. PELICAN reaches an accuracy of 86.8% with a training database composed of simulated data and a fraction of 10% of real data. The ability of PELICAN to deal with the different causes of non-representativeness between the training and test databases, and its robustness against survey properties and observational conditions, put it on the forefront of the light curves classification tools for the LSST era.


title: Failure of the $L^1$ pointwise ergodic theorem for $\mathrm{PSL}_2(\mathbb{R})$

abstract: Amos Nevo established the pointwise ergodic theorem in $L^p$ for measure-preserving actions of $\mathrm{PSL}_2(\mathbb{R})$ on probability spaces with respect to ball averages and every $p>1$. This paper shows by explicit example that Nevo's Theorem cannot be extended to $p=1$.


title: Radial velocities in the outermost disk toward the anticenter

abstract: We measure the mean Galactocentric radial component of the velocity of stars ($v_R$) in the disk at 8 kpc$<R<28$ kpc in the direction of the anticenter. For this, we use the Apache Point Galactic Evolution Experiment (APOGEE). Furthermore, we compare the result with HI maps along the same line of sight. We find an increase in positive (expansion) $v_R$ at $R\approx 9-13$ kpc, reaching a maximum of $\approx 6$ km/s, and a decrease at large values of $R$ reaching a negative (contraction) value of $\approx -10$ km/s for $R>17$ kpc. Negative velocities are also observed in 21 cm HI maps, possibly dominated by local gas emission. Among the possible dynamical causes for these non-zero $v_R$, factors such as the effect of the Galactic bar, streams, or mergers do not seem appropriate to explain our observations. An explanation might be the gravitational attraction of overdensities in a spiral arm. As a matter of fact, we see a change of regime from positive to negative velocities around $R\approx 15$ kpc, in the position where we cross the Outer spiral arm in the anticenter. The mass in spiral arms necessary to produce these velocities would be about 3\% of the mass of the disk, consistent with our knowledge of the spiral arms. Another scenario that we explore is a simple class of out-of-equilibrium systems in which radial motions are generally created by the monolithic collapse of isolated self-gravitating overdensities.


title: Second bounded cohomology and WWPD

abstract: Given a group acting on a Gromov hyperbolic space, Bestvina and Fujiwara introduced the WPD property --- weak proper discontinuity --- for studying the 2nd bounded cohomology of the group. We carry out a more general study of second bounded cohomology using a 'really' weak property discontinuity property known as WWPD that was introduced by Bestvina, Bromberg, and Fujiwara.


title: Risk Neutral Reformulation Approach to Risk Averse Stochastic Programming

abstract: The aim of this paper is to show that in some cases risk averse multistage stochastic programming problems can be reformulated in a form of risk neutral setting. This is achieved by a change of the reference probability measure making ``bad" (extreme) scenarios more frequent. As a numerical example we demonstrate advantages of such change-of-measure approach applied to the Brazilian Interconnected Power System operation planning problem.


title: The i3+3 Design for Phase I Clinical Trials

abstract: Purpose: The 3+3 design has been shown to be less likely to achieve the objectives of phase I dose-finding trials when compared with more advanced model-based designs. One major criticism of the 3+3 design is that it is based on simple rules, does not depend on statistical models for inference, and leads to unsafe and unreliable operating characteristics. On the other hand, being rule-based allows 3+3 to be easily understood and implemented in practice, making it the first choice among clinicians. Is it possible to have a rule-based design with great performance? Methods: We propose a new rule-based design called i3+3, where the letter "i" represents the word "interval". The i3+3 design is based on simple but more advanced rules that account for the variabilities in the observed data. We compare the operating characteristics for the proposed i3+3 design with other popular phase I designs by simulation. Results: The i3+3 design is far superior than the 3+3 design in trial safety and the ability to identify the true MTD. Compared with model-based phase I designs, i3+3 also demonstrates comparable performances. In other words, the i3+3 design possesses both the simplicity and transparency of the rule-based approaches, and the superior operating characteristics seen in model-based approaches. An online R Shiny tool (https://i3design.shinyapps.io/i3plus3/) is provided to illustrate the i3+3 design, although in practice it requires no software to design or conduct a dose-finding trial. Conclusion: The i3+3 design could be a practice-altering method for the clinical community.


title: A Private SUSY 4HDM with FCNC in the Up-sector

abstract: We present a SUSY model with four Higgs doublets of the ``private type", where each fermion type (up, down, and charged leptons) obtain their masses from a different Higgs doublet H_f ( f = u_1, d, e). The conditions for anomaly cancellation, imply that the remaining Higgs doublet of the model (H_{u_2}), must have the same hypercharge as H_{u_1}, and thus can only couple to up-type quarks, which opens the possibility to have FCNC's only in this sector. We study the Lagrangian of the model, and in particular the Higgs potential, in order to identify the Higgs mass-eigenstates and their interactions; for the Yukawa matrices we consider the 4-texture case. We obtain constraints on the model parameters by using LHC measurements on the properties of the 125 GeV Higgs boson ($h$), and identify viable regions of parameter space. These constraints are then used to evaluate the prospects to detect the FCNC decay mode t -> ch at the future high-luminosity (HL) option for the LHC, which are compared with current limits from LHC-run2. Then, we also evaluate the FCNC decay of the next heavier Higgs boson H_2 -> tc, which can reach typically BR(H_2 -> tc) = O(10^{-4}-10^{-5} ); the search for the signal at HL-LHC is also studied, finding that it could be detectable for some specific regions of the model parameter space.


title: Secular period variations of RR Lyrae stars in NGC 4147

abstract: A compilation of photometric time-series data of NGC 4147 spanning 58 years, enabled to study the secular behaviour of the pulsation period in the population of RR Lyrae stars in the cluster. The traditional O-C diagram approach was employed. Three stars display a significant period decrease and two a comparable period increase. Nevertheless, the median value of the whole sample is nearly zero, which seems consistent with theoretical expectations. The limitations of empirical determinations of secular period changes in globular clusters to corroborate evolutionary theoretical predictions are discussed.


title: The effect of horizontal plasma inhomogeneities in 3D NLTE radiation transfer in stellar atmospheres

abstract: {We aim to demonstrate the effect of atmospheric inhomogeneities on the emergent specific intensity and radiation flux of a spectral line radiation.} {We self-consistently solve the NLTE problem for a two-level atom in a 3D atmosphere using the Cartesian grid. For that purpose, we use the 3D radiative transfer code PORTA. By examining simple examples, we study cases where the temperature inhomogeneities in the atmosphere models lead to the modification of the emergent radiation.} {We show that specific temperature inhomogeneities in the model atmospheres influence the emerging radiation, and that interpretation of the stellar spectra based on a plane-parallel atmosphere models can lead to erroneous conclusions about the atmospheric structure.}


title: Modeling of Radiation-Induced Defect Recovery in 3C-SiC Under High Field Bias Conditions

abstract: In this work, the implications of high field bias conditions in radiation-induced defect recovery in 3C-SiC crystals is studied. It is well known that transient heating effects (or thermal spikes) occur when energetic swift heavy ions (SHIs) deposit energy to the surrounding medium via ionization. Here, we explore the dynamics of this transient event under high background electric fields in 3C-SiC, which is what occurs when an ion strike coincides with field-sensitive volumes. In this study, we use the Ensemble Monte Carlo method to quantify how the energy deposition of the ionized regions change in response to high background electric fields. Subsequently, we study the relationship between the radiation-induced thermal spike and defect recovery using molecular dynamics simulations. We find that field strengths below the critical breakdown of wide bandgap devices are sufficient to exacerbate the localized heating, which subsequently enhances the defect recovery effect. This work is beneficial for 3C-SiC electronics and materials used in high radiation environments.


title: Integrating Phase 2 into Phase 3 based on an Intermediate Endpoint While Accounting for a Cure Proportion -- with an Application to the Design of a Clinical Trial in Acute Myeloid Leukemia

abstract: For a trial with primary endpoint overall survival for a molecule with curative potential, statistical methods that rely on the proportional hazards assumption may underestimate the power and the time to final analysis. We show how a cure proportion model can be used to get the necessary number of events and appropriate timing via simulation. If Phase 1 results for the new drug are exceptional and/or the medical need in the target population is high, a Phase 3 trial might be initiated after Phase 1. Building in a futility interim analysis into such a pivotal trial may mitigate the uncertainty of moving directly to Phase 3. However, if cure is possible, overall survival might not be mature enough at the interim to support a futility decision. We propose to base this decision on an intermediate endpoint that is sufficiently associated with survival. Planning for such an interim can be interpreted as making a randomized Phase 2 trial a part of the pivotal trial: if stopped at the interim, the trial data would be analyzed and a decision on a subsequent Phase 3 trial would be made. If the trial continues at the interim then the Phase 3 trial is already underway. To select a futility boundary, a mechanistic simulation model that connects the intermediate endpoint and survival is proposed. We illustrate how this approach was used to design a pivotal randomized trial in acute myeloid leukemia, discuss historical data that informed the simulation model, and operational challenges when implementing it.


title: Adaptive Internal Models: Explaining the Oculomotor System and the Cerebellum

abstract: We propose a new model of the oculomotor system, particularly the vestibulo-ocular reflex, gaze fixation, and smooth pursuit. Our key insight is to exploit recent developments on adaptive internal models. The outcome is a simple model that includes the interactions between the brainstem and the cerebellum and that recovers behaviors from more than 15 oculomotor experiments. In addition, we put forward a thesis that the cerebellum embodies internal models of all persistent, exogenous reference and disturbance signals acting on the body and observable through the error signals it receives.


title: Oxygen-induced in-situ manipulation of the interlayer coupling and exciton recombination in Bi2Se3/MoS2 2D heterostructures

abstract: 2D heterostructures are more than a sum of the parent 2D materials, but are also a product of the interlayer coupling, which can induce new properties. In this paper we present a method to tune the interlayer coupling in Bi2Se3/MoS2 2D heterostructures by regulating the oxygen presence in the atmosphere, while applying laser or thermal energy. Our data suggests the interlayer coupling is tuned through the diffusive intercalation and de-intercalation of oxygen molecules. When one layer of Bi2Se3 is grown on monolayer MoS2, an influential interlayer coupling is formed that quenches the signature photoluminescence (PL) peaks. However, thermally annealing in the presence of oxygen disrupts the interlayer coupling, facilitating the emergence of the MoS2 PL peak. DFT calculations predict intercalated oxygen increases the interlayer separation ~17%, disrupting the interlayer coupling and inducing the layers to behave more electronically independent. The interlayer coupling can then be restored by thermally annealing in N2 or Ar, where the peaks will re-quench. Hence, this is an interesting oxygen-induced switching between "non-radiative" and "radiative" exciton recombination. This switching can also be accomplished locally, controllably, and reversibly using a low-power focused laser, while changing the environment from pure N2 to air. This allows for the interlayer coupling to be precisely manipulated with submicron spatial resolution, facilitating site-programmable 2D light-emitting pixels whose emission intensity could be precisely varied by a factor exceeding 200x. Our results show that these atomically-thin 2D heterostructures may be excellent candidates for oxygen sensing.


title: Cosmological Implications of the Generalized Entropy Based Holographic Dark Energy Models in Dynamical Chern-Simons Modified Gravity

abstract: Recently, Tsallis, R\'{e}nyi and Sharma-Mitall and entropies have widely been used to study the gravitational and cosmological setups. We consider a flat FRW universe with linear interaction between dark energy and dark matter. We discuss the dark energy models using Tsallis, R\'{e}nyi and Sharma-Mitall entropies in the framework of Chern-Simons modified gravity. We explore various cosmological parameters (equation of state parameter, squared sound of speed ) and cosmological plane ($\omega_{d}-\omega_{d}'$, where $\omega_{d}'$, is the evolutionary equation of state parameter). It is observed that the equation of state parameter gives quintessence-like nature of the universe in most of the cases. Also, the squared speed of sound shows stability of the models for Tsallis, R\'{e}nyi dark energy model while unstable behavior for Sharma-Mitall dark energy model. The $\omega_{d}-\omega_{d}'$ plane represents the thawing region for all dark energy models.


title: Braided fermions from Hurwitz algebras

abstract: Some curious structural similarities between a recent braid- and Hurwitz algebraic description of the unbroken internal symmetries for a single generations of Standard Model fermions were recently identified. The non-trivial braid groups that can be represented using the four normed division algebras are $B_2$ and $B_3^c$, exactly those required to represent a single generation of fermions in terms of simple three strand ribbon braids. These braided fermion states can be identified with the basis states of the minimal left ideals of the Clifford algebra $C\ell(6)$, generated from the nested left actions of the complex octonions $\mathbb{C}\otimes\mathbb{O}$ on itself. That is, the ribbon spectrum can be related to octonion algebras. Some speculative ideas relating to ongoing research that attempts to construct a unified theory based on braid groups and Hurwitz algebras are discussed.


title: Steinberg groups for Jordan pairs - an introduction with open problems

abstract: The paper gives an introduction to Steinberg groups for root graded Jordan pairs, a theory developed in the recent book "Steinberg groups for Jordan pairs" by Ottmar Loos and the author.


title: Fractional Burgers models in creep and stress relaxation tests

abstract: Classical and thermodynamically consistent fractional Burgers models are examined in creep and stress relaxation tests. Using the Laplace transform method, the creep compliance and relaxation modulus are obtained in integral form, that yielded, when compared to the thermodynamical requirements, the narrower range of model parameters in which the creep compliance is a Bernstein function while the relaxation modulus is completely monotonic. Moreover, the relaxation modulus may even be oscillatory function with decreasing amplitude. The asymptotic analysis of the creep compliance and relaxation modulus is performed near the initial time-instant and for large time as well.


title: Leading higher-derivative corrections to Kerr geometry

abstract: We compute the most general leading-order correction to Kerr solution when the Einstein-Hilbert action is supplemented with higher-derivative terms, including the possibility of dynamical couplings controlled by scalars. The model we present depends on five parameters and it contains, as particular cases, Einstein-dilaton-Gauss-Bonnet gravity, dynamical Chern-Simons gravity and the effective action coming from Heterotic Superstring theory. By solving the corrected field equations, we find the modified Kerr metric that describes rotating black holes in these theories. We express the solution as a series in the spin parameter $\chi$, and we show that including enough terms in the expansion we are able to describe black holes with large spin. For the computations in the text we use an expansion up to order $\chi^{14}$, which is accurate for $\chi<0.7$, but we provide as well a Mathematica notebook that computes the solution at any given order. We study several properties of the corrected black holes, such as geometry of the horizon, ergosphere, light rings and scalar hair. Some of the corrections violate parity, and we highlight in those cases plots of horizons and ergospheres without $\mathbb{Z}_{2}$ symmetry.


title: On the strong convergence of partial sums with respect to bounded Vilenkin systems

abstract: In this paper we investigate some strong convergence theorems for partial sums with respect to Vilenkin system.


title: Dark energy as an emergent phenomenon

abstract: The cosmological constant, which was introduced by Einstein a century ago to allow for a static universe, experienced a revival two decades ago under the label dark energy as a parameter to model the observed accelerated expansion of the universe. Its physical nature has however remained enigmatic. Here we use the Einstein equations without cosmological constant to show that the origin of the accelerated expansion is not in the equations but has to do with the boundary conditions related to the causal horizon, which exists because the age of the universe is finite. Via transformations to conformal coordinates and Euclidian spacetime we find a resonance condition that uniquely determines the dimensionless parameter $\Omega_\Lambda$ that governs the observed cosmic acceleration: $\Omega_\Lambda =\textstyle{2\over 3}\,(\pi\, t_H/t_c)^2$, where $t_H$ is the Hubble time and $t_c$ is the conformal age of the universe. This explanation leads to a somewhat modified cosmology, in which the expansion rate of the early universe is 2.1 times faster than in the standard model. We show that Big Bang nucleosynthesis calculations with the faster expansion rate requires the mean baryon density to be raised to the level of the total matter density to agree with the observed deuterium abundance. This appears to eliminate the need to invoke the existence of some yet to be discovered exotic particles to explain dark matter, since all of it may be baryonic while still remaining consistent with the observed abundances of the light elements.


title: Generalized Kasha's Scheme for Classifying Two-Dimensional Excitonic Molecular Aggregates: Temperature Dependent Absorption Peak Frequency Shift

abstract: We propose a generalized theoretical framework for classifying two-dimensional (2D) excitonic molecular aggregates based on an analysis of temperature dependent spectra. In addition to the monomer-aggregate absorption peak shift, which defines the conventional J- and H-aggregates, we incorporate the peak shift associated with increasing temperature as a measure to characterize the exciton band structure. First we show that there is a one-to-one correspondence between the monomer-aggregate and the T-dependent peak shifts for Kasha's well-established model of 1D aggregates, where J-aggregates exhibit further redshift upon increasing temperature and H-aggregates exhibit further blueshift. On the contrary, 2D aggregate structures are capable of supporting the two other combinations: blueshifting J-aggregates and redshifting H-aggregates, owing to their more complex exciton band structures. Secondly, using spectral lineshape theory, the T-dependent shift is associated with the relative abundance of states on each side of the bright state. We further establish that the density of states can be connected to the microscopic packing condition leading to these four classes of aggregates by separately considering the short and long-range contribution to the excitonic couplings. In particular the T-dependent shift is shown to be an unambiguous signature for the sign of net short-range couplings: Aggregates with net negative (positive) short-range couplings redshift (blueshift) with increasing temperature. Lastly, comparison with experiments shows that our theory can be utilized to quantitatively account for the observed but previously unexplained T-dependent absorption lineshapes. Thus, our work provides a firm ground for elucidating the structure-function relationships for molecular aggregates and is fully compatible with existing experimental and theoretical structure characterization tools.


title: Shrinking the Quark Gluon Plasma

abstract: In recent years the understanding on the limits of the smallest possible droplet of the Quark Gluon Plasma has been called into question. Experimental results from both the Large Hadron Collider and the Relativistic Heavy Ion Collider have provided hints that the Quark Gluon Plasma may be produced in systems as small as that formed in pPb or dAu collisions. Yet alternative explanations still exist from correlations arising from quarks and gluons in a color glass condensate picture. In order to resolve these two scenarios, a system size scan has been proposed at the Large Hadron Collider for collisions of ArAr and OO. Here we make predictions for a possible future run of ArAr and OO collisions at the Large Hadron Collider and study the system size dependence of a variety of flow observables. We find that linear response (from the initial conditions to the final flow harmonics) becomes more dominant in smaller systems whereas linear+cubic response can accurately predict multi-particle cumulants for a wide range of centralities in large systems.


title: Structure-Property Relationship in Layered BaMn2Sb2 and Ba2Mn3Sb2O2

abstract: Layered transition-metal compounds have received great attention owing to their novel physical properties. Here, we present the structural, electronic, thermal, and magnetic properties of BaMn2Sb2 and Ba2Mn3Sb2O2 single crystals, both with the layered structure analogous to high-temperature superconductors. While the Mn moment in the MnSb4 tetrahedral environment forms G-type antiferromagnetic (AFM) ordering in both BaMn2Sb2 (TN1~443 K) and Ba2Mn3Sb2O2 (TN1~314 K), a short-range AFM order is found in the intercalated MnO2 layer at a much lower temperature (TN2~60 K) in Ba2Mn3Sb2O2. The directions of the ordered moments in these two magnetic sub-lattices of Ba2Mn3Sb2O2 are perpendicular to each other, even though the system is electrically conductive. This indicates that the large magnetic moments in these compounds are highly localized, leading to negligible coupling between MnSb4 and MnO2 layers in Ba2Mn3Sb2O2. These findings provide an insight into the structure-magnetism-based design principle for new superconductors.


title: Diverging exchange force and form of the exact density matrix functional

abstract: For translationally invariant one-band lattice models, we exploit the ab initio knowledge of the natural orbitals to simplify reduced density matrix functional theory (RDMFT). Striking underlying features are discovered: First, within each symmetry sector, the interaction functional $\mathcal{F}$ depends only on the natural occupation numbers $\bf{n}$. The respective sets $\mathcal{P}^1_N$ and $\mathcal{E}^1_N$ of pure and ensemble $N$-representable one-matrices coincide. Second, and most importantly, the exact functional is strongly shaped by the geometry of the polytope $\mathcal{E}^1_N \equiv \mathcal{P}^1_N $, described by linear constraints $D^{(j)}(\bf{n})\geq 0$. For smaller systems, it follows as $\mathcal{F}[\bf{n}]=\sum_{i,i'} \overline{V}_{i,i'} \sqrt{D^{(i)}(\bf{n})D^{(i')}(\bf{n})}$. This generalizes to systems of arbitrary size by replacing each $D^{(i)}$ by a linear combination of $\{D^{(j)}(\bf{n})\}$ and adding a non-analytical term involving the interaction $\hat{V}$. Third, the gradient $\mathrm{d}\mathcal{F}/\mathrm{d}\bf{n}$ is shown to diverge on the boundary $\partial\mathcal{E}^1_N$, suggesting that the fermionic exchange symmetry manifests itself within RDMFT in the form of an "exchange force". All findings hold for systems with non-fixed particle number as well and $\hat{V}$ can be any $p$-particle interaction. As an illustration, we derive the exact functional for the Hubbard square.


title: Transformed Snapshot Interpolation with High Resolution Transforms

abstract: In the last few years, several methods have been developed to deal with jump singularities in parametric or stochastic hyperbolic PDEs. They typically use some alignment of the jump-sets in physical space before performing well established reduced order modelling techniques such as reduced basis methods, POD or simply interpolation. In the current literature, the transforms are typically of low resolution in space, mostly low order polynomials, Fourier modes or constant shifts. In this paper, we discuss higher resolution transforms in one of the recent methods, the transformed snapshot interpolation (TSI). We introduce a new discretization of the transforms with an appropriate behaviour near singularities and consider their numerical computation via an optimization procedure.


title: Rotational Light Curves of Jupiter from UV to Mid-Infrared and Implications for Brown Dwarfs and Exoplanets

abstract: Rotational modulations are observed on brown dwarfs and directly imaged exoplanets, but the underlying mechanism is not well understood. Here, we analyze Jupiter's rotational light curves at 12 wavelengths from the ultraviolet (UV) to the mid-infrared (mid-IR). Peak-to-peak amplitudes of Jupiter's light curves range from sub percent to 4% at most wavelengths, but the amplitude exceeds 20% at 5 $\rm \mu m$, a wavelength sensing Jupiter's deep troposphere. Jupiter's rotational modulations are primarily caused by discrete patterns in the cloudless belts instead of the cloudy zones. The light-curve amplitude is dominated by the sizes and brightness contrasts of the Great Red Spot (GRS), expansions of the North Equatorial Belt (NEB), patchy clouds in the North Temperate Belt (NTB) and a train of hot spots in the NEB. In reflection, the contrast is controlled by upper tropospheric and stratospheric hazes, clouds, and chromophores in the clouds. In thermal emission, the small rotational variability is caused by the spatial distribution of temperature and opacities of gas and aerosols; the large variation is caused by the $\rm NH_{3}$ cloud holes and thin-thick clouds. The methane-band light curves exhibit opposite-shape behavior compared with the UV and visible wavelengths, caused by wavelength-dependent brightness change of the GRS. Light-curve evolution is induced by periodic events in the belts and longitudinal drifting of the GRS and patchy clouds in the NTB. This study suggests several interesting mechanisms related to distributions of temperature, gas, hazes, and clouds for understanding the observed rotational modulations on brown dwarfs and exoplanets.


title: A comment on ergodic theorem for amenable groups

abstract: We prove a version of ergodic theorem for an action of an amenable group, where a F{\o} lner sequence needs not to be tempered. Instead, it is assumed that a function satisfies certain mixing condition.


title: Optimal Decision-Making in Mixed-Agent Partially Observable Stochastic Environments via Reinforcement Learning

abstract: Optimal decision making with limited or no information in stochastic environments where multiple agents interact is a challenging topic in the realm of artificial intelligence. Reinforcement learning (RL) is a popular approach for arriving at optimal strategies by predicating stimuli, such as the reward for following a strategy, on experience. RL is heavily explored in the single-agent context, but is a nascent concept in multiagent problems. To this end, I propose several principled model-free and partially model-based reinforcement learning approaches for several multiagent settings. In the realm of normative reinforcement learning, I introduce scalable extensions to Monte Carlo exploring starts for partially observable Markov Decision Processes (POMDP), dubbed MCES-P, where I expand the theory and algorithm to the multiagent setting. I first examine MCES-P with probably approximately correct (PAC) bounds in the context of multiagent setting, showing MCESP+PAC holds in the presence of other agents. I then propose a more sample-efficient methodology for antagonistic settings, MCESIP+PAC. For cooperative settings, I extend MCES-P to the Multiagent POMDP, dubbed MCESMP+PAC. I then explore the use of reinforcement learning as a methodology in searching for optima in realistic and latent model environments. First, I explore a parameterized Q-learning approach in modeling humans learning to reason in an uncertain, multiagent environment. Next, I propose an implementation of MCES-P, along with image segmentation, to create an adaptive team-based reinforcement learning technique to positively identify the presence of phenotypically-expressed water and pathogen stress in crop fields.


title: Electrochemical kinetics of SEI growth on carbon black, II: Modeling

abstract: Mathematical models of capacity fade can reduce the time and cost of lithium-ion battery development and deployment, and growth of the solid-electrolyte interphase (SEI) is a major source of capacity fade. Experiments in Part I reveal nonlinear voltage dependence and strong charge-discharge asymmetry in SEI growth on carbon black negative electrodes, which is not captured by previous models. Here, we present a theoretical model for the electrochemical kinetics of SEI growth coupled to lithium intercalation, which accurately predicts experimental results with few adjustable parameters. The key hypothesis is that the initial SEI is a mixed ion-electron conductor, and its electronic conductivity varies approximately with the square of the local lithium concentration, consistent with hopping conduction of electrons along percolating networks. By including a lithium-ion concentration dependence for the electronic conductivity in the SEI, the bulk SEI thus modulates the overpotential and exchange current of the electrolyte reduction reaction. As a result, SEI growth is promoted during lithiation but suppressed during delithiation. This new insight establishes the fundamental electrochemistry of SEI growth kinetics. Our model improves upon existing models by introducing the effects of electrochemical SEI growth and its dependence on potential, current magnitude, and current direction in predicting capacity fade.


title: How Reliable is the Crowdsourced Knowledge of Security Implementation?

abstract: Stack Overflow (SO) is the most popular online Q&A site for developers to share their expertise in solving programming issues. Given multiple answers to certain questions, developers may take the accepted answer, the answer from a person with high reputation, or the one frequently suggested. However, researchers recently observed exploitable security vulnerabilities in popular SO answers. This observation inspires us to explore the following questions: How much can we trust the security implementation suggestions on SO? If suggested answers are vulnerable, can developers rely on the community's dynamics to infer the vulnerability and identify a secure counterpart? To answer these highly important questions, we conducted a study on SO posts by contrasting secure and insecure advices with the community-given content evaluation. We investigated whether SO incentive mechanism is effective in improving security properties of distributed code examples. Moreover, we also traced duplicated answers to assess whether the community behavior facilitates propagation of secure and insecure code suggestions. We compiled 953 different groups of similar security-related code examples and labeled their security, identifying 785 secure answer posts and 644 insecure ones. Compared with secure suggestions, insecure ones had higher view counts (36,508 vs. 18,713), received a higher score (14 vs. 5), and had significantly more duplicates (3.8 vs. 3.0) on average. 34% of the posts provided by highly reputable so-called trusted users were insecure. Our findings show that there are lots of insecure snippets on SO, while the community-given feedback does not allow differentiating secure from insecure choices. Moreover, the reputation mechanism fails in indicating trustworthy users with respect to security questions, ultimately leaving other users wandering around alone in a software security minefield.


title: StreamBox-HBM: Stream Analytics on High Bandwidth Hybrid Memory

abstract: Stream analytics have an insatiable demand for memory and performance. Emerging hybrid memories combine commodity DDR4 DRAM with 3D-stacked High Bandwidth Memory (HBM) DRAM to meet such demands. However, achieving this promise is challenging because (1) HBM is capacity-limited and (2) HBM boosts performance best for sequential access and high parallelism workloads. At first glance, stream analytics appear a particularly poor match for HBM because they have high capacity demands and data grouping operations, their most demanding computations, use random access. This paper presents the design and implementation of StreamBox-HBM, a stream analytics engine that exploits hybrid memories to achieve scalable high performance. StreamBox-HBM performs data grouping with sequential access sorting algorithms in HBM, in contrast to random access hashing algorithms commonly used in DRAM. StreamBox-HBM solely uses HBM to store Key Pointer Array (KPA) data structures that contain only partial records (keys and pointers to full records) for grouping operations. It dynamically creates and manages prodigious data and pipeline parallelism, choosing when to allocate KPAs in HBM. It dynamically optimizes for both the high bandwidth and limited capacity of HBM, and the limited bandwidth and high capacity of standard DRAM. StreamBox-HBM achieves 110 million records per second and 238 GB/s memory bandwidth while effectively utilizing all 64 cores of Intel's Knights Landing, a commercial server with hybrid memory. It outperforms stream engines with sequential access algorithms without KPAs by 7x and stream engines with random access algorithms by an order of magnitude in throughput. To the best of our knowledge, StreamBox-HBM is the first stream engine optimized for hybrid memories.


title: Exploding Markov operators

abstract: A special class of doubly stochastic (Markov) operators is constructed. These operators come from measure preserving transformations and inherit some of their properties, namely ergodicity and positivity of entropy, yet they may have no pointwise factors.


title: Spatiotemporal Calibration of Atmospheric Nitrogen Dioxide Concentration Estimates From an Air Quality Model for Connecticut

abstract: A spatiotemporal calibration and resolution refinement model was fitted to calibrate nitrogen dioxide (NO$_2$) concentration estimates from the Community Multiscale Air Quality (CMAQ) model, using two sources of observed data on NO$_2$ that differed in their spatial and temporal resolutions. To refine the spatial resolution of the CMAQ model estimates, we leveraged information using additional local covariates including total traffic volume within 2 km, population density, elevation, and land use characteristics. Predictions from this model greatly improved the bias in the CMAQ estimates, as observed by the much lower mean squared error (MSE) at the NO$_2$ monitor sites. The final model was used to predict the daily concentration of ambient NO$_2$ over the entire state of Connecticut on a grid with pixels of size 300 x 300 m. A comparison of the prediction map with a similar map for the CMAQ estimates showed marked improvement in the spatial resolution. The effect of local covariates was evident in the finer spatial resolution map, where the contribution of traffic on major highways to ambient NO$_2$ concentration stands out. An animation was also provided to show the change in the concentration of ambient NO$_2$ over space and time for 1994 and 1995.


title: The ISTI Rapid Response on Exploring Cloud Computing 2018

abstract: This report describes eighteen projects that explored how commercial cloud computing services can be utilized for scientific computation at national laboratories. These demonstrations ranged from deploying proprietary software in a cloud environment to leveraging established cloud-based analytics workflows for processing scientific datasets. By and large, the projects were successful and collectively they suggest that cloud computing can be a valuable computational resource for scientific computation at national laboratories.


title: Spin-orbit driven band inversion in bilayer graphene by van der Waals proximity effect

abstract: Spin orbit coupling (SOC) is the key to realizing time-reversal invariant topological phases of matter. Famously, SOC was predicted by Kane and Mele to stabilize a quantum spin Hall insulator; however, the weak intrinsic SOC in monolayer graphene has precluded experimental observation. Here, we exploit a layer-selective proximity effect---achieved via van der Waals contact to a semiconducting transition metal dichalcogenide--to engineer Kane-Mele SOC in ultra-clean \textit{bilayer} graphene. Using high-resolution capacitance measurements to probe the bulk electronic compressibility, we find that SOC leads to the formation of a distinct incompressible, gapped phase at charge neutrality. The experimental data agrees quantitatively with a simple theoretical model in which the new phase results from SOC-driven band inversion. In contrast to Kane-Mele SOC in monolayer graphene, the inverted phase is not expected to be a time reversal invariant topological insulator, despite being separated from conventional band insulators by electric field tuned phase transitions where crystal symmetry mandates that the bulk gap must close. Electrical transport measurements, conspicuously, reveal that the inverted phase has a conductivity $\sim e^2/h$, which is suppressed by exceptionally small in-plane magnetic fields. The high conductivity and anomalous magnetoresistance are consistent with theoretical models that predict helical edge states within the inversted phase, that are protected from backscattering by an emergent spin symmetry that remains robust even for large Rashba SOC. Our results pave the way for proximity engineering of strong topological insulators as well as correlated quantum phases in the strong spin-orbit regime in graphene heterostructures.


title: Site Specific Knight Shift Measurements of the Dilute Kondo lattice System Ce$_{1-x}$La$_x$CoIn$_5$

abstract: $^{115}$In Nuclear magnetic resonance data are presented for a series of Ce$_{1-x}$La$_x$CoIn$_5$ crystals with different La dilutions, $x$. Multiple In(1) sites associated with different numbers of nearest-neighbor cerium atoms exhibit different Knight shifts and spin lattice relaxation rates. Analysis of the temperature dependence of these sites reveals both an evolution of the heavy electron coherence as a function of dilution, as well as spatial inhomogeneity associated with a complete suppression of antiferromagnetic fluctuations in the vicinity of the La sites. Quantum critical fluctuations persist within disconnected Ce clusters with dilution levels up to 75%, despite the fact that specific heat shows Fermi liquid behavior in dilute samples.


title: An Adaptive Weighted Deep Forest Classifier

abstract: A modification of the confidence screening mechanism based on adaptive weighing of every training instance at each cascade level of the Deep Forest is proposed. The idea underlying the modification is very simple and stems from the confidence screening mechanism idea proposed by Pang et al. to simplify the Deep Forest classifier by means of updating the training set at each level in accordance with the classification accuracy of every training instance. However, if the confidence screening mechanism just removes instances from training and testing processes, then the proposed modification is more flexible and assigns weights by taking into account the classification accuracy. The modification is similar to the AdaBoost to some extent. Numerical experiments illustrate good performance of the proposed modification in comparison with the original Deep Forest proposed by Zhou and Feng.


title: Stochastic phase field $\alpha$-Navier-Stokes vesicle-fluid interaction model

abstract: We consider a stochastic perturbation of the phase field alpha-Navier-Stokes model with vesicle-fluid interaction. It consists in a system of nonlinear evolution partial differential equations modeling the fluid-structure interaction associated to the dynamics of an elastic vesicle immersed in a moving incompressible viscous fluid. This system of equations couples a phase-field equation -- for the interface between the fluid and the vesicle -- to the alpha-Navier-Stokes equation -- for the viscous fluid -- with an extra nonlinear interaction term, namely the bending energy. The stochastic perturbation is an additive space-time noise of trace class on each equation of the system. We prove the existence and uniqueness of solution in classical spaces of $L^{2}$ functions with estimates of non-linear terms and bending energy. It is based on a priori estimate about the regularity of solutions of finite dimensional systems, and tightness of the approximated solution.


title: Projective Decomposition and Matrix Equivalence up to Scale

abstract: A data matrix may be seen simply as a means of organizing observations into rows ( e.g., by measured object) and into columns ( e.g., by measured variable) so that the observations can be analyzed with mathematical tools. As a mathematical object, a matrix defines a linear mapping between points representing weighted combinations of its rows (the row vector space) and points representing weighted combinations of its columns (the column vector space). From this perspective, a data matrix defines a relationship between the information that labels its rows and the information that labels its columns, and numerical methods are used to analyze this relationship. A first step is to normalize the data, transforming each observation from scales convenient for measurement to a common scale, on which addition and multiplication can meaningfully combine the different observations. For example, z-transformation rescales every variable to the same scale, standardized variation from an expected value, but ignores scale differences between measured objects. Here we develop the concepts and properties of projective decomposition, which applies the same normalization strategy to both rows and columns by separating the matrix into row- and column-scaling factors and a scale-normalized matrix. We show that different scalings of the same scale-normalized matrix form an equivalence class, and call the scale-normalized, canonical member of the class its scale-invariant form that preserves all pairwise relative ratios. Projective decomposition therefore provides a means of normalizing the broad class of ratio-scale data, in which relative ratios are of primary interest, onto a common scale without altering the ratios of interest, and simultaneously accounting for scale effects for both organizations of the matrix values. Both of these properties distinguish it from z-transformation.


title: BitCracker: BitLocker meets GPUs

abstract: BitLocker is a full-disk encryption feature available in recent Windows versions. It is designed to protect data by providing encryption for entire volumes and it makes use of a number of different authentication methods. In this paper we present a solution, named BitCracker, to attempt the decryption, by means of a dictionary attack, of memory units encrypted by BitLocker with a user supplied password or the recovery password. To that purpose, we resort to GPU (Graphics Processing Units) that are, by now, widely used as general-purpose coprocessors in high performance computing applications. BitLocker decryption process requires the computation of a very large number of SHA- 256 hashes and also AES, so we propose a very fast solution, highly tuned for Nvidia GPU, for both of them. We analyze the performance of our CUDA implementation on several Nvidia GPUs and we carry out a comparison of our SHA-256 hash with the Hashcat password cracker tool. Finally, we present our OpenCL version, recently released as a plugin of the John The Ripper tool.


title: Distributions of Human Exposure to Ozone During Commuting Hours in Connecticut using the Cellular Device Network

abstract: Epidemiologic studies have established associations between various air pollutants and adverse health outcomes for adults and children. Due to high costs of monitoring air pollutant concentrations for subjects enrolled in a study, statisticians predict exposure concentrations from spatial models that are developed using concentrations monitored at a few sites. In the absence of detailed information on when and where subjects move during the study window, researchers typically assume that the subjects spend their entire day at home, school or work. This assumption can potentially lead to large exposure assignment bias. In this study, we aim to determine the distribution of the exposure assignment bias for an air pollutant (ozone) when subjects are assumed to be static as compared to accounting for individual mobility. To achieve this goal, we use cell-phone mobility data on approximately 400,000 users in the state of Connecticut during a week in July, 2016, in conjunction with an ozone pollution model, and compare individual ozone exposure assuming static versus mobile scenarios. Our results show that exposure models not taking mobility into account often provide poor estimates of individuals commuting into and out of urban areas: the average 8-hour maximum difference between these estimates can exceed 80 parts per billion (ppb). However, for most of the population, the difference in exposure assignment between the two models is small, thereby validating many current epidemiologic studies focusing on exposure to ozone.


title: A function algebra providing new Mergelyan type theorems in several complex variables

abstract: For compact sets $K\subset \mathbb C^{d}$, we introduce a subalgebra $A_{D}(K)$ of $A(K)$, which allows us to obtain Mergelyan type theorems for products of planar compact sets as well as for graphs of functions.


title: File System in Data-Centric Computing

abstract: The moving computation on the edge or near to data is the new trend that can break the bandwidth wall and to unleash the power of next generation NVM or SCM memory. File system is the important OS subsystem that plays the role of mediator between the user-space application and storage device. The key goal of the file system is to represent the file abstraction and to build the files' namespace. In the current paradigm the file system needs to copy the metadata and user data in the DRAM of the host with the goal to access and to modify the user data on the host side. The DAX approach doesn't change the concept but to build the way to bypass the page cache via the direct access to file's content in persistent memory. Generally speaking, for the case of data-centric computing, the file system needs to solve the opposite task not to copy data into page cache but to deliver the processing activity near data on the storage device side.


title: Sheaves: A Topological Approach to Big Data

abstract: This document develops general concepts useful for extracting knowledge embedded in large graphs or datasets that have pair-wise relationships, such as cause-effect-type relations. Almost no underlying assumptions are made, other than that the data can be presented in terms of pair-wise relationships between objects/events. This assumption is used to mine for patterns in the dataset, defining a reduced graph or dataset that boils-down or concentrates information into a more compact form. The resulting extracted structure or set of patterns are manifestly symbolic in nature, as they capture and encode the graph structure of the dataset in terms of a (generative) grammar. This structure is identified as having the formal mathematical structure of a sheaf. In essence, this paper introduces the basic concepts of sheaf theory into the domain of graphical datasets.


title: AVA-ActiveSpeaker: An Audio-Visual Dataset for Active Speaker Detection

abstract: Active speaker detection is an important component in video analysis algorithms for applications such as speaker diarization, video re-targeting for meetings, speech enhancement, and human-robot interaction. The absence of a large, carefully labeled audio-visual dataset for this task has constrained algorithm evaluations with respect to data diversity, environments, and accuracy. This has made comparisons and improvements difficult. In this paper, we present the AVA Active Speaker detection dataset (AVA-ActiveSpeaker) that will be released publicly to facilitate algorithm development and enable comparisons. The dataset contains temporally labeled face tracks in video, where each face instance is labeled as speaking or not, and whether the speech is audible. This dataset contains about 3.65 million human labeled frames or about 38.5 hours of face tracks, and the corresponding audio. We also present a new audio-visual approach for active speaker detection, and analyze its performance, demonstrating both its strength and the contributions of the dataset.


title: Graph Neural Networks with convolutional ARMA filters

abstract: Popular graph neural networks implement convolution operations on graphs based on polynomial spectral filters. In this paper, we propose a novel graph convolutional layer inspired by the auto-regressive moving average (ARMA) filter that, compared to polynomial ones, provides a more flexible frequency response, is more robust to noise, and better captures the global graph structure. We propose a graph neural network implementation of the ARMA filter with a recursive and distributed formulation, obtaining a convolutional layer that is efficient to train, localized in the node space, and can be transferred to new graphs at test time. We perform a spectral analysis to study the filtering effect of the proposed ARMA layer and report experiments on four downstream tasks: semi-supervised node classification, graph signal classification, graph classification, and graph regression. Results show that the proposed ARMA layer brings significant improvements over graph neural networks based on polynomial filters.


title: ECrits - Visualizing Support Ticket Escalation Risk

abstract: Managing support tickets in large, multi-product organizations is difficult. Failure to meet the expectations of customers can lead to the escalation of support tickets, which is costly for IBM in terms of customer relationships and resources spent addressing the escalation. Keeping the customer happy is an important task in requirements engineering, which often comes in the form of handling their problems brought forth in support tickets. Proper attention to customers, their issues, and the bottom-up requirements that surface through bug reports can be difficult when the support process involves spending a lot of time managing customers to prevent escalations. For any given support analyst, understanding the customer is achievable through time spent looking through past and present support tickets within their organization; however, this solution does not scale up to account for all support tickets across all product teams. ECrits is a tool developed to help mitigate information overload by selectively mining customer information from support ticket repositories, displaying that data to support analysts, and doing predictive modelling on that data to suggest which support tickets are likely to escalate.


title: A squeezing invariant measurement to test displacement of quantum Gaussian states

abstract: We consider a hypothesis testing problem for displacement parameters of n independent copies of an m-mode squeezed quantum Gaussian state whose mixture parameter is known. Given n>1, we construct a quantum measurement as a test using an observable which is invariant by n-fold tensor product of any m-mode squeezing operator. For a pure state case, we calculate the type II error probability of this test. We compare this test with a Hotelling's T-squared test which is based on heterodyne measurements.


title: Efficient Representation Learning Using Random Walks for Dynamic Graphs

abstract: An important part of many machine learning workflows on graphs is vertex representation learning, i.e., learning a low-dimensional vector representation for each vertex in the graph. Recently, several powerful techniques for unsupervised representation learning have been demonstrated to give the state-of-the-art performance in downstream tasks such as vertex classification and edge prediction. These techniques rely on random walks performed on the graph in order to capture its structural properties. These structural properties are then encoded in the vector representation space. However, most contemporary representation learning methods only apply to static graphs while real-world graphs are often dynamic and change over time. Static representation learning methods are not able to update the vector representations when the graph changes; therefore, they must re-generate the vector representations on an updated static snapshot of the graph regardless of the extent of the change in the graph. In this work, we propose computationally efficient algorithms for vertex representation learning that extend random walk based methods to dynamic graphs. The computation complexity of our algorithms depends upon the extent and rate of changes (the number of edges changed per update) and on the density of the graph. We empirically evaluate our algorithms on real world datasets for downstream machine learning tasks of multi-class and multi-label vertex classification. The results show that our algorithms can achieve competitive results to the state-of-the-art methods while being computationally efficient.


title: Learning to Remember More with Less Memorization

abstract: Memory-augmented neural networks consisting of a neural controller and an external memory have shown potentials in long-term sequential learning. Current RAM-like memory models maintain memory accessing every timesteps, thus they do not effectively leverage the short-term memory held in the controller. We hypothesize that this scheme of writing is suboptimal in memory utilization and introduces redundant computation. To validate our hypothesis, we derive a theoretical bound on the amount of information stored in a RAM-like system and formulate an optimization problem that maximizes the bound. The proposed solution dubbed Uniform Writing is proved to be optimal under the assumption of equal timestep contributions. To relax this assumption, we introduce modifications to the original solution, resulting in a solution termed Cached Uniform Writing. This method aims to balance between maximizing memorization and forgetting via overwriting mechanisms. Through an extensive set of experiments, we empirically demonstrate the advantages of our solutions over other recurrent architectures, claiming the state-of-the-arts in various sequential modeling tasks.


title: Coding Techniques for Future Wi-Fi

abstract: Future wireless applications such as high definition video streaming, wireless cloud radio access networks, and cellular data offload are bandwidth-hungry, which highlights the need for Wi-Fi links exceeding 1 Gb/s. This article discusses coding schemes and, in particular, the use of Low-Density Parity-Check (LDPC) codes as Forward Error Correction (FEC) in future Wi-Fi standards. Moreover, we consider advanced design strategies such as root-check LDPC structures, computer-aided design for short blocks and high-performance decoding algorithms with low latency. We then conclude the article with a discussion of FEC challenges for future Wi-FI applications.


title: Tunable Magnetic Skyrmions in Ultrathin Magnetic Nanostructures for Cellular-Level Neurostimulation

abstract: In 2016, the Global Burden of Disease reported that neurological disorders were the principal cause of disability-adjusted life years (DALYs) and the second leading cause of deaths. Research in the last decade has pushed neuroscience to design and implement low-cost, efficient, implantable, flexible electrodes/probes and 3D arrays for neuron stimulation and sensing. Electrical arrays used in current CMOS-based technologies can be affected by the migration of cells (such as astrocytes) that attempt to seal off the electronic devices, causing increased impedance and alternations in the electric field. In this regard, magnetic nanodevices can be better candidates. A wide assortment of magnetic skyrmion-based device ideas and models have as of late been proposed featuring their potential applications. In this paper we propose a highly tunable skyrmion-based spintronic nanodevice for neuron stimulation. The effects of tunable material and magnetic properties specifically Dzyaloshinskii-Moriya interaction (DMI), perpendicular magnetic anisotropy (PMA) constant, number of skyrmions and device dimensions on stable skyrmion nucleation and smooth skyrmion dynamics in a magnetic ultra-thin film have been extensively studied. The aim of this study was to meet the standard therapeutic specifications of neuron stimulation, which is an electric field of about 10 mV/mm for a duration of 50 {\mu}s. From Faraday's Laws of Induction, skyrmion dynamics that generates an alternating magnetic flux density induces an electric field for a certain time duration in the ultra-thin film. The results of this work show that on tuning the skyrmion dynamics, the induced electric field can reach the standard value required for neurostimulation, thereby providing a strong futuristic possibility to exploit skyrmion-based spintronic nanodevices for neuron stimulation.


title: Viscosity independent diffusion mediated by death and reproduction in biofilms

abstract: Bacterial biofilms, surface-attached communities of cells, are in some respects similar to colloidal solids; both are densely packed with non-zero yield stresses. However, unlike non-living materials, bacteria reproduce and die, breaking mechanical equilibrium and inducing collective dynamic responses. We report experiments and theory investigating the motion of immotile Vibrio cholerae, which can kill each other and reproduce in biofilms. We vary viscosity by using bacterial variants that secrete different amounts of extracellular matrix polymers, but are otherwise identical. Unlike thermally-driven diffusion, in which diffusivity decreases with increased viscosity, we find that cellular motion mediated by death and reproduction is independent of viscosity over timescales relevant to bacterial reproduction. To understand this surprising result, we use two separate modeling approaches. First we perform explicitly mechanical simulations of one-dimensional chains of Voigt-Kelvin elements that can die and reproduce. Next, we perform an independent statistical approach, modeling Brownian motion with the classic Langevin equation under an effective temperature that depends on cellular division rate. The diffusion of cells in both approaches agrees quite well, supporting a kinetic interpretation for the effective temperature used here and developed in previous work. As the viscoelastic behavior of biofilms is believed to play a large role in their anomalous biological properties, such as antibiotic resistance, the independence of cellular diffusive motion --- important for biofilm growth and remodeling --- on viscoelastic properties likely holds ecological, medical, and industrial relevance.


title: A surface in odd characteristic with discrete and non-finitely generated automorphism group

abstract: It was proved by Tien-Cuong Dinh and me that there is a smooth complex projective surface whose automorphism group is discrete and not finitely generated. In this paper, we will show that there is a smooth projective surface, birational to some K3 surface, such that the automorphism group is discrete and not finitely generated, over any algebraically closed field of odd characteristic except precisely an algebraic closure of the prime field.


title: A duality formula between elliptic determinants

abstract: We prove a duality formula between two elliptic determinants. We present a proof which is a variant of the Izergin-Korepin method which is a method originally introduced to analyze and compute partition functions of integrable lattice models.


title: Observations of A Fast-Expanding and UV-Bright Type Ia Supernova SN 2013gs

abstract: In this paper, we present extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2013gs discovered during the Tsinghua-NAOC Transient Survey. The photometric observations in the optical show that the light curves of SN 2013gs is similar to that of normal SNe Ia, with an absolute peak magnitude of $M_{B}$ = $-$19.25 $\pm$ 0.15 mag and a post-maximum decline rate $\Delta$m$_{15}$(B) = 1.00 $ \pm $ 0.05 mag. \emph{Gehrels Swift} UVOT observations indicate that SN 2013gs shows unusually strong UV emission (especially in the $uvw1$ band) at around the maximum light (M$_{uvw1}$ $\sim$ $-$18.9 mag). The SN is characterized by relatively weak Fe~{\sc ii} {\sc iii} absorptions at $\sim$ 5000{\AA} in the early spectra and a larger expansion velocity ($v_{Si}$ $\sim$ 13,000 km s$^{-1}$ around the maximum light) than the normal-velocity SNe Ia. We discuss the relation between the $uvw1-v$ color and some observables, including Si~{\sc ii} velocity, line strength of Si~{\sc ii} $\lambda$6355, Fe~{\sc ii}/{\sc iii} lines and $\Delta m_{15}$(B). Compared to other fast-expanding SNe Ia, SN 2013gs exhibits Si and Fe absorption lines with similar strength and bluer $uvw1-v$ color. We briefly discussed the origin of the observed UV dispersion of SNe Ia.


title: Metrics matter in community detection

abstract: We present a critical evaluation of normalized mutual information (NMI) as an evaluation metric for community detection. NMI exaggerates the leximin method's performance on weak communities: Does leximin, in finding the trivial singletons clustering, truly outperform eight other community detection methods? Three NMI improvements from the literature are AMI, rrNMI, and cNMI. We show equivalences under relevant random models, and for evaluating community detection, we advise one-sided AMI under the $\mathbb{M}_{\mathrm{all}}$ model (all partitions of $n$ nodes). This work seeks (1) to start a conversation on robust measurements, and (2) to advocate evaluations which do not give "free lunch".


title: Investigating the Effect of Magnetic Dipole-Dipole Interaction on Magnetic Particle Spectroscopy (MPS): Implications for Magnetic Nanoparticle-based Bioassays and Magnetic Particle Imaging (MPI)

abstract: Superparamagnetic iron oxide nanoparticles (SPIONs), with comparable size to biomolecules (such as proteins, nucleic acids, etc.) and unique magnetic properties, good biocompatibility, low toxicity, potent catalytic behavior, are promising candidates for many biomedical applications. There is one property present in most SPION systems, yet it has not been fully exploited, which is the dipole-dipole interaction (also called dipolar interaction) between the SPIONs. It is known that the magnetic dynamics of an ensemble of SPIONs are substantially influenced by the dipolar interactions. However, the exact way it affects the performance of magnetic particle-based bioassays and magnetic particle imaging (MPI) is still an open question. The purpose of this paper is to give a partial answer to this question. This is accomplished by numerical simulations on the dipolar interactions between two nearby SPIONs and experimental measurements on an ensemble of SPIONs using our lab-based magnetic particle spectroscopy (MPS) system. Our results show that even moderate changes in the SPION concentration may have substantial effects on the magnetic dynamics of the SPION system and the harmonic signal magnitudes can be increased or decreased by 60%, depending on the values of MPS system parameters.


title: Exponential Strong Converse for Successive Refinement with Causal Decoder Side Information

abstract: We consider the $k$-user successive refinement problem with causal decoder side information and derive an exponential strong converse theorem. The rate-distortion region for the problem can be derived as a straightforward extension of the two-user case by Maor and Merhav (2008). We show that for any rate-distortion tuple outside the rate-distortion region of the $k$-user successive refinement problem with causal decoder side information, the joint excess-distortion probability approaches one exponentially fast. Our proof follows by judiciously adapting the recently proposed strong converse technique by Oohama using the information spectrum method, the variational form of the rate-distortion region and H\"older's inequality. The lossy source coding problem with causal decoder side information considered by El Gamal and Weissman is a special case ($k=1$) of the current problem. Therefore, the exponential strong converse theorem for the El Gamal and Weissman problem follows as a corollary of our result.


title: Connected sum of CR manifolds with positive CR Yamabe constant

abstract: Suppose $M_{1}$ and $M_{2}$ are $3$-dimensional closed (compact without boundary) CR manifolds with positive CR Yamabe constant. In this note, we show that the connected sum of $M_{1}$ and $M_{2}$ also admits a CR structure with positive CR Yamabe constant.


title: Stochastic Modeling of 3-D Compositional Distribution in the Crust with Bayesian Inference and Application to Geoneutrino Observation in Japan

abstract: Geoneutrino observations, first achieved by KamLAND in 2005 and followed by Borexino in 2010, have accumulated statistics and improved sensitivity for more than ten years. The uncertainty of the geoneutrino flux at the surface is now reduced to a level small enough to set useful constraints on U and Th abundances in the bulk silicate earth (BSE). However, in order to make inferences on earth's compositional model, the contributions from the local crust need to be understood within a similar uncertainty. Here we develop a new method to construct a stochastic crustal composition model utilizing Bayesian inference. While the methodology has general applicability, it incorporates all the local uniqueness in its probabilistic framework. Unlike common approaches for this type of problem, our method does not depend on crustal segmentation into upper, (middle) and lower, whose classification and boundaries are not always well defined. We also develop a new modeling method to infer rock composition distributions that conserve mass balance and therefore do not bias the results. Combined with a new vast collection of geochemical data for rock samples in the Japan arc, we apply this method to geoneutrino observation at Kamioka, Japan. Currently a difficulty remains in the handling of correlations in the flux integration; we conservatively assume maximum correlation, which leads to large flux estimation errors of 60~70%. Despite the large errors, this is the first local crustal model for geoneutrino flux prediction with probabilistic error estimation in a reproducible way.


title: Data-Driven Reconstruction of Gravitationally Lensed Galaxies using Recurrent Inference Machines

abstract: We present a machine learning method for the reconstruction of the undistorted images of background sources in strongly lensed systems. This method treats the source as a pixelated image and utilizes the Recurrent Inference Machine (RIM) to iteratively reconstruct the background source given a lens model. Our architecture learns to minimize the likelihood of the model parameters (source pixels) given the data using the physical forward model (ray tracing simulations) while implicitly learning the prior of the source structure from the training data. This results in better performance compared to linear inversion methods, where the prior information is limited to the 2-point covariance of the source pixels approximated with a Gaussian form, and often specified in a relatively arbitrary manner. We combine our source reconstruction network with a convolutional neural network that predicts the parameters of the mass distribution in the lensing galaxies directly from telescope images, allowing a fully automated reconstruction of the background source images and the foreground mass distribution.


title: On weaving frames in Hilbert spaces

abstract: In this paper, we obtain some new properties of weaving frames and present some conditions under which a family of frames is woven in Hilbert spaces. Some characterizations of weaving frames in terms of operators are given. We also give a condition associated with synthesis operators of frames such that the sequence of frames is woven. Finally, for a family of woven frames, we show that they are stable under invertible operators and small perturbations.


title: Estimation of Planetary Photometric Emissions for Extremely Close-in Exoplanets

abstract: The fine precision of photometric data available from missions like Kepler provide researchers with the ability to measure changes in light on the order of tens of parts per million (ppm). This level of precision allows researchers to measure the loss of light due to exoplanet transits as well as the light emitted by an exoplanet, or planetary photometric emissions. The planetary photometric emissions are due to the thermal emissions of the exoplanet, and reflected stellar light. In many cases it is assumed that the incident stellar light may be modeled as plane parallel rays. For extremely close-in exoplanets the finite angular size of the host star must be taken into account and the plane parallel ray model breaks down. One consequence of modeling the incident stellar radiation in this manner is the creation of three distinct zones as opposed to the two zones present in the plane parallel ray model. The three zones are the fully illuminated, penumbral, and un-illuminated zones. The existence of the penumbral zone means that more than half of the exoplanet will be at least partially illuminated by the host star. In this work we will present a complete derivation of the reflective luminosity of the fully illuminated zone. In addition, we will present an outline for the derivation of the reflected intensity distribution of the penumbral zone. Within this work we will also derive a new expression for the thermal luminosity of an exoplanet by treating each of the three zones as a blackbody emitting at a constant temperature. Finally, it will be shown that an estimation of the radius of the exoplanet requires proper accounting of light from the penumbral zone during the primary transit. Not doing so risks underestimating the radius of the exoplanet and an overestimation of its geometric albedo and nightside temperature.


title: Hacking the Brain: Triggering Neuroplasticity for Enhancing Musical Talent: A study on Monkey and Human behavior after Exposure to Videogames and Visual/Auditory Stimuli to Increase Musical Abilities through Neuroplasticity

abstract: In this paper, we analyze the cognitive improvements that can be achieved through hacking the brain through the use of multiple methods to enhance neuroplasticity. Exposure to gaming, for example, has proven conducive for learning real-world abilities through auditory and visual stimuli. We will discuss cortical magnification and receptive field sizes, as well as topographic brain maps in the context of neuroplasticity. We finally propose more studies to be performed to improve musical talent and musical abilities.


title: $C^3$ matching for asymptotically flat spacetimes

abstract: We propose a criterion for finding the minimum distance at which an interior solution of Einstein's equations can be matched with an exterior asymptotically flat solution. It is based upon the analysis of the eigenvalues of the Riemann curvature tensor and their first derivatives, implying $C^3$ differentiability conditions. The matching itself is performed by demanding continuity of the curvature eigenvalues across the matching surface. We apply the $C^3$ matching approach to spherically symmetric perfect fluid spacetimes and obtain the physically meaningful condition that density and pressure should vanish on the matching surface. Several perfect fluid solutions in Newton and Einstein gravity are tested.


title: A Divergence-Conforming Hybridized Discontinuous Galerkin Method for the Incompressible Reynolds Averaged Navier-Stokes Equations

abstract: We introduce a hybridized discontinuous Galerkin method for the incompressible Reynolds Averaged Navier-Stokes equations coupled with the Spalart-Allmaras one equation turbulence model. With a special choice of velocity and pressure spaces for both element and trace degrees of freedom, we arrive at a method which returns point-wise divergence-free mean velocity fields and properly balances momentum and energy. We further examine the use of different polynomial degrees and meshes to see how the order of the scalar eddy viscosity affects the convergence of the mean velocity and pressure fields, specifically for the method of manufactured solutions. As is standard with hybridized discontinuous Galerkin methods, static condensation can be employed to remove the element degrees of freedom and thus dramatically reduce the global number of degrees of freedom. Numerical results illustrate the effectiveness of the proposed methodology.


title: Hierarchical Reinforcement Learning via Advantage-Weighted Information Maximization

abstract: Real-world tasks are often highly structured. Hierarchical reinforcement learning (HRL) has attracted research interest as an approach for leveraging the hierarchical structure of a given task in reinforcement learning (RL). However, identifying the hierarchical policy structure that enhances the performance of RL is not a trivial task. In this paper, we propose an HRL method that learns a latent variable of a hierarchical policy using mutual information maximization. Our approach can be interpreted as a way to learn a discrete and latent representation of the state-action space. To learn option policies that correspond to modes of the advantage function, we introduce advantage-weighted importance sampling. In our HRL method, the gating policy learns to select option policies based on an option-value function, and these option policies are optimized based on the deterministic policy gradient method. This framework is derived by leveraging the analogy between a monolithic policy in standard RL and a hierarchical policy in HRL by using a deterministic option policy. Experimental results indicate that our HRL approach can learn a diversity of options and that it can enhance the performance of RL in continuous control tasks.


title: Nonlinear excitation of geodesic acoustic mode by toroidal Alfv\'en eigenmodes and impact on plasma performance

abstract: Spontaneous nonlinear excitation of geodesic acoustic mode (GAM) by toroidal Alfv\'en eigenmode (TAE) is investigated using nonlinear gyrokinetic theory. It is found that, the nonlinear decay process depends on thermal ion \beta value. Here, \beta is the plasma thermal to magnetic pressure ratio. In the low-\beta limit, TAE decays into a GAM and a lower TAE sideband in the toroidicity induced shear Alfv\'en wave continuous spectrum gap; while in the high-\beta limit, TAE decays into a GAM and a propagating kinetic TAE in the continuum. Both cases are investigated for the spontaneous decay conditions. The nonlinear saturation levels of both GAM and daughter wave are derived. The corresponding power balance and wave particle power transfer to thermal plasma are computed. Implications on thermal plasma heating are also discussed.


title: Instability of unidirectional flows for the 2D $\alpha$-Euler equations

abstract: We study stability of unidirectional flows for the linearized 2D $\alpha$-Euler equations on the torus. The unidirectional flows are steady states whose vorticity is given by Fourier modes corresponding to a vector $\mathbf p \in \mathbb Z^{2}$. We linearize the $\alpha$-Euler equation and write the linearized operator $L_{B} $ in $\ell^{2}(\mathbb Z^{2})$ as a direct sum of one-dimensional difference operators $L_{B,\mathbf q}$ in $\ell^{2}(\mathbb Z)$ parametrized by some vectors $\mathbf q\in\mathbb Z^2$ such that the set $\{\mathbf q +n \mathbf p:n \in \mathbb Z\}$ covers the entire grid $\mathbb Z^{2}$. The set $\{\mathbf q +n \mathbf p:n \in \mathbb Z\}$ can have zero, one, or two points inside the disk of radius $\|\mathbf p\|$. We consider the case where the set $\{\mathbf q +n \mathbf p:n \in \mathbb Z\}$ has exactly one point in the open disc of radius $\mathbf p$. We show that unidirectional flows that satisfy this condition are linearly unstable. Our main result is an instability theorem that provides a necessary and sufficient condition for the existence of a positive eigenvalue to the operator $L_{B,\mathbf q}$ in terms of equations involving certain continued fractions. Moreover, we are also able to provide a complete characterization of the corresponding eigenvector. The proof is based on the use of continued fractions techniques expanding upon the ideas of Friedlander and Howard.


title: Demonstration of a strain-mediated magnetoelectric write and read unit in a Co60Fe20B20/ Pb(Mg1/3Nb2/3)0.7Ti0.3O3 heterostructure

abstract: Taking advantage of the Magnetoelectric (ME) and its inverse effect, this article demonstrates strain-mediated magnetoelectric write and read operations simultaneously in Co60Fe20B20/ Pb(Mg1/3Nb2/3)0.7Ti0.3O3 heterostructures without using any symmetry breaking magnetic field at room temperature. By applying an external DC-voltage across a (011)-cut PMN-PT substrate, the ferroelectric polarization is re-oriented, which results in an anisotropic in-plane strain that transfers to the CoFeB thin film and changes its magnetic anisotropy Hk. The change in Hk in-turn results in a 90o rotation of the magnetic easy axis for sufficiently high voltages. Simultaneously, the inverse effect is employed to read changes of the magnetic properties. Because the Piezoelectric (PE)/FerroMagnetic (FM) system is fully coupled, the change of magnetization in FM induces an elastic stress in the PE layer, which generates a piezoelectric potential in the system that can be used to readout the magnetic state of the FM layer. Our experimental results are in excellent qualitative agreement with a recently proposed, experimentally benchmarked equivalent circuit model that considers how magnetic properties are electrically controlled in such ME/PE heterostructure and how a back-voltage is generated due to changing magnetic properties in a self-consistent model.


title: Adaptive Fusion for RGB-D Salient Object Detection

abstract: RGB-D salient object detection aims to identify the most visually distinctive objects in a pair of color and depth images. Based upon an observation that most of the salient objects may stand out at least in one modality, this paper proposes an adaptive fusion scheme to fuse saliency predictions generated from two modalities. Specifically, we design a two-streamed convolutional neural network (CNN), each of which extracts features and predicts a saliency map from either RGB or depth modality. Then, a saliency fusion module learns a switch map that is used to adaptively fuse the predicted saliency maps. A loss function composed of saliency supervision, switch map supervision, and edge-preserving constraints is designed to make full supervision, and the entire network is trained in an end-to-end manner. Benefited from the adaptive fusion strategy and the edge-preserving constraint, our approach outperforms state-of-the-art methods on three publicly available datasets.


title: Stereoscopic Dark Flash for Low-light Photography

abstract: In this work, we present a camera configuration for acquiring "stereoscopic dark flash" images: a simultaneous stereo pair in which one camera is a conventional RGB sensor, but the other camera is sensitive to near-infrared and near-ultraviolet instead of R and B. When paired with a "dark" flash (i.e., one having near-infrared and near-ultraviolet light, but no visible light) this camera allows us to capture the two images in a flash/no-flash image pair at the same time, all while not disturbing any human subjects or onlookers with a dazzling visible flash. We present a hardware prototype of this camera that approximates an idealized camera, and we present an imaging procedure that let us acquire dark flash stereo pairs that closely resemble those we would get from that idealized camera. We then present a technique for fusing these stereo pairs, first by performing registration and warping, and then by using recent advances in hyperspectral image fusion and deep learning to produce a final image. Because our camera configuration and our data acquisition process allow us to capture true low-noise long exposure RGB images alongside our dark flash stereo pairs, our learned model can be trained end-to-end to produce a fused image that retains the color and tone of a real RGB image while having the low-noise properties of a flash image.


title: A Non-Linear Roth Theorem for Sets of Positive Density

abstract: Suppose that $A \subset \mathbb{R}$ has positive upper density, \[ \limsup_{|I| \to \infty} \frac{|A \cap I|}{|I|} = \delta > 0,\] and $P(t) \in \mathbb{R}[t]$ is a polynomial with no constant or linear term, or more generally a non-flat curve: a locally differentiable curve which doesn't "resemble a line" near $0$ or $\infty$. Then for any $R_0 \leq R$ sufficiently large, there exists some $x_R \in A$ so that \[ \inf_{R_0 \leq T \leq R} \frac{|\{ 0 \leq t < T : x_R - t \in A, \ x_R - P(t) \in A \}|}{T} \geq c_P \cdot \delta^2 \] for some absolute constant $c_P > 0$, that depends only on $P$.


title: Monochromatic disconnection of graphs

abstract: For an edge-colored graph $G$, we call an edge-cut $M$ of $G$ monochromatic if the edges of $M$ are colored with a same color. The graph $G$ is called monochromatically disconnected if any two distinct vertices of $G$ are separated by a monochromatic edge-cut. For a connected graph $G$, the monochromatic disconnection number, denoted by $md(G)$, of $G$ is the maximum number of colors that are needed in order to make $G$ monochromatically disconnected. We will show that almost all graphs have monochromatic disconnection numbers equal to 1. We also obtain the Nordhaus-Gaddum-type results for $md(G)$.


title: An arbitrary two-particle high-dimensional Bell state measurement by auxiliary entanglement

abstract: Bell state measurement (BSM) plays a vital role in quantum information. There are many researches on BSM of qubit Bell state, however, there is no definite solution of how to realize high-dimensional Bell state measurement (HDBSM). In this paper, We present a scheme for realizing arbitrary high-dimensional two-particle Bell state measurement by auxiliary entanglement. In our scheme, the maximal entangled states with different degrees of freedom (DOF) are used as auxiliary states, and then the HDBSM can be achieved by implementing Bell measurement with different degrees of freedom. We present the detailed HDBSM scheme for three and four-dimensional Bell state, and the general formula for arbitrary high-dimensional Bell state. Furthermore, We have designed the experimental scheme for three-dimensional BSM, and given the method of arbitrary dimension. This method can promote the high-dimensional quantum information task.


title: An Irregular-shaped Ring-Pair Magnet Array with a Monotonic Field Gradient for 2D Head Imaging in Low-field Portable MRI

abstract: We present a design and the optimization of an irregular-shaped ring-pair magnet array that generates a 1D monotonic field pattern for 2D head imaging in a low-field portable MRI system. The magnet rings are discretized into fan-shaped ring segments with varying outer diameters for the design and optimization. Besides, the inner radii of ring-pairs are tapered from outside in to provide the controlled field inhomogeneity. Genetic algorithm (GA) was used, and a current model for a fan-shaped ring segment was derived to have a fast forward calculation in the optimization. A monotonic field pattern is successfully obtained along the x-direction in a cylindrical field of view (FoV), with a relatively strong magnetic field (132.98 mT) and the homogeneity of 151840 ppm. The proposed array was further evaluated by applying its field as a spatial encoding magnetic field (SEM) for imaging numerically. Due to the field monotonicity, the reconstructed image by applying the fields of the proposed array shows clearer features (a higher structural similarity index) with a reduced error rate compared to that using a sparse dipolar Halbach array. The proposed magnet array is a promising alternative to supply SEM for imaging in a permanent-magnet-based low-field portable MRI system.


title: Analysis of a Two-Layer Neural Network via Displacement Convexity

abstract: Fitting a function by using linear combinations of a large number $N$ of `simple' components is one of the most fruitful ideas in statistical learning. This idea lies at the core of a variety of methods, from two-layer neural networks to kernel regression, to boosting. In general, the resulting risk minimization problem is non-convex and is solved by gradient descent or its variants. Unfortunately, little is known about global convergence properties of these approaches. Here we consider the problem of learning a concave function $f$ on a compact convex domain $\Omega\subseteq {\mathbb R}^d$, using linear combinations of `bump-like' components (neurons). The parameters to be fitted are the centers of $N$ bumps, and the resulting empirical risk minimization problem is highly non-convex. We prove that, in the limit in which the number of neurons diverges, the evolution of gradient descent converges to a Wasserstein gradient flow in the space of probability distributions over $\Omega$. Further, when the bump width $\delta$ tends to $0$, this gradient flow has a limit which is a viscous porous medium equation. Remarkably, the cost function optimized by this gradient flow exhibits a special property known as displacement convexity, which implies exponential convergence rates for $N\to\infty$, $\delta\to 0$. Surprisingly, this asymptotic theory appears to capture well the behavior for moderate values of $\delta, N$. Explaining this phenomenon, and understanding the dependence on $\delta,N$ in a quantitative manner remains an outstanding challenge.


title: An Empirical Study of Speculative Concurrency in Ethereum Smart Contracts

abstract: We use historical data to estimate the potential benefit of speculative techniques for executing Ethereum smart contracts in parallel. We replay transaction traces of sampled blocks from the Ethereum blockchain over time, using a simple speculative execution engine. In this engine, miners attempt to execute all transactions in a block in parallel, rolling back those that cause data conflicts. Aborted transactions are then executed sequentially. Validators execute the same schedule as miners. We find that our speculative technique yields estimated speed-ups starting at about 8-fold in 2016, declining to about 2-fold at the end of 2017, where speed-up is measured using either gas costs or instruction counts. We also observe that a small set of contracts are responsible for many data conflicts resulting from speculative concurrent execution.


title: Population-Guided Large Margin Classifier for High-Dimension Low -Sample-Size Problems

abstract: Various applications in different fields, such as gene expression analysis or computer vision, suffer from data sets with high-dimensional low-sample-size (HDLSS), which has posed significant challenges for standard statistical and modern machine learning methods. In this paper, we propose a novel linear binary classifier, denoted by population-guided large margin classifier (PGLMC), which is applicable to any sorts of data, including HDLSS. PGLMC is conceived with a projecting direction w given by the comprehensive consideration of local structural information of the hyperplane and the statistics of the training samples. Our proposed model has several advantages compared to those widely used approaches. First, it is not sensitive to the intercept term b. Second, it operates well with imbalanced data. Third, it is relatively simple to be implemented based on Quadratic Programming. Fourth, it is robust to the model specification for various real applications. The theoretical properties of PGLMC are proven. We conduct a series of evaluations on two simulated and six real-world benchmark data sets, including DNA classification, digit recognition, medical image analysis, and face recognition. PGLMC outperforms the state-of-the-art classification methods in most cases, or at least obtains comparable results.


title: Matrix versions of the Hellinger distance

abstract: On the space of positive definite matrices we consider distance functions of the form $d(A,B)=\left[\tr\mathcal{A}(A,B)-\tr\mathcal{G}(A,B)\right]^{1/2},$ where $\mathcal{A}(A,B)$ is the arithmetic mean and $\mathcal{G}(A,B)$ is one of the different versions of the geometric mean. When $\mathcal{G}(A,B)=A^{1/2}B^{1/2}$ this distance is $\|A^{1/2}-B^{1/2}\|_2,$ and when $\mathcal{G}(A,B)=(A^{1/2}BA^{1/2})^{1/2}$ it is the Bures-Wasserstein metric. We study two other cases: $\mathcal{G}(A,B)=A^{1/2}(A^{-1/2}BA^{-1/2})^{1/2}A^{1/2},$ the Pusz-Woronowicz geometric mean, and $\mathcal{G}(A,B)=\exp\big(\frac{\log A+\log B}{2}\big),$ the log Euclidean mean. With these choices $d(A,B)$ is no longer a metric, but it turns out that $d^2(A,B)$ is a divergence. We establish some (strict) convexity properties of these divergences. We obtain characterisations of barycentres of $m$ positive definite matrices with respect to these distance measures.