title: Consistency vs. Availability in Distributed Real-Time Systems

abstract: In distributed applications, Brewer's CAP theorem tells us that when networks become partitioned (P), one must give up either consistency (C) or availability (A). Consistency is agreement on the values of shared variables; availability is the ability to respond to reads and writes accessing those shared variables. Availability is a real-time property whereas consistency is a logical property. We have extended the CAP theorem to relate quantitative measures of these two properties to quantitative measures of communication and computation latency (L), obtaining a relation called the CAL theorem that is linear in a max-plus algebra. This paper shows how to use the CAL theorem in various ways to help design real-time systems. We develop a methodology for systematically trading off availability and consistency in application-specific ways and to guide the system designer when putting functionality in end devices, in edge computers, or in the cloud. We build on the Lingua Franca coordination language to provide system designers with concrete analysis and design tools to make the required tradeoffs in deployable software.


title: Corporate Culture and Organizational Fragility

abstract: Complex organizations accomplish tasks through many steps of collaboration among workers. Corporate culture supports collaborations by establishing norms and reducing misunderstandings. Because a strong corporate culture relies on costly, voluntary investments by many workers, we model it as an organizational public good, subject to standard free-riding problems, which become severe in large organizations. Our main finding is that voluntary contributions to culture can nevertheless be sustained, because an organization's equilibrium productivity is endogenously highly sensitive to individual contributions. However, the completion of complex tasks is then necessarily fragile to small shocks that damage the organization's culture.


title: On efficient quantum block encoding of pseudo-differential operators

abstract: Block encoding lies at the core of many existing quantum algorithms. Meanwhile, efficient and explicit block encodings of dense operators are commonly acknowledged as a challenging problem. This paper presents a comprehensive study of the block encoding of a rich family of dense operators: the pseudo-differential operators (PDOs). First, a block encoding scheme for generic PDOs is developed. Then we propose a more efficient scheme for PDOs with a separable structure. Finally, we demonstrate an explicit and efficient block encoding algorithm for PDOs with a dimension-wise fully separable structure. Complexity analysis is provided for all block encoding algorithms presented. The application of theoretical results is illustrated with worked examples, including the representation of variable coefficient elliptic operators and the computation of the inverse of elliptic operators without invoking quantum linear system algorithms (QLSAs).


title: Dal tramonto all'alba. Breve ricognizione sulla fisica italiana nel primo trentennio del XIX secolo

abstract: This paper presents an overview of the state of scientific research in physics in the Italian peninsula for the first thirty years of XIX century. In doing so, we focus geographically on the Kingdom of Sardinia and the Lombardo - Veneto territories, because of the role played in the development of Italian physics and the political future of the peninsula; we line out a quantitative analysis of the scientific community using basic tool of graph theory; we somehow try to avoid the major personalities of the century (e.g. Volta, Avogadro, Lagrange) to highlight the development of the local communities. The language of the paper is Italian.


title: Capacity-CRB Tradeoff in OFDM Integrated Sensing and Communication Systems

abstract: Integrated sensing and communication (ISAC) has emerged as a key technology for future communication systems. In this paper, we provide a general framework to reveal the fundamental tradeoff between sensing and communication in OFDM systems, where a unified ISAC waveform is exploited to perform both tasks. In particular, we define the Capacity-Bayesian Cramer Rao Bound (BCRB) region in the asymptotically case when the number of subcarriers is large. Specifically, we show that the asymptotically optimal input distribution that achieves the Pareto boundary point of the Capacity-BCRB region is Gaussian and the entire Pareto boundary can be obtained by solving a convex power allocation problem. Moreover, we characterize the structure of the sensing-optimal power allocation in the asymptotically case. Finally, numerical simulations are conducted to verify the theoretical analysis and provide useful insights.


title: An Optimized, Easy-to-use, Open-source GPU Solver for Large-scale Inverse Homogenization Problems

abstract: We propose a high-performance GPU solver for inverse homogenization problems to design high-resolution 3D microstructures. Central to our solver is a favorable combination of data structures and algorithms, making full use of the parallel computation power of today's GPUs through a software-level design space exploration. This solver is demonstrated to optimize homogenized stiffness tensors, such as bulk modulus, shear modulus, and Poisson's ratio, under the constraint of bounded material volume. Practical high-resolution examples with 512^3=134.2 million finite elements run in less than 40 seconds per iteration with a peak GPU memory of 9 GB on an NVIDIA GeForce GTX 1080Ti GPU. Besides, our GPU implementation is equipped with an easy-to-use framework with less than 20 lines of code to support various objective functions defined by the homogenized stiffness tensors. Our open-source high-performance implementation is publicly accessible at https://github.com/lavenklau/homo3d.


title: Rationalization for Explainable NLP: A Survey

abstract: Recent advances in deep learning have improved the performance of many Natural Language Processing (NLP) tasks such as translation, question-answering, and text classification. However, this improvement comes at the expense of model explainability. Black-box models make it difficult to understand the internals of a system and the process it takes to arrive at an output. Numerical (LIME, Shapley) and visualization (saliency heatmap) explainability techniques are helpful; however, they are insufficient because they require specialized knowledge. These factors led rationalization to emerge as a more accessible explainable technique in NLP. Rationalization justifies a model's output by providing a natural language explanation (rationale). Recent improvements in natural language generation have made rationalization an attractive technique because it is intuitive, human-comprehensible, and accessible to non-technical users. Since rationalization is a relatively new field, it is disorganized. As the first survey, rationalization literature in NLP from 2007-2022 is analyzed. This survey presents available methods, explainable evaluations, code, and datasets used across various NLP tasks that use rationalization. Further, a new subfield in Explainable AI (XAI), namely, Rational AI (RAI), is introduced to advance the current state of rationalization. A discussion on observed insights, challenges, and future directions is provided to point to promising research opportunities.


title: Unifying Structure Reasoning and Language Model Pre-training for Complex Reasoning

abstract: Recent pre-trained language models (PLMs) equipped with foundation reasoning skills have shown remarkable performance on downstream complex tasks. However, the significant structure reasoning skill has been rarely studied, which involves modeling implicit structure information within the text and performing explicit logical reasoning over them to deduce the conclusion. This paper proposes a unified learning framework that combines explicit structure reasoning and language pre-training to endow PLMs with the structure reasoning skill. It first identifies several elementary structures within contexts to construct structured queries and performs step-by-step reasoning along the queries to identify the answer entity. The fusion of textual semantics and structure reasoning is achieved by using contextual representations learned by PLMs to initialize the representation space of structures, and performing stepwise reasoning on this semantic representation space. Experimental results on four datasets demonstrate that the proposed model achieves significant improvements in complex reasoning tasks involving diverse structures, and shows transferability to downstream tasks with limited training data and effectiveness for complex reasoning of KGs modality.


title: ExClaim: Explainable Neural Claim Verification Using Rationalization

abstract: With the advent of deep learning, text generation language models have improved dramatically, with text at a similar level as human-written text. This can lead to rampant misinformation because content can now be created cheaply and distributed quickly. Automated claim verification methods exist to validate claims, but they lack foundational data and often use mainstream news as evidence sources that are strongly biased towards a specific agenda. Current claim verification methods use deep neural network models and complex algorithms for a high classification accuracy but it is at the expense of model explainability. The models are black-boxes and their decision-making process and the steps it took to arrive at a final prediction are obfuscated from the user. We introduce a novel claim verification approach, namely: ExClaim, that attempts to provide an explainable claim verification system with foundational evidence. Inspired by the legal system, ExClaim leverages rationalization to provide a verdict for the claim and justifies the verdict through a natural language explanation (rationale) to describe the model's decision-making process. ExClaim treats the verdict classification task as a question-answer problem and achieves a performance of 0.93 F1 score. It provides subtasks explanations to also justify the intermediate outcomes. Statistical and Explainable AI (XAI) evaluations are conducted to ensure valid and trustworthy outcomes. Ensuring claim verification systems are assured, rational, and explainable is an essential step toward improving Human-AI trust and the accessibility of black-box systems.


title: Improving Deep Regression with Ordinal Entropy

abstract: In computer vision, it is often observed that formulating regression problems as a classification task often yields better performance. We investigate this curious phenomenon and provide a derivation to show that classification, with the cross-entropy loss, outperforms regression with a mean squared error loss in its ability to learn high-entropy feature representations. Based on the analysis, we propose an ordinal entropy loss to encourage higher-entropy feature spaces while maintaining ordinal relationships to improve the performance of regression tasks. Experiments on synthetic and real-world regression tasks demonstrate the importance and benefits of increasing entropy for regression.


title: On the status of the star Schulte 12 in the association Cyg OB2

abstract: High-resolution spectra of the LBV candidate Schulte 12 in the Cyg OB2 association were obtained at the 6-meter BTA telescope with the NES echelle spectrograph on the arbitrary dates in 2001-2022. Variability of the emission profile of H$\alpha$ and absorptions of HeI, SiII with time was found. Based on the radial velocity measurements at 10 observation dates, radial velocity variability with an amplitude of $\Delta$Vr$\approx$8 km/s relative to the average value of the heliocentric velocity Vr=$-15.6\pm2.6$ km/s was revealed. This indicates the presence of a companion in the system. Based on the reliable intensity measurements of a sample of DIBs, color excess E(B-V)=1.74$\pm0.03^{m}$ was determined. This results in the interstellar extinction value Av$\approx 5.6^m$ that is only about half of the total extinction. Taking current Schulte 12 parameters, including Gaia~EDR3 parallax, we estimated its absolute magnitude as Mv$\approx -9.2^m$ and luminosity log$(L/L_{\odot})\approx$5.5, which does not exceed the Humphreys-Davidson limit.


title: A fast and flexible machine learning approach to data quality monitoring

abstract: We present a machine learning based approach for real-time monitoring of particle detectors. The proposed strategy evaluates the compatibility between incoming batches of experimental data and a reference sample representing the data behavior in normal conditions by implementing a likelihood-ratio hypothesis test. The core model is powered by recent large-scale implementations of kernel methods, nonparametric learning algorithms that can approximate any continuous function given enough data. The resulting algorithm is fast, efficient and agnostic about the type of potential anomaly in the data. We show the performance of the model on multivariate data from a drift tube chambers muon detector.


