title: Numerical Relativity Simulations of the Neutron Star Merger GW170817: Long-Term Remnant Evolutions, Winds, Remnant Disks, and Nucleosynthesis

abstract: We present a systematic numerical-relativity study of the dynamical ejecta, winds and nucleosynthesis in neutron star merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOS) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamics turbulence up to 100 milliseconds postmerger. Spiral density waves propagating from the neutron star remnant to the disk trigger a wind with mass flux ${\sim}0.1{-}0.5\,{\rm M_\odot/s}$ persisting for the entire simulation as long as the remnant does not collapse to black hole. This wind has average electron fraction $\gtrsim 0.3$ and average velocity ${\sim}0.1-0.17\,$c and thus is a site for the production of weak $r$-process elements (mass number $A<195$). Disks around long-lived remnants have masses ${\sim}0.1{-}0.2\,{\rm M_\odot}$, temperatures peaking at $\lesssim10\,$MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multi-component kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the comparable-mass long-lived mergers robustly accounts for all the $r$-process peaks, from mass number ${\sim}75$ to actinides in terms of solar abundances. Total abundandes are weakly dependent on the EOS, while the mass ratio affect the production of first peak elements.


title: Debiased Galaxy Cluster Pressure Profiles from X-ray Observations and Simulations

abstract: We present an updated model for the average cluster pressure profile, adjusted for hydrostatic mass bias by combining results from X-ray observations with cosmological simulations. Our model estimates this bias by fitting a power-law to the relation between the "true" halo mass and X-ray cluster mass in hydrodynamic simulations (IllustrisTNG, BAHAMAS, and MACSIS). As an example application, we consider the REXCESS X-ray cluster sample and the Universal Pressure Profile (UPP) derived from scaled and stacked pressure profiles. We find adjusted masses, $M_\mathrm{500c},$ that are $\lesssim$15% higher and scaled pressures $P/P_\mathrm{500c}$ that have $\lesssim$35% lower normalization than previously inferred. Our Debiased Pressure Profile (DPP) is well-fit by a Generalized Navarro-Frenk-White (GNFW) function, with parameters $[P_0,c_{500},\alpha,\beta,\gamma]=[5.048,1.217,1.192,5.490,0.433]$ and does not require a mass-dependent correction term. When the DPP is used to model the Sunyaev-Zel'dovich (SZ) effect, we find that the integrated Compton $Y-M$ relation has only minor deviations from self-similar scaling. The thermal SZ angular power spectrum is lower in amplitude by approximately 30%, assuming nominal cosmological parameters (e.g. $\Omega_\text{m}=0.3$, $\sigma_8 = 0.8$), and is broadly consistent with recent Planck results without requiring additional bias corrections.


title: Mildly-Hierarchical triple dynamics and applications to the outer solar system

abstract: Three-body interactions are ubiquitous in astrophysics. For instance, Kozai-Lidov oscillations in hierarchical triple systems have been studied extensively and applied to a wide range of astrophysical systems. However, mildly-hierarchical triples also play an important role, but they are less explored. In this work we consider the secular dynamics of a test particle in a mildly-hierarchical configuration. We find the limit within which the secular approximation is reliable, present resonances and chaotic regions using surface of sections, and characterize regions of phase space that allow large eccentricity and inclination variations. Finally, we apply the secular results to the outer solar system. We focus on the distribution of extreme trans-neptunian objects (eTNOs) under the perturbation of a possible outer planet (Planet-9), and find that in addition to a low inclination Planet-9, a polar or a counter-orbiting one could also produce pericenter clustering of eTNOs, while the polar one leads to a wider spread of eTNO inclinations.


title: QQ-system and Weyl-type transfer matrices in integrable SO(2r) spin chains

abstract: We propose the full system of Baxter Q-functions (QQ-system) for the integrable spin chains with the symmetry of the $D_r$ Lie algebra. We use this QQ-system to derive new Weyl-type formulas expressing transfer matrices in all symmetric and antisymmetric (fundamental) representations through $r+1$ basic Q-functions. Our functional relations are consistent with the Q-operators proposed recently by one of the authors and verified explicitly on the level of operators at small finite length.


title: A Striped Electron Fluid on (111) KTaO$_3$

abstract: A recent study has revealed that the low carrier density electron gas (2DEG) induced at the interface of EuO and (111) KTaO$_3$ exhibits a broken symmetry phase with a strong in-plane anisotropy of the resistivity. We present a minimal tight binding model of this (111) 2DEG, including the large spin-orbit coupling from the Ta ions, which reveals a hexagonal Fermi surface with a highly enhanced 2$k_F$ electronic susceptibility. We argue that repulsive electronic interactions, together with a ferromagnetic EuO substrate, favor a magnetic stripe instability leading to a partially gapped Fermi surface. Such a stripe state, or its vestigial nematicity, could explain the observed transport anisotropy. We propose a $k\cdot p$ theory for the low energy $j=3/2$ states, which captures the key results from our tight-binding study, and further reveals the intertwined dipolar and octupolar modulations underlying this magnetic stripe order. We conclude by speculating on the relation of this stripe order to the superconductivity seen in this material.


title: Self-accelerating root search and optimisation methods based on rational interpolation

abstract: Iteration methods based on barycentric rational interpolation are derived that exhibit accelerating orders of convergence. For univariate root search, the derivative-free methods approach quadratic convergence and the first-derivative methods approach cubic convergence. For univariate optimisation, the convergence order of the derivative-free methods approaches 1.62 and the order of the first-derivative methods approaches 2.42. Generally, performance advantages are found with respect to low-memory iteration methods. In optimisation problems where the objective function and gradient is calculated at each step, the full-memory iteration methods converge asymptotically 1.8 times faster than the secant method. Frameworks for multivariate root search and optimisation are also proposed, though without discovery of practical parameter choices.


title: A Minimal Approach to Baryogenesis via Affleck-Dine and Inflaton Mass Terms

abstract: We present a minimal approach to the generation of the baryon ($B$) asymmetry of the Universe, in which the asymmetry is generated in a complex inflaton condensate via $B$-violating quadratic inflaton potential terms using the Affleck-Dine (AD) mechanism. We show that the $B$-violating quadratic mass terms create an oscillating asymmetry in the complex inflaton condensate at late times. The final asymmetry transferred to the Standard Model sector at reheating is naturally reduced to the magnitude of the observed $B$ asymmetry by the effect of averaging over the $B$ oscillations. This approach to baryogenesis can easily be realised in a wide range of inflation models. We also show that AD baryogenesis via inflaton mass terms is free of the tuning of field dynamics required by the conventional AD mechanism.


title: Exciting Prospects for Detecting Late-Time Neutrinos from Core-Collapse Supernovae

abstract: The importance of detecting neutrinos from a Milky Way core-collapse supernova is well known. An under-studied phase is proto-neutron star cooling. For SN 1987A, this seemingly began at about 2 s, and is thus probed by only 6 of the 19 events (and only the $\bar{\nu}_e$ flavor) in the Kamiokande-II and IMB detectors. With the higher statistics expected for present and near-future detectors, it should be possible to measure detailed neutrino signals out to very late times. We present the first comprehensive study of neutrino detection during the proto-neutron star cooling phase, considering a variety of outcomes, using all flavors, and employing detailed detector physics. For our nominal model, the event yields (at 10 kpc) after 10 s -- the approximate duration of the SN 1987A signal -- far exceed the entire SN 1987A yield, with $\simeq$250 $\bar{\nu}_e$ events (to 50 s) in Super-Kamiokande, $\simeq$110 $\nu_e$ events (to 40 s) in DUNE, and $\simeq$10 $\nu_\mu, \nu_\tau, \bar{\nu}_\mu, \bar{\nu}_\tau$ events (to 20 s) in JUNO. These data would allow unprecedented probes of the proto-neutron star, including the onset of neutrino transparency and hence its transition to a neutron star. If a black hole forms, even at very late times, this can be clearly identified. But will the detectors fulfill their potential for this perhaps once-ever opportunity for an all-flavor, high-statistics detection of a core collapse? Maybe. Further work is urgently needed, especially for DUNE to thoroughly investigate and improve its MeV capabilities.


title: Two-dipole and three-dipole dispersion coefficients for interaction of alkaline-earth atoms with alkaline-earth atoms and alkaline-earth ions

abstract: Apropos to the growing interest in the study of long-range interactions which for their applications in cold atom physics, we have performed theoretical calculation for the two-dipole $C_6$ and three-dipole $C_9$ dispersion coefficients involving alkaline-earth atoms with alkaline-earth atoms and alkaline-earth ions. The $C_6$ and $C_9$ coefficients are expressed in terms of the dynamic dipole polarizabilities, which are calculated using relativistic methods. Thereafter, the calculated $C_6$ coefficients for the considered alkaline-earth atoms among themselves are compared with the previously reported values. Due to unavailability of any other earlier theoretical or experimental results, for the $C_6$ coefficients for alkaline-earth atoms with alkaline-earth ions and the $C_9$ coefficients, we have performed separate fitting calculations and compared. Our calculations match in an excellent manner with the fitting calculations. We have also reported the oscillator strengths for the leading transitions and static dipole polarizabilities for the ground states of the alkaline-earth ions, i.e., Mg$^+$, Ca$^+$, Sr$^+$, and Ba$^+$ as well as the alkaline-earth atoms, i.e., Mg, Ca, Sr, and Ba. These, when compared with the available experimental results, show good agreement.


title: Merging of unequal mass binary black holes in non-axisymmetric galactic nuclei

