title: Optomechanical parametric oscillation of a quantum light-fluid lattice

abstract: Two-photon coherent states are one of the main building pillars of non-linear and quantum optics. It is the basis for the generation of minimum-uncertainty quantum states and entangled photon pairs, applications not obtainable from standard coherent states or one-photon lasers. Here we describe a fully-resonant optomechanical parametric amplifier involving a polariton condensate in a trap lattice quadratically coupled to mechanical modes. The quadratic coupling derives from non-resonant virtual transitions to extended discrete excited states induced by the optomechanical coupling. Non-resonant continuous wave (cw) laser excitation leads to striking experimental consequences, including the emergence of optomechanically induced inter-site parametric oscillations and inter-site tunneling of polaritons at discrete inter-trap detunings corresponding to sums of energies of the two involved mechanical oscillations (20 and 60 GHz confined vibrations). We show that the coherent mechanical oscillations correspond to parametric resonances with threshold condition different to that of standard linear optomechanical self-oscillation. The associated Arnold tongues display a complex scenario of states within the instability region. The observed new phenomena can have applications for the generation of entangled phonon pairs, squeezed mechanical states relevant in sensing and quantum computation, and for the bidirectional frequency conversion of signals in a technologically relevant range.

categories: cond-mat.mes-hall cond-mat.quant-gas physics.optics quant-ph


title: When are Iterative Gaussian Processes Reliably Accurate?

abstract: While recent work on conjugate gradient methods and Lanczos decompositions have achieved scalable Gaussian process inference with highly accurate point predictions, in several implementations these iterative methods appear to struggle with numerical instabilities in learning kernel hyperparameters, and poor test likelihoods. By investigating CG tolerance, preconditioner rank, and Lanczos decomposition rank, we provide a particularly simple prescription to correct these issues: we recommend that one should use a small CG tolerance ($\epsilon \leq 0.01$) and a large root decomposition size ($r \geq 5000$). Moreover, we show that L-BFGS-B is a compelling optimizer for Iterative GPs, achieving convergence with fewer gradient updates.


title: Investigating Galactic binary cluster candidates with Gaia EDR3

abstract: A number of stellar open cluster (OC) pairs in the Milky Way occupy similar positions in the phase space (coordinates, parallax and proper motions) and therefore may constitute physically interacting systems. The characterization of such objects based on observational data is a fundamental step towards a proper understanding of their physical status and to investigate cluster pair formation in the Galaxy. In this work, we employed Gaia EDR3 data to investigate a set of 16 OCs distributed as seven stellar aggregates. We determined structural parameters and applied a decontamination technique that allowed to obtain unambiguous lists of member stars. The studied OCs span Galactocentric distances and ages in the ranges ~7 < $R_G$ (kpc) < ~11 and 7.3 <= log t <= 9.2. Eight OCs were found to constitute 4 gravitationally bound pairs (NGC5617-Trumpler22, Collinder394-NGC6716, Ruprecht100-Ruprecht101, NGC659-NGC663, the latter being a dynamically unevolved binary) and other 4 clusters constitute 2 interacting, but gravitationally unbound, pairs (King16-Berkeley4, NGC2383-NGC2384, the latter being a dissolving OC). Other 4 OCs (Dias1, Pismis19, Czernik20, NGC1857) seem not associated to any stellar aggregates. Apparently, clusters within bound and dynamically evolved pairs tend to present ratios of half-light to tidal radius larger than single clusters located at similar $R_G$, suggesting that mutual tidal interactions may possibly affect their structural parameters. Unbound or dynamically unevolved systems seems to present less noticeable signature of tidal forces on their structure. Moreover, the core radius seems more importantly correlated with the clusters' internal dynamical relaxation process.


title: Homomorphisms of distributive lattices as restrictions of congruences. III. Rectangular lattices and two convex sublattices

abstract: Let $L$ be a finite lattice and let $I$ be an ideal of $L$. Then the restriction map is a bounded lattice homomorphism of the congruence lattice of~$L$ into the congruence lattice of $I$. In a 2009 paper, the authors proved the converse. In a 2012 paper, G. Cz\'edli proved an analogous result for rectangular lattices. In this paper, we prove a stronger form of Cz\'edli's result and provide a short, elementary, and direct proof.


title: Observational Constraints on $f(T)$ Gravity from Model-Independent Data

abstract: We establish new constraints on $f(T)$ gravity models by using cosmological data. In particular, we investigate the restrictions given by the gas mass fraction measurements of galaxy clusters and transversal BAO data. Both data sets are regarded as weakly dependent on a fiducial cosmology. In addition, we also include a CMB measurement of the temperature power spectrum first peak, along with $H(z)$ values from cosmic chronometers and supernovae data from the Pantheon data set. We also perform a forecast for future constraints on the deviation of $f(T)$ models from the $\Lambda$CDM scenario by following the specifications of the J-PAS and Euclid surveys and find significant improvements on the constraints of the $b$-parameter, when compared to the results of the statistical analysis.


title: Benign Overfitting in Adversarially Robust Linear Classification

abstract: "Benign overfitting", where classifiers memorize noisy training data yet still achieve a good generalization performance, has drawn great attention in the machine learning community. To explain this surprising phenomenon, a series of works have provided theoretical justification in over-parameterized linear regression, classification, and kernel methods. However, it is not clear if benign overfitting still occurs in the presence of adversarial examples, i.e., examples with tiny and intentional perturbations to fool the classifiers. In this paper, we show that benign overfitting indeed occurs in adversarial training, a principled approach to defend against adversarial examples. In detail, we prove the risk bounds of the adversarially trained linear classifier on the mixture of sub-Gaussian data under $\ell_p$ adversarial perturbations. Our result suggests that under moderate perturbations, adversarially trained linear classifiers can achieve the near-optimal standard and adversarial risks, despite overfitting the noisy training data. Numerical experiments validate our theoretical findings.


title: The thermodynamic trends of intrinsic defects in primary halide perovskites: A first-principles study

abstract: Defects in halide perovskites play an essential role in determining the efficiency and stability of the resulting optoelectronic devices. Here, we present a systematic study of intrinsic point defects in six primary metal halide perovskites, MAPbI$_3$, MAPbBr$_3$, MAPbCl$_3$, FAPbI$_3$, CsPbI$_3$ and MASnI$_3$, using density functional theory calculations with the SCAN+rVV10 functional. We analyse the impact of changing anions and cations on the defect formation energies and the charge state transitions levels and identify the physical origins underlying the observed trends. Dominant defects in the lead-iodide compounds are the A$^+$ cation interstitials (A = Cs, MA, FA), charge-compensated by I$^-$ interstitials or lead $({2-})$ vacancies. In the lead-bromide and -chloride compounds, halide vacancies become relatively more prominent, and for MAPbBr$_3$, the Pb$^{2+}$ interstitial also becomes important. The trends can be explained in terms of the changes in electrostatic interactions and chemical bonding upon replacing cations and anions. Defect physics in MASnI$_3$ is strongly dominated by tin $({2-})$ vacancies, promoted by the easy oxidation of the tin perovskite. Intrinsically, all compounds are mildly p-doped, except for MASnI$_3$, which is strongly p-doped. All acceptor levels created by defects in the six perovskites are shallow. Some defects, halide vacancies and Pb or Sn interstitials in particular, create deep donor traps. Although these traps might hamper the electronic behavior of MAPbBr$_3$ and MAPbCl$_3$, in iodine-based perovskites their equilibrium concentrations are too small to affect the materials' properties.


title: Modeling COVID-19 Transmission using IDSIM, an Epidemiological-Modelling Desktop App with Multi-Level Immunization Capabilities

abstract: The COVID-19 pandemic has placed unprecedented demands on local public health units in Ontario, Canada, one of which was the need for in-house epidemiological-modelling capabilities. To address this need, Ontario Tech University and the Durham Region Health Department developed a native Windows desktop app that performs epidemiological modelling of infectious diseases. The app is an implementation of a multi-stratified compartmental epidemiological model that can accommodate multiple virus variants and levels of vaccination, as well as public health measures such as physical distancing, contact tracing followed by quarantine, and testing followed by isolation. This article presents the epidemiological model and epidemiological-simulation results obtained using the developed app. The simulations investigate the effects of different factors on COVID-19 transmission in Durham Region, including vaccination coverage, vaccine effectiveness, waning of vaccine-induced immunity, advent of the Omicron variant and effect of COVID-19 booster vaccines in reducing the number of infections and severe cases. Results indicate that, for the Delta variant, natural immunity, in addition to vaccination-induced immunity, is necessary to achieve herd immunity and that waning of vaccine-induced immunity lengthens the time necessary to reach herd immunity. In the absence of additional public health measures, a wave driven by the Omicron variant is predicted to pose significant public health challenges with infections predicted to peak in approximately two to three months, depending on the rate of administration of booster doses.


title: First order linear logic and tensor type calculus for categorial grammars

abstract: We study relationship between first order multiplicative linear logic (MLL1), which has been known to provide representations to different categorial grammars, and the recently introduced extended tensor type calculus (ETTC). We identify a fragment of MLL1, which seems sufficient for many grammar representations, and establish a correspondence between ETTC and this fragment. The system ETTC, thus, can be seen as an alternative syntax and intrinsic deductive system together with a geometric representation for the latter. We also give a natural deduction formulation of ETTC, which might be convenient.


title: Two-Stage and Sequential Unbiased Estimation of N in Binomial Trials, when the Probability of Success p is Unknown

abstract: We propose two-stage and sequential procedures to estimate the unknown parameter N of a binomial distribution with unknown parameter p, when we reinforce data with an independent sample of a negative-binomial experiment having the same p.