title: Improving Signed Propagation of Graph Neural Network Under Multiple Classes

abstract: Message-passing Graph Neural Networks (GNNs), which collect information from adjacent nodes achieve dismal performance on heterophilic graphs. Various schemes have been proposed to solve this problem, and propagating signed information on heterophilic edges has gained great attention. Recently, some works provided theoretical analysis that signed propagation always leads to performance improvement under a binary class scenario. However, we notice that prior analyses do not align well with multi-class benchmark datasets. This paper provides a new understanding of signed propagation for multi-class scenarios and points out two drawbacks in terms of message-passing and parameter update: (1) Message-passing: if two nodes belong to different classes but have a high similarity, signed propagation can decrease the separability. (2) Parameter update: the prediction uncertainty (e.g., conflict evidence) of signed neighbors increases during training, which can impede the stability of the algorithm. Based on the observation, we introduce two novel strategies for improving signed propagation under multi-class graphs. The proposed scheme combines calibration to secure robustness while reducing uncertainty. We show the efficacy of our theorem through extensive experiments on six benchmark graph datasets.


title: Ultra-stable and versatile high-energy resolution setup for attosecond photoelectron spectroscopy

abstract: Attosecond photoelectron spectroscopy is often performed with interferometric experimental setups that require outstanding stability. We demonstrate and characterize in detail an actively stabilized, versatile, high spectral resolution attosecond beamline. The active-stabilization system can remain ultra-stable for several hours with an RMS stability of 13 as and a total pump-probe delay scanning range of \sim 400 fs. A tunable femtosecond laser source to drive high-order harmonic generation allows for precisely addressing atomic and molecular resonances. Furthermore, the interferometer includes a spectral shaper in 4f-geometry in the probe arm as well as a tunable bandpass filter in the pump arm, which offer additional high flexibility in terms of tunability as well as narrowband or polychromatic probe pulses. We show that spectral phase measurements of photoelectron wavepackets with the rainbow RABBIT technique (reconstruction of attosecond beating by two photon transitions) with narrowband probe pulses can significantly improve the photoelectron energy resolution. In this setup, the temporal-spectral resolution of photoelectron spectroscopy can reach a new level of accuracy and precision.


title: Submodular Hypergraph Partitioning: Metric Relaxations and Fast Algorithms via an Improved Cut-Matching Game

abstract: Despite there being significant work on developing spectral, and metric embedding based approximation algorithms for hypergraph generalizations of conductance, little is known regarding the approximability of hypergraph partitioning objectives beyond this. This work proposes algorithms for a general model of hypergraph partitioning that unifies both undirected and directed versions of many well-studied partitioning objectives. The first contribution of this paper introduces polymatroidal cut functions, a large class of cut functions amenable to approximation algorithms via metric embeddings and routing multicommodity flows. We demonstrate an $O(\sqrt{\log n})$-approximation, where $n$ is the number of vertices in the hypergraph, for these problems by rounding relaxations to metrics of negative-type. The second contribution of this paper generalizes the cut-matching game framework of Khandekar et. al. to tackle polymatroidal cut functions. This yields the first almost-linear time $O(\log n)$-approximation algorithm for standard versions of undirected and directed hypergraph partitioning. A technical consequence of our construction is that a cut-matching game which greatly relaxes the set of allowed actions for both players can be used to partition hypergraphs with negligible impact on the approximation ratio. We believe this to be of independent interest.


title: Topological semiinfinite tensor (super)modules

abstract: We construct universal monoidal categories of topological tensor supermodules over the Lie superalgebras $\mathfrak{gl}(V\oplus \Pi V)$ and $\mathfrak{osp}(V\oplus \Pi V)$ associated with a Tate space $V$. Here $V\oplus \Pi V$ is a $\mathbb{Z}/2\mathbb{Z}$-graded topological vector space whose even and odd parts are isomorphic to $V$. We discuss the purely even case first, by introducing monoidal categories$\widehat{\mathbf{T}}_{\mathfrak{gl}(V)}$, $\widehat{\mathbf{T}}_{\mathfrak{o}(V)}$ and $\widehat{\mathbf{T}}_{\mathfrak{sp}(V)}$, and show that these categories are anti-equivalent to respective previously studied categories $\mathbb{T}_{\mathfrak{gl}(V)}$, $\mathbb{T}_{\mathfrak{o}(V)}$, $\mathbb{T}_{\mathfrak{sp}(V)}$. These latter categories have certain universality properties as monoidal categories, which consequently carry over to $\widehat{\mathbf{T}}_{\mathfrak{gl}(V)}$, $\widehat{\mathbf{T}}_{\mathfrak{o}(V)}$ and $\widehat{\mathbf{T}}_{\mathfrak{sp}(V)}$. Moreover, the categories $\mathbb{T}_{\mathfrak{o}(V)}$ and $\mathbb{T}_{\mathfrak{sp}(V)}$ are known to be equivalent, and this implies the equivalence of the categories $\widehat{\mathbf{T}}_{\mathfrak{o}(V)}$ and $\widehat{\mathbf{T}}_{\mathfrak{sp}(V)}$. After introducing a supersymmetric setting, we establish the equivalence of the category $\widehat{\mathbf{T}}_{\mathfrak{gl}(V)}$ with the category $\widehat{\mathbf{T}}_{\mathfrak{gl}(V\oplus \Pi V)}$, and the equivalence of both categories $\widehat{\mathbf{T}}_{\mathfrak{o}(V)}$ and $\widehat{\mathbf{T}}_{\mathfrak{sp}(V)}$ with $\widehat{\mathbf{T}}_{\mathfrak{osp}(V\oplus \Pi V)}$.


title: Strongest coronal magnetic fields in solar cycles 23-24: probing, statistics, and implications

abstract: Strong coronal magnetic field, when present, manifests itself as bright microwave sources at high frequencies produced by gyroresonant (GR) emission mechanism in thermal coronal plasma. The highest frequency at which this emission is observed is proportional to the absolute value of the strongest coronal magnetic field on the line of sight. Although no coronal magnetic field larger than roughly 2,000 G was expected, recently the field at least twice larger has been reported. Here, we report a search for and statistical study of such strong coronal magnetic fields using high-frequency GR emission. A historic record of spatially resolved microwave observations at high frequencies, 17 and 34 GHz, is available from Nobeyama RadioHeliograph for more than 20 years (1995-2018). Here we employ this data set to identify sources of bright GR emission at 34 GHz and perform a statistical analysis of the identified GR cases to quantify the strongest coronal magnetic fields during two solar cycles. We found that although active regions with the strong magnetic field are relatively rare (less than 1% of all active regions), they appear regularly on the Sun. These active regions are associated with prominent manifestations of solar activity.


title: The Risk-Taking Software Engineer: A Framed Portrait

abstract: Background: Risk-taking is prevalent in a host of activities performed by software engineers on a daily basis, yet there is scant research on it. Aims and Method: We study if software engineers' risk-taking is affected by framing effects and by software engineers' personality. To this end, we perform a survey experiment with 124 software engineers. Results: We find that framing substantially affects their risk-taking. None of the "Big Five" personality traits are related to risk-taking in software engineers after correcting for multiple testing. Conclusions: Software engineers and their managers must be aware of framing effects and account for them properly.


title: On Abelian Groups Having All Proper Characteristic Subgroups Isomorphic

abstract: We consider two variants of those Abelian groups with all proper characteristic subgroups isomorphic and give an in-depth study of their basic and specific properties in either parallel or contrast to the Abelian groups with all proper fully invariant subgroups isomorphic, which are studied in details by the current authors in Commun. Algebra (2015). In addition, we also examine those Abelian groups having at least one proper characteristic subgroup isomorphic to the whole group. The established by us results somewhat extend those obtained by Grinshpon-Nikolskaya in Tomsk State Univ. J. Math. & Mech. (2011, 2012) and in Commun. Algebra (2011), respectively.


title: New Challenges for Content Privacy in Speech and Audio

abstract: Privacy in speech and audio has many facets. A particularly under-developed area of privacy in this domain involves consideration for information related to content and context. Speech content can include words and their meaning or even stylistic markers, pathological speech, intonation patterns, or emotion. More generally, audio captured in-the-wild may contain background speech or reveal contextual information such as markers of location, room characteristics, paralinguistic sounds, or other audible events. Audio recording devices and speech technologies are becoming increasingly commonplace in everyday life. At the same time, commercialised speech and audio technologies do not provide consumers with a range of privacy choices. Even where privacy is regulated or protected by law, technical solutions to privacy assurance and enforcement fall short. This position paper introduces three important and timely research challenges for content privacy in speech and audio. We highlight current gaps and opportunities, and identify focus areas, that could have significant implications for developing ethical and safer speech technologies.


title: Sublattice extraordinary-log phase and new special point of the antiferromagnetic Potts model

abstract: We study the surface criticality of a three-dimensional classical antiferromagnetic Potts model, whose bulk critical behaviors belongs to the XY model because of emergent O(2) symmetry. We find that the surface antiferromagnetic next-nearest neighboring interactions can drive the extraordinary-log phase to the ordinary phase, the transition between the two phases belongs to the universality class of the well-known special transition of the XY model. Further strengthening the surface next-nearest neighboring interactions, the extraordinary-log phase reappears, but the main critical behaviors are dominated on the sublattices of the model; the special point between the ordinary phase and the sublattice extraordinary-log phase belongs to a new universality class.


title: Discovery of a magnetic Dirac system with large intrinsic non-linear Hall effect

abstract: Magnetic materials exhibiting topological Dirac fermions are attracting significant attention for their promising technological potential in spintronics. In these systems, the combined effect of the spin-orbit coupling and magnetic order enables the realization of novel topological phases with exotic transport properties, including the anomalous Hall effect and magneto-chiral phenomena. Herein, we report experimental signature of topological Dirac antiferromagnetism in TaCoTe2 via angle-resolved photoelectron spectroscopy (ARPES) and first-principles density functional theory (DFT) calculations. In particular, we find the existence of spin-orbit coupling-induced gaps at the Fermi level, consistent with the manifestation of a large intrinsic non-linear Hall conductivity. Remarkably, we find that the latter is extremely sensitive to the orientation of the N\'eel vector, suggesting TaCoTe2 a suitable candidate for the realization of non-volatile spintronic devices with an unprecedented level of intrinsic tunability.


