title: A relation of thermodynamic relevance between the superadditivity, concavity and homogeneity properties of real-valued functions

abstract: We provide a necessary and sufficient condition for the validity of the following Landsberg-Thirring theorem: for a real-valued function on a convex set, any two of the properties of superadditivity, concavity and homogeneity implies the third. Applications to statistical thermodynamics, following Thirring and Landsberg, are briefly reviewed.


title: Controllability of LTI Networked Systems with Heterogeneous Dynamics

abstract: In this paper, a necessary and sufficient condition for the controllability of networked systems with heterogeneous dynamics is established where the nodes are higher dimensional linear time invariant systems and the network topology is directed and weighted. The controllability of networked system over some specific topologies are also examined and some non-controllability results are obtained. The theoretical results are demonstrated with examples.


title: Sharp conditions for the linearization of finite elasticity

abstract: We consider the topic of linearization of finite elasticity for pure traction problems. We characterize the variational limit for the approximating sequence of rescaled nonlinear elastic energies. We show that the limiting minimal value can be strictly lower than the minimal value of the standard linear elastic energy if a strict compatibility condition for external loads does not hold. The results are provided for both the compressible and the incompressible case.


title: Planar narrow-band-pass filter based on Si resonant metasurface

abstract: Optically resonant dielectric metasurfaces offer unique capability to fully control the wavefront, polarisation, intensity or spectral content of light based on the excitation and interference of different electric and magnetic Mie multipolar resonances. Recent advances of the wide accessibility in the nanofabrication and nanotechnologies have led to a surge in the research field of high-quality functional optical metasurfaces which can potentially replace or even outperform conventional optical components with ultra-thin feature. Replacing conventional optical filtering components with metasurface technology offers remarkable advantages including lower integration cost, ultra-thin compact configuration, easy combination with multiple functions and less restriction on materials. Here we propose and experimentally demonstrate a planar narrow-band-pass filter based on the optical dielectric metasurface composed of Si nanoresonators in array. A broadband transmission spectral valley (around 200~nm) has been realised by combining electric and magnetic dipole resonances adjacent to each other. Meanwhile, we obtain a narrow-band transmission peak by exciting a high-quality leaky mode which is formed by partially breaking a bound state in the continuum generated by the collective longitudinal magnetic dipole resonances in the metasurface. Our proposed metasurface-based filter shows a stable performance for oblique light incidence with small angles (within 10 deg). Our work imply many potential applications of nanoscale photonics devices such as displays, spectroscopy, etc.


title: FraCaS: Temporal Analysis

abstract: In this paper, we propose an implementation of temporal semantics which is suitable for inference problems. This implementation translates syntax trees to logical formulas, suitable for consumption by the Coq proof assistant. We support several phenomena including: temporal references, temporal adverbs, aspectual classes and progressives. We apply these semantics to the complete FraCaS testsuite. We obtain an accuracy of 81 percent overall and 73 percent for problems explicitly marked as related to temporal reference.


title: Testing General Relativity with NuSTAR data of Galactic Black Holes

abstract: Einstein's theory of General Relativity predicts that the spacetime metric around astrophysical black holes is described by the Kerr solution. In this work, we employ state-of-the-art in relativistic reflection modeling to analyze a selected set of NuSTAR spectra of Galactic black holes to obtain the most robust and precise constraints on the Kerr black hole hypothesis possible today. Our constraints are much more stringent than those from other electromagnetic techniques, and with some sources we get stronger constraints than those currently available from gravitational waves.


title: Asymptotic behavior of Ext for pairs of modules of large complexity over graded complete intersections

abstract: Let $M$ and $N$ be finitely generated graded modules over a graded complete intersection $R$ such that $\operatorname{Ext}_R^i(M,N)$ has finite length for all $i\gg 0$. We show that the even and odd Hilbert polynomials, which give the lengths of $\operatorname{Ext}^i_R(M,N)$ for all large even $i$ and all large odd $i$, have the same degree and leading coefficient whenever the highest degree of these polynomials is at least the dimension of $M$ or $N$. Refinements of this result are given when $R$ is regular in small codimensions.


title: SMART Frame Selection for Action Recognition

abstract: Action recognition is computationally expensive. In this paper, we address the problem of frame selection to improve the accuracy of action recognition. In particular, we show that selecting good frames helps in action recognition performance even in the trimmed videos domain. Recent work has successfully leveraged frame selection for long, untrimmed videos, where much of the content is not relevant, and easy to discard. In this work, however, we focus on the more standard short, trimmed action recognition problem. We argue that good frame selection can not only reduce the computational cost of action recognition but also increase the accuracy by getting rid of frames that are hard to classify. In contrast to previous work, we propose a method that instead of selecting frames by considering one at a time, considers them jointly. This results in a more efficient selection, where good frames are more effectively distributed over the video, like snapshots that tell a story. We call the proposed frame selection SMART and we test it in combination with different backbone architectures and on multiple benchmarks (Kinetics, Something-something, UCF101). We show that the SMART frame selection consistently improves the accuracy compared to other frame selection strategies while reducing the computational cost by a factor of 4 to 10 times. Additionally, we show that when the primary goal is recognition performance, our selection strategy can improve over recent state-of-the-art models and frame selection strategies on various benchmarks (UCF101, HMDB51, FCVID, and ActivityNet).


title: A Declarative Metamorphic Testing Framework for Autonomous Driving

abstract: Autonomous driving has gained much attention from both industry and academia. Currently, Deep Neural Networks (DNNs) are widely used for perception and control in autonomous driving. However, several fatal accidents caused by autonomous vehicles have raised serious safety concerns about autonomous driving models. Some recent studies have successfully used the metamorphic testing technique to detect thousands of potential issues in some popularly used autonomous driving models. However, prior study is limited to a small set of metamorphic relations, which do not reflect rich, real-world traffic scenarios and are also not customizable. This paper presents a novel declarative rule-based metamorphic testing framework called RMT. RMT provides a rule template with natural language syntax, allowing users to flexibly specify an enriched set of testing scenarios based on real-world traffic rules and domain knowledge. RMT automatically parses human-written rules to metamorphic relations using an NLP-based rule parser referring to an ontology list and generates test cases with a variety of image transformation engines. We evaluated RMT on three autonomous driving models. With an enriched set of metamorphic relations, RMT detected a significant number of abnormal model predictions that were not detected by prior work. Through a large-scale human study on Amazon Mechanical Turk, we further confirmed the authenticity of test cases generated by RMT and the validity of detected abnormal model predictions.


title: Solving Problems with Inconsistent Constraints with a Modified Augmented Lagrangian Method

abstract: We present a numerical method for the minimization of constrained optimization problems where the objective is augmented with large quadratic penalties of inconsistent equality constraints. Such objectives arise from quadratic integral penalty methods for the direct transcription of optimal control problems. The Augmented Lagrangian Method (ALM) has a number of advantages over the Quadratic Penalty Method (QPM). However, if the equality constraints are inconsistent, then ALM might not converge to a point that minimizes the bias of the objective and penalty term. Therefore, we present a modification of ALM that fits our purpose. We prove convergence of the modified method and bound its local convergence rate by that of the unmodified method. Numerical experiments demonstrate that the modified ALM can minimize certain quadratic penalty-augmented functions faster than QPM, whereas the unmodified ALM converges to a minimizer of a significantly different problem.


title: Camera-aware Proxies for Unsupervised Person Re-Identification

