title: Leveraging Synthetic Data in Object Detection on Unmanned Aerial Vehicles

abstract: Acquiring data to train deep learning-based object detectors on Unmanned Aerial Vehicles (UAVs) is expensive, time-consuming and may even be prohibited by law in specific environments. On the other hand, synthetic data is fast and cheap to access. In this work, we explore the potential use of synthetic data in object detection from UAVs across various application environments. For that, we extend the open-source framework DeepGTAV to work for UAV scenarios. We capture various large-scale high-resolution synthetic data sets in several domains to demonstrate their use in real-world object detection from UAVs by analyzing multiple training strategies across several models. Furthermore, we analyze several different data generation and sampling parameters to provide actionable engineering advice for further scientific research. The DeepGTAV framework is available at https://git.io/Jyf5j.


title: Rapid Transition to Remote Instruction of Physics Labs During Spring 2020: Instructor Perspectives

abstract: Laboratory courses are an important part of the undergraduate physics curriculum. During physics labs, students can engage in authentic, hands-on experimental practices, which can prepare them for graduate school, research laboratories, and jobs in industry. Due to the COVID-19 pandemic in spring 2020, colleges and universities across the world rapidly transitioned to teaching labs remotely. In this work, we report results from a survey of physics lab instructors on how they adapted their courses in the transition to emergency remote teaching. We found that the instructors who responded to the survey faced numerous challenges when transitioning their classes to remote instruction, particularly in providing students with a similar experience to the in-person labs. In addition, we identified common themes in the instructors' responses including changing learning goals of the courses to be more concept-focused, reducing group work due to equity and technological concerns, and using a variety simulation tools, as well as report on factors that the instructors hoped to continue once they have returned to in-person instruction.


title: Dominant Kitaev interactions in the honeycomb materials Na$_3$Co$_2$SbO$_6$ and Na$_2$Co$_2$TeO$_6$

abstract: Cobaltates with 3$d$ based layered honeycomb structure were recently proposed as Kitaev magnets and putative candidates to host the long-sought Kitaev spin liquid. Here we present inelastic neutron scattering results down to 50 mK for powder samples of Na$_3$Co$_2$SbO$_6$ and Na$_2$Co$_2$TeO$_6$, with high resolution in regions of low momentum and energy transfers. We compare the experimental data below the antiferromagnetic zigzag ordering temperature with dynamical structure factors obtained within spin wave theory. We search the wide parameter range of a $K$-$J_1$-$\Gamma$-$\Gamma^{\prime}$-$J_3$ spin 1/2 model and identify the best fits to constant momentum cuts of the inelastic neutron data. The powder average limits the selection of a unique parameter set for each material, but we see clear trends towards ferromagnetic Kitaev exchange in both compounds, and ratios $|K/J_1|$ could be as large as $5\ldots 25$; in contrast, antiferromagnetic Kitaev exchange cannot be ruled out, but requires a fine-tuning of all involved spin exchange couplings with only moderate ratio $|K/J_1| \sim 1$. Our experimental data are incompatible with a purely isotropic Heisenberg model.


title: Entropy-Regularized Partially Observed Markov Decision Processes

abstract: We investigate partially observed Markov decision processes (POMDPs) with cost functions regularized by entropy terms describing state, observation, and control uncertainty. Standard POMDP techniques are shown to offer bounded-error solutions to these entropy-regularized POMDPs, with exact solutions possible when the regularization involves the joint entropy of the state, observation, and control trajectories. Our joint-entropy result is particularly surprising since it constitutes a novel, tractable formulation of active state estimation.


title: Mod $p$ local-global compatibility for $\mathrm{GSp}_4(\mathbb{Q}_p)$ in the ordinary case

abstract: Let $F$ be a totally real field of even degree in which $p$ splits completely. Let $\overline{r}:G_F \rightarrow \mathrm{GSp}_4(\overline{\mathbb{F}}_p)$ be a modular Galois representation unramified at all finite places away from $p$ and upper-triangular, maximally nonsplit, and of parallel weight at places dividing $p$. Fix a place $w$ dividing $p$. Assuming certain genericity conditions and Taylor--Wiles assumptions, we prove that the $\mathrm{GSp}_4(F_w)$-action on the corresponding Hecke-isotypic part of the space of mod $p$ automorphic forms on a compact mod center form of $\mathrm{GSp}_4$ with infinite level at $w$ determines $\overline{r}|_{G_{F_w}}$.


title: Fully coupled drift-less Forward and backward Stochastic Differential Equations in a degenerate case

abstract: Existence and uniqueness results of fully coupled forward stochastic differential equations without drifts and backward stochastic differential equations in a degenerate case are obtained for an arbitrarily large time duration.


title: The Wavelength-shifting Optical Module

abstract: The Wavelength-shifting Optical Module (WOM) is a novel photosensor concept for the instrumentation of large detector volumes with single-photon sensitivity. The key objective is to improve the signal-to-noise ratio which is achieved by decoupling the photosensitive area of a sensor from the cathode area of its photomultiplier tube (PMT). The WOM consists of a transparent tube with two PMTs attached to its ends. The tube is coated with wavelength-shifting paint absorbing ultra-violet photons with nearly $100\,\%$ efficiency. Depending on the environment, e.g. air (ice), up to $73\,\%$ $(41\,\%)$ of the subsequently emitted optical photons can be captured by total internal reflection and propagate towards the PMTs where they are recorded. The optical properties of the paint, the geometry of the tube and the coupling of the tube to the PMTs have been optimized for maximal sensitivity based on theoretical derivations and experimental evaluations. Prototypes were built to demonstrate the technique and to develop a reproducible construction process. Important measurable characteristics of the WOM are the wavelength dependent effective area, the transit time spread of detected photons and the signal-to-noise ratio. The WOM outperforms bare PMTs especially with respect to the low signal-to-noise ratio with an increase of a factor up to 8.9 in air (5.2 in ice). Since the gain in sensitivity is mostly in the UV-regime, the WOM is an ideal sensor for Cherenkov and scintillation detectors.


title: Density Regression with Bayesian Additive Regression Trees

abstract: Flexibly modeling how an entire density changes with covariates is an important but challenging generalization of mean and quantile regression. While existing methods for density regression primarily consist of covariate-dependent discrete mixture models, we consider a continuous latent variable model in general covariate spaces, which we call DR-BART. The prior mapping the latent variable to the observed data is constructed via a novel application of Bayesian Additive Regression Trees (BART). We prove that the posterior induced by our model concentrates quickly around true generative functions that are sufficiently smooth. We also analyze the performance of DR-BART on a set of challenging simulated examples, where it outperforms various other methods for Bayesian density regression. Lastly, we apply DR-BART to two real datasets from educational testing and economics, to study student growth and predict returns to education. Our proposed sampler is efficient and allows one to take advantage of BART's flexibility in many applied settings where the entire distribution of the response is of primary interest. Furthermore, our scheme for splitting on latent variables within BART facilitates its future application to other classes of models that can be described via latent variables, such as those involving hierarchical or time series data.


title: Quantitative assessment of fitting errors associated with streak camera noise in Thomson scattering data analysis

abstract: Thomson scattering measurements in High Energy Density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominate source of measurement uncertainty. In this paper we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data. These methods will have general applicability to other measurements made using optical streak cameras.


title: What Does Nature Minimize In Every Incompressible Flow?

abstract: In this paper, we discover the fundamental quantity that Nature minimizes in almost all flows encountered in everyday life: river, rain, flow in a pipe, blood flow, airflow over an airplane, etc. We show that the norm of the pressure gradient over the field is minimum at every instant of time! We call it the principle of minimum pressure gradient (PMPG). The principle is deeply rooted in classical mechanics via Gauss' principle of least constraint. Therefore, while we prove mathematically that Navier-Stokes' equation represents the necessary condition for minimization of the pressure gradient, the PMPG stands on its own philosophy independent of Navier-Stokes'. It turns any fluid mechanics problem into a minimization one. We demonstrate this intriguing property by solving three of the classical problems in fluid mechanics using the PMPG without resorting to Navier-Stokes' equation. In fact, the inviscid version of the PMPG allowed solving the long-standing problem of the aerohydrodynamic lift over smooth cylindrical shapes where Euler's equation fails to provide a unique answer. Moreover, the result challenges the accepted wisdom about lift generation on an airfoil, which has prevailed over a century. The PMPG is expected to be transformative for theoretical modeling of fluid mechanics as it encodes a complicated nonlinear partial differential equation into a simple minimization problem. The principle even transcends Navier-Stokes' equations in its applicability to non-Newtonian fluids with arbitrary constitutive relations and fluids subject to arbitrary forcing (e.g. electric or magnetic).


title: Morphological classifiers

abstract: This work proposes a new type of classifier called Morphological Classifier (MC). MCs aggregate concepts from mathematical morphology and supervised learning. The outcomes of this aggregation are classifiers that may preserve shape characteristics of classes, subject to the choice of a stopping criterion and structuring element. MCs are fundamentally based on set theory, and their classification model can be a mathematical set itself. Two types of morphological classifiers are proposed in the current work, namely, Morphological k-NN (MkNN) and Morphological Dilation Classifier (MDC), which demonstrate the feasibility of the approach. This work provides evidence regarding the advantages of MCs, e.g., very fast classification times as well as competitive accuracy rates. The performance of MkNN and MDC was tested using p -dimensional datasets. MCs tied or outperformed 14 well established classifiers in 5 out of 8 datasets. In all occasions, the obtained accuracies were higher than the average accuracy obtained with all classifiers. Moreover, the proposed implementations utilize the power of the Graphics Processing Units (GPUs) to speed up processing.


title: Crash Data Augmentation Using Conditional Generative Adversarial Networks (CGAN) for Improving Safety Performance Functions