abstract: In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $\varphi-$GPU $N$-body code, which enabled us to fully exploit supercomputers equipped with graphic processing units (GPUs). As a model for the galactic nuclei, we adopted initial axisymmetric, rotating models, aimed at reproducing the properties of a galactic nucleus emerging from a galaxy merger event, containing two SMBHs which were unbound initially. We found no 'final-parsec problem', as our SMBHs tend to pair and shrink without showing significant signs of stalling. This confirms earlier results and extends them to large particle numbers and rotating systems. We find that the SMBHB hardening depends on the binary-reduced mass ratio via a single parameter function. Our results suggest that, at a fixed value for the SMBHB primary mass, the merger time of highly asymmetric binaries is up to four order of magnitudes smaller than the equal-mass binaries. This can significantly affect the population of SMBHs potentially detectable as gravitational wave sources.


title: Reduced equations for an active model of the hydroelastic waves in the cochlea

abstract: Building upon our earlier passive models for the cochlea, here we enhance the model with an active mechanism. Starting with a one-chamber simplification leading to a system of a time-dependent PDE in two spatial variables for the pressure coupled to a PDE in one spatial variable for the oscillation of the basilar membrane, we rigorously establish the validity of a dimension reduction to a system to two ODE's. We then present numerical simulations demonstrating the ability of this reduced active system to distinguish and amplify multi-frequency input signals.


title: van der Waals coefficients of the multi-layered MoS$_2$ with alkali metals

abstract: The van der Waals coefficients and the separation dependent retardation functions of the interactions between the atomically thin films of the multi-layered transition metal molybdenum disulfide (MoS$_2$) dichalcogenides with the alkali atoms are investigated. First, we determine the frequency-dependent dielectric permittivity and intrinsic carrier density values for different layers of MoS$_2$ by adopting various fitting models to the recently measured optical data reported by Yu and co-workers [Sci. Rep. {\bf 5}, 16996 (2015)] using spectroscopy ellipsometry. Then, dynamic electric dipole polarizabilities of the alkali atoms are evaluated very accurately by employing the relativistic coupled-cluster theory. We also demonstrate the explicit change in the above coefficients for different number of layers. These studies are highly useful for the optoelectronics, sensing and storage applications using layered MoS$_2$.


title: An Assessment of the Radio Frequency Electromagnetic Field Exposure from A Massive MIMO 5G Testbed

abstract: Current radiofrequency electromagnetic field (RF-EMF) exposure limits have become a critical concern for fifth-generation (5G) mobile network deployment. Regulation is not harmonized and in certain countries and regions it goes beyond the guidelines set out by the International Commission on Non-Ionizing Radiation Protection (ICNIRP). Using a massive multiple-input-multiple-output (mMIMO) testbed with beamforming capabilities that is capable of mimicking realistic 5G base station (BS) performance, this paper presents an experimental and statistical assessment of its associated RF-EMF exposure within a real-world indoor environment. The mMIMO testbed has up to 128 channels with user-programmable software defined radio (SDR) capability. It could perform zero-forcing precoding after channel state information (CSI) acquisition for different beamforming scenarios with respect to the associated user terminal antenna setups and positions. With 64 active mMIMO transmit antennas, 8 beamforming scenarios have been considered for single-user (SU) and multi-user (MU) downlink communications at different locations. Using a calibrated triaxial isotropic field-probe, the received channel power heat map for each beamforming scenario was acquired and then converted into an RF-EMF heat map. The relevant RF-EMF statistics was evaluated based on the variations of beam profiles and number of users.


title: Optically induced topological superconductivity via Floquet interaction engineering

abstract: We propose a mechanism for light-induced unconventional superconductivity in a two-valley semiconductor with a massive Dirac type band structure. The superconducting phase results from the out-of-equilibrium excitation of carriers in the presence of Coulomb repulsion and is stabilized by coupling the driven semiconductor to a bosonic or fermionic thermal bath. We consider a circularly-polarized light pump and show that by controlling the detuning of the pump frequency relative to the band gap, different types of chiral superconductivity would be induced. The emergence of novel superconducting states, such as the chiral $p$-wave pairing, results from the Floquet engineering of the interaction. This is realized by modifying the form of the Coulomb interaction by projecting it into the states that are resonant with the pump frequency. We show that the resulting unconventional pairing in our system can host topologically protected chiral bound states. We discuss a promising experimental platform to realize our proposal and detect the signatures of the emergent superconducting state.


title: COVID-19 in South Asia: Real-time monitoring of reproduction and case fatality rate

abstract: As the ravages caused by COVID-19 pandemic are becoming inevitable with every moment, monitoring and understanding of transmission and fatality rate has become even more paramount for containing its spread. The key purpose of this analysis is to report the real-time effective reproduction rate ($R_t$ ) and case fatality rates (CFR) of COVID-19 in South Asia region. Data for this study are extracted from JHU CSSE COVID-19 Data source up to July 31, 2020. $R_t$ is estimated using exponential growth and time-dependent methods. R0 package in R-language is employed to estimate $R_t$ by fitting the existing epidemic curve. Case fatality rate is estimated by using Naive and Kaplan-Meier methods. Owing to exponential increase in cases of COVID-19, the pandemic will ensue in India, Maldives and in Nepal as $R_t$ was estimated greater than 1 for these countries. Although case fatality rates are found lesser as compared to other highly affected regions in the world, strict monitoring of deaths for better health facilities and care of patients is emphasized. More regional level cooperation and efforts are the need of time to minimize the detrimental effects of the virus.


title: Design based incomplete U-statistics

abstract: U-statistics are widely used in fields such as economics, machine learning, and statistics. However, while they enjoy desirable statistical properties, they have an obvious drawback in that the computation becomes impractical as the data size $n$ increases. Specifically, the number of combinations, say $m$, that a U-statistic of order $d$ has to evaluate is $O(n^d)$. Many efforts have been made to approximate the original U-statistic using a small subset of combinations since Blom (1976), who referred to such an approximation as an incomplete U-statistic. To the best of our knowledge, all existing methods require $m$ to grow at least faster than $n$, albeit more slowly than $n^d$, in order for the corresponding incomplete U-statistic to be asymptotically efficient in terms of the mean squared error. In this paper, we introduce a new type of incomplete U-statistic that can be asymptotically efficient, even when $m$ grows more slowly than $n$. In some cases, $m$ is only required to grow faster than $\sqrt{n}$. Our theoretical and empirical results both show significant improvements in the statistical efficiency of the new incomplete U-statistic.


title: Transverse spin sum rule of the proton

abstract: The transversely-polarized state of a proton with arbitrary momentum is not an eigenstate of transverse angular momentum operator. The latter does not commute with the QCD Hamiltonian. However, the expectation value of the transverse angular momentum in the state is well-defined and grows proportionally to the energy of the particle. The transverse spin content of the proton is analyzed in terms of the QCD angular momentum structure. In particular, we reconfirm that the generalized parton distributions $H+E$ provide transverse orbital angular momentum densities of quarks and gluons in the infinite momentum frame.


title: Achieving large zeta-potentials with charged porous surfaces

abstract: We discuss an electro-osmotic flow near charged porous coatings of a finite hydrodynamic permeability, impregnated with an outer electrolyte solution. It is shown that their electrokinetic (zeta) potential is generally augmented compared to the surface electrostatic potential, thanks to a large liquid slip at their surface emerging due to an electro-osmotic flow in the enriched by counter-ions porous films. The inner flow shows a very rich behavior controlled by the volume charge density of the coating, its Brinkman length, and concentration of added salt. Interestingly, even for relatively small Brinkman length the zeta-potential can, in some cases, become huge, providing a very fast outer flow in the bulk electrolyte. When the Brinkman length is large enough, the zeta-potential could be extremely high even at practically vanishing surface potential. To describe the slip velocity in a simple manner, we introduce a concept of an electro-osmotic slip length and demonstrate that the latter is always defined by the hydrodynamic permeability of the porous film, and also, depending on the regime, either by its volume charge density or the salt concentration. These results provide a framework for the rational design of porous coatings to enhance electrokineic phenomena, and for tuning their properties by adjusting bulk electrolyte concentrations, with direct applications in microfluidics.


title: Leveraging Vehicle Connectivity and Autonomy to Stabilize Flow in Mixed Traffic Conditions: Accounting for Human-driven Vehicle Driver Behavioral Heterogeneity and Perception-reaction Time Delay

abstract: The erratic nature of human driving tends to trigger undesired waves that amplify as successive driver reactions propagate from the errant vehicle to vehicles upstream. Known as phantom jams, this phenomenon has been identified in the literature as one of the main causes of traffic congestion. This paper is based on the premise that vehicle automation and connectivity can help mitigate such jams. In the paper, we design a controller for use in a connected and autonomous vehicle (CAV) to stabilize the flow of human-driven vehicles (HDVs) that are upstream of the CAV, and consequently to lower collision risk in the upstream traffic environment. In modeling the HDV dynamics in the mixed traffic stream, we duly consider HDV driver heterogeneity and the time delays associated with their perception reaction time. We can find that the maximum number of HDVs that a CAV can stabilize is lower when human drivers potential time delay and heterogeneity are considered, compared to the scenario where such are not considered. This result suggests that heterogeneity and time delay in HDV behavior impairs the CAVs capability to stabilize traffic. Therefore, in designing CAV controllers for traffic stabilization, it is essential to consider such uncertainty-related conditions. In our demonstration, we also show that the designed controller can significantly improve both the stability of the mixed traffic stream and the safety of both CAVs and HDVs in the stream. The results are useful for real-time calibration of the model parameters that characterize HDV movements in the mixed stream.


title: Secure Simultaneous Information and Power Transfer for Downlink Multi-user Massive MIMO

abstract: In this paper, downlink secure transmission in simultaneous information and power transfer (SWIPT) system enabled with massive multiple-input multiple-output (MIMO) is studied. A base station (BS) with a large number of antennas transmits energy and information signals to its intended users, but these signals are also received by an active eavesdropper. The users and eavesdropper employ a power splitting technique to simultaneously decode information and harvest energy. Massive MIMO helps the BS to focus energy to the users and prevent information leakage to the eavesdropper. The harvested energy by each user is employed for decoding information and transmitting uplink pilot signals for channel estimation. It is assumed that the active eavesdropper also harvests energy in the downlink and then contributes during the uplink training phase. Achievable secrecy rate is considered as the performance criterion and a closed-form lower bound for it is derived. To provide secure transmission, the achievable secrecy rate is then maximized through an optimization problem with constraints on the minimum harvested energy by the user and the maximum harvested energy by the eavesdropper. Numerical results show the effectiveness of using massive MIMO in providing physical layer security in SWIPT systems and also show that our closed-form expressions for the secrecy rate are accurate.