title: Non-isothermal multicomponent flows with mass diffusion and heat conduction

abstract: A type-I model of non-isothermal multicomponent systems of gases describing mass diffusive and heat conductive phenomena is presented. The derivation of the model and a convergence result among thermomechanical theories in the smooth regime are discussed. Furthermore, the global-in-time existence of weak solutions and the weak-strong uniqueness property are established for the corresponding system with zero barycentric velocity.


title: A quantum model of charge capture and release onto/from deep traps

abstract: The rapid development of optical technologies and applications revealed the critical role of point defects affecting device performance. One of the powerful tools to study influence of defects on charge capture and recombination processes is thermoluminescence. The popular models behind thermoluminescence and carrier capture processes are semi-classic though. They offer good qualitative description, but implicitly exclude quantum nature of the accompanying parameters, such as frequency factors and capture cross sections. As a consequence, results obtained for a specific host material cannot be successfully extrapolated to other materials. Thus, the main purpose of our work is to introduce a reliable analytical model that describes non-radiative capture and release of electrons from/to the conduction band (CB). The proposed model is governed by Bose-Einstein statistics (for phonon occupation) and Fermi's golden rule (for resonant charge transfer between the trap and the CB). The constructed model offers a physical interpretation of the capture coefficients and frequency factors, and seamlessly includes the Coulomb neutral/attractive nature of traps. It connects the frequency factor to the overlap of wavefunctions of the delocalized CB and trap states, and suggests a strong dependence on the density of charge distribution, i.e. the ionicity/covalency of the chemical bonds within the host. Separation of the resonance condition from the accumulation/dissipation of phonons on the site leads to the conclusion that the capture cross-section does not necessarily depend on the trap depth. The model is verified by comparison to reported experimental data, showing good agreement. As such, the model generates reliable information about trap states whose exact nature is not completely understood and allows to do materials research in a more systematic way.


title: Dense RGB SLAM with Neural Implicit Maps

abstract: There is an emerging trend of using neural implicit functions for map representation in Simultaneous Localization and Mapping (SLAM). Some pioneer works have achieved encouraging results on RGB-D SLAM. In this paper, we present a dense RGB SLAM method with neural implicit map representation. To reach this challenging goal without depth input, we introduce a hierarchical feature volume to facilitate the implicit map decoder. This design effectively fuses shape cues across different scales to facilitate map reconstruction. Our method simultaneously solves the camera motion and the neural implicit map by matching the rendered and input video frames. To facilitate optimization, we further propose a photometric warping loss in the spirit of multi-view stereo to better constrain the camera pose and scene geometry. We evaluate our method on commonly used benchmarks and compare it with modern RGB and RGB-D SLAM systems. Our method achieves favorable results than previous methods and even surpasses some recent RGB-D SLAM methods.The code is at poptree.github.io/DIM-SLAM/.


title: Exponential sum approximations of finite completely monotonic functions

abstract: Bernstein's theorem (also called Hausdorff--Bernstein--Widder theorem) enables the integral representation of a completely monotonic function. We introduce a finite completely monotonic function, which is a completely monotonic function with a finite positive integral interval of the integral representation. We consider the exponential sum approximation of a finite completely monotonic function based on the Gaussian quadrature with a variable transformation. If the variable transformation is analytic on an open Bernstein ellipse, the maximum absolute error decreases at least geometrically with respect to the number of exponential functions. The maximization of the decreasing rate of the error bound can be achieved by using a variable transformation represented by Jacobi's delta amplitude function (also called dn function). The error curve is expanded by introducing basis functions, which are eigenfunctions of a fourth order differential operator, satisfy orthogonality conditions, and have the interlacing property of zeros by Kellogg's theorem.


title: On constructing benchmark quantum circuits with known near-optimal transformation cost

abstract: Current quantum devices impose strict connectivity constraints on quantum circuits, making circuit transformation necessary before running logical circuits on real quantum devices. Many quantum circuit transformation (QCT) algorithms have been proposed in the past several years. This paper proposes a novel method for constructing benchmark circuits and uses these benchmark circuits to evaluate state-of-the-art QCT algorithms, including TKET from Cambridge Quantum Computing, Qiskit from IBM, and three academic algorithms SABRE, SAHS, and MCTS. These benchmarks have known near-optimal transformation costs and thus are called QUEKNO (for quantum examples with known near-optimality). Compared with QUEKO benchmarks designed by Tan and Cong (2021), which all have zero optimal transformation costs, QUEKNO benchmarks are more general and can provide a more faithful evaluation for QCT algorithms (like TKET) which use subgraph isomorphism to find the initial mapping. Our evaluation results show that SABRE can generate transformations with conspicuously low average costs on the 53-qubit IBM Q Rochester and Google's Sycamore in both gate size and depth objectives.


title: A combinatorial formula for LLT cumulants of melting lollipops in terms of spanning trees

abstract: We prove a combinatorial formula for LLT cumulants of melting lollipops as a positive combination of LLT polynomials indexed by spanning trees. The result gives an affirmative answer to a general positivity question for this class of unicellular LLT cumulants, and gives an independent proof of their Schur-positivity. In the special case of the complete graph, we also express the formula in terms of parking functions.


title: Data-driven reduced order modeling for parametric PDE eigenvalue problems using Gaussian process regression

abstract: In this article, we propose a data-driven reduced basis (RB) method for the approximation of parametric eigenvalue problems. The method is based on the offline and online paradigms. In the offline stage, we generate snapshots and construct the basis of the reduced space, using a POD approach. Gaussian process regressions (GPR) are used for approximating the eigenvalues and projection coefficients of the eigenvectors in the reduced space. All the GPR corresponding to the eigenvalues and projection coefficients are trained in the offline stage, using the data generated in the offline stage. The output corresponding to new parameters can be obtained in the online stage using the trained GPR. The proposed algorithm is used to solve affine and non-affine parameter-dependent eigenvalue problems. The numerical results demonstrate the robustness of the proposed non-intrusive method.


title: Phonon-mediated dark to bright plasmon conversion

abstract: The optical response of a matter excitation embedded in nanophotonic devices is commonly described by the Drude-Lorentz model. Here, we demonstrate that this widely used approach fails in the case where quantum-confined plasmons of a two-dimensional electron gas interact strongly with optical phonons. We propose a new quantum model which contains the semiclassical Drude-Lorentz one for simple electronic potentials, but predicts very different results in symmetry-broken potentials. We unveil a new mechanism for the oscillator strength transfer between bright phonon-polariton and dark plasmon modes, enabling thus new quantum degrees of freedom for designing the optical response of nanostructures.


title: Lead-free, luminescent perovskite nanocrystals obtained through ambient condition synthesis

abstract: Heterovalent substitution of toxic lead is an increasingly popular design strategy to obtain environmentally sustainable variants of the exciting material class of halide perovskites. Perovskite nanocrystals (NCs) obtained through solution-based methods exhibit exceedingly high optical quality. Unfortunately, most of these synthesis routes still require reaction under inert gas and at very high temperatures. Herein we present a novel synthesis routine for lead-free double perovskite NCs. We combine hot injection and ligand-assisted reprecipitation (LARP) methods to achieve a low-temperature and ambient atmosphere-based synthesis for manganese-doped Cs_{2}NaBiCl_{6} NCs. Mn incorporation is critical for the otherwise non-emissive material, with a 9:1 Bi:Mn precursor ratio maximizing the bright orange photoluminescence (PL) and quantum yield (QY). Higher temperatures slightly increased the material's performance, yet NCs synthesized at room temperature were still emissive, highlighting the versatility of the synthetic approach. Furthermore, the NCs show excellent long-term stability in ambient conditions, facilitating additional investigations and energy-related applications.


title: Exploring Methods for Building Dialects-Mandarin Code-Mixing Corpora: A Case Study in Taiwanese Hokkien

abstract: In natural language processing (NLP), code-mixing (CM) is a challenging task, especially when the mixed languages include dialects. In Southeast Asian countries such as Singapore, Indonesia, and Malaysia, Hokkien-Mandarin is the most widespread code-mixed language pair among Chinese immigrants, and it is also common in Taiwan. However, dialects such as Hokkien often have a scarcity of resources and the lack of an official writing system, limiting the development of dialect CM research. In this paper, we propose a method to construct a Hokkien-Mandarin CM dataset to mitigate the limitation, overcome the morphological issue under the Sino-Tibetan language family, and offer an efficient Hokkien word segmentation method through a linguistics-based toolkit. Furthermore, we use our proposed dataset and employ transfer learning to train the XLM (cross-lingual language model) for translation tasks. To fit the code-mixing scenario, we adapt XLM slightly. We found that by using linguistic knowledge, rules, and language tags, the model produces good results on CM data translation while maintaining monolingual translation quality.


title: Validation of non-equilibrium kinetics in CO2-N2 plasmas

abstract: This work explores the effect of N2 addition on CO2 dissociation and on the vibrational kinetics of CO2 and CO under various non-equilibrium plasma conditions. A self-consistent kinetic model, previously validated for pure CO2 and CO2-O2 discharges, is further extended by adding the kinetics of N2. The vibrational kinetics considered include levels up to v = 10 for CO, v = 59 for N2 and up to v1 = 2 and v2 = v3 = 5, respectively for the symmetric stretch, bending and asymmetric stretch modes of CO2, and account for electron-impact excitation and de-excitation (e-V), vibration-to-translation (V-T) and vibration-to-vibration energy exchange (V-V) processes. The kinetic scheme is validated by comparing the model predictions with recent experimental data measured in a DC glow discharge operating in pure CO2 and in CO2-N2 mixtures, at pressures in the range 0.6 - 4 Torr (80.00 - 533.33 Pa) and a current of 50 mA. The experimental results show a higher vibrational temperature of the different modes of CO2 and CO and an increased dissociation fraction of CO2, that can reach values as high as 70%, when N2 is added to the plasma. On the one hand, the simulations suggest that the former effect is the result of the CO2-N2 and CO-N2 V-V transfers and the reduction of quenching due to the decrease of atomic oxygen concentration; on the other hand, the dilution of CO2 and dissociation products, CO and O2, reduces the importance of back reactions and contributes to the higher CO2 dissociation fraction with increased N2 content in the mixture, while the N2(B3Pg) electronically excited state further enhances the CO2 dissociation.