abstract: This paper tackles the purely unsupervised person re-identification (Re-ID) problem that requires no annotations. Some previous methods adopt clustering techniques to generate pseudo labels and use the produced labels to train Re-ID models progressively. These methods are relatively simple but effective. However, most clustering-based methods take each cluster as a pseudo identity class, neglecting the large intra-ID variance caused mainly by the change of camera views. To address this issue, we propose to split each single cluster into multiple proxies and each proxy represents the instances coming from the same camera. These camera-aware proxies enable us to deal with large intra-ID variance and generate more reliable pseudo labels for learning. Based on the camera-aware proxies, we design both intra- and inter-camera contrastive learning components for our Re-ID model to effectively learn the ID discrimination ability within and across cameras. Meanwhile, a proxy-balanced sampling strategy is also designed, which facilitates our learning further. Extensive experiments on three large-scale Re-ID datasets show that our proposed approach outperforms most unsupervised methods by a significant margin. Especially, on the challenging MSMT17 dataset, we gain $14.3\%$ Rank-1 and $10.2\%$ mAP improvements when compared to the second place. Code is available at: \texttt{https://github.com/Terminator8758/CAP-master}.


title: Atmospheric Pressure Mass Spectrometry of Single Viruses and Nanoparticles by Nanoelectromechanical Systems

abstract: Mass spectrometry of intact nanoparticles and viruses can serve as a potent characterization tool for material science and biophysics. Inaccessible by widespread commercial techniques, the mass of single nanoparticles and viruses (>10MDa) can be readily measured by NEMS (Nanoelectromechanical Systems) based Mass Spectrometry, where charged and isolated analyte particles are generated by Electrospray Ionization (ESI) in air and transported onto the NEMS resonator for capture and detection. However, the applicability of NEMS as a practical solution is hindered by their miniscule surface area, which results in poor limit-of-detection and low capture efficiency values. Another hindrance is the necessity to house the NEMS inside complex vacuum systems, which is required in part to focus analytes towards the miniscule detection surface of the NEMS. Here, we overcome both limitations by integrating an ion lens onto the NEMS chip. The ion lens is composed of a polymer layer, which charges up by receiving part of the ions incoming from the ESI tip and consequently starts to focus the analytes towards an open window aligned with the active area of the NEMS electrostatically. With this integrated system, we have detected the mass of gold and polystyrene nanoparticles under ambient conditions and with two orders-of-magnitude improvement in capture efficiency compared to the state-of-the-art. We then applied this technology to obtain the mass spectrum of SARS-CoV-2 and BoHV-1 virions. With the increase in analytical throughput, the simplicity of the overall setup and the operation capability under ambient conditions, the technique demonstrates that NEMS Mass Spectrometry can be deployed for mass detection of engineered nanoparticles and biological samples efficiently.


title: On $\ell_\infty$-Grothendieck subspaces

abstract: A closed subspace $S$ of $\ell_\infty$ is said to be a \emph{$\ell_\infty$-Grothendieck subspace} if $c_0\subset S$ (hence $\ell_\infty\subset S^{**}$) and every $\sigma(S^*,S)$-convergent sequence in $S^*$ is $\sigma(S^*,\ell_\infty)$-convergent. Here we give examples of closed subspaces of $\ell_\infty$ containing $c_0$ which are or fail to be $\ell_\infty$-Grothendieck.


title: How does neural activity encode spontaneous motor behavior in zebrafish larvae ?

abstract: The origins of spontaneous movements have been investigated in human as well as in other vertebrates. Studies have reported an increase in neuronal activity one second before the onset of a given movement: this is known as readiness potential. The mechanisms underlying this increase are still unclear. Zebrafish larva is an ideal animal model to study the neuronal basis of spontaneous movements. Because of its small size and transparency, this vertebrate is an ideal candidate to apply optical recording methods. In order to understand what neuronal activity causes the execution of a specific tail movement at a given time, we will mainly use a prediction approach.


title: A multiple-relaxation-time lattice Boltzmann model based four-level finite-difference scheme for one-dimensional diffusion equation

abstract: In this paper, we first present a multiple-relaxation-time lattice Boltzmann (MRT-LB) model for one-dimensional diffusion equation where the D1Q3 (three discrete velocities in one-dimensional space) lattice structure is considered. Then through the theoretical analysis, we derive an explicit four-level finite-difference scheme from this MRT-LB model. The results show that the four-level finite-difference scheme is unconditionally stable, and through adjusting the weight coefficient $\omega_{0}$ and the relaxation parameters $s_1$ and $s_2$ corresponding to the first and second moments, it can also have a sixth-order accuracy in space. Finally, we also test the four-level finite-difference scheme through some numerical simulations, and find that the numerical results are consistent with our theoretical analysis.


title: Intelligent Reflecting Surfaces: Sum-Rate Optimization Based on Statistical CSI

abstract: In this paper, we consider a multi-user multiple-input multiple-output (MIMO) system aided by multiple intelligent reflecting surfaces (IRSs) that are deployed to increase the coverage and, possibly, the rank of the channel. We propose an optimization algorithm to configure the IRSs, which is aimed at maximizing the network sum-rate by exploiting only the statistical characterization of the environment, such as the distribution of the locations of the users and the distribution of the multipath channels. As a consequence, the proposed approach does not require the estimation of the instantaneous channel state information (CSI) for system optimization, thus significantly relaxing (or even avoiding) the need of frequently reconfiguring the IRSs, which constitutes one of the most critical issues in IRS-assisted systems. Numerical results confirm the validity of the proposed approach. It is shown, in particular, that IRS-assisted wireless systems that are optimized based on statistical CSI still provide large performance gains as compared to the baseline scenarios in which no IRSs are deployed.


title: Multi-kT/m Focusing Gradient in a Linear Active Plasma Lens

abstract: Active plasma lenses are compact devices developed as a promising beam-focusing alternative for charged particle beams, capable of short focal lengths for high-energy beams. We have previously shown that linear magnetic fields with gradients of around 0.3 kT/m can be achieved in argon-filled plasma lenses that preserve beam emittance [C.A. Lindstr{\o}m et al., Phys. Rev. Lett. 121, 194801 (2018)]. Here we show that with argon in a 500 {\mu}m diameter capillary, the fields are still linear with a focusing gradient of 3.6 kT/m, which is an order of magnitude higher than the gradients of quadrupole magnets. The current pulses that generate the magnetic field are provided by compact Marx banks, and are highly repeatable. These results establish active plasma lenses as an ideal device for pulsed particle beam applications requiring very high focusing gradients that are uniform throughout the lens aperture.


title: Blockchain-Based Approach for Securing Spectrum Trading in Multibeam Satellite Systems

abstract: This paper presents a blockchain-based approach for securing spectrum sharing in multi-beam satellite systems. Satellite spectrum is a scarce resource that requires highly efficient management schemes for optimized sharing by network users. However, spectrum sharing is vulnerable to attacks by malicious protocol participants. In order to ensure efficient spectrum management in the face of dishonest satellite users or cyber attackers, it is important for spectrum sharing mechanism to provide transparency and traceability of the trading process so as to enable the system to detect, and hence eliminate, unauthorized access by malicious users. We address these requirements by proposing the use of blockchain which, apart from its ability to provide transparency and traceability, ensures an immutable means for keeping track of user trading reputation. Besides, in order to address the practical constraints of heterogeneous user nodes, we also propose the use of edge computing to support users with limited computing power. In this paper, we propose a blockchain-based spectrum trading framework and, based on which, a multibeam satellite spectrum sharing algorithm for interference pricing and heterogeneous spectrum demands is devised to improve the efficiency of satellite spectrum. By leveraging on the system characteristics of blockchain, a dynamic spectrum sharing mechanism with traceability, openness and transparency for whole trading process is presented. Numerical results are also provided to evaluate the system benefits and spectrum pricing of the proposed mechanism.


title: Deep Reinforcement Learning for Joint Spectrum and Power Allocation in Cellular Networks

abstract: A wireless network operator typically divides the radio spectrum it possesses into a number of subbands. In a cellular network those subbands are then reused in many cells. To mitigate co-channel interference, a joint spectrum and power allocation problem is often formulated to maximize a sum-rate objective. The best known algorithms for solving such problems generally require instantaneous global channel state information and a centralized optimizer. In fact those algorithms have not been implemented in practice in large networks with time-varying subbands. Deep reinforcement learning algorithms are promising tools for solving complex resource management problems. A major challenge here is that spectrum allocation involves discrete subband selection, whereas power allocation involves continuous variables. In this paper, a learning framework is proposed to optimize both discrete and continuous decision variables. Specifically, two separate deep reinforcement learning algorithms are designed to be executed and trained simultaneously to maximize a joint objective. Simulation results show that the proposed scheme outperforms both the state-of-the-art fractional programming algorithm and a previous solution based on deep reinforcement learning.


title: Disjointly non-singular operators on Banach lattices

abstract: An operator $T$ from a Banach lattice $E$ into a Banach space is disjointly non-singular ($DN$-$S$, for short) if no restriction of $T$ to a subspace generated by a disjoint sequence is strictly singular. We obtain several results for $DN$-$S$ operators, including a perturbative characterization. For $E=L_p$ ($1< p<\infty$) we improve the results, and we show that the $DN$-$S$ operators have a different behavior in the cases $p=2$ and $p\neq 2$. As an application we prove that the strongly embedded subspaces of $L_p$ form an open subset in the set of all closed subspaces.


title: GPU acceleration of a patient-specific airway image segmentation and its assessment

abstract: Image segmentation plays an important role in computer vision, object detection, traffic control, and video surveillance. Typically, it is a critical step in the 3D reconstruction of a specific organ in medical image processing which unveils the detailed tomography of organ, tumor, and nerve, and thus helping to improve the quality of surgical pathology. However, there may be high computational requirements in it. With the advent of GPUs, more complex and realistic models can be simulated, but the deployment of these facilities also requires a huge amount of capital. As a consequence, how to make good use of these computational resource is essential to GPU computing. This study discusses the image segmentation of 3D airway reconstruction, identifies the computing-intensive task, and parallelizes the algorithm of image segmentation in order to obtain theoretical maximum speedup in terms of the benchmark ratio of GPU to CPU. There are five steps involved, which are the image acquisition, pre-processing, segmentation, reconstruction, and object recognition. It is worth to note that it takes 85\% of time on segmentation. This study successfully accelerates the image segmentation of 3D airway reconstruction by optimizing the memory usage, grid and block setting and multiple GPUs communication, thereby gaining a total speedup of 61.8 on two GPUs (Nvidia K40).


title: Unsupervised Scale-Invariant Multispectral Shape Matching

abstract: Alignment between non-rigid stretchable structures is one of the most challenging tasks in computer vision, as the invariant properties are hard to define, and there is no labeled data for real datasets. We present unsupervised neural network architecture based upon the spectral domain of scale-invariant geometry. We build on top of the functional maps architecture, but show that learning local features, as done until now, is not enough once the isometry assumption breaks. We demonstrate the use of multiple scale-invariant geometries for solving this problem. Our method is agnostic to local-scale deformations and shows superior performance for matching shapes from different domains when compared to existing spectral state-of-the-art solutions.


title: Inference in experiments conditional on observed imbalances in covariates

abstract: Double blind randomized controlled trials are traditionally seen as the gold standard for causal inferences as the difference-in-means estimator is an unbiased estimator of the average treatment effect in the experiment. The fact that this estimator is unbiased over all possible randomizations does not, however, mean that any given estimate is close to the true treatment effect. Similarly, while pre-determined covariates will be balanced between treatment and control groups on average, large imbalances may be observed in a given experiment and the researcher may therefore want to condition on such covariates using linear regression. This paper studies the theoretical properties of both the difference-in-means and OLS estimators \emph{conditional} on observed differences in covariates. By deriving the statistical properties of the conditional estimators, we can establish guidance for how to deal with covariate imbalances. We study both inference with OLS, as well as with a new version of Fisher's exact test, where the randomization distribution comes from a small subset of all possible assignment vectors.


title: Sensitivity of metabolic networks: perturbing metabolite concentrations

abstract: Sensitivity studies the network response to perturbations. We consider local perturbations of the concentrations of metabolites at an equilibrium. We investigate the responses in the network, both of the metabolite concentrations and of the reaction fluxes. Our approach is purely qualitative, rather than quantitative. In fact, our analysis is based, solely, on the stoichiometry of the reaction network and we do not require any quantitative information on the reaction rates. Instead, the description is done only in algebraic terms, and the only data required is the network structure. Biological applications include metabolic control and network reconstruction.


title: Top-$k$ Ranking Bayesian Optimization

abstract: This paper presents a novel approach to top-$k$ ranking Bayesian optimization (top-$k$ ranking BO) which is a practical and significant generalization of preferential BO to handle top-$k$ ranking and tie/indifference observations. We first design a surrogate model that is not only capable of catering to the above observations, but is also supported by a classic random utility model. Another equally important contribution is the introduction of the first information-theoretic acquisition function in BO with preferential observation called multinomial predictive entropy search (MPES) which is flexible in handling these observations and optimized for all inputs of a query jointly. MPES possesses superior performance compared with existing acquisition functions that select the inputs of a query one at a time greedily. We empirically evaluate the performance of MPES using several synthetic benchmark functions, CIFAR-$10$ dataset, and SUSHI preference dataset.


title: Fiducial inference then and now

abstract: We conduct a review of the fiducial approach to statistical inference, following its journey from its initiation by R. A. Fisher, through various problems and criticisms, on to its general neglect, and then to its more recent resurgence. Emphasis is laid on the functional model formulation, which helps clarify the very limited conditions under which fiducial inference can be conducted in an unambiguous and self-consistent way.


title: GrADyS: Exploring movement awareness for efficient routing in Ground-and-Air Dynamic Sensor Networks

abstract: Several situations exist where a geographic region of some size needs to be scanned or monitored through many sensors. Still, it is either absolutely impossible or prohibitively expensive to deploy and maintain wireless communication infrastructure for the distributed sensors. Either because the region is hidden behind walls, not easily accessible, hard to get through, or infected with some lethal bacteria or virus transmitter. In this case, the best is to scatter (disposable) sensors in the region and let them transmit the collected sensor data by wireless means to an overflying UAV/drone. Which then physically hauls the collected data from the monitored area to a central base station that functions as a gateway to the Internet. The project GrADyS aims to research two sets of problems regarding such data collection. The former aims to coordinate several autonomous UAVs in a distributed manner to collect the generated data while relying only on ad-hoc communication. The latter aims to develop routing protocols to mesh networks Bluetooth Mesh's Low Power Nodes. Both research lines already present preliminary results that are presented in this paper.


title: Rectifiable curves in proximally smooth sets

abstract: We provide an algorithm of constructing a rectifiable curve between two sufficiently close points of a proximally smooth set in a uniformly convex and uniformly smooth Banach space. Our algorithm returns a reasonably short curve between two sufficiently close points of a proximally smooth set, is iterative and uses a certain modification of the metric projection. We estimate the length of a constructed curve and its deviation from the segment with the same endpoints. These estimates coincide up to a constant factor with those for the geodesics in a proximally smooth set in a Hilbert space.


title: Confused Modulo Projection based Somewhat Homomorphic Encryption -- Cryptosystem, Library and Applications on Secure Smart Cities

abstract: With the development of cloud computing, the storage and processing of massive visual media data has gradually transferred to the cloud server. For example, if the intelligent video monitoring system cannot process a large amount of data locally, the data will be uploaded to the cloud. Therefore, how to process data in the cloud without exposing the original data has become an important research topic. We propose a single-server version of somewhat homomorphic encryption cryptosystem based on confused modulo projection theorem named CMP-SWHE, which allows the server to complete blind data processing without \emph{seeing} the effective information of user data. On the client side, the original data is encrypted by amplification, randomization, and setting confusing redundancy. Operating on the encrypted data on the server side is equivalent to operating on the original data. As an extension, we designed and implemented a blind computing scheme of accelerated version based on batch processing technology to improve efficiency. To make this algorithm easy to use, we also designed and implemented an efficient general blind computing library based on CMP-SWHE. We have applied this library to foreground extraction, optical flow tracking and object detection with satisfactory results, which are helpful for building smart cities. We also discuss how to extend the algorithm to deep learning applications. Compared with other homomorphic encryption cryptosystems and libraries, the results show that our method has obvious advantages in computing efficiency. Although our algorithm has some tiny errors ($10^{-6}$) when the data is too large, it is very efficient and practical, especially suitable for blind image and video processing.


title: Determination of the helpfulness of physics exam study methods

abstract: Studying for physics exams can be difficult and stressful, especially during a student's introductory year in physics. For students who do not plan to major in physics, the desire to do well is based less on understanding concepts and more on achieving a better grade. For this reason, students want to study as efficiently as possible by using the most optimal study methods. We have taken surveys over the past three years to determine how students study for exams and compared that to their exam grades. We found that students who studied using methods that they rated as more helpful did better on the exams. By utilizing the study results, we are able to present our current and future students with study methods that have been rated as being more helpful, and give them advice on ways to optimize their study time for exams.


title: Bayesian unsupervised learning reveals hidden structure in concentrated electrolytes

abstract: Electrolytes play an important role in a plethora of applications ranging from energy storage to biomaterials. Notwithstanding this, the structure of concentrated electrolytes remains enigmatic. Many theoretical approaches attempt to model the concentrated electrolytes by introducing the idea of ion pairs, with ions either being tightly `paired' with a counter-ion, or `free' to screen charge. In this study we reframe the problem into the language of computational statistics, and test the null hypothesis that all ions share the same local environment. Applying the framework to molecular dynamics simulations, we show that this null hypothesis is not supported by data. Our statistical technique suggests the presence of distinct local ionic environments; surprisingly, these differences arise in like charge correlations rather than unlike charge attraction. The resulting fraction of particles in non-aggregated environments shows a universal scaling behaviour across different background dielectric constants and ionic concentrations.


title: An Information-Theoretic Framework for Unifying Active Learning Problems

abstract: This paper presents an information-theoretic framework for unifying active learning problems: level set estimation (LSE), Bayesian optimization (BO), and their generalized variant. We first introduce a novel active learning criterion that subsumes an existing LSE algorithm and achieves state-of-the-art performance in LSE problems with a continuous input domain. Then, by exploiting the relationship between LSE and BO, we design a competitive information-theoretic acquisition function for BO that has interesting connections to upper confidence bound and max-value entropy search (MES). The latter connection reveals a drawback of MES which has important implications on not only MES but also on other MES-based acquisition functions. Finally, our unifying information-theoretic framework can be applied to solve a generalized problem of LSE and BO involving multiple level sets in a data-efficient manner. We empirically evaluate the performance of our proposed algorithms using synthetic benchmark functions, a real-world dataset, and in hyperparameter tuning of machine learning models.


title: A deep learning method for solving Fokker-Planck equations

abstract: The time evolution of the probability distribution of a stochastic differential equation follows the Fokker-Planck equation, which usually has an unbounded, high-dimensional domain. Inspired by our early study in \cite{li2018data}, we propose a mesh-free Fokker-Planck solver, in which the solution to the Fokker-Planck equation is now represented by a neural network. The presence of the differential operator in the loss function improves the accuracy of the neural network representation and reduces the the demand of data in the training process. Several high dimensional numerical examples are demonstrated.


title: On some $q$-series identities related to a generalized divisor function and their implications

abstract: In this article, a $q$-series examined by Kluyver and Uchimura is generalized. This allows us to find generalization of the identities in the random acyclic digraph studied by Simon, Crippa, and Collenberg in 1993. As one of the corollaries of our main theorem, we get results of Dilcher and Andrews, Crippa, and Simon. This main theorem involves a surprising new generalization of the divisor function $\sigma_s(n)$, which we denote by $\sigma_{s,z}(n)$. Analytic properties of $\sigma_{s,z}(n)$ are also studied. As a special case of one of our theorem we obtain result from a recent paper of Bringmann and Jennings-Shaffer.


title: Three-phase coexistence in binary charged lipid membranes in hypotonic solution

abstract: We investigated the phase separation of dioleoylphosphatidylserine (DOPS) and dipalmitoylphosphatidylcholine (DPPC) in giant unilamellar vesicles in hypotonic solution using fluorescence and confocal laser scanning microscopy. Although phase separation in charged lipid membranes is generally suppressed by the electrostatic repulsion between the charged headgroups, osmotic stress can promote the formation of charged lipid domains. Interestingly, we observed three-phase coexistence even in DOPS/DPPC binary lipid mixtures. The three phases were DPPC-rich, dissociated DOPS-rich, and nondissociated DOPS-rich phases. The two forms of DOPS were found to coexist owing to the ionization of the DOPS headgroup, such that the system could be regarded as quasi-ternary. The three formed phases with differently ionized DOPS domains were successfully identified experimentally by monitoring the adsorption of positively charged particles. In addition, coarse-grained molecular dynamics simulations confirmed the stability of the three-phase coexistence. Attraction mediated by hydrogen bonding between protonated DOPS molecules and reduction of the electrostatic interactions at the domain boundaries stabilized the three-phase coexistence.


title: The effects of beta-cell mass and function, intercellular coupling, and islet synchrony on $\textrm{Ca}^{2+}$ dynamics

abstract: Type 2 diabetes (T2D) is a challenging metabolic disorder characterized by a substantial loss of $\beta$-cell mass and alteration of $\beta$-cell function in the islets of Langerhans, disrupting insulin secretion and glucose homeostasis. The mechanisms for deficiency in $\beta$-cell mass and function during the hyperglycemia development and T2D pathogenesis are complex. To study the relative contribution of $\beta$-cell mass to $\beta$-cell function in T2D, we make use of a comprehensive electrophysiological model of human $\beta$-cell clusters. We find that defect in $\beta$-cell mass causes a functional decline in single $\beta$-cell, impairment in intra-islet synchrony, and changes in the form of oscillatory patterns of membrane potential and intracellular $\textrm{Ca}^{2+}$ concentration, which can lead to changes in insulin secretion dynamics and in insulin levels. The model demonstrates a good correspondence between suppression of synchronizing electrical activity and published experimental measurements. We then compare the role of gap junction-mediated electrical coupling with both $\beta$-cell synchronization and metabolic coupling in the behavior of $\textrm{Ca}^{2+}$ concentration dynamics within human islets. Our results indicate that inter-$\beta$-cellular electrical coupling depicts a more important factor in shaping the physiological regulation of islet function and in human T2D. We further predict that varying the whole-cell conductance of delayed rectifier $\textrm{K}^{+}$ channels modifies oscillatory activity patterns of $\beta$-cell population lacking intercellular coupling, which significantly affect $\textrm{Ca}^{2+}$ concentration and insulin secretion.


title: Minimax Strikes Back

abstract: Deep Reinforcement Learning (DRL) reaches a superhuman level of play in many complete information games. The state of the art search algorithm used in combination with DRL is Monte Carlo Tree Search (MCTS). We take another approach to DRL using a Minimax algorithm instead of MCTS and learning only the evaluation of states, not the policy. We show that for multiple games it is competitive with the state of the art DRL for the learning performances and for the confrontations.


title: Entropic-Wasserstein barycenters: PDE characterization, regularity and CLT

abstract: In this paper, we investigate properties of entropy-penalized Wasserstein barycenters introduced by Bigot, Cazelles and Papadakis (2019) as a regularization of Wasserstein barycenters first presented by Agueh and Carlier (2011). After characterizing these barycenters in terms of a system of Monge-Amp\`ere equations, we prove some global moment and Sobolev bounds as well as higher regularity properties. We finally establish a central limit theorem for entropic-Wasserstein barycenters.


title: Searches for Compact Binary Coalescence Events using Neural Networks in LIGO/Virgo Second Observation Period

abstract: We present results on the search for the coalescence of compact binary mergers using convolutional neural networks and the LIGO/Virgo data, corresponding to the O2 observation period. Two-dimensional images in time and frequency are used as input, and two sets of neural networks are trained separately for low mass (0.2 - 2.0 Msun) and high mass (25 - 100 Msun) compact binary coalescence events. We explored neural networks trained with input information from a single or a pair of interferometers, indicating that the use of information from pairs leads to an improved performance. A scan over the full O2 data set using the convolutional neural networks for detection demonstrates that the performance is compatible with that from canonical pipelines using matched filtering techniques. No additional events with significant signal-to-noise ratio are found in the O2 data.


title: A Generally Covariant Theory of Quantized Dirac Field in de Sitter Spacetime

abstract: As a sequel to our previous work\cite{Feng2020}, we propose in this paper a quantization scheme for Dirac field in de Sitter spacetime. Our scheme is covariant under both general transformations and Lorentz transformations. We first present a Hamiltonian structure, then quantize the field following the standard approach of constrained systems. For the free field, the time-dependent quantized Hamiltonian is diagonalized by Bogliubov transformation and the eigen-states at each instant are interpreted as the observed particle states at that instant. The measurable energy-momentum of observed particle/antiparticles are the same as obtained for Klein-Gordon field. Moreover, the energy-momentum also satisfies geodesic equation, a feature justifying its measurability. As in \cite{Feng2020}, though the mathematics is carried out in terms of conformal coordinates for the sake of convenience, the whole theory can be transformed into any other coordinates based on general covariance. It is concluded that particle/antiparticle states, such as vacuum states in particular are time-dependent and vacuum states at one time evolves into non-vacuum states at later times. Formalism of perturbational calculation is provided with an extended Dirac picture.


title: Self-supervised monocular depth estimation from oblique UAV videos

abstract: UAVs have become an essential photogrammetric measurement as they are affordable, easily accessible and versatile. Aerial images captured from UAVs have applications in small and large scale texture mapping, 3D modelling, object detection tasks, DTM and DSM generation etc. Photogrammetric techniques are routinely used for 3D reconstruction from UAV images where multiple images of the same scene are acquired. Developments in computer vision and deep learning techniques have made Single Image Depth Estimation (SIDE) a field of intense research. Using SIDE techniques on UAV images can overcome the need for multiple images for 3D reconstruction. This paper aims to estimate depth from a single UAV aerial image using deep learning. We follow a self-supervised learning approach, Self-Supervised Monocular Depth Estimation (SMDE), which does not need ground truth depth or any extra information other than images for learning to estimate depth. Monocular video frames are used for training the deep learning model which learns depth and pose information jointly through two different networks, one each for depth and pose. The predicted depth and pose are used to reconstruct one image from the viewpoint of another image utilising the temporal information from videos. We propose a novel architecture with two 2D CNN encoders and a 3D CNN decoder for extracting information from consecutive temporal frames. A contrastive loss term is introduced for improving the quality of image generation. Our experiments are carried out on the public UAVid video dataset. The experimental results demonstrate that our model outperforms the state-of-the-art methods in estimating the depths.


title: Photometric and spectroscopic analysis of Comet 29P/Schwassmann-Wachmann 1 activity

abstract: We carried out photometric and spectroscopic observations of comet 29P/Schwassmann-Wachmann 1 at the 6-m BTA telescope (SAO RAS, Russia) and the 1.6-m telescope of the National Laboratory for Astro-physics (LNA, Brazil) on February 20, 2012, and on May 31, 2011, respectively. The spectra revealed the presence of CO+ and N2+ emissions in the cometary coma at a distance of 5.25 AU from the Sun. The ratio [N2+]/[CO+] within the projected slit is 0.013. The images obtained through BVR filters showed a bright, dust coma, indicating a high level of activity. We estimated a colour index and a colour excess for the comet. The parameter Af\r{ho}, which is used as an indicator of a cometary activity, was measured to be 2584 cm in the reference optical aperture of \r{ho}=10^4 km. The dust production constituted 33 kg/s and 9.3*10^3 kg/s, it was obtained using different methods. We also investigated the morphology of the comet using image enhancement techniques and found two jets in the coma.


title: Siamese Anchor Proposal Network for High-Speed Aerial Tracking

abstract: In the domain of visual tracking, most deep learning-based trackers highlight the accuracy but casting aside efficiency. Therefore, their real-world deployment on mobile platforms like the unmanned aerial vehicle (UAV) is impeded. In this work, a novel two-stage Siamese network-based method is proposed for aerial tracking, \textit{i.e.}, stage-1 for high-quality anchor proposal generation, stage-2 for refining the anchor proposal. Different from anchor-based methods with numerous pre-defined fixed-sized anchors, our no-prior method can 1) increase the robustness and generalization to different objects with various sizes, especially to small, occluded, and fast-moving objects, under complex scenarios in light of the adaptive anchor generation, 2) make calculation feasible due to the substantial decrease of anchor numbers. In addition, compared to anchor-free methods, our framework has better performance owing to refinement at stage-2. Comprehensive experiments on three benchmarks have proven the superior performance of our approach, with a speed of around 200 frames/s.


title: Optimal Control of Diffusion Processes with Terminal Constraint in Law

abstract: Stochastic optimal control problems with constraints on the probability distribution of the final output are considered. Necessary conditions for optimality in the form of a coupled system of partial differential equations involving a forward Fokker-Planck equation and a backward Hamilton-Jacobi-Bellman equation are proved using convex duality techniques.


title: Static object detection and segmentation in videos based on dual foregrounds difference with noise filtering

abstract: This paper presents static object detection and segmentation method in videos from cluttered scenes. Robust static object detection is still challenging task due to presence of moving objects in many surveillance applications. The level of difficulty is extremely influenced by on how you label the object to be identified as static that do not establish the original background but appeared in the video at different time. In this context, background subtraction technique based on the frame difference concept is applied to the identification of static objects. Firstly, we estimate a frame differencing foreground mask image by computing the difference of each frame with respect to a static reference frame. The Mixture of Gaussian MOG method is applied to detect the moving particles and then outcome foreground mask is subtracted from frame differencing foreground mask. Pre-processing techniques, illumination equalization and de-hazing methods are applied to handle low contrast and to reduce the noise from scattered materials in the air e.g. water droplets and dust particles. Finally, a set of mathematical morphological operation and largest connected-component analysis is applied to segment the object and suppress the noise. The proposed method was built for rock breaker station application and effectively validated with real, synthetic and two public data sets. The results demonstrate the proposed approach can robustly detect, segmented the static objects without any prior information of tracking.


title: Spatial and Spectral Mode-Selection Effects in Topological Lasers with Frequency-Dependent Gain

abstract: We develop a semiclassical theory of laser oscillation into a chiral edge state of a topological photonic system endowed with a frequency-dependent gain. As an archetypal model of this physics, we consider a Harper-Hofstadter lattice embedding population-inverted two-level atoms as gain material. We show that a suitable design of the spatial distribution of gain and of its spectral shape provides flexible mode selection mechanisms that can stabilize single-mode lasing into an edge state. Implications of our results for recent experiments are outlined.


title: Visuo-Locomotive Complexity as a Component of Parametric Systems for Architecture Design

abstract: A people-centred approach for designing large-scale built-up spaces necessitates systematic anticipation of user's embodied visuo-locomotive experience from the viewpoint of human-environment interaction factors pertaining to aspects such as navigation, wayfinding, usability. In this context, we develop a behaviour-based visuo-locomotive complexity model that functions as a key correlate of cognitive performance vis-a-vis internal navigation in built-up spaces. We also demonstrate the model's implementation and application as a parametric tool for the identification and manipulation of the architectural morphology along a navigation path as per the parameters of the proposed visuospatial complexity model. We present examples based on an empirical study in two healthcare buildings, and showcase the manner in which a dynamic and interactive parametric (complexity) model can promote behaviour-based decision-making throughout the design process to maintain desired levels of visuospatial complexity as part of a navigation or wayfinding experience.


title: Quantum reinforcement learning in continuous action space

abstract: Quantum reinforcement learning (QRL) is one promising algorithm proposed for near-term quantum devices. Early QRL proposals are effective at solving problems in discrete action space, but often suffer from the curse of dimensionality in the continuous domain due to discretization. To address this problem, we propose a quantum Deep Deterministic Policy Gradient algorithm that is efficient at solving both classical and quantum sequential decision problems in the continuous domain. As an application, our method can solve the quantum state-generation problem in a single shot: it only requires a one-shot optimization to generate a model that outputs the desired control sequence for arbitrary target state. In comparison, the standard quantum control method requires optimizing for each target state. Moreover, our method can also be used to physically reconstruct an unknown quantum state.


title: Markov-modulated generalized Ornstein-Uhlenbeck processes and an application in risk theory

abstract: We derive the Markov-modulated generalized Ornstein-Uhlenbeck process by embedding a Markov-modulated random recurrence equation in continuous time. The obtained process turns out to be the unique solution of a certain stochastic differential equation driven by a bivariate Markov-additive process. We present this stochastic differential equation as well as its solution explicitely in terms of the driving Markov-additive process. Moreover, we give necessary and sufficient conditions for strict stationarity of the Markov-modulated generalized Ornstein-Uhlenbeck process, and prove that its stationary distribution is given by the distribution of a specific exponential functional of Markov-additive processes. Finally we propose an application of the Markov-modulated generalized Ornstein-Uhlenbeck process as Markov-modulated risk model with stochastic investment. This generalizes Paulsen's risk process to a Markov-switching environment. We derive a formula in this risk model that expresses the ruin probability in terms of the distribution of an exponential functional of a Markov-additive process.


title: Fundamental Limits and Tradeoffs in Invariant Representation Learning

abstract: A wide range of machine learning applications such as privacy-preserving learning, algorithmic fairness, and domain adaptation/generalization among others, involve learning invariant representations of the data that aim to achieve two competing goals: (a) maximize information or accuracy with respect to a target response, and (b) maximize invariance or independence with respect to a set of protected features (e.g., for fairness, privacy, etc). Despite their wide applicability, theoretical understanding of the optimal tradeoffs -- with respect to accuracy, and invariance -- achievable by invariant representations is still severely lacking. In this paper, we provide an information theoretic analysis of such tradeoffs under both classification and regression settings. More precisely, we provide a geometric characterization of the accuracy and invariance achievable by any representation of the data; we term this feasible region the information plane. We provide an inner bound for this feasible region for the classification case, and an exact characterization for the regression case, which allows us to either bound or exactly characterize the Pareto optimal frontier between accuracy and invariance. Although our contributions are mainly theoretical, a key practical application of our results is in certifying the potential sub-optimality of any given representation learning algorithm for either classification or regression tasks. Our results shed new light on the fundamental interplay between accuracy and invariance, and may be useful in guiding the design of future representation learning algorithms.


title: A Light Field Front-end for Robust SLAM in Dynamic Environments

abstract: There is a general expectation that robots should operate in urban environments often consisting of potentially dynamic entities including people, furniture and automobiles. Dynamic objects pose challenges to visual SLAM algorithms by introducing errors into the front-end. This paper presents a Light Field SLAM front-end which is robust to dynamic environments. A Light Field captures a bundle of light rays emerging from a single point in space, allowing us to see through dynamic objects occluding the static background via Synthetic Aperture Imaging(SAI). We detect apriori dynamic objects using semantic segmentation and perform semantic guided SAI on the Light Field acquired from a linear camera array. We simultaneously estimate both the depth map and the refocused image of the static background in a single step eliminating the need for static scene initialization. The GPU implementation of the algorithm facilitates running at close to real time speeds of 4 fps. We demonstrate that our method results in improved robustness and accuracy of pose estimation in dynamic environments by comparing it with state of the art SLAM algorithms.


title: Multi-Label Noise Robust Collaborative Learning for Remote Sensing Image Classification

abstract: The development of accurate methods for multi-label classification (MLC) of remote sensing (RS) images is one of the most important research topics in RS. The MLC methods based on convolutional neural networks (CNNs) have shown strong performance gains in RS. However, they usually require a high number of reliable training images annotated with multiple land-cover class labels. Collecting such data is time-consuming and costly. To address this problem, the publicly available thematic products, which can include noisy labels, can be used to annotate RS images with zero-labeling cost. However, multi-label noise (which can be associated with wrong and missing label annotations) can distort the learning process of the MLC methods. To address this problem, we propose a novel multi-label noise robust collaborative learning (RCML) method to alleviate the negative effects of multi-label noise during the training phase of a CNN model. RCML identifies, ranks and excludes noisy multi-labels in RS images based on three main modules: 1) the discrepancy module; 2) the group lasso module; and 3) the swap module. The discrepancy module ensures that the two networks learn diverse features, while producing the same predictions. The task of the group lasso module is to detect the potentially noisy labels assigned to multi-labeled training images, while the swap module is devoted to exchange the ranking information between two networks. Unlike the existing methods that make assumptions about noise distribution, our proposed RCML does not make any prior assumption about the type of noise in the training set. The experiments conducted on two multi-label RS image archives confirm the robustness of the proposed RCML under extreme multi-label noise rates. Our code is publicly available at: http://www.noisy-labels-in-rs.org


title: Model-Based Actor-Critic with Chance Constraint for Stochastic System

abstract: Safety is essential for reinforcement learning (RL) applied in real-world situations. Chance constraints are suitable to represent the safety requirements in stochastic systems. Previous chance-constrained RL methods usually have a low convergence rate, or only learn a conservative policy. In this paper, we propose a model-based chance constrained actor-critic (CCAC) algorithm which can efficiently learn a safe and non-conservative policy. Different from existing methods that optimize a conservative lower bound, CCAC directly solves the original chance constrained problems, where the objective function and safe probability is simultaneously optimized with adaptive weights. In order to improve the convergence rate, CCAC utilizes the gradient of dynamic model to accelerate policy optimization. The effectiveness of CCAC is demonstrated by a stochastic car-following task. Experiments indicate that compared with previous RL methods, CCAC improves the performance while guaranteeing safety, with a five times faster convergence rate. It also has 100 times higher online computation efficiency than traditional safety techniques such as stochastic model predictive control.


title: Co-evolution between Codon Usage and Protein-Protein Interaction in Bacteria

abstract: We study the correlation between the codon usage bias of genetic sequences and the network features of protein-protein interaction (PPI) in bacterial species. We use PCA techniques in the space of codon bias indices to show that genes with similar patterns of codon usage have a significantly higher probability that their encoded proteins are functionally connected and interacting. Importantly, this signal emerges when multiple aspects of codon bias are taken into account at the same time. The present study extends our previous observations on E.Coli over a wide set of 34 bacteria. These findings could allow for future investigations on the possible effects of codon bias on the topology of the PPI network, with the aim of improving existing bioinformatics methods for predicting protein interactions.


title: Comet 29P/Schwassmann-Wachmann 1 dust environment from photometric observation at the SOAR Telescope

abstract: We report photometric observations of comet 29P/Schwassmann-Wachmann 1 made on August 12, 2016 with the broadband B, V, R and I filters and the SOAR 4.1-meter telescope (Chile). We find the comet active at that time. Enhanced images obtained in all filters reveal three radial features in the 29P/ Schwassmann-Wachmann 1 coma, regardless of the image-processing algorithm. Using a high-resolution spectrum of comet 29P/Schwassmann-Wachmann 1 reported by Ivanova et al. (2018) on the same date, we estimate the relative contribution of the gaseous emission and the continuum to the total response measured with our broadband B and V filters. The gaseous-emission contribution appears to be very small 2.5%. We compute the dust production Af\r{ho} in 29P/ Schwassmann-Wachmann 1 for the four filters and find its growth with the wavelength, from 3,393 cm in the B filter to 8,561 cm in the I filter. We model the color slope of dust in Comet 29P/Schwassmann-Wachmann 1 using agglomerated debris particles. Simultaneous analysis of the color slope in the B-R and R-I pairs suggests a single dominant chemical species of 29P/Schwassmann-Wachmann 1 dust particles consisting of Fe-Mg silicates and obeying a power-law size distribution with index n=2.55. This conclusion is consistent with the previous thermal-emission study of 29P/Schwassmann-Wachmann 1.


title: On non-locality in the Calculus of Variations

abstract: Non-locality is being intensively studied in various PDE-contexts and in variational problems. The numerical approximation also looks challenging, as well as the application of these models to Continuum Mechanics and Image Analysis, among other areas. Even though there is a growing body of deep and fundamental knowledge about non-locality, for variational principles there are still very basic questions that have not been addressed so far. Taking some of these as a motivation, we describe a general perspective on distinct classes of non-local variational principles setting a program for the analysis of this kind of problems. We start such program with the simplest problem possible: that of scalar, uni-dimensional cases, under a particular class of non-locality. Even in this simple initial scenario, one finds quite unexpected facts to the point that our intuition about local, classic problems can no longer guide us for these new problems. There are three main issues worth highlighting, in the particular situation treated: $\bullet$ natural underlying spaces involve different non-local types of derivatives as, for instance, fractional Sobolev spaces; $\bullet$ no convexity of integrands is required for existence of minimizers; $\bullet$ optimality is formulated in terms of quite special integral equations rather than differential equations. We are thus able to provide some specific answers to the initial questions that motivated our investigation. In subsequent papers, we will move on to consider the higher dimensional situation driven by the possibility that no convexity or quasiconvexity might be involved in weak lower semicontinuity in a full vector, higher dimensional situation.


title: The Confirmation of Three Faint Variable Stars and the Observation of Eleven Others in the Vicinity of Kepler-8b by the Lookout Observatory

abstract: A commissioning survey of the Lookout Observatory has observed fourteen faint (V ~ 13 to 17) variables in the region of the exoplanet Kepler-8b. Three of these are variable star candidates discovered by the Asteroid Terrestrial-Impact Last Alert System (ATLAS) and confirmed here. The ATLAS survey identified 315,000 probable variables within its wide-field survey in 2018. The faintness (down to r ~ 18) and small amplitudes (down to 0.02 mag) included in these candidates makes external validation difficult. The confirmation of this handful of variable stars lends credibility to the ATLAS catalog. Lastly, the agreement between various surveys and this new one validates the use of this instrument for variable star and exoplanet research.


title: The Complex-Scaled Half-Space Matching Method

abstract: The Half-Space Matching (HSM) method has recently been developed as a new method for the solution of 2D scattering problems with complex backgrounds, providing an alternative to Perfectly Matched Layers (PML) or other artificial boundary conditions. Based on half-plane representations for the solution, the scattering problem is rewritten as a system of integral equations in which the unknowns are restrictions of the solution to the boundaries of a finite number of overlapping half-planes contained in the domain: this integral equation system is coupled to a standard finite element discretisation localised around the scatterer. While satisfactory numerical results have been obtained for real wavenumbers, wellposedness and equivalence to the original scattering problem have been established only for complex wavenumbers. In the present paper, by combining the HSM framework with a complex-scaling technique, we provide a new formulation for real wavenumbers which is provably well-posed and has the attraction for computation that the complex-scaled solutions of the integral equation system decay exponentially at infinity. The analysis requires the study of double-layer potential integral operators on intersecting infinite lines, and their analytic continuations. The effectiveness of the method is validated by preliminary numerical results.


title: Lasing condition for trapped modes in subwavelength--wired PT--symmetric resonators

abstract: The ability to control the laser modes within a subwavelength resonator is of key relevance in modern optoelectronics. This work deals with the theoretical research on optical properties of a PT--symmetric nano--scaled dimer formed by two dielectric wires, one is with loss and the other with gain, wrapped with graphene sheets. We show the existence of two non--radiating trapped modes which transform into radiating modes by increasing the gain--loss parameter. Moreover, these modes reach the lasing condition for suitable values of this parameter, a fact that makes these modes to achieve an ultra high quality factor that is manifested on the response of the structure when it is excited by a plane wave. Unlike other mechanism that transform trapped modes into radiating modes, we show that the variation of gain--loss parameter in the balanced loss--gain structure here studied leads to a variation in the phase difference between induced dipole moments on each wires, without appreciable variation in the modulus of these dipole moments. We provide an approximated method that reproduces the main results provided by the rigorous calculation. Our theoretical findings reveal the possibility to develop unconventional optical devices and structures with enhanced functionality.


title: SymFields: An Open Source Symbolic Fields Analysis Tool for General Curvilinear Coordinates in Python

abstract: An open source symbolic tool for vector fields analysis 'SymFields' is developed in Python. The SymFields module is constructed upon Python symbolic module sympy, which could only conduct scaler field analysis. With SymFields module, you can conduct vector analysis for general curvilinear coordinates regardless whether it is orthogonal or not. In SymFields, the differential operators based on metric tensor are normalized to real physical values, which means your can use real physical value of the vector fields as inputs. This could greatly free the physicists from the tedious calculation under complicated coordinates.


title: Extremely low linear polarization of comet C/2018 V1 (Machholz-Fujikawa-Iwamoto)

abstract: We measured the degree of linear polarization P of comet C/2018 V1 (Machholz-Fujikawa-Iwamoto) with the broadband Johnson V filter in mid-November of 2018. Within a radius of \r{ho}=17,000 km of the inner coma, we detected an extremely low linear polarization at phase angles from 83 to 91.2 degree and constrained the polarization maximum to Pmax = (6.8 +/- 1.8)%. This is the lowest Pmax ever measured in a comet. Using model agglomerated debris particles, we reproduced the polarimetric response of comet C/2018 V1. Four retrieved refractive indices closely match what was experimentally found in Mg-rich silicates with little or no iron content. Moreover, the size distribution of the agglomerated debris particles appears in good quantitative agreement with the in situ findings of comet 1P/Halley. The dust model of polarization of comet C/2018 V1 suggests a strongly negative polarization with amplitude |Pmin| = 5%-7%; whereas, an interpretation based on gaseous emission requires no negative polarization at small phase angles. This dramatic difference could be used to discriminate gaseous-emission and dust explanations in low-Pmax comets in future.


title: Analysis of NARXNN for State of Charge Estimation for Li-ion Batteries on various Drive Cycles

abstract: Electric Vehicles (EVs) are rapidly increasing in popularity as they are environment friendly. Lithium Ion batteries are at the heart of EV technology and contribute to most of the weight and cost of an EV. State of Charge (SOC) is a very important metric which helps to predict the range of an EV. There is a need to accurately estimate available battery capacity in a battery pack such that the available range in a vehicle can be determined. There are various techniques available to estimate SOC. In this paper, a data driven approach is selected and a Nonlinear Autoregressive Network with Exogenous Inputs Neural Network (NARXNN) is explored to accurately estimate SOC. NARXNN has been shown to be superior to conventional Machine Learning techniques available in the literature. The NARXNN model is developed and tested on various EV Drive Cycles like LA92, US06, UDDS and HWFET to test its performance on real world scenarios. The model is shown to outperform conventional statistical machine learning methods and achieve a Mean Squared Error (MSE) in the 1e-5 range.


title: Stability and oscillation of linear delay differential equations

abstract: There is a close connection between stability and oscillation of delay differential equations. For the first-order equation $$ x^{\prime}(t)+c(t)x(\tau(t))=0,~~t\geq 0, $$ where $c$ is locally integrable of any sign, $\tau(t)\leq t$ is Lebesgue measurable, $\lim_{t\rightarrow\infty}\tau(t)=\infty$, we obtain sharp results, relating the speed of oscillation and stability. We thus unify the classical results of Myshkis and Lillo. We also generalise the $3/2-$stability criterion to the case of measurable parameters, improving $1+1/e$ to the sharp $3/2$ constant.


title: Interaction rates in cosmology: heavy particle production and scattering

abstract: We study transition rates and cross sections from first principles in a spatially flat radiation dominated cosmology. We consider a model of scalar particles to study scattering and heavy particle production from pair annihilation, drawing more general conclusions. The S-matrix formulation is ill suited to study these ubiquitous processes in a rapidly expanding cosmology. We introduce a physically motivated adiabatic expansion that relies on wavelengths much smaller than the particle horizon at a given time. The leading order in this expansion dominates the transition rates and cross sections. Several important and general results are direct consequences of the cosmological redshift and a finite particle horizon: i) a violation of local Lorentz invariance, ii) freeze-out of the production cross section at a finite time, iii) sub-threshold production of heavier particles as a consequence of the uncertainty in the local energy from a finite particle horizon, a manifestation of the \emph{antizeno} effect. If heavy dark matter is produced via annihilation of a lighter species, sub-threshold production yields an enhanced abundance. We discuss several possible consequences of these effects.


title: Political Posters Identification with Appearance-Text Fusion

abstract: In this paper, we propose a method that efficiently utilizes appearance features and text vectors to accurately classify political posters from other similar political images. The majority of this work focuses on political posters that are designed to serve as a promotion of a certain political event, and the automated identification of which can lead to the generation of detailed statistics and meets the judgment needs in a variety of areas. Starting with a comprehensive keyword list for politicians and political events, we curate for the first time an effective and practical political poster dataset containing 13K human-labeled political images, including 3K political posters that explicitly support a movement or a campaign. Second, we make a thorough case study for this dataset and analyze common patterns and outliers of political posters. Finally, we propose a model that combines the power of both appearance and text information to classify political posters with significantly high accuracy.


title: Universality in microdroplet nucleation during solvent exchange in Hele-Shaw like channels

abstract: Micro and nanodroplets have many important applications such as in drug delivery, liquid-liquid extraction, nanomaterial synthesis and cosmetics. A commonly used method to generate a large number of micro or nanodroplets in one simple step is solvent exchange (also called nanoprecipitation), in which a good solvent of the droplet phase is displaced by a poor one, generating an oversaturation pulse that leads to droplet nucleation. Despite its crucial importance, the droplet growth resulting from the oversaturation pulse in this ternary system is still poorly understood. We experimentally and theoretically study this growth in Hele-Shaw like channels by measuring the total volume of the oil droplets that nucleates out of it. In order to prevent the oversaturated oil from exiting the channel, we decorated some of the channels with a porous region in the middle. Solvent exchange is performed with various solution compositions, flow rates and channel geometries, and the measured droplets volume is found to increase with the P\'eclet number $Pe$ with an approximate effective power law $V\propto Pe^{0.50}$. A theoretical model is developed to account for this finding. With this model we can indeed explain the $V\propto Pe^{1/2}$ scaling, including the prefactor, which can collapse all data of the "porous" channels onto one universal curve, irrespective of channel geometry and composition of the mixtures. Our work provides a macroscopic approach to this bottom-up method of droplet generation and may guide further studies on oversaturation and nucleation in ternary systems.


title: Light Higgs searches in $t \bar t \phi$ production at the LHC

abstract: In this paper we propose a new reconstruction method to explore the low mass region in the associated production of top-quark pairs ($t\bar{t}$) with a generic scalar boson ($\phi$) at the LHC. The new method of mass reconstruction shows an improved resolution of at least a factor of two in the low mass region when compared to previous methods, without the loss of sensitivity of previous analyses. It turns out that it also leads to an improvement of the mass reconstruction of the 125 GeV Higgs for the same production process. We use an effective Lagrangian to describe a scalar with a generic Yukawa coupling to the top quarks. A full phenomenological analysis was performed, using Standard Model background and signal events generated with MadGraph5\_aMC@NLO and reconstructed using a kinematic fit. The use of CP-sensitive variables allows then to maximize the distinction between CP-even and CP-odd components of the Yukawa couplings. Confidence Levels (CLs) for the exclusion of $\phi$ bosons with mixed CP (both CP-even and CP-odd components) were determined as a function of the top Yukawa couplings to the $\phi$ boson. The mass range analysed starts slightly above the $\Upsilon$ mass up to 40 GeV, although the analysis can be used for an arbitrary mass. %We focus on dileptonic final states of the $t\bar{t}\phi$ system, with $\phi\rightarrow b\bar{b}$. If no new light scalar is found, exclusion limits at 95\% CL for the absolute value of the CP-even and CP-odd Yukawa are derived. %couplings are expected to be, approximately, as low as 0.10 and 0.50, respectively, at the end of the High Luminosity phase of the LHC (HL-LHC). Finally, we analyse how these limits constrain the parameter space of the complex two-Higgs doublet model (C2HDM).


title: Stability from graph symmetrisation arguments with applications to inducibility

abstract: We present a sufficient condition for the stability property of extremal graph problems that can be solved via Zykov's symmetrisation. Our criterion is stated in terms of an analytic limit version of the problem. We show that, for example, it applies to the inducibility problem for an arbitrary complete bipartite graph $B$, which asks for the maximum number of induced copies of $B$ in an $n$-vertex graph, and to the inducibility problem for $K_{2,1,1,1}$ and $K_{3,1,1}$, the only complete partite graphs on at most five vertices for which the problem was previously open.


title: DCCRGAN: Deep Complex Convolution Recurrent Generator Adversarial Network for Speech Enhancement

abstract: Generative adversarial network (GAN) still exists some problems in dealing with speech enhancement (SE) task. Some GAN-based systems adopt the same structure from Pixel-to-Pixel directly without special optimization. The importance of the generator network has not been fully explored. Other related researches change the generator network but operate in the time-frequency domain, which ignores the phase mismatch problem. In order to solve these problems, a deep complex convolution recurrent GAN (DCCRGAN) structure is proposed in this paper. The complex module builds the correlation between magnitude and phase of the waveform and has been proved to be effective. The proposed structure is trained in an end-to-end way. Different LSTM layers are used in the generator network to sufficiently explore the speech enhancement performance of DCCRGAN. The experimental results confirm that the proposed DCCRGAN outperforms the state-of-the-art GAN-based SE systems.


title: Attraction and condensation of driven tracers in a narrow channel

abstract: Emergent bath-mediated attraction and condensation arise when multiple particles are simultaneously driven through an equilibrated bath under geometric constraints. While such scenarios are observed in a variety of non-equilibrium phenomena, with an abundance of experimental and numerical evidence, little quantitative understanding of how these interactions arise is currently available. Here we approach the problem by studying the behavior of two driven "tracer" particles, propagating through a bath in a 1D lattice with excluded-volume interactions. We analytically explore the mechanism responsible for the tracers' emergent interactions and compute the resulting effective attractive potential. This mechanism is then numerically shown to extend to a realistic model of hard driven Brownian discs confined to a narrow 2D channel.


title: Analysis of a two-fluid Taylor-Couette flow with one non-Newtonian fluid

abstract: We study the dynamic behaviour of two viscous fluid films confined between two concentric cylinders rotating at a small relative velocity. It is assumed that the fluids are immiscible and that the volume of the outer fluid film is large compared to the volume of the inner one. Moreover, while the outer fluid is considered to have constant viscosity, the rheological behaviour of the inner thin film is determined by a strain-dependent power-law. Starting from a Navier--Stokes system, we formally derive evolution equations for the interface separating the two fluids. Two competing effects drive the dynamics of the interface, namely, the surface tension and the shear stresses induced by the rotation of the cylinders. When the two effects are comparable, the solutions behave, for large times, as in the Newtonian regime. We also study the regime in which the surface tension effects dominate the stresses induced by the rotation of the cylinders. In this case, we prove local existence of positive weak solutions both for shear-thinning and shear-thickening fluids. In the latter case, we show that interfaces which are initially close to a circle converge to a circle in finite time and keep that shape for later times.


title: The role of time estimation in decreased impatience in Intertemporal Choice

abstract: The role of specific cognitive processes in deviations from constant discounting in intertemporal choice is not well understood. We evaluated decreased impatience in intertemporal choice tasks independent of discounting rate and non-linearity in long-scale time representation; nonlinear time representation was expected to explain inconsistencies in discounting rate. Participants performed temporal magnitude estimation and intertemporal choice tasks. Psychophysical functions for time intervals were estimated by fitting linear and power functions, while discounting functions were estimated by fitting exponential and hyperbolic functions. The temporal magnitude estimates of 65% of the participants were better fit with power functions (mostly compression). 63% of the participants had intertemporal choice patterns corresponding best to hyperbolic functions. Even when the perceptual bias in the temporal magnitude estimations was compensated in the discounting rate computation, the data of 8 out of 14 participants continued exhibiting temporal inconsistency. The results suggest that temporal inconsistency in discounting rate can be explained to different degrees by the bias in temporal representations. Non-linearity in temporal representation and discounting rate should be evaluated on an individual basis. Keywords: Intertemporal choice, temporal magnitude, model comparison, impatience, time inconsistency


title: Reconfigurable Intelligent Surface assisted Multi-user Communications: How Many Reflective Elements Do We Need?

abstract: Reconfigurable intelligent surfaces (RISs) consisting of multiple reflective elements are a promising technique to enhance communication quality as they can create favorable propagation conditions. In this letter, we characterize the fundamental relations between the number of reflective elements and the system sum-rate in RIS-assisted multi-user communications. It is known from previous works that the received signal-to-noise ratio~(SNR) can linearly increase with the squared number of RIS reflective elements, but how many elements are sufficient to provide an acceptable system sum-rate still remains an open problem. To this end, we derive the asymptotic capacity with zero-forcing (ZF) precoding, and then discuss how many reflective elements are required so that the ratio of the system sum-rate to the capacity can exceed a predefined threshold. Numerical results verify our analysis.


title: (Decision and regression) tree ensemble based kernels for regression and classification

abstract: Tree based ensembles such as Breiman's random forest (RF) and Gradient Boosted Trees (GBT) can be interpreted as implicit kernel generators, where the ensuing proximity matrix represents the data-driven tree ensemble kernel. Kernel perspective on the RF has been used to develop a principled framework for theoretical investigation of its statistical properties. Recently, it has been shown that the kernel interpretation is germane to other tree-based ensembles e.g. GBTs. However, practical utility of the links between kernels and the tree ensembles has not been widely explored and systematically evaluated. Focus of our work is investigation of the interplay between kernel methods and the tree based ensembles including the RF and GBT. We elucidate the performance and properties of the RF and GBT based kernels in a comprehensive simulation study comprising of continuous and binary targets. We show that for continuous targets, the RF/GBT kernels are competitive to their respective ensembles in higher dimensional scenarios, particularly in cases with larger number of noisy features. For the binary target, the RF/GBT kernels and their respective ensembles exhibit comparable performance. We provide the results from real life data sets for regression and classification to show how these insights may be leveraged in practice. Overall, our results support the tree ensemble based kernels as a valuable addition to the practitioner's toolbox. Finally, we discuss extensions of the tree ensemble based kernels for survival targets, interpretable prototype and landmarking classification and regression. We outline future line of research for kernels furnished by Bayesian counterparts of the frequentist tree ensembles.


title: Phase retrieval by projections in R^n requires 2n-2 projections (This paper has a gap)

abstract: We will answer the most significant open problem in real phase retrieval by projections by showing it requires at least $2n-2$ projections to do phase retrieval in $\RR^n$.


title: Digital Reconstruction of Elmina Castle for Mobile Virtual Reality via Point-based Detail Transfer

abstract: Reconstructing 3D models from large, dense point clouds is critical to enable Virtual Reality (VR) as a platform for entertainment, education, and heritage preservation. Existing 3D reconstruction systems inevitably make trade-offs between three conflicting goals: the efficiency of reconstruction (e.g., time and memory requirements), the visual quality of the constructed scene, and the rendering speed on the VR device. This paper proposes a reconstruction system that simultaneously meets all three goals. The key idea is to avoid the resource-demanding process of reconstructing a high-polygon mesh altogether. Instead, we propose to directly transfer details from the original point cloud to a low polygon mesh, which significantly reduces the reconstruction time and cost, preserves the scene details, and enables real-time rendering on mobile VR devices. While our technique is general, we demonstrate it in reconstructing cultural heritage sites. We for the first time digitally reconstruct the Elmina Castle, a UNESCO world heritage site at Ghana, from billions of laser-scanned points. The reconstruction process executes on low-end desktop systems without requiring high processing power, making it accessible to the broad community. The reconstructed scenes render on Oculus Go in 60 FPS, providing a real-time VR experience with high visual quality. Our project is part of the Digital Elmina effort (http://digitalelmina.org/) between University of Rochester and University of Ghana.


title: An energy stable and maximum bound preserving scheme with variable time steps for time fractional Allen-Cahn equation

abstract: In this work, we propose a Crank-Nicolson-type scheme with variable steps for the time fractional Allen-Cahn equation. The proposed scheme is shown to be unconditionally stable (in a variational energy sense), and is maximum bound preserving. Interestingly, the discrete energy stability result obtained in this paper can recover the classical energy dissipation law when the fractional order $\alpha \rightarrow 1.$ That is, our scheme can asymptotically preserve the energy dissipation law in the $\alpha \rightarrow 1$ limit. This seems to be the first work on variable time-stepping scheme that can preserve both the energy stability and the maximum bound principle. Our Crank-Nicolson scheme is build upon a reformulated problem associated with the Riemann-Liouville derivative. As a by product, we build up a reversible transformation between the L1-type formula of the Riemann-Liouville derivative and a new L1-type formula of the Caputo derivative, with the help of a class of discrete orthogonal convolution kernels. This is the first time such a \textit{discrete} transformation is established between two discrete fractional derivatives. We finally present several numerical examples with an adaptive time-stepping strategy to show the effectiveness of the proposed scheme.


title: Tree approximation in quasi-trees

abstract: In this paper we investigate the geometric properties of quasi-trees, and prove some equivalent criteria. We give a general construction of a tree that approximates the ends of a geodesic space, and use this to prove that every quasi-tree is $(1,C)$-quasi-isometric to a simplicial tree. As a consequence, we show that Gromov's tree approximation lemma for hyperbolic spaces can be improved in the case of quasi-trees to give a uniform approximation for any set of points, independent of cardinality. From this we show that having uniform tree approximation for finite subsets is equivalent to being able to uniformly approximate the entire space by a tree. As another consequence, we note that the boundary of a quasi-tree is isometric to the boundary of its approximating tree under a certain choice of visual metric, and that this gives a natural extension of the standard metric on the boundary of a tree.


title: Arithmetic statistics of Galois groups

abstract: We develop a computational framework for the statistical characterization of Galois characters with finite image, with application to characterizing Galois groups and establishing equivalence of characters of finite images of $\mathrm{Gal}(\overline{\mathbb{Q}}/\mathbb{Q})$.


title: Perturbed Adaptive Belief Propagation Decoding for High-Density Parity-Check Codes

abstract: Algebraic codes such as BCH code are receiving renewed interest as their short block lengths and low/no error floors make them attractive for ultra-reliable low-latency communications (URLLC) in 5G wireless networks. This paper aims at enhancing the traditional adaptive belief propagation (ABP) decoding, which is a soft-in-soft-out (SISO) decoding for high-density parity-check (HDPC) algebraic codes, such as Reed-Solomon (RS) codes, Bose-Chaudhuri-Hocquenghem (BCH) codes, and product codes. The key idea of traditional ABP is to sparsify certain columns of the parity-check matrix corresponding to the least reliable bits with small log-likelihood-ratio (LLR) values. This sparsification strategy may not be optimal when some bits have large LLR magnitudes but wrong signs. Motivated by this observation, we propose a Perturbed ABP (P-ABP) to incorporate a small number of unstable bits with large LLRs into the sparsification operation of the parity-check matrix. In addition, we propose to apply partial layered scheduling or hybrid dynamic scheduling to further enhance the performance of P-ABP. Simulation results show that our proposed decoding algorithms lead to improved error correction performances and faster convergence rates than the prior-art ABP variants.


title: GlocalNet: Class-aware Long-term Human Motion Synthesis

abstract: Synthesis of long-term human motion skeleton sequences is essential to aid human-centric video generation with potential applications in Augmented Reality, 3D character animations, pedestrian trajectory prediction, etc. Long-term human motion synthesis is a challenging task due to multiple factors like, long-term temporal dependencies among poses, cyclic repetition across poses, bi-directional and multi-scale dependencies among poses, variable speed of actions, and a large as well as partially overlapping space of temporal pose variations across multiple class/types of human activities. This paper aims to address these challenges to synthesize a long-term (> 6000 ms) human motion trajectory across a large variety of human activity classes (>50). We propose a two-stage activity generation method to achieve this goal, where the first stage deals with learning the long-term global pose dependencies in activity sequences by learning to synthesize a sparse motion trajectory while the second stage addresses the generation of dense motion trajectories taking the output of the first stage. We demonstrate the superiority of the proposed method over SOTA methods using various quantitative evaluation metrics on publicly available datasets.


title: Prevalence Estimation from Random Samples and Census Data with Participation Bias

abstract: Countries officially record the number of COVID-19 cases based on medical tests of a subset of the population with unknown participation bias. For prevalence estimation, the official information is typically discarded and, instead, small random survey samples are taken. We derive (maximum likelihood and method of moment) prevalence estimators, based on a survey sample, that additionally utilize the official information, and that are substantially more accurate than the simple sample proportion of positive cases. Put differently, using our estimators, the same level of precision can be obtained with substantially smaller survey samples. We take into account the possibility of measurement errors due to the sensitivity and specificity of the medical testing procedure. The proposed estimators and associated confidence intervals are implemented in the companion open source R package cape.


title: A Spatiotemporal Functional Model for Bike-Sharing Systems -- An Example based on the City of Helsinki

abstract: Understanding the usage patterns for bike-sharing systems is essential in terms of supporting and enhancing operational planning for such schemes. Studies have demonstrated how factors such as weather conditions influence the number of bikes that should be available at bike-sharing stations at certain times during the day. However, the influences of these factors usually vary over the course of a day, and if there is good temporal resolution, there could also be significant effects only for some hours/minutes (rush hours, the hours when shops are open, and so forth). Thus, in this paper, an analysis of Helsinki's bike-sharing data from 2017 is conducted that considers full temporal and spatial resolutions. Moreover, the data are available at a very high frequency. Hence, the station hire data is analysed in a spatiotemporal functional setting, where the number of bikes at a station is defined as a continuous function of the time of day. For this completely novel approach, we apply a functional spatiotemporal hierarchical model to investigate the effect of environmental factors and the magnitude of the spatial and temporal dependence. Challenges in computational complexity are faced using a bootstrapping approach. The results show the necessity of splitting the bike-sharing stations into two clusters based on the similarity of their spatiotemporal functional observations in order to model the station hire data of Helsinki's bike-sharing system effectively. The estimated functional influences of the proposed factors are different for the two clusters. Moreover, the estimated parameters reveal high random effects in the data that are not explained by the mean of the process. In this random-effects model, the temporal autoregressive parameter dominates the spatial dependence.


title: Two-sided inequalities for the density function's maximum of weighted sum of chi-square variables

abstract: Two--sided bounds are constructed for a probability density function of a weighted sum of chi-square variables. Both cases of central and non-central chi-square variables are considered. The upper and lower bounds have the same dependence on the parameters of the sum and differ only in absolute constants. The estimates obtained will be useful, in particular, when comparing two Gaussian random elements in a Hilbert space and in multidimensional central limit theorems, including the infinite-dimensional case.


title: Sequential detection of Replay attacks

abstract: One of the most studied forms of attacks on the cyber-physical systems is the replay attack. The statistical similarities of the replay signal and the true observations make the replay attack difficult to detect. In this paper, we have addressed the problem of replay attack detection by adding watermarking to the control inputs and then performed resilient detection using cumulative sum (CUSUM) test on the joint statistics of the innovation signal and the watermarking signal. We derive the expression of the Kullback-Liebler divergence (KLD) between the two joint distributions before and after the replay attack, which is asymptotically inversely proportional to the detection delay. We perform structural analysis of the derived KLD expression and suggest a technique to improve the KLD for the systems with relative degree greater than one. A scheme to find the optimal watermarking signal variance for a fixed increase in the control cost to maximize the KLD under the CUSUM test is presented. We provide various numerical simulation results to support our theory. The proposed method is also compared with a state-of-the-art method.


title: Percolation Transitions in Growing Networks Under Achlioptas Processes: Analytic Solutions

abstract: Networks are ubiquitous in diverse real-world systems. Many empirical networks grow as the number of nodes increases with time. Percolation transitions in growing random networks can be of infinite order. However, when the growth of large clusters is suppressed under some effects, e.g., the Achlioptas process, the transition type changes to the second order. However, analytical results for the critical behavior, such as the transition point, critical exponents, and scaling relations are rare. Here, we derived them explicitly as a function of a control parameter $m$ representing the suppression strength using the scaling ansatz. We then confirmed the results by solving the rate equation and performing numerical simulations. Our results clearly show that the transition point approaches unity and the order-parameter exponent $\beta$ approaches zero algebraically as $m \to \infty$, whereas they approach these values exponentially for a static network. Moreover, the upper critical dimension becomes $d_u=4$ for growing networks, whereas it is $d_u=2$ for static ones.


title: Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

abstract: Understanding and predicting divertor heat-load width ${\lambda}_q$ is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for ${\lambda}_q$ using the extreme-scale edge gyrokinetic code XGC1 in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C.S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013); R.J. Goldston, Nucl. Fusion 52, 013009 (2012)], and furthermore predicted over six times wider ${\lambda}_q$ than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for ${\lambda}_q$ as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide ${\lambda}_q^X{GC}$ for the full-current Q = 10 ITER plasma. The new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider ${\lambda}q_^{XGC} is identified to be the weakly-collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.