title: Co-design and Co-simulation for Engineering Systems: Insights from the Sustainable Infrastructure Planning Game

abstract: This paper draws on perspectives from co-design as an integrative and collaborative design activity and co-simulation as a supporting information system to advance engineering design methods for problems of societal significance. Design and implementation of the Sustainable Infrastructure Planning Game provides a prototypical co-design artifact that leverages the High Level Architecture co-simulation standard. Three role players create a strategic infrastructure plan for agriculture, water, and energy sectors to meet sustainability objectives for a growing and urbaninzing population in a fictional desert nation. An observational study conducts 15 co-design sessions to understand underlying dynamics between actors and how co-simulation capabilities influence design outcomes. Results characterize the dependencies and conflicts between player roles based on technical exchange of resource flows, identifying tension between agriculture and water roles based on water demands for irrigation. Analysis shows a correlation between data exchange, facilitated by synchronous co-simulation, and highly-ranked achievement of joint sustainability outcomes. Conclusions reflect on the opportunities and challenges presented by co-simulation in co-design settings to address engineering systems problems.


title: A characterization of the uniform convergence points set of some convergent sequence of functions

abstract: We characterize the uniform convergence points set of a pointwisely convergent sequence of real-valued functions defined on a perfectly normal space. We prove that if $X$ is a perfectly normal space which can be covered by a disjoint sequence of dense subsets and $A\subseteq X$, then $A$ is the set of points of the uniform convergence for some convergent sequence $(f_n)_{n\in\omega}$ of functions $f_n:X\to \mathbb R$ if and only if $A$ is $G_\delta$-set which contains all isolated points of $X$. This result generalizes a theorem of J\'{a}n Bors\'{i}k published in 2019.


title: A Two-Stage Metaheuristic Algorithm for the Dynamic Vehicle Routing Problem in Industry 4.0 approach

abstract: Industry 4.0 is a concept that assists companies in developing a modern supply chain (MSC) system when they are faced with a dynamic process. Because Industry 4.0 focuses on mobility and real-time integration, it is a good framework for a dynamic vehicle routing problem (DVRP). This research works on DVRP. The aim of this research is to minimize transportation cost without exceeding the capacity constraint of each vehicle while serving customer demands from a common depot. Meanwhile, new orders arrive at a specific time into the system while the vehicles are executing the delivery of existing orders. This paper presents a two-stage hybrid algorithm for solving the DVRP. In the first stage, construction algorithms are applied to develop the initial route. In the second stage, improvement algorithms are applied. Experimental results were designed for different sizes of problems. Analysis results show the effectiveness of the proposed algorithm.


title: An Efficient Sliding Mesh Interface Method for High-Order Discontinuous Galerkin Schemes

abstract: Sliding meshes are a powerful method to treat deformed domains in computational fluid dynamics, where different parts of the domain are in relative motion. In this paper, we present an efficient implementation of a sliding mesh method into a discontinuous Galerkin compressible Navier-Stokes solver and its application to a large eddy simulation of a 1-1/2 stage turbine. The method is based on the mortar method and is high-order accurate. It can handle three-dimensional sliding mesh interfaces with various interface shapes. For plane interfaces, which are the most common case, conservativity and free-stream preservation are ensured. We put an emphasis on efficient parallel implementation. Our implementation generates little computational and storage overhead. Inter-node communication via MPI in a dynamically changing mesh topology is reduced to a bare minimum by ensuring a priori information about communication partners and data sorting. We provide performance and scaling results showing the capability of the implementation strategy. Apart from analytical validation computations and convergence results, we present a wall-resolved implicit LES of the 1-1/2 stage Aachen turbine test case as a large scale practical application example.


title: Directional Laplacian Centrality for Cyber Situational Awareness

abstract: Cyber operations is drowning in diverse, high-volume, multi-source data. In order to get a full picture of current operations and identify malicious events and actors analysts must see through data generated by a mix of human activity and benign automated processes. Although many monitoring and alert systems exist, they typically use signature-based detection methods. We introduce a general method rooted in spectral graph theory to discover patterns and anomalies without a priori knowledge of signatures. We derive and propose a new graph-theoretic centrality measure based on the derivative of the graph Laplacian matrix in the direction of a vertex. To build intuition about our measure we show how it identifies the most central vertices in standard network data sets and compare to other graph centrality measures. Finally, we focus our attention on studying its effectiveness in identifying important IP addresses in network flow data. Using both real and synthetic network flow data, we conduct several experiments to test our measure's sensitivity to two types of injected attack profiles, and show that vertices participating in injected attack profiles exhibit noticeable changes in our centrality measures, even when the injected anomalies are relatively small, and in the presence of simulated network dynamics.


title: On the Characterization of Generalized Derivatives for the Solution Operator of the Bilateral Obstacle Problem

abstract: We consider optimal control problems for a wide class of bilateral obstacle problems where the control appears in a possibly nonlinear source term. The non-differentiability of the solution operator poses the main challenge for the application of efficient optimization methods and the characterization of Bouligand generalized derivatives of the solution operator is essential for their theoretical foundation and numerical realization. In this paper, we derive specific elements of the Bouligand generalized differential if the control operator satisfies natural monotonicity properties. We construct monotone sequences of controls where the solution operator is G\^ateaux differentiable and characterize the corresponding limit element of the Bouligand generalized differential as being the solution operator of a Dirichlet problem on a quasi-open domain. In contrast to a similar recent result for the unilateral obstacle problem [RU19], we have to deal with an opposite monotonic behavior of the active and strictly active sets corresponding to the upper and lower obstacle. Moreover, the residual is no longer a nonnegative functional on H^{-1} and its representation as the difference of two nonnegative Radon measures requires special care. This necessitates new proof techniques that yield two elements of the Bouligand generalized differential. Also for the unilateral case we obtain an additional element to that derived in [RU19].


title: Quantum memory enhanced dissipative entanglement creation in non-equilibrium steady states

abstract: This Article investigates dissipative preparation of entangled non-equilibrium steady states (NESS). We construct a collision model where the open system consists of two qubits which are coupled to heat reservoirs with different temperatures. The baths are modeled by sequences of qubits interacting with the open system. The model can be studied in different dynamical regimes: with and without environmental memory effects. We report that only a certain bath temperature range allows for entangled NESS. Furthermore, we obtain minimal and maximal critical values for the heat current through the system. Surprisingly, quantum memory effects play a crucial role in the long time limit. First, memory effects broaden the parameter region where quantum correlated NESS may be dissipatively prepared and, secondly, they increase the attainable concurrence. Most remarkably, we find a heat current range that does not only allow but guarantees that the NESS is entangled. Thus, the heat current can witness entanglement of non-equilibrium steady states.


title: A two-carrier scheme: evading the 3dB quantum penalty of heterodyne readout in gravitational-wave detectors

abstract: Precision measurements using traditional heterodyne readout suffer a 3dB quantum noise penalty compared with homodyne readout. The extra noise is caused by the quantum fluctuations in the image vacuum. We propose a two-carrier gravitational-wave detector design that evades the 3dB quantum penalty of heterodyne readout. We further propose a new way of realising frequency-dependent squeezing utilising two-mode squeezing in our scheme. It naturally achieves more precise audio frequency signal measurements with radio frequency squeezing. In addition, the detector is compatible with other quantum nondemolition techniques.


title: Gentle introduction to rigorous Renormalization Group: a worked fermionic example

abstract: Much of our understanding of critical phenomena is based on the notion of Renormalization Group (RG), but the actual determination of its fixed points is usually based on approximations and truncations, and predictions of physical quantities are often of limited accuracy. The RG fixed points can be however given a fully rigorous and non-perturbative characterization, and this is what is presented here in a model of symplectic fermions with a nonlocal ("long-range") kinetic term depending on a parameter $\varepsilon$ and a quartic interaction. We identify the Banach space of interactions, which the fixed point belongs to, and we determine it via a convergent approximation scheme. The Banach space is not limited to relevant interactions, but it contains all possible irrelevant terms with short-ranged kernels, decaying like a stretched exponential at large distances. As the model shares a number of features in common with $\phi^4$ or Ising models, the result can be used as a benchmark to test the validity of truncations and approximations in RG studies. The analysis is based on results coming from Constructive RG to which we provide a tutorial and self-contained introduction. In addition, we prove that the fixed point is analytic in $\varepsilon$, a somewhat surprising fact relying on the fermionic nature of the problem.

categories: hep-th cond-mat.stat-mech math-ph math.MP math.PR


title: Coherence measures induced by norm functions

abstract: Which matrix norms induce proper measures for quantifying quantum coherence? We study this problem for two important classes of norms and show that (i) coherence measures cannot be induced by any unitary similarity invariant norm, and (ii) the $\ell_{q,p}$-norm induces a coherence measure if and only if $q=1$ and $1 \leq p \leq 2$, thus giving a new class of coherence measures with simple closed forms that are easy to compute. These results extend and unify previously known facts about norm-induced coherence measures, and lead to a broader framework for understanding what functionals can be coherence measures.


title: Distributed Personalized Gradient Tracking with Convex Parametric Models