title: Counterfactual Explanation and Instance-Generation using Cycle-Consistent Generative Adversarial Networks

abstract: The image-based diagnosis is now a vital aspect of modern automation assisted diagnosis. To enable models to produce pixel-level diagnosis, pixel-level ground-truth labels are essentially required. However, since it is often not straight forward to obtain the labels in many application domains such as in medical image, classification-based approaches have become the de facto standard to perform the diagnosis. Though they can identify class-salient regions, they may not be useful for diagnosis where capturing all of the evidences is important requirement. Alternatively, a counterfactual explanation (CX) aims at providing explanations using a casual reasoning process of form "If X has not happend, Y would not heppend". Existing CX approaches, however, use classifier to explain features that can change its predictions. Thus, they can only explain class-salient features, rather than entire object of interest. This hence motivates us to propose a novel CX strategy that is not reliant on image classification. This work is inspired from the recent developments in generative adversarial networks (GANs) based image-to-image domain translation, and leverages to translate an abnormal image to counterpart normal image (i.e. counterfactual instance CI) to find discrepancy maps between the two. Since it is generally not possible to obtain abnormal and normal image pairs, we leverage Cycle-Consistency principle (a.k.a CycleGAN) to perform the translation in unsupervised way. We formulate CX in terms of a discrepancy map that, when added from the abnormal image, will make it indistinguishable from the CI. We evaluate our method on three datasets including a synthetic, tuberculosis and BraTS dataset. All these experiments confirm the supremacy of propose method in generating accurate CX and CI.


title: Quasi-optimal Reinforcement Learning with Continuous Actions

abstract: Many real-world applications of reinforcement learning (RL) require making decisions in continuous action environments. In particular, determining the optimal dose level plays a vital role in developing medical treatment regimes. One challenge in adapting existing RL algorithms to medical applications, however, is that the popular infinite support stochastic policies, e.g., Gaussian policy, may assign riskily high dosages and harm patients seriously. Hence, it is important to induce a policy class whose support only contains near-optimal actions, and shrink the action-searching area for effectiveness and reliability. To achieve this, we develop a novel \emph{quasi-optimal learning algorithm}, which can be easily optimized in off-policy settings with guaranteed convergence under general function approximations. Theoretically, we analyze the consistency, sample complexity, adaptability, and convergence of the proposed algorithm. We evaluate our algorithm with comprehensive simulated experiments and a dose suggestion real application to Ohio Type 1 diabetes dataset.


title: On the Algebraic Properties of Flame Graphs

abstract: Flame graphs are a popular way of representing profiling data. In this paper we propose a possible mathematical definition of flame graphs. In doing so, we gain some interesting algebraic properties almost for free, which in turn allow us to define some operations that can allow to perform an in-depth performance regression analysis. The typical documented use of a flame graph is via its graphical representation, whereby one scans the picture for the largest plateaux. Whilst this method is effective at finding the main sources of performance issues, it leaves quite a large amount of data potentially unused. By combining a mathematical precise definition of flame graphs with some statistical methods we show how to generalise this visual procedure and make the best of the full set of collected profiling data.


title: A central limit theorem for a sequence of conditionally centered and $\alpha$-mixing random fields

abstract: A central limit theorem is established for a sum of random variables belonging to a sequence of random fields. The fields are assumed to have zero mean conditional on the past history and to satisfy certain conditional $\alpha$-mixing conditions in space or time. The limiting normal distribution is obtained for increasing spatial domain or increasing length of the sequence. The applicability of the theorem is demonstrated by examples regarding estimating functions for a space-time point process and a space-time Markov process.


title: Tuning electronic properties and contact type in van der Waals heterostructures of bilayer SnS and graphene

abstract: Using first-principles calculations, we study the structural and electronic properties of the bilayer SnS/graphene, bilayer SnS/bilayer graphene (AA-stacked), bilayer SnS/bilayer graphene (AB-stacked) and monolayer SnS/graphene/monolayer SnS van der Waals (vdW) heterostructures. Electronic properties of all components of the vdW heterostructures are well preserved, which reflects the weakness of the vdW interaction. In the cases of bilayer SnS/graphene and bilayer SnS/bilayer graphene (AA-stacked), an Ohmic contact is formed which can be turned first into p-type and then into n-type Schottky contacts via application of an external electric field. Calculations show that an Ohmic contact is also formed at the interface of bilayer SnS/bilayer graphene (AB-stacked) heterostructure, but interestingly, by applying the perpendicular electric field a transition from semimetal/semiconductor contact to semiconductor/semiconductor one occurs which can enhance its optical properties. Alternatively, in the monolayer SnS/graphene/monolayer SnS vdW heterosructure, a p-type Schottky contact is established that changes into Ohmic contact under an applied electric field. Our results clearly indicate that the electronic properties of the vdW heterostructures can be tuned efficiently by external electric field, which is important in designing of new nanoelectronic devices.


title: A high order unfitted finite element method for time-Harmonic Maxwell interface problems

abstract: In this paper, we propose a high order unfitted finite element method for solving time-harmonic Maxwell interface problems. The unfitted finite element method is based on a mixed formulation in the discontinuous Galerkin framework on a Cartesian mesh with possible hanging nodes. The $H^2$ regularity of the solution to Maxwell interface problems with $C^2$ interfaces in each subdomain is proved. Practical interface resolving mesh conditions are introduced under which the hp inverse estimates on three-dimensional curved domains are proved. Stability and hp a priori error estimate of the unfitted finite element method are proved. Numerical results are included to illustrate the performance of the method.


title: Molecular hydrogen and its proxies HCO$^+$ and CO in the diffuse interstellar medium

abstract: There is a robust polyatomic chemistry in diffuse, partially-molecular interstellar gas that is readily accessible in absorption at radio/mm/sub-mm wavelengths. Accurate column densities are derived owing to the weak internal excitation, so relative molecular abundances are well known with respect to each other but not with respect to H2. Here we consider the use of proxies for hydrogen column densities N(H2) and N(H) = N(HI)+2N(H2) based on measurements of HCO+ absorption and CO emission and absorption, and we compare these with results obtained by others when observing HI, H2 and CO toward stars and AGN. We consider the use of HCO+ as a proxy for H2 and show that the assumption of a relative abundance N(H2) = N(HCO+)/3x10^{-9} gives the same view of the atomic-molecular hydrogen transition that is seen in UV absorption toward stars. CO on the other hand shows differences between the radio and optical regimes because emission is always detected when N(\hcop) > 6x10^{11}\pcc or N(H2) > 2x10^20\pcc. Wide variations in the integrated CO {J=1-0} brightness W_CO and N(CO)/N(H2) imply equivalent variations in the CO-H2 conversion factor even while the ensemble mean is near the usual Galactic values. Gas/reddening ratios found in absorption toward stars, N(H)/E(B-V) = 6.2x10^21 H \pcc/mag overall or 6.8x10^21 H \pcc/mag for sightlines at E(B-V) <= 0.08 mag lacking H2 are well below the Galactic mean measured at low reddening and high Galactic latitude, 8.3x10^21 H \pcc/mag.


title: Heating or Cooling: Study of Advective Heat Transport in the Inflow and the Outflow of Optically Thin Advection-dominated Accretion Flows