abstract: In this paper, we present a crash frequency data augmentation method based on Conditional Generative Adversarial Networks to improve crash frequency models. The proposed method is evaluated by comparing the performance of Base SPFs (developed using original data) and Augmented SPFs (developed using original data plus synthesised data) in terms of hotspot identification performance, model prediction accuracy, and dispersion parameter estimation accuracy. The experiments are conducted using simulated and real-world crash data sets. The results indicate that the synthesised crash data by CGAN have the same distribution as the original data and the Augmented SPFs outperforms Base SPFs in almost all aspects especially when the dispersion parameter is low.


title: Co-interacting dark matter and conformally coupled light scalars

abstract: The increasing observational pressure on the standard cosmological model motivates analyses going beyond the paradigm of the collision-less cold dark matter (DM). Since the only clear evidence for the existence of DM is based on gravitational interactions, it seems particularly fitting to study them in this sector where extensions to the standard model can be naturally introduced. A promising avenue can be obtained using modifications of the space-time metric coupled to DM and induced by the presence of a new ultra-light scalar field $\phi$. The $\phi$ field can contribute to the DM density and can couple all the matter species universally, including additional heavy DM particles. We present a simple two-component DM model employing derivative conformal interactions between the two DM species. This can simultaneously: 1) guarantee the necessary symmetries to stabilize the dark species, 2) predict a subdominant thermal relic density of the heavy DM component, and 3) alleviate small-scale structure tensions of the cold DM scenario due to the possible co-interactions in the dark sector. The scenario is highly predictive with future observational prospects ranging from the Large Hadron Collider (LHC) to gravitational-wave searches, and can be generalized to more rich and realistic dark sector models.


title: Gravitational waveforms for compact binaries from second-order self-force theory

abstract: We produce gravitational waveforms for nonspinning compact binaries undergoing a quasicircular inspiral. Our approach is based on a two-timescale expansion of the Einstein equations in second-order self-force theory, which allows first-principles waveform production in tens of milliseconds. Although the approach is designed for extreme mass ratios, our waveforms agree remarkably well with those from full numerical relativity, even for comparable-mass systems. Our results will be invaluable in accurately modelling extreme-mass-ratio inspirals for the LISA mission and intermediate-mass-ratio systems currently being observed by the LIGO-Virgo-KAGRA Collaboration.


title: Quantum vortex instability and black hole superradiance

abstract: Vortices and black holes set the scene for many interesting dynamical processes in physics. Here, we study the dynamical instability of quantised vortices and rotational superradiance around rotating black holes, illustrating in the process that the same physics is at play in these two seemingly disparate phenomena. We also compare the instability of the vortex to the black hole bomb instability, which occurs for massive scalar fields in the Kerr spacetime. Taking inspiration from the analogy between black hole bomb modes and the hydrogen spectrum, the vortex instability is compared with nuclear resonances involved in $\alpha$-decay.


title: Symmetry of magnetic correlations in spin-triplet superconductor UTe2

abstract: The temperature dependence of the low-energy magnetic excitations in the spin-triplet superconductor UTe$_2$ was measured via inelastic neutron scattering in the normal and superconducting states. The imaginary part of the dynamic susceptibility follows the behavior of interband correlations in a hybridized Kondo lattice with an appropriate characteristic energy. These excitations are a lower-dimensional analogue of phenomena observed in other Kondo lattice materials, such that their presence is not necessarily due to dominance of ferromagnetic or antiferromagnetic correlations. The onset of superconductivity alters the magnetic excitations noticeably on the same energy scales, suggesting that these changes originate from additional electronic structure modification.


title: An algebraic attack on stream ciphers with application to nonlinear filter generators and WG-PRNG

abstract: In this paper, we propose a new algebraic attack on stream ciphers. Starting from the well-known attack due to Courtois and Meier, we design an attack especially effective against nonlinear filter generators. We test it on two toy stream ciphers and we show that the level of security of one of stream ciphers submitted to the NIST competition on Lightweight Cryptography, WG-PRNG, is less than that stated before now.


title: Angular Momentum Eigenstates of the Isotropic 3-D Harmonic Oscillator: Phase-Space Distributions and Coalescence Probabilities

abstract: The isotropic 3-dimensional harmonic oscillator potential can serve as an approximate description of many systems in atomic, solid state, nuclear, and particle physics. In particular, the question of 2 particles binding (or coalescing) into angular momentum eigenstates in such a potential has interesting applications. We compute the probabilities for coalescence of two distinguishable, non-relativistic particles into such a bound state, where the initial particles are represented by generic wave packets of given average positions and momenta. We use a phase-space formulation and hence need the Wigner distribution functions of angular momentum eigenstates in isotropic 3-dimensional harmonic oscillators. These distribution functions have been discussed in the literature before but we utilize an alternative approach to obtain these functions. Along the way, we derive a general formula that expands angular momentum eigenstates in terms of products of 1-dimensional harmonic oscillator eigenstates.


title: High-resolution, quantitative signal subspace imaging for synthetic aperture radar

abstract: We consider synthetic aperture radar imaging of a region containing point-like targets. Measurements are the set of frequency responses to scattering by the targets taken over a collection of individual spatial locations along the flight path making up the synthetic aperture. Because signal subspace imaging methods do not work on these measurements directly, we rearrange the frequency response at each spatial location using the Prony method and obtain a matrix that is suitable for these methods. We arrange the set of these Prony matrices as one block-diagonal matrix and introduce a signal subspace imaging method for it. We show that this signal subspace method yields high-resolution and quantitative images provided that the signal-to-noise ratio is sufficiently high. We give a resolution analysis for this imaging method and validate this theory using numerical simulations. Additionally, we show that this imaging method is stable to random perturbations to the travel times and validate this theory with numerical simulations using the random travel time model for random media.


title: Efficient $d$-dimensional molecular dynamics simulations for studies of the glass-jamming transition

abstract: We develop an algorithm suitable for parallel molecular dynamics simulations in $d$ spatial dimensions and describe its implementation in C++. All routines work in arbitrary $d$; the maximum simulated $d$ is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass/jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. Scaling of simulation runtimes with the number of particles $N$ and number of simulation threads $n_{\rm threads}$ is comparable to popular MD codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled $d = 6$ liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.


title: Human Activity Recognition on wrist-worn accelerometers using self-supervised neural networks

abstract: Measures of Activity of Daily Living (ADL) are an important indicator of overall health but difficult to measure in-clinic. Automated and accurate human activity recognition (HAR) using wrist-worn accelerometers enables practical and cost efficient remote monitoring of ADL. Key obstacles in developing high quality HAR is the lack of large labeled datasets and the performance loss when applying models trained on small curated datasets to the continuous stream of heterogeneous data in real-life. In this work we design a self-supervised learning paradigm to create a robust representation of accelerometer data that can generalize across devices and subjects. We demonstrate that this representation can separate activities of daily living and achieve strong HAR accuracy (on multiple benchmark datasets) using very few labels. We also propose a segmentation algorithm which can identify segments of salient activity and boost HAR accuracy on continuous real-life data.


title: Perceptual Evaluation of 360 Audiovisual Quality and Machine Learning Predictions

abstract: In an earlier study, we gathered perceptual evaluations of the audio, video, and audiovisual quality for 360 audiovisual content. This paper investigates perceived audiovisual quality prediction based on objective quality metrics and subjective scores of 360 video and spatial audio content. Thirteen objective video quality metrics and three objective audio quality metrics were evaluated for five stimuli for each coding parameter. Four regression-based machine learning models were trained and tested here, i.e., multiple linear regression, decision tree, random forest, and support vector machine. Each model was constructed using a combination of audio and video quality metrics and two cross-validation methods (k-Fold and Leave-One-Out) were investigated and produced 312 predictive models. The results indicate that the model based on the evaluation of VMAF and AMBIQUAL is better than other combinations of audio-video quality metric. In this study, support vector machine provides higher performance using k-Fold (PCC = 0.909, SROCC = 0.914, and RMSE = 0.416). These results can provide insights for the design of multimedia quality metrics and the development of predictive models for audiovisual omnidirectional media.


title: Projection theorems for linear-fractional families of projections

abstract: Marstrand's theorem states that applying a generic rotation to a planar set $A$ before projecting it orthogonally to the $x$-axis almost surely gives an image with the maximal possible dimension $\min(1, \dim A)$. We first prove, using the transversality theory of Peres-Schlag locally, that the same result holds when applying a generic complex linear-fractional transformation in $PSL(2,\C)$ or a generic real linear-fractional transformation in $PGL(3,\R)$. We next show that, under some necessary technical assumptions, transversality locally holds for restricted families of projections corresponding to one-dimensional subgroups of $PSL(2,\C)$ or $PGL(3,\R)$. Third, we demonstrate, in any dimension, local transversality and resulting projection statements for the families of closest-point projections to totally-geodesic subspaces of hyperbolic and spherical geometries.


title: A Simplicity Bubble Problem in Formal-Theoretic Learning Systems

abstract: When mining large datasets in order to predict new data, limitations of the principles behind statistical machine learning pose a serious challenge not only to the Big Data deluge, but also to the traditional assumptions that data generating processes are biased toward low algorithmic complexity. Even when one assumes an underlying algorithmic-informational bias toward simplicity in finite dataset generators, we show that current approaches to machine learning (including deep learning, or any formal-theoretic hybrid mix of top-down AI and statistical machine learning approaches), can always be deceived, naturally or artificially, by sufficiently large datasets. In particular, we demonstrate that, for every learning algorithm (with or without access to a formal theory), there is a sufficiently large dataset size above which the algorithmic probability of an unpredictable deceiver is an upper bound (up to a multiplicative constant that only depends on the learning algorithm) for the algorithmic probability of any other larger dataset. In other words, very large and complex datasets can deceive learning algorithms into a ``simplicity bubble'' as likely as any other particular non-deceiving dataset. These deceiving datasets guarantee that any prediction effected by the learning algorithm will unpredictably diverge from the high-algorithmic-complexity globally optimal solution while converging toward the low-algorithmic-complexity locally optimal solution, although the latter is deemed a global one by the learning algorithm. We discuss the framework and additional empirical conditions to be met in order to circumvent this deceptive phenomenon, moving away from statistical machine learning towards a stronger type of machine learning based on, and motivated by, the intrinsic power of algorithmic information theory and computability theory.