title: Antimatter in the Milky Way

abstract: Recent astronomical observations indicating a strikingly abundant presence of antimatter in the Galaxy, in particular, of anti-stars are reviewed. Long-time earlier theoretical predictions are briefly discussed.


title: Kinetic-scale current sheets in the solar wind at 1 AU: Scale-dependent properties and critical current density

abstract: We present analysis of 17,043 proton kinetic-scale current sheets collected over 124 days of Wind spacecraft measurements in the solar wind at 11 Samples/s magnetic field resolution. The current sheets have thickness $\lambda$ from a few tens to one thousand kilometers with typical value around 100 km or from about 0.1 to 10$\lambda_{p}$ in terms of local proton inertial length $\lambda_{p}$. We found that the current density is larger for smaller scale current sheets, $J_0\approx 6\; {\rm nA/m^2} \cdot (\lambda/100\;{\rm km})^{-0.56}$ , but does not statistically exceed critical value $J_A$ corresponding to the drift between ions and electrons of local Alv\'{e}n speed. The observed trend holds in normalized units, $J_0/J_{A}\approx 0.17\cdot (\lambda/\lambda_{p})^{-0.51}$. The current sheets are statistically force-free with magnetic shear angle correlated with current sheet spatial scale, $\Delta \theta\approx 19^{\circ}\cdot (\lambda/\lambda_{p})^{0.5}$. The observed correlations are consistent with local turbulence being the source of proton kinetic-scale current sheets in the solar wind, while mechanisms limiting the current density remain to be understood.


title: Hydrodynamic Simulations of Electron-capture Supernovae: Progenitor and Dimension Dependence

abstract: We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in \texttt{FLASH} with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with 6 progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same setup using the canonical, spherically-symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy ($\sim1.36$-$1.48\times10^{50}$ erg) and the mass ($\sim0.017$-$0.018M_\odot$) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount ($\lesssim1.8\times10^{-3}~M_\odot$) of low-$Y_e$ ($0.35<Y_e<0.4$) component. The baryonic mass of the protoneutron star is $\sim1.34~M_\odot$ ($\sim1.357~M_\odot$) with the new (n8.8) progenitor models when simulations end at $\sim400$ ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate ONeMg core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.


title: Random cohort effects and age groups dependency structure for mortality modelling and forecasting: Mixed-effects time-series model approach

abstract: There have been significant efforts devoted to solving the longevity risk given that a continuous growth in population ageing has become a severe issue for many developed countries over the past few decades. The Cairns-Blake-Dowd (CBD) model, which incorporates cohort effects parameters in its parsimonious design, is one of the most well-known approaches for mortality modelling at higher ages and longevity risk. This article proposes a novel mixed-effects time-series approach for mortality modelling and forecasting with considerations of age groups dependence and random cohort effects parameters. The proposed model can disclose more mortality data information and provide a natural quantification of the model parameters uncertainties with no pre-specified constraint required for estimating the cohort effects parameters. The abilities of the proposed approach are demonstrated through two applications with empirical male and female mortality data. The proposed approach shows remarkable improvements in terms of forecast accuracy compared to the CBD model in the short-, mid-and long-term forecasting using mortality data of several developed countries in the numerical examples.


title: Elimination (a,b)-trees with fast, durable updates

abstract: Many concurrent dictionary implementations are designed and optimized for read-mostly workloads with uniformly distributed keys, and often perform poorly on update-heavy workloads. In this work, we first present a concurrent (a,b)-tree, the OCC-ABtree, which outperforms its fastest competitor by up to 2x on uniform update-heavy workloads, and is competitive on other workloads. We then turn our attention to skewed update-heavy workloads (which feature many inserts/deletes on the same key) and introduce the Elim-ABtree, which uses a new optimization called publishing elimination. In publishing elimination, concurrent inserts and deletes to a key are reordered to eliminate them. This reduces the number of writes in the data structure. The Elim-ABtree achieves up to 2.5x the performance of its fastest competitor (including the OCC-ABtree). The OCC-ABtree and Elim-ABtree are linearizable. We also introduce durable linearizable versions (for systems with Intel Optane DCPMM non-volatile main memory) that are nearly as fast.


title: Chudnovsky's Conjecture and the stable Harbourne-Huneke containment for general points

