title: Some problems on ruled hypersurfaces in nonflat complex space forms

abstract: We study ruled real hypersurfaces whose shape operators have constant squared norm in nonflat complex space forms. In particular, we prove the nonexistence of such hypersurfaces in the projective case. We also show that biharmonic ruled real hypersurfaces in nonflat complex space forms are minimal, which provides their classification due to a known result of Lohnherr and Reckziegel.


title: Improving the Drell--Yan probe of small x partons at the LHC via an azimuthal angle cut

abstract: Predictions for Drell--Yan lepton pair production at low dilepton mass and small x at the LHC usually have a large scale dependence. This can be decreased by determining an "optimal factorization scale". In this paper, we reduce this scale by imposing a cutoff in azimuthal angle between the transverse momentum of the leptons, properly taking into account Sudakov effects. This allows one to probe the parton distributions at smaller scales eliminating most of the current theoretical uncertainty.


title: Magnetic fluctuations and superconducting pairing in $\varepsilon$-iron

abstract: We study Coulomb correlation effects and their role in superconductivity of $\varepsilon$-iron under pressure from 12 to 33 GPa by using a combination of density functional and dynamical mean-field theory. Our results indicate a persistence of the Fermi-liquid behavior below the temperature $\sim$1000 K. The Coulomb correlations are found to substantially renormalize the density of states, reducing the distance from the peak to the Fermi level to 0.4 eV compared to 0.75 eV obtained in DFT calculations. We find significant antiferromagnetic correlations, which are accompanied by the formation of short-lived local magnetic moments. We use the obtained results as a starting point for construction of the multi-band Bethe-Salpeter equation, which eigenvalues indicate that antiferromagnetic spin fluctuations may result in the superconducting pairing in $\varepsilon$-Fe. Moreover, the tendency to superconducting instability becomes weaker with the increase of pressure, which may explain the disappearance of superconductivity at $\sim$30 GPa.


title: FastGAE: Scalable Graph Autoencoders with Stochastic Subgraph Decoding

abstract: Graph autoencoders (AE) and variational autoencoders (VAE) are powerful node embedding methods, but suffer from scalability issues. In this paper, we introduce FastGAE, a general framework to scale graph AE and VAE to large graphs with millions of nodes and edges. Our strategy, based on an effective stochastic subgraph decoding scheme, significantly speeds up the training of graph AE and VAE while preserving or even improving performances. We demonstrate the effectiveness of FastGAE on various real-world graphs, outperforming the few existing approaches to scale graph AE and VAE by a wide margin.


title: From Chiral Kinetic Theory To Spin Hydrodynamics

abstract: The total angular momentum is conserved in the evolution of the Quark-Gluon Plasma (QGP) created in heavy-ion collision, and consists of two sectors: the orbital angular momentum (OAM) caused by kinetic motion, and the spin, an intrinsic degree of freedom of quarks and gluons. Microscopic scattering processes allow the conversion between these two components, hence the spin density could eventually have non-trivial influence on the QGP evolution. A hydrodynamic theory, with the spin polarization effect properly taken into account, is required to quantitatively study the polarization rate for observed hadrons, e.g. the $\Lambda$-hyperon. In this work, we start with chiral kinetic theory and construct the spin hydrodynamic framework for a chiral spinor system. We obtain the equations of motion of second-order dissipative relativistic fluid dynamics with non-trivial spin polarization density.


title: Phonon-induced dephasing in quantum dot-cavity QED: Limitations of the polaron master equation

abstract: A semiconductor quantum dot (QD) embedded within an optical microcavity is a system of fundamental importance within quantum information processing. The optimization of quantum coherence is crucial in such applications, requiring an in-depth understanding of the relevant decoherence mechanisms. We provide herein a critical review of prevalent theoretical treatments of the QD-cavity system coupled to longitudinal acoustic phonons, comparing predictions against a recently obtained exact solution. Within this review we consider a range of temperatures and exciton-cavity coupling strengths. Predictions of the polaron Nakajima-Zwanzig (NZ) and time-convolutionless (TCL) master equations, as well as a variation of the former adapted for adiabatic continuous wave excitation (CWE), are compared against an asymptotically exact solution based upon Trotter's decomposition (TD) theorem. The NZ and TCL implementations, which apply a polaron transformation to the Hamiltonian and subsequently treat the exciton-cavity coupling to second order, do not offer a significant improvement accuracy relative to the polaron transformation alone. The CWE adaptation provides a marked improvement, capturing the broadband features of the absorption spectrum (not present in NZ and TCL implementations). We attribute this difference to the effect of the Markov approximation, and particularly its unsuitability in pulsed excitation regime. Even the CWE adaptation, however, breaks down in the regime of high temperature ($50K$) and strong exciton-cavity coupling ($g \gtrsim 0.2$ meV). The TD solution is of comparable computational complexity to the above-mentioned master equation approaches, yet remains accurate at higher temperatures and across a broad range of exciton-cavity coupling strengths (at least up to $g=1.5$ meV).


title: Vehicle Ego-Lane Estimation with Sensor Failure Modeling

abstract: We present a probabilistic ego-lane estimation algorithm for highway-like scenarios that is designed to increase the accuracy of the ego-lane estimate, which can be obtained relying only on a noisy line detector and tracker. The contribution relies on a Hidden Markov Model (HMM) with a transient failure model. The proposed algorithm exploits the OpenStreetMap (or other cartographic services) road property lane number as the expected number of lanes and leverages consecutive, possibly incomplete, observations. The algorithm effectiveness is proven by employing different line detectors and showing we could achieve much more usable, i.e. stable and reliable, ego-lane estimates over more than 100 Km of highway scenarios, recorded both in Italy and Spain. Moreover, as we could not find a suitable dataset for a quantitative comparison with other approaches, we collected datasets and manually annotated the Ground Truth about the vehicle ego-lane. Such datasets are made publicly available for usage from the scientific community.


title: Topological spin torque emerging in classical-spin systems with different time scales

abstract: In classical spin systems with two largely different inherent time scales, the configuration of the fast spins almost instantaneously follows the slow-spin dynamics. We develop the emergent effective theory for the slow-spin degrees of freedom and demonstrate that this generally includes a topological spin torque. This torque gives rise to anomalous real-time dynamics. It derives from the holonomic constraints defining the fast-spin configuration space and is given in terms of a topological charge density which becomes a quantized homotopy invariant when integrated.

categories: cond-mat.mes-hall cond-mat.str-el physics.atom-ph physics.chem-ph


title: Non-abelian tri-vector deformations in d=11 supergravity

abstract: A truncation of the SL(5) Exceptional Field Theory that allows to describe spacetimes of the form $M_4 \times M_7$ with the 4-form flux on $M_4$ is constructed. The resulting theory is used to test the recently proposed tri-vector generalisation of Yang-Baxter deformations applied to the AdS${}_4 \times \mathbb{S}^7$ solution in $d=11$ supergravity. We present two new supergravity solutions corresponding to non-abelian non-unimodular tri-vector deformations of AdS${}_4 \times \mathbb{S}^7$.


title: Status Updates with Priorities: Lexicographic Optimality

abstract: In this paper, we consider a transmission scheduling problem, in which several streams of status update packets with diverse priority levels are sent through a shared channel to their destinations. We introduce a notion of Lexicographic age optimality, or simply lex-age-optimality, to evaluate the performance of multi-class status update policies. In particular, a lex-age-optimal scheduling policy first minimizes the Age of Information (AoI) metrics for high-priority streams, and then, within the set of optimal policies for high-priority streams, achieves the minimum AoI metrics for low-priority streams. We propose a new scheduling policy named Preemptive Priority, Maximum Age First, Last-Generated, First-Served (PP-MAF-LGFS), and prove that the PP-MAF-LGFS scheduling policy is lex-age-optimal. This result holds (i) for minimizing any time-dependent, symmetric, and non-decreasing age penalty function; (ii) for minimizing any non-decreasing functional of the stochastic process formed by the age penalty function; and (iii) for the cases where different priority classes have distinct arrival traffic patterns, age penalty functions, and age penalty functionals. For example, the PP-MAF-LGFS scheduling policy is lex-age-optimal for minimizing the mean peak age of a high-priority stream and the time-average age of a low-priority stream. Numerical results are provided to illustrate our theoretical findings.


title: $L^p$ regularity for a class of averaging operators on the Heisenberg group

abstract: We prove $L^p_{comp}\to L^p_{s}$ boundedness for averaging operators associated to a class of curves in the Heisenberg group $\mathbb{H}^1$ via $L^2$ estimates for related oscillatory integrals and Bourgain-Demeter decoupling inequalities on the cone. We also construct a Sobolev space adapted to translations on the Heisenberg group to which these averaging operators map all $L^p$ functions boundedly.


title: Design and Experiments with a Low-Cost Single-Motor Modular Aquatic Robot

abstract: We present a novel design for a low-cost robotic boat powered by a single actuator, useful for both modular and swarming applications. The boat uses the conservation of angular momentum and passive flippers to convert the motion of a single motor into an adjustable paddling motion for propulsion and steering. We develop design criteria for modularity and swarming and present a prototype implementing these criteria. We identify significant mechanical sensitivities with the presented design, theorize about the cause of the sensitivities, and present an improved design for future work.


title: Pure Differentially Private Summation from Anonymous Messages