abstract: We present a distributed optimization algorithm for solving online personalized optimization problems over a network of computing and communicating nodes, each of which linked to a specific user. The local objective functions are assumed to have a composite structure and to consist of a known time-varying (engineering) part and an unknown (user-specific) part. Regarding the unknown part, it is assumed to have a known parametric (e.g., quadratic) structure a priori, whose parameters are to be learned along with the evolution of the algorithm. The algorithm is composed of two intertwined components: (i) a dynamic gradient tracking scheme for finding local solution estimates and (ii) a recursive least squares scheme for estimating the unknown parameters via user's noisy feedback on the local solution estimates. The algorithm is shown to exhibit a bounded regret under suitable assumptions. Finally, a numerical example corroborates the theoretical analysis.


title: Entanglement and Impropriety

abstract: The relationship between quantum entanglement and classical impropriety is considered in the context of multi-modal squeezed states of light. Replacing operators with complex Gaussian random variables in the Bogoliubov transformations for squeezed states, we find that the resulting transformed variables are not only correlated but also improper. A simple threshold exceedance model of photon detection is considered and used to demonstrate how the behavior of improper Gaussian random variables can mimic that of entangled photon pairs when coincidence post-selection is performed.


title: Horizon Hair from Inversion Symmetry

abstract: The extreme Reissner-Nordstr\"om solution has a discrete conformal isometry that maps the future event horizon to future null infinity and vice versa, the Couch-Torrence inversion isometry. We study the dynamics of a probe Maxwell field on the extreme Reissner-Nordstr\"om solution in light of this symmetry. We present a gauge fixing that is compatible with the inversion symmetry. The gauge fixing allows us to relate the gauge parameter at the future horizon to future null infinity, which further allows us to study global charges for large gauge symmetries in the exterior of the extreme Reissner-Nordstr\"om black hole. Along the way, we construct Newman-Penrose and Aretakis like conserved quantities along future null infinity and the future event horizon, respectively, and relate them via the Couch-Torrence inversion symmetry.


title: Large-Scale Lipschitz Estimates for Elliptic Systems with Periodic High-Contrast Coefficients

abstract: This paper is concerned with the large-scale regularity in the homogenization of elliptic systems of elasticity with periodic high-contrast coefficients. We obtain the large-scale Lipschitz estimate that is uniform with respect to the contrast ratio $\delta^2$ for $0<\delta<\infty$. Our study also covers the case of soft inclusions ($\delta=0$) as well as the case of stiff inclusions ($\delta=\infty$). The large-scale Lipschitz estimate, together with classical local estimates, allows us to establish explicit bounds for the matrix of fundamental solutions and its derivatives.


title: Deep Detail Enhancement for Any Garment

abstract: Creating fine garment details requires significant efforts and huge computational resources. In contrast, a coarse shape may be easy to acquire in many scenarios (e.g., via low-resolution physically-based simulation, linear blend skinning driven by skeletal motion, portable scanners). In this paper, we show how to enhance, in a data-driven manner, rich yet plausible details starting from a coarse garment geometry. Once the parameterization of the garment is given, we formulate the task as a style transfer problem over the space of associated normal maps. In order to facilitate generalization across garment types and character motions, we introduce a patch-based formulation, that produces high-resolution details by matching a Gram matrix based style loss, to hallucinate geometric details (i.e., wrinkle density and shape). We extensively evaluate our method on a variety of production scenarios and show that our method is simple, light-weight, efficient, and generalizes across underlying garment types, sewing patterns, and body motion.


title: Multiple metastable states in an off-lattice Potts model

abstract: The interactions between a group of components are commonly studied in several areas of science (social science, biology, material science, complex dynamical systems, among others) using the methods of thermodynamics and statistical mechanics. In this work we study the properties of the recently proposed off-lattice, two-dimensional Potts model [Eur. Phys. J. B 87, 78 (2014)], originally motivated by the dynamics of agent opinions, and which is described by a Hamiltonian obtained by a maximum entropy inference procedure. We performed microcanonical and canonical Monte Carlo simulations of the first-order phase transition in the model, revealing a caloric curve with metastable regions. Furthermore, we report a "switching" behavior between multiple metastable states. We also note that the thermodynamics of the model has striking similarities with systems having long-range interactions, even though the interactions are short-ranged.


title: Cross-correlation between Subaru Hyper Suprime-Cam Galaxy Weak Lensing and Planck Cosmic Microwave Background Lensing

abstract: Cross-correlations between galaxy weak lensing (WL) and Cosmic Microwave Background (CMB) lensing are a powerful tool to probe matter fluctuations at intermediate redshifts and to detect residual systematics in either probe. In this paper, we study the cross-correlation of galaxy WL from the Hyper Suprime-Cam Subaru Strategic Program (HSC) first data release and CMB lensing from the final Planck data release, for HSC source galaxies at 0.3< z < 1.5. HSC is the deepest Stage-III galaxy WL survey, and provides both a great opportunity to study the high-redshift universe and new challenges related to its exceptionally high source density, such as source blending. The cross-correlation signal is measured at a significance level of 3.1$\sigma$. The amplitude of our best-fit model with respect to the best-fit 2018 Planck cosmology is $A = 0.81\pm 0.25$, consistent with $A=1$. Our result is also consistent with previous CMB lensing and galaxy WL cross-correlation studies using different surveys. We perform tests with respect to the WL $B$-modes, the point-spread-function, photometric redshift errors, and thermal Sunyaev-Zel'dovich leakage, and find no significant evidence of residual systematics.


title: Predicting Risk of Developing Diabetic Retinopathy using Deep Learning

abstract: Diabetic retinopathy (DR) screening is instrumental in preventing blindness, but faces a scaling challenge as the number of diabetic patients rises. Risk stratification for the development of DR may help optimize screening intervals to reduce costs while improving vision-related outcomes. We created and validated two versions of a deep learning system (DLS) to predict the development of mild-or-worse ("Mild+") DR in diabetic patients undergoing DR screening. The two versions used either three-fields or a single field of color fundus photographs (CFPs) as input. The training set was derived from 575,431 eyes, of which 28,899 had known 2-year outcome, and the remaining were used to augment the training process via multi-task learning. Validation was performed on both an internal validation set (set A; 7,976 eyes; 3,678 with known outcome) and an external validation set (set B; 4,762 eyes; 2,345 with known outcome). For predicting 2-year development of DR, the 3-field DLS had an area under the receiver operating characteristic curve (AUC) of 0.79 (95%CI, 0.78-0.81) on validation set A. On validation set B (which contained only a single field), the 1-field DLS's AUC was 0.70 (95%CI, 0.67-0.74). The DLS was prognostic even after adjusting for available risk factors (p<0.001). When added to the risk factors, the 3-field DLS improved the AUC from 0.72 (95%CI, 0.68-0.76) to 0.81 (95%CI, 0.77-0.84) in validation set A, and the 1-field DLS improved the AUC from 0.62 (95%CI, 0.58-0.66) to 0.71 (95%CI, 0.68-0.75) in validation set B. The DLSs in this study identified prognostic information for DR development from CFPs. This information is independent of and more informative than the available risk factors.


title: Investigation of temperature dependence of magnetic properties of Cr$_2$O$_3$ thin film structure using a magnetic field imaging technique based on Nitrogen-Vacancy centres in diamond crystal

abstract: This work presents a magnetic field imaging method based on color centres in diamond crystal applied to thin film structure. To demonstrate the capacity of our device we have used it for characterization of magnetic properties in microscopic scale of Cr$_2$O$_3$ thin film structure above and below N\'eel temperature. The obtained measurement results clearly identify the detection of the magnetic phase transition of Cr$_2$O$_3$ thin film with an unexpected diamagnetic like behaviour at 19$^{\circ}$C (below the N\'eel temperature of Cr$_2$O$_3$). To have better insights in the magnetic fields created by the thin films we present simulations of the magnetic fields near the thin film surface. We also analysed the optically detected magnetic resonance (ODMR) profiles to be sure that the measured property is related only to the shift of the optically detected magnetic resonance. We demonstrate how Nitrogen-Vacancy centre based magnetometry can deliver a convenient platform for research of phase transition dynamics and spatial magnetic field distributions of thin films, that can rival other widely used measurement techniques.


title: Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small $x$

abstract: We compute the differential cross-section for inclusive prompt photon$+$quark production in deeply inelastic scattering of electrons off nuclei at small $x$ ($e+A$ DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling $\alpha_{\mathrm{s}}$) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale $Q_{s,A}(x)$. We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive $e+A$ DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small $x$ in EIC kinematics require jet reconstruction at small $k_\perp$, it will be important to study their feasibility relative to photon-hadron correlations.


title: Exploring Navigation Styles in a FutureLearn MOOC

abstract: This paper presents for the first time a detailed analysis of fine-grained navigation style identification in MOOCs backed by a large number of active learners. The result shows 1) whilst the sequential style is clearly in evidence, the global style is less prominent; 2) the majority of the learners do not belong to either category; 3) navigation styles are not as stable as believed in the literature; and 4) learners can, and do, swap between navigation styles with detrimental effects. The approach is promising, as it provides insight into online learners' temporal engagement, as well as a tool to identify vulnerable learners, which potentially benefit personalised interventions (from teachers or automatic help) in Intelligent Tutoring Systems (ITS).


title: Can We Spot the "Fake News" Before It Was Even Written?

abstract: Given the recent proliferation of disinformation online, there has been also growing research interest in automatically debunking rumors, false claims, and "fake news." A number of fact-checking initiatives have been launched so far, both manual and automatic, but the whole enterprise remains in a state of crisis: by the time a claim is finally fact-checked, it could have reached millions of users, and the harm caused could hardly be undone. An arguably more promising direction is to focus on fact-checking entire news outlets, which can be done in advance. Then, we could fact-check the news before it was even written: by checking how trustworthy the outlets that published it is. We describe how we do this in the Tanbih news aggregator, which makes readers aware of what they are reading. In particular, we develop media profiles that show the general factuality of reporting, the degree of propagandistic content, hyper-partisanship, leading political ideology, general frame of reporting, and stance with respect to various claims and topics.


title: Triplet-blockaded Josephson supercurrent in double quantum dots