categories: cs.IT cs.AI cs.LG cs.LO math.IT math.ST stat.TH


title: K-polystability of 3-dimensional log Fano pairs of Maeda type

abstract: Using the Abban-Zhuang theory and the classification of three-dimensional log smooth log Fano pairs due to Maeda, we prove that threefold log Fano pairs $(X, D)$ of Maeda type with reducible boundary $D$ are K-unstable, with four exceptions. We also correct several inaccuracies in Maeda's classification.


title: High time resolution broad-band polarimetry: technique, calibration and standards

abstract: Regular large-scale polarimetric observations in Crimean astrophysical observatory began in the early 1960s. In 2002 - 2017 the single-channel aperture photometer-polarimeter with a quarter-wave plate at the 2.6-m Shajn mirror telescope (SMT) was used. We accumulated a large homogeneous data set of polarimetric observations of different types of objects that are to be published separately. Correct polarimetric data processing requires high polarization standards and zero-polarization stars. We aim to improve the data reduction and calibration process to obtain further results with highest possible accuracy. High time resolution broad-band (WR, R, V, B, U) polarization observations are made of 98 known standard stars (527 time series with total duration about 184 hours). We determined values of linear and circular polarization for 98 nearby Northern bright stars. This catalogue is not compilative, but obtained using the same instrument and technique during large time interval. It will be used for our future research and it may be used by other authors. We implemented the least squares approach for determination of the Stokes parameters. It allowed us to obtain results with the accuracy better then we obtained using previously used methods. We report suspicious or variable stars that are not suitable as standards for high precision polarimetry.


title: A rodent paradigm for studying perceptual decisions under asymmetric reward

abstract: Many real-life decisions involve both perceptual processes and weighing the consequences of different actions. However, the neural mechanisms underlying perceptual decisions have typically been examined separately from those underlying economic decisions. Here, we trained rats to make choices informed by both perceptual and value cues on a trial-by-trial basis. As in typical perceptual tasks, subjects were rewarded for correctly categorizing a tone relative to a learned threshold. To add an economic component, a light indicated, on each trial, whether correct responses to one side gave higher rewards than correct responses to the other side. As such, on trials with some perceptual uncertainty, it could be worthwhile to choose the unlikely option, if it had higher expected value. We found that, despite subjects sensitivity to the frequency of the cue and the reward sizes, their behavior was not optimal: subjects tended to shift their choices in a stimulus-independent way following light flashes. Moreover, subjects tended to under-shift, which could be interpreted as being over-confident in their perceptual beliefs or as being risk-averse.


title: Randomize the Future: Asymptotically Optimal Locally Private Frequency Estimation Protocol for Longitudinal Data

abstract: Longitudinal data tracking under Local Differential Privacy (LDP) is a challenging task. Baseline solutions that repeatedly invoke a protocol designed for one-time computation lead to linear decay in the privacy or utility guarantee with respect to the number of computations. To avoid this, the recent approach of Erlingsson et al. (2020) exploits the potential sparsity of user data that changes only infrequently. Their protocol targets the fundamental problem of frequency estimation protocol for longitudinal binary data, with $\ell_\infty$ error of $O ( (1 / \epsilon) \cdot (\log d)^{3 / 2} \cdot k \cdot \sqrt{ n \cdot \log ( d / \beta ) } )$, where $\epsilon$ is the privacy budget, $d$ is the number of time periods, $k$ is the maximum number of changes of user data, and $\beta$ is the failure probability. Notably, the error bound scales polylogarithmically with $d$, but linearly with $k$. In this paper, we break through the linear dependence on $k$ in the estimation error. Our new protocol has error $O ( (1 / \epsilon) \cdot (\log d) \cdot \sqrt{ k \cdot n \cdot \log ( d / \beta ) } )$, matching the lower bound up to a logarithmic factor. The protocol is an online one, that outputs an estimate at each time period. The key breakthrough is a new randomizer for sequential data, FutureRand, with two key features. The first is a composition strategy that correlates the noise across the non-zero elements of the sequence. The second is a pre-computation technique which, by exploiting the symmetry of input space, enables the randomizer to output the results on the fly, without knowing future inputs. Our protocol closes the error gap between existing online and offline algorithms.


title: Nonnegative OPLS for Supervised Design of Filter Banks: Application to Image and Audio Feature Extraction

abstract: Audio or visual data analysis tasks usually have to deal with high-dimensional and nonnegative signals. However, most data analysis methods suffer from overfitting and numerical problems when data have more than a few dimensions needing a dimensionality reduction preprocessing. Moreover, interpretability about how and why filters work for audio or visual applications is a desired property, especially when energy or spectral signals are involved. In these cases, due to the nature of these signals, the nonnegativity of the filter weights is a desired property to better understand its working. Because of these two necessities, we propose different methods to reduce the dimensionality of data while the nonnegativity and interpretability of the solution are assured. In particular, we propose a generalized methodology to design filter banks in a supervised way for applications dealing with nonnegative data, and we explore different ways of solving the proposed objective function consisting of a nonnegative version of the orthonormalized partial least-squares method. We analyze the discriminative power of the features obtained with the proposed methods for two different and widely studied applications: texture and music genre classification. Furthermore, we compare the filter banks achieved by our methods with other state-of-the-art methods specifically designed for feature extraction.


title: Improving the Efficiency of Off-Policy Reinforcement Learning by Accounting for Past Decisions

abstract: Off-policy learning from multistep returns is crucial for sample-efficient reinforcement learning, particularly in the experience replay setting now commonly used with deep neural networks. Classically, off-policy estimation bias is corrected in a per-decision manner: past temporal-difference errors are re-weighted by the instantaneous Importance Sampling (IS) ratio (via eligibility traces) after each action. Many important off-policy algorithms such as Tree Backup and Retrace rely on this mechanism along with differing protocols for truncating ("cutting") the ratios ("traces") to counteract the excessive variance of the IS estimator. Unfortunately, cutting traces on a per-decision basis is not necessarily efficient; once a trace has been cut according to local information, the effect cannot be reversed later, potentially resulting in the premature truncation of estimated returns and slower learning. In the interest of motivating efficient off-policy algorithms, we propose a multistep operator that permits arbitrary past-dependent traces. We prove that our operator is convergent for policy evaluation, and for optimal control when targeting greedy-in-the-limit policies. Our theorems establish the first convergence guarantees for many existing algorithms including Truncated IS, Non-Markov Retrace, and history-dependent TD($\lambda$). Our theoretical results also provide guidance for the development of new algorithms that jointly consider multiple past decisions for better credit assignment and faster learning.


title: Tolman-Ehrenfest-Klein Law in non-Riemannian geometries

abstract: Heat always flows from hotter to a colder temperature until thermal equilibrium be finally restored in agreement with the usual (zeroth, first and second) laws of thermodynamics. However, Tolman and Ehrenfest demonstrated that the relation between inertia and weight uniting all forms of energy in the framework of general relativity implies that the standard equilibrium condition is violated in order to maintain the validity of the first and second law of thermodynamics. Here we demonstrate that the thermal equilibrium condition for a static self-gravitating fluid, besides being violated, is also heavily dependent on the underlying spacetime geometry (whether Riemannian or non-Riemannian). As a particular example, a new equilibrium condition is deduced for a large class of Weyl and f(R) type gravity theories. Such results suggest that experiments based on the foundations of the heat theory (thermal sector) may also be used for confronting gravity theories and prospect the intrinsic geometric nature of the spacetime structure.


title: A combinatorial proof of the Gaussian product inequality beyond the MTP${}_2$ case

abstract: A combinatorial proof of the Gaussian product inequality (GPI) is given under the assumption that each component of a centered Gaussian random vector $\boldsymbol{X} = (X_1, \ldots, X_d)$ of arbitrary length can be written as a linear combination, with coefficients of identical sign, of the components of a standard Gaussian random vector. This condition on $\boldsymbol{X}$ is shown to be strictly weaker than the assumption that the density of the random vector $(|X_1|, \ldots, |X_d|)$ is multivariate totally positive of order $2$, abbreviated MTP${}_2$, for which the GPI is already known to hold. Under this condition, the paper highlights a new link between the GPI and the monotonicity of a certain ratio of gamma functions.


title: A Survey on Perceptually Optimized Video Coding

abstract: To provide users with more realistic visual experiences, videos are developing in the trends of Ultra High Definition (UHD), High Frame Rate (HFR), High Dynamic Range (HDR), Wide Color Gammut (WCG) and high clarity. However, the data amount of videos increases exponentially, which requires high efficiency video compression for storage and network transmission. Perceptually optimized video coding aims to maximize compression efficiency by exploiting visual redundancies. In this paper, we present a broad and systematic survey on perceptually optimized video coding. Firstly, we present problem formulation and framework of the perceptually optimized video coding, which includes visual perception modelling, visual quality assessment and perceptual video coding optimization. Secondly, recent advances on visual factors, computational perceptual models and quality assessment models are presented. Thirdly, we review perceptual video coding optimizations from four key aspects, including perceptually optimized bit allocation, rate-distortion optimization, transform and quantization, filtering and enhancement. In each part, problem formulation, working flow, recent advances, advantages and challenges are presented. Fourthly, perceptual coding performances of the latest coding standards and tools are experimentally analyzed. Finally, challenging issues and future opportunities are identified.