abstract: The shuffled (aka anonymous) model has recently generated significant interest as a candidate distributed privacy framework with trust assumptions better than the central model but with achievable errors smaller than the local model. We study pure differentially private (DP) protocols in the shuffled model for summation, a basic and widely used primitive: - For binary summation where each of n users holds a bit as an input, we give a pure $\epsilon$-DP protocol for estimating the number of ones held by the users up to an error of $O_\epsilon(1)$, and each user sends $O_\epsilon(\log n)$ messages each of 1 bit. This is the first pure protocol in the shuffled model with error $o(\sqrt{n})$ for constant $\epsilon$. Using this protocol, we give a pure $\epsilon$-DP protocol that performs summation of real numbers in $[0, 1]$ up to an error of $O_{\epsilon}(1)$, and where each user sends $O_{\epsilon}(\log^3 n)$ messages each of $O(\log\log n)$ bits. - In contrast, we show that for any pure $\epsilon$-DP protocol for binary summation in the shuffled model having absolute error $n^{0.5-\Omega(1)}$, the per user communication has to be at least $\Omega_{\epsilon}(\sqrt{\log n})$ bits. This implies the first separation between the (bounded-communication) multi-message shuffled model and the central model, and the first separation between pure and approximate DP protocols in the shuffled model. To prove our lower bound, we consider (a generalization of) the following question: given $\gamma$ in $(0, 1)$, what is the smallest m for which there are two random variables $X^0, X^1$ supported on $\{0, \dots ,m\}$ such that (i) the total variation distance between $X^0$ and $X^1$ is at least $1-\gamma$, and (ii) the moment generating functions of $X^0$ and $X^1$ are within a constant factor of each other everywhere? We show that the answer is $m = \Theta(\sqrt{\log(1/\gamma)})$.


title: FRB-periodicity: mild pulsars in tight O/B-star binaries

abstract: Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $\eta = L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $\eta \leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.


title: Autonomous Navigation in Unknown Environments using Sparse Kernel-based Occupancy Mapping

abstract: This paper focuses on real-time occupancy mapping and collision checking onboard an autonomous robot navigating in an unknown environment. We propose a new map representation, in which occupied and free space are separated by the decision boundary of a kernel perceptron classifier. We develop an online training algorithm that maintains a very sparse set of support vectors to represent obstacle boundaries in configuration space. We also derive conditions that allow complete (without sampling) collision-checking for piecewise-linear and piecewise-polynomial robot trajectories. We demonstrate the effectiveness of our mapping and collision checking algorithms for autonomous navigation of an Ackermann-drive robot in unknown environments.


title: The space of almost calibrated $(1,1)$ forms on a compact K\"ahler manifold

abstract: The space $\mathcal{H}$ of "almost calibrated" $(1,1)$ forms on a compact K\"ahler manifold plays an important role in the study of the deformed Hermitian-Yang-Mills equation of mirror symmetry as emphasized by recent work of the second author and Yau, and is related by mirror symmetry to the space of positive Lagrangians studied by Solomon. This paper initiates the study of the geometry of $\mathcal{H}$. We show that $\mathcal{H}$ is an infinite dimensional Riemannian manifold with non-positive sectional curvature. In the hypercritical phase case we show that $\mathcal{H}$ has a well-defined metric structure, and that its completion is a ${\rm CAT}(0)$ geodesic metric space, and hence has an intrinsically defined ideal boundary. Finally, we show that in the hypercritical phase case $\mathcal{H}$ admits $C^{1,1}$ geodesics, improving a result of the second author and Yau. Using results of Darvas-Lempert we show that this result is sharp.


title: Junctions of mass-deformed nonlinear sigma models on $SO(2N)/U(N)$ and $Sp(N)/U(N)$ II

abstract: We study vacua, walls and three-pronged junctions of mass-deformed nonlinear sigma models on $SO(2N)/U(N)$ and $Sp(N)/U(N)$ for generic $N$. We review and discuss the on-shell component Lagrangians of the ${\mathcal{N}}=2$ nonlinear sigma model on the Grassmann manifold, which are obtained in the ${\mathcal{N}}=1$ superspace formalism and in the harmonic superspace formalism. We also show that the K\"{a}hler potential of the ${\mathcal{N}}=2$ nonlinear sigma model on the complex projective space, which is obtained in the projective superspace formalism, is equivalent to the K\"{a}hler potential of the ${\mathcal{N}}=2$ nonlinear sigma model with the Fayet-Iliopoulos parameters $c^a=(0,0,c=1)$ on the complex projective space, which is obtained in the ${\mathcal{N}}=1$ superspace formalism.


title: Explicit Wiretap Channel Codes via Source Coding, Universal Hashing, and Distribution Approximation, When the Channels' Statistics are Uncertain

abstract: We consider wiretap channels with uncertainty on the eavesdropper channel under (i) noisy blockwise type II, (ii) compound, or (iii) arbitrarily varying models. We present explicit wiretap codes that can handle these models in a unified manner and only rely on three primitives, namely source coding with side information, universal hashing, and distribution approximation. Our explicit wiretap codes achieve the best known single-letter achievable rates, previously obtained non-constructively, for the models considered. Our results are obtained for strong secrecy, do not require a pre-shared secret between the legitimate users, and do not require any symmetry properties on the channel. An extension of our results to compound main channels is also derived via new capacity-achieving polar coding schemes for compound settings.


title: Machine Learning for Predicting Epileptic Seizures Using EEG Signals: A Review

abstract: With the advancement in artificial intelligence (AI) and machine learning (ML) techniques, researchers are striving towards employing these techniques for advancing clinical practice. One of the key objectives in healthcare is the early detection and prediction of disease to timely provide preventive interventions. This is especially the case for epilepsy, which is characterized by recurrent and unpredictable seizures. Patients can be relieved from the adverse consequences of epileptic seizures if it could somehow be predicted in advance. Despite decades of research, seizure prediction remains an unsolved problem. This is likely to remain at least partly because of the inadequate amount of data to resolve the problem. There have been exciting new developments in ML-based algorithms that have the potential to deliver a paradigm shift in the early and accurate prediction of epileptic seizures. Here we provide a comprehensive review of state-of-the-art ML techniques in early prediction of seizures using EEG signals. We will identify the gaps, challenges, and pitfalls in the current research and recommend future directions.


title: Curvature of magnetic field lines in compressible magnetized turbulence: Statistics, magnetization predictions, gradient curvature, modes and self-gravitating media

abstract: Magnetic field lines in interstellar media have a rich morphology, which could be characterized by geometrical parameters such as curvature and torsion. In this paper, we explore the statistical properties of magnetic field line curvature $\kappa$ in compressible magnetized turbulence. We see that both the mean and standard deviation of magnetic field line curvature obey power-law relations to the magnetization. Moreover, the power-law tail of the curvature probability distribution function is also proportional to the Alfvenic Mach number. We also explore whether the curvature method could be used in the field-tracing Velocity Gradient Technique. In particular, we observe that there is a relation between the mean and standard deviation of the curvature probed by velocity gradients to $M_A$. Finally we discuss how curvature is contributed by different MHD modes in interstellar turbulence, and suggests that the eigenvectors of MHD modes could be possibly represented by the natural Fernet-Serrat frame of the magnetic field lines. We discuss possible theoretical and observational applications of the curvature technique, including the extended understanding on a special length scale that characterize the importance of magnetic field curvature in driving MHD turbulence, and how it could be potentially used to study self-gravitating system.


title: Self-Directed Online Machine Learning for Topology Optimization

abstract: Topology optimization by optimally distributing materials in a given domain requires non-gradient optimizers to solve highly complicated problems. However, with hundreds of design variables or more involved, solving such problems would require millions of Finite Element Method (FEM) calculations whose computational cost is huge and impractical. Here we report Self-directed Online Learning Optimization (SOLO) which integrates Deep Neural Network (DNN) with FEM calculations. A DNN learns and substitutes the objective as a function of design variables. A small number of training data is generated dynamically based on the DNN's prediction of the optimum. The DNN adapts to the new training data and gives better prediction in the region of interest until convergence. The optimum predicted by the DNN is proved to converge to the true global optimum through iterations. Our algorithm was tested by four types of problems including compliance minimization, fluid-structure optimization, heat transfer enhancement and truss optimization. It reduced the computational time by 2 ~ 5 orders of magnitude compared with directly using heuristic methods, and outperformed all state-of-the-art algorithms tested in our experiments. This approach enables solving large multi-dimensional optimization problems.


title: A metric space with transfinite asymptotic dimension $2\omega$

abstract: We build an example of a metric space with transfinite asymptotic dimension $2\omega$.


title: Fermion masses and mixing in modular A4 Symmetry

abstract: We consider a flavor model based on A4 modular group to account for both lepton and quark parameters (masses and mixing). The inverse seesaw mechanism is considered to produce the light neutrino masses. Lepton masses and mixing are obtained in terms of Yukawa coupling ratios and values of the modulus {\tau} nearby some fixed points for inverted neutrino mass hierarchy. The quark masses and mixing are arisen at the same {\tau} values used in inverted neutrino mass hierarchy and are in agreements with the recent data.


title: On the computation of intersection numbers for twisted cocycles

abstract: Intersection numbers of twisted cocycles arise in mathematics in the field of algebraic geometry. Quite recently, they appeared in physics: Intersection numbers of twisted cocycles define a scalar product on the vector space of Feynman integrals. With this application, the practical and efficient computation of intersection numbers of twisted cocycles becomes a topic of interest. An existing algorithm for the computation of intersection numbers of twisted cocycles requires in intermediate steps the introduction of algebraic extensions (for example square roots), although the final result may be expressed without algebraic extensions. In this article I present an improvement of this algorithm, which avoids algebraic extensions.


title: Superheavy dark matter in $R+R^2$ cosmology with conformal anomaly

abstract: Cosmological evolution and particle creation in $R^2$-modified gravity are considered for the case of the dominant decay of the scalaron into a pair of gauge bosons due to conformal anomaly. It is shown that in the process of thermalization superheavy dark matter with the coupling strength typical for the GUT SUSY can be created. Such dark matter would have the proper cosmological density if the particle mass is close to $10^{12}$ GeV.


title: Hybrid Pass Transistor Logic with Ambipolar Transistors