abstract: Serial double quantum dots created in semiconductor nanostructures provide a versatile platform for investigating two-electron spin quantum states, which can be tuned by electrostatic gating and an external magnetic field. In this work, we directly measure the supercurrent reversal between adjacent charge states of an InAs nanowire double quantum dot with superconducting leads, in good agreement with theoretical models. In the even charge parity sector, we observe a supercurrent blockade with increasing magnetic field, corresponding to the spin singlet to triplet transition. Our results demonstrate a direct spin to supercurrent conversion, the superconducting equivalent of the Pauli spin blockade. This effect can be exploited in hybrid quantum architectures coupling the quantum states of spin systems and superconducting circuits.


title: Classical Sub-subleading Soft Photon and Soft Graviton Theorems in Four Spacetime Dimensions

abstract: Classical soft photon and soft graviton theorems determine long wavelength electromagnetic and gravitational waveforms for a general classical scattering process in terms of the electric charges and asymptotic momenta of the ingoing and outgoing macroscopic objects. Performing Fourier transformation of the electromagnetic and gravitational waveforms in the frequency variable one finds electromagnetic and gravitational waveforms at late and early retarded time. Here extending the formalism developed in \cite{1912.06413}, we derive sub-subleading electromagnetic and gravitational waveforms which behave like $u^{-2}(\ln u)$ at early and late retarded time $u$ in four spacetime dimensions. We also have derived the sub-subleading soft photon theorem analyzing two loop amplitudes in scalar QED. Finally, we give the structure of leading non-analytic contribution to (sub)$^{n}$-leading classical soft photon and graviton theorems which behave like $u^{-n}(\ln u)^{n-1}$ for early and late retarded time $u$.


title: An Automated, End-to-End Framework for Modeling Attacks From Vulnerability Descriptions

abstract: Attack graphs are one of the main techniques used to automate the risk assessment process. In order to derive a relevant attack graph, up-to-date information on known attack techniques should be represented as interaction rules. Designing and creating new interaction rules is not a trivial task and currently performed manually by security experts. However, since the number of new security vulnerabilities and attack techniques continuously and rapidly grows, there is a need to frequently update the rule set of attack graph tools with new attack techniques to ensure that the set of interaction rules is always up-to-date. We present a novel, end-to-end, automated framework for modeling new attack techniques from textual description of a security vulnerability. Given a description of a security vulnerability, the proposed framework first extracts the relevant attack entities required to model the attack, completes missing information on the vulnerability, and derives a new interaction rule that models the attack; this new rule is integrated within MulVAL attack graph tool. The proposed framework implements a novel pipeline that includes a dedicated cybersecurity linguistic model trained on the the NVD repository, a recurrent neural network model used for attack entity extraction, a logistic regression model used for completing the missing information, and a novel machine learning-based approach for automatically modeling the attacks as MulVAL's interaction rule. We evaluated the performance of each of the individual algorithms, as well as the complete framework and demonstrated its effectiveness.


title: Unsupervised Deep Metric Learning with Transformed Attention Consistency and Contrastive Clustering Loss

abstract: Existing approaches for unsupervised metric learning focus on exploring self-supervision information within the input image itself. We observe that, when analyzing images, human eyes often compare images against each other instead of examining images individually. In addition, they often pay attention to certain keypoints, image regions, or objects which are discriminative between image classes but highly consistent within classes. Even if the image is being transformed, the attention pattern will be consistent. Motivated by this observation, we develop a new approach to unsupervised deep metric learning where the network is learned based on self-supervision information across images instead of within one single image. To characterize the consistent pattern of human attention during image comparisons, we introduce the idea of transformed attention consistency. It assumes that visually similar images, even undergoing different image transforms, should share the same consistent visual attention map. This consistency leads to a pairwise self-supervision loss, allowing us to learn a Siamese deep neural network to encode and compare images against their transformed or matched pairs. To further enhance the inter-class discriminative power of the feature generated by this network, we adapt the concept of triplet loss from supervised metric learning to our unsupervised case and introduce the contrastive clustering loss. Our extensive experimental results on benchmark datasets demonstrate that our proposed method outperforms current state-of-the-art methods for unsupervised metric learning by a large margin.


title: A Survey and Insights on Deployments of the Connected and Autonomous Vehicles in US

abstract: CV/ITS (Connected Vehicle, Intelligent Transportation System) and AV/ADS (Autonomous Vehicle, Automated Driving System) have been emerging for the sake of saving people lives, improving traffic efficiency and helping the environment for decades. There are separate efforts led respectively by USDOT with state DOTs for CV, and private sectors through market driven approach from start-ups and technology companies for AV. By CV/ITS effort there are 97 deployments of V2X communications utilizing the 5.9 GHz band, 18,877 vehicles with aftermarket V2X communications devices, and 8,098 infrastructure V2X devices installed at the roadsides. However, CV/ITS still cannot be massively deployed in US markets due to lack of regulations, dedicated wireless spectrum bands, sustainable financial & business models with mature supply chain, etc. In the other side, technology-driven AV market has been much slower than expected mainly because of immaturity of AI technology to handle different complex driving scenarios in a cost effective way. In this paper, we first present these two parallel journeys focusing on the deployments including operating models, scenarios and applications, evaluations and lessons learning. Then, come up with recommendations to a cooperative CAV approach driving a more feasible, safer, affordable and cost effective transportation, but require a great industry collaboration from Automotive, Transportation. ICT and Cloud.


title: When do solar erupting hot magnetic flux ropes form?

abstract: We investigate the formation times of eruptive magnetic flux ropes relative to the onset of solar eruptions, which is important for constraining models of coronal mass ejection (CME) initiation. We inspected uninterrupted sequences of 131 \AA\ images that spanned more than eight hours and were obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) to identify the formation times of hot flux ropes that erupted in CMEs from locations close to the limb. The appearance of the flux ropes as well as their evolution toward eruptions were determined using morphological criteria. Two-thirds (20/30) of the flux ropes were formed well before the onset of the eruption (from 51 minutes to more than eight hours), and their formation was associated with the occurrence of a confined flare. We also found four events with preexisting hot flux ropes whose formations occurred a matter of minutes (from three to 39) prior to the eruptions without any association with distinct confined flare activity. Six flux ropes were formed once the eruptions were underway. However, in three of them, prominence material could be seen in 131 \AA\ images, which may indicate the presence of preexisting flux ropes that were not hot. The formation patterns of the last three groups of hot flux ropes did not show significant differences. For the whole population of events, the mean and median values of the time difference between the onset of the eruptive flare and the appearance of the hot flux rope were 151 and 98 minutes, respectively. Our results provide, on average, indirect support for CME models that involve preexisting flux ropes; on the other hand, for a third of the events, models in which the ejected flux rope is formed during the eruption appear more appropriate.


title: Bipartite Graph Reasoning GANs for Person Image Generation

abstract: We present a novel Bipartite Graph Reasoning GAN (BiGraphGAN) for the challenging person image generation task. The proposed graph generator mainly consists of two novel blocks that aim to model the pose-to-pose and pose-to-image relations, respectively. Specifically, the proposed Bipartite Graph Reasoning (BGR) block aims to reason the crossing long-range relations between the source pose and the target pose in a bipartite graph, which mitigates some challenges caused by pose deformation. Moreover, we propose a new Interaction-and-Aggregation (IA) block to effectively update and enhance the feature representation capability of both person's shape and appearance in an interactive way. Experiments on two challenging and public datasets, i.e., Market-1501 and DeepFashion, show the effectiveness of the proposed BiGraphGAN in terms of objective quantitative scores and subjective visual realness. The source code and trained models are available at https://github.com/Ha0Tang/BiGraphGAN.


title: Uncertainty Quantification of Structural Systems with Subset of Data

abstract: Quantification of the impact of uncertainty in material properties as well as the input ground motion on structural responses is an important step in implementing a performance-based earthquake engineering (PBEE) framework. Among various sources of uncertainty, the variability in the input ground motions, a.k.a. record-to-record, greatly affects the assessment results. The objective of this paper is to quantify the uncertainty in structural response with hybrid uncertainty sources. In this paper, multiple matrix completion methods are proposed and applied on a case study structure. The matrix completion method is a means to estimate the analyses results for the entire set of input parameters by conducting analysis for only a small subset of analyses. The main algorithmic contributions of our proposed method are twofold. First, we develop a sampling technique for choosing a subset of representative simulations, which allows improving the accuracy of the estimated response. An unsupervised machine learning technique is used for this purpose. Next, the proposed matrix completion method for uncertainty quantification is further refined by incorporating a regression model that is trained on the available partial simulations. The regression model improves the initial sampling as it provides a rough estimation of the structural responses. Finally, the proposed algorithm is applied to a multi-degree-of-freedom system, and the structural responses (i.e., displacements and base shear) are estimated. Results show that the proposed algorithm can effectively estimate the response from a full set of nonlinear simulations by conducting analyses only on a small portion of the set.


title: Influence Spread in the Heterogeneous Multiplex Linear Threshold Model

abstract: The linear threshold model (LTM) has been used to study spread on single-layer networks defined by one inter-agent sensing modality and agents homogeneous in protocol. We define and analyze the heterogeneous multiplex LTM to study spread on multi-layer networks with each layer representing a different sensing modality and agents heterogeneous in protocol. Protocols are designed to distinguish signals from different layers: an agent becomes active if a sufficient number of its neighbors in each of any $a$ of the $m$ layers is active. We focus on Protocol OR, when $a=1$, and Protocol AND, when $a=m$, which model agents that are most and least readily activated, respectively. We develop theory and algorithms to compute the size of the spread at steady state for any set of initially active agents and to analyze the role of distinguished sensing modalities, network structure, and heterogeneity. We show how heterogeneity manages the tension in spreading dynamics between sensitivity to inputs and robustness to disturbances.


title: Seeing Relativity -- III. Journeying within the Kerr metric toward the negative gravity region