abstract: Advection is believed to be the dominant cooling mechanism in optically thin advection-dominated accretion flows (ADAF's). When outflow is considered, however, the first impression is that advection should be of opposite sign in the inflow and the outflow, due to the opposite direction of radial motion. Then how is the energy balance achieved simultaneously? We investigate the problem in this paper, analysing the profiles of different components of advection with self-similar solutions of ADAF's in spherical coordinates ($r\theta\phi$). We find that for $n < 3\gamma/2-1$, where $n$ is the density index in $\rho \propto r^{-n}$ and $\gamma$ is the heat capacity ratio, the radial advection is a heating mechanism in the inflow and a cooling mechanism in the outflow. It becomes 0 for $n = 3\gamma/2-1$, and turns to a cooling mechanism in the inflow and a heating mechanism in the outflow for $n > 3\gamma/2-1$. The energy conservation is only achieved when the latitudinal ($\theta$-direction) advection is considered, which takes an appropriate value to maintain energy balance, so that the overall effect of advection, no matter the parameter choices, is always a cooling mechanism that cancels out the viscous heating everywhere. For the extreme case of $n=3/2$, latitudinal motion stops, viscous heating is balanced solely by radial advection, and no outflow is developed.


title: An encounter-based approach to the escape problem

abstract: We revise the encounter-based approach to imperfect diffusion-controlled reactions, which employs the statistics of encounters between a diffusing particle and the reactive region to implement surface reactions. We extend this approach to deal with a more general setting, in which the reactive region is surrounded by a reflecting boundary with an escape region. We derive a spectral expansion for the full propagator and investigate the behavior and probabilistic interpretations of the associated probability flux density. In particular, we obtain the joint probability density of the escape time and the number of encounters with the reactive region before escape, and the probability density of the first-crossing time of a prescribed number of encounters. We briefly discuss generalizations of the conventional Poissonian-type surface reaction mechanism described by Robin boundary condition and potential applications of this formalism in chemistry and biophysics.

categories: cond-mat.stat-mech math-ph math.MP physics.chem-ph


title: A Perron-Frobenius analysis of wall-bounded turbulence

abstract: The Perron-Frobenius operator (PFO) is adapted from dynamical-system theory to the study of turbulent channel flow. It is shown that, as long as the analysis is restricted to the system attractor, the PFO can be used to differentiate causality and coherence from simple correlation without performing interventional experiments, and that the key difficulty remains collecting enough data to populate the operator matrix. This is alleviated by limiting the analysis to two-dimensional projections of the phase space, and developing a series of indicators to choose the best parameter pairs from a large number of possibilities. The techniques thus developed are applied to the study of bursting in the inertial layer of the channel, with emphasis on the process by which bursts are reinitiated after they have decayed. Conditional averaging over phase-space trajectories suggested by the PFO shows, somewhat counter-intuitively, that a key ingredient for the burst recovery is the development of a low-shear region near the wall, overlaid by a lifted shear layer. This is confirmed by a computational experiment in which the control of the mean velocity profile by the turbulence fluctuations is artificially relaxed. The behaviour of the mean velocity profile is thus modified, but the association of low wall shear with the initiation of the bursts is maintained.


title: Estimation of Sea State Parameters from Ship Motion Responses Using Attention-based Neural Networks

abstract: On-site estimation of sea state parameters is crucial for ship navigation systems' accuracy, stability, and efficiency. Extensive research has been conducted on model-based estimating methods utilizing only ship motion responses. Model-free approaches based on machine learning (ML) have recently gained popularity, and estimation from time-series of ship motion responses using deep learning (DL) methods has given promising results. Accordingly, in this study, we apply the novel, attention-based neural network (AT-NN) for estimating sea state parameters (wave height, zero-crossing period, and relative wave direction) from raw time-series data of ship pitch, heave, and roll motions. Despite using reduced input data, it has been successfully demonstrated that the proposed approaches by modified state-of-the-art techniques (based on convolutional neural networks (CNN) for regression, multivariate long short-term memory CNN, and sliding puzzle neural network) reduced estimation MSE by 23% and MAE by 16% compared to the original methods. Furthermore, the proposed technique based on AT-NN outperformed all tested methods (original and enhanced), reducing estimation MSE by up to 94% and MAE by up to 70%. Finally, we also proposed a novel approach for interpreting the uncertainty estimation of neural network outputs based on the Monte-Carlo dropout method to enhance the model's trustworthiness.


title: Genetically Modified Wolf Optimization with Stochastic Gradient Descent for Optimising Deep Neural Networks

abstract: When training Convolutional Neural Networks (CNNs) there is a large emphasis on creating efficient optimization algorithms and highly accurate networks. The state-of-the-art method of optimizing the networks is done by using gradient descent algorithms, such as Stochastic Gradient Descent (SGD). However, there are some limitations presented when using gradient descent methods. The major drawback is the lack of exploration, and over-reliance on exploitation. Hence, this research aims to analyze an alternative approach to optimizing neural network (NN) weights, with the use of population-based metaheuristic algorithms. A hybrid between Grey Wolf Optimizer (GWO) and Genetic Algorithms (GA) is explored, in conjunction with SGD; producing a Genetically Modified Wolf optimization algorithm boosted with SGD (GMW-SGD). This algorithm allows for a combination between exploitation and exploration, whilst also tackling the issue of high-dimensionality, affecting the performance of standard metaheuristic algorithms. The proposed algorithm was trained and tested on CIFAR-10 where it performs comparably to the SGD algorithm, reaching high test accuracy, and significantly outperforms standard metaheuristic algorithms.


title: Time-Conditioned Generative Modeling of Object-Centric Representations for Video Decomposition and Prediction

abstract: When perceiving the world from multiple viewpoints, humans have the ability to reason about the complete objects in a compositional manner even when an object is completely occluded from certain viewpoints. Meanwhile, humans are able to imagine novel views after observing multiple viewpoints. Recent remarkable advances in multi-view object-centric learning still leaves some unresolved problems: 1) The shapes of partially or completely occluded objects can not be well reconstructed. 2) The novel viewpoint prediction depends on expensive viewpoint annotations rather than implicit rules in view representations. In this paper, we introduce a time-conditioned generative model for videos. To reconstruct the complete shape of an object accurately, we enhance the disentanglement between the latent representations of objects and views, where the latent representations of time-conditioned views are jointly inferred with a Transformer and then are input to a sequential extension of Slot Attention to learn object-centric representations. In addition, Gaussian processes are employed as priors of view latent variables for video generation and novel-view prediction without viewpoint annotations. Experiments on multiple datasets demonstrate that the proposed model can make object-centric video decomposition, reconstruct the complete shapes of occluded objects, and make novel-view predictions.


title: T-symmetry in String Geometry Theory

abstract: String geometry theory is one of the candidates of non-perturbative formulation of string theory. In this paper, we have shown that dimensionally reduced string geometry theories have what we call T-symmetry. In case of the dimensional reduction in space-like directions, the T-symmetry transformation gives the T-dual transformation between the type IIA and IIB perturbative vacua. In case of the dimensional reduction in the direction of string geometry time $\bar{\tau}$, the T-symmetry transformation is independent of the T-dual transformation, and gives a symmetry that cannot be seen in the perturbative string theories.


title: Second-Order Coverage Control for Multi-Agent UAV Photogrammetry

abstract: Unmanned Aerial Vehicles equipped with cameras can be used to automate image capture for generating 3D models via photogrammetry. Current methods rely on a single vehicle to capture images sequentially, or use pre-planned and heuristic imaging configurations. We seek to provide a multi-agent control approach to capturing the images required to 3D map a region. A photogrammetry cost function is formulated that captures the importance of sharing feature-dense areas across multiple images for successful photogrammetry reconstruction. A distributed second-order coverage controller is used to minimise this cost and move agents to an imaging configuration. This approach prioritises high quality images that are simultaneously captured, leading to efficient and scalable 3D mapping of a region. We demonstrate our approach with a hardware experiment, generating and comparing 3D reconstructions from image sets captured using our approach to those captured using traditional methods.


title: First post-Newtonian $N$-body problem in Einstein-Cartan theory with the Weyssenhoff fluid: Lagrangian and first integrals

abstract: The rotational dynamics of an $N$-body system at the first post-Newtonian order in Einstein-Cartan theory is derived. This result is achieved by performing the point-particle limit of the equations of motion of the Weyssenhoff fluid, which models the quantum spin effects residing inside the bodies. For the special case of binary systems, we determine the Lagrangian function and the resulting first integrals underlying the translational dynamics and the spin precession.


title: Determination of nonthermal bonding origin of a novel photoexcited lattice instability in SnSe

abstract: Interatomic forces that bind materials are largely determined by an often complex interplay between the electronic band-structure and the atomic arrangements to form its equilibrium structure and dynamics. As these forces also determine the phonon dispersion, lattice dynamics measurements are often crucial tools for understanding how materials transform between different structures. This is the case for the mono-chalcogenides which feature a number of lattice instabilities associated with their network of resonant bonds and a large tunability in their functional properties. SnSe hosts a novel lattice instability upon above-bandgap photoexcitation that is distinct from the distortions associated with its high temperature phase transition, demonstrating that photoexcitation can alter the interatomic forces significantly different than thermal excitation. Here we report decisive time-resolved X-ray scattering-based measurements of the nonequlibrium lattice dynamics in SnSe. By fitting interatomic force models to the excited-state dispersion, we determine this instability as being primarily due to changes in the fourth-nearest neighbor bonds that connect bilayers, with relatively little change to the intralayer resonant bonds. In addition to providing critical insight into the nonthermal bonding origin of the instability in SnSe, such measurements will be crucial for understanding and controlling materials properties under non-equilibrium conditions.


title: Using deterministic tourist walk as a small-world metric on Watts-Strogatz networks

abstract: The Watts-Strogatz model (WS) has been demonstrated to effectively describe real-world networks due to its ability to reproduce the small-world properties commonly observed in a variety of systems, including social networks, computer networks, biochemical reactions, and neural networks. As the presence of small-world properties is a prevalent characteristic in many real-world networks, the measurement of "small-worldness" has become a crucial metric in the field of network science, leading to the development of various methods for its assessment over the past two decades. In contrast, the deterministic tourist walk (DTW) method has emerged as a prominent technique for texture analysis and network classification. In this paper, we propose the use of a modified version of the DTW method to classify networks into three categories: regular networks, random networks, and small-world networks. Additionally, we construct a small-world metric, denoted by the coefficient $\chi$, from the DTW method. Results indicate that the proposed method demonstrates excellent performance in the task of network classification, achieving over $90\%$ accuracy. Furthermore, the results obtained using the coefficient $\chi$ on real-world networks provide evidence that the proposed method effectively serves as a satisfactory small-world metric.


title: Slice Transformer and Self-supervised Learning for 6DoF Localization in 3D Point Cloud Maps

abstract: Precise localization is critical for autonomous vehicles. We present a self-supervised learning method that employs Transformers for the first time for the task of outdoor localization using LiDAR data. We propose a pre-text task that reorganizes the slices of a $360^\circ$ LiDAR scan to leverage its axial properties. Our model, called Slice Transformer, employs multi-head attention while systematically processing the slices. To the best of our knowledge, this is the first instance of leveraging multi-head attention for outdoor point clouds. We additionally introduce the Perth-WA dataset, which provides a large-scale LiDAR map of Perth city in Western Australia, covering $\sim$4km$^2$ area. Localization annotations are provided for Perth-WA. The proposed localization method is thoroughly evaluated on Perth-WA and Appollo-SouthBay datasets. We also establish the efficacy of our self-supervised learning approach for the common downstream task of object classification using ModelNet40 and ScanNN datasets. The code and Perth-WA data will be publicly released.


title: A Practical Introduction to Regression Discontinuity Designs: Extensions

abstract: This monograph, together with its accompanying first part Cattaneo, Idrobo and Titiunik (2020), collects and expands the instructional materials we prepared for more than $50$ short courses and workshops on Regression Discontinuity (RD) methodology that we taught between 2014 and 2023. In this second monograph, we discuss several topics in RD methodology that build on and extend the analysis of RD designs introduced in Cattaneo, Idrobo and Titiunik (2020). Our first goal is to present an alternative RD conceptual framework based on local randomization ideas. This methodological approach can be useful in RD designs with discretely-valued scores, and can also be used more broadly as a complement to the continuity-based approach in other settings. Then, employing both continuity-based and local randomization approaches, we extend the canonical Sharp RD design in multiple directions: fuzzy RD designs, RD designs with discrete scores, and multi-dimensional RD designs. The goal of our two-part monograph is purposely practical and hence we focus on the empirical analysis of RD designs.


title: Successive Subspace Learning for Cardiac Disease Classification with Two-phase Deformation Fields from Cine MRI