abstract: In our previous work with Grifo and H\`a, we showed the stable Harbourne-Huneke containment and Chudnovsky's conjecture for the defining ideal of sufficiently many general points in $\mathbb{P}^N$. In this paper, we establish the conjectures for all remaining cases, and hence, give the affirmative answer to Harbourne-Huneke containment and Chudnovsky's conjecture for any number of general points in $\mathbb{P}^N$ for all $N$. Our new technique is to develop the Cremona reduction process that provides effective lower bounds for the Waldschmidt constant of the defining ideals of generic points in projective spaces.


title: Discrete LQR and ILQR methods based on high order Runge-Kutta methods

abstract: In this paper, discrete linear quadratic regulator (DLQR) and iterative linear quadratic regulator (ILQR) methods based on high-order Runge-Kutta (RK) discretization are proposed for solving linear and nonlinear quadratic optimal control problems respectively. As discovered in [W. Hager, Runge-Kutta method in optimal control and the discrete adjoint system, Numer. Math.,2000, pp. 247-282], direct approach with RK discretization is equivalent with indirect approach based on symplectic partitioned Runge-Kutta (SPRK) integration. In this paper, we will reconstruct this equivalence by the analogue of continuous and discrete dynamic programming. Then, based on the equivalence, we discuss the issue that the internal-stage controls produced by direct approach may have lower order accuracy than the RK method used. We propose order conditions for internal-stage controls and then demonstrate that third or fourth order explicit RK discretization cannot avoid the order reduction phenomenon. To overcome this obstacle, we calculate node control instead of internal-stage controls in DLQR and ILQR methods. And numerical examples will illustrate the validity of our methods. Another advantage of our methods is high computational efficiency which comes from the usage of feedback technique. In this paper, we also demonstrate that ILQR is essentially a quasi-Newton method with linear convergence rate.


title: On gamma matrices of local zeta functions associated with homogeneous cones

abstract: The purpose of this paper is to investigate coefficient matrices of functional equations of zeta functions associated with homogeneous cones, which are given explicitly in the previous paper, in detail. We prove that the coefficient matrix can be decomposed into variable-wise matrices regardless of the choice of homogeneous cones. Moreover, under a suitable condition, we show that the associated zeta functions admit a kind of completion forms.


title: Unknotting twisted knots with Gauss diagram forbidden moves

abstract: Twisted knot theory, introduced by M.O.Bourgoin, is a generalization of virtual knot theory. It is easily shown that any virtual knot can be deformed into a trivial knot by a finite sequence of generalized Reidemeister moves and two "forbidden moves" $F_{1}$ and $F_{2}$. Similarly, we show that any twisted knot also can be deformed into a trivial knot or a trivial knot with a bar by a finite sequence of extended Reidemeister moves and three "forbidden moves" $T_{4}$, $F_{1}$ (or $F_{2}$) and $F_{3}$ (or $F_{4}$) .


title: Higher Frobenius-Schur indicators for semisimple Hopf algebras in positive characteristic

abstract: Let $H$ be a semisimple Hopf algebra over an algebraically closed field $\mathbbm{k}$ of characteristic $p>\dim_{\mathbbm{k}}(H)^{1/2}$. We show that the antipode $S$ of $H$ satisfies the equality $S^2(h)=\mathbf{u}h\mathbf{u}^{-1}$, where $h\in H$, $\mathbf{u}=S(\Lambda_{(2)})\Lambda_{(1)}$ and $\Lambda$ is a nonzero integral of $H$. The formula of $S^2$ enables us to define higher Frobenius-Schur indicators for the Hopf algebra $H$. This generalizes the notions of higher Frobenius-Schur indicators from the case of characteristic 0 to the case of characteristic $p>\dim_{\mathbbm{k}}(H)^{1/2}$. These indicators defined here share some properties with the ones defined over a field of characteristic 0. Especially, all these indicators are gauge invariants for the tensor category Rep$(H)$ of finite dimensional representations of $H$.


title: Entropy Regularized Optimal Transport Independence Criterion

abstract: We introduce an independence criterion based on entropy regularized optimal transport. Our criterion can be used to test for independence between two samples. We establish non-asymptotic bounds for our test statistic and study its statistical behavior under both the null hypothesis and the alternative hypothesis. The theoretical results involve tools from U-process theory and optimal transport theory. We also offer a random feature type approximation for large-scale problems, as well as a differentiable program implementation for deep learning applications. We present experimental results on existing benchmarks for independence testing, illustrating the interest of the proposed criterion to capture both linear and nonlinear dependencies in synthetic data and real data.


title: Shock Reflection in Plane Symmetry

abstract: We consider the problem of shock reflection on a solid wall in plane symmetry for a barotropic fluid. We establish a local in time solution after the point of reflection, thereby determining the state behind the reflected shock. The location of the reflected shock in space time being unknown the problem constitutes a free boundary problem. Furthermore, the quantity describing the thermodynamic state of the fluid along the wall being also unknown, complete data for the problem is only given in a single point, the reflection point.


title: Closed-form meromorphic solutions of some third order boundary layer ordinary differential equations

abstract: We introduce a general third order non-linear autonomous ODE which covers many ODEs coming from boundary layer problems, like the Falkner-Skan equation and the Cheng-Minkowycz equation. Using Wiman-Valiron theory and complex analytic methods recently developed, for the generic cases, it is shown that all their meromorphic solutions must be rational, or rational in one exponential, and then we find all of them explicitly. For a few non-generic cases, some solutions, which are meromorphic or singlevalued, are also obtained. Our results also explain why it is so difficult to obtain new closed-form solutions of the Falkner-Skan equation.


title: Lower Bounds for Regulators of Number Fields in terms of their Discriminants

abstract: We prove inequalities that compare the regulator of a number field with its absolute discriminant. We refine some ideas in Silverman's work in 1984 where such general inequalities are first proven. In order to prove our main theorems, we combine these refinements with the authors' recent results on bounding the product of heights of relative units in a number field extension.


title: Random Anisotropy Magnet at Finite Temperature

abstract: Abstract We present finite-temperature Monte Carlo studies of a 2D random-anisotropy magnet on lattices containing one million spins. The correlated spin-glass state predicted by analytical theories is reproduced in simulations, as are the field-cooled and zero-field-cooled magnetization curves observed in experiments. The orientations of lattice spins begin to freeze when the temperature is lowered. The freezing transition is due to the energy barriers generated by the random anisotropy rather than due to random interactions in conventional spin-glasses. We describe freezing by introducing the time-dependent spin-glass order parameter $q$ and the spin-melting time $\tau_{M}$ defined via $q=\tau_{M}/t$ above freezing, where t is the time of the experiment represented by the number of Monte Carlo steps.


title: Echo state graph neural networks with analogue random resistor arrays

abstract: Recent years have witnessed an unprecedented surge of interest, from social networks to drug discovery, in learning representations of graph-structured data. However, graph neural networks, the machine learning models for handling graph-structured data, face significant challenges when running on conventional digital hardware, including von Neumann bottleneck incurred by physically separated memory and processing units, slowdown of Moore's law due to transistor scaling limit, and expensive training cost. Here we present a novel hardware-software co-design, the random resistor array-based echo state graph neural network, which addresses these challenges. The random resistor arrays not only harness low-cost, nanoscale and stackable resistors for highly efficient in-memory computing using simple physical laws, but also leverage the intrinsic stochasticity of dielectric breakdown to implement random projections in hardware for an echo state network that effectively minimizes the training cost thanks to its fixed and random weights. The system demonstrates state-of-the-art performance on both graph classification using the MUTAG and COLLAB datasets and node classification using the CORA dataset, achieving 34.2x, 93.2x, and 570.4x improvement of energy efficiency and 98.27%, 99.46%, and 95.12% reduction of training cost compared to conventional graph learning on digital hardware, respectively, which may pave the way for the next generation AI system for graph learning.


title: BP-Net: Cuff-less, Calibration-free, and Non-invasive Blood Pressure Estimation via a Generic Deep Convolutional Architecture

abstract: Objective: The paper focuses on development of robust and accurate processing solutions for continuous and cuff-less blood pressure (BP) monitoring. In this regard, a robust deep learning-based framework is proposed for computation of low latency, continuous, and calibration-free upper and lower bounds on the systolic and diastolic BP. Method: Referred to as the BP-Net, the proposed framework is a novel convolutional architecture that provides longer effective memory while achieving superior performance due to incorporation of casual dialated convolutions and residual connections. To utilize the real potential of deep learning in extraction of intrinsic features (deep features) and enhance the long-term robustness, the BP-Net uses raw Electrocardiograph (ECG) and Photoplethysmograph (PPG) signals without extraction of any form of hand-crafted features as it is common in existing solutions. Results: By capitalizing on the fact that datasets used in recent literature are not unified and properly defined, a benchmark dataset is constructed from the MIMIC-I and MIMIC-III databases obtained from PhysioNet. The proposed BP-Net is evaluated based on this benchmark dataset demonstrating promising performance and shows superior generalizable capacity. Conclusion: The proposed BP-Net architecture is more accurate than canonical recurrent networks and enhances the long-term robustness of the BP estimation task. Significance: The proposed BP-Net architecture addresses key drawbacks of existing BP estimation solutions, i.e., relying heavily on extraction of hand-crafted features, such as pulse arrival time (PAT), and; Lack of robustness. Finally, the constructed BP-Net dataset provides a unified base for evaluation and comparison of deep learning-based BP estimation algorithms.


title: ViNMT: Neural Machine Translation Toolkit

abstract: We present an open-source toolkit for neural machine translation (NMT). The new toolkit is mainly based on vaulted Transformer (Vaswani et al., 2017) along with many other improvements detailed below, in order to create a self-contained, simple to use, consistent and comprehensive framework for Machine Translation tasks of various domains. It is tooled to support both bilingual and multilingual translation tasks, starting from building the model from respective corpora, to inferring new predictions or packaging the model to serving-capable JIT format.


title: Power Under Multiplicity Project (PUMP): Estimating Power, Minimum Detectable Effect Size, and Sample Size When Adjusting for Multiple Outcomes in Multi-level Experiments

abstract: For randomized controlled trials (RCTs) with a single intervention being measured on multiple outcomes, researchers often apply a multiple testing procedure (such as Bonferroni or Benjamini-Hochberg) to adjust $p$-values. Such an adjustment reduces the likelihood of spurious findings, but also changes the statistical power, sometimes substantially, which reduces the probability of detecting effects when they do exist. However, this consideration is frequently ignored in typical power analyses, as existing tools do not easily accommodate the use of multiple testing procedures. We introduce the PUMP R package as a tool for analysts to estimate statistical power, minimum detectable effect size, and sample size requirements for multi-level RCTs with multiple outcomes. Multiple outcomes are accounted for in two ways. First, power estimates from PUMP properly account for the adjustment in $p$-values from applying a multiple testing procedure. Second, as researchers change their focus from one outcome to multiple outcomes, different definitions of statistical power emerge. PUMP allows researchers to consider a variety of definitions of power, as some may be more appropriate for the goals of their study. The package estimates power for frequentist multi-level mixed effects models, and supports a variety of commonly-used RCT designs and models and multiple testing procedures. In addition to the main functionality of estimating power, minimum detectable effect size, and sample size requirements, the package allows the user to easily explore sensitivity of these quantities to changes in underlying assumptions.


title: Has the impact flux of small and large asteroids varied through time on Mars, the Earth and the Moon?

abstract: The impact flux over the last 3 Ga in the inner Solar System is commonly assumed to be constant through time. However, asteroid break-up events in the main belt may have been responsible for cratering spikes over the last ~2 Ga on the Earth-Moon system. We investigate here the possible variations of the size frequency distributions of impactors from the record of small craters of 521 martian impact craters larger than 20 km in diameter. We show that 49 craters (out of the 521) correspond to the complete crater population of this size formed over the last 600 Ma. Our results on Mars show that the flux of both small (> 5 m) and large asteroids (> 1 km) are coupled, does not vary between each other over the last 600 Ma. Existing data sets for large craters on the Earth and the Moon are analyzed and compared to our results on Mars. On Earth, we infer the formation location of a set of impact craters thanks to plate tectonic reconstruction and show that a cluster of craters formed during the Ordovician period, about 470 Ma ago, appears to be a preservation bias. On the Moon, the late increase seen in the crater age signal can be due to the uncertain calibration method used to date those impacts (i.e. rock abundance in lunar impact ejecta), and other calibrations are consistent with a constant crater production rate. We conclude to a coupling of the crater production rate between kilometer-size craters and down to ~100 m in diameter in the inner Solar System. This is consistent with the traditional model for delivering asteroids to planet-crossing orbits: the Yarkovsky effect slowly pushes the large debris from asteroid break-ups towards orbital resonances while smaller debris are grinded through collisional cascades. This suggests that the influence of past asteroid break-ups in the cratering rate for D > 100 m is limited or inexistent.


title: Learned Coarse Models for Efficient Turbulence Simulation

abstract: Turbulence simulation with classical numerical solvers requires high-resolution grids to accurately resolve dynamics. Here we train learned simulators at low spatial and temporal resolutions to capture turbulent dynamics generated at high resolution. We show that our proposed model can simulate turbulent dynamics more accurately than classical numerical solvers at the comparably low resolutions across various scientifically relevant metrics. Our model is trained end-to-end from data and is capable of learning a range of challenging chaotic and turbulent dynamics at low resolution, including trajectories generated by the state-of-the-art Athena++ engine. We show that our simpler, general-purpose architecture outperforms various more specialized, turbulence-specific architectures from the learned turbulence simulation literature. In general, we see that learned simulators yield unstable trajectories; however, we show that tuning training noise and temporal downsampling solves this problem. We also find that while generalization beyond the training distribution is a challenge for learned models, training noise, added loss constraints, and dataset augmentation can help. Broadly, we conclude that our learned simulator outperforms traditional solvers run on coarser grids, and emphasize that simple design choices can offer stability and robust generalization.


title: The lower moments of nucleon structure functions in lattice QCD with physical quark masses

abstract: We present results for the nucleon structure functions and form factors obtained from 2+1 flavor lattice QCD with physical light quark masses ($m_{\pi}=135$ MeV) in a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared ${\mathscr{O}}(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $\beta=1.82$ and $2.00$ corresponding to lattice spacings of $0.085$ fm and $0.064$ fm respectively. The lower moments of structure functions, $\langle x \rangle_{u-d}$ and $\langle x \rangle_{\Delta u - \Delta d}$ given by the twist-2 operators being properly renormalized, are evaluated in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV only at $\beta=1.82$, since the renormalization factors at $\beta=2.00$ have not yet determined nonperturbatively in the RI/MOM scheme. Instead, at two lattice spacings, we evaluate appropriate ratios of $g_{A}/g_{V}$ and $\langle x \rangle_{u-d}/\langle x \rangle_{\Delta u -\Delta d}$, which are not renormalized in the continuum limit. These quantities thus can be directly compared with the experimental data without the renormalization.


title: Machine Learning Application Development: Practitioners' Insights

abstract: Nowadays, intelligent systems and services are getting increasingly popular as they provide data-driven solutions to diverse real-world problems, thanks to recent breakthroughs in Artificial Intelligence (AI) and Machine Learning (ML). However, machine learning meets software engineering not only with promising potentials but also with some inherent challenges. Despite some recent research efforts, we still do not have a clear understanding of the challenges of developing ML-based applications and the current industry practices. Moreover, it is unclear where software engineering researchers should focus their efforts to better support ML application developers. In this paper, we report about a survey that aimed to understand the challenges and best practices of ML application development. We synthesize the results obtained from 80 practitioners (with diverse skills, experience, and application domains) into 17 findings; outlining challenges and best practices for ML application development. Practitioners involved in the development of ML-based software systems can leverage the summarized best practices to improve the quality of their system. We hope that the reported challenges will inform the research community about topics that need to be investigated to improve the engineering process and the quality of ML-based applications.


title: Data-Free Knowledge Transfer: A Survey

abstract: In the last decade, many deep learning models have been well trained and made a great success in various fields of machine intelligence, especially for computer vision and natural language processing. To better leverage the potential of these well-trained models in intra-domain or cross-domain transfer learning situations, knowledge distillation (KD) and domain adaptation (DA) are proposed and become research highlights. They both aim to transfer useful information from a well-trained model with original training data. However, the original data is not always available in many cases due to privacy, copyright or confidentiality. Recently, the data-free knowledge transfer paradigm has attracted appealing attention as it deals with distilling valuable knowledge from well-trained models without requiring to access to the training data. In particular, it mainly consists of the data-free knowledge distillation (DFKD) and source data-free domain adaptation (SFDA). On the one hand, DFKD aims to transfer the intra-domain knowledge of original data from a cumbersome teacher network to a compact student network for model compression and efficient inference. On the other hand, the goal of SFDA is to reuse the cross-domain knowledge stored in a well-trained source model and adapt it to a target domain. In this paper, we provide a comprehensive survey on data-free knowledge transfer from the perspectives of knowledge distillation and unsupervised domain adaptation, to help readers have a better understanding of the current research status and ideas. Applications and challenges of the two areas are briefly reviewed, respectively. Furthermore, we provide some insights to the subject of future research.


title: Counting substructures and eigenvalues II: quadrilaterals

abstract: Let $G$ be a graph and $\lambda(G)$ be the spectral radius of $G$. A previous result due to Nikiforov [Linear Algebra Appl., 2009] in spectral graph theory asserted that every graph $G$ on $m\geq 10$ edges contains a 4-cycle if $\lambda(G)>\sqrt{m}$. Define $f(m)$ to be the minimum number of copies of 4-cycles in such a graph. A consequence of a recent theorem due to Zhai et al. [European J. Combin., 2021] shows that $f(m)=\Omega(m)$. In this article, by somewhat different techniques, we prove that $f(m)=\Theta(m^2)$. We left the solution to $\lim\limits_{m\rightarrow \infty} \frac{f(m)}{m^2}$ as a problem, and also mention other ones for further study.


title: What is Event Knowledge Graph: A Survey

abstract: Besides entity-centric knowledge, usually organized as Knowledge Graph (KG), events are also an essential kind of knowledge in the world, which trigger the spring up of event-centric knowledge representation form like Event KG (EKG). It plays an increasingly important role in many downstream applications, such as search, question-answering, recommendation, financial quantitative investments, and text generation. This paper provides a comprehensive survey of EKG from history, ontology, instance, and application views. Specifically, to characterize EKG thoroughly, we focus on its history, definition, schema induction, acquisition, related representative graphs/systems, and applications. The development processes and trends are studied therein. We further summarize prospective directions to facilitate future research on EKG.


title: Efficient Channel Estimation for RIS-Aided MIMO Communications with Unitary Approximate Message Passing

abstract: Reconfigurable intelligent surface (RIS) is very promising for wireless networks to achieve high energy efficiency, extended coverage, improved capacity, massive connectivity, etc. To unleash the full potentials of RIS-aided communications, acquiring accurate channel state information is crucial, which however is very challenging. For RIS-aided multiple-input and multiple-output (MIMO) communications, the existing channel estimation methods have computational complexity growing rapidly with the number of RIS units $N$ (e.g., in the order of $N^2$ or $N^3$) and/or have special requirements on the matrices involved (e.g., the matrices need to be sparse for algorithm convergence to achieve satisfactory performance), which hinder their applications. In this work, instead of using the conventional signal model in the literature, we derive a new signal model obtained through proper vectorization and reduction operations. Then, leveraging the unitary approximate message passing (UAMP), we develop a more efficient channel estimator that has complexity linear with $N$ and does not have special requirements on the relevant matrices, thanks to the robustness of UAMP. These facilitate the applications of the proposed algorithm to a general RIS-aided MIMO system with a larger $N$. Moreover, extensive numerical results show that the proposed estimator delivers much better performance and/or requires significantly less number of training symbols, thereby leading to notable reductions in both training overhead and latency.


title: Energy-Aware Multi-Robot Task Allocation in Persistent Tasks

abstract: The applicability of the swarm robots to perform foraging tasks is inspired by their compact size and cost. A considerable amount of energy is required to perform such tasks, especially if the tasks are continuous and/or repetitive. Real-world situations in which robots perform tasks continuously while staying alive (survivability) and maximizing production (performance) require energy awareness. This paper proposes an energy-conscious distributed task allocation algorithm to solve continuous tasks (e.g., unlimited foraging) for cooperative robots to achieve highly effective missions. We consider efficiency as a function of the energy consumed by the robot during exploration and collection when food is returned to the collection bin. Finally, the proposed energy-efficient algorithm minimizes the total transit time to the charging station and time consumed while recharging and maximizes the robot's lifetime to perform maximum tasks to enhance the overall efficiency of collaborative robots. We evaluated the proposed solution against a typical greedy benchmarking strategy (assigning the closest collection bin to the available robot and recharging the robot at maximum) for efficiency and performance in various scenarios. The proposed approach significantly improved performance and efficiency over the baseline approach.


title: ERNIE-ViLG: Unified Generative Pre-training for Bidirectional Vision-Language Generation

abstract: Conventional methods for the image-text generation tasks mainly tackle the naturally bidirectional generation tasks separately, focusing on designing task-specific frameworks to improve the quality and fidelity of the generated samples. Recently, Vision-Language Pre-training models have greatly improved the performance of the image-to-text generation tasks, but large-scale pre-training models for text-to-image synthesis task are still under-developed. In this paper, we propose ERNIE-ViLG, a unified generative pre-training framework for bidirectional image-text generation with transformer model. Based on the image quantization models, we formulate both image generation and text generation as autoregressive generative tasks conditioned on the text/image input. The bidirectional image-text generative modeling eases the semantic alignments across vision and language. For the text-to-image generation process, we further propose an end-to-end training method to jointly learn the visual sequence generator and the image reconstructor. To explore the landscape of large-scale pre-training for bidirectional text-image generation, we train a 10-billion parameter ERNIE-ViLG model on a large-scale dataset of 145 million (Chinese) image-text pairs which achieves state-of-the-art performance for both text-to-image and image-to-text tasks, obtaining an FID of 7.9 on MS-COCO for text-to-image synthesis and best results on COCO-CN and AIC-ICC for image captioning.


title: A simple method for measuring inequality

abstract: To simultaneously overcome the limitation of the Gini index in that it is less sensitive to inequality at the tails of income distribution and the limitation of the inter-decile ratios that ignore inequality in the middle of income distribution, an inequality index is introduced. It comprises three indicators, namely, the Gini index, the income share held by the top 10%, and the income share held by the bottom 10%. The data from the World Bank database and the Organization for Economic Co-operation and Development Income Distribution Database between 2005 and 2015 are used to demonstrate how the inequality index works. The results show that it can distinguish income inequality among countries that share the same Gini index but have different income gaps between the top 10% and the bottom 10%. It could also distinguish income inequality among countries that have the same ratio of income share held by the top 10% to income share held by the bottom 10% but differ in the values of the Gini index. In addition, the inequality index could capture the dynamics where the Gini index of a country is stable over time but the ratio of income share of the top 10% to income share of the bottom 10% is increasing. Furthermore, the inequality index could be applied to other scientific disciplines as a measure of statistical heterogeneity and for size distributions of any non-negative quantities.


title: Modelling of Bi-directional Spatio-Temporal Dependence and Users' Dynamic Preferences for Missing POI Check-in Identification

abstract: Human mobility data accumulated from Point-of-Interest (POI) check-ins provides great opportunity for user behavior understanding. However, data quality issues (e.g., geolocation information missing, unreal check-ins, data sparsity) in real-life mobility data limit the effectiveness of existing POI-oriented studies, e.g., POI recommendation and location prediction, when applied to real applications. To this end, in this paper, we develop a model, named Bi-STDDP, which can integrate bi-directional spatio-temporal dependence and users' dynamic preferences, to identify the missing POI check-in where a user has visited at a specific time. Specifically, we first utilize bi-directional global spatial and local temporal information of POIs to capture the complex dependence relationships. Then, target temporal pattern in combination with user and POI information are fed into a multi-layer network to capture users' dynamic preferences. Moreover, the dynamic preferences are transformed into the same space as the dependence relationships to form the final model. Finally, the proposed model is evaluated on three large-scale real-world datasets and the results demonstrate significant improvements of our model compared with state-of-the-art methods. Also, it is worth noting that the proposed model can be naturally extended to address POI recommendation and location prediction tasks with competitive performances.


title: A new class of differential quasivariational inequalities with an application to a quasistatic viscoelastic frictional contact problem

abstract: The overarching goal of this paper is to introduce and investigate a new nonlinear system driven by a nonlinear differential equation, a history-dependent quasivariational inequality, and a parabolic variational inequality in Banach spaces. Such a system can be used to model quasistatic frictional contact problems for viscoelastic materials with long memory, damage and wear. By using the Banach fixed point theorem, we prove an existence and uniqueness theorem of solution for such a system under some mild conditions. As a novel application, we obtain a unique solvability of a quasistatic viscoelastic frictional contact problem with long memory, damage and wear.


title: Distributed Random Reshuffling over Networks

abstract: In this paper, we consider distributed optimization problems where $n$ agents, each possessing a local cost function, collaboratively minimize the average of the local cost functions over a connected network. To solve the problem, we propose a distributed random reshuffling (D-RR) algorithm that invokes the random reshuffling (RR) update in each agent. We show that D-RR inherits favorable characteristics of RR for both smooth strongly convex and smooth nonconvex objective functions. In particular, for smooth strongly convex objective functions, D-RR achieves $\mathcal{O}(1/T^2)$ rate of convergence (where $T$ counts epoch number) in terms of the squared distance between the iterate and the global minimizer. When the objective function is assumed to be smooth nonconvex, we show that D-RR drives the squared norm of gradient to $0$ at a rate of $\mathcal{O}(1/T^{2/3})$. These convergence results match those of centralized RR (up to constant factors) and outperform the distributed stochastic gradient descent (DSGD) algorithm if we run a relatively large number of epochs. Finally, we conduct a set of numerical experiments to illustrate the efficiency of the proposed D-RR method on both strongly convex and nonconvex distributed optimization problems.


title: Tri-Hexagonal charge order in kagome metal CsV$_{3}$Sb$_{5}$ revealed by $^{121}$Sb NQR

abstract: We report $^{121}$Sb nuclear quadrupole resonance (NQR) measurements on kagome superconductor CsV$_3$Sb$_5$ with $T_{\rm c}=2.5$ K. $^{121}$Sb NQR spectra split after a charge density wave (CDW) transition at $94$ K, which demonstrates a commensurate CDW state. The coexistence of the high temperature phase and the CDW phase between $91$ K and $94$ K manifests that it is a first order phase transition. The CDW order exhibits Tri-Hexagonal deformation with a lateral shift between the adjacent kagome layers, which is consistent with $2 \times 2 \times 2$ superlattice modulation. The superconducting state coexists with CDW order and shows a conventional s-wave behavior in the bulk state.


title: Homogenization for polynomial optimization with unbounded sets

abstract: This paper considers polynomial optimization with unbounded sets. We give a homogenization formulation and propose a hierarchy of Moment-SOS relaxations to solve it. Under the assumptions that the feasible set is closed at infinity and the ideal of homogenized equality constraining polynomials is real radical, we show that this hierarchy of Moment-SOS relaxations has finite convergence, if some optimality conditions (i.e., the linear independence constraint qualification, strict complementarity and second order sufficiency conditions) hold at every minimizer, including the one at infinity. Moreover, we prove extended versions of Putinar-Vasilescu type Positivstellensatz for polynomials that are nonnegative on unbounded sets. The classical Moment-SOS hierarchy with denominators is also studied. In particular, we give a positive answer to a conjecture of Mai, Lasserre and Magron in their recent work.


title: Domain Adaptation with Category Attention Network for Deep Sentiment Analysis

abstract: Domain adaptation tasks such as cross-domain sentiment classification aim to utilize existing labeled data in the source domain and unlabeled or few labeled data in the target domain to improve the performance in the target domain via reducing the shift between the data distributions. Existing cross-domain sentiment classification methods need to distinguish pivots, i.e., the domain-shared sentiment words, and non-pivots, i.e., the domain-specific sentiment words, for excellent adaptation performance. In this paper, we first design a Category Attention Network (CAN), and then propose a model named CAN-CNN to integrate CAN and a Convolutional Neural Network (CNN). On the one hand, the model regards pivots and non-pivots as unified category attribute words and can automatically capture them to improve the domain adaptation performance; on the other hand, the model makes an attempt at interpretability to learn the transferred category attribute words. Specifically, the optimization objective of our model has three different components: 1) the supervised classification loss; 2) the distributions loss of category feature weights; 3) the domain invariance loss. Finally, the proposed model is evaluated on three public sentiment analysis datasets and the results demonstrate that CAN-CNN can outperform other various baseline methods.


title: A simple method for estimating the Lorenz curve

abstract: Given many popular functional forms for the Lorenz curve do not have a closed-form expression for the Gini index and no study has utilized the observed Gini index to estimate parameter(s) associated with the corresponding parametric functional form, a simple method for estimating the Lorenz curve is introduced. It utilizes 3 indicators, namely, the Gini index and the income shares of the bottom and the top in order to calculate the values of parameters associated with the specified functional form which has a closed-form expression for the Gini index. No error minimization technique is required in order to estimate the Lorenz curve. The data on the Gini index and the income shares of 4 countries that have different level of income inequality, economic, sociological, and regional backgrounds from the United Nations University-World Income Inequality Database are used to illustrate how the simple method works. The overall results indicate that the estimated Lorenz curves fit the actual observations practically well. This simple method could be useful in the situation where the availability of data on income distribution is low. However, if more data on income distribution are available, this study shows that the specified functional form could be used to directly estimate the Lorenz curve. Moreover, the estimated values of the Gini index calculated based on the specified functional form are virtually identical to their actual observations.


title: Neural Hierarchical Factorization Machines for User's Event Sequence Analysis

abstract: Many prediction tasks of real-world applications need to model multi-order feature interactions in user's event sequence for better detection performance. However, existing popular solutions usually suffer two key issues: 1) only focusing on feature interactions and failing to capture the sequence influence; 2) only focusing on sequence information, but ignoring internal feature relations of each event, thus failing to extract a better event representation. In this paper, we consider a two-level structure for capturing the hierarchical information over user's event sequence: 1) learning effective feature interactions based event representation; 2) modeling the sequence representation of user's historical events. Experimental results on both industrial and public datasets clearly demonstrate that our model achieves significantly better performance compared with state-of-the-art baselines.


title: Controlling thermoelectric, heat, and energy currents through a quantum dot in sequential and cotunneling Coulomb-blockade regimes

abstract: Thermal transport through a Coulomb-blockade quantum dot (QD) coupled to two metallic leads is studied using five different approaches to the master equation in which sequential and coutuneling terms are taken into account. In the presence of intradot Coulomb interactions, a plateau in the thermo-particle, the heat, and the energy currents is seen. The current plateau diminishes at a high thermal bias between the leads. It is shown that the Pauli, the Redfield, the Lindblad-type equation with first order tunneling rates, and first-order von-Neumann master equations give very similar thermal transport indicating the conservation of coherency in the electron transport in sequential tunneling between the QD and leads. In contrast, the thermal transport is suppressed when coutuneling processes are taken into account via a second-order von-Neumann master equation. The consideration of second order effects with respect to the QD-leads coupling brings in a wealth of virtual processes at the contact to each lead. These virtual processes directly weaken the effects of the contribution of the first order direct processes to the overall transport, and introduce important other aspects of the transport, as level broadening, energy shifts, and lifetimes in the time-domain. As a result the current plateau formed via the Coulomb interaction diminishes, when second order and cotunneling processes are considered.


title: Macroeconomic and financial management in an uncertain world: What can we learn from complexity science?

abstract: This paper discusses serious drawbacks of existing knowledge in macroeconomics and finance in explaining and predicting economic and financial phenomena. Complexity science is proposed as an alternative approach to be used in order to better understand how economy and financial market work. This paper argues that understanding characteristics of complex system could greatly benefit financial analysts, financial regulators, as well as macroeconomic policy makers.

categories: econ.GN nlin.AO q-fin.EC q-fin.GN stat.ME


title: High Precision Calculations of the Higgs Boson Mass

abstract: In this paper we are going to review the status of the computations of the perturbative quantum corrections to the Higgs boson mass in the Standard Model and in its supersymmetric extensions. In particular, supersymmetric theories require a very accurate computation of the Higgs boson mass, which includes corrections up to even three-loop level, since their predictions are limited by theoretical uncertainties. A discussion about these uncertainties in the context of the minimal and next to minimal supersymmetric standard models is included.


title: Uncovering migration systems through spatio-temporal tensor co-clustering

abstract: A central problem in the study of human mobility is that of migration systems. Typically, migration systems are defined as a set of relatively stable movements of people between two or more locations over time. While these emergent systems are expected to vary over time, they ideally contain a stable underlying structure that could be discovered empirically. There have been some notable attempts to formally or informally define migration systems, however they have been limited by being hard to operationalize, and by defining migration systems in ways that ignore origin/destination aspects and/or fail to account for migration dynamics. In this work we propose a novel method, spatio-temporal (ST) tensor co-clustering, stemming from signal processing and machine learning theory. To demonstrate its effectiveness for describing stable migration systems we focus on domestic migration between counties in the US from 1990-2018. Relevant data for this period has been made available through the US Internal Revenue Service. Specifically, we concentrate on three illustrative case studies: (i) US Metropolitan Areas, (ii) the state of California, and (iii) Louisiana, focusing on detecting exogenous events such as Hurricane Katrina in 2005. Finally, we conclude with discussion and limitations of this approach.


title: On the three graph invariants related to matching of finite simple graphs

abstract: Let $G$ be a finite simple graph on the vertex set $V(G)$ and let $\text{ind-match}(G)$, $\text{min-match}(G)$ and $\text{match}(G)$ denote the induced matching number, the minimum matching number and the matching number of $G$, respectively. It is known that the inequalities $\text{ind-match}(G) \leq \text{min-match}(G) \leq \text{match}(G) \leq 2\text{min-match}(G)$ and $\text{match}(G) \leq \left\lfloor |V(G)|/2 \right\rfloor$ hold in general. In the present paper, we determine the possible tuples $(p, q, r, n)$ with $\text{ind-match}(G) = p$, $\text{min-match}(G) = q$, $\text{match}(G) = r$ and $|V(G)| = n$ arising from connected simple graphs. As an application of this result, we also determine the possible tuples $(p', q, r, n)$ with ${\rm{reg}}(G) = p'$, $\text{min-match}(G) = q$, $\text{match}(G) = r$ and $|V(G)| = n$ arising from connected simple graphs, where $I(G)$ is the edge ideal of $G$ and ${\rm{reg}}(G) = {\rm{reg}}(K[V(G)]/I(G))$ is the Castelnuovo--Mumford regularity of the quotient ring $K[V(G)]/I(G)$.


title: Energetic Formulation of Large-Deformation Poroelasticity

abstract: The modeling of coupled fluid transport and deformation in a porous medium is essential to predict the various geomechanical process such as CO2 sequestration, hydraulic fracturing, and so on. Current applications of interest, for instance, that include fracturing or damage of the solid phase, require a nonlinear description of the large deformations that can occur. This paper presents a variational energy-based continuum mechanics framework to model large-deformation poroelasticity. The approach begins from the total free energy density that is additively composed of the free energy of the components. A variational procedure then provides the balance of momentum, fluid transport balance, and pressure relations. A numerical approach based on finite elements is applied to analyze the behavior of saturated and unsaturated porous media using a nonlinear constitutive model for the solid skeleton. Examples studied include the Terzaghi and Mandel problems; a gas-liquid phase-changing fluid; multiple immiscible gases; and unsaturated systems where we model injection of fluid into soil. The proposed variational approach can potentially have advantages for numerical methods as well as for combining with data-driven models in a Bayesian framework.


title: CSformer: Bridging Convolution and Transformer for Compressive Sensing

abstract: Convolution neural networks (CNNs) have succeeded in compressive image sensing. However, due to the inductive bias of locality and weight sharing, the convolution operations demonstrate the intrinsic limitations in modeling the long-range dependency. Transformer, designed initially as a sequence-to-sequence model, excels at capturing global contexts due to the self-attention-based architectures even though it may be equipped with limited localization abilities. This paper proposes CSformer, a hybrid framework that integrates the advantages of leveraging both detailed spatial information from CNN and the global context provided by transformer for enhanced representation learning. The proposed approach is an end-to-end compressive image sensing method, composed of adaptive sampling and recovery. In the sampling module, images are measured block-by-block by the learned sampling matrix. In the reconstruction stage, the measurement is projected into dual stems. One is the CNN stem for modeling the neighborhood relationships by convolution, and the other is the transformer stem for adopting global self-attention mechanism. The dual branches structure is concurrent, and the local features and global representations are fused under different resolutions to maximize the complementary of features. Furthermore, we explore a progressive strategy and window-based transformer block to reduce the parameter and computational complexity. The experimental results demonstrate the effectiveness of the dedicated transformer-based architecture for compressive sensing, which achieves superior performance compared to state-of-the-art methods on different datasets.


title: BatchLens: A Visualization Approach for Analyzing Batch Jobs in Cloud Systems

abstract: Cloud systems are becoming increasingly powerful and complex. It is highly challenging to identify anomalous execution behaviors and pinpoint problems by examining the overwhelming intermediate results/states in complex application workflows. Domain scientists urgently need a friendly and functional interface to understand the quality of the computing services and the performance of their applications in real time. To meet these needs, we explore data generated by job schedulers and investigate general performance metrics (e.g., utilization of CPU, memory and disk I/O). Specifically, we propose an interactive visual analytics approach, BatchLens, to provide both providers and users of cloud service with an intuitive and effective way to explore the status of system batch jobs and help them conduct root-cause analysis of anomalous behaviors in batch jobs. We demonstrate the effectiveness of BatchLens through a case study on the public Alibaba bench workload trace datasets.


title: Spin-injection-generated shock waves and solitons in a ferromagnetic thin film: the spin piston problem

abstract: The unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction (expansion) waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with "contact" spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wavetrain. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.


title: Systematic dispersion compensation for spectral domain optical coherence tomography using time-frequency analysis and iterative optimization for iridocorneal angle imaging

abstract: Dispersion is a common phenomenon in optics due to the frequency dependence of the refractive index in polychromatic light. This issue, if left untreated in optical coherence tomography (OCT) imaging, leads to signal broadening of the coherence length and deterioration of the axial resolution. We report a new numeric method for the systematic dispersion compensation in a spectral-domain (SD) OCT for imaging the iridocorneal angle of human cadaver eyes. The dispersion compensation for our OCT system is calculated by an automated iterative process that minimizes the wavenumber-dependent variance of the ridge extracted from the energy distribution of a mirror's spectral interferogram using Short-Time Fourier Transform (STFT) Time-Frequency Analysis (TFA). The average axial resolution of 2.7 um in air was achieved at a range of depths up to 2 mm. Compensated OCT images of the iridocorneal angle in human cadaver eyes were much clearer than non-compensated images. We demonstrate the feasibility, effectiveness, and robustness of the proposed method for dispersion compensation in an SD-OCT by evaluating both the mirror and human cadaver eye measurements. We also verified that our imaging system is able to visualize the iridocorneal angle details, such as trabecular meshwork (TM), Schlemm's canal (SC), and collector channels (CCs), which are important ocular outflow structures and play a crucial role in glaucoma managements.


title: SimSR: Simple Distance-based State Representation for Deep Reinforcement Learning

abstract: This work explores how to learn robust and generalizable state representation from image-based observations with deep reinforcement learning methods. Addressing the computational complexity, stringent assumptions and representation collapse challenges in existing work of bisimulation metric, we devise Simple State Representation (SimSR) operator. SimSR enables us to design a stochastic approximation method that can practically learn the mapping functions (encoders) from observations to latent representation space. In addition to the theoretical analysis and comparison with the existing work, we experimented and compared our work with recent state-of-the-art solutions in visual MuJoCo tasks. The results shows that our model generally achieves better performance and has better robustness and good generalization.


title: An Intelligent Self-driving Truck System For Highway Transportation

abstract: Recently, there have been many advances in autonomous driving society, attracting a lot of attention from academia and industry. However, existing works mainly focus on cars, extra development is still required for self-driving truck algorithms and models. In this paper, we introduce an intelligent self-driving truck system. Our presented system consists of three main components, 1) a realistic traffic simulation module for generating realistic traffic flow in testing scenarios, 2) a high-fidelity truck model which is designed and evaluated for mimicking real truck response in real-world deployment, 3) an intelligent planning module with learning-based decision making algorithm and multi-mode trajectory planner, taking into account the truck's constraints, road slope changes, and the surrounding traffic flow. We provide quantitative evaluations for each component individually to demonstrate the fidelity and performance of each part. We also deploy our proposed system on a real truck and conduct real world experiments which shows our system's capacity of mitigating sim-to-real gap. Our code is available at https://github.com/InceptioResearch/IITS


title: Properties beyond mass for unresolved haloes across redshift and cosmology using correlations with local halo environment

abstract: The structural and dynamic properties of the dark matter halos, though an important ingredient in understanding large-scale structure formation, require more conservative particle resolution than those required by halo mass alone in a simulation. This reduces the parameter space of the simulations, more severely for high-redshift and large-volume mocks which are required by the next-generation large sky surveys. Here, we incorporate redshift and cosmology dependence into an algorithm that assigns accurate halo properties such as concentration, spin, velocity, and spatial distribution to the sub-resolution haloes in a simulation. By focusing on getting the right correlations with halo mass and local tidal anisotropy $\alpha$ measured at $4 \times$ halo radius, our method will also recover the correlations of these small scale structural properties with the large-scale environment, i.e., the halo assembly bias at all scales greater than $5 \times$ halo radius. We find that the distribution of halo properties is universal with redshift and cosmology. By applying the algorithm to a large volume simulation $(600 h^{-1}{\rm Mpc})$, we can access the $30-500$ particle haloes, thus gaining an order of magnitude in halo mass and two to three orders of magnitude in number density at $z=2-4$. This technique reduces the cost of mocks required for the estimation of covariance matrices, weak lensing studies, or any large-scale clustering analysis with less massive haloes.


title: Space Charge Driven Emittance Growth and the Effect of Octupoles in IOTA

abstract: The Integrable Optics Test Accelerator (IOTA) at Fermilab is a small machine dedicated to a broad frontier accelerator physics program. An important aspect of this program is to investigate the potential benefits of the resonance free tune spread achievable with integrable optics to store and accelerate high intensity proton beams for which space charge is significant. In this context, a good understanding of proton beam emittance growth and particle loss mechanisms is essential. Assuming nominal design parameters, simulations show that for a bunched beam, the bulk of emittance growth takes place immediately following injection, typically within tens of turns. We attempt to account for this growth using a simplified RMS mismatch theory; some of its limitations and possible improvements are briefly discussed. We then compare theoretical predictions to simulations performed using the PIC code pyORBIT. Further exploring ways to mitigate emittance growth and reduce particle loss, we compare two beam matching strategies: (1) matching at the injection point (2) matching at the center of the nonlinear (octupole) insertion region where $\beta_x = \beta_y$. To observe how nonlinearity affects emittance growth and whether it dominates growth due to mismatch, we track two different distributions. Finally, we explore the potential of using octupoles in a quasi-integrable configuration to mitigate growth using a variety of initial distributions both at reduced and full intensities.


title: Ferromagnetic Enhancement in LaMnO3 Films with Release and Flexure

abstract: A variety of novel phenomena and functionalities emerge from lowering the dimensionality of materials and enriching the degrees of freedom in modulation. In this work, it is found that the saturation magnetization of LaMnO3 (LMO) films is largely enhanced by 56% after releasing from a brand-new phase of tetragonal strontium aluminate buffer layer, and is significantly increased by 92% with bending films to a curvature of 1 mm-1 using a water-assisted direct-transferring method. Meanwhile, the Curie temperature of LMO films has been improved by 13 K. High-resolution spherical aberration-corrected scanning transmission electron microscopy and first-principles calculations unambiguously demonstrate that the enhanced ferromagnetism is attributed to the strengthened Mn-O-Mn super-exchange interactions from the augmented characteristics of the unconventional P21/n structure caused by the out-of-plane lattice shrinking after strain releasing and increased flexure degree of freestanding LMO films. This work paves a way to achieve large-scale and crack-and-wrinkle-free freestanding films of oxides with largely improved functionalities.


title: Gorenstein braid cones and crepant resolutions

abstract: To any poset $P$, we associate a convex cone called a braid cone. We also associate a fan and study the toric varieties the cone and fan define. The fan always defines a smooth toric variety $X_P$, while the toric variety $U_P$ of the cone may be singular. We show that $X_{P} \dashrightarrow U_{P}$ is a crepant resolution of singularities if and only if $P$ is bounded. Next, we aim to determine when $U_P$ is Gorenstein or $\mathbb{Q}$-Gorenstein. We prove that whether or not $U_P$ is ($\mathbb{Q}$)-Gorenstein depends only on the biconnected components of the Hasse diagram of $P$. In the case that $P$ has a minimum or maximum element, we show that the Gorenstein property of $U_P$ is completely determined by the M\"obius function of $P$. We also provide a recursive method that determines if $U_P$ is ($\mathbb{Q}$)-Gorenstein in this case. We conjecture that $U_P$ is Gorenstein if and only if it is $\mathbb{Q}$-Gorenstein. We verify this conjecture for posets of length $1$ and also for posets with a minimum or maximum element.


title: An Interface-Driven Adaptive Variational Procedure for Fully Eulerian Fluid-Structure Interaction via Phase-field Modeling

abstract: In this paper, we present a novel interface-driven adaptive variational procedure using a fully Eulerian description of fluid-structure interaction. The proposed fully-Eulerian procedure involves a fixed background unstructured mesh on which the fluid-structure interface is treated implicitly. We model the fluid-structure interaction by the phase-field finite element formulation relying on a partitioned staggered integration of the convective Allen-Cahn equation with the unified momentum equation for both solid and fluid dynamics. We employ the positivity preserving variational scheme for a bounded and stable solution of the Allen-Cahn phase-field equation. To evaluate the solid stresses, the left Cauchy-Green deformation tensor is convected at each time step to trace the evolution of the solid strain in the Eulerian reference frame. We utilize the residual based error indicators and the newest vertex bisection algorithm for the adaptive refinement/coarsening of the unstructured mesh. The proposed nonlinear adaptive partitioned procedure restricts the coarsening step to the last non-linear iteration while simultaneously ensuring convergence properties of the coupled governing equations. We perform a detailed convergence and accuracy analysis via two benchmark problems namely, the pure solid system and a coupled fluid-solid system with an interface in a rectangle domain. We next systematically assess the performance of the adaptive procedure in terms of conservation properties for the increasing complexity of problems. Finally, we demonstrate our fully-Eulerian interface-driven adaptive FSI model to simulate the contact and bouncing phenomenon between an elastic solid and a rigid wall.


title: Cameron's operator in terms of determinants, and hypergeometric numbers

abstract: By studying Cameron's operator in terms of determinants, two kinds of "integer" sequences of incomplete numbers were introduced. One was the sequence of restricted numbers, including $s$-step Fibonacci sequences. Another was the sequence of associated numbers, including Lam\'e sequences of higher order. By the classical Trudi's formula and the inverse relation, more expressions were able to be obtained. These relations and identities can be extended to those of sequence of negative integers or rational numbers. As applications, we consider hypergeometric Bernoulli, Cauchy and Euler numbers with some modifications.


title: Bayesian Optimization of Function Networks

abstract: We consider Bayesian optimization of the output of a network of functions, where each function takes as input the output of its parent nodes, and where the network takes significant time to evaluate. Such problems arise, for example, in reinforcement learning, engineering design, and manufacturing. While the standard Bayesian optimization approach observes only the final output, our approach delivers greater query efficiency by leveraging information that the former ignores: intermediate output within the network. This is achieved by modeling the nodes of the network using Gaussian processes and choosing the points to evaluate using, as our acquisition function, the expected improvement computed with respect to the implied posterior on the objective. Although the non-Gaussian nature of this posterior prevents computing our acquisition function in closed form, we show that it can be efficiently maximized via sample average approximation. In addition, we prove that our method is asymptotically consistent, meaning that it finds a globally optimal solution as the number of evaluations grows to infinity, thus generalizing previously known convergence results for the expected improvement. Notably, this holds even though our method might not evaluate the domain densely, instead leveraging problem structure to leave regions unexplored. Finally, we show that our approach dramatically outperforms standard Bayesian optimization methods in several synthetic and real-world problems.


title: Improving heavy Dirac neutrino prospects at future hadron colliders using machine learning

abstract: In this work, by using the machine learning methods, we study the sensitivities of heavy pseudo-Dirac neutrino $N$ in the inverse seesaw at the high-energy hadron colliders. The production process for the signal is $pp \to \ell N \to 3 \ell + E_T^{\rm miss}$, while the dominant background is $p p \to W Z \to 3 \ell + E_T^{\rm miss}$. We use either the Multi-Layer Perceptron or the Boosted Decision Tree with Gradient Boosting to analyse the kinematic observables and optimize the discrimination of background and signal events. It is found that the reconstructed $Z$ boson mass and heavy neutrino mass from the charged leptons and missing transverse energy play crucial roles in separating the signal from backgrounds. The prospects of heavy-light neutrino mixing $|V_{\ell N}|^2$ (with $\ell = e,\,\mu$) are estimated by using machine learning at the hadron colliders with $\sqrt{s}=14$ TeV, 27 TeV, and 100 TeV, and it is found that $|V_{\ell N}|^2$ can be improved up to ${\cal O} (10^{-6})$ for heavy neutrino mass $m_N = 100$ GeV and ${\cal O} (10^{-4})$ for $m_N = 1$ TeV.


title: On Orbifold Gromov-Witten Classes

abstract: We discuss some questions about Gromov-Witten classes of target stacks.


title: Quantum computational study of chloride ion attack on chloromethane for chemical accuracy and quantum noise effects with UCCSD and k-UpCCGSD ansatzes

abstract: Quantum computing is expected to play an important role in solving the problem of huge computational costs in various applications by utilizing the collective properties of quantum states, including superposition, interference, and entanglement, to perform computations. Quantum mechanical (QM) methods are candidates for various applications and can provide accurate absolute energy calculations in structure-based methods. QM methods are powerful tools for describing reaction pathways and their potential energy surfaces (PESs). In this study, we applied quantum computing to describe the PES of the bimolecular nucleophilic substitution (SN2) reaction between chloromethane and chloride ions. We performed noiseless and noise simulations using quantum algorithms and compared the accuracy and noise effects of the ansatzes. In noiseless simulations, the results from UCCSD and k-UpCCGSD are similar to those of full configurational interaction (FCI) with the same active space, which indicates that quantum algorithms can describe the PES of the SN2 reaction. In noise simulations, UCCSD is more susceptible to quantum noise than k-UpCCGSD. Therefore, k-UpCCGSD can serve as an alternative to UCCSD to reduce quantum noisy effects in the noisy intermediate-scale quantum era, and k-UpCCGSD is sufficient to describe the PES of the SN2 reaction in this work. The results showed the applicability of quantum computing to the SN2 reaction pathway and provided valuable information for structure-based molecular simulations with quantum computing.


title: Bayesian Testing Of Granger Causality In Functional Time Series

abstract: We develop a multivariate functional autoregressive model (MFAR), which captures the cross-correlation among multiple functional time series and thus improves forecast accuracy. We estimate the parameters under the Bayesian dynamic linear models (DLM) framework. In order to capture Granger causality from one FAR series to another we employ Bayes Factor. Motivated by the broad application of functional data in finance, we investigate the causality between the yield curves of two countries. Furthermore, we illustrate a climatology example, examining whether the weather conditions Granger cause pollutant daily levels in a city.


title: Homogeneous Formulation of Convex Quadratic Programs for Infeasibility Detection

abstract: Convex Quadratic Programs (QPs) have come to play a central role in the computation of control action for constrained dynamical systems. In this paper, we present a novel Homogeneous QP (HQP) formulation which is obtained by embedding the original QP in a larger space. The key properties of the HQP are: (i) is always feasible, (ii) an optimal solution to QP can be readily obtained from a solution to HQP, and (iii) infeasibility of QP corresponds to a particular solution of HQP. An immediate consequence is that all the existing algorithms for QP are now also capable of robustly detecting infeasibility. In particular, we present an Infeasible Interior Point Method (IIPM) for the HQP and show polynomial iteration complexity when applied to HQP. A key distinction with prior IPM approaches is that we do not need to solve second-order cone programs. Numerical experiments on the formulation are provided using existing codes.


title: SplitBrain: Hybrid Data and Model Parallel Deep Learning

abstract: The recent success of deep learning applications has coincided with those widely available powerful computational resources for training sophisticated machine learning models with huge datasets. Nonetheless, training large models such as convolutional neural networks using model parallelism (as opposed to data parallelism) is challenging because the complex nature of communication between model shards makes it difficult to partition the computation efficiently across multiple machines with an acceptable trade-off. This paper presents SplitBrain, a high performance distributed deep learning framework supporting hybrid data and model parallelism. Specifically, SplitBrain provides layer-specific partitioning that co-locates compute intensive convolutional layers while sharding memory demanding layers. A novel scalable group communication is proposed to further improve the training throughput with reduced communication overhead. The results show that SplitBrain can achieve nearly linear speedup while saving up to 67\% of memory consumption for data and model parallel VGG over CIFAR-10.


title: Higher-order social-ecological network as a simplicial complex

abstract: A social-ecological network is a formal representation of a corresponding social-ecological system, and encodes a relation within a given system as an interaction. Conventionally, such networks have been defined as encoding and representing pairwise interactions among the fundamental units of the system. This work proposes a combinatorial definition of social-ecological network by means of its structure as a simplicial complex. The proposed definition is a comprehensive one that takes into account the heterogeneity of interactions within a given SES, and the higher-order social-ecological network modelled using this definition is able to represent the modelled SES by capturing all orders of interactions within the system. Such a social-ecological network consequently, is better equipped to capture and represent the structural details of the real-world SES, and is thus capable of facilitating a deeper insight into the complex behaviour of the represented SES emergent through the higher-order interactions within the system, as compared to the conventional graph-theoretic network that exclusively models pairwise interactions.


title: A Critical Review of Inductive Logic Programming Techniques for Explainable AI

abstract: Despite recent advances in modern machine learning algorithms, the opaqueness of their underlying mechanisms continues to be an obstacle in adoption. To instill confidence and trust in artificial intelligence systems, Explainable Artificial Intelligence has emerged as a response to improving modern machine learning algorithms' explainability. Inductive Logic Programming (ILP), a subfield of symbolic artificial intelligence, plays a promising role in generating interpretable explanations because of its intuitive logic-driven framework. ILP effectively leverages abductive reasoning to generate explainable first-order clausal theories from examples and background knowledge. However, several challenges in developing methods inspired by ILP need to be addressed for their successful application in practice. For example, existing ILP systems often have a vast solution space, and the induced solutions are very sensitive to noises and disturbances. This survey paper summarizes the recent advances in ILP and a discussion of statistical relational learning and neural-symbolic algorithms, which offer synergistic views to ILP. Following a critical review of the recent advances, we delineate observed challenges and highlight potential avenues of further ILP-motivated research toward developing self-explanatory artificial intelligence systems.


title: InverseMV: Composing Piano Scores with a Convolutional Video-Music Transformer

abstract: Many social media users prefer consuming content in the form of videos rather than text. However, in order for content creators to produce videos with a high click-through rate, much editing is needed to match the footage to the music. This posts additional challenges for more amateur video makers. Therefore, we propose a novel attention-based model VMT (Video-Music Transformer) that automatically generates piano scores from video frames. Using music generated from models also prevent potential copyright infringements that often come with using existing music. To the best of our knowledge, there is no work besides the proposed VMT that aims to compose music for video. Additionally, there lacks a dataset with aligned video and symbolic music. We release a new dataset composed of over 7 hours of piano scores with fine alignment between pop music videos and MIDI files. We conduct experiments with human evaluation on VMT, SeqSeq model (our baseline), and the original piano version soundtrack. VMT achieves consistent improvements over the baseline on music smoothness and video relevance. In particular, with the relevance scores and our case study, our model has shown the capability of multimodality on frame-level actors' movement for music generation. Our VMT model, along with the new dataset, presents a promising research direction toward composing the matching soundtrack for videos. We have released our code at https://github.com/linchintung/VMT


title: Evolutionary correlation, regime switching, spectral dynamics and optimal trading strategies for cryptocurrencies and equities

abstract: This paper uses new and recently established methodologies to study the evolutionary dynamics of the cryptocurrency market, and compares the findings with that of the equity market. We begin by applying random matrix theory and principal components analysis (PCA) to correlation matrices of both collections, highlighting clear differences in the eigenspectra exhibited. We then explore the heterogeneity of both asset classes, studying the time-varying dynamics of underlying sector behaviours, and determine the collective similarity within each collection. We then turn to a study of structural break dynamics and evolutionary power spectra, where we quantify the collective affinity in structural breaks and evolutionary behaviours of underlying sector time series. Finally, we implement two algorithms simulating `portfolio choice' dynamics to compare the effectiveness of stock selection and sector allocation in cryptocurrency portfolios. There, we highlight the importance of both endeavours and comment on noteworthy implications for cryptocurrency portfolio management.