abstract: In this paper we study some features of the Kerr metric both from an analytic and a visual point of view by performing accurate raytracing in various situations. We focus on features that are unique to the maximal analytic extension of the Kerr metric as compared to that of the Schwarzschild or even the Reissner-Nordstr\"om one. A large number of new, yet underexplored phenomena appear, especially regarding the structure of bounded null geodesics and the aspect of the negative gravity regions whose visual characteristics are shown both from outside and inside it.


title: Gravitational formfactors of the $ \rho $ meson in QCD sum rules

abstract: By using the quark part of the energy-momentum tensor current, the gravitational formfactors of the $ \rho $ meson are calculated within the light-cone sum rules method. In the considered version, the energy-momentum tensor current is not conserved and as a result, there appear nine formfactors, six (three) of which correspond to the conservation (nonconservation) of the energy-momentum tensor current. We also compare our results with the one existing in the literature.


title: Obnoxious facility location: the case of weighted demand points

abstract: The problem considered in this paper is the weighted obnoxious facility location in the convex hull of demand points. The objective function is to maximize the smallest weighted distance between a facility and a set of demand points. Three new optimal solution approaches are proposed. Two variants of the "Big Triangle Small Triangle" global optimization method, and a procedure based on intersection points between Apollonius circles. We also compared the results with a multi-start approach using the non-linear multi-purpose software SNOPT. Problems with 1,000 demand points are optimally solved in a fraction of a second of computer time.


title: Probability Link Models with Symmetric Information Divergence

abstract: This paper introduces link functions for transforming one probability distribution to another such that the Kullback-Leibler and R\'enyi divergences between the two distributions are symmetric. Two general classes of link models are proposed. The first model links two survival functions and is applicable to models such as the proportional odds and change point, which are used in survival analysis and reliability modeling. A prototype application involving the proportional odds model demonstrates advantages of symmetric divergence measures over asymmetric measures for assessing the efficacy of features and for model averaging purposes. The advantages include providing unique ranks for models and unique information weights for model averaging with one-half as much computation requirement of asymmetric divergences. The second model links two cumulative probability distribution functions. This model produces a generalized location model which are continuous counterparts of the binary probability models such as probit and logit models. Examples include the generalized probit and logit models which have appeared in the survival analysis literature, and a generalized Laplace model and a generalized Student-$t$ model, which are survival time models corresponding to the respective binary probability models. Lastly, extensions to symmetric divergence between survival functions and conditions for copula dependence information are presented.

categories: stat.ML cs.IT cs.LG math.IT math.ST stat.TH


title: GRIMGEP: Learning Progress for Robust Goal Sampling in Visual Deep Reinforcement Learning

abstract: Designing agents, capable of learning autonomously a wide range of skills is critical in order to increase the scope of reinforcement learning. It will both increase the diversity of learned skills and reduce the burden of manually designing reward functions for each skill. Self-supervised agents, setting their own goals, and trying to maximize the diversity of those goals have shown great promise towards this end. However, a currently known limitation of agents trying to maximize the diversity of sampled goals is that they tend to get attracted to noise or more generally to parts of the environments that cannot be controlled (distractors). When agents have access to predefined goal features or expert knowledge, absolute Learning Progress (ALP) provides a way to distinguish between regions that can be controlled and those that cannot. However, those methods often fall short when the agents are only provided with raw sensory inputs such as images. In this work we extend those concepts to unsupervised image-based goal exploration. We propose a framework that allows agents to autonomously identify and ignore noisy distracting regions while searching for novelty in the learnable regions to both improve overall performance and avoid catastrophic forgetting. Our framework can be combined with any state-of-the-art novelty seeking goal exploration approaches. We construct a rich 3D image based environment with distractors. Experiments on this environment show that agents using our framework successfully identify interesting regions of the environment, resulting in drastically improved performances. The source code is available at https://sites.google.com/view/grimgep.


title: Central Limit Theorem in High Dimensions : The Optimal Bound on Dimension Growth Rate

abstract: In this article, we try to give an answer to the simple question: ``\textit{What is the critical growth rate of the dimension $p$ as a function of the sample size $n$ for which the Central Limit Theorem holds uniformly over the collection of $p$-dimensional hyper-rectangles ?''}. Specifically, we are interested in the normal approximation of suitably scaled versions of the sum $\sum_{i=1}^{n}X_i$ in $\mathcal{R}^p$ uniformly over the class of hyper-rectangles $\mathcal{A}^{re}=\{\prod_{j=1}^{p}[a_j,b_j]\cap\mathcal{R}:-\infty\leq a_j\leq b_j \leq \infty, j=1,\ldots,p\}$, where $X_1,\dots,X_n$ are independent $p-$dimensional random vectors with each having independent and identically distributed (iid) components. We investigate the critical cut-off rate of $\log p$ below which the uniform central limit theorem (CLT) holds and above which it fails. According to some recent results of Chernozukov et al. (2017), it is well known that the CLT holds uniformly over $\mathcal{A}^{re}$ if $\log p=o\big(n^{1/7}\big)$. They also conjectured that for CLT to hold uniformly over $\mathcal{A}^{re}$, the optimal rate is $\log p = o\big(n^{1/3}\big)$. We show instead that under some conditions, the CLT holds uniformly over $\mathcal{A}^{re}$, when $\log p=o\big(n^{1/2}\big)$. More precisely, we show that if $\log p =\epsilon \sqrt{n}$ for some sufficiently small $\epsilon>0$, the normal approximation is valid with an error $\epsilon$, uniformly over $\mathcal{A}^{re}$. Further, we show by an example that the uniform CLT over $\mathcal{A}^{re}$ fails if $\limsup_{t\rightarrow \infty} n^{-(1/2+\delta)} \log p >0$ for some $\delta>0$. Hence the critical rate of the growth of $p$ for the validity of the CLT is given by $\log p=o\big(n^{1/2}\big)$.


title: Almost Hermitian Identities

abstract: We study the local commutation relation between the Lefschetz operator and the exterior differential on an almost complex manifold with a compatible metric. The identity that we obtain generalizes the backbone of the local K\"ahler identities to the setting of almost Hermitian manifolds, allowing for new global results for such manifolds.


title: DeepDrummer : Generating Drum Loops using Deep Learning and a Human in the Loop

abstract: DeepDrummer is a drum loop generation tool that uses active learning to learn the preferences (or current artistic intentions) of a human user from a small number of interactions. The principal goal of this tool is to enable an efficient exploration of new musical ideas. We train a deep neural network classifier on audio data and show how it can be used as the core component of a system that generates drum loops based on few prior beliefs as to how these loops should be structured. We aim to build a system that can converge to meaningful results even with a limited number of interactions with the user. This property enables our method to be used from a cold start situation (no pre-existing dataset), or starting from a collection of audio samples provided by the user. In a proof of concept study with 25 participants, we empirically demonstrate that DeepDrummer is able to converge towards the preference of our subjects after a small number of interactions.


title: Similarities and differences between nickelate and cuprate films grown on a SrTiO$_3$ substrate

abstract: The recent discovery of superconductivity in Sr-doped NdNiO$_2$ films grown on SrTiO$_3$ started a novel field within unconventional superconductivity. To understand the similarities and differences between nickelate and cuprate layers on the same SrTiO$_3$ substrate, here based on the density functional theory we have systematically investigated the structural, electronic, and magnetic properties of NdNiO$_2$/SrTiO$_3$ and CaCuO$_2$/SrTiO$_3$ systems. Our results revealed a strong lattice reconstruction in the case of NdNiO$_2$/SrTiO$_3$, resulting in a polar film, with the surface and interfacial NiO$_2$ layers presenting opposite displacements. However, for CaCuO$_2$/SrTiO$_3$, the distortions of those same two CuO$_2$ layers were in the same direction. In addition, we found this distortion to be approximately independent of the studied range of film thickness for both the nickelate and cuprates films. Furthermore, we also observed a two-dimensional electron gas at the interface between NdNiO$_2$ and SrTiO$_3$, caused by the polar discontinuity, in agreement with recent literature. For NdNiO$_2$/SrTiO$_3$ the two-dimensional electron gas extends over several layers, while for CaCuO$_2$/SrTiO$_3$ this electronic rearrangement is very localized at the interface between CaCuO$_2$ and SrTiO$_3$. The electronic reconstruction found at the interface involves a strong occupation of the Ti $3d_{xy}$ state. In both cases, there is a significant electronic charge transfer from the surface Ni or Cu layers to the Ti interface layer. The interfacial Ni and Cu layer is hole and electron doped, respectively. By introducing magnetism and electronic correlation, we observed that the $d_{3z^2-r^2}$ orbital of Ni becomes itinerant while the same orbital for Cu remains doubly occupied, establishing a clear two- vs one-orbital active framework for the description of these systems.


title: GANBERT: Generative Adversarial Networks with Bidirectional Encoder Representations from Transformers for MRI to PET synthesis

abstract: Synthesizing medical images, such as PET, is a challenging task due to the fact that the intensity range is much wider and denser than those in photographs and digital renderings and are often heavily biased toward zero. Above all, intensity values in PET have absolute significance, and are used to compute parameters that are reproducible across the population. Yet, usually much manual adjustment has to be made in pre-/post- processing when synthesizing PET images, because its intensity ranges can vary a lot, e.g., between -100 to 1000 in floating point values. To overcome these challenges, we adopt the Bidirectional Encoder Representations from Transformers (BERT) algorithm that has had great success in natural language processing (NLP), where wide-range floating point intensity values are represented as integers ranging between 0 to 10000 that resemble a dictionary of natural language vocabularies. BERT is then trained to predict a proportion of masked values images, where its "next sentence prediction (NSP)" acts as GAN discriminator. Our proposed approach, is able to generate PET images from MRI images in wide intensity range, with no manual adjustments in pre-/post- processing. It is a method that can scale and ready to deploy.