title: Transition pathways in Cylinder-Gyroid interface

abstract: When two distinct ordered phases contact, the interface may exhibit rich and fascinating structures. Focusing on the Cylinder-Gyroid interface system, transition pathways connecting various interface morphologies are studied armed with the Landau--Brazovskii model. Specifically, minimum energy paths are obtained by computing transition states with the saddle dynamics. We present four primary transition pathways connecting different local minima, representing four different mechanisms of the formation of the Cylinder-Gyroid interface. The connection of Cylinder and Gyroid can be either direct or indirect via Fddd with three different orientations. Under different displacements, each of the four pathways may have the lowest energy.


title: Dynamic friction unraveled using an atomically defined model system

abstract: The pervasive phenomenon of friction has been studied at the nanoscale by controlled manipulation of single atoms and molecules, which permitted a precise determination of the static friction force necessary to initiate motion. However, much less is known about the atomistic dynamics during manipulation. Here we reveal the complete manipulation process of a carbon monoxide molecule on a copper surface at low temperatures using a combination of atomic force microscopy, vibrational spectroscopy for different isotope molecules and density functional theory. We measured the energy dissipation associated with manipulation and relate its origin to hysteresis involving an intermediate state, which enables an atomistic interpretation of dynamic friction. Our results show how friction forces can be controlled and optimized, facilitating new fundamental insights for tribology.


title: Motility and Swimming: Universal Description and Generic Trajectories

abstract: Autonomous locomotion is a ubiquitous phenomenon in biology and in physics of active systems at microscopic scale. This includes prokaryotic, eukaryotic cells (crawling and swimming) and artificial swimmers. An outstanding feature is the ability of these entities to follow complex trajectories, ranging from straight, curved (circular, helical...), to random-like ones. The non-straight nature of these trajectories is often explained as a consequence of the asymmetry of the particle or the medium in which it moves, or due to the presence of bounding walls, etc... Here, we show that straight, circular and helical trajectories emerge naturally in the absence of asymmetry of the swimmer or that of suspending medium. Our first proof is based on general considerations, without referring to an explicit form of a model. We show that these three trajectories correspond to self-congruent solutions. Self-congruency means that the states of the system at different moments of time can be made identical by an appropriate combination of rotation and translation of the coordinate space. We show that these solutions are exhibited by spherically symmetric particles as a result of a series of pitchfork bifurcations as the activity is increased. Self-congruent dynamics in one and two dimensions are analyzed as well. Finally, we present a simple explicit nonlinear exactly solvable model of fully isotropic phoretic particle that shows the transitions from a non-motile state to straight motion to circular motion to helical motion as a series of spontaneous symmetry-breaking bifurcations. Whether a system exhibits or not a given trajectory only depends on the numerical values of parameters entering the model, while asymmetry of swimmer shape, or anisotropy of the suspending medium , or influence of bounding walls are not necessary.


title: Safety and Liveness Guarantees through Reach-Avoid Reinforcement Learning

abstract: Reach-avoid optimal control problems, in which the system must reach certain goal conditions while staying clear of unacceptable failure modes, are central to safety and liveness assurance for autonomous robotic systems, but their exact solutions are intractable for complex dynamics and environments. Recent successes in reinforcement learning methods to approximately solve optimal control problems with performance objectives make their application to certification problems attractive; however, the Lagrange-type objective used in reinforcement learning is not suitable to encode temporal logic requirements. Recent work has shown promise in extending the reinforcement learning machinery to safety-type problems, whose objective is not a sum, but a minimum (or maximum) over time. In this work, we generalize the reinforcement learning formulation to handle all optimal control problems in the reach-avoid category. We derive a time-discounted reach-avoid Bellman backup with contraction mapping properties and prove that the resulting reach-avoid Q-learning algorithm converges under analogous conditions to the traditional Lagrange-type problem, yielding an arbitrarily tight conservative approximation to the reach-avoid set. We further demonstrate the use of this formulation with deep reinforcement learning methods, retaining zero-violation guarantees by treating the approximate solutions as untrusted oracles in a model-predictive supervisory control framework. We evaluate our proposed framework on a range of nonlinear systems, validating the results against analytic and numerical solutions, and through Monte Carlo simulation in previously intractable problems. Our results open the door to a range of learning-based methods for safe-and-live autonomous behavior, with applications across robotics and automation. See https://github.com/SafeRoboticsLab/safety_rl for code and supplementary material.


title: Photoinduced charge density wave transition like a puppet on a string

abstract: Charge density wave (CDW) materials can undergo an ultrafast phase transition after an ultrashort laser pulse excitation, and the suggested underlying mechanisms have always been associated with two main features: excitonic interaction-induced CDW charge order and electron-phonon coupling-induced periodic lattice distortion (PLD). Here, beyond these two mechanisms, we reveal that photoexcitation induced CDW phase transition in the prototypic CDW example 1T-TiSe2 is similar to a puppet on a string: six Ti-Se bonds connected to each distorted Ti atom acting as six strings controlling the PLD and in turn the CDW orders. The photoexcitation induced modulation on charge population of the Ti-Se bonds generates a laser-fluence-dependent interatomic-repulsive force along each Ti-Se bond. The nonequal length of these six Ti-Se bonds gives rise to a net force exerted on the central distorted Ti atom to push it toward the suppressing of the PLD and thus the CDW orders. We further illustrate that the dynamics of each distorted Ti atom behaves as though it attached to a spring in a simple harmonic motion with a fluence-dependent oscillation frequency, uniting two previously reported scaling laws for phase transition time. These findings significantly advance the understanding of CDW instability and provide new insights into how photoexcitation induced modulation on charge population may lead to phase transitions by directly connecting interatomic forces with reaction pathways.


title: Insula Interoception, Active Inference and Feeling Representation

abstract: The body sends interoceptive visceral information through deep brain structures to the cerebral cortex. The insula cortex, organized in hierarchical modules, is the major cortical region receiving interoceptive afferents and contains visceral topographic maps. Yet, the biological significance of the insula's modular architecture in relation to deep brain regions remains unsolved. In this opinion, we propose the Insula Hierarchical Modular Adaptive Interoception Control (IMAC) model to suggest that insula modules (granular, dysgranular and agranular subregions), forming networks with prefrontal (supplementary motor area, dorsolateral and ventromedial cortices) and striatum (posterior, dorsomedial and ventromedial) subregions, are specialized for higher-order interoceptive representations, recruited in a context-dependent manner to support habitual, model-based and exploratory adaptive behavior. We then discuss how insula interoceptive representations, or metaceptions, could give rise to conscious interoceptive feelings built up from low-order visceral representations and associated basic emotions located in deep interoceptive brain structures.


title: A Geometric Classification of the Holomorphic Vertex Operator Algebras of Central Charge 24

abstract: We associate with a generalised deep hole of the Leech lattice vertex operator algebra a generalised hole diagram. We show that this Dynkin diagram determines the generalised deep hole up to conjugacy and that there are exactly 70 such diagrams. In an earlier work we proved a bijection between the generalised deep holes and the strongly rational, holomorphic vertex operator algebras of central charge 24 with non-trivial weight-1 space. Hence, we obtain a new, geometric classification of these vertex operator algebras.


title: Long-term secure distributed storage using quantum key distribution network with third-party verification

abstract: The quantum key distribution network with Vernam's One Time Pad encryption and secret sharing are powerful security tools to realize an information theoretically secure distributed storage system. In our previous work, a single-password-authenticated secret sharing scheme based on the QKD network and Shamir's secret sharing was experimentally demonstrated; it confirmed ITS data transmission, storage, authentication, and integrity. To achieve data integrity, an ITS message authentication code tag is employed and a data owner of the secret sharing performs both the MAC tag generation and verification. However, for a scenario in which the data owner and end users are different entities, the above approach may not work since the data owner can cheat the end users. In this paper, we resolve this problem by proposing an ITS integrity protection scheme employing a third-party verification with time-stamp.


title: Extremal isolated horizons with $\Lambda$ and the related unique type D black holes

abstract: We extend our previous work in which we derived the most general form of an induced metric describing the geometry of an axially symmetric extremal isolated horizon (EIH) in asymptotically flat spacetime. Here we generalize it to EIHs in asymptotically (anti-)de Sitter spacetime. The resulting metric conveniently forms a 6-parameter family which, in addition to a cosmological constant $\Lambda$, depends on the area of the horizon, total electric and magnetic charges, and two deficit angles representing conical singularities at poles. Such a metric is consistent with results obtained in the context of near-horizon geometries. Moreover, we study extremal horizons of all black holes within the class of Plebanski-Demianski exact (electro)vacuum spacetimes of the algebraic type D. In an important special case of non-accelerating black holes, that is the famous Kerr-Newman-NUT-(A)dS metric, we were able to identify the corresponding extremal horizons, including their position and geometry, and find explicit relations between the physical parameters of the metric and the geometrical parameters of the EIHs.


title: Run-of-Mine Stockyard Recovery Scheduling and Optimisation for Multiple Reclaimers

abstract: Stockpiles are essential in the mining value chain, assisting in maximising value and production. Quality control of taken minerals from the stockpiles is a major concern for stockpile managers where failure to meet some requirements can lead to losing money. This problem was recently investigated using a single reclaimer, and basic assumptions. This study extends the approach to consider multiple reclaimers in preparing for short and long-term deliveries. The engagement of multiple reclaimers complicates the problem in terms of their interaction in preparing a delivery simultaneously and safety distancing of reclaimers. We also consider more realistic settings, such as handling different minerals with different types of reclaimers. We propose methods that construct a solution step by step to meet precedence constraints for all reclaimers in the stockyard. We study various instances of the problem using greedy algorithms, Ant Colony Optimisation (ACO), and propose an integrated local search method determining an efficient schedule. We fine-tune and compare the algorithms and show that the ACO combined with local search can yield efficient solutions.