title: An Efficient and Robust Committee Structure for Sharding Blockchain

abstract: Nowadays, sharding is deemed as a promising way to save traditional blockchain protocols from their low scalability. However, such technique also brings several potential risks and huge communication overheads. An improper design may give rise to the inconsistent state among different committees. Further, the communication overheads arising from cross-shard transactions unfortunately reduce the system's performance. In this paper, we first summarize five essential issues that all sharding blockchain designers face. For each issue, we discuss its key challenge and propose our suggested solutions. In order to break the performance bottlenecks, we propose a reputation mechanism for selecting leaders. The term of reputation in our design reflects each node's honest computation resources. In addition, we introduce a referee committee and partial sets in each committee, and design a recovery procedure in case the leader is malicious. Under the design, we prove that malicious leaders will not hurt the system and will be evicted. Furthermore, we conduct a series of simulations to evaluate our design. The results show that selecting leaders by the reputation can dramatically improve the system performance.


title: Inhomogeneous disordering at a photo-induced charge density wave transition

abstract: Using ultrashort laser pulses, it has become possible to probe the dynamics of long-range order in solids on microscopic timescales. In the conventional description of symmetry-broken phases within time-dependent Ginzburg-Landau theory, the order parameter evolves coherently, with small fluctuations along an average trajectory. Recent experiments, however, indicate that some systems can support a different scenario, named ultrafast inhomogeneous disordering, where the average order parameter is no longer representative of the state on the atomic scale. Here we theoretically show that ultrafast disordering can occur in a minimal, yet paradigmatic, model for a Peierls instability if atomic scale inhomogeneities of both the electronic structure and the charge density wave order parameter are taken into account. The latter is achieved using a non-equilibrium generalization of statistical dynamical mean-field theory, coupled to stochastic differential equations for the order parameter.