title: Varying impacts of letters of recommendation on college admissions: Approximate balancing weights for subgroup effects in observational studies

abstract: In a pilot program during the 2016-17 admissions cycle, the University of California, Berkeley invited many applicants for freshman admission to submit letters of recommendation. We use this pilot as the basis for an observational study of the impact of submitting letters of recommendation on subsequent admission, with the goal of estimating how impacts vary across pre-defined subgroups. Understanding this variation is challenging in observational studies, however, because estimated impacts reflect both actual treatment effect variation and differences in covariate balance across groups. To address this, we develop balancing weights that directly optimize for ``local balance'' within subgroups while maintaining global covariate balance between treated and control units. We then show that this approach has a dual representation as a form of inverse propensity score weighting with a hierarchical propensity score model. In the UC Berkeley pilot study, our proposed approach yields excellent local and global balance, unlike more traditional weighting methods, which fail to balance covariates within subgroups. We find that the impact of letters of recommendation increases with the predicted probability of admission, with mixed evidence of differences for under-represented minority applicants.


title: tf-Darshan: Understanding Fine-grained I/O Performance in Machine Learning Workloads

abstract: Machine Learning applications on HPC systems have been gaining popularity in recent years. The upcoming large scale systems will offer tremendous parallelism for training through GPUs. However, another heavy aspect of Machine Learning is I/O, and this can potentially be a performance bottleneck. TensorFlow, one of the most popular Deep-Learning platforms, now offers a new profiler interface and allows instrumentation of TensorFlow operations. However, the current profiler only enables analysis at the TensorFlow platform level and does not provide system-level information. In this paper, we extend TensorFlow Profiler and introduce tf-Darshan, both a profiler and tracer, that performs instrumentation through Darshan. We use the same Darshan shared instrumentation library and implement a runtime attachment without using a system preload. We can extract Darshan profiling data structures during TensorFlow execution to enable analysis through the TensorFlow profiler. We visualize the performance results through TensorBoard, the web-based TensorFlow visualization tool. At the same time, we do not alter Darshan's existing implementation. We illustrate tf-Darshan by performing two case studies on ImageNet image and Malware classification. We show that by guiding optimization using data from tf-Darshan, we increase POSIX I/O bandwidth by up to 19% by selecting data for staging on fast tier storage. We also show that Darshan has the potential of being used as a runtime library for profiling and providing information for future optimization.


title: GANDALF: Generative Adversarial Networks with Discriminator-Adaptive Loss Fine-tuning for Alzheimer's Disease Diagnosis from MRI

abstract: Positron Emission Tomography (PET) is now regarded as the gold standard for the diagnosis of Alzheimer's Disease (AD). However, PET imaging can be prohibitive in terms of cost and planning, and is also among the imaging techniques with the highest dosage of radiation. Magnetic Resonance Imaging (MRI), in contrast, is more widely available and provides more flexibility when setting the desired image resolution. Unfortunately, the diagnosis of AD using MRI is difficult due to the very subtle physiological differences between healthy and AD subjects visible on MRI. As a result, many attempts have been made to synthesize PET images from MR images using generative adversarial networks (GANs) in the interest of enabling the diagnosis of AD from MR. Existing work on PET synthesis from MRI has largely focused on Conditional GANs, where MR images are used to generate PET images and subsequently used for AD diagnosis. There is no end-to-end training goal. This paper proposes an alternative approach to the aforementioned, where AD diagnosis is incorporated in the GAN training objective to achieve the best AD classification performance. Different GAN lossesare fine-tuned based on the discriminator performance, and the overall training is stabilized. The proposed network architecture and training regime show state-of-the-art performance for three- and four- class AD classification tasks.


title: sputniPIC: an Implicit Particle-in-Cell Code for Multi-GPU Systems

abstract: Large-scale simulations of plasmas are essential for advancing our understanding of fusion devices, space, and astrophysical systems. Particle-in-Cell (PIC) codes have demonstrated their success in simulating numerous plasma phenomena on HPC systems. Today, flagship supercomputers feature multiple GPUs per compute node to achieve unprecedented computing power at high power efficiency. PIC codes require new algorithm design and implementation for exploiting such accelerated platforms. In this work, we design and optimize a three-dimensional implicit PIC code, called sputniPIC, to run on a general multi-GPU compute node. We introduce a particle decomposition data layout, in contrast to domain decomposition on CPU-based implementations, to use particle batches for overlapping communication and computation on GPUs. sputniPIC also natively supports different precision representations to achieve speed up on hardware that supports reduced precision. We validate sputniPIC through the well-known GEM challenge and provide performance analysis. We test sputniPIC on three multi-GPU platforms and report a 200-800x performance improvement with respect to the sputniPIC CPU OpenMP version performance. We show that reduced precision could further improve performance by 45% to 80% on the three platforms. Because of these performance improvements, on a single node with multiple GPUs, sputniPIC enables large-scale three-dimensional PIC simulations that were only possible using clusters.


title: Matching for random systems with an application to minimal weight expansions

abstract: We extend the notion of matching for one-dimensional dynamical systems to random matching for random dynamical systems on an interval. We prove that for a large family of piecewise affine random systems of the interval the property of random matching implies that any invariant density is piecewise constant. We further introduce a one-parameter family of random dynamical systems that produce signed binary expansions of numbers in the interval [-1,1]. This family has random matching for Lebesgue almost every parameter. We use this to prove that the frequency of the digit 0 in the associated signed binary expansions never exceeds 1/2.


title: The VISCACHA survey -- II. Structure of star clusters in the Magellanic Clouds periphery

abstract: We provide a homogeneous set of structural parameters of 83 star clusters located at the periphery of the Small Magellanic Cloud (SMC) and the Large Magellanic Cloud (LMC). The clusters' stellar density and surface brightness profiles were built from deep, AO assisted optical images, and uniform analysis techniques. The structural parameters were obtained from King and Elson et al. model fittings. Integrated magnitudes and masses (for a subsample) are also provided. The sample contains mostly low surface brightness clusters with distances between 4.5 and 6.5 kpc and between 1 and 6.5 kpc from the LMC and SMC centres, respectively. We analysed their spatial distribution and structural properties, comparing them with those of inner clusters. Half-light and Jacobi radii were estimated, allowing an evaluation of the Roche volume tidal filling. We found that: (i) for our sample of LMC clusters, the tidal radii are, on average, larger than those of inner clusters from previous studies; (ii) the core radii dispersion tends to be greater for LMC clusters located towards the southwest, with position angles of $\sim$200 degrees and about $\sim$5 degrees from the LMC centre, i.e., those LMC clusters nearer to the SMC; (iii) the core radius evolution for clusters with known age is similar to that of inner clusters; (iv) SMC clusters with galactocentric distances closer than 4 kpc are overfilling; (v) the recent Clouds collision did not leave marks on the LMC clusters' structure that our analysis could reveal.


title: $\alpha$-decay properties of even-even superheavy nuclei

abstract: In this paper, we have calculated the $\alpha$-decay half-lives of superheavy nuclei with $106 \leq Z \leq 126$ and a neutron number of $150 \leq N \leq 200$ within proximity potentials and deformed-spherical Coulomb potentials by using $WS4$ $\alpha$-decay energy and the semi-classical Wentzel-Kramers-Brillouin approximation for penetration probability. Besides, we have included the preformation factor within the cluster-formation model. We have investigated magic numbers and submagic numbers in the mentioned region; with high probability 162, 178, and 184 are predicted as neutron magic numbers. We also have confirmed that there is good agreement between our predicted half-lives and the ones obtained from semiempirical relationships such as Royer, VSS, UDL, and SemFIS2.


title: Connected Incomplete Preferences

abstract: The standard model of choice in economics is the maximization of a complete and transitive preference relation over a fixed set of alternatives. While completeness of preferences is usually regarded as a strong assumption, weakening it requires care to ensure that the resulting model still has enough structure to yield interesting results. This paper takes a step in this direction by studying the class of "connected preferences", that is, preferences that may fail to be complete but have connected maximal domains of comparability. We offer four new results. Theorem 1 identifies a basic necessary condition for a continuous preference to be connected in the sense above, while Theorem 2 provides sufficient conditions. Building on the latter, Theorem 3 characterizes the maximal domains of comparability. Finally, Theorem 4 presents conditions that ensure that maximal domains are arc-connected.


title: Ultraviolet Line Identifications and Spectral Formation Near Max-Light in Type Ia Supernovae 2011fe

abstract: We present ultraviolet line identifications of near maximum-light HST observations of SN 2011fe using synthetic spectra generated from both SYNOW and $\texttt{PHOENIX}$. We find the spectrum to be dominated by blends of iron group elements Fe, Co, and Ni (as expected due to heavy line blanketing by these elements in the UV) and for the first time identify lines from C IV and Si IV in a supernova spectrum. We also find that classical delayed detonation models of Type Ia supernovae are able to accurately reproduce the flux levels of SN 2011fe in the UV. Further analysis reveals that photionization edges play an important role in feature formation in the far-UV, and that temperature variations in the outer layers of the ejecta significantly alter the Fe III/Fe II ratio producing large flux changes in the far-UV and velocity shifts in mid-UV features. SN 2011fe is the best observed core-normal SNe Ia, therefore analysis its of UV spectra shows the power of UV spectra in discriminating between different metalicities and progenitor scenarios of Type Ia supernovae, due to the fact that the UV probes the outermost layers of the Type Ia supernova, which are most sensitive to metalicity and progenitor variations.


title: Comparison of Model Predictive and Reinforcement Learning Methods for Fault Tolerant Control

abstract: A desirable property in fault-tolerant controllers is adaptability to system changes as they evolve during systems operations. An adaptive controller does not require optimal control policies to be enumerated for possible faults. Instead it can approximate one in real-time. We present two adaptive fault-tolerant control schemes for a discrete time system based on hierarchical reinforcement learning. We compare their performance against a model predictive controller in presence of sensor noise and persistent faults. The controllers are tested on a fuel tank model of a C-130 plane. Our experiments demonstrate that reinforcement learning-based controllers perform more robustly than model predictive controllers under faults, partially observable system models, and varying sensor noise levels.