abstract: Cardiac cine magnetic resonance imaging (MRI) has been used to characterize cardiovascular diseases (CVD), often providing a noninvasive phenotyping tool.~While recently flourished deep learning based approaches using cine MRI yield accurate characterization results, the performance is often degraded by small training samples. In addition, many deep learning models are deemed a ``black box," for which models remain largely elusive in how models yield a prediction and how reliable they are. To alleviate this, this work proposes a lightweight successive subspace learning (SSL) framework for CVD classification, based on an interpretable feedforward design, in conjunction with a cardiac atlas. Specifically, our hierarchical SSL model is based on (i) neighborhood voxel expansion, (ii) unsupervised subspace approximation, (iii) supervised regression, and (iv) multi-level feature integration. In addition, using two-phase 3D deformation fields, including end-diastolic and end-systolic phases, derived between the atlas and individual subjects as input offers objective means of assessing CVD, even with small training samples. We evaluate our framework on the ACDC2017 database, comprising one healthy group and four disease groups. Compared with 3D CNN-based approaches, our framework achieves superior classification performance with 140$\times$ fewer parameters, which supports its potential value in clinical use.


title: New symmetries for Dyson's rank function

abstract: At the 1987 Ramanujan Centenary meeting Dyson asked for a coherent group-theoretical structure for Ramanujan's mock theta functions analogous to Hecke's theory of modular forms. Many of Ramanujan's mock theta functions can be written in terms of $R(\zeta_p,q)$, where $R(z,q)$ is the two-variable generating function of Dyson's rank function and $\zeta_p$ is a primitive $p$-th root of unity. In his lost notebook Ramanujan gives the $5$-dissection of $R(\zeta_5,q)$. This result is related to Dyson's famous rank conjecture which was proved by Atkin and Swinnerton-Dyer. In 2016 the first author showed that there is an analogous result for the $p$-dissection of $R(\zeta_p,q)$ when $p$ is any prime greater than $3$, by extending work of Bringmann and Ono, and Ahlgren and Treneer. It was also shown how the group $\Gamma_1(p)$ acts on the elements of the $p$-dissection of $R(\zeta_p,q)$. We extend this to the group $\Gamma_0(p)$, thus revealing new and surprising symmetries for Dyson's rank function.


title: On x-independence of RQ = FLQ/F2Q ratio at low x

abstract: We give predictions for the ratio RQ(x,Q2)=FLQ(x,Q2)/F2Q(x,Q2) at small values of Bjorken variable x in the first three orders of perturbation theory. We show an approximate x-independence of RQ(x,Q2) at low x and non-large Q2 values (Q2 < 8--10 m2Q), irrespectively on the gluon density in a proton used in the calculations. This observation could be useful in subsequent phenomenological studies of the heavy flavor production at future lepton-hadron and hadron-hadron colliders.


title: Leveraging Spatial and Temporal Correlations for Network Traffic Compression

abstract: The deployment of modern network applications is increasing the network size and traffic volumes at an unprecedented pace. Storing network-related information (e.g., traffic traces) is key to enable efficient network management. However, this task is becoming more challenging due to the ever-increasing data transmission rates and traffic volumes. In this paper, we present a novel method for network traffic compression that exploits spatial and temporal patterns naturally present in network traffic. We consider a realistic scenario where traffic measurements are performed at multiple links of a network topology using tools like SNMP or NetFlow. Such measurements can be seen as multiple time series that exhibit spatial and temporal correlations induced by the network topology, routing or user behavior. Our method leverages graph learning methods to effectively exploit both types of correlations for traffic compression. The experimental results show that our solution is able to outperform GZIP, the \textit{de facto} traffic compression method, improving by 50\%-65\% the compression ratio on three real-world networks.


title: Study of Dynamic Interaction Between Low Re Aerodynamic Load and Flexible-Biomimetic Wings with Tailorable Stiffness by FSI Modeling

abstract: In the present work, we investigate dynamic interaction and response of flexible bio-inspired morphing wing structure to a low Reynolds aerodynamic load. The aspects of inspiration are as follows. First, the segmentation of the wing into rigid and flexible segments. Considering a leading edge constitution of bone and muscle. In addition to a flexible trailing edge composed of feathers. Second, the material properties provided by experimental biology in literature are adopted such as the bending stiffness and Young's modulus. The development of numerical models allowing non uniform distribution of properties are developed and implemented into an OpenFoam finite volume solver that couples fluid dynamics to a structural solid dynamics solver through the FSI interface. In the course of this work, the validation is performed for a NACA6409 airfoil considering a rigid segment of 40\% and flexible segment 60\% chord length in order to test the aero-structure behavior for an aerodynamic load of air flow at low Reynolds number of $\boldsymbol{5\times 10^5}$ for the fluid and feather inspired material properties. The results suggest that bio-inspired techniques can be reproduced in engineering configurations.


title: Boundedness of the dyadic maximal function on graded Lie groups

abstract: Let $1<p\leq \infty$ and let $n\geq 2.$ It was proved independently by C. Calder\'on, R. Coifman and G. Weiss that the dyadic maximal function \begin{equation*} \mathcal{M}^{d\sigma}_Df(x)=\sup_{j\in\mathbb{Z}}\left|\smallint\limits_{\mathbb{S}^{n-1}}f(x-2^jy)d\sigma(y)\right| \end{equation*} is a bounded operator on $L^p(\mathbb{R}^n)$ where $d\sigma(y)$ is the surface measure on $\mathbb{S}^{n-1}.$ In this paper we prove an analogue of this result on arbitrary graded Lie groups. More precisely, to any finite Borel measure $d\sigma$ with compact support on a graded Lie group $G,$ we associate the corresponding dyadic maximal function $\mathcal{M}_D^{d\sigma}$ using the homogeneous structure of the group. Then, we prove a criterion in terms of the order (at zero and at infinity) of the group Fourier transform $\widehat{d\sigma}$ of $d\sigma$ with respect to a fixed Rockland operator $\mathcal{R}$ on $G$ that assures the boundedness of $\mathcal{M}_D^{d\sigma}$ on $L^p(G)$ for all $1<p\leq \infty.$


title: Raw or Cooked? Object Detection on RAW Images

abstract: Images fed to a deep neural network have in general undergone several handcrafted image signal processing (ISP) operations, all of which have been optimized to produce visually pleasing images. In this work, we investigate the hypothesis that the intermediate representation of visually pleasing images is sub-optimal for downstream computer vision tasks compared to the RAW image representation. We suggest that the operations of the ISP instead should be optimized towards the end task, by learning the parameters of the operations jointly during training. We extend previous works on this topic and propose a new learnable operation that enables an object detector to achieve superior performance when compared to both previous works and traditional RGB images. In experiments on the open PASCALRAW dataset, we empirically confirm our hypothesis.


title: The Euler Characteristic of Finite Categories

abstract: We associate a rational number $\chi(\mathcal{A})$ to every category $\mathcal{A}$ whose object and morphism sets are finite. We show that the assignment $\chi$ is additive under disjoint union and it preserves products. Hence we consider $\chi$ as an Euler measure on a family of categories where this assignment also obeys a version of the inclusion--exclusion principle.


title: Infinite-dimensional port-Hamiltonian systems with a stationary interface

abstract: We consider two systems of two conservation laws that are defined on complementary spatial intervals and coupled by an interface as a single port-Hamiltonian system. In case of a fixed interface position, we characterize the boundary and interface conditions for which the associated port-Hamiltonian operator generates a contraction semigroup. Furthermore, we present sufficient conditions for the exponential stability of the generated $C_0$-semigroup. The results are illustrated by the example of two acoustic waveguides coupled by an interface consisting of some membrane.


title: The Best of Both Worlds: Accurate Global and Personalized Models through Federated Learning with Data-Free Hyper-Knowledge Distillation

abstract: Heterogeneity of data distributed across clients limits the performance of global models trained through federated learning, especially in the settings with highly imbalanced class distributions of local datasets. In recent years, personalized federated learning (pFL) has emerged as a potential solution to the challenges presented by heterogeneous data. However, existing pFL methods typically enhance performance of local models at the expense of the global model's accuracy. We propose FedHKD (Federated Hyper-Knowledge Distillation), a novel FL algorithm in which clients rely on knowledge distillation (KD) to train local models. In particular, each client extracts and sends to the server the means of local data representations and the corresponding soft predictions -- information that we refer to as ``hyper-knowledge". The server aggregates this information and broadcasts it to the clients in support of local training. Notably, unlike other KD-based pFL methods, FedHKD does not rely on a public dataset nor it deploys a generative model at the server. We analyze convergence of FedHKD and conduct extensive experiments on visual datasets in a variety of scenarios, demonstrating that FedHKD provides significant improvement in both personalized as well as global model performance compared to state-of-the-art FL methods designed for heterogeneous data settings.


title: Parikh Automata on Infinite Words

abstract: Parikh automata on finite words were first introduced by Klaedtke and Rue{\ss} [Automata, Languages and Programming, 2003]. In this paper, we introduce several variants of Parikh automata on infinite words and study their expressiveness. We show that one of our new models is equivalent to synchronous blind counter machines introduced by Fernau and Stiebe [Fundamenta Informaticae, 2008]. All our models admit {\epsilon}-elimination, which to the best of our knowledge is an open question for blind counter automata. We then study the classical decision problems of the new automata models.


title: The Conditional Cauchy-Schwarz Divergence with Applications to Time-Series Data and Sequential Decision Making

abstract: The Cauchy-Schwarz (CS) divergence was developed by Pr\'{i}ncipe et al. in 2000. In this paper, we extend the classic CS divergence to quantify the closeness between two conditional distributions and show that the developed conditional CS divergence can be simply estimated by a kernel density estimator from given samples. We illustrate the advantages (e.g., rigorous faithfulness guarantee, lower computational complexity, higher statistical power, and much more flexibility in a wide range of applications) of our conditional CS divergence over previous proposals, such as the conditional KL divergence and the conditional maximum mean discrepancy. We also demonstrate the compelling performance of conditional CS divergence in two machine learning tasks related to time series data and sequential inference, namely time series clustering and uncertainty-guided exploration for sequential decision making.


title: An Experimental Investigation of Cavitation Bulk Nanobubbles Characteristics: Effects of pH and Surface-active Agents