title: Deconfounded Visual Grounding

abstract: We focus on the confounding bias between language and location in the visual grounding pipeline, where we find that the bias is the major visual reasoning bottleneck. For example, the grounding process is usually a trivial language-location association without visual reasoning, e.g., grounding any language query containing sheep to the nearly central regions, due to that most queries about sheep have ground-truth locations at the image center. First, we frame the visual grounding pipeline into a causal graph, which shows the causalities among image, query, target location and underlying confounder. Through the causal graph, we know how to break the grounding bottleneck: deconfounded visual grounding. Second, to tackle the challenge that the confounder is unobserved in general, we propose a confounder-agnostic approach called: Referring Expression Deconfounder (RED), to remove the confounding bias. Third, we implement RED as a simple language attention, which can be applied in any grounding method. On popular benchmarks, RED improves various state-of-the-art grounding methods by a significant margin. Code will soon be available at: https://github.com/JianqiangH/Deconfounded_VG.


title: Hamiltonian Monodromy via spectral Lax pairs

abstract: Hamiltonian Monodromy is the simplest topological obstruction to the existence of global action-angle coordinates in a completely integrable system. We show that this property can be studied in a neighborhood of a focus-focus singularity by a spectral Lax pair approach. From the Lax pair, we derive a Riemann surface which allows us to compute in a straightforward way the corresponding Monodromy matrix. The general results are applied to the Jaynes-Cummings model and the spherical pendulum.


title: An empirical Bayes approach to estimating dynamic models of co-regulated gene expression

abstract: Time-course gene expression datasets provide insight into the dynamics of complex biological processes, such as immune response and organ development. It is of interest to identify genes with similar temporal expression patterns because such genes are often biologically related. However, this task is challenging due to the high dimensionality of these datasets and the nonlinearity of gene expression time dynamics. We propose an empirical Bayes approach to estimating ordinary differential equation (ODE) models of gene expression, from which we derive a similarity metric between genes called the Bayesian lead-lag $R^2$ (LLR2). Importantly, the calculation of the LLR2 leverages biological databases that document known interactions amongst genes; this information is automatically used to define informative prior distributions on the ODE model's parameters. As a result, the LLR2 is a biologically-informed metric that can be used to identify clusters or networks of functionally-related genes with co-moving or time-delayed expression patterns. We then derive data-driven shrinkage parameters from Stein's unbiased risk estimate that optimally balance the ODE model's fit to both data and external biological information. Using real gene expression data, we demonstrate that our methodology allows us to recover interpretable gene clusters and sparse networks. These results reveal new insights about the dynamics of biological systems.