abstract: In comparison to the conventional complementary pull-up and pull-down logic structure, the pass transistor logic (PTL) family reduces the number of transistors required to perform logic functions, thereby reducing both area and power consumption. However, this logic family requires inter-stage inverters to ensure signal integrity in cascaded logic circuits, and inverters must be used to provide each logical input signal in its complementary form. These inverters and complementary signals increase the device count and significantly degrade overall system efficiency. Dual-gate ambipolar field-effect transistors natively provide a single-transistor XNOR operation and permit highly-efficient and compact circuits due to their ambipolar capabilities. Similar to PTL, logic circuits based on ambipolar field-effect transistors require complementary signals. Therefore, numerous inverters are required, with significant energy and area costs. Ambipolar field-effect transistors are a natural match for PTL, as hybrid ambipolar-PTL circuits can simultaneously use these inverters to satisfy their necessity in both PTL and ambipolar circuits. We therefore propose a new hybrid ambipolar-PTL logic family that exploits the compact logic of PTL and the ambipolar capabilities of ambipolar field-effect transistors. Novel hybrid ambipolar-PTL circuits were designed and simulated in SPICE, demonstrating strong signal integrity along with the efficiency advantages of using the required inverters to simultaneously satisfy the requirements of PTL and ambipolar circuits. In comparison to the ambipolar field-effect transistors in the conventional CMOS logic structure, our hybrid full adder circuit can reduce propagation delay by 47%, energy consumption by 88%, energy-delay product by a factor of 9, and area-energy-delay product by a factor of 20.


title: Dust entrainment in galactic winds

abstract: Winds driven by stellar feedback are an essential part of the galactic ecosystem and are the main mechanism through which low-mass galaxies regulate their star formation. These winds are generally observed to be multi-phase with detections of entrained neutral and molecular gas. They are also thought to enrich the circum-galactic medium around galaxies with metals and dust. This ejected dust encodes information about the integrated star formation and outflow history of the galaxy. It is therefore, important to understand how much dust is entrained and driven out of the disc by galactic winds. Here we demonstrate that stellar feedback is efficient in driving dust-enriched winds and eject enough material to account for the amount of extraplanar dust observed in nearby galaxies. The amount of dust in the wind depends on the sites from where they are launched, with dustier galaxies launching more dust enriched outflows. Moreover, the outflowing cold-dense gas is significantly more dust-enriched than the volume filling hot tenuous material, naturally reproducing the complex multiphase structure of the outflowing wind observed in nearby galaxies. These results provide an important new insight into the dynamics, structure, and composition of galactic winds and their role in determining the dust content of the extragalactic gas in galaxies.


title: Scatter in Sunyaev--Zel'dovich effect scaling relations explained by inter-cluster variance in mass accretion histories

abstract: X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev--Zel'dovich effect observable-mass scaling relation ($Y_\mathrm{SZ}-M$) due solely to inter-cluster variation in mass accretion histories. We study the dependence of $Y_\mathrm{SZ}-M$ on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9\%$ scatter in $Y_\mathrm{SZ}-M$ at $z=0$, increasing as the aperture used to compute $Y_\mathrm{SZ}$ increases from $R_\mathrm{500c}$ to $5R_\mathrm{500c}$. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in $Y_\mathrm{SZ}-M$. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, $Y_\mathrm{SZ}-M$ deviates more from self-similarity and scatter increases. We show that the $Y_\mathrm{SZ}-M$ residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.


title: Hyper-optimized tensor network contraction