title: Modified Gravity (MOG) and Heavy Neutron Star in Mass Gap

abstract: The modified gravity (MOG) theory is applied to the gravitational wave binary merger GW190814 to demonstrate that the modified Tolman-Oppenheimer-Volkoff equation for a neutron star can produce a mass $M=2.6 -2.7 M_\odot$, allowing for the binary secondary component to be identified as a heavy neutron star in the hypothesized mass gap $2.5 - 5 M_\odot$.


title: Compact Spaces with a $P$-base

abstract: In the paper, we investigate (scattered) compact spaces with a $P$-base for some poset $P$. More specifically, we prove that, under the assumption $\omega_1<\mathfrak{b}$, any compact space with an $\omega^\omega$-base is first-countable and any scattered compact space with an $\omega^\omega$-base is countable. These give positive solutions to Problems 8.6.9 and 8.7.7 in \cite{Banakh2019}. Using forcing, we also prove that in a model of $\omega_1<\mathfrak{b}$, there is a non-first-countable compact space with a $P$-base for some poset $P$ with calibre~$\omega_1$.


title: Reduction and Coherent States

abstract: We apply a quantum version of dimensional reduction to Gaussian coherent states in Bargmann space to obtain squeezed states on complex projective spaces. This leads to a definition of a family of squeezed spin states with excellent semi-classical properties, governed by a symbol calculus. We prove semiclassical norm estimates and a propagation result.


title: Fault-Tolerant Control of Degrading Systems with On-Policy Reinforcement Learning

abstract: We propose a novel adaptive reinforcement learning control approach for fault tolerant control of degrading systems that is not preceded by a fault detection and diagnosis step. Therefore, \textit{a priori} knowledge of faults that may occur in the system is not required. The adaptive scheme combines online and offline learning of the on-policy control method to improve exploration and sample efficiency, while guaranteeing stable learning. The offline learning phase is performed using a data-driven model of the system, which is frequently updated to track the system's operating conditions. We conduct experiments on an aircraft fuel transfer system to demonstrate the effectiveness of our approach.


title: Optical Properties of Dilute CuAl$_{1-x}$Fe$_x$O$_2$ Delafossite Alloys

abstract: For powder samples of CuAl$_{1-x}$Fe$_x$O$_2$ ($x =$ 0, 0.01, 0.05, and 0.1), measured optical properties are compared with model simulations and phonon spectra are compared with simulations based on weighted dynamical matrix approach.


title: A brief introduction to observational entropy

abstract: In the past several years, observational entropy has been developed as both a (time-dependent) quantum generalization of Boltzmann entropy, and as a rather general framework to encompass classical and quantum equilibrium and non-equilibrium coarse-grained entropy. In this paper we review the construction, interpretation, most important properties, and some applications of this framework. The treatment is self-contained and relatively pedagogical, aimed at a broad class of researchers.


title: Characterizing viscoelastic materials via ensemble-based data assimilation of bubble collapse observations

abstract: Viscoelastic material properties at high strain rates are needed to model many biological and medical systems. Bubble cavitation can induce such strain rates, and the resulting bubble dynamics are sensitive to the material properties. Thus, in principle, these properties can be inferred via measurements of the bubble dynamics. Estrada et al. (2018) demonstrated such bubble-dynamic high-strain-rate rheometry by using least-squares shooting to minimize the difference between simulated and experimental bubble radius histories. We generalize their technique to account for additional uncertainties in the model, initial conditions, and material properties needed to uniquely simulate the bubble dynamics. Ensemble-based data assimilation minimizes the computational expense associated with the bubble cavitation model. We test an ensemble Kalman filter (EnKF), an iterative ensemble Kalman smoother (IEnKS), and a hybrid ensemble-based 4D--Var method (En4D--Var) on synthetic data, assessing their estimations of the viscosity and shear modulus of a Kelvin--Voigt material. Results show that En4D--Var and IEnKS provide better moduli estimates than EnKF. Applying these methods to the experimental data of Estrada et al. (2018) yields similar material property estimates to those they obtained, but provides additional information about uncertainties. In particular, the En4D--Var yields lower viscosity estimates for some experiments, and the dynamic estimators reveal a potential mechanism that is unaccounted for in the model, whereby the viscosity is reduced in some cases due to material damage occurring at bubble collapse.


title: Meshless Approximation and Helmholtz-Hodge Decomposition of Vector Fields

abstract: The analysis of vector fields is crucial for the understanding of several physical phenomena, such as natural events (e.g., analysis of waves), diffusive processes, electric and electromagnetic fields. While previous work has been focused mainly on the analysis of 2D or 3D vector fields on volumes or surfaces, we address the meshless analysis of a vector field defined on an arbitrary domain, without assumptions on its dimension and discretisation. The meshless approximation of the Helmholtz-Hodge decomposition of a vector field is achieved by expressing the potential of its components as a linear combination of radial basis functions and by computing the corresponding conservative, irrotational, and harmonic components as solution to a least-squares or to a differential problem. To this end, we identify the conditions on the kernel of the radial basis functions that guarantee the existence of their derivatives. Finally, we demonstrate our approach on 2D and 3D vector fields measured by sensors or generated through simulation.


title: Unitary representations of the Cherednik algebra: $V^*$-homology

abstract: We give a non-negative combinatorial formula, in terms of Littlewood-Richardson numbers, for the homology of the unitary representations of the cyclotomic rational Cherednik algebra, and as a consequence, for the graded Betti numbers for the ideals of a class of subspace arrangements arising from the reflection arrangements of complex reflection groups.


title: On optical transmittance of ultra diluted gas

abstract: The paper proposes a model of optical transmittance of ultra diluted gas taking into account gas particles non-locality, the quantum effect of wave function spreading derived from solving the Schr\"odinger equation for a free particle. A significant increase in the transmittance of such gas is envisaged as compared to the classical predictions. Some quantitative and qualitative consequences of the model are indicated and falsifying experiments are proposed. The classic Beer-Lambert law equation within range of its applicability is derived from the model. Remarks to some astrophysical phenomena and possible interpretations of Quantum Mechanics are made. An experiment consistent with the predictions of this model is referenced.


title: A Framework of Hierarchical Attacks to Network Controllability

abstract: Network controllability robustness reflects how well a networked dynamical system can maintain its controllability against destructive attacks. This paper investigates the network controllability robustness from the perspective of a malicious attack. A framework of hierarchical attack is proposed, by means of edge- or node-removal attacks. Edges (or nodes) in a target network are classified hierarchically into categories, with different priorities to attack. The category of critical edges (or nodes) has the highest priority to be selected for attack. Extensive experiments on nine synthetic networks and nine real-world networks show the effectiveness of the proposed hierarchical attack strategies for destructing the network controllability. From the protection point of view, this study suggests that the critical edges and nodes should be hidden from the attackers. This finding helps better understand the network controllability and better design robust networks.


title: Adaptive music: Automated music composition and distribution

abstract: Creativity, or the ability to produce new useful ideas, is commonly associated to the human being; but there are many other examples in nature where this phenomenon can be observed. Inspired by this fact, in engineering and particularly in computational sciences, many different models have been developed to tackle a number of problems. Composing music, a form of art broadly present along the human history, is the main topic addressed in this thesis. Taking advantage of the kind of ideas that bring diversity and creativity to nature and computation, we present Melomics: an algorithmic composition method based on evolutionary search. The solutions have a genetic encoding based on formal grammars and these are interpreted in a complex developmental process followed by a fitness assessment, to produce valid music compositions in standard formats. The system has exhibited a high creative power and versatility to produce music of different types and it has been tested, proving on many occasions the outcome to be indistinguishable from the music made by human composers. The system has also enabled the emergence of a set of completely novel applications: from effective tools to help anyone to easily obtain the precise music that they need, to radically new uses, such as adaptive music for therapy, exercise, amusement and many others. It seems clear that automated composition is an active research area and that countless new uses will be discovered.


title: Approximation in (Poly-) Logarithmic Space

abstract: We develop new approximation algorithms for classical graph and set problems in the RAM model under space constraints. As one of our main results, we devise an algorithm for d-Hitting Set that runs in time n^{O(d^2 + d/\epsilon})}, uses O((d^2 + d/\epsilon) log n) bits of space, and achieves an approximation ratio of O((d/{\epsilon}) n^{\epsilon}) for any positive \epsilon \leq 1 and any natural number d. In particular, this yields a factor-O(log n) approximation algorithm which runs in time n^{O(log n)} and uses O(log^2 n) bits of space (for constant d). As a corollary, we obtain similar bounds for Vertex Cover and several graph deletion problems. For bounded-multiplicity problem instances, one can do better. We devise a factor-2 approximation algorithm for Vertex Cover on graphs with maximum degree \Delta, and an algorithm for computing maximal independent sets which both run in time n^{O(\Delta)} and use O(\Delta log n) bits of space. For the more general d-Hitting Set problem, we devise a factor-d approximation algorithm which runs in time n^{O(d {\delta}^2)} and uses O(d {\delta}^2 log n) bits of space on set families where each element appears in at most \delta sets. For Independent Set restricted to graphs with average degree d, we give a factor-(2d) approximation algorithm which runs in polynomial time and uses O(log n) bits of space. We also devise a factor-O(d^2) approximation algorithm for Dominating Set on d-degenerate graphs which runs in time n^{O(log n)} and uses O(log^2 n) bits of space. For d-regular graphs, we show how a known randomized factor-O(log d) approximation algorithm can be derandomized to run in time n^{O(1)} and use O(log n) bits of space. Our results use a combination of ideas from the theory of kernelization, distributed algorithms and randomized algorithms.