abstract: Understanding the behavior of nanobubbles (NBs) in various aqueous solutions is a challenging task. The present work investigates the effects of various surfactants (i.e., anionic, cationic, and nonionic) and pH medium on bulk NBs formation, size, concentration, bubble size distribution (BSD), zeta potential, and stability. The effect of surfactant was investigated at various concentrations above and below critical micelle concentrations. NBs were created in DI water using a piezoelectric transducer. The stability of NBs was assessed by tracking the change in size and concentration over time. NBs size is small in the neutral medium compared to the other surfactant or pH mediums. The size, concentration, BSD, and stability of NBs are strongly influenced by the zeta potential rather than the solution medium. BSD curve shifts to lower bubble sizes when the magnitude of zeta potential is high in any solution. NBs were observed to exist for a long time, either in pure water, surfactant, or pH solutions. The longevity of NBs is shortened in environments with pH less than 3. Surfactant adsorption on the NBs surface increases with surfactant concentration up to a certain limit, beyond which it declines considerably. The Derjaguin-LandauVerwey-Overbeek (DLVO) theory was used to interpret the NBs stability, which resulted in a total potential energy barrier that is positive and greater than 43.90kBT for pH ranging from 6.0 to 11.0, whereas, for pH below 6, the potential energy barrier essentially vanishes. Moreover, an effort has also been made to elucidate the plausible prospect of ion distribution and its alignment surrounding NBs in cationic and anionic surfactants. The present research will extend the in-depth investigation of NBs for industrial applications involving NBs.


title: A Conditional Abundance Matching Method of Extending Simulated Halo Merger Trees to Resolve Low-Mass Progenitors and Sub-halos

abstract: We present an algorithm to extend subhalo merger trees in a low-resolution dark-matter-only simulation by conditionally matching them to those in a high-resolution simulation. The algorithm is general and can be applied to simulation data with different resolutions using different target variables. We instantiate the algorithm by a case in which trees from ELUCID, a constrained simulation of $(500h^{-1}{\rm Mpc})^3$ volume of the local universe, are extended by matching trees from TNGDark, a simulation with much higher resolution. Our tests show that the extended trees are statistically equivalent to the high-resolution trees in the joint distribution of subhalo quantities and in important summary statistics relevant to modeling galaxy formation and evolution in halos. The extended trees preserve certain information of individual systems in the target simulation, including properties of resolved satellite subhalos, and shapes and orientations of their host halos. With the extension, subhalo merger trees in a cosmological scale simulation are extrapolated to a mass resolution comparable to that in a higher-resolution simulation carried out in a smaller volume, which can be used as the input for (sub)halo-based models of galaxy formation. The source code of the algorithm, and halo merger trees extended to a mass resolution of $\sim 2 \times 10^8 h^{-1}M_\odot$ in the entire ELUCID simulation, are available.


title: Vision Aided Environment Semantics Extraction and Its Application in mmWave Beam Selection

abstract: In this letter, we propose a novel mmWave beam selection method based on the environment semantics extracted from user-side camera images. Specifically, we first define the environment semantics as the spatial distribution of the scatterers that affect the wireless propagation channels and utilize the keypoint detection technique to extract them from the input images. Then, we design a deep neural network with the environment semantics as the input that can output the optimal beam pairs at the mobile station (MS) and the base station (BS). Compared with the existing beam selection approaches that directly use images as the input, the proposed semantic-based method can explicitly obtain the environmental features that account for the propagation of wireless signals, thus reducing the storage and computational burden. Simulation results show that the proposed method can precisely estimate the location of the scatterers and outperform the existing works based on computer vision or light detection and ranging (LIDAR).


title: Soft Sensing Regression Model: from Sensor to Wafer Metrology Forecasting

abstract: The semiconductor industry is one of the most technology-evolving and capital-intensive market sectors. Effective inspection and metrology are necessary to improve product yield, increase product quality and reduce costs. In recent years, many semiconductor manufacturing equipments are equipped with sensors to facilitate real-time monitoring of the production process. These production-state and equipment-state sensor data provide an opportunity to practice machine-learning technologies in various domains, such as anomaly/fault detection, maintenance scheduling, quality prediction, etc. In this work, we focus on the task of soft sensing regression, which uses sensor data to predict impending inspection measurements that used to be measured in wafer inspection and metrology systems. We proposed an LSTM-based regressor and designed two loss functions for model training. Although engineers may look at our prediction errors in a subjective manner, a new piece-wise evaluation metric was proposed for assessing model accuracy in a mathematical way. The experimental results demonstrated that the proposed model can achieve accurate and early prediction of various types of inspections in complicated manufacturing processes.


title: Expectations for time-delay measurements in active galactic nuclei with the Vera Rubin Observatory

abstract: The Vera Rubin Observatory will provide an unprecedented set of time-dependent observations of the sky. The planned Legacy Survey of Space and Time (LSST) operating for 10 years will provide dense lightcurves for thousands of active galactic nuclei (AGN) in Deep Drilling Fields (DDFs) and less dense lightcurves for millions of AGN. We model the prospects for measuring time delays for emission lines with respect to the continuum, using these data. We model the artificial lightcurves using Timmer-Koenig algorithm, we use the exemplary cadence to sample them, we supplement lightcurves with the expected contamination by the strong emission lines (Hbeta, Mg II and CIV as well as with Fe II pseudo-continuum and the starlight). We choose the suitable photometric bands appropriate for the redshift and compare the assumed line time delay with the recovered time delay for 100 statistical realizations of the light curves. We show that time delays for emission lines can be well measured from the Main Survey for the bright tail of the quasar distribution (about 15% of all sources) with the accuracy within 1 sigma error, for DDFs results for fainter quasars are also reliable when all 10 years of data are used. There are also some prospects to measure the time delays for the faintest quasars at the smallest redshifts from the first two years of data, and eventually even from the first season. The entire quasar population will allow obtaining results of apparently high accuracy but in our simulations, we see a systematic offset between the assumed and recovered time delay depending on the redshift and source luminosity which will not disappear even in the case of large statistics. Such a problem might affect the slope of the radius-luminosity relation and cosmological applications of quasars if simulations correcting for such effects are not performed.


title: Modified Erd\H{o}s-Ginzburg-Ziv constants for $\mathbb{Z}_2^d$

abstract: Let $G$ be a finite abelian group written additively, and let $r$ be a multiple of its exponent. The modified Erd\H{o}s-Ginzburg-Ziv constant $\mathsf{s}_r'(G)$ is the smallest integer $s$ such that every zero-sum sequence of length $s$ over $G$ has a zero-sum subsequence of length $r$. We find exact values of $\mathsf{s}_{2k}'(\mathbb{Z}_2^d)$ for $d \leq 2k+1$.


title: Modular algorithms for Gross-Stark units and Stark-Heegner points

abstract: In recent work, Darmon, Pozzi and Vonk explicitly construct a modular form whose spectral coefficients are $p$-adic logarithms of Gross-Stark units and Stark-Heegner points. Here we describe how this construction gives rise to a practical algorithm for explicitly computing these logarithms to specified precision, and how to recover the exact values of the Gross-Stark units and Stark-Heegner points from them. Key tools are overconvergent modular forms, reduction theory of quadratic forms and Newton polygons. As an application, we tabulate Brumer-Stark units in narrow Hilbert class fields of real quadratic fields with discriminants up to $10000$, for primes less than $20$, as well as Stark-Heegner points on elliptic curves.


title: Leveraging driver vehicle and environment interaction: Machine learning using driver monitoring cameras to detect drunk driving

abstract: Excessive alcohol consumption causes disability and death. Digital interventions are promising means to promote behavioral change and thus prevent alcohol-related harm, especially in critical moments such as driving. This requires real-time information on a person's blood alcohol concentration (BAC). Here, we develop an in-vehicle machine learning system to predict critical BAC levels. Our system leverages driver monitoring cameras mandated in numerous countries worldwide. We evaluate our system with n=30 participants in an interventional simulator study. Our system reliably detects driving under any alcohol influence (area under the receiver operating characteristic curve [AUROC] 0.88) and driving above the WHO recommended limit of 0.05g/dL BAC (AUROC 0.79). Model inspection reveals reliance on pathophysiological effects associated with alcohol consumption. To our knowledge, we are the first to rigorously evaluate the use of driver monitoring cameras for detecting drunk driving. Our results highlight the potential of driver monitoring cameras and enable next-generation drunk driver interaction preventing alcohol-related harm.


title: Informing policy via dynamic models: Cholera in Haiti

abstract: Public health decisions must be made about when and how to implement interventions to control an infectious disease epidemic. These decisions should be informed by data on the epidemic as well as current understanding about the transmission dynamics. Such decisions can be posed as statistical questions about scientifically motivated dynamic models. Thus, we encounter the methodological task of building credible, data-informed decisions based on stochastic, partially observed, nonlinear dynamic models. This necessitates addressing the tradeoff between biological fidelity and model simplicity, and the reality of misspecification for models at all levels of complexity. We assess current methodological approaches to these issues via a case study of the 2010-2019 cholera epidemic in Haiti. We consider three dynamic models developed by expert teams to advise on vaccination policies. We evaluate previous methods used for fitting these models, and we demonstrate modified data analysis strategies leading to improved statistical fit. Specifically, we present approaches for diagnosing model misspecification and the consequent development of improved models. Additionally, we demonstrate the utility of recent advances in likelihood maximization for high-dimensional nonlinear dynamic models, enabling likelihood-based inference for spatiotemporal incidence data using this class of models. Our workflow is reproducible and extendable, facilitating future investigations of this disease system.


title: Towards Quantification of Assurance for Learning-enabled Components

abstract: Perception, localization, planning, and control, high-level functions often organized in a so-called pipeline, are amongst the core building blocks of modern autonomous (ground, air, and underwater) vehicle architectures. These functions are increasingly being implemented using learning-enabled components (LECs), i.e., (software) components leveraging knowledge acquisition and learning processes such as deep learning. Providing quantified component-level assurance as part of a wider (dynamic) assurance case can be useful in supporting both pre-operational approval of LECs (e.g., by regulators), and runtime hazard mitigation, e.g., using assurance-based failover configurations. This paper develops a notion of assurance for LECs based on i) identifying the relevant dependability attributes, and ii) quantifying those attributes and the associated uncertainty, using probabilistic techniques. We give a practical grounding for our work using an example from the aviation domain: an autonomous taxiing capability for an unmanned aircraft system (UAS), focusing on the application of LECs as sensors in the perception function. We identify the applicable quantitative measures of assurance, and characterize the associated uncertainty using a non-parametric Bayesian approach, namely Gaussian process regression. We additionally discuss the relevance and contribution of LEC assurance to system-level assurance, the generalizability of our approach, and the associated challenges.