abstract: Tensor networks represent the state-of-the-art in computational methods across many disciplines, including the classical simulation of quantum many-body systems and quantum circuits. Several applications of current interest give rise to tensor networks with irregular geometries. Finding the best possible contraction path for such networks is a central problem, with an exponential effect on computation time and memory footprint. In this work, we implement new randomized protocols that find very high quality contraction paths for arbitrary and large tensor networks. We test our methods on a variety of benchmarks, including the random quantum circuit instances recently implemented on Google quantum chips. We find that the paths obtained can be very close to optimal, and often many orders or magnitude better than the most established approaches. As different underlying geometries suit different methods, we also introduce a hyper-optimization approach, where both the method applied and its algorithmic parameters are tuned during the path finding. The increase in quality of contraction schemes found has significant practical implications for the simulation of quantum many-body systems and particularly for the benchmarking of new quantum chips. Concretely, we estimate a speed-up of over 10,000$\times$ compared to the original expectation for the classical simulation of the Sycamore `supremacy' circuits.


title: The frequency of gaseous debris discs around white dwarfs

abstract: 1-3 per cent of white dwarfs are orbited by planetary dusty debris detectable as infrared emission in excess above the white dwarf flux. In a rare subset of these systems, a gaseous disc component is also detected via emission lines of the Ca II 8600\r{A} triplet, broadened by the Keplerian velocity of the disc. We present the first statistical study of the fraction of debris discs containing detectable amounts of gas in emission at white dwarfs within a magnitude and signal-to-noise limited sample. We select 7705 single white dwarfs spectroscopically observed by the Sloan Digital Sky Survey (SDSS) and $Gaia$ with magnitudes $g$ $\leq$ 19. We identify five gaseous disc hosts, all of which have been previously discovered. We calculate the occurrence rate of a white dwarf hosting a debris disc detectable via Ca II emission lines as 0.067$\pm$$_{0.025}^{0.042}$ per cent. This corresponds to an occurrence rate for a dusty debris disc to have an observable gaseous component in emission as 4$\pm$$_{2}^{4}$ per cent. Given that variability is a common feature of the emission profiles of gaseous debris discs, and the recent detection of a planetesimal orbiting within the disc of SDSSJ122859.93+104032.9, we propose that gaseous components are tracers for the presence of planetesimals embedded in the discs and outline a qualitative model. We also present spectroscopy of the Ca II triplet 8600\r{A} region for 20 white dwarfs hosting dusty debris discs in an attempt to identify gaseous emission. We do not detect any gaseous components in these 20 systems, consistent with the occurrence rate that we calculated.


title: Universality and Non-Universality in Distributed Nuclear Burning in Homogeneous Isotropic Turbulence

abstract: Nuclear burning plays a key role in a wide range of astrophysical stellar transients, including thermonuclear, pair instability, and core collapse supernovae, as well as kilonovae and collapsars. Turbulence is now understood to also play a key role in these astrophysical transients. Here we demonstrate that turbulent nuclear burning may lead to large enhancements above the uniform background burning rate, since turbulent dissipation gives rise to temperature fluctuations, and in general the nuclear burning rates are highly sensitive to temperature. We derive results for the turbulent enhancement of the nuclear burning rate under the influence of strong turbulence in the distributed burning regime in homogeneous isotropic turbulence, using probability distribution function (PDF) methods. We demonstrate that the turbulent enhancement obeys a universal scaling law in the limit of weak turbulence. We further demonstrate that, for a wide range of key nuclear reactions, such as C$^{12}$(O$^{16}$, $\alpha$)Mg$^{24}$ and triple-$\alpha$, even relatively modest temperature fluctuations, of the order ten percent, can lead to enhancements of 1 - 3 orders of magnitude in the turbulent nuclear burning rate. We verify the predicted turbulent enhancement directly against numerical simulations, and find very good agreement. We also present an estimation for the onset of turbulent detonation initiation, and discuss implications of our results for the modeling of stellar transients.


title: First Results on Dark Matter Substructure from Astrometric Weak Lensing

abstract: Low-mass structures of dark matter (DM) are expected to be entirely devoid of light-emitting regions and baryons. Precisely because of this lack of baryonic feedback, small-scale substructures of the Milky Way are a relatively pristine testing ground for discovering aspects of DM microphysics and primordial fluctuations on subgalactic scales. In this work, we report results from the first search for Galactic DM subhalos with time-domain astrometric weak gravitational lensing. The analysis is based on a matched-filter template of local lensing corrections to the proper motion of stars in the Magellanic Clouds. We describe a data analysis pipeline detailing sample selection, background subtraction, and handling outliers and other systematics. For tentative candidate lenses, we identify a signature based on an anomalous parallax template that can unequivocally confirm the presence of a DM lens, opening up prospects for robust discovery potential with full time-series data. We present our constraints on substructure fraction $f_l \lesssim 5$ at 90% CL (and $f_l \lesssim 2$ at 50% CL) for compact lenses with radii $r_l < 1\,\mathrm{pc}$, with best sensitivity reached for lens masses $M_l$ around $10^7$-$10^8\,M_\odot$. Parametric improvements are expected with future astrometric data sets; by end of mission, $Gaia$ could reach $f_l \lesssim 10^{-3}$ for these massive point-like objects, and be sensitive to lighter and/or more extended subhalos for $\mathcal{O}(1)$ substructure fractions.


title: Advection by large-scale spiral flows in galaxy clusters

abstract: The intracluster medium of a galaxy cluster often shows an extended quasi-spiral structure, accentuated by tangential discontinuities known as cold fronts (CFs). These discontinuities are thought to isolate between low-entropy, high-metallicity gas inside (i.e., below) the CF that was lifted from the center of the cluster by some radial factor $f_i$, and high-entropy, low-metallicity gas outside the CF that was pushed inward by a factor $f_o$. We find broad support for such a picture, by comparing the entropy and metallicity discontinuities with the respective azimuthal averages, using newly deprojected thermal profiles in clusters A2029, A2142, A2204 and Centaurus, supplemented by deprojected CFs from the literature. In particular, the mean advection factors $f_K$ and $f_Z$, inferred respectively from entropy and metallicity, strongly correlate ($\mathcal{R}=0.7^{+0.2}_{-0.3}$) with each other, consistent with large-scale advection. However, unlike sloshing simulations, in which the inside/outside phases are an inflow/outflow settling back to equilibrium after a violent perturbation, our results are more consistent with an outflow/inflow, with the fast, Mach $\mathcal{M}_i\sim0.8$ gas inside the CF being a rapidly heated or mixed outflow, probably originating from the cD galaxy, and gas outside the CF being a $\mathcal{M}_o\sim0.03$, slowly-cooling inflow. In particular, entropy indicates an outside advection factor $f_{Ko}\simeq 1.33\pm0.04$ that is approximately constant in all CFs, gauging the distance traversed by inflowing gas within a cooling time. In contrast, $1.1\lesssim f_{Ki}\lesssim 2.5$ and $1\lesssim f_Z\lesssim 17$ vary considerably among clusters, and strongly correlate ($3.1\sigma\text{-}4.2\sigma$) with the virial mass, $f_{Ki}\propto M_{200}^{0.14\pm0.07}$ and $f_Z\propto M_{200}^{1.4\pm0.4}$, suggesting that each cluster sustains a quasi-steady spiral flow.


title: A robust two-parameter description of the stellar profile of elliptical galaxies

abstract: The stellar density profile a galaxy is typically summarised with two numbers: total stellar mass and half-light radius. The total mass of a galaxy, however, is not a well-defined quantity, due to the finite depth of photometric observations and the arbitrariness of the distinction between galaxy and diffuse intra-group light. This limits our ability to make accurate comparisons between models and observations. I wish to provide a more robust two-parameter description of the stellar density distribution of elliptical galaxies, in terms of quantities that can be measured unambiguously. I propose to use the stellar mass enclosed within 10 kpc in projection, $M_{*,10}$, and the mass-weighted stellar density slope within the same aperture, $\Gamma_{*,10}$, for this purpose. I measured the distribution in $M_{*,10}$ and $\Gamma_{*,10}$ of a sample of elliptical galaxies from the Sloan Digital Sky Survey and the Galaxy And Mass Assembly survey, using photometry from the Hyper Suprime-Cam survey. The pair of values of $(M_{*,10},\Gamma_{*,10})$ can be used to predict the stellar density profile in the inner 10 kpc of a galaxy with better than 20% accuracy. Similarly, $M_{*,10}$ and $\Gamma_{*,10}$ can be combined to obtain a proxy for stellar velocity dispersion at least as good as the stellar mass fundamental plane. As a first application, I then compared the observed $M_{*,10}-\Gamma_{*,10}$ relation of elliptical galaxies with that of similarly selected galaxies in the EAGLE Reference simulation. Observed and simulated galaxies match at $M_{*,10}=10^{11}M_\odot$, but the EAGLE $M_{*,10}-\Gamma_{*,10}$ relation is shallower and has a larger intrinsic scatter compared to observations. This new parameterisation of the stellar density profile of massive elliptical galaxies provides a more robust way of comparing results from different photometric surveys and from hydrodynamical simulations.


title: The formation of exponential disk galaxies in MOND

abstract: The formation and evolution of galaxies is highly dependent on the dynamics of stars and gas, which is governed by the underlying law of gravity. To investigate how the formation and evolution of galaxies takes place in Milgromian gravity (MOND), we present full hydrodynamical simulations with the Phantom of Ramses (POR) code. These are the first-ever galaxy formation simulations done in MOND with detailed hydrodynamics, including star formation, stellar feedback, radiative transfer and supernovae. These models start from simplified initial conditions, in the form of isolated, rotating gas spheres in the early Universe. These collapse and form late-type galaxies obeying several scaling relations, which was not a priori expected. The formed galaxies have a compact bulge and a disk with exponentially decreasing surface mass density profiles and scale lengths consistent with observed galaxies, and vertical stellar mass distributions with distinct exponential profiles (thin and thick disk). This work thus shows for the first time that disk galaxies with exponential profiles in both gas and stars are a generic outcome of collapsing gas clouds in MOND. These models have a slight lack of stellar angular momentum because of their somewhat compact stellar bulge, which is connected to the simple initial conditions and the negligible later gas accretion. We also analyse how the addition of more complex baryonic physics changes the main resulting properties of the models and find this to be negligibly so in the Milgromian framework.


title: Sound propagation in a two-dimensional Bose gas across the superfluid transition

abstract: Motivated by recent experiments in Phys. Rev. Lett. 121, 145301 (2018), we study sound propagation in a two-dimensional (2D) Bose gas across the superfluid-thermal transition using classical field dynamics. Below the transition temperature we find a Bogoliubov and a non-Bogoliubov mode, above it we find the normal sound mode and the diffusive mode, as we determine from the dynamical structure factor. Our simulations of the experimental procedure agree with the measured velocities, and show that below the transition temperature the measurements detect the Bogoliubov mode. Above the transition, they either detect the normal sound mode for low densities or weak interactions, or the diffusive mode for high densities or strong interactions. As a key observation, we discuss the weak coupling regime in which the non-Bogoliubov mode has a higher velocity than the Bogoliubov mode, in contrast to a hydrodynamic scenario. We propose to detect this regime via step-pulse density perturbation, which simultaneously detects both sound modes


title: Two intermediate-mass transiting brown dwarfs from the TESS mission

abstract: We report the discovery of two intermediate-mass brown dwarfs (BDs), TOI-569b and TOI-1406b, from NASA's Transiting Exoplanet Survey Satellite mission. TOI-569b has an orbital period of $P = 6.55604 \pm 0.00016$ days, a mass of $M_b = 64.1 \pm 1.9 M_J$, and a radius of $R_b = 0.75 \pm 0.02 R_J$. Its host star, TOI-569, has a mass of $M_\star = 1.21 \pm 0.03 M_\odot$, a radius of $R_\star = 1.47 \pm 0.03 R_\odot$, $\rm [Fe/H] = +0.29 \pm 0.09$ dex, and an effective temperature of $T_{\rm eff} = 5768 \pm 110K$. TOI-1406b has an orbital period of $P = 10.57415 \pm 0.00063$ days, a mass of $M_b =46.0 \pm 2.7 M_J$, and a radius of $R_b = 0.86 \pm 0.03 R_J$. The host star for this BD has a mass of $M_\star =1 .18 \pm 0.09 M_\odot$, a radius of $R_\star = 1.35 \pm 0.03 R_\odot$, $ \rm [Fe/H] = -0.08 \pm 0.09$ dex and an effective temperature of $T_{\rm eff} = 6290 \pm 100K$. Both BDs are in circular orbits around their host stars and are older than 3 Gyr based on stellar isochrone models of the stars. TOI-569 is one of two slightly evolved stars known to host a transiting BD (the other being KOI-415). TOI-1406b is one of three known transiting BDs to occupy the mass range of $40-50 M_J$ and one of two to have a circular orbit at a period near 10 days (with the first being KOI-205b).Both BDs have reliable ages from stellar isochrones in addition to their well-constrained masses and radii, making them particularly valuable as tests for substellar isochrones in the BD mass-radius diagram.


title: Electrical probes of the non-Abelian spin liquid in Kitaev materials

abstract: Recent thermal-conductivity measurements evidence a magnetic-field-induced non-Abelian spin liquid phase in the Kitaev material $\alpha$-$\mathrm{RuCl}_{3}$. Although the platform is a good Mott insulator, we propose experiments that electrically probe the spin liquid's hallmark chiral Majorana edge state and bulk anyons, including their exotic exchange statistics. We specifically introduce circuits that exploit interfaces between electrically active systems and Kitaev materials to `perfectly' convert electrons from the former into emergent fermions in the latter---thereby enabling variations of transport probes invented for topological superconductors and fractional quantum Hall states. Along the way we resolve puzzles in the literature concerning interacting Majorana fermions, and also develop an anyon-interferometry framework that incorporates nontrivial energy-partitioning effects. Our results illuminate a partial pathway towards topological quantum computation with Kitaev materials.


title: Evolution of the grain size distribution in galactic discs

abstract: Dust is formed out of stellar material and is constantly affected by different mechanisms occurring in the ISM. Dust grains behave differently under these mechanisms depending on their sizes, and therefore the dust grain size distribution also evolves as part of the dust evolution itself. Following how the grain size distribution evolves is a difficult computing task that is just recently being overtaking. Smoothed particle hydrodynamic (SPH) simulations of a single galaxy as well as cosmological simulations are producing the first predictions of the evolution of the dust grain size distribution. We compare for the first time the evolution of the dust grain size distribution predicted by the SPH simulations with the results provided by the observations. We analyse how the radial distribution of the small to large grain mass ratio (D(S)/D(L)) changes over the whole discs in three galaxies: M 101, NGC 628 and M 33. We find good agreement between the observed radial distribution of D(S)/D(L) and what is obtained from the SPH simulations of a single galaxy. The central parts of NGC 628, at high metallicity and with a high molecular gas fraction, are mainly affected not only by accretion but also by coagulation of dust grains. The centre of M 33, having lower metallicity and lower molecular gas fraction, presents an increase of D(S)/D(L), showing that shattering is very effective in creating a large fraction of small grains. Observational results provided by our galaxies confirm the general relations predicted by the cosmological simulations based on the two grain size approximation. However, we present evidence that the simulations could be overestimating the amount of large grains in high massive galaxies.


title: Discovery of a low-mass companion embedded in the disk of the young massive star MWC 297 with VLT/SPHERE

abstract: We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics.


title: Gauge Theories for the Thermal Hall Effect

abstract: We consider the thermal Hall effect of fermionic matter coupled to emergent gauge fields in 2+1 dimensions. While the low-temperature thermal Hall conductivity of bulk topological phases can be connected to chiral edge states and a gravitational anomaly, there is no such interpretation at nonzero temperatures above 2+1 dimensional quantum critical points. In the limit of a large number of matter flavors, the leading contribution to the thermal Hall conductivity is that from the fermionic matter. The next-to-leading contribution is from the gauge fluctuations, and this has a sign which is opposite to that of the matter contribution. We illustrate this by computations on a Dirac Chern-Simons theory of the quantum phase transition in a square-lattice antiferromagnet involving the onset of semion topological order. We find similar results for a model of the pseudogap metal with Fermi pockets coupled to an emergent U(1) gauge field. We note connections to recent observations on the hole-doped cuprates: our theory captures the main trends, but the overall magnitude of the effect is smaller than that observed.


title: Spontaneous scalarization of charged Reissner-Nordstr\"om black holes: Analytic treatment along the existence line