title: From Finite Vector Field Data to Combinatorial Dynamical Systems in the Sense of Forman

abstract: The main goal is to construct a combinatorial dynamical system in the sense of Forman from finite vector field data. We use a linear minimization problem with binary variables and linear equality constraints. The solution of the minimization problem induces an admissible matching for the combinatorial dynamical system. They are three main steps for the method: Construct a simplicial complex, compute a vector for each simplex, solve the minimization problem, and apply the induced matching. We argue the effectiveness of the method by testing it on the Lotka-Volterra model and the Lorenz attractor model. We are able to retrieve the cyclic behaviour in the Lotka-Volterra model and the chaotic behaviour of the Lorenz attractor. Two extensions to the algorithm are shown. We use the barycentric subdivision to obtain a better resolution. We add conditions on the minimization problem to obtain a solution which induce a gradient matching.


title: Sub-Chain Beam for mmWave Devices: A Trade-off between Power Saving and Beam Correspondence

abstract: Beam correspondence, or downlink-uplink (DL-UL) beam reciprocity, refers to the assumption that the best beams in the DL are also the best beams in the UL. This is an important assumption that allows the existing beam management framework in 5G to rely heavily on DL beam sweeping and avoid UL beam sweeping: UL beams are inferred from the measurements of the DL reference signals. Beam correspondence holds when the radio configurations are symmetric in the DL and UL. However, as mmWave technology matures, the DL and the UL face different constraints often breaking the beam correspondence. For example, power constraints may require a UE to activate only a portion of its antenna array for UL transmission, while still activating the full array for DL reception. Meanwhile, if the UL beam with sub-array, named as sub-chain beam in this paper, has a similar radiation pattern as the DL beam, the beam correspondence can still hold. This paper proposes methods for sub-chain beam codebook design to achieve a trade-off between the power saving and beam correspondence.


title: High Throughput Multi-Channel Parallelized Diffraction Convolutional Neural Network Accelerator

abstract: Convolutional neural networks are paramount in image and signal processing including the relevant classification and training tasks alike and constitute for the majority of machine learning compute demand today. With convolution operations being computationally intensive, next generation hardware accelerators need to offer parallelization and algorithmic-hardware homomorphism. Fortunately, diffractive display optics is capable of million-channel parallel data processing at low latency, however, thus far only showed tens of Hertz slow single image and kernel capability, thereby significantly underdelivering from its performance potential. Here, we demonstrate an operation-parallelized high-throughput Fourier optic convolutional neural network accelerator. For the first time simultaneously processing of multiple kernels in Fourier domain enabled by optical diffraction has been achieved alongside with already conventional in the field input parallelism. Additionally, we show an about one hundred times system speed up over existing optical diffraction-based processors and this demonstration rivals performance of modern electronic solutions. Therefore, this system is capable of processing large-scale matrices about ten times faster than state of art electronic systems.


title: Analysis of ECG data to detect Atrial Fibrillation

abstract: Atrial fibrillation(termed as AF/Afib henceforth) is a discrete and often rapid heart rhythm that can lead to clots near the heart. We can detect Afib by ECG signal by the absence of p and inconsistent intervals between R waves as shown in fig(1). Existing methods revolve around CNN that are used to detect afib but most of them work with 12 point lead ECG data where in our case the health gauge watch deals with single-point ECG data. Twelve-point lead ECG data is more accurate than a single point. Furthermore, the health gauge watch data is much noisier. Implementing a model to detect Afib for the watch is a test of how the CNN is changed/modified to work with real life data


title: A Robust Initialization of Residual Blocks for Effective ResNet Training without Batch Normalization

abstract: Batch Normalization is an essential component of all state-of-the-art neural networks architectures. However, since it introduces many practical issues, much recent research has been devoted to designing normalization-free architectures. In this paper, we show that weights initialization is key to train ResNet-like normalization-free networks. In particular, we propose a slight modification to the summation operation of a block output to the skip-connection branch, so that the whole network is correctly initialized. We show that this modified architecture achieves competitive results on CIFAR-10, CIFAR-100 and ImageNet without further regularization nor algorithmic modifications.


title: Maxwell's Demon must remain sebservient to Clausius's statement

abstract: Using classical thermodynamics, we argue that Maxwell's demon loses its battle against Clausius as any temperature difference or other thermodynamic forces it creates is immediately compensated by spontaneous counterbalancing flows that bring about equilibration by slower particles in principle. Being constrained by these spontaneously generated equilibration processes in which he actively but unwittingly participates, the demon is incapable of destroying equilibrium and violating the second law. In fact, our investigation shows that he is unintentionally designed to support it, and does not alter the temperature.


title: New Insights into the Criteria of Fast Radio Burst in the Light of FRB 20121102A

abstract: The total event number of fast radio bursts (FRBs) is accumulating rapidly with the improvement of existing radio telescopes and the completion of new facilities. Especially, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) Collaboration has just reported more than one thousand bursts in a short observing period of 47 days \citep{LiD2021}. The interesting bimodal distribution in their work motivates us to revisit the definition of FRBs. In this work, we ascribe the bimodal distribution to two physical kinds of radio bursts, which may have different radiation mechanisms. We propose to use brightness temperature to separate two subtypes. For FRB 20121102A, the critical brightness temperature is $T_{\rm B,cri}\simeq10^{33}\,\rm K$. Bursts with $T_{\rm B}\geq T_{\rm B,cri}$ are denoted as "classical" FRBs, and further we find a tight pulse width-fluence relation ($T\propto\mathcal{F_\nu}^{0.306}$) for them. On the contrary, the other bursts are considered as "atypical" bursts that may originate from a different physical process. We suggest that for each FRB event, a similar dividing line should exist but $T_{\rm B,cri}$ is not necessarily the same. Its exact value depends on FRB radiation mechanism and properties of the source.


title: Coherent reaction between molecular and atomic Bose-Einstein condensates: integrable model

abstract: We solve a model that describes a stimulated conversion between ultracold bosonic atoms and molecules. The reaction is triggered by a linearly time-dependent transition throughout the Feshbach resonance. Our solution predicts a nonexponential dependence, with a dynamic phase transition, of the reaction efficiency on the transition rate. We find that the emerging phase can have a thermalized energy distribution with the temperature defined by the rate of the transition. This phase, however, has strong purely quantum correlations.


title: Learning with Proper Partial Labels

abstract: Partial-label learning is a kind of weakly-supervised learning with inexact labels, where for each training example, we are given a set of candidate labels instead of only one true label. Recently, various approaches on partial-label learning have been proposed under different generation models of candidate label sets. However, these methods require relatively strong distributional assumptions on the generation models. When the assumptions do not hold, the performance of the methods is not guaranteed theoretically. In this paper, we propose the notion of properness on partial labels. We show that this proper partial-label learning framework requires a weaker distributional assumption and includes many previous partial-label learning settings as special cases. We then derive a unified unbiased estimator of the classification risk. We prove that our estimator is risk-consistent, and we also establish an estimation error bound. Finally, we validate the effectiveness of our algorithm through experiments.


title: Analysis and modelling of macroscopic and microscopic dynamics of a pedestrian cross-flow

abstract: In this work we investigate the behaviour of a human crowd in a cross-flow. We first analyse the results of a set of controlled experiments in which subjects were divided into two groups, in such a way to explore different density settings, and asked to walk through the crossing area. We study the results of the experiment by analysing, along with traditional indicators such as density and velocity, also walking and body orientation, studying how these microscopic observables are influenced by density. Furthermore, we report a preliminary but quantitative analysis on the emergence of self-organising stripes in the crossing area. We also try to reproduce the empirical results using a hierarchy of models, which differ in the details of the body shape (using a disk-shaped body vs a more realistic elliptical shape) and in how collision avoiding is performed (using only information regarding "centre of mass" distance and velocity, or actually introducing body shape information). We verified that the most detailed model (i.e., using body shape information and an elliptical body) outperforms in a significant way the simplest one (using only centre of mass distance and velocity, and disk-shaped bodies). Furthermore, we observed that if elliptical bodies are introduced without using such information in collision avoidance, the performance of the model is relatively poor. Nevertheless, the difference between the different models is relevant only in describing the "tails" of the observable distributions, suggesting that the more complex models could be of practical use only for describing high density settings. We also verified that "stripe formation" emerges in all models.


title: Regular Sequences On Squares Of Monomial Ideals

abstract: We use initially regular sequences that consist of linear sums to explore the depth of $R/I^2$, when $I$ is a monomial ideal in a polynomial ring $R$. We give conditions under which these linear sums form regular or initially regular sequences on $R/I^2$. We then obtain a criterion for when $\depth R/I^2>1$ and a lower bound on $\depth R/I^2$.


title: Selective Multiple Power Iteration: from Tensor PCA to gradient-based exploration of landscapes