title: Analytic lattice cohomology of isolated curve singularities

abstract: We construct a lattice cohomology ${\mathbb H}^*(C,o)=\oplus_{q\geq 0}{\mathbb H}^q(C,o)$ and a graded root ${\mathfrak R}(C,o)$ to any complex isolated curve singularity $(C,o)$. Each ${\mathbb H}^q(C,o)$ is a ${\mathbb Z}$-graded ${\mathbb Z}[U]$-module. The Euler characteristic of ${\mathbb H}^*(C,o)$ is the delta-invariant of $(C,o)$. The construction is based on the multivariable Hilbert series of the multifiltration provided by valuations of the normalization. Several examples are discussed, e.g. Gorenstein curves (where an additional symmetry is established), plane curves (in particular, Newton non-degenerate ones), ordinary $r$-tuples. We also prove that a flat deformation $(C_t,o)_{t\in ({\mathbb C},0)}$ of isolated curve singularities induces an explicit degree zero graded ${\mathbb Z}[U]$-module morphism ${\mathbb H}^0(C_{t=0},o)\to {\mathbb H}^0(C_{t\not=0},o)$, and a graded (graph) map of degree zero at the level of graded roots ${\mathfrak R}(C_{t\not=0},o)\to{\mathfrak R}(C_{t=0},o)$. In the treatment of the deformation functor we need a second construction of the lattice cohomology in terms of the system of linear subspace arrangements associated with the above filtration.


title: Soft X-ray spectro-ptychography on boron nitride nanotubes, carbon nanotubes and permalloy nanorods

abstract: Spectro-ptychography offers improved spatial resolution and additional phase spectral information relative to that provided by scanning transmission X-ray microscopes (STXM). However, carrying out ptychography at the lower range of soft X-ray energies (e.g., below 200 eV to 600 eV) on samples with weakly scattering signals can be challenging. We present soft X-ray ptychography results at energies as low as 180 eV and illustrate the capabilities with results from permalloy nanorods (Fe 2p), carbon nanotubes (C 1s), and boron nitride bamboo nanostructures (B 1s, N1s). We describe optimization of low energy X-ray spectro-ptychography and discuss important challenges associated with measurement approaches, reconstruction algorithms, and their effects on the reconstructed images. A method for evaluating the increase in radiation dose when using overlapping sampling is presented.


title: Closed formula for the transport of micro-nano-particle across model porous media

abstract: In the last decade the Fick-Jacobs approximation has been exploited to capture the transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that for fore-aft symmetric channels the flux of non-interacting systems is fully captured by its linear response regime. For this case we derive a very simple formula that captures the correct trends and that can be exploited as a simple tool to design experiments or simulations. Finally, we show that higher order corrections in the flux may appear for non-symmetric channels.


title: SuperScaler: Supporting Flexible DNN Parallelization via a Unified Abstraction

abstract: With the growing model size, deep neural networks (DNN) are increasingly trained over massive GPU accelerators, which demands a proper parallelization plan that transforms a DNN model into fine-grained tasks and then schedules them to GPUs for execution. Due to the large search space, the contemporary parallelization plan generators often rely on empirical rules that couple transformation and scheduling, and fall short in exploring more flexible schedules that yield better memory usage and compute efficiency. This tension can be exacerbated by the emerging models with increasing complexity in their structure and model size. SuperScaler is a system that facilitates the design and generation of highly flexible parallelization plans. It formulates the plan design and generation into three sequential phases explicitly: model transformation, space-time scheduling, and data dependency preserving. Such a principled approach decouples multiple seemingly intertwined factors and enables the composition of highly flexible parallelization plans. As a result, SuperScaler can not only generate empirical parallelization plans, but also construct new plans that achieve up to 3.5X speedup compared to state-of-the-art solutions like DeepSpeed, Megatron and Alpa, for emerging DNN models like Swin-Transformer and AlphaFold2, as well as well-optimized models like GPT-3.


title: Advances in the Implementation of the Exactly Energy Conserving Semi-Implicit (ECsim) Particle in Cell Method

abstract: The Energy Conserving semi-implicit method (ECsim), presented by Lapenta in 2017, is a Particle in Cell (PIC) algorithm for the simulation of plasmas. Energy conservation is achieved within a semi-implicit formulation that does not require any non-linear solver. A mass matrix is introduced to express linearly the particle-field coupling. With the mass matrix the algorithm preserves energy conservation to machine precision. The construction of the mass matrix is the central nature of the method and also the main cost of the computational cycle. We analyze here three methods that modify the construction of the mass matrix. First, we consider how the sub-cycling of the particle motion modifies the mass matrix. Second, we introduce a form of smoothing that reduces the noise while retaining exact energy conservation. Finally, we discuss an approximation of the mass matrix that transform the ECsim scheme to the implicit moment method.


title: Adapting a Language Model While Preserving its General Knowledge

abstract: Domain-adaptive pre-training (or DA-training for short), also known as post-training, aims to train a pre-trained general-purpose language model (LM) using an unlabeled corpus of a particular domain to adapt the LM so that end-tasks in the domain can give improved performances. However, existing DA-training methods are in some sense blind as they do not explicitly identify what knowledge in the LM should be preserved and what should be changed by the domain corpus. This paper shows that the existing methods are suboptimal and proposes a novel method to perform a more informed adaptation of the knowledge in the LM by (1) soft-masking the attention heads based on their importance to best preserve the general knowledge in the LM and (2) contrasting the representations of the general and the full (both general and domain knowledge) to learn an integrated representation with both general and domain-specific knowledge. Experimental results will demonstrate the effectiveness of the proposed approach.


title: Tier Balancing: Towards Dynamic Fairness over Underlying Causal Factors

abstract: The pursuit of long-term fairness involves the interplay between decision-making and the underlying data generating process. In this paper, through causal modeling with a directed acyclic graph (DAG) on the decision-distribution interplay, we investigate the possibility of achieving long-term fairness from a dynamic perspective. We propose Tier Balancing, a technically more challenging but more natural notion to achieve in the context of long-term, dynamic fairness analysis. Different from previous fairness notions that are defined purely on observed variables, our notion goes one step further, capturing behind-the-scenes situation changes on the unobserved latent causal factors that directly carry out the influence from the current decision to the future data distribution. Under the specified dynamics, we prove that in general one cannot achieve the long-term fairness goal only through one-step interventions. Furthermore, in the effort of approaching long-term fairness, we consider the mission of "getting closer to" the long-term fairness goal and present possibility and impossibility results accordingly.


title: Matchings under distance constraints II

abstract: This paper introduces the \emph{$d$-distance $b$-matching problem}, in which we are given a bipartite graph $G=(S,T;E)$ with $S=\{s_1,\dots,s_n\}$, a weight function on the edges, an integer $d\in\mathbb{Z}_+$ and a degree bound function $b:S\cup T\to\mathbb{Z}_+$. The goal is to find a maximum-weight subset $M\subseteq E$ of the edges satisfying the following two conditions: 1) the degree of each node $v\in S\cup T$ is at most $b(v)$ in $M$, 2) if $s_it,s_jt\in M$, then $|i-j|\geq d$. In the cyclic version of the problem, the nodes in $S$ are considered to be in cyclic order. We get back the \emph{(cyclic) $d$-distance matching problem} when $b(s) = 1$ for $s\in S$ and $b(t) = \infty$ for $t\in T$. We prove that the $d$-distance matching problem is APX-hard even in the unweighted case. We show that $(2-\frac{1}{d})$ is a tight upper bound on the integrality gap of the natural integer programming model for the cyclic $d$-distance $b$-matching problem provided that $(2d-1)$ divides the size of $S$. For the non-cyclic case, the integrality gap is shown to be at most $(2-\frac{2}{d})$. The proofs give approximation algorithms with guarantees matching these bounds, and also improve the best known algorithms for the (cyclic) $d$-distance matching problem. In a related problem, our goal is to find a permutation of $S$ maximizing the weight of the optimal $d$-distance $b$-matching. This problem can be solved in polynomial time for the (cyclic) $d$-distance matching problem -- even though the (cyclic) $d$-distance matching problem itself is NP-hard and also hard to approximate arbitrarily. For (cyclic) $d$-distance $b$-matchings, however, we prove that finding the best permutation is NP-hard even if $b\equiv 2$ or $d=2$, and we give $e$-approximation algorithms.


title: Note on Milnor numbers of irreducible germs

abstract: Let $(\bf {V,0})\subset (\mathbb{C}^n,0)$ be a germ of a complex hypersurface and let $f: (\mathbb{C}^n,0)\to(\mathbb{C}^n,0)$ be a germ of a finite holomorphic mapping. If germs $(\bf {V,0})$ and ${\bf W}:=(F^{-1}(\bf{ V})),0)$ are irreducible and with isolated singularities, then $$\mu(F^{-1}(\bf{ V}))\ge \mu(\bf {V}),$$ where $\mu$ denotes the Milnor number.


title: Monitoring of the heart movements using a FMCW radar and correlation with an ECG

abstract: Monitoring the activity of the heart is important for diagnosing and preventing cardiovascular diseases. The electrocardiogram (ECG) is the gold standard for diagnosing such diseases. It monitors the heart's electrical activity, and while this is highly correlated with the cardiac mechanical activity, it does not provide all the information. Other sensors such as the echocardiograph are able to monitor the heart's movements, but such tools are expensive and hard to operate. Therefore, contactless monitoring of the heart using RF sensing has gained interest in recent years. In this paper, we describe a process to extract the movements of the heart from millimeter wave radar with high accuracy, and thus we provide a noninvasive and affordable way to monitor cardiac movements. We then demonstrate the correlation between the observed movements and the ECG. Furthermore, we propose an algorithm to synchronize the ECG signal and the processed signal from the radar sensor. The results we obtained provide insights on the mechanical activity of the heart, which could assist cardiologists in their diagnoses.