abstract: It has recently been demonstrated that charged black holes can support spatially regular matter configurations made of massless scalar fields which are non-minimally coupled to the electromagnetic field of the charged spacetime. Intriguingly, using numerical techniques, it has been revealed that the resonant spectra of the composed charged-black-hole-nonminimally-coupled-scalar-field configurations are characterized by charge-dependent discrete scalarization bands $\alpha\in\{[\alpha^{-}_{n}({\bar Q}),\alpha^{+}_{n}({\bar Q}]\}_{n=0}^{n=\infty}$, where $\alpha$ is the dimensionless coupling constant of the theory and ${\bar Q}\equiv Q/M$ is the dimensionless charge-to-mass ratio of the central supporting black hole. In the present paper we use {\it analytical} techniques in order to study the physical and mathematical properties of the spatially regular non-minimally coupled scalar field configurations (linearized scalar clouds) which are supported by the central charged Reissner-Nordstr\"om black holes. In particular, we derive a remarkably compact formula for the discrete resonant spectrum $\{\alpha^-_n({\bar Q})\}_{n=0}^{n=\infty}$ which characterizes the composed black-hole-linearized-field configurations along the {\it existence-line} of the system, the critical line which separates bare Reissner-Nordstr\"om black holes from hairy scalarized black-hole configurations. The analytical results are confirmed by direct numerical computations.


title: Rational CFT With Three Characters: The Quasi-Character Approach

abstract: Quasi-characters are vector-valued modular functions having an integral, but not necessarily positive, q-expansion. Using modular differential equations, a complete classification has been provided in arXiv:1810.09472 for the case of two characters. These in turn generate all possible admissible characters, of arbitrary Wronskian index, in order two. Here we initiate a study of the three-character case. We conjecture several infinite families of quasi-characters and show in examples that their linear combinations an generate admissible characters with arbitrarily large Wronskian index. The structure is completely different from the order two case, and the novel coset construction of arXiv:1602.01022 plays a key role in discovering the appropriate families. Using even unimodular lattices, we construct some explicit three-character CFT corresponding to the new admissible characters.

categories: hep-th cond-mat.str-el math-ph math.MP math.NT


title: Observational constraints on the optical and near-infrared emission from the neutron star-black hole binary merger S190814bv

abstract: On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.


title: Synthesizing three-body interaction of spin chirality with superconducting qubits

abstract: Superconducting qubits provide a competitive platform for quantum simulation of complex dynamics that lies at the heart of quantum many-body systems, because of the flexibility and scalability afforded by the nature of microfabrication. However, in a multiqubit device, the physical form of couplings between qubits is either an electric (capacitor) or magnetic field (inductor), and the associated quadratic field energy determines that only two-body interaction in the Hamiltonian can be directly realized. Here we propose and experimentally synthesize the three-body spin-chirality interaction in a superconducting circuit based on Floquet engineering. By periodically modulating the resonant frequencies of the qubits connected with each other via capacitors, we can dynamically turn on and off qubit-qubit couplings, and further create chiral flows of the excitations in the three-qubit circular loop. Our result is a step toward engineering dynamical and many-body interactions in multiqubit superconducting devices, which potentially expands the degree of freedom in quantum simulation tasks.


title: Impact of Residual Carbon Impurities and Gallium Vacancies on Trapping Effects in AlGaN/GaN MIS-HEMTs

abstract: Effects of residual C impurities and Ga vacancies on the dynamic instabilities of AlN/AlGaN/GaN metal insulator semiconductor high electron mobility transistors are investigated. Secondary ion mass spectroscopy, positron annihilation spectroscopy, steady state and time-resolved photoluminescence (PL) measurements have been performed in conjunction with electrical characterization and current transient analyses. The correlation between yellow luminescence (YL), C- and Ga vacancy concentration is investigated. Time-resolved PL indicating the C$_{\mathrm{N}}$O$_{\mathrm{N}}$ complex as the main source of the YL, while Ga vacancies or related complexes with C seem not to play a major role. The device dynamic performance is found to be significantly dependent on the C concentration close to the channel of the transistor. Additionally, the magnitude of the YL is found to be in agreement with the threshold voltage shift and with the on-resistance degradation. Trap analysis of the GaN buffer shows an apparent activation energy of $\sim$0.8eV for all samples, pointing to a common dominating trapping process and that the growth parameters affect solely the density of trap centres. It is inferred that the trapping process is likely to be directly related to C based defects.


title: BOFFIN TTS: Few-Shot Speaker Adaptation by Bayesian Optimization

abstract: We present BOFFIN TTS (Bayesian Optimization For FIne-tuning Neural Text To Speech), a novel approach for few-shot speaker adaptation. Here, the task is to fine-tune a pre-trained TTS model to mimic a new speaker using a small corpus of target utterances. We demonstrate that there does not exist a one-size-fits-all adaptation strategy, with convincing synthesis requiring a corpus-specific configuration of the hyper-parameters that control fine-tuning. By using Bayesian optimization to efficiently optimize these hyper-parameter values for a target speaker, we are able to perform adaptation with an average 30% improvement in speaker similarity over standard techniques. Results indicate, across multiple corpora, that BOFFIN TTS can learn to synthesize new speakers using less than ten minutes of audio, achieving the same naturalness as produced for the speakers used to train the base model.


title: Review of Pulsar Timing Array for Gravitational Wave Research

abstract: Ongoing research on Pulsar Timing Array (PTA) to detect gravitational radiation is reviewed. Here, we discuss the use of millisecond pulsars as a gravitational wave detector, the sources of gravitational radiation detectable by PTAs and the current status of PTA experiments pointing out the future possibilities.


title: Dropout Prediction over Weeks in MOOCs by Learning Representations of Clicks and Videos

abstract: This paper addresses a key challenge in MOOC dropout prediction, namely to build meaningful representations from clickstream data. While a variety of feature extraction techniques have been explored extensively for such purposes, to our knowledge, no prior works have explored modeling of educational content (e.g. video) and their correlation with the learner's behavior (e.g. clickstream) in this context. We bridge this gap by devising a method to learn representation for videos and the correlation between videos and clicks. The results indicate that modeling videos and their correlation with clicks bring statistically significant improvements in predicting dropout.


title: Constraining the Mass of the Emerging Galaxy Cluster SpARCS1049+56 at z=1.71 with Infrared Weak Lensing

abstract: In the hierarchical structure formation model of the universe, galaxy clusters are assembled through a series of mergers. Accordingly, it is expected that galaxy clusters in the early universe are actively forming and dynamically young. Located at a high redshift of z=1.71, SpARCS1049+56 offers a unique look into the galaxy cluster formation process. This cluster has been shown to be rich in cluster galaxies and to have intense star formation. Its high redshift pushes a weak-lensing analysis beyond the regime of the optical spectrum into that of the infrared. Equipped with deep Hubble Space Telescope Wide Field Camera 3 UVIS and IR observations, we present a weak-lensing characterization of SpARCS1049+56. As few IR weak-lensing studies have been performed, we discuss the details of PSF modeling and galaxy shape measurement for an IR weak-lensing procedure and the systematics that come with the territory. It will be critical to understand these systematics in future weak-lensing studies in the IR with the next generation space telescopes such as JWST, Euclid, and WFIRST. Through a careful analysis, the mass distribution of this young galaxy cluster is mapped and the convergence peak is detected at a 3.3 sigma level. The weak-lensing mass of the cluster is estimated to be $3.5\pm1.2\times10^{14}\ \text{M}_\odot$ and is consistent with the mass derived from a mass-richness scaling relation. This mass is extreme for a cluster at such a high redshift and suggests that SpARCS1049+56 is rare in the standard $\Lambda$CDM universe.


title: On indicated coloring of lexicographic product of graphs

abstract: Indicated coloring is a graph coloring game in which two players collectively color the vertices of a graph in the following way. In each round the first player (Ann) selects a vertex, and then the second player (Ben) colors it properly, using a fixed set of colors. The goal of Ann is to achieve a proper coloring of the whole graph, while Ben is trying to prevent the realization of this project. The smallest number of colors necessary for Ann to win the game on a graph $G$ (regardless of Ben's strategy) is called the indicated chromatic number of $G$, denoted by $\chi_i(G)$. In this paper, we have shown that for any graphs $G$ and $H$, $G[H]$ is $k$-indicated colorable for all $k\geq\mathrm{col}(G)\mathrm{col}(H)$. Also, we have shown that for any graph $G$ and for some classes of graphs $H$ with $\chi(H)=\chi_i(H)=\ell$, $G[H]$ is $k$-indicated colorable if and only if $G[K_\ell]$ is $k$-indicated colorable. As a consequence of this result we have shown that for some particular families of graphs $G$ and $H$, $G[H]$ is $k$-indicated colorable for every $k\geq \chi(G[H])$. This serves as a partial answer to one of the questions raised by A. Grzesik in \cite{and}. In addition, if $G$ is a Bipartite graph or a $\{P_5,K_3\}$-free graph (or) a $\{P_5,Paw\}$-free graph and if $H$ is from the same families of graphs, then we have shown that $\chi_i(G[H])=\chi(G[H])$.


title: Magnetically controlled vector based on E coli Nissle 1917

abstract: In this paper, it is first discovered that EcN has natural magnetically controlled properties, that is, E coli Nissle 1917 cells move in a gradient magnetic field of permanent magnets without artificial magnetic labeling. It is also shown in the paper that the cultivation of EcN in a medium enriched with iron chelates and under the influence of an external magnetic field increases several times the magnetophoretic mobility of E coli Nissle 1917 cells compared to cultivation under standard conditions. Thus, the speed of movement of E coli Nissle 1917 cells in a gradient magnetic field of laboratory magnets is achieved of the order of several mm/s. It has also been shown for the first time that the cultivation of E coli Nissle 1917 biomass is accelerated by cultivation in an external magnetic field, while the change in the concentration of iron chelates in the medium has no effect on the cultivation dynamics of the E coli Nissle 1917 culture.


title: On the inverse best approximation property of systems of subspaces of a Hilbert space

abstract: Let $H$ be a Hilbert space and $H_1,...,H_n$ be closed subspaces of $H$. Denote by $P_k$ the orthogonal projection onto $H_k$, $k=1,2,...,n$. Following Patrick L. Combettes and Noli N. Reyes, we will say that the system of subspaces $H_1,...,H_n$ possesses the inverse best approximation property (IBAP) if for arbitrary elements $x_1\in H_1,...,x_n\in H_n$ there exists an element $x\in H$ such that $P_k x=x_k$ for all $k=1,2,...,n$. We provide various new necessary and sufficient conditions for a system of $n$ subspaces to possess the IBAP. Using the main characterization theorem, we study properties of the systems of subspaces which possess the IBAP, obtain a sufficient condition for a system of subspaces to possess the IBAP, and provide examples of systems of subspaces which possess the IBAP. These results are applied to a problem of probability theory. Let $(\Omega,\mathcal{F},\mu)$ be a probability space and $\mathcal{F}_1,...,\mathcal{F}_n$ be sub-$\sigma$-algebras of $\mathcal{F}$. We will say that the collection $\mathcal{F}_1,...,\mathcal{F}_n$ possesses the inverse marginal property (IMP) if for arbitrary random variables $\xi_1,...,\xi_n$ such that (1) $\xi_k$ is $\mathcal{F}_k$-measurable, $k=1,2,...,n$; (2) $E|\xi_k|^2<\infty$, $k=1,2,...,n$; (3) $E\xi_1=E\xi_2=...=E\xi_n$, there exists a random variable $\xi$ such that $E|\xi|^2<\infty$ and $E(\xi|\mathcal{F}_k)=\xi_k$ for all $k=1,2,...,n$. We will show that a collection of sub-$\sigma$-algebras possesses the IMP if and only if the system of corresponding marginal subspaces possesses the IBAP. We consider two examples; in the first example $\Omega=\mathbb{N}$, in the second example $\Omega=[a,b)$. For these examples we establish relations between the IMP, the IBAP, closedness of the sum of marginal subspaces and "fast decreasing" of tails of the measure $\mu$.


title: A Domain Semantics for Higher-Order Recursive Processes

abstract: The polarized SILL programming language uniformly integrates functional programming and session-typed message-passing concurrency. It supports general recursion, asynchronous and synchronous communication, and higher-order programs that communicate channels and processes. We give polarized SILL a domain-theoretic semantics---the first denotational semantics for a language with this combination of features. Session types in polarized SILL denote pairs of domains of unidirectional communications. Processes denote continuous functions between these domains, and process composition is interpreted by a trace operator. We illustrate our semantics by validating expected program equivalences.


title: Size of bulk fermions in the SYK model

abstract: The study of quantum gravity in the form of the holographic duality has uncovered and motivated the detailed investigation of various diagnostics of quantum chaos. One such measure is the operator size distribution, which characterizes the size of the support region of an operator and its evolution under Heisenberg evolution. In this work, we examine the role of the operator size distribution in holographic duality for the Sachdev-Ye-Kitaev (SYK) model. Using an explicit construction of AdS$_2$ bulk fermion operators in a putative dual of the low temperature SYK model, we study the operator size distribution of the boundary and bulk fermions. Our result provides a direct derivation of the relationship between (effective) operator size of both the boundary and bulk fermions and bulk $\text{SL}(2; \mathbb{R})$ generators.


title: A Posteriori Error Estimates for Self-Similar Solutions to the Euler Equations

abstract: The main goal of this paper is to analyze a family of "simplest possible" initial data for which, as shown by numerical simulations, the incompressible Euler equations have multiple solutions. We take here a first step toward a rigorous validation of these numerical results. Namely, we consider the system of equations corresponding to a self-similar solution, restricted to a bounded domain with smooth boundary. Given an approximate solution obtained via a finite dimensional Galerkin method, we establish a posteriori error bounds on the distance between the numerical approximation and the exact solution having the same boundary data.


title: Mutual Information-based State-Control for Intrinsically Motivated Reinforcement Learning

abstract: In reinforcement learning, an agent learns to reach a set of goals by means of an external reward signal. In the natural world, intelligent organisms learn from internal drives, bypassing the need for external signals, which is beneficial for a wide range of tasks. Motivated by this observation, we propose to formulate an intrinsic objective as the mutual information between the goal states and the controllable states. This objective encourages the agent to take control of its environment. Subsequently, we derive a surrogate objective of the proposed reward function, which can be optimized efficiently. Lastly, we evaluate the developed framework in different robotic manipulation and navigation tasks and demonstrate the efficacy of our approach. A video showing experimental results is available at https://youtu.be/CT4CKMWBYz0


title: Extending Machine Learning to Predict Unbalanced Physics Course Outcomes

abstract: Machine learning algorithms have recently been used to classify students as those likely to receive an A or B or students likely to receive a C, D, or F in a physics class. The performance metrics used in that study become unreliable when the outcome variable is substantially unbalanced. This study seeks to further explored the classification of students who will receive a C, D, and F and extend those methods to predicting whether a student will receive a D or F. The sample used for this work ($N=7184$) is substantially unbalanced with only 12\% of the students receiving a D or F. Applying the same methods as the previous study produced a classifier that was very inaccurate, classifying only 20\% of the D or F cases correctly. This study will focus on the random forest machine learning algorithm. By adjusting the random forest decision threshold, the correct classification rate of the D or F outcome rose to 46\%. This study also investigated the previous finding that demographic variables such as gender, underrepresented minority status, and first generation status had low variable importance for predicting class outcomes. Downsampling revealed that this was not the result of the underrepresentation of these students. An optimized classification model was constructed which predicted the D and F outcome with 46\% accuracy and C, D, and F outcome with 69\% accuracy; the accuracy of prediction of these outcomes is called "sensitivity" in the machine learning literature. Substantial variation was detected when this classification model was applied to predict the C, D, or F outcome for underrepresented demographic groups with 61\% sensitivity for women, 67\% for underrepresented minority students, and 78\% for first-generation students. Similar variation was observed for the D and F outcome.


title: Learning Probabilistic Intersection Traffic Models for Trajectory Prediction

abstract: Autonomous agents must be able to safely interact with other vehicles to integrate into urban environments. The safety of these agents is dependent on their ability to predict collisions with other vehicles' future trajectories for replanning and collision avoidance. The information needed to predict collisions can be learned from previously observed vehicle trajectories in a specific environment, generating a traffic model. The learned traffic model can then be incorporated as prior knowledge into any trajectory estimation method being used in this environment. This work presents a Gaussian process based probabilistic traffic model that is used to quantify vehicle behaviors in an intersection. The Gaussian process model provides estimates for the average vehicle trajectory, while also capturing the variance between the different paths a vehicle may take in the intersection. The method is demonstrated on a set of time-series position trajectories. These trajectories are reconstructed by removing object recognition errors and missed frames that may occur due to data source processing. To create the intersection traffic model, the reconstructed trajectories are clustered based on their source and destination lanes. For each cluster, a Gaussian process model is created to capture the average behavior and the variance of the cluster. To show the applicability of the Gaussian model, the test trajectories are classified with only partial observations. Performance is quantified by the number of observations required to correctly classify the vehicle trajectory. Both the intersection traffic modeling computations and the classification procedure are timed. These times are presented as results and demonstrate that the model can be constructed in a reasonable amount of time and the classification procedure can be used for online applications.


title: Decoupling of epitaxy related trapping effects in AlGaN/GaN metal-insulator semiconductor high electron mobility transistors

abstract: The decoupling of epitaxial factors influencing on the dynamic instabilities of AlGaN/GaN metal-insulator semiconductor high electron mobility transistors is investigated. Three different sets of samples have been analyzed by means of dynamic instabilities in the threshold voltage (V$_{\mathrm{th}}$ shift). Secondary ion mass spectroscopy, steady-state photoluminescence (PL) measurements have been performed in conjunction with electrical characterization. The device dynamic performance is found to be significantly dependent on both the C concentration close to the channel as well as on the distance between the channel and the higher doped C region. Additionally, we note that experiments studying trapping should avoid large variations in the sheet carrier density (N$_{\mathrm{s}}$). This change in the N$_{\mathrm{s}}$ itself has a significant impact on the V$_{\mathrm{th}}$ shift. This experimental trends are also supported by a basic model and device simulation. Finally, the relationship between the yellow luminescence (YL) and the band edge (BE) ratio and the V$_{\mathrm{th}}$ shift is investigated. As long as the basic layer structure is not changed, the YL/BE ratio obtained from steady-state PL is demonstrated to be a valid method in predicting trap concentrations in the GaN channel layer.


title: A note on Lie algebra cohomology

abstract: Given a finite dimensional Lie algebra $L$ let $I$ be the augmentation ideal in the universal enveloping algebra $U(L)$. We study the conditions on $L$ under which the $Ext$-groups $Ext (k,k)$ for the trivial $L$-module $k$ are the same when computed in the category of all $U(L)$-modules or in the category of $I$-torsion $U(L)$-modules. We also prove that the Rees algebra $\oplus _{n\geq 0}I^n$ is Noetherian if and only if $L$ is nilpotent. An application to cohomology of equivariant sheaves is given.


title: Avoidance of split overlaps

abstract: We generalize Axel Thue's familiar definition of overlaps in words, and show that there are no infinite words containing split occurrences of these generalized overlaps. Along the way we prove a useful theorem about repeated disjoint occurrences in words -- an interesting natural variation on the classical de Bruijn sequences.


title: Partially Observable Games for Secure Autonomy

abstract: Technology development efforts in autonomy and cyber-defense have been evolving independently of each other, over the past decade. In this paper, we report our ongoing effort to integrate these two presently distinct areas into a single framework. To this end, we propose the two-player partially observable stochastic game formalism to capture both high-level autonomous mission planning under uncertainty and adversarial decision making subject to imperfect information. We show that synthesizing sub-optimal strategies for such games is possible under finite-memory assumptions for both the autonomous decision maker and the cyber-adversary. We then describe an experimental testbed to evaluate the efficacy of the proposed framework.


title: Reply to: Overconfidence in Bayesian analyses of galaxy rotation curves

abstract: Cameron {\it et al.}~2019 (hereafter C19) recommends a more cautious and rigorous approach to statistical analysis in astronomy. We welcome this particular side of their communication as it helps stimulating the effort towards the adoption of better statistical methods in galaxy rotation curves, an effort to which we contributed with Rodrigues {\it et al.}~2018 (hereafter R18). Indeed, R18 was the first work that, in order to conclude on the universality of the acceleration scale $a_0$, studied the posterior distributions on $a_0$ of a large set of galaxies. As C19 agrees, the credible intervals were found within the Bayesian framework, that is, the marginalized posteriors on $a_0$ were found using Bayes' theorem to update the priors in light of the observational data; this process was done without introducing any approximation. Considering R18, C19 also remarks that: i) better methods to select the nuisance parameters and the corresponding priors could be used; ii) a quality cut based on $\chi^2$ values should not be used, and iii) the compatibility of the posteriors should be assessed in a more robust way. In the following, after first clarifying the context of our work, we address these criticisms.


title: Impossibility of convergence of a Heun function on the boundary of the disc of convergence

abstract: The Heun's equation is the Fuchsian equation of second order with four regular singularities. Heun functions generalize well-known special functions such as Spheroidal Wave, Lam\'{e}, Mathieu, hypergeometric-type functions, etc. The recursive relation of coefficients starts to appear by putting a power series into the Heun equation. A local Heun function consists of the three term recurrence relation in its power series, and we prove that the function is not convergent on the boundary of the disc of convergence.


title: Consistent Blandford-Znajek Expansion

abstract: The Blandford-Znajek mechanism is the continuous extraction of energy from a rotating black hole via plasma currents flowing on magnetic field lines threading the horizon. In the discovery paper, Blandford and Znajek demonstrated the mechanism by solving the equations of force-free electrodynamics in a perturbative expansion valid at small black hole spin. Attempts to extend this perturbation analysis to higher order have encountered inconsistencies.We overcome this problem using the method of matched asymptotic expansions, taking care to resolve all of the singular surfaces (light surfaces) in the problem. Working with the monopole field configuration, we show explicitly how the inconsistencies are resolved in this framework and calculate the field configuration to one order higher than previously known. However, there is no correction to the energy extraction rate at this order. These results confirm the basic consistency of the split monopole at small spin and lay a foundation for further perturbative studies of the Blandford-Znajek mechanism.


title: Exploring Maximum Entropy Distributions with Evolutionary Algorithms

abstract: This paper shows how to evolve numerically the maximum entropy probability distributions for a given set of constraints, which is a variational calculus problem. An evolutionary algorithm can obtain approximations to some well-known analytical results, but is even more flexible and can find distributions for which a closed formula cannot be readily stated. The numerical approach handles distributions over finite intervals. We show that there are two ways of conducting the procedure: by direct optimization of the Lagrangian of the constrained problem, or by optimizing the entropy among the subset of distributions which fulfill the constraints. An incremental evolutionary strategy easily obtains the uniform, the exponential, the Gaussian, the log-normal, the Laplace, among other distributions, once the constrained problem is solved with any of the two methods. Solutions for mixed ("chimera") distributions can be also found. We explain why many of the distributions are symmetrical and continuous, but some are not.

categories: stat.ME cs.IT cs.LG cs.NE math.IT stat.CO


title: The dynamics of the early universe in a model with radiation and a generalized Chaplygin gas

abstract: The early universe is modeled through the quantization of a Friedmann-Robertson-Walker model with positive curvature of the spatial hypersurfaces. In this model, the universe is filled by two fluids: radiation and a generalized Chaplygin gas. The quantization of this model is made following the prescriptions due to J. A. Wheeler and B. DeWitt. Using the Schutz's formalism, the time notion is recovered and the Wheeler-DeWitt equation transforms into a time dependent Schr\"{o}dinger equation, which rules the dynamics of the early universe, under the action of an effective potential $V_{eff}$. That potential, depends on three parameters. Depending on the values of these parameters, $V_{eff}$ may have two different shapes. $V_{eff}(a)$ may have the shape of a barrier or the shape of a well followed by a barrier. We solve, numerically, the appropriate time dependent Schr\"{o}dinger equation and obtain the time evolution of an initial wave function, for both cases. These wave functions satisfy suitable boundary conditions. For both shapes of $V_{eff}$, we compute the tunneling probability, which is a function of the mean kinetic energy associated to the radiation energy $E_m$ and of the three parameters of the generalized Chaplygin gas: $\alpha$, $A$ and $B$. The tunneling probabilities, for both shapes of $V_{eff}$, indicates that the universe should nucleate with the highest possible values of $E_m$, $\alpha$, $A$ and $B$. Finally, we study the classical universe evolution after the wavefunction has tunneled $V_{eff}$. The calculations show that the universe may emerge from the Planck era in an inflationary phase.


title: Brain Tumor Segmentation by Cascaded Deep Neural Networks Using Multiple Image Scales

abstract: Intracranial tumors are groups of cells that usually grow uncontrollably. One out of four cancer deaths is due to brain tumors. Early detection and evaluation of brain tumors is an essential preventive medical step that is performed by magnetic resonance imaging (MRI). Many segmentation techniques exist for this purpose. Low segmentation accuracy is the main drawback of existing methods. In this paper, we use a deep learning method to boost the accuracy of tumor segmentation in MR images. Cascade approach is used with multiple scales of images to induce both local and global views and help the network to reach higher accuracies. Our experimental results show that using multiple scales and the utilization of two cascade networks is advantageous.


title: Initial estimate for minimum energy pathways and transition states using velocities in internal coordinates

abstract: Many algorithms for finding reaction pathways require an initial estimate of the minimum energy path (MEP). Most estimation methods use a variational approach and thus must be seeded from an even simpler path, such as one generated by Cartesian interpolation. Often, care must be taken to avoid atomic intersections in this seed path, and the estimator itself may potentially converge to multiple undesirable local minima. As an alternative we form an initial estimate by numerically integrating a velocity vector field that is projected from redundant internal coordinates into the Cartesian manifold. We compare this method to the image dependent pair potential, the geodesic method, and linear Cartesian interpolation using three test cases: the rotation of a methyl group in ethane, HCN$\to$HNC tautomerization, and HONO elimination from dimethylnitramine. In the first test case, a zero-temperature string calculation seeded with our method converges to the MEP in significantly fewer geometry and SCF cycles than any of the others, while in the second, only the geodesic method slightly outperformed ours. In the third test case, we used the midpoint of each path as an initial guess for a transition state calculation. Our midpoint was geometrically the closest to the true transition state and converged in the fewest geometry and SCF cycles.


title: Time-dependence in non-Hermitian quantum systems

abstract: We present a coherent and consistent framework for explicit time-dependence in non-Hermitian quantum mechanics. The area of non-Hermitian quantum mechanics has been growing rapidly over the past twenty years. This has been driven by the fact that $\mathcal{PT}$-symmetric non-Hermitian systems exhibit real energy eigenvalues and unitary time evolution. Historically, the introduction of time into the world of non-Hermitian quantum mechanics has been a conceptually difficult problem to address, as it requires the Hamiltonian to become unobservable. We solve this issue with the introduction of a new observable energy operator and explain why its instigation is a necessary and natural progression in this setting. For the first time, the introduction of time has allowed us to make sense of the parameter regime in which the $\mathcal{PT}$-symmetry is spontaneously broken. Ordinarily, in the time-independent setting, the energy eigenvalues become complex and the wave functions are asymptotically unbounded. We demonstrate that in the time-dependent setting this broken symmetry can be mended and analysis on the spontaneously broken $\mathcal{PT}$ regime is indeed possible. We provide many examples of this mending on a wide range of different systems, beginning with a $2\times2$ matrix model and extending to higher dimensional matrix models and coupled harmonic oscillator systems with infinite Hilbert space. Furthermore, we use the framework to perform analysis on time-dependent quasi-exactly solvable models. We present the "eternal life" of entropy in this thesis. Ordinarily, for entangled quantum systems coupled to the environments, the entropy decays rapidly to zero. However, in the spontaneously broken regime, we find the entropy decays asymptotically to a non-zero value. We create an elegant framework for Darboux and Darboux/Crum transformations for time-dependent non-Hermitian Hamiltonians.


title: CP Violation in Same-sign Dilepton Production at the LHC

abstract: If the neutrino is a Majorana particle, low-energy lepton-number-violating (LNV) processes, such as neutrinoless double-beta ($0\nu\beta\beta$) decay, are possible. It may also be possible to observe high-energy $0\nu\beta\beta$-like LNV processes at the LHC. These are distinguished by the presence of same-sign dileptons in the final state (e.g., ${\bar u} d \to t {\bar b} \, e^- \mu^-$). In this paper, we show that CP-violating triple products (TPs) may be present in the process, and may be measurable at the LHC. If a nonzero TP were observed, it would give us much information about the underlying new physics (NP). We would know that there are (at least) two interfering NP amplitudes, with different weak phases and different Lorentz structures. And if we had some knowledge of the NP, e.g., by direct production of NP particles, we could get information about the magnitudes and relative phases of its couplings.


title: Natural orbitals in multiconfiguration calculations of hyperfine structure parameters

abstract: We are reinvestigating the hyperfine structure of sodium using a fully relativistic multiconfiguration approach. In the fully relativistic approach, the computational strategy somewhat differs from the original nonrelativistic counterpart used in J\"onsson et al. (Phys. Rev. A 53 (1996) 4021). Numerical instabilities force us to use a layer-by-layer approach that has some broad unexpected effects. Core correlation is found to be significant and therefore requires to be described in an adequate orbital basis. The natural-orbital basis provides an interesting alternative to the orbital basis from the layer-by-layer approach, allowing us to overcome some deficits of the latter, giving rise to magnetic dipole hyperfine structure constant values in excellent agreement with observations. Effort is made to assess the reliability of the natural-orbital bases and to illustrate their efficiency.


title: Fast Probabilistic Hosting Capacity Analysis for Active Distribution Systems

abstract: Interconnection studies for distributed energy resources (DERs) can currently take months since they entail simulating a large number of power flow scenarios. If DERs are to be actively controlled, probabilistic hosting capacity analysis (PHCA) studies become more time-consuming since they require solving multiple optimal power flow (OPF) tasks. PHCA is expedited here by leveraging the powerful tool of multiparametric programming (MPP). Using an approximate grid model, optimal DER setpoints are decided by a quadratic program, which depends on analysis and uncertain parameters in a possibly nonlinear fashion. By reformulating this program, feasible and infeasible OPF instances alike are handled in a unified way to uniquely reveal the location, frequency, and severity of feeder constraint violations. The effect of voltage regulators is also captured by novel approximate models. Upon properly extending MPP to PHCA, we were able to find the exact minimizers for 518,400 OPF instances on the IEEE 123-bus feeder by solving only 6,905 of them, and 86,400 instances on a 1,160-bus feeder by solving only 2,111 instances. This accelerated PHCA by a factor of 10. Thus, a utility can promptly infer grid statistics using real-world data without a probabilistic characterization of uncertain parameters.


title: Parallel 3DPIFCM Algorithm for Noisy Brain MRI Images

abstract: In this paper we implemented the algorithm we developed in [1] called 3DPIFCM in a parallel environment by using CUDA on a GPU. In our previous work we introduced 3DPIFCM which performs segmentation of images in noisy conditions and uses particle swarm optimization for finding the optimal algorithm parameters to account for noise. This algorithm achieved state of the art segmentation accuracy when compared to FCM (Fuzzy C-Means), IFCMPSO (Improved Fuzzy C-Means with Particle Swarm Optimization), GAIFCM (Genetic Algorithm Improved Fuzzy C-Means) on noisy MRI images of an adult Brain. When using a genetic algorithm or PSO (Particle Swarm Optimization) on a single machine for optimization we witnessed long execution times for practical clinical usage. Therefore, in the current paper our goal was to speed up the execution of 3DPIFCM by taking out parts of the algorithm and executing them as kernels on a GPU. The algorithm was implemented using the CUDA [13] framework from NVIDIA and experiments where performed on a server containing 64GB RAM , 8 cores and a TITAN X GPU with 3072 SP cores and 12GB of GPU memory. Our results show that the parallel version of the algorithm performs up to 27x faster than the original sequential version and 68x faster than GAIFCM algorithm. We show that the speedup of the parallel version increases as we increase the size of the image due to better utilization of cores in the GPU. Also, we show a speedup of up to 5x in our Brainweb experiment compared to other generic variants such as IFCMPSO and GAIFCM.


title: Multimodal fusion of imaging and genomics for lung cancer recurrence prediction

abstract: Lung cancer has a high rate of recurrence in early-stage patients. Predicting the post-surgical recurrence in lung cancer patients has traditionally been approached using single modality information of genomics or radiology images. We investigate the potential of multimodal fusion for this task. By combining computed tomography (CT) images and genomics, we demonstrate improved prediction of recurrence using linear Cox proportional hazards models with elastic net regularization. We work on a recent non-small cell lung cancer (NSCLC) radiogenomics dataset of 130 patients and observe an increase in concordance-index values of up to 10%. Employing non-linear methods from the neural network literature, such as multi-layer perceptrons and visual-question answering fusion modules, did not improve performance consistently. This indicates the need for larger multimodal datasets and fusion techniques better adapted to this biological setting.


title: Microscopic Origin of Spin-Orbit Torque in Ferromagnetic Heterostructures: A First Principles Approach

abstract: We present an {\it ab initio}-based theoretical framework which elucidates the origin of the spin-orbit torque (SOT) in Normal-Metal(NM)/Ferromagnet(FM) heterostructures. The SOT is decomposed into two contributions, namely, {\it spin-Hall} and the {\it spin-orbital} components. We find that {\it (i)} the Field-Like (FL) SOT is dominated by the spin-orbital component and {\it (ii)} both components contribute to the damping-like torque with comparable magnitude in the limit of thick Pt film. The contribution of the spin-orbital component to the DL-SOT is present only for NMs with strong SOC coupling strength. We demonstrate that the FL-SOT can be expressed in terms of the non-equilibrium spin-resolved orbital moment accumulation. The calculations reveal that the experimentally reported oxygen-induced sign-reversal of the FL-SOT in Pt/Co bilayers is due to the significant reduction of the majority-spin orbital moment accumulation on the interfacial NM atoms.


title: UNCC Biomedical Semantic Question Answering Systems. BioASQ: Task-7B, Phase-B

abstract: In this paper, we detail our submission to the 2019, 7th year, BioASQ competition. We present our approach for Task-7b, Phase B, Exact Answering Task. These Question Answering (QA) tasks include Factoid, Yes/No, List Type Question answering. Our system is based on a contextual word embedding model. We have used a Bidirectional Encoder Representations from Transformers(BERT) based system, fined tuned for biomedical question answering task using BioBERT. In the third test batch set, our system achieved the highest MRR score for Factoid Question Answering task. Also, for List type question answering task our system achieved the highest recall score in the fourth test batch set. Along with our detailed approach, we present the results for our submissions, and also highlight identified downsides for our current approach and ways to improve them in our future experiments.


title: 3DPIFCM Novel Algorithm for Segmentation of Noisy Brain MRI Images

abstract: We present a novel algorithm named 3DPIFCM, for automatic segmentation of noisy MRI Brain images. The algorithm is an extension of a well-known IFCM (Improved Fuzzy C-Means) algorithm. It performs fuzzy segmentation and introduces a fitness function that is affected by proximity of the voxels and by the color intensity in 3D images. The 3DPIFCM algorithm uses PSO (Particle Swarm Optimization) in order to optimize the fitness function. In addition, the 3DPIFCM uses 3D features of near voxels to better adjust the noisy artifacts. In our experiments, we evaluate 3DPIFCM on T1 Brainweb dataset with noise levels ranging from 1% to 20% and on a synthetic dataset with ground truth both in 3D. The analysis of the segmentation results shows a significant improvement in the segmentation quality of up to 28% compared to two generic variants in noisy images and up to 60% when compared to the original FCM (Fuzzy C-Means).


title: On the validation of Newcomb-Benford law and Weibull distribution in neuromuscular transmission

abstract: The neuromuscular junction represents a relevant substrate for revealing important biophysical mechanisms of synaptic transmission. In this context, calcium ions are important in the synapse machinery, providing the nervous impulse transmission to the muscle fiber. In this work, we carefully investigated whether intervals of spontaneous electrical activity, recorded in seven different calcium concentrations, conform Newcomb-Benford law. Our analysis revealed that electrical discharge of neuromuscular junction obeys the expected values for Newcomb-Benford law for first and second digits, while first-two digits do not perfectly follows the law. We next examined previous theoretical studies, establishing a relation between the law and lognormal and Weibull distributions. We showed that Weibull distribution is more appropriate to fit the intervals as compared to lognormal distribution. Altogether, the present findings strongly suggest that spontaneous activity is a base-scale invariant phenomenon.


title: Function approximation by neural nets in the mean-field regime: Entropic regularization and controlled McKean-Vlasov dynamics

abstract: We consider the problem of function approximation by two-layer neural nets with random weights that are "nearly Gaussian" in the sense of Kullback-Leibler divergence. Our setting is the mean-field limit, where the finite population of neurons in the hidden layer is replaced by a continuous ensemble. We show that the problem can be phrased as global minimization of a free energy functional on the space of (finite-length) paths over probability measures on the weights. This functional trades off the $L^2$ approximation risk of the terminal measure against the KL divergence of the path with respect to an isotropic Brownian motion prior. We characterize the unique global minimizer and examine the dynamics in the space of probability measures over weights that can achieve it. In particular, we show that the optimal path-space measure corresponds to the F\"ollmer drift, the solution to a McKean-Vlasov optimal control problem closely related to the classic Schr\"odinger bridge problem. While the F\"ollmer drift cannot in general be obtained in closed form, thus limiting its potential algorithmic utility, we illustrate the viability of the mean-field Langevin diffusion as a finite-time approximation under various conditions on entropic regularization. Specifically, we show that it closely tracks the F\"ollmer drift when the regularization is such that the minimizing density is log-concave.


title: Lozenge Tiling Function Ratios for Hexagons with Dents on Two Sides

abstract: We give a formula for the number of lozenge tilings of a hexagon on the triangular lattice with unit triangles removed from arbitrary positions along two non-adjacent, non-opposite sides. Our formula implies that for certain families of such regions, the ratios of their numbers of tilings are given by simple product formulas.


title: Quasi-one-dimensional system as a high-temperature superconductor

abstract: It is well-known that quasi-one-dimensional superconductors suffer from the pairing fluctuations that significantly reduce the superconducting temperature or even completely suppress any coherent behavior. Here we demonstrate that a coupling to a robust pair condensate changes the situation dramatically. In this case the quasi-one-dimensional system can be a high temperature superconductor governed by the proximity to the Lifshitz transition at which the Fermi level approaches the lower edge of the single-particle spectrum.


title: Coherent electronic transport in periodic crystals

abstract: We consider independent electrons in a periodic crystal in their ground state, and turn on a uniform electric field at some prescribed time. We rigorously define the current per unit volume and study its properties using both linear response and adiabatic theory. Our results provide a unified framework for various phenomena such as the quantification of Hall conductivity of insulators with broken time-reversibility, the ballistic regime of electrons in metals, Bloch oscillations in the long-time response of metals, and the static conductivity of graphene. We identify explicitly the regime in which each holds.