abstract: We propose Selective Multiple Power Iterations (SMPI), a new algorithm to address the important Tensor PCA problem that consists in recovering a spike $\bf{v_0}^{\otimes k}$ corrupted by a Gaussian noise tensor $\bf{Z} \in (\mathbb{R}^n)^{\otimes k}$ such that $\bf{T}=\sqrt{n} \beta \bf{v_0}^{\otimes k} + \bf{Z}$ where $\beta$ is the signal-to-noise ratio (SNR). SMPI consists in generating a polynomial number of random initializations, performing a polynomial number of symmetrized tensor power iterations on each initialization, then selecting the one that maximizes $\langle \bf{T}, \bf{v}^{\otimes k} \rangle$. Various numerical simulations for $k=3$ in the conventionally considered range $n \leq 1000$ show that the experimental performances of SMPI improve drastically upon existent algorithms and becomes comparable to the theoretical optimal recovery. We show that these unexpected performances are due to a powerful mechanism in which the noise plays a key role for the signal recovery and that takes place at low $\beta$. Furthermore, this mechanism results from five essential features of SMPI that distinguish it from previous algorithms based on power iteration. These remarkable results may have strong impact on both practical and theoretical applications of Tensor PCA. (i) We provide a variant of this algorithm to tackle low-rank CP tensor decomposition. These proposed algorithms also outperforms existent methods even on real data which shows a huge potential impact for practical applications. (ii) We present new theoretical insights on the behavior of SMPI and gradient descent methods for the optimization in high-dimensional non-convex landscapes that are present in various machine learning problems. (iii) We expect that these results may help the discussion concerning the existence of the conjectured statistical-algorithmic gap.


title: On nonlinear transformations in quantum computation

abstract: While quantum computers are naturally well-suited to implementing linear operations, it is less clear how to implement nonlinear operations on quantum computers. However, nonlinear subroutines may prove key to a range of applications of quantum computing from solving nonlinear equations to data processing and quantum machine learning. Here we develop a series of basic subroutines for implementing nonlinear transformations of input quantum states. Our algorithms are framed around the concept of a weighted state, a mathematical entity describing the output of an operational procedure involving both quantum circuits and classical post-processing.


title: Electric-field-induced parametric excitation of exchange magnons in a CoFeB/MgO junction

abstract: Electric-field controlled magnetization dynamics is an important integrant in low power spintronic devices. In this letter, we demonstrate electric-field induced parametric excitation for CoFeB/MgO junctions by using interfacial in-plane magnetic anisotropy. When the in-plane magnetic anisotropy and the external magnetic field are parallel to each other, magnons are efficiently excited by using electric-field induced parametric excitation. Its wavelength and wavenumber are tuned by changing input power and frequency of the applied voltage. A generalized phenomenological model is developed to explain the underlying role of the electric-field torque. Electrical excitation with no Joule heating offers a good opportunity for developing magnonic devices and exploring various nonlinear dynamics in magnetic systems.


title: Warm inflation with bulk viscous pressure for different solutions of an anisotropic universe

abstract: We study a warm inflationary model for different expansions assuming an anisotropic universe described by Bianchi I metric. The universe is filled with a scalar field or inflaton, radiation, and bulk viscous pressure. We carry out the inflationary analysis for different solutions of such universe in two different cases of the bulk viscosity coefficient $\xi$ and the dissipation coefficient $\Gamma$ as constant and variable parameters, respectively. We compare the obtained results with the recent observations, in order to find the observational constraints on the parameters space of the models. Moreover, we attempt to present a better judgment among the considered models by calculation of the non-linear parameter $f_{NL}$ describing the non-Gaussianity property of the models. Additionally, we investigate the warm inflationary models with viscous pressure from the Weak Gravity Conjecture approach, considering the swampland criteria.


title: Adversarial Attacks against Windows PE Malware Detection: A Survey of the State-of-the-Art

abstract: Malware has been one of the most damaging threats to computers that span across multiple operating systems and various file formats. To defend against ever-increasing and ever-evolving malware, tremendous efforts have been made to propose a variety of malware detection that attempt to effectively and efficiently detect malware so as to mitigate possible damages as early as possible. Recent studies have shown that, on the one hand, existing ML and DL techniques enable superior solutions in detecting newly emerging and previously unseen malware. However, on the other hand, ML and DL models are inherently vulnerable to adversarial attacks in the form of adversarial examples. In this paper, we focus on malware with the file format of portable executable (PE) in the family of Windows operating systems, namely Windows PE malware, as a representative case to study the adversarial attack methods in such adversarial settings. To be specific, we start by first outlining the general learning framework of Windows PE malware detection based on ML/DL and subsequently highlighting three unique challenges of performing adversarial attacks in the context of Windows PE malware. Then, we conduct a comprehensive and systematic review to categorize the state-of-the-art adversarial attacks against PE malware detection, as well as corresponding defenses to increase the robustness of Windows PE malware detection. Finally, we conclude the paper by first presenting other related attacks against Windows PE malware detection beyond the adversarial attacks and then shedding light on future research directions and opportunities. In addition, a curated resource list of adversarial attacks and defenses for Windows PE malware detection is also available at https://github.com/ryderling/adversarial-attacks-and-defenses-for-windows-pe-malware-detection.


title: Measuring the Non-Axially-Symmetric Surface Temperature Distribution of the Central Compact Object in Puppis A

abstract: The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822$-$4300, the CCO in the Puppis A supernova remnant (SNR), using $471$ ks of XMM-Newton data. We compute the energy-dependent pulse profiles in sixteen energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: $kT_{\rm warm} = (1+z) \times 0.222_{-0.019}^{+0.018}$ keV and $kT_{\rm hot} = (1+z) \times 0.411\pm0.011$ keV. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a minimum value of $1.\!^{\circ}1 \pm 0.\!^{\circ}2$, and an expectation value of $9.\!^{\circ}35 \pm 0.\!^{\circ}17$ among the most probable geometries. The discovery of this asymmetry, along with the factor of $\approx2$ temperature difference between the two emitting regions, may indicate that RX J0822$-$4300 was born with a strong, tangled crustal magnetic field.


title: Massive On-shell Recursion Relations for n-point Amplitudes

abstract: We construct two and three-line shifts for tree-level amplitude with massless and/or massive particles, and provide a method to construct general multi-line shifts for all masses. We choose the massless-massive BCFW shift from these shifts and examine its validity in renormalizable theories. Using such a shift, we find that amplitudes with at least one massless vector boson are constructible. This reveals the importance of gauge theory in the construction of amplitudes with massive particles. We also find that this kind of amplitudes have a cancellation related to group structure among different channels, which is essential for constructibility. Furthermore, we show that in the limit of large shift parameter $z$, the amplitude with four massive vector bosons, which can include transverse massive vector particles, have structures proportional to the amplitude with shifted vector particles replaced by Goldstone bosons in the leading order. This is responsible for the failure of massive-massive BCFW recursion relations in the amplitudes with four massive vector bosons.


title: Mean field game for modeling of COVID-19 spread

abstract: The paper presents the one of possible approach to model the epidemic propagation. The proposed model is based on the mean-field control inside separate groups of population, namely, suspectable (S), infected (I), removed (R) and cross-immune (C). In the paper the numerical algorithm to solve such a problem is presented, which ensures the conservation the total mass of population during timeline. Numerical experiments demonstrate the result of modelling the propagation of COVID-19 virus during two 100 day periods in Novosibirsk (Russia).


title: The Lichtenbaum conjecture for abelian extension of imaginary quadratic fields

abstract: In this paper, we prove the cohomological Lichtenbaum conjecture of abelian extensions of imaginary quadratic fields up to a finite set of bad primes.


title: A Multi-Objective Degree-Based Network Anonymization Approach

abstract: Enormous amounts of data collected from social networks or other online platforms are being published for the sake of statistics, marketing, and research, among other objectives. The consequent privacy and data security concerns have motivated the work on degree-based data anonymization. We propose and study a new multi-objective anonymization approach that generalizes the known degree anonymization problem and attempts at improving it as a more realistic model for data security/privacy. Our suggested model guarantees a convenient privacy level based on modifying the degrees in a way that respects some given local restrictions, per node, such that the total modifications at the global level (in the whole graph/network) are bounded by some given value. The corresponding multi-objective graph realization approach is solved using Integer Linear Programming to obtain the best possible solutions. Our experimental studies provide empirical evidence of the effectiveness of the new approach; by specifically showing that the introduced anonymization algorithm has a negligible effect on the way nodes are clustered, thereby preserving valuable network information while significantly improving data privacy.


title: Signed and Unsigned Partial Information Decompositions of Continuous Network Interactions

abstract: We investigate the partial information decomposition (PID) framework as a tool for edge nomination. We consider both the $I_{\cap}^{\text{min}}$ and $I_{\cap}^{\text{PM}}$ PIDs, from arXiv:1004.2515 and arXiv:1801.09010 respectively, and we both numerically and analytically investigate the utility of these frameworks for discovering significant edge interactions. In the course of our work, we extend both the $I_{\cap}^{\text{min}}$ and $I_{\cap}^{\text{PM}}$ PIDs to a general class of continuous trivariate systems. Moreover, we examine how each PID apportions information into redundant, synergistic, and unique information atoms within the source-bivariate PID framework. Both our simulation experiments and analytic inquiry indicate that the atoms of the $I_{\cap}^{\text{PM}}$ PID have a non-specific sensitivity to high predictor-target mutual information, regardless of whether or not the predictors are truly interacting. By contrast, the $I_{\cap}^{\text{min}}$ PID is quite specific, although simulations suggest that it lacks sensitivity.


title: Asymptotic analysis of ML-covariance parameter estimators based on covariance approximations

abstract: Given a zero-mean Gaussian random field with a covariance function that belongs to a parametric family of covariance functions, we introduce a new notion of likelihood approximations, termed truncated-likelihood functions. Truncated-likelihood functions are based on direct functional approximations of the presumed family of covariance functions. For compactly supported covariance functions, within an increasing-domain asymptotic framework, we provide sufficient conditions under which consistency and asymptotic normality of estimators based on truncated-likelihood functions are preserved. We apply our result to the family of generalized Wendland covariance functions and discuss several examples of Wendland approximations. For families of covariance functions that are not compactly supported, we combine our results with the covariance tapering approach and show that ML estimators, based on truncated-tapered likelihood functions, asymptotically minimize the Kullback-Leibler divergence, when the taper range is fixed.


title: Investigating Effect of Dialogue History in Multilingual Task Oriented Dialogue Systems

abstract: While the English virtual assistants have achieved exciting performance with an enormous amount of training resources, the needs of non-English-speakers have not been satisfied well. Up to Dec 2021, Alexa, one of the most popular smart speakers around the world, is able to support 9 different languages [1], while there are thousands of languages in the world, 91 of which are spoken by more than 10 million people according to statistics published in 2019 [2]. However, training a virtual assistant in other languages than English is often more difficult, especially for those low-resource languages. The lack of high-quality training data restricts the performance of models, resulting in poor user satisfaction. Therefore, we devise an efficient and effective training solution for multilingual task-orientated dialogue systems, using the same dataset generation pipeline and end-to-end dialogue system architecture as BiToD[5], which adopted some key design choices for a minimalistic natural language design where formal dialogue states are used in place of natural language inputs. This reduces the room for error brought by weaker natural language models, and ensures the model can correctly extract the essential slot values needed to perform dialogue state tracking (DST). Our goal is to reduce the amount of natural language encoded at each turn, and the key parameter we investigate is the number of turns (H) to feed as history to model. We first explore the turning point where increasing H begins to yield limiting returns on the overall performance. Then we examine whether the examples a model with small H gets wrong can be categorized in a way for the model to do few-shot finetuning on. Lastly, will explore the limitations of this approach, and whether there is a certain type of examples that this approach will not be able to resolve.


title: A One-Dimensional Model for Star Formation Near a "Leaky" Black Hole

abstract: In the presence of a Weyl scaling invariant cosmological action, black holes no longer have an event horizon and an apparent horizon. So they have no trapped surfaces, and may be "leaky", emitting a "black hole wind" which could lead to star formation in the neighborhood of the hole. In this paper we formulate and analyze a one-dimensional model for star formation resulting from a postulated outgoing black hole particle flux, incident on a distant spherical surface modeled as a set of planar disks surrounding the hole. Using the Toomre analysis of the Jeans instability of a disk, we compute conditions for a disk collapse instability and estimate the collapse time. We suggest a mechanism for giving the disk angular momentum around the central black hole. This gives a possible explanation for the puzzling observation of young stars forming in the vicinity of the black hole Sagittarius A* central to the Milky Way galaxy.


title: Model Selection in Batch Policy Optimization

abstract: We study the problem of model selection in batch policy optimization: given a fixed, partial-feedback dataset and $M$ model classes, learn a policy with performance that is competitive with the policy derived from the best model class. We formalize the problem in the contextual bandit setting with linear model classes by identifying three sources of error that any model selection algorithm should optimally trade-off in order to be competitive: (1) approximation error, (2) statistical complexity, and (3) coverage. The first two sources are common in model selection for supervised learning, where optimally trading-off these properties is well-studied. In contrast, the third source is unique to batch policy optimization and is due to dataset shift inherent to the setting. We first show that no batch policy optimization algorithm can achieve a guarantee addressing all three simultaneously, revealing a stark contrast between difficulties in batch policy optimization and the positive results available in supervised learning. Despite this negative result, we show that relaxing any one of the three error sources enables the design of algorithms achieving near-oracle inequalities for the remaining two. We conclude with experiments demonstrating the efficacy of these algorithms.


title: Physics Constrained Flow Neural Network for Short-Timescale Predictions in Data Communications Networks

abstract: Machine learning is gaining growing momentum in various recent models for the dynamic analysis of information flows in data communications networks. These preliminary models often rely on off-the-shelf learning models to predict from historical statistics while disregarding the physics governing the generating behaviors of these flows. This paper instead introduces Flow Neural Network (FlowNN) to improve the feature representation with learned physical bias. This is implemented by an induction layer, working upon the embedding layer, to impose the physics connected data correlations, and a self-supervised learning strategy with stop-gradient to make the learned physics universal. For the short-timescale network prediction tasks, FlowNN achieves 17% - 71% of loss decrease than the state-of-the-art baselines on both synthetic and real-world networking datasets, which shows the strength of this new approach.


title: High-order tensor flow processing using integrated photonic circuits

abstract: Tensor analytics lays mathematical basis for the prosperous promotion of multiway signal processing. To increase computing throughput, mainstream processors transform tensor convolutions to matrix multiplications to enhance parallelism of computing. However, such order-reducing transformation produces data duplicates and consumes additional memory. Here, we demonstrate an integrated photonic tensor flow processor without tensor-matrix transformation, which outputs the convolved tensor as the input tensor 'flows' through the processor. The hybrid manipulation of optical dimensions of wavelength, time, and space enables the direct representation and processing of high-order tensors in optical domain. In the proof-of-concept experiment, processing of multi-channel images and videos is accomplished at the frequency of 20 GHz. A convolutional neural network is demonstrated on the processor, which achieves an accuracy of 97.9 percent on action recognition.


title: How do the finite-size particles modify the drag in Taylor-Couette turbulent flow

abstract: We experimentally investigate the drag modification by neutrally buoyant finite-size particles with various aspect ratios in a Taylor-Couette (TC) turbulent flow. The current Reynolds number, $Re$, ranges from $6.5\times10^3$ to $2.6\times10^4$, and the particle volume fraction, $\Phi$, is up to $10\%$. Particles with three kinds of aspect ratio, $\lambda$, are used: $\lambda=1/3$ (oblate), $\lambda=1$ (spherical) and $\lambda=3$ (prolate). Unlike the case of bubbly TC flow, we find that the suspended finite-size particles increase the drag of the TC system regardless of their aspect ratios. In addition, the normalized friction coefficient, $c_{f,\Phi}/c_{f,\Phi=0}$, decreases with increasing $Re$, the reason could be that in the current low volume fractions the turbulent stress becomes dominant at higher $Re$. As $Re$ increases, the particles distribute more evenly in the entire system, which results from both the greater turbulence intensity and the more pronounced finite-size effects of the particles at higher $Re$. Moreover, it is found that the variation of the particle aspect ratios leads to different particle collective effects. The suspended spherical particles, which tend to cluster near the walls and form a particle layer, significantly affect the boundary layer and result in maximum drag modification. The minimal drag modification is found in the oblate case, where the particles preferentially cluster in the bulk region, and thus the particle layer is absent. Based on the optical measurement results, it can be concluded that, in the low volume fraction ranges ($\Phi=0.5\%$ and $\Phi = 2\%$ here), the larger drag modification is connected to the near-wall particle clustering. The present findings suggest that the particle shape plays a significant role in drag modification.


title: Second-Order Real Nodal-Line Semimetal in Three-Dimensional Graphdiyne

abstract: Real topological phases featuring real Chern numbers and second-order boundary modes have been a focus of current research, but finding their material realization remains a challenge. Here, based on first-principles calculations and theoretical analysis, we reveal the already experimentally synthesized three-dimensional (3D) graphdiyne as the first realistic example of the recently proposed second-order real nodal-line semimetal. We show that the material hosts a pair of real nodal rings, each protected by two topological charges: a real Chern number and a 1D winding number. The two charges generate distinct topological boundary modes at distinct boundaries. The real Chern number leads to a pair of hinge Fermi arcs, whereas the winding number protects a double drumhead surface bands. We develop a low-energy model for 3D graphdiyne which captures the essential topological physics. Experimental aspects and possible topological transition to a 3D real Chern insulator phase are discussed.


title: Globally convergent visual-feature range estimation with biased inertial measurements

abstract: The design of a globally convergent position observer for feature points from visual information is a challenging problem, especially for the case with only inertial measurements and without assumptions of uniform observability, which remained open for a long time. We give a solution to the problem in this paper assuming that only the bearing of a feature point, and biased linear acceleration and rotational velocity of a robot -- all in the body-fixed frame -- are available. Further, in contrast to existing related results, we do not need the value of the gravitational constant either. The proposed approach builds upon the parameter estimation-based observer recently developed in (Ortega et al., Syst. Control Lett., vol.85, 2015) and its extension to matrix Lie groups in our previous work. Conditions on the robot trajectory under which the observer converges are given, and these are strictly weaker than the standard persistency of excitation and uniform complete observability conditions. Finally, as an illustration, we apply the proposed design to the visual inertial navigation problem.


title: Age of Information in Energy Harvesting Aided Massive Multiple Access Networks

abstract: Given the proliferation of the massive machine type communication devices (MTCDs) in beyond 5G (B5G) wireless networks, energy harvesting (EH) aided next generation multiple access (NGMA) systems have drawn substantial attention in the context of energy-efficient data sensing and transmission. However, without adaptive time slot (TS) and power allocation schemes, NGMA systems relying on stochastic sampling instants might lead to tardy actions associated both with high age of information (AoI) as well as high power consumption. For mitigating the energy consumption, we exploit a pair of sleep-scheduling policies, namely the multiple vacation (MV) policy and start-up threshold (ST) policy, which are characterized in the context of three typical multiple access protocols, including time-division multiple access (TDMA), frequency-division multiple access (FDMA) and non-orthogonal multiple access (NOMA). Furthermore, we derive closed-form expressions for the MTCD system's peak AoI, which are formulated as the optimization objective under the constraints of EH power, status update rate and stability conditions. An exact linear search based algorithm is proposed for finding the optimal solution by fixing the status update rate. As a design alternative, a low complexity concave-convex procedure (CCP) is also formulated for finding a near-optimal solution relying on the original problem's transformation into a form represented by the difference of two convex problems. Our simulation results show that the proposed algorithms are beneficial in terms of yielding a lower peak AoI at a low power consumption in the context of the multiple access protocols considered.


title: Making sense of electrical vehicle discussions using sentiment analysis on closely related news and user comments

abstract: We used a token-wise and document-wise sentiment analysis using both unsupervised and supervised models applied to both news and user reviews dataset. And our token-wise sentiment analysis found a statistically significant difference in sentiment between the two groups (both of which were very large N), our document-wise supervised sentiment analysis found no significant difference in sentiment.


title: Robust and Precise Facial Landmark Detection by Self-Calibrated Pose Attention Network

abstract: Current fully-supervised facial landmark detection methods have progressed rapidly and achieved remarkable performance. However, they still suffer when coping with faces under large poses and heavy occlusions for inaccurate facial shape constraints and insufficient labeled training samples. In this paper, we propose a semi-supervised framework, i.e., a Self-Calibrated Pose Attention Network (SCPAN) to achieve more robust and precise facial landmark detection in challenging scenarios. To be specific, a Boundary-Aware Landmark Intensity (BALI) field is proposed to model more effective facial shape constraints by fusing boundary and landmark intensity field information. Moreover, a Self-Calibrated Pose Attention (SCPA) model is designed to provide a self-learned objective function that enforces intermediate supervision without label information by introducing a self-calibrated mechanism and a pose attention mask. We show that by integrating the BALI fields and SCPA model into a novel self-calibrated pose attention network, more facial prior knowledge can be learned and the detection accuracy and robustness of our method for faces with large poses and heavy occlusions have been improved. The experimental results obtained for challenging benchmark datasets demonstrate that our approach outperforms state-of-the-art methods in the literature.


title: MVDG: A Unified Multi-view Framework for Domain Generalization

abstract: To generalize the model trained in source domains to unseen target domains, domain generalization (DG) has recently attracted lots of attention. Since target domains can not be involved in training, overfitting source domains is inevitable. As a popular regularization technique, the meta-learning training scheme has shown its ability to resist overfitting. However, in the training stage, current meta-learning-based methods utilize only one task along a single optimization trajectory, which might produce a biased and noisy optimization direction. Beyond the training stage, overfitting could also cause unstable prediction in the test stage. In this paper, we propose a novel multi-view DG framework to effectively reduce the overfitting in both the training and test stage. Specifically, in the training stage, we develop a multi-view regularized meta-learning algorithm that employs multiple optimization trajectories to produce a suitable optimization direction for model updating. We also theoretically show that the generalization bound could be reduced by increasing the number of tasks in each trajectory. In the test stage, we utilize multiple augmented images to yield a multi-view prediction to alleviate unstable prediction, which significantly promotes model reliability. Extensive experiments on three benchmark datasets validate that our method can find a flat minimum to enhance generalization and outperform several state-of-the-art approaches.


title: Precision identified-hadron spectrum analysis for 5 TeV $\bf p$-$\bf Pb$ collisions -- Part II

abstract: This is the second part of a two-part article on precision modeling of identified-hadron (PID) $p_t$ spectra from 5 TeV $p$-Pb collisions. In Part I a revised two-component (soft + hard) model (TCM) of PID spectra was introduced, an apparent detection inefficiency for protons was corrected and jet-related spectrum hard components isolated via an improved method were compared to a fixed TCM serving as reference. In the present article (Part II) the TCM is further elaborated to describe centrality variation of spectrum hard components (and therefore entire spectra) within data statistical uncertainties as determined by a standard statistical measure (Z-scores). The completed PID TCM is then used to investigate properties of spectrum and yield ratios (e.g. $p/\pi$ ratios) and ensemble-mean $\bar p_t$ data. Systematic differences are observed between meson and baryon hard components. With increasing $p$-Pb centrality meson hard components shift to lower $p_t$ while baryon hard components shift to higher $p_t$. Those trends explain quantitatively the main features and centrality variation of spectrum ratios. Ensemble-mean $\bar p_t$ trends are also predicted at the level of data statistical uncertainties, and details of the $\bar p_t$ centrality trends are explained quantitatively.


title: Flakify: A Black-Box, Language Model-based Predictor for Flaky Tests

abstract: Software testing assures that code changes do not adversely affect existing functionality. However, a test case can be flaky, i.e., passing and failing across executions, even for the same version of the source code. Flaky test cases introduce overhead to software development as they can lead to unnecessary attempts to debug production or testing code. The state-of-the-art ML-based flaky test case predictors rely on pre-defined sets of features that are either project-specific, require access to production code, which is not always available to software test engineers. Therefore, in this paper, we propose Flakify, a black-box, language model-based predictor for flaky test cases. Flakify relies exclusively on the source code of test cases, thus not requiring to (a) access to production code (black-box), (b) rerun test cases, (c) pre-define features. To this end, we employed CodeBERT, a pre-trained language model, and fine-tuned it to predict flaky test cases using the source code of test cases. We evaluated Flakify on two publicly available datasets (FlakeFlagger and IDoFT) for flaky test cases and compared our technique with the FlakeFlagger approach using two different evaluation procedures: cross-validation and per-project validation. Flakify achieved high F1-scores on both datasets using cross-validation and per-project validation, and surpassed FlakeFlagger by 10 and 18 percentage points in terms of precision and recall, respectively, when evaluated on the FlakeFlagger dataset, thus reducing the cost bound to be wasted on unnecessarily debugging test cases and production code by the same percentages. Flakify also achieved significantly higher prediction results when used to predict test cases on new projects, suggesting better generalizability over FlakeFlagger. Our results further show that a black-box version of FlakeFlagger is not a viable option for predicting flaky test cases.


title: Two-mode light states before and after delocalized single-photon addition

abstract: We studied the effect of delocalized single-photon addition (DPA) on two input modes containing four cases: two independent coherent states (CSs), two independent thermal states (TSs), two independent single-mode squeezed vacuums (SVs), and an entangled two-mode squeezed vacuum (TMSV). In essence, four types of new non-Gaussian entangled light states are generated. We studied three different resources (including entanglement, discorrelation and Wigner negativity) for each two-mode light state. The output states after DPA are entangled, with more parameters and complex structures, characterizing more Wigner negativity or even discorrelation. In contrast, the CSs case is the most tunable protocol, because its negativity under partial transposition, discorrelation, and Wigner logarithmic negativity are more sensitive to superposition phase than those in TSs, SVs and TMSV cases.


title: Asymptotic normality of least squares estimators to stochastic differential equations driven by fractional Brownian motions

abstract: We will consider the following stochastic differential equation (SDE): \begin{equation} X_t=X_0+\int_0^tb(X_s,\theta_0)ds+\sigma B_t,~~~t\in(0,T], \end{equation} where $\{B_t\}_{t\ge 0}$ is a fractional Brownian motion with Hurst index $H\in(1/2,1)$, $\theta_0$ is a parameter that contains a bounded and open convex subset $\Theta\subset\mathbb{R}^d$, $\{b(\cdot,\theta),\theta\in\Theta\}$ is a family of drift coefficients with $b(\cdot,\theta):\mathbb{R}\rightarrow\mathbb{R}$, and $\sigma\in\mathbb{R}$ is assumed to be the known diffusion coefficient.


title: The Color Gradients of the Globular Cluster Systems in M87 and M49

abstract: Combining data from the ACS Virgo Cluster Survey (ACSVCS) and the Next Generation Virgo cluster Survey (NGVS), we extend previous studies of color gradients of the globular cluster (GC) systems of the two most massive galaxies in the Virgo cluster, M87 and M49, to radii of $\sim 15~R_e$ ($\sim 200$ kpc for M87 and $\sim 250$ kpc for M49). We find significant negative color gradients, i.e., becoming bluer with increasing distance, out to these large radii. The gradients are driven mainly by the outwards decrease of the ratio of red to blue GC numbers. The color gradients are also detected out to $\sim 15~R_e$ in the red and blue sub-populations of GCs taken separately. In addition, we find a negative color gradient when we consider the satellite low-mass elliptical galaxies as a system, i.e., the satellite galaxies closer to the center of the host galaxy usually have redder color indices, both for their stars and GCs. According to the "two phase" formation scenario of massive early-type galaxies, the host galaxy accretes stars and GCs from low-mass satellite galaxies in the second phase. So the accreted GC system naturally inherits the negative color gradient present in the satellite population. This can explain why the color gradient of the GC system can still be observed at large radii after multiple minor mergers.


title: Primordial black hole formation from massless scalar isocurvature

abstract: We numerically study the primordial black hole (PBH) formation by an isocurvature perturbation of a massless scalar field on super Hubble scales in the radiation-dominated universe. As a first step we perform simulations of spherically symmetric configurations. For the initial condition, we employ the spatial gradient expansion and provide the general form of the growing mode solutions valid up through the second order in this expansion. The initial scalar field profile is assumed to be Gaussian with a characteristic comoving wavenumber $k$; $\sim\exp(-k^2R^2)$, where $R$ is the radial coordinate. We find that a PBH is formed for a sufficiently large amplitude of the scalar field profile. Nevertheless, we find that the late time behavior of the gravitational collapse is dominated by the dynamics of the fluid but not by the scalar field, which is analogous to the PBH formation from an adiabatic perturbation in the radiation-dominated universe.


title: Deterministic generation of a perfect soliton crystal with a saturable absorber

abstract: We numerically investigate the deterministic generation of a perfect soliton crystal (PSC) in an optical microresonator functionalized with a saturable absorber (SA). The SA allows the direct formation of a PSC from an initial, periodic Turing roll. It prevents passage through a chaotic state, which induces a stochastic nature as regards the number of generated dissipative Kerr solitons. We show that PSCs form deterministically, and the number is controlled by adjusting the input power and SA parameter. Our work provides a simple approach for obtaining a stable PSC that offers an ultra-high repetition rate and a high comb output power.