title: The Structure of Radiatively Inefficient Black Hole Accretion Flows

abstract: We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting.


title: PIRK: Scalable Interval Reachability Analysis for High-Dimensional Nonlinear Systems

abstract: Reachability analysis is a critical tool for the formal verification of dynamical systems and the synthesis of controllers for them. Due to their computational complexity, many reachability analysis methods are restricted to systems with relatively small dimensions. One significant reason for such limitation is that those approaches, and their implementations, are not designed to leverage parallelism. They use algorithms that are designed to run serially within one compute unit and they can not utilize widely-available high-performance computing (HPC) platforms such as many-core CPUs, GPUs and Cloud-computing services. This paper presents PIRK, a tool to efficiently compute reachable sets for general nonlinear systems of extremely high dimensions. PIRK has been tested on several systems, with state dimensions ranging from ten up to 4 billion. The scalability of PIRK's parallel implementations is found to be highly favorable.


title: Entropy numbers and Marcinkiewicz-type discretization theorem

abstract: This paper studies the behavior of the entropy numbers of classes of functions with bounded integral norms from a given finite dimensional linear subspace. Upper bounds of these entropy numbers in the uniform norm are obtained and applied to establish a Marcinkiewicz type discretization theorem for integral norms of functions from a given finite dimensional subspace.


title: Flexible Multiple Access Enabling Low-Latency Communications: Introducing NOMA-R

abstract: Various verticals in 5G and beyond (B5G) networks require very stringent latency guarantees, while at the same time envisioning massive connectivity. As a result, choosing the optimal multiple access (MA) technique to achieve low latency is a key enabler of B5G. In particular, this issue is more acute in uplink transmissions due to the potentially high number of collisions. On this premise, in the present contribution we discuss the issue of delay-sensitive uplink connectivity using optimized MA techniques; to this end, we perform a comparative analysis of various MA approaches with respect to the achievable effective capacity (EC). As opposed to standard rate (PHY) or throughput (MAC) analyses, we propose the concept of the effective capacity as a suitable metric for characterizing jointly PHY-MAC layer delays. The palette of investigated MA approaches includes standard orthogonal MA (OMA) and power domain non orthogonal MA (NOMA) in uplink scenarios, both considering random pairing and optimized pairing alternatives. It further extends to encompass a recently proposed third alternative, referred to as NOMA-Relevant (NOMA-R), which extends OMA and NOMA approaches by flexibly selecting the MA technique. We show that optimizing both user pairing and MA selection increases the network EC, especially when stringent delay constraints are in place; thus a flexible MA is a potentially preferable strategy for future low latency applications


title: Gevrey regularity for the Vlasov-Poisson system

abstract: We prove propagation of $\frac{1}{s}$-Gevrey regularity $(s\in(0,1))$ for the Vlasov-Poisson system on $\mathbb{T}^d$ using a Fourier space method in analogy to the results proved for the 2D Euler system in \cite{KV} and \cite{LO}. More precisely, we give a quantitative estimate for the growth in time of the $\frac{1}{s}$-Gevrey norm for the solution of the system in terms of the force field and the gradient in the velocity variable of the distribution of matter. As an application, we show global existence of $\frac{1}{s}$-Gevrey solutions ($s\in (0,1)$) for the Vlasov-Poisson system in $\mathbb{T}^3$. Furthermore, the propagation of Gevrey regularity can be easily modified to obtain the same result in $\mathbb{R}^d$. In particular, this implies global existence of analytic $(s=1)$ and $\frac{1}{s}$-Gevrey solutions ($s\in (0,1)$) for the Vlasov-Poisson system in $\mathbb{R}^3$.


title: Automatic weak imposition of free slip boundary conditions via Nitsche's method: application to nonlinear problems in geodynamics

abstract: Imposition of free slip boundary conditions in science and engineering simulations presents a challenge when the simulation domain is non-trivial. Inspired by recent progress in symbolic computation of discontinuous Galerkin finite element methods, we present a symmetric interior penalty form of Nitsche's method to weakly impose these slip boundary conditions and present examples of its use in the Stokes subsystem motivated by problems in geodynamics. We compare numerical results with well established benchmark problems. We also examine performance of the method with iterative solvers.


title: Bounded Incentives in Manipulating the Probabilistic Serial Rule

abstract: The Probabilistic Serial mechanism is well-known for its desirable fairness and efficiency properties. It is one of the most prominent protocols for the random assignment problem. However, Probabilistic Serial is not incentive-compatible, thereby these desirable properties only hold for the agents' declared preferences, rather than their genuine preferences. A substantial utility gain through strategic behaviors would trigger self-interested agents to manipulate the mechanism and would subvert the very foundation of adopting the mechanism in practice. In this paper, we characterize the extent to which an individual agent can increase its utility by strategic manipulation. We show that the incentive ratio of the mechanism is $\frac{3}{2}$. That is, no agent can misreport its preferences such that its utility becomes more than 1.5 times of what it is when reports truthfully. This ratio is a worst-case guarantee by allowing an agent to have complete information about other agents' reports and to figure out the best response strategy even if it is computationally intractable in general. To complement this worst-case study, we further evaluate an agent's utility gain on average by experiments. The experiments show that an agent' incentive in manipulating the rule is very limited. These results shed some light on the robustness of Probabilistic Serial against strategic manipulation, which is one step further than knowing that it is not incentive-compatible.


title: The Rockerverse: Packages and Applications for Containerization with R

abstract: The Rocker Project provides widely used Docker images for R across different application scenarios. This article surveys downstream projects that build upon the Rocker Project images and presents the current state of R packages for managing Docker images and controlling containers. These use cases cover diverse topics such as package development, reproducible research, collaborative work, cloud-based data processing, and production deployment of services. The variety of applications demonstrates the power of the Rocker Project specifically and containerisation in general. Across the diverse ways to use containers, we identified common themes: reproducible environments, scalability and efficiency, and portability across clouds. We conclude that the current growth and diversification of use cases is likely to continue its positive impact, but see the need for consolidating the Rockerverse ecosystem of packages, developing common practices for applications, and exploring alternative containerisation software.


title: Binary Classification from Positive Data with Skewed Confidence

abstract: Positive-confidence (Pconf) classification [Ishida et al., 2018] is a promising weakly-supervised learning method which trains a binary classifier only from positive data equipped with confidence. However, in practice, the confidence may be skewed by bias arising in an annotation process. The Pconf classifier cannot be properly learned with skewed confidence, and consequently, the classification performance might be deteriorated. In this paper, we introduce the parameterized model of the skewed confidence, and propose the method for selecting the hyperparameter which cancels out the negative impact of skewed confidence under the assumption that we have the misclassification rate of positive samples as a prior knowledge. We demonstrate the effectiveness of the proposed method through a synthetic experiment with simple linear models and benchmark problems with neural network models. We also apply our method to drivers' drowsiness prediction to show that it works well with a real-world problem where confidence is obtained based on manual annotation.


title: Inverse FIP Effects in Giant Solar Flares Found from Earth X-Ray Albedo with Suzaku/XIS

abstract: We report X-ray spectroscopic results for four giant solar flares occurred on 2005 September 7 (X17.0), 2005 September 8 (X5.4), 2005 September 9 (X6.2), and 2006 December 5 (X9.0), obtained from Earth albedo data with the X-ray imaging spectrometer (XIS) onboard Suzaku. The good energy resolution of the XIS (FWHM$\sim$100 eV) enables us to separate a number of line-like features and detect the underlying continuum emission. These features include Si He$\alpha$, Si Ly$\alpha$, S He$\alpha$, S Ly$\alpha$, Ar He$\alpha$, and Ca He$\alpha$ originating from solar flares as well as fluorescent Ar K$\alpha$ and Ar K$\beta$ from the Earth atmosphere. Absolute elemental abundances (X/H) averaged over the four flares are obtained to be $\sim$2.0 (Ca), $\sim$0.7 (Si), $\sim$0.3 (S), and $\sim$0.9 (Ar) at around flare peaks. This abundance pattern is similar to those of active stars' coronae showing inverse first ionization potential (i-FIP) effects, i.e., elemental abundances decrease with decreasing FIP with a turnover at the low end of FIP. The abundances are almost constant during the flares, with an exception of Si which increases by a factor of $\sim$2 in the decay phase. The evolution of the Si abundance is consistent with a picture that the i-FIP plasma originates from the chromosphere evaporation and then mixes with the surrounding low-FIP biased materials. Flare-to-flare abundance varied by a factor of 2, agreeing with past observations of solar flares. Finally, we emphasize that Earth albedo data acquired by X-ray astronomy satellites like Suzaku and XRISM can significantly contribute to studies of solar physics.


title: Closed Form Solutions For The Quantum LRS Bianchi IX and VIII Wheeler DeWitt Equations For An Arbitrary Hartle-Hawking Ordering Parameter

abstract: In this paper we present the quantum LRS Bianchi IX models or, as they may be better known, the Taub-NUT models in addition to the quantum LRS Bianchi VIII models, both coupled to a stiff matter source. We solve their Wheeler DeWitt equations for an arbitrary Hartle-Hawking ordering parameter using separation of variables. Afterwards we construct a superposition of their respective wave functions and discuss some of the interesting qualitative properties they possess.


title: The Formation Height of Millimeter-wavelength Emission in the Solar Chromosphere

abstract: In the past few years, the ALMA radio telescope has become available for solar observations. ALMA diagnostics of the solar atmosphere are of high interest because of the theoretically expected linear relationship between the brightness temperature at mm wavelengths and the local gas temperature in the solar atmosphere. Key for the interpretation of solar ALMA observations is understanding where in the solar atmosphere the ALMA emission originates. Recent theoretical studies have suggested that ALMA bands at 1.2 (band 6) and 3 mm (band 3) form in the middle and upper chromosphere at significantly different heights. We study the formation of ALMA diagnostics using a 2.5D radiative MHD model that includes the effects of ion-neutral interactions (ambipolar diffusion) and non-equilibrium ionization of hydrogen and helium. Our results suggest that in active regions and network regions, observations at both wavelengths most often originate from similar heights in the upper chromosphere, contrary to previous results. Non-equilibrium ionization increases the opacity in the chromosphere so that ALMA mostly observe spicules and fibrils along the canopy fields. We combine these modeling results with observations from IRIS, SDO and ALMA to suggest a new interpretation for the recently reported "dark chromospheric holes", regions of very low temperatures in the chromosphere.


title: Green equivalences in equivariant mathematics

abstract: We establish Green equivalences for all Mackey 2-functors, without assuming Krull-Schmidt. By running through the examples of Mackey 2-functors, we recover all variants of the Green equivalence and Green correspondence known in representation theory and obtain new ones in several other contexts. Such applications include equivariant stable homotopy theory in topology and equivariant sheaves in geometry.


title: Caustics of weakly Lagrangian distributions

abstract: We study semiclassical sequences of distributions $u_h$ associated to a Lagrangian submanifold of phase space $\lag \subset T^*X$. If $u_h$ is a semiclassical Lagrangian distribution, which concentrates at a maximal rate on $\lag,$ then the asymptotics of $u_h$ are well-understood by work of Arnol'd, provided $\lag$ projects to $X$ with a stable Lagrangian singularity. We establish sup-norm estimates on $u_h$ under much more general hypotheses on the rate at which it is concentrating on $\lag$ (again assuming a stable projection). These estimates apply to sequences of eigenfunctions of integrable and KAM Hamiltonians.


title: Modelling and Quantifying Membership Information Leakage in Machine Learning

abstract: Machine learning models have been shown to be vulnerable to membership inference attacks, i.e., inferring whether individuals' data have been used for training models. The lack of understanding about factors contributing success of these attacks motivates the need for modelling membership information leakage using information theory and for investigating properties of machine learning models and training algorithms that can reduce membership information leakage. We use conditional mutual information leakage to measure the amount of information leakage from the trained machine learning model about the presence of an individual in the training dataset. We devise an upper bound for this measure of information leakage using Kullback--Leibler divergence that is more amenable to numerical computation. We prove a direct relationship between the Kullback--Leibler membership information leakage and the probability of success for a hypothesis-testing adversary examining whether a particular data record belongs to the training dataset of a machine learning model. We show that the mutual information leakage is a decreasing function of the training dataset size and the regularization weight. We also prove that, if the sensitivity of the machine learning model (defined in terms of the derivatives of the fitness with respect to model parameters) is high, more membership information is potentially leaked. This illustrates that complex models, such as deep neural networks, are more susceptible to membership inference attacks in comparison to simpler models with fewer degrees of freedom. We show that the amount of the membership information leakage is reduced by $\mathcal{O}(\log^{1/2}(\delta^{-1})\epsilon^{-1})$ when using Gaussian $(\epsilon,\delta)$-differentially-private additive noises.

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


title: On highly accurate computations of the Coulomb two-center systems with unit charges

abstract: Results of our recent highly accurate computations of the Coulomb two-center systems with the unit electrical charges $X^{+} X^{+} e^{-}$ and $X^{+} Y^{+} e^{-}$ are discussed. In particular, we have determined (to very high accuracy) the total energies of the ground $1 s \sigma-$states in the two-center adiabatic (or molecular) H$_{2}^{+}$, D$_2^{+}$, HD$^{+}$, HT$^{+}$, T$^{+}_{2}$ and DT$^{+}$ ions. In these computations we applied the new masses of hydrogen isotopes, which have been measured in the recent high-energy experiments. We also derived (and tested) our accurate mass-interpolation formula for the total energies of the model two-center $X^{+} X^{+} e^{-}$ ions with very heavy masses of the point $X^{+}$ particles ($M_X \ge$ 100,000 $m_e$). By using this formula we analyze the overall accuracy and validity of the adiabatic two-center approximation in application to the three-body atomic (or Coulomb) systems.


title: Connection coefficients for ultraspherical polynomials with argument doubling and generalized bispectrality

abstract: We start by presenting a generalization of a discrete wave equation that is particularly satisfied by the entries of the matrix coefficients of the refinement equation corresponding to the multiresolution analysis of Alpert. The entries are in fact functions of two discrete variables and they can be expressed in terms of the Legendre polynomials. Next, we generalize these functions to the case of the ultraspherical polynomials and show that these new functions obey two generalized eigenvalue problems in each of the two discrete variables, which constitute a generalized bispectral problem. At the end, we make some connections to other problems.


title: The image of a point charge in an infinite conducting cylinder

abstract: The electrostatics problem of a point charge next to a conducting plane is best solved by placing an image charge placed on the opposite side. For a charge between two parallel planes this can be solved with image charges outside the planes at evenly spaced intervals moving out to infinity. What is the corresponding image of a point charge is when placed on the axis of a cylinder?. The potential of a point charge in a cylinder is well known and may expressed in many forms involving integrals or series of Bessel functions, but none of which elude to an image. In fact the image is complex (in both senses), consisting of infinitely many rings on a disk with some surface charge distribution. This manuscript attempts to describe the image as accurately as possible, and in doing so finds simple accurate approximations for the potential.


title: Treatment effect estimation with disentangled latent factors

abstract: Much research has been devoted to the problem of estimating treatment effects from observational data; however, most methods assume that the observed variables only contain confounders, i.e., variables that affect both the treatment and the outcome. Unfortunately, this assumption is frequently violated in real-world applications, since some variables only affect the treatment but not the outcome, and vice versa. Moreover, in many cases only the proxy variables of the underlying confounding factors can be observed. In this work, we first show the importance of differentiating confounding factors from instrumental and risk factors for both average and conditional average treatment effect estimation, and then we propose a variational inference approach to simultaneously infer latent factors from the observed variables, disentangle the factors into three disjoint sets corresponding to the instrumental, confounding, and risk factors, and use the disentangled factors for treatment effect estimation. Experimental results demonstrate the effectiveness of the proposed method on a wide range of synthetic, benchmark, and real-world datasets.


title: Dark-soliton molecules in an exciton-polariton superfluid

abstract: The general theory of dark solitons relies on repulsive interactions and therefore predicts the impossibility to form dark-soliton bound states. One important exception to this prediction is the observation of bound solitons in non-local nonlinear media. Here, we report that exciton-polariton superfluids can also sustain dark-soliton molecules although the interactions are fully local. With a novel all optical technique, we create two dark solitons and bind them to each other to form an unconventional dark-soliton molecule. We demonstrate that the stability of this structure and the separation distance between two dark-solitons is tightly connected to the driven-dissipative nature of the polariton fluid.


title: On the dual positive cones and the algebraicity of a compact K\"ahler manifold

abstract: We investigate the algebraicity of compact K\"ahler manifolds admitting a positive rational Hodge class of bidimension $(1,1)$. We prove that if the dual K\"ahler cone of a compact K\"ahler manifold $X$ contains a rational class as an interior point, then its Albanese variety is projective. As a consequence, we answer the Oguiso--Peternell problem for Ricci-flat compact K\"ahler manifolds. We also study related algebraicity problems for threefolds.


title: Regularization Helps with Mitigating Poisoning Attacks: Distributionally-Robust Machine Learning Using the Wasserstein Distance

abstract: We use distributionally-robust optimization for machine learning to mitigate the effect of data poisoning attacks. We provide performance guarantees for the trained model on the original data (not including the poison records) by training the model for the worst-case distribution on a neighbourhood around the empirical distribution (extracted from the training dataset corrupted by a poisoning attack) defined using the Wasserstein distance. We relax the distributionally-robust machine learning problem by finding an upper bound for the worst-case fitness based on the empirical sampled-averaged fitness and the Lipschitz-constant of the fitness function (on the data for given model parameters) as regularizer. For regression models, we prove that this regularizer is equal to the dual norm of the model parameters. We use the Wine Quality dataset, the Boston Housing Market dataset, and the Adult dataset for demonstrating the results of this paper.

categories: cs.LG cs.CR eess.SP math.OC math.ST stat.ML stat.TH


title: Constructing a variational quasi-reversibility method for a Cauchy problem for elliptic equations

abstract: In the recent developments of regularization theory for inverse and ill-posed problems, a variational quasi-reversibility (QR) method has been designed to solve a class of time-reversed quasi-linear parabolic problems. Known as a PDE-based approach, this method relies on adding a suitable perturbing operator to the original problem and consequently, on gaining the corresponding fine stabilized operator, which leads us to a forward-like problem. In this work, we establish new conditional estimates for such operators to solve a prototypical Cauchy problem for elliptic equations. This problem is based on the stationary case of the inverse heat conduction problem, where one wants to identify the heat distribution in a certain medium, given the partial boundary data. Using the new QR method, we obtain a second-order initial value problem for a wave-type equation, whose weak solvability can be deduced using a priori estimates and compactness arguments. Weighted by a Carleman-like function, a new type of energy estimates is explored in a variational setting when we investigate the H\"older convergence rate of the proposed scheme. Besides, a linearized version of this scheme is analyzed. Numerical examples are provided to corroborate our theoretical analysis.


title: The Indian Chefs Process

abstract: This paper introduces the Indian Chefs Process (ICP), a Bayesian nonparametric prior on the joint space of infinite directed acyclic graphs (DAGs) and orders that generalizes Indian Buffet Processes. As our construction shows, the proposed distribution relies on a latent Beta Process controlling both the orders and outgoing connection probabilities of the nodes, and yields a probability distribution on sparse infinite graphs. The main advantage of the ICP over previously proposed Bayesian nonparametric priors for DAG structures is its greater flexibility. To the best of our knowledge, the ICP is the first Bayesian nonparametric model supporting every possible DAG. We demonstrate the usefulness of the ICP on learning the structure of deep generative sigmoid networks as well as convolutional neural networks.


title: A Sequential Constraint Method for Solving Variational Inequality over the Intersection of Fixed Point Sets

abstract: We consider the variational inequality problem over the intersection of fixed point sets of firmly nonexpansive operators. In order to solve the problem, we present an algorithm and subsequently show the strong convergence of the generated sequence to the solution of the considered problem.


title: Creating High Quality All-Sky Visualizations of Astronomy Image Data Sets: HiPS and Montage

abstract: We describe a case study to use the Montage image mosaic engine to create maps of the ALLWISE image data set in the Hierarchical Progressive Survey (HiPS) sky-tesselation scheme. Our approach demonstrates that Montage reveals the science content of infrared images in greater detail than has hitherto been possible in HiPS maps. The approach exploits two unique (to our knowledge) characteristics of the Montage image mosaic engine: background modeling to rectify the time variable image backgrounds to common levels; and an adaptive image stretch to present images for visualization. The creation of the maps is supported by the development of four new tools that when fully tested will become part of the Montage distribution. The compute intensive part of the processing lies in the reprojection of the images, and we show how we optimized the processing for efficient creation of mosaics that are used in turn to create maps in the HiPS tiling scheme. We plan to apply our methodology to infrared image data sets such a those delivered by Spitzer, 2MASS, IRAS and Planck.


title: Faint Standards for ZYJHK from the UKIDSS and VISTA Surveys

abstract: The currently defined "UKIRT Faint Standards" have JHK magnitudes between 10 and 15, with K_median=11.2. These stars will be too bright for the next generation of large telescopes. We have used multi-epoch observations taken as part of the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Visible and Infrared Survey Telescope for Astronomy (VISTA) surveys to identify non-variable stars with JHK magnitudes in the range 16-19. The stars were selected from the UKIDSS Deep Extragalactic Survey (DXS) and Ultra Deep Survey (UDS), the WFCAM calibration data (WFCAMCAL08B), the VISTA Deep Extragalactic Observations (VIDEO) and UltraVISTA. Sources selected from the near-infrared databases were paired with the Pan-STARRS Data Release 2 of optical to near-infrared photometry and the Gaia astrometric Data Release 2. Colour indices and other measurements were used to exclude sources that did not appear to be simple single stars. From an initial selection of 169 sources, we present a final sample of 81 standard stars with ZYJHK magnitudes, or a subset, each with 20 to 600 observations in each filter. The new standards have Ks_median=17.5. The relative photometric uncertainty for the sample is <0.006 mag and the absolute uncertainty is estimated to be <~0.02 mag. The sources are distributed equatorially and are accessible from both hemispheres.


title: BUDD: Multi-modal Bayesian Updating Deforestation Detections

abstract: The global phenomenon of forest degradation is a pressing issue with severe implications for climate stability and biodiversity protection. In this work we generate Bayesian updating deforestation detection (BUDD) algorithms by incorporating Sentinel-1 backscatter and interferometric coherence with Sentinel-2 normalized vegetation index data. We show that the algorithm provides good performance in validation AOIs. We compare the effectiveness of different combinations of the three data modalities as inputs into the BUDD algorithm and compare against existing benchmarks based on optical imagery.


title: Abnormally low thermal conductivity of 2D selenene: An ab initio study

abstract: The lattice thermal conductivity and thermal transport properties of 2D $\alpha$-selenene are investigated based on the first-principles calculations. The isotropic in-plane thermal conductivity is as low as 3.04 W m$^{-1}$ K$^{-1}$ at room temperature, even abnormally lower than $\alpha$-tellurene which processes analogous configuration and lower Debye temperature. We find this abnormal phenomenon reasonably stems from the larger anharmonicity of the acoustic phonon branch. Moreover, the phonon spectra, elastic properties, and related thermal properties are also exhibited. Acoustic phonons contribute mainly to the total thermal conductivity. Furthermore, size effect, boundary effect, the total phase space for three-phonon processes, phonon group velocity and relaxation time are further investigated, and the last one is unveiled to be the key ingredient of thermal transport in 2D selenene.


title: High Speed Mid-Wave Infrared Uni-traveling Carrier Photodetector

abstract: Mid-wave infrared (MWIR) frequency comb is expected to dramatically improve the precision and sensitivity of molecular spectroscopy. For high resolution application, high speed MWIR photodetector is one of the key components, however, the commercially available high speed MWIR photodetector only has sub-GHz bandwidth currently. In this paper, we demonstrate, for the first time to our knowledge, a high speed mid-wave infrared (MWIR) uni-traveling carrier photodetector based on InAs/GaSb type-II superlattice (T2SL) at room temperature. The device exhibits a cutoff wavelength of 5.6{\mu}m, and 3dB bandwidth of 6.58 GHz for a 20{\mu}m diameter device at 300K. These promising results show the device has potential to be utilized in high speed applications such as frequency comb spectroscopy, free space communication and others. The limitations on the high frequency performance of the photodetectors are also discussed.


title: Quantifying Observed Prior Impact

abstract: We distinguish two questions (i) how much information does the prior contain? and (ii) what is the effect of the prior? Several measures have been proposed for quantifying effective prior sample size, for example Clarke [1996] and Morita et al. [2008]. However, these measures typically ignore the likelihood for the inference currently at hand, and therefore address (i) rather than (ii). Since in practice (ii) is of great concern, Reimherr et al. [2014] introduced a new class of effective prior sample size measures based on prior-likelihood discordance. We take this idea further towards its natural Bayesian conclusion by proposing measures of effective prior sample size that not only incorporate the general mathematical form of the likelihood but also the specific data at hand. Thus, our measures do not average across datasets from the working model, but condition on the current observed data. Consequently, our measures can be highly variable, but we demonstrate that this is because the impact of a prior can be highly variable. Our measures are Bayes estimates of meaningful quantities and well communicate the extent to which inference is determined by the prior, or framed differently, the amount of effort saved due to having prior information. We illustrate our ideas through a number of examples including a Gaussian conjugate model (continuous observations), a Beta-Binomial model (discrete observations), and a linear regression model (two unknown parameters). Future work on further developments of the methodology and an application to astronomy are discussed at the end.


title: Global simulations of stellar dynamos

abstract: The dynamo mechanism, responsible for the solar magnetic activity, is still an open problem in astrophysics. Different theories proposed to explain such phenomena have failed in reproducing the observational properties of the solar magnetism. Thus, ab-initio computational modeling of the convective dynamo in a spherical shell turns out as the best alternative to tackle this problem. In this work we review the efforts performed in global simulations over the past decades. Regarding the development and sustain of mean-flows, as well as mean magnetic field, we discuss the points of agreement and divergence between the different modeling strategies. Special attention is given to the implicit large-eddy simulations performed with the EULAG-MHD code.


title: Detecting Light Boson Dark Matter through Conversion into Magnon

abstract: Light boson dark matter such as axion or hidden photon can be resonantly converted into a magnon in a magnetic insulator under the magnetic field, which can be detected experimentally. We provide a quantum mechanical formulation for the magnon event rate and show that the result is consistent with that obtained by a classical calculation. Besides, it is pointed out that the experimental setup of the QUAX proposal for the axion detection also works as a detector of hidden photon dark matter. It has good sensitivity in the mass range around 1 meV, which is beyond astrophysical constraints.


title: Interpretable Rumor Detection in Microblogs by Attending to User Interactions

abstract: We address rumor detection by learning to differentiate between the community's response to real and fake claims in microblogs. Existing state-of-the-art models are based on tree models that model conversational trees. However, in social media, a user posting a reply might be replying to the entire thread rather than to a specific user. We propose a post-level attention model (PLAN) to model long distance interactions between tweets with the multi-head attention mechanism in a transformer network. We investigated variants of this model: (1) a structure aware self-attention model (StA-PLAN) that incorporates tree structure information in the transformer network, and (2) a hierarchical token and post-level attention model (StA-HiTPLAN) that learns a sentence representation with token-level self-attention. To the best of our knowledge, we are the first to evaluate our models on two rumor detection data sets: the PHEME data set as well as the Twitter15 and Twitter16 data sets. We show that our best models outperform current state-of-the-art models for both data sets. Moreover, the attention mechanism allows us to explain rumor detection predictions at both token-level and post-level.


title: The impact of isomers on a kilonova associated with neutron star mergers

abstract: We investigate the significance of isomers on a kilonova associated with neutron star mergers (NSMs) for the first time. We calculate the evolution of abundances and energy generation rates ($\dot{\epsilon}_{\rm nuc}$) of ejecta from NSMs, taking into account $\beta^-$ decay through isomers. We find that for ejecta with electron fraction ($Y_e$) of $0.2-0.3$, $\dot{\epsilon}_{\rm nuc}$ is appreciably changed from those without isomers, due to the large change in timing of $\beta^-$ decay through the isomeric states of ${}^{123, 125, 127} \rm Sn$ and ${}^{128} \rm Sb$. In particular, the effects of the isomers on $\dot{\epsilon}_{\rm nuc}$ are prominent for ejecta of $Y_e \sim 0.25$, which could emit a fraction of an early, blue component of a kilonova observed in GW170817. When the excitation from a ground state to its isomeric state is important, the isomers of ${}^{129} \rm Sb$ and ${}^{129, 131} \rm Te$ also cause additional and appreciable change in $\dot{\epsilon}_{\rm nuc}$. Furthermore, we demonstrate that larger amounts of lanthanide-free ejecta result in a better fit of an observed light curve of the kilonova in GW170817 if the isomers are taken into account.


title: A Stochastic Subgradient Method for Nonsmooth Nonconvex Multi-Level Composition Optimization

abstract: We propose a single time-scale stochastic subgradient method for constrained optimization of a composition of several nonsmooth and nonconvex functions. The functions are assumed to be locally Lipschitz and differentiable in a generalized sense. Only stochastic estimates of the values and generalized derivatives of the functions are used. The method is parameter-free. We prove convergence with probability one of the method, by associating with it a system of differential inclusions and devising a nondifferentiable Lyapunov function for this system. For problems with functions having Lipschitz continuous derivatives, the method finds a point satisfying an optimality measure with error of order $1/\sqrt{N}$, after executing $N$ iterations with constant stepsize.


title: Constructing a subgradient from directional derivatives for functions of two variables

abstract: For any scalar-valued bivariate function that is locally Lipschitz continuous and directionally differentiable, it is shown that a subgradient may always be constructed from the function's directional derivatives in the four compass directions, arranged in a so-called "compass difference". When the original function is nonconvex, the obtained subgradient is an element of Clarke's generalized gradient, but the result appears to be novel even for convex functions. The function is not required to be represented in any particular form, and no further assumptions are required, though the result is strengthened when the function is additionally L-smooth in the sense of Nesterov. For certain optimal-value functions and certain parametric solutions of differential equation systems, these new results appear to provide the only known way to compute a subgradient. These results also imply that centered finite differences will converge to a subgradient for bivariate nonsmooth functions. As a dual result, we find that any compact convex set in two dimensions contains the midpoint of its interval hull. Examples are included for illustration, and it is demonstrated that these results do not extend directly to functions of more than two variables or sets in higher dimensions.


title: Second order subexponentiality and infinite divisibility

abstract: We characterize the second order subexponentiality of an infinitely divisible distribution on the real line under an exponential moment assumption. We investigate the asymptotic behaviour of the difference between the tails of an infinitely divisible distribution and its L\'evy measure. Moreover, we study the second order asymptotic behaviour of the tail of the $t$-th convolution power of an infinitely divisible distribution. The density version for a self-decomposable distribution on the real line without an exponential moment assumption is also given. Finally, the regularly varying case for a self-decomposable distribution on the half line is discussed.


title: Clues to the origin and properties of magnetic white dwarfs

abstract: A significant fraction of white dwarfs possess a magnetic field with strengths ranging from a few kG up to about 1000 MG. However, the incidence of magnetism varies when the white dwarf population is broken down into different spectral types providing clues on the formation of magnetic fields in white dwarfs. Several scenarios for the origin of magnetic fields have been proposed from a fossil field origin to dynamo generation at various stages of evolution. Offset dipoles are often assumed sufficient to model the field structure, however time-resolved spectropolarimetric observations have revealed more complex structures such as magnetic spots or multipoles. Surface mapping of these field structures combined with measured rotation rates help distinguish scenarios involving single star evolution from other scenarios involving binary interactions. I describe key observational properties of magnetic white dwarfs such as age, mass, and field strength, and confront proposed formation scenarios with these properties.


title: Assistive Relative Pose Estimation for On-orbit Assembly using Convolutional Neural Networks

abstract: Accurate real-time pose estimation of spacecraft or object in space is a key capability necessary for on-orbit spacecraft servicing and assembly tasks. Pose estimation of objects in space is more challenging than for objects on Earth due to space images containing widely varying illumination conditions, high contrast, and poor resolution in addition to power and mass constraints. In this paper, a convolutional neural network is leveraged to uniquely determine the translation and rotation of an object of interest relative to the camera. The main idea of using CNN model is to assist object tracker used in on space assembly tasks where only feature based method is always not sufficient. The simulation framework designed for assembly task is used to generate dataset for training the modified CNN models and, then results of different models are compared with measure of how accurately models are predicting the pose. Unlike many current approaches for spacecraft or object in space pose estimation, the model does not rely on hand-crafted object-specific features which makes this model more robust and easier to apply to other types of spacecraft. It is shown that the model performs comparable to the current feature-selection methods and can therefore be used in conjunction with them to provide more reliable estimates.


title: Even cycles and perfect matchings in claw-free plane graphs

abstract: Lov{\'a}sz showed that a matching covered graph $G$ has an ear decomposition starting with an arbitrary edge of $G$. Let $G$ be a graph which has a perfect matching. We call $G$ cycle-nice if for each even cycle $C$ of $G$, $G-V(C)$ has a perfect matching. If $G$ is a cycle-nice matching covered graph, then $G$ has ear decompositions starting with an arbitrary even cycle of $G$. In this paper, we characterize cycle-nice claw-free plane graphs. We show that the only cycle-nice simple 3-connected claw-free plane graphs are $K_4$, $W_5$ and $\overline C_6$. Furthermore, every cycle-nice 2-connected claw-free plane graph can be obtained from a graph in the family ${\cal F}$ by a sequence of three types of operations, where ${\cal F}$ consists of even cycles, a diamond, $K_4$, and $\overline C_6$.


title: Machine Learning in Thermodynamics: Prediction of Activity Coefficients by Matrix Completion

abstract: Activity coefficients, which are a measure of the non-ideality of liquid mixtures, are a key property in chemical engineering with relevance to modeling chemical and phase equilibria as well as transport processes. Although experimental data on thousands of binary mixtures are available, prediction methods are needed to calculate the activity coefficients in many relevant mixtures that have not been explored to-date. In this report, we propose a probabilistic matrix factorization model for predicting the activity coefficients in arbitrary binary mixtures. Although no physical descriptors for the considered components were used, our method outperforms the state-of-the-art method that has been refined over three decades while requiring much less training effort. This opens perspectives to novel methods for predicting physico-chemical properties of binary mixtures with the potential to revolutionize modeling and simulation in chemical engineering.


title: Low-rank quaternion tensor completion for recovering color videos and images

abstract: Low-rank quaternion tensor completion method, a novel approach to recovery color videos and images is proposed in this paper. We respectively reconstruct a color image and a color video as a quaternion matrix (second-order tensor) and a third-order quaternion tensor by encoding the red, green, and blue channel pixel values on the three imaginary parts of a quaternion. Different from some traditional models which treat color pixel as a scalar and represent color channels separately, whereas, during the quaternion-based reconstruction, it is significant that the inherent color structures of color images and color videos can be completely preserved. Under the definition of Tucker rank, the global low-rank prior to quaternion tensor is encoded as the nuclear norm of unfolding quaternion matrices. Then, by applying the ADMM framework, we provide the tensor completion algorithm for any order quaternion tensors, which theoretically can be well used to recover missing entries of any multidimensional data with color structures. Simulation results for color videos and color images recovery show the superior performance and efficiency of the proposed method over some state-of-the-art existing ones.


title: Quaternion matrix regression for color face recognition

abstract: Regression analysis-based approaches have been widely studied for face recognition (FR) in the past several years. More recently, to better deal with some difficult conditions such as occlusions and illumination, nuclear norm based matrix regression methods have been proposed to characterize the low-rank structure of the error image, which generalize the one-dimensional, pixel-based error model to the two-dimensional structure. These methods, however, are inherently devised for grayscale image based FR and without exploiting the color information which is proved beneficial for FR of color face images. Benefiting from quaternion representation, which is capable of encoding the cross-channel correlation of color images, we propose a novel color FR method by formulating the color FR problem as a nuclear norm based quaternion matrix regression (NQMR). We further develop a more robust model called R-NQMR by using the logarithm of the nuclear norm, instead of the original nuclear norm, which adaptively assigns weights on different singular values, and then extend it to deal with the mixed noise. The proposed models, then, are solved using the effective alternating direction multiplier method (ADMM). Experiments on several public face databases demonstrate the superior performance and efficacy of the proposed approaches for color FR, especially for some difficult conditions (occlusion, illumination and mixed noise) over some state-of-the-art regression analysis-based approaches.


title: An Improved Quadrature Voltage-Controlled Oscillator with Through-Silicon-Via Inductor in Three-dimensional Integrated Circuits

abstract: Low-power quadrature voltage-controlled oscillator (QVCO) design utilizing transformer-feedback and current-reuse techniques with increased frequency range is proposed in this paper. With increasing demand for QVCOs in on-chip applications, the conventional spiral inductor based approaches for QVCOs has become a major bottleneck due to their large size. To address this concern, we propose to replace the conventional spiral inductor based approaches with through-silicon-via (TSV) inductor based approach in three-dimensional integrated circuits (3D ICs). In addition, the proposed QVCO circuit can provide higher frequency range of operation compared with conventional designs. Experimental results show by replacing conventional spiral transformers with TSV transformers, up to 3.9x reduction in metal resource consumption. The proposed QVCOs achieves a phase noise of -114 $dBc/Hz$@1 $MHz$ and -111.2 $dBc/Hz$@1 $MHz$ at the carrier of 2.5 $GHz$ for toroidal TSV transformed based-QVCO and vertical spiral transformer based-QVCO respectively. The power consumption is only 1.5 $mW$ and 1.7 $mW$ for toroidal TSV transformed based-QVCO and vertical spiral transformer based-QVCO respectively, under the supply voltage of 0.7 $V$.


title: Adaptive Estimation and Statistical Inference for High-Dimensional Graph-Based Linear Models

abstract: We consider adaptive estimation and statistical inference for high-dimensional graph-based linear models. In our model, the coordinates of regression coefficients correspond to an underlying undirected graph. Furthermore, the given graph governs the piecewise polynomial structure of the regression vector. In the adaptive estimation part, we apply graph-based regularization techniques and propose a family of locally adaptive estimators called the Graph-Piecewise-Polynomial-Lasso. We further study a one-step update of the Graph-Piecewise-Polynomial-Lasso for the problem of statistical inference. We develop the corresponding theory, which includes the fixed design and the sub-Gaussian random design. Finally, we illustrate the superior performance of our approaches by extensive simulation studies and conclude with an application to an Arabidopsis thaliana microarray dataset.


title: Widely tunable cavity-enhanced frequency combs

abstract: We describe the cavity-enhancement of frequency combs over a wide tuning range of 450-700 nm (> 7900 cm$^{-1}$), covering nearly the entire visible spectrum. Tunable visible frequency combs from a synchronously-pumped optical parametric oscillator are coupled into a 4-mirror, dispersion-managed cavity with a finesse of 600 to 1400. An intracavity absorption path length enhancement greater than 190 is obtained over the entire tuning range, while preserving intracavity spectral bandwidths capable of supporting sub-200 fs pulse durations. These tunable cavity-enhanced frequency combs can find many applications in nonlinear optics and spectroscopy.


title: Kalibre: Knowledge-based Neural Surrogate Model Calibration for Data Center Digital Twins

abstract: Computational fluid dynamics (CFD) model has been widely used for prototyping data centers. Evolving it to high-fidelity {\em digital twin} is desirable for the management and operations of large-scale data centers. Manually calibrating CFD model parameters to achieve twin-class fidelity by specially trained domain expert is tedious and labor-intensive. To reduce manual efforts, existing automatic calibration approaches developed for various computational models apply heuristics to search model configurations within an empirically defined parameter bound. However, in the context of CFD, each search step requires long-lasting CFD model's iterated solving, rendering these approaches impractical with increased model complexity. This paper presents Kalibre, a knowledge-based neural surrogate approach that performs CFD model calibration by iterating four key steps of i) training a neural surrogate model based on CFD-generated data, ii) finding the optimal parameters at the moment through neural surrogate retraining based on sensor-measured data, iii) configuring the found parameters back to the CFD model, and iv) validating the CFD model using sensor-measured data as the ground truth. Thus, the parameter search is offloaded to the neural surrogate which is ultra-faster than CFD model's iterated solving. To speed up the convergence of Kalibre, we integrate prior knowledge of the twinned data center's thermophysics into the neural surrogate design to improve its learning efficiency. With about five hours computation on a 32-core processor, Kalibre achieves mean absolute errors (MAEs) of $0.81^o$C and $0.75^o$C in calibrating two CFD models for two production data halls hosting thousands of servers each while requires fewer CFD solving processes than existing baseline approaches.


title: Large time asymptotics for a cubic nonlinear Schr\"odinger system in one space dimension, II

abstract: This is a sequel to the paper "Large time asymptotics for a cubic nonlinear Schr\"odinger system in one space dimension" by the same authors. We continue to study the Cauchy problem for the two-component system of cubic nonlinear Schr\"odinger equations in one space dimension. We provide criteria for large time decay or non-decay in $L^2$ of the small amplitude solutions in terms of the Fourier transforms of the initial data.


title: On perturbative constraints for vacuum f(R) gravity

abstract: Perturbative techniques are important for modified theories of gravity since they allow to calculate deviations from General Relativity without recurring to exact solutions, which can be difficult to find. When applied to models such as $f(R)$ gravity, these techniques introduce corrections in the field equations that involve higher order derivatives. Such corrections must be handled carefully to have a well defined perturbative scheme, and this can be achieved through the method of perturbative constraints, where the coefficient of the additional term in the action is used as expansion parameter for the quantities of interest. In this work, we implement a perturbative framework that compares solutions in modified theories of gravity with solutions of the Einstein field equations, by following the guidelines of perturbation theory constructed in General Relativity together with the perturbative constraints rationale. By using this formalism, we demonstrate that a consistent $f(R)$ perturbation theory in vacuum, for an important class of $f(R)$ functions, produces no additional effects with respect to what is expected from the perturbation theory of General Relativity. From this result, we argue that there are fundamental limitations that explain why the solutions of some $f(R)$ models can be disconnected from their general relativistic counterparts, in the sense that the limit that leads from the $f(R)$ action to General Relativity does not transform the solutions accordingly.


title: Theoretical prediction of a low-energy Stone-Wales graphene with intrinsic type-III Dirac-cone

abstract: Based on first-principles method we predict a new low-energy Stone-Wales graphene SW40, which has an orthorhombic lattice with Pbam symmetry and 40 carbon atoms in its crystalline cell forming well-arranged Stone-Wales patterns. The calculated total energy of SW40 is just about 133 meV higher than that of graphene, indicating its excellent stability exceeds all the previously proposed graphene allotropes. We find that SW40 processes intrinsic Type-III Dirac-cone (Phys. Rev. Lett., 120, 237403, 2018) formed by band-crossing of a local linear-band and a local flat-band, which can result in highly anisotropic Fermions in the system. Interestingly, such intrinsic type-III Dirac-cone can be effectively tuned by inner-layer strains and it will be transferred into Type-II and Type-I Dirac-cones under tensile and compressed strains, respectively. Finally, a general tight-binding model was constructed to understand the electronic properties nearby the Fermi-level in SW40. The results show that type-III Dirac-cone feature can be well understood by the $\pi$-electron interactions between adjacent Stone-Wales defects.


title: PulseSatellite: A tool using human-AI feedback loops for satellite image analysis in humanitarian contexts

abstract: Humanitarian response to natural disasters and conflicts can be assisted by satellite image analysis. In a humanitarian context, very specific satellite image analysis tasks must be done accurately and in a timely manner to provide operational support. We present PulseSatellite, a collaborative satellite image analysis tool which leverages neural network models that can be retrained on-the fly and adapted to specific humanitarian contexts and geographies. We present two case studies, in mapping shelters and floods respectively, that illustrate the capabilities of PulseSatellite.


title: A Chandra X-ray Survey of Optically Selected AGN Pairs

abstract: We present a Chandra archival study of optically selected AGN pairs at a median redshift $\bar{z} \sim$ 0.1. Out of 1286 AGN pairs (with projected separations $r_p < 100$ kpc and velocity offsets $\Delta v < 600$ km s$^{-1}$) optically identified from the Sloan Digital Sky Survey Seventh Data Release, we find 67 systems with archival Chandra observations, which represent the largest sample of optically selected AGN pairs studied in the X-ray. Among the 67 AGN pairs, 21 systems have both nuclei detected in the X-rays, 36 have one nucleus detected in the X-rays, whereas 10 have no X-ray detection. The X-ray detection rate, 78/134 = 58% ($\pm 7$% 1$\sigma$ Poisson errors), is significantly higher than that (23/134 = $17\% \pm 4\%$) of a comparison sample of star-forming galaxy pairs, lending support to the optical AGN classification. In the conservative case where X-ray contamination from star formation is removed, the X-ray detection rate becomes $27\% \pm 4\%$, consistent with predictions from latest galaxy merger simulations. The 2-10 keV X-ray luminosity $L_{2-10~{\rm keV}}$ increases with decreasing projected separation in AGN pairs for $r_p \gtrsim 15$ kpc, suggesting an enhancement of black hole accretion even in early-stage mergers. On the other hand, $L_{2-10~{\rm keV}}$ appears to decrease with decreasing projected separation at $r_p \lesssim 15$ kpc, which is contradictory to predictions from merger simulations. The apparent decrease in $L_{2-10~{\rm keV}}$ of AGN pairs at $r_p \lesssim 15$ kpc may be caused by: (i) enhanced absorbing columns from merger-induced gas inflows, (ii) feedback effects from early-stage mergers, and/or (iii) small number statistics. Future X-ray studies with larger samples are needed to put our results on firmer statistical ground.


title: A regularity theory for stochastic partial differential equations with a super-linear diffusion coefficient and a spatially homogeneous colored noise

abstract: Existence, uniqueness, and regularity of a strong solution are obtained for stochastic PDEs with a colored noise $F$ and its super-linear diffusion coefficient: $$ du=(a^{ij}u_{x^ix^j}+b^iu_{x^i}+cu)dt+\xi|u|^{1+\lambda}dF, \quad (t,x)\in(0,\infty)\times\mathbb{R}^d, $$ where $\lambda \geq 0$ and the coefficients depend on $(\omega,t,x)$. The strategy of handling nonlinearity of the diffusion coefficient is to find a sharp estimation for a general Lipschitz case, and apply it to the super-linear case. Moreover, investigation for the estimate provides a range of $\lambda$, a sufficient condition for the unique solvability, where the range depends on the spatial covariance of $F$ and the spatial dimension $d$.


title: Rough differential equations with path-dependent coefficients

abstract: We establish the existence of solutions to path-dependent rough differential equations with non-anticipative coefficients. Regularity assumptions on the coefficients are formulated in terms of horizontal and vertical derivatives.


title: Number of constituent quark scaling of elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV

abstract: We briefly summarize our recent study on the number of constituent quark (NCQ) scaling of hadron elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV. With the inclusion of hadron production via the quark coalescence model at intermediate $p_T$, the viscous hydrodynamics at low $p_T$, and jet fragmentation at high $p_T$, our $Hydro-Coal-Frag$ model provides a nice description of the $p_T$-spectra and differential elliptic flow $v_2(p_T)$ of pions, kaons and protons over the $p_T$ range from 0 to 6 GeV. Our results demonstrate that including the quark coalescence is essential for reproducing the observed approximate NCQ scaling of hadron $v_2$ at intermediate $p_T$ in experiments, indicating strongly the existence of partonic degrees of freedom and the formation of quark-gluon plasma in high multiplicity p--Pb collisions at the LHC.


title: Merger rate density of binary black holes formed in open clusters

abstract: Several binary black holes (BBHs) have been observed using gravitational wave detectors. For the formation mechanism of BBHs, two main mechanisms, isolated binary evolution and dynamical formation in dense star clusters, have been suggested. Future observations are expected to provide more information about BBH distributions, and it will help us to distinguish the two formation mechanisms. For the star cluster channel, globular clusters have mainly been investigated. However, recent simulations have suggested that BBH formation in open clusters is not negligible. We estimate a local merger rate density of BBHs originated from open clusters using the results of our $N$-body simulations of open clusters with four different metallicities. We find that the merger rate per cluster is the highest for our 0.1 solar metallicity model. Assuming a cosmic star formation history and a metallicity evolution with dispersion, we estimate the local merger rate density of BBHs originated from open clusters to be $\sim 70~{\rm yr}^{-1} {\rm Gpc}^{-3}$. This value is comparable to the merger rate density expected from the first and second observation runs of LIGO and Virgo. In addition, we find that BBH mergers obtained from our simulations can reproduce the distribution of primary mass and mass ratio of merging BBHs estimated from the LIGO and Virgo observations.


title: The $4 \times 4$ orthostochastic variety

abstract: Orthostochastic matrices are the entrywise squares of orthogonal matrices, and naturally arise in various contexts, including notably definite symmetric determinantal representations of real polynomials. However, defining equations for the real variety were previously known only for $3 \times 3$ matrices. We study the real variety of $4 \times 4$ orthostochastic matrices, and find a minimal defining set of equations consisting of 6 quintics and 3 octics. The techniques used here involve a wide range of both symbolic and computational methods, in computer algebra and numerical algebraic geometry.


title: ImVoteNet: Boosting 3D Object Detection in Point Clouds with Image Votes

abstract: 3D object detection has seen quick progress thanks to advances in deep learning on point clouds. A few recent works have even shown state-of-the-art performance with just point clouds input (e.g. VoteNet). However, point cloud data have inherent limitations. They are sparse, lack color information and often suffer from sensor noise. Images, on the other hand, have high resolution and rich texture. Thus they can complement the 3D geometry provided by point clouds. Yet how to effectively use image information to assist point cloud based detection is still an open question. In this work, we build on top of VoteNet and propose a 3D detection architecture called ImVoteNet specialized for RGB-D scenes. ImVoteNet is based on fusing 2D votes in images and 3D votes in point clouds. Compared to prior work on multi-modal detection, we explicitly extract both geometric and semantic features from the 2D images. We leverage camera parameters to lift these features to 3D. To improve the synergy of 2D-3D feature fusion, we also propose a multi-tower training scheme. We validate our model on the challenging SUN RGB-D dataset, advancing state-of-the-art results by 5.7 mAP. We also provide rich ablation studies to analyze the contribution of each design choice.


title: FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45 m telescope (FUGIN). VI. Dense gas and mini-starbursts in the W43 giant molecular cloud complex

abstract: We performed new large-scale $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=$1--0 observations of the W43 giant molecular cloud complex in the tangential direction of the Scutum arm ($l\sim {30^\circ}$) as a part of the FUGIN project. The low-density gas traced by $^{12}$CO is distributed over 150 pc $\times$ 100 pc ($l \times b$), and has a large velocity dispersion (20-30 km s$^{-1}$). However, the dense gas traced by C$^{18}$O is localized in the W43 Main, G30.5, and W43 South (G29.96-0.02) high-mass star-forming regions in the W43 GMC complex, which have clumpy structures. We found at least two clouds with a velocity difference of $\sim$ 10-20 km s$^{-1}$, both of which are likely to be physically associated with these high-mass star-forming regions based on the results of high $^{13}$CO $J=$ 3-2 to $J =$ 1-0 intensity ratio and morphological correspondence with the infrared dust emission. The velocity separation of these clouds in W43 Main, G30.5, and W43 South is too large for each cloud to be gravitationally bound. We also revealed that the dense gas in the W43 GMC has a high local column density, while "the current SFE" of entire the GMC is low ($\sim 4\%$) compared with the W51 and M17 GMC. We argue that the supersonic cloud-cloud collision hypothesis can explain the origin of the local mini-starbursts and dense gas formation in the W43 GMC complex.


title: Finite and symmetric Mordell-Tornheim multiple zeta values

abstract: We introduce finite and symmetric Mordell-Tornheim type of multiple zeta values and give a new approach to the Kaneko-Zagier conjecture stating that the finite and symmetric multiple zeta values satisfy the same relations.


title: Early-detection and classification of live bacteria using time-lapse coherent imaging and deep learning

abstract: We present a computational live bacteria detection system that periodically captures coherent microscopy images of bacterial growth inside a 60 mm diameter agar-plate and analyzes these time-lapsed holograms using deep neural networks for rapid detection of bacterial growth and classification of the corresponding species. The performance of our system was demonstrated by rapid detection of Escherichia coli and total coliform bacteria (i.e., Klebsiella aerogenes and Klebsiella pneumoniae subsp. pneumoniae) in water samples. These results were confirmed against gold-standard culture-based results, shortening the detection time of bacterial growth by >12 h as compared to the Environmental Protection Agency (EPA)-approved analytical methods. Our experiments further confirmed that this method successfully detects 90% of bacterial colonies within 7-10 h (and >95% within 12 h) with a precision of 99.2-100%, and correctly identifies their species in 7.6-12 h with 80% accuracy. Using pre-incubation of samples in growth media, our system achieved a limit of detection (LOD) of ~1 colony forming unit (CFU)/L within 9 h of total test time. This computational bacteria detection and classification platform is highly cost-effective (~$0.6 per test) and high-throughput with a scanning speed of 24 cm2/min over the entire plate surface, making it highly suitable for integration with the existing analytical methods currently used for bacteria detection on agar plates. Powered by deep learning, this automated and cost-effective live bacteria detection platform can be transformative for a wide range of applications in microbiology by significantly reducing the detection time, also automating the identification of colonies, without labeling or the need for an expert.


title: Spiking Inception Module for Multi-layer Unsupervised Spiking Neural Networks

abstract: Spiking Neural Network (SNN), as a brain-inspired approach, is attracting attention due to its potential to produce ultra-high-energy-efficient hardware. Competitive learning based on Spike-Timing-Dependent Plasticity (STDP) is a popular method to train an unsupervised SNN. However, previous unsupervised SNNs trained through this method are limited to a shallow network with only one learnable layer and cannot achieve satisfactory results when compared with multi-layer SNNs. In this paper, we eased this limitation by: 1)We proposed a Spiking Inception (Sp-Inception) module, inspired by the Inception module in the Artificial Neural Network (ANN) literature. This module is trained through STDP-based competitive learning and outperforms the baseline modules on learning capability, learning efficiency, and robustness. 2)We proposed a Pooling-Reshape-Activate (PRA) layer to make the Sp-Inception module stackable. 3)We stacked multiple Sp-Inception modules to construct multi-layer SNNs. Our algorithm outperforms the baseline algorithms on the hand-written digit classification task, and reaches state-of-the-art results on the MNIST dataset among the existing unsupervised SNNs.


title: Estimating the Jet Power of Mrk\,231 During the 2017-2018 Flare

abstract: Long-term 17.6~GHz radio monitoring of the broad absorption line quasar, Mrk\,231, detected a strong flare in late 2017. This triggered four epochs of Very Long Baseline Array (VLBA) observations from 8.4~GHz to 43~GHz over a 10-week period as well as an X-ray observation with NuSTAR. This was the third campaign of VLBA monitoring that we have obtained. The 43~GHz VLBA was degraded in all epochs with only 7 of 10 antennas available in three epochs and 8 in the first epoch. However, useful results were obtained due to a fortuitous capturing of a complete short 100~mJy flare at 17.6~GHz: growth and decay. This provided useful constraints on the physical model of the ejected plasma that were not available in previous campaigns. We consider four classes of models, discrete ejections (both protonic and positronic) and jetted (protonic and positronic). The most viable model is a "dissipative bright knot" in a faint background leptonic jet with an energy flux $\sim10^{43}$ ergs/sec. Inverse Compton scattering calculations (based on these models) in the ambient quasar photon field explains the lack of a detectable increase in X-ray luminosity measured by NuSTAR. We show that the core (the bright knot) moves towards a nearby secondary at $\approx 0.97$c. The background jet is much fainter. Evidently, the high frequency VLBA core does not represent the point of origin of blazar jets, in general, and optical depth "core shift" estimates of jet points of origin can be misleading.


title: Continuous gravitational wave from magnetized white dwarfs

abstract: Recent evidence of super-Chandrasekhar white dwarfs (WDs), from the observations of over-luminous type Ia supernovae (SNeIa), has been a great astrophysical discovery. However, no such massive WDs have so far been observed directly as their luminosities are generally quite low. Hence it immediately raises the question of whether there is any possibility of detecting them directly. The search for super-Chandrasekhar WDs is very important as SNeIa are used as standard candles in cosmology. In this article, we show that continuous gravitational wave can allow us to detect such super-Chandrasekhar WDs directly.


title: Circle actions on 8-dimensional almost complex manifolds with 4 fixed points

abstract: Consider a circle action on an 8-dimensional compact almost complex manifold with 4 fixed points. To the author's knowledge, $S^2 \times S^6$ is the only known example of such a manifold. In this paper, we prove that if the circle acts on an 8-dimensional compact almost complex manifold $M$ with 4 fixed points, all the Chern numbers and the Hirzebruch $\chi_y$-genus of $M$ agree with those of $S^2 \times S^6$. In particular, $M$ is unitary cobordant to $S^2 \times S^6$.


title: On the Steady Prandtl Boundary Layer Expansions

abstract: In this paper, we consider the zero-viscosity limit of the 2D steady Navier-Stokes equations in $(0,L)\times\mathbb{R}^+$ with non-slip boundary conditions. By estimating the stream-function of the remainder, we justify the validity of the Prandtl boundary layer expansions.


title: Neutron stars with a generalized Proca hair and spontaneous vectorization

abstract: In a class of generalized Proca theories, we study the existence of neutron star solutions with a nonvanishing temporal component of the vector field $A_\mu$ approaching 0 toward spatial infinity, as they may be the endpoints of tachyonic instabilities of neutron star solutions in general relativity with $A_{\mu}=0$. Such a phenomenon is called spontaneous vectorization, which is analogous to spontaneous scalarization in scalar-tensor theories with nonminimal couplings to the curvature or matter. For the nonminimal coupling $\beta X R$, where $\beta$ is a coupling constant and $X=-A_{\mu}A^{\mu}/2$, we show that there exist both 0-node and 1-node vector-field solutions, irrespective of the choice of the equations of state of nuclear matter. The 0-node solution, which is present only for $\beta=-{\cal O}(0.1)$, may be induced by some nonlinear effects such as the selected choice of initial conditions. The 1-node solution exists for $\beta=-{\cal O}(1)$, which suddenly emerges above a critical central density of star and approaches the general relativistic branch with the increasing central density. We compute the mass $M$ and radius $r_s$ of neutron stars for some realistic equations of state and show that the $M$-$r_s$ relations of 0-node and 1-node solutions exhibit notable difference from those of scalarized solutions in scalar-tensor theories. Finally, we discuss the possible endpoints of tachyonic instabilities.


title: Doubly regular Diophantine quadruples

abstract: For a nonzero integer n, a set of m distinct nonzero integers {a_1,a_2,...,a_m} such that a_i a_j + n is a perfect square for all 1 <= i < j <= m, is called a D(n)-m-tuple. In this paper, by using properties of so-called regular Diophantine m-tuples and certain family of elliptic curves, we show that there are infinitely many essentially different sets consisting of perfect squares which are simultaneously D(n_1)-quadruples and D(n_2)-quadruples with distinct non-zero squares n_1 and n_2.


title: Low Complexity Iterative Rake Detector for Orthogonal Time Frequency Space Modulation

abstract: This paper presents a linear complexity iterative rake detector for the recently proposed orthogonal time frequency space (OTFS) modulation scheme. The basic idea is to extract and combine the received multipath components of the transmitted symbols in the delay-Doppler grid using linear diversity combining schemes like maximal ratio combining (MRC), equal gain combining and selection combining to improve the SNR of the combined signal. We reformulate the OTFS input-output relation in the vector form by placing some null symbols in the delay-Doppler grid thereby exploiting the block circulant property of the channel matrix. Using the new input-output relation we propose a low complexity iterative detector based on the MRC scheme. The bit error rate (BER) performance of the proposed detector will be compared with the state of the art message passing detector and orthogonal frequency division multiplexing (OFDM) scheme employing a single tap minimum mean square error (MMSE) equalizer. We also show that the frame error rate (FER) performance of the MRC detector can be improved by employing error correcting codes operating in the form of a turbo decision feedback equalizer (DFE).


title: Matching numbers and dimension of edge ideals

abstract: Let $G$ be a finite simple graph on the vertex set $V(G) = \{x_{1}, \ldots, x_{n}\}$ and match$(G)$, min-match$(G)$ and ind-match$(G)$ the matching number, minimum matching number and induced matching number of $G$, respectively. Let $K[V(G)] = K[x_{1}, \ldots, x_{n}]$ denote the polynomial ring over a field $K$ and $I(G) \subset K[V(G)]$ the edge ideal of $G$. The relationship between these graph-theoretic invariants and ring-theoretic invariants of the quotient ring $K[V(G)]/I(G)$ has been studied. In the present paper, we study the relationship between match$(G)$, min-match$(G)$, ind-match$(G)$ and $\dim K[V(G)]/I(G)$.


title: Review of proton and nuclear shape fluctuations at high energy

abstract: Determining the inner structure of protons and nuclei in terms of their fundamental constituents has been one of the main tasks of high energy nuclear and particle physics experiments. This quest started as a mapping of the (average) parton densities as a function of longitudinal momentum fraction and resolution scale. Recently, the field has progressed to more differential imaging, where one important development is the description of the event-by-event quantum fluctuations in the wave function of the colliding hadron. In this Review, recent developments on the extraction of proton and nuclear transverse geometry with event-by-event fluctuations from collider experiments at high energy is presented. The importance of this fundamentally interesting physics in other collider experiments like in studies of the properties of the Quark Gluon Plasma is also illustrated.


title: Strengthened inequalities for the mean width and the $\ell$-norm

abstract: Barthe proved that the regular simplex maximizes the mean width of convex bodies whose John ellipsoid (maximal volume ellipsoid contained in the body) is the Euclidean unit ball; or equivalently, the regular simplex maximizes the $\ell$-norm of convex bodies whose L\"owner ellipsoid (minimal volume ellipsoid containing the body) is the Euclidean unit ball. Schmuckenschl\"ager verified the reverse statement; namely, the regular simplex minimizes the mean width of convex bodies whose L\"owner ellipsoid is the Euclidean unit ball. Here we prove stronger stability versions of these results. We also consider related stability results for the mean width and the $\ell$-norm of the convex hull of the support of centered isotropic measures on the unit sphere.


title: Inequality for local energy of Ising models with quenched randomness and its application

abstract: In this study, we extend the lower bound on the average of the local energy of the Ising model with quenched randomness [J. Phys. Soc. Jpn. 76, 074711 (2007)] obtained for a symmetric distribution to an asymmetric one. Compared with the case of symmetric distribution, our bound has a non-trivial term. By applying the acquired bound to a Gaussian distribution, we obtain the lower bounds on the expectation of the square of the correlation function. Thus, we demonstrate that in the Ising model in a Gaussian random field, the spin-glass order parameter generally has a finite value at any temperature, regardless of the forms of the other interactions.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: Superconductivity and charge density wave under a time-dependent periodic field in the one-dimensional attractive Hubbard model

abstract: We investigate the competition between superconductivity (SC) and charge density wave (CDW) under a time-dependent periodic field in the attractive Hubbard model. By employing the time-dependent exact diagonalization method, we show that the driving frequency and amplitude of the external field can control the enhancement of either the superconducting pair or the CDW correlations, which are degenerate in the ground state of the half-filled attractive Hubbard model in the absence of the field. In the strong-coupling limit of the attractive Hubbard interaction, the controllability is characterized by the anisotropic interaction of the effective model. The anisotropy is induced by the external field and lifts the degeneracy of SC and CDW. We find that the enhancement or suppression of the superconducting pair and CDW correlations in the periodically-driven attractive Hubbard model can be well interpreted by the quench dynamics of the effective model derived in the strong-coupling limit.


title: Depth Based Semantic Scene Completion with Position Importance Aware Loss

abstract: Semantic Scene Completion (SSC) refers to the task of inferring the 3D semantic segmentation of a scene while simultaneously completing the 3D shapes. We propose PALNet, a novel hybrid network for SSC based on single depth. PALNet utilizes a two-stream network to extract both 2D and 3D features from multi-stages using fine-grained depth information to efficiently captures the context, as well as the geometric cues of the scene. Current methods for SSC treat all parts of the scene equally causing unnecessary attention to the interior of objects. To address this problem, we propose Position Aware Loss(PA-Loss) which is position importance aware while training the network. Specifically, PA-Loss considers Local Geometric Anisotropy to determine the importance of different positions within the scene. It is beneficial for recovering key details like the boundaries of objects and the corners of the scene. Comprehensive experiments on two benchmark datasets demonstrate the effectiveness of the proposed method and its superior performance. Models and Video demo can be found at: https://github.com/UniLauX/PALNet.


title: Pre-defined Sparsity for Low-Complexity Convolutional Neural Networks

abstract: The high energy cost of processing deep convolutional neural networks impedes their ubiquitous deployment in energy-constrained platforms such as embedded systems and IoT devices. This work introduces convolutional layers with pre-defined sparse 2D kernels that have support sets that repeat periodically within and across filters. Due to the efficient storage of our periodic sparse kernels, the parameter savings can translate into considerable improvements in energy efficiency due to reduced DRAM accesses, thus promising significant improvements in the trade-off between energy consumption and accuracy for both training and inference. To evaluate this approach, we performed experiments with two widely accepted datasets, CIFAR-10 and Tiny ImageNet in sparse variants of the ResNet18 and VGG16 architectures. Compared to baseline models, our proposed sparse variants require up to 82% fewer model parameters with 5.6times fewer FLOPs with negligible loss in accuracy for ResNet18 on CIFAR-10. For VGG16 trained on Tiny ImageNet, our approach requires 5.8times fewer FLOPs and up to 83.3% fewer model parameters with a drop in top-5 (top-1) accuracy of only 1.2% (2.1%). We also compared the performance of our proposed architectures with that of ShuffleNet andMobileNetV2. Using similar hyperparameters and FLOPs, our ResNet18 variants yield an average accuracy improvement of 2.8%.


title: Schwinger Boson mean field theory of kagome Heisenberg antiferromagnet with Dzyaloshinskii-Moriya interaction

abstract: We have studied the effect of the Dzyaloshinskii-Moriya interaction on kagome Heisenberg anti-ferromagnet using Schwinger boson mean field theory (SBMFT). Within SBMFT framework, Messio et al had argued that the ground state of kagome antiferromagnet is possibly a chiral topological spin liquid (Phys. Rev. Lett. 108, 207204 (2012)). Thus, we have computed zero-temperature ground state phase diagram considering the time-reversal breaking states as well as fully symmetric Ans\"atze. We discuss the relevance of these results in experiments and other studies. Finally, we have computed the static and dynamic spin structure factors in relevant phases.


title: New biharmonic bases in commutative algebras of the second rank and monogenic functions related to the biharmonic equation

abstract: Among all two-dimensional commutative algebras of the second rank a totally of all their biharmonic bases $\{e_1,e_2\}$, satisfying conditions $\left(e_1^2+ e_2^2\right)^{2} = 0$, $e_1^2 + e_2^2 \ne 0$, is found in an explicit form. A set of "analytic" (monogenic) functions satisfying the biharmonic equation and defined in the real planes generated by the biharmonic bases is built. A characterization of biharmonic functions in bounded simply connected domains by real components of some monogenic functions is found.


title: Multimode interference induced optical nonreciprocity and routing in an optical microcavity

abstract: Optical nonreciprocity and routing using optocal microcavities draw much atttention in recent years. Here, we report the results of the study on the nonreciprocity and routing using optomechanical multimode interference in an optical microcavity. The optomechanical system used here possesses multi-optical modes and a mechanical mode. Optomechanical induced transparency and absorption, appear in the system due to the interference between different paths. The system can present significant nonreciprocity and routing properties when appropriate parameters of the system are set. We design quantum devices, such as diode, circulator and router, which are important applications. Our work shows that optomechanical multimode system can be used as a promising platform for buliding photonic and quantum network.


title: Distribution of quantum coherence and quantum phase transition in the Ising system

abstract: Quantifying of quantum coherence of a given system not only plays an important role in quantum information science but also promote our understanding on some basic problems, such as quantum phase transition. Conventional quantum coherence measurements, such as $l_1$ norm of coherence and relative entropy of coherence, has been widely used to study quantum phase transition, which usually are basis-dependent. The recent quantum version of the Jensen-Shannon divergence meet all the requirements of a good coherence measure. It is not only a metric but also can be basis-independent. Here, based on the quantum renormalization group method we propose an analysis on the critical behavior of two types Ising systems when distribution of quantum coherence. We directly obtain the trade-off relation, critical phenomena, singular behavior, and scaling behavior for both quantum block spin system. Furthermore, the monogamy relation in the multipartite system is also studied in detail. These new result expand the result that quantum coherence can decompose into various contributions as well as enlarge the applications in using basis-independent quantum coherence to reflect quantum critical phenomena.


title: qBSA: Logic Design of a 32-bit Block-Skewed RSFQ Arithmetic Logic Unit

abstract: Single flux quantum (SFQ) circuits are an attractive beyond-CMOS technology because they promise two orders of magnitude lower power at clock frequencies exceeding 25 GHz.However, every SFQ gate is clocked creating very deep gate-level pipelines that are difficult to keep full, particularly for sequences that include data-dependent operations. This paper proposes to increase the throughput of SFQ pipelines by re-designing the datapath to accept and operate on least-significant bits (LSBs) clock cycles earlier than more significant bits. This skewed datapath approach reduces the latency of the LSB side which can be feedback earlier for use in subsequent data-dependent operations increasing their throughput. In particular,we propose to group the bits into 4-bit blocks that are operatedon concurrently and create block-skewed datapath units for 32-bit operation. This skewed approach allows a subsequent data-dependent operation to start evaluating as soon as the first 4-bit block completes. Using this general approach, we developa block-skewed MIPS-compatible 32-bit ALU. Our gate-level Verilog design improves the throughput of 32-bit data dependent operations by 2x and 1.5x compared to previously proposed 4-bit bit-slice and 32-bit Ladner-Fischer ALUs respectively.


title: On-chip deterministic operation of quantum dots in dual-mode waveguides for a plug-and-play single-photon source

abstract: A deterministic source of coherent single photons is an enabling device of quantum-information processing for quantum simulators, and ultimately a full-fledged quantum internet. Quantum dots (QDs) in nanophotonic structures have been employed as excellent sources of single photons, and planar waveguides are well suited for scaling up to multiple photons and emitters exploring near-unity photon-emitter coupling and advanced active on-chip functionalities. An ideal single-photon source requires suppressing noise and decoherence, which notably has been demonstrated in electrically-contacted heterostructures. It remains a challenge to implement deterministic resonant excitation of the QD required for generating coherent single photons, since residual light from the excitation laser should be suppressed without compromising source efficiency and scalability. Here, we present the design and realization of a novel planar nanophotonic device that enables deterministic pulsed resonant excitation of QDs through the waveguide. Through nanostructure engineering, the excitation light and collected photons are guided in two orthogonal waveguide modes enabling deterministic operation. We demonstrate a coherent single-photon source that simultaneously achieves high-purity ($g^{(2)}(0)$ = 0.020 $\pm$ 0.005), high-indistinguishability ($V$ = 96 $\pm$ 2 %), and $>$80 % coupling efficiency into the waveguide. The novel `plug-and-play' coherent single-photon source could be operated unmanned for several days and will find immediate applications, e.g., for constructing heralded multi-photon entanglement sources for photonic quantum computing or sensing.


title: Patient Specific Biomechanics Are Clinically Significant In Accurate Computer Aided Surgical Image Guidance

abstract: Augmented Reality is used in Image Guided surgery (AR IG) to fuse surgical landmarks from preoperative images into a video overlay. Physical simulation is essential to maintaining accurate position of the landmarks as surgery progresses and ensuring patient safety by avoiding accidental damage to vessels etc. In liver procedures, AR IG simulation accuracy is hampered by an inability to model stiffness variations unique to the patients disease. We introduce a novel method to account for patient specific stiffness variation based on Magnetic Resonance Elastography (MRE) data. To the best of our knowledge we are the first to demonstrate the use of in-vivo biomechanical data for AR IG landmark placement. In this early work, a comparative evaluation of our MRE data driven simulation and the traditional method shows clinically significant differences in accuracy during landmark placement and motivates further animal model trials.


title: Environment-aware Reconfigurable Noise Suppression

abstract: The paper proposes an efficient, robust, and reconfigurable technique to suppress various types of noises for any sampling rate. The theoretical analyses, subjective and objective test results show that the proposed noise suppression (NS) solution significantly enhances the speech transmission index (STI), speech intelligibility (SI), signal-to-noise ratio (SNR), and subjective listening experience. The STI and SI consists of 5 levels, i.e., bad, poor, fair, good, and excellent. The most common noisy condition is of SNR ranging from -5 to 8 dB. For the input SNR between -5 and 2.5 dB, the proposed NS improves the STI and SI from "fair" to "good". For the input SNR between 2.5 to 8 dB, the STI and SI are improved from "good" to "excellent" by the proposed NS. The proposed NS can be adopted in both capture and playback paths for voice over internet protocol, voice-trigger, and automatic speech recognition applications.


title: Query-Sequence Optimization on a Reconfigurable Hardware-Accelerated System

abstract: Hardware acceleration of database query processing can be done with the help of FPGAs. In particular, they are partially reconfigurable during runtime, which allows for the runtime adaption of the hardware to a variety of queries. Reconfiguration itself, however, takes some time. As the affected area of the FPGA is not available for computations during the reconfiguration, avoiding some of the reconfigurations can improve overall performance. This paper presents optimizations based on query sequences, which reduces the impact of the reconfigurations. Knowledge of coming queries is used to (I) speculatively start reconfiguration already when a query is still running and (II) avoid overwriting of reconfigurable regions that will be used again in subsequent queries. We evaluate our optimizations with a calibrated model and measurements for various parameter values. Improvements in execution time of up to 21% can be obtained even with sequences of only two queries


title: Turbulent Schmidt Number Measurements Over Three-Dimensional Cubic Arrays

abstract: We present turbulent Schmidt number Sc_t estimations above three-dimensional urban canopies, where Sc_t is a property of the flow defined as the ratio of the eddy diffusivity of momentum (K_M) to the eddy diffusivity of mass (D_t). Despite the fact that Sc_t modelling is of great interest, inter alia, for pollutant dispersion simulations conducted via computational fluid dynamics, no universal value is known. Simultaneous measurements of fluid velocity and mass concentration are carried out in a water channel for three staggered arrays of cubical obstacles corresponding to isolated flow, wake-interference, and skimming-flow regimes. A passive tracer is released from a continuous point source located at a height z=1.67H, where H is the obstacle height. The results show an increase of Sc_t with height above the canopy for all three arrays, with values at z=2H (Sc_t = 0.6) about double compared to that at z=H. The observed Sc_t agrees well with that modelled by using a simple formulation for Sc_t based on expressions for K_M and D_t published in previous studies. Comparisons with other Sc_t models found in the literature are also presented and discussed.


title: Investigation of Numerical Dispersion with Time Step of The FDTD Methods: Avoiding Erroneous Conclusions

abstract: It is widely thought that small time steps lead to small numerical errors in the finite-difference time-domain (FDTD) simulations. In this paper, we investigated how time steps impact on numerical dispersion of two FDTD methods including the FDTD(2,2) method and the FDTD(2,4) method. Through rigorously analytical and numerical analysis, it is found that small time steps of the FDTD methods do not always have small numerical errors. Our findings reveal that these two FDTD methods present different behaviors with respect to time steps: (1) for the FDTD(2,2) method, smaller time steps limited by the Courant-Friedrichs-Lewy (CFL) condition increase numerical dispersion and lead to larger simulation errors; (2) for the FDTD(2,4) method, as time step increases, numerical dispersion errors first decrease and then increase. Our findings are also comprehensively validated from one- to three-dimensional cases through several numerical examples including wave propagation, resonant frequencies of cavities and a practical electromagnetic compatibility (EMC) problem.


title: Stationary and dynamical properties of two harmonically trapped bosons in the crossover from two dimensions to one

abstract: We unravel the stationary properties and the interaction quench dynamics of two bosons, confined in a two-dimensional anisotropic harmonic trap. A transcendental equation is derived giving access to the energy spectrum and revealing the dependence of the energy gaps on the anisotropy parameter. The relation between the two and the one dimensional scattering lengths as well as the Tan contacts is established. The contact, capturing the two-body short range correlations, shows an increasing tendency for a larger anisotropy. Subsequently, the interaction quench dynamics from attractive to repulsive values and vice versa is investigated for various anisotropies. A closed analytical form of the expansion coefficients of the two-body wavefunction, during the time evolution is constructed. The response of the system is studied by means of the time-averaged fidelity, the spectra of the spatial extent of the cloud in each direction and the one-body density. It is found that as the anisotropy increases, the system becomes less perturbed independently of the interactions while for fixed anisotropy quenches towards the non-interacting regime perturb the system in the most efficient manner. Furthermore, we identify that in the tightly confined direction more frequencies are involved in the dynamics stemming from higher-lying excited states.


title: Concise Read-Only Specifications for Better Synthesis of Programs with Pointers -- Extended Version

abstract: In program synthesis there is a well-known trade-off between concise and strong specifications: if a specification is too verbose, it might be harder to write than the program; if it is too weak, the synthesised program might not match the user's intent. In this work we explore the use of annotations for restricting memory access permissions in program synthesis, and show that they can make specifications much stronger while remaining surprisingly concise. Specifically, we enhance Synthetic Separation Logic (SSL), a framework for synthesis of heap-manipulating programs, with the logical mechanism of read-only borrows. We observe that this minimalistic and conservative SSL extension benefits the synthesis in several ways, making it more (a) expressive (stronger correctness guarantees are achieved with a modest annotation overhead), (b) effective (it produces more concise and easier-to-read programs), (c) efficient (faster synthesis), and (d) robust (synthesis efficiency is less affected by the choice of the search heuristic). We explain the intuition and provide formal treatment for read-only borrows. We substantiate the claims (a)--(d) by describing our quantitative evaluation of the borrowing-aware synthesis implementation on a series of standard benchmark specifications for various heap-manipulating programs.


title: Strong Shape Theory of Continuous Maps

abstract: The work is motivated by the papers [Ba1], [Ba2], [Ba7], [Ba11], [Be] and [Be-Tu]. In particular, the strong homology groups of continuous maps were defined and studied in [Be] and [Be-Tu]. To show that given groups are homology type functor, it was required to construct the corresponding shape category. In this paper, we study this very problem. In particular, using the methods developed in [Ba7], [Ba-Ts]. the strong shape theory of continuous maps of compact metric spaces, so-called strong fiber shape theory is constructed.