title: Direct sampling for recovering sound soft scatterers from point source measurements

abstract: In this paper, we consider the inverse problem of recovering a sound soft scatterer from the measured scattered field. The scattered field is assumed to be induced by a point source on a curve/surface that is known. Here we will propose and analyze new direct sampling methods for this problem. The first method we consider uses a far-field transformation of the near-field data which will allow us to derive explicit bounds in the resolution analysis for the direct sampling method's imaging functional. Two direct sampling methods will be studied using the far-field transformation. For these imaging functionals, we will use the Funk-Hecke identities to study the resolution analysis. We will also study a direct sampling method for the case of the given Cauchy data. Numerical examples are given to show the applicability of the new imaging functionals for recovering a sound soft scatterer in 2D.


title: Chern Invariants of Topological Continua; a Self-Consistent Nonlocal Hydrodynamic Model

abstract: The bulk-edge correspondence is a fundamental principle of topological wave physics, which states that the difference in gap Chern numbers between the interfaced materials is equal to the net number of topological edge modes. Although this principle works well for periodic photonic topological insulators, difficulties arise in the case of continuous systems with no intrinsic periodicity, due to the absence of a finite Brillouin zone, which may lead to an ill-behaved response for diverging wavenumbers. This problem has been solved previously by introducing an ad hoc material model including a spatial cutoff wavenumber. However, this method introduces other difficulties, such as an unphysical response at large wavenumbers, and the need to interpolate the interfaced materials permittivity functions. In this work, we show that the inclusion of realistic spatial dispersion (e.g., hydrodynamic nonlocality) formally restores the bulk-edge correspondence by avoiding an unphysical response at large wavenumbers, forming a complete, self-consistent model. This model clarifies the subtle and rich topological wave physics of continuous media.


title: RegGuard: Leveraging CPU Registers for Mitigation of Control- and Data-Oriented Attacks

abstract: CPU registers are small discrete storage units, used to hold temporary data and instructions within the CPU. Registers are not addressable in the same way memory is, which makes them immune from memory attacks and manipulation by other means. In this paper, we take advantage of this to provide a protection mechanism for critical program data; both active local variables and control objects on the stack. This protection effectively eliminates the threat of control- and data-oriented attacks, even by adversaries with full knowledge of the active stack. Our solution RegGuard, is a compiler register allocation strategy that utilises the available CPU registers to hold critical variables during execution. Unlike conventional allocations schemes, RegGuard prioritises the security significance of a program variable over its expected performance gain. Our scheme can deal effectively with saved registers to the stack, i.e., when the compiler needs to free up registers to make room for the variables of a new function call. With RegGuard, critical data objects anywhere on the entire stack are effectively protected from corruption, even by adversaries with arbitrary read and write access. While our primary design focus is on security, performance is very important for a scheme to be adopted in practice. RegGuard is still benefiting from the performance gain normally associated with register allocations, and the overhead is within a few percent of other unsecured register allocation schemes for most cases. We present detailed experiments that showcase the performance of RegGuard using different benchmark programs and the C library on ARM64 platform.


title: Pick-and-Mix Information Operators for Probabilistic ODE Solvers

abstract: Probabilistic numerical solvers for ordinary differential equations compute posterior distributions over the solution of an initial value problem via Bayesian inference. In this paper, we leverage their probabilistic formulation to seamlessly include additional information as general likelihood terms. We show that second-order differential equations should be directly provided to the solver, instead of transforming the problem to first order. Additionally, by including higher-order information or physical conservation laws in the model, solutions become more accurate and more physically meaningful. Lastly, we demonstrate the utility of flexible information operators by solving differential-algebraic equations. In conclusion, the probabilistic formulation of numerical solvers offers a flexible way to incorporate various types of information, thus improving the resulting solutions.


title: Gapless Dirac magnons in CrCl$_{3}$

abstract: Bosonic Dirac materials are testbeds for dissipationless spin-based electronics. In the quasi two-dimensional honeycomb lattice of CrX$_{3}$ (X=Cl, Br, I), Dirac magnons have been predicted at the crossing of acoustical and optical spin waves, analogous to Dirac fermions in graphene. Here we show that, distinct from CrBr$_{3}$ and CrI$_{3}$, gapless Dirac magnons are present in bulk CrCl$_{3}$, with inelastic neutron scattering intensity at low temperatures approaching zero at the Dirac $K$ point. Upon warming, magnon-magnon interactions induce strong renormalization and decreased lifetimes, with a ~25% softening of the upper magnon branch intensity from 5 to 50 K, though magnon features persist well above T$_{N}$. Moreover, an unusual negative thermal expansion (NTE) of the $a$-axis lattice constant and anomalous phonon behavior are observed below 50 K, indicating magnetoelastic and spin-phonon coupling arising from an increase in the in-plane spin correlations that begins tens of Kelvin above T$_{N}$.


title: Closed-loop Feedback Registration for Consecutive Images of Moving Flexible Targets

abstract: Advancement of imaging techniques enables consecutive image sequences to be acquired for quality monitoring of manufacturing production lines. Registration for these image sequences is essential for in-line pattern inspection and metrology, e.g., in the printing process of flexible electronics. However, conventional image registration algorithms cannot produce accurate results when the images contain many similar and deformable patterns in the manufacturing process. Such a failure originates from a fact that the conventional algorithms only use the spatial and pixel intensity information for registration. Considering the nature of temporal continuity and consecution of the product images, in this paper, we propose a closed-loop feedback registration algorithm for matching and stitching the deformable printed patterns on a moving flexible substrate. The algorithm leverages the temporal and spatial relationships of the consecutive images and the continuity of the image sequence for fast, accurate, and robust point matching. Our experimental results show that our algorithm can find more matching point pairs with a lower root mean squared error (RMSE) compared to other state-of-the-art algorithms while offering significant improvements to running time.


title: Almost coherent modules and almost coherent sheaves

abstract: We review the theory of almost coherent modules that was introduced in "Almost Ring Theory" by Gabber and Ramero. Then we globalize it by developing a new theory of almost coherent sheaves on schemes and on a class of "nice" formal schemes. We show that these sheaves satisfy many properties similar to usual coherent sheaves, i.e. the Almost Proper Mapping Theorem, the Formal GAGA, etc. We also construct an almost version of the Grothendieck twisted image functor $f^!$ and verify its properties. Lastly, we study sheaves of $p$-adic nearby cycles on admissible formal models of rigid spaces and show that these sheaves provide examples of almost coherent sheaves. This gives a new proof of the finiteness result for \'etale cohomology of proper rigid spaces obtained before in the work of Peter Scholze "$p$-adic Hodge Theory For Rigid-Analytic Varieties".


title: SciXGen: A Scientific Paper Dataset for Context-Aware Text Generation

abstract: Generating texts in scientific papers requires not only capturing the content contained within the given input but also frequently acquiring the external information called \textit{context}. We push forward the scientific text generation by proposing a new task, namely \textbf{context-aware text generation} in the scientific domain, aiming at exploiting the contributions of context in generated texts. To this end, we present a novel challenging large-scale \textbf{Sci}entific Paper Dataset for Conte\textbf{X}t-Aware Text \textbf{Gen}eration (SciXGen), consisting of well-annotated 205,304 papers with full references to widely-used objects (e.g., tables, figures, algorithms) in a paper. We comprehensively benchmark, using state-of-the-arts, the efficacy of our newly constructed SciXGen dataset in generating description and paragraph. Our dataset and benchmarks will be made publicly available to hopefully facilitate the scientific text generation research.


title: Reduced Basis Approximations of Parameterized Dynamical Partial Differential Equations via Neural Networks

abstract: Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a time-stepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. We demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.


title: Air, helium and water leakage in rubber O-ring seals with application to syringes

abstract: We study the leakage of fluids (liquids or gases) in syringes with glass barrel, steel plunger and rubber O-ring stopper. The leakrate depends on the interfacial surface roughness and on the viscoelastic properties of the rubber. Random surface roughness is produced by sandblasting the rubber O-rings. We present a very simple theory for gas flow which takes into account both the diffusive and ballistic flow. The theory shows that the interfacial fluid flow (leakage) channels are so narrow that the gas flow is mainly ballistic (the so called Knudsen limit). We compare the leakrate obtained using air and helium. For barrels filled with water we observe no leakage even if leakage occurs for gases. We interpret this as resulting from capillary (Laplace pressure or surface energy) effects.


title: Learning controllers for performance through LMI regions

abstract: In an open-loop experiment, an input sequence is applied to an unknown linear time-invariant system (in continuous or discrete time) affected also by an unknown-but-bounded disturbance sequence (with an energy or instantaneous bound); the corresponding state sequence is measured. The goal is to design directly from the input and state sequences a controller that enforces a certain performance specification on the transient behaviour of the unknown system. The performance specification is expressed through a subset of the complex plane where closed-loop eigenvalues need to belong, a so called LMI region. For this control design problem, we provide here convex programs to enforce the performance specification from data in the form of linear matrix inequalities (LMI). For generic LMI regions, these are sufficient conditions to assign the eigenvalues within the LMI region for all possible dynamics consistent with data, and become necessary and sufficient conditions for special LMI regions. In this way, we extend classical model-based conditions from a seminal work in the literature to the setting of data-driven control from noisy data. Through two numerical examples, we investigate how these data-based conditions compare with each other.


title: Contrastive Document Representation Learning with Graph Attention Networks

abstract: Recent progress in pretrained Transformer-based language models has shown great success in learning contextual representation of text. However, due to the quadratic self-attention complexity, most of the pretrained Transformers models can only handle relatively short text. It is still a challenge when it comes to modeling very long documents. In this work, we propose to use a graph attention network on top of the available pretrained Transformers model to learn document embeddings. This graph attention network allows us to leverage the high-level semantic structure of the document. In addition, based on our graph document model, we design a simple contrastive learning strategy to pretrain our models on a large amount of unlabeled corpus. Empirically, we demonstrate the effectiveness of our approaches in document classification and document retrieval tasks.


title: The circumbinary rings of GG Carinae: indications of disc eccentricity growth in the B[e] supergiant's atomic emission lines

abstract: B[e] supergiants have unusual circumstellar environments which may include thin, concentric rings displaying atomic line emission. GG Carinae is a B[e] supergiant binary which exhibits such a geometry in its circumbinary environment. We study atomic emission lines arising from GG Carinae's circumbinary disc in FEROS spectra collected between 1998 and 2015. We find that semi-forbidden Fe\,II] and permitted Ca\,II emission are formed in the same thin circumbinary ring previously reported to have forbidden [O\,I] and [Ca\,II] emission. We find that there are two circumbinary rings orbiting with projected velocities of $84.6\pm1.0$\,km\,s$^{-1}$ and $27.3\pm0.6$\,km\,s$^{-1}$. Deprojecting these velocities from the line-of-sight, and using updated binary masses presented by \cite{Porter2021GGPhotometry}, we find that the radii of the circumbinary rings are $2.8^{+0.9}_{-1.1}$\,AU and $27^{+9}_{-10}$\,AU for the inner ring and outer ring respectively. We find evidence of subtle dynamical change in the inner circumbinary ring over the 17 years spanned by the data, manifesting in variability in the ratio of the intensity of the blueshifted peak to the redshifted peak of its emission lines and the central velocity becoming more blueshifted. We perform smoothed-particle hydrodynamic simulations of the system which suggest that these observed changes are consistent with pumping of the eccentricity of a radially thin circumbinary ring by the inner binary. We find a systemic velocity of the GG Carinae system of $-23.2 \pm 0.4$\,km\,s$^{-1}$.


title: An Open Natural Language Processing Development Framework for EHR-based Clinical Research: A case demonstration using the National COVID Cohort Collaborative (N3C)

abstract: While we pay attention to the latest advances in clinical natural language processing (NLP), we can notice some resistance in the clinical and translational research community to adopt NLP models due to limited transparency, interpretability, and usability. In this study, we proposed an open natural language processing development framework. We evaluated it through the implementation of NLP algorithms for the National COVID Cohort Collaborative (N3C). Based on the interests in information extraction from COVID-19 related clinical notes, our work includes 1) an open data annotation process using COVID-19 signs and symptoms as the use case, 2) a community-driven ruleset composing platform, and 3) a synthetic text data generation workflow to generate texts for information extraction tasks without involving human subjects. The corpora were derived from texts from three different institutions (Mayo Clinic, University of Kentucky, University of Minnesota). The gold standard annotations were tested with a single institution's (Mayo) ruleset. This resulted in performances of 0.876, 0.706, and 0.694 in F-scores for Mayo, Minnesota, and Kentucky test datasets, respectively. The study as a consortium effort of the N3C NLP subgroup demonstrates the feasibility of creating a federated NLP algorithm development and benchmarking platform to enhance multi-institution clinical NLP study and adoption. Although we use COVID-19 as a use case in this effort, our framework is general enough to be applied to other domains of interest in clinical NLP.


title: Marital Stability With Committed Couples: A Revealed Preference Analysis

abstract: We present a revealed preference characterization of marital stability where some couples are committed. A couple is committed if they can only divorce upon mutual consent. We provide theoretical insights into the potential of the characterization for identifying intrahousehold consumption patterns. We show that when there is no price variation for private goods between potential couples, it is only possible to identify intrahousehold resource allocations for non-committed couples. Simulation exercises using household data drawn from the Longitudinal Internet Studies for the Social Sciences (LISS) panel support our theoretical findings. Our results show that in the presence of price variation, the empirical implications of marital stability can be used for identifying household consumption allocations for both committed and non-committed couples.


title: Symmetry breaking for ground states of biharmonic NLS via Fourier extension estimates

abstract: We consider ground states solutions $u \in H^2(\mathbb{R}^N)$ of biharmonic (fourth-order) nonlinear Schr\"odinger equations of the form $$ \Delta^2 u + 2a \Delta u + b u - |u|^{p-2} u = 0 \quad \mbox{in $\mathbb{R}^N$} $$ with positive constants $a, b > 0$ and exponents $2 < p < 2^*$, where $2^* = \frac{2N}{N-4}$ if $N > 4$ and $2^* = \infty$ if $N \leq 4$. By exploiting a connection to the adjoint Stein--Tomas inequality on the unit sphere and by using trial functions due to Knapp, we prove a general symmetry breaking result by showing that all ground states $u\in H^2(\mathbb{R}^N)$ in dimension $N \geq 2$ fail to be radially symmetric for all exponents $2 < p < \frac{2N+2}{N-1}$ in a suitable regime of $a,b>0$. As applications of our main result, we also prove symmetry breaking for a minimization problem with constrained $L^2$-mass and for a related problem on the unit ball in $\mathbb{R}^N$ subject to Dirichlet boundary conditions.


title: Adversarial attacks against Bayesian forecasting dynamic models

abstract: The last decade has seen the rise of Adversarial Machine Learning (AML). This discipline studies how to manipulate data to fool inference engines, and how to protect those systems against such manipulation attacks. Extensive work on attacks against regression and classification systems is available, while little attention has been paid to attacks against time series forecasting systems. In this paper, we propose a decision analysis based attacking strategy that could be utilized against Bayesian forecasting dynamic models.


title: Style Agnostic 3D Reconstruction via Adversarial Style Transfer

abstract: Reconstructing the 3D geometry of an object from an image is a major challenge in computer vision. Recently introduced differentiable renderers can be leveraged to learn the 3D geometry of objects from 2D images, but those approaches require additional supervision to enable the renderer to produce an output that can be compared to the input image. This can be scene information or constraints such as object silhouettes, uniform backgrounds, material, texture, and lighting. In this paper, we propose an approach that enables a differentiable rendering-based learning of 3D objects from images with backgrounds without the need for silhouette supervision. Instead of trying to render an image close to the input, we propose an adversarial style-transfer and domain adaptation pipeline that allows to translate the input image domain to the rendered image domain. This allows us to directly compare between a translated image and the differentiable rendering of a 3D object reconstruction in order to train the 3D object reconstruction network. We show that the approach learns 3D geometry from images with backgrounds and provides a better performance than constrained methods for single-view 3D object reconstruction on this task.


title: Random geometric graphs and the spherical Wishart matrix

abstract: We consider the random geometric graph on $n$ vertices drawn uniformly from a $d$--dimensional sphere. We focus on the sparse regime, when the expected degree is constant independent of $d$ and $n$. We show that, when $d$ is larger than $n$ by logarithmic factors, this graph is comparable to the Erd\H{o}s--R\'enyi random graph of the same edge density in the \emph{inclusion divergence} between the graph laws. This divergence functions in certain ways like a relaxation of the total variation distance, but is strong enough to distinguish Erd\H{o}s--R\'enyi graphs of different densities with a higher resolution than the total variation distance. To do the analysis, we derive some exact statistics of the \emph{spherical Wishart matrix}, the Gram matrix of $n$ independent uniformly random $d$--dimensional spherical vectors. In particular we give expressions for the characteristic function of the spherical Wishart matrix which are well--approximated using steepest descent.


title: Possible X-ray Quasi-Periodic Eruptions in a Tidal Disruption Event Candidate

abstract: X-ray Quasi-Periodic Eruptions (QPEs) are a recently discovered phenomenon associated with supermassive black holes at the centers of galaxies. They are high amplitude soft X-ray flares that recur on timescales of hours, but what causes these flares remains uncertain. In the two years since their original discovery, four known QPE-hosting galaxies have been found, with varying properties and levels of activity. We have conducted a blind algorithm-assisted search of the XMM-Newton Source Catalog and found a fifth QPE candidate, XMMSL1 J024916.6-041244. This is a star-forming galaxy hosting a relatively low-mass nuclear black hole, and has previously been identified as a Tidal Disruption Event candidate. An XMM-Newton pointed observation of the source in 2006 exhibited nearly two QPE-like flares in soft X-rays, and, unlike in other QPE sources, there are hints of corresponding dips in the UV light curves. Afterwards, a series of Swift observations observed the rapid dimming of the source; thereafter, in August 2021, we triggered a second XMM-Newton observation, which revealed that the source is detected, but the QPEs are no longer present. Here we report on (I) the strategy we used to systematically search through XMM-Newton archival data; (II) the properties of J0249 and its QPE flares; and (III) the relative behaviors and properties of the QPE sample to date, now 5 members large.


title: Biracks and Switch Braid Quivers

abstract: We consider birack and switch colorings of braids. We define a switch structure on the set of permutation representations of the braid group and consider when such a representation is a switch automorphism. We define quiver-valued invariants of braids using finite switches and biracks and use these to categorify the birack 2-cocycle invariant for braids. We obtain new polynomial invariants of braids via decategorification of these quivers.


title: Noisy intermediate-scale quantum computing algorithm for solving an $n$-vertex MaxCut problem with log($n$) qubits

abstract: Quantum computers are devices, which allow more efficient solutions of problems as compared to their classical counterparts. As the timeline to developing a quantum-error corrected computer is unclear, the quantum computing community has dedicated much attention to developing algorithms for currently available noisy intermediate-scale quantum computers (NISQ). Thus far, within NISQ, optimization problems are one of the most commonly studied and are quite often tackled with the quantum approximate optimization algorithm (QAOA). This algorithm is best known for computing graph partitions with a maximal separation of edges (MaxCut), but can easily calculate other problems related to graphs. Here, I present a novel quantum optimization algorithm, which uses exponentially less qubits as compared to the QAOA while requiring a significantly reduced number of quantum operations to solve the MaxCut problem. Such an improved performance allowed me to partition graphs with 32 nodes on publicly available 5 qubit gate-based quantum computers without any preprocessing such as division of the graph into smaller subgraphs. These results represent a 40% increase in graph size as compared to state-of-art experiments on gate-based quantum computers such as Google Sycamore. The obtained lower bound is 54.9% on the solution for actual hardware benchmarks and 77.6% on ideal simulators of quantum computers. Furthermore, large-scale optimization problems represented by graphs of a 128 nodes are tackled with simulators of quantum computers, again without any predivision into smaller subproblems and a lower solution bound of 67.9% is achieved. The study presented here paves way to using powerful genetic optimizer in synergy with quantum computers


title: The Galois module structure of holomorphic poly-differentials and Riemann-Roch spaces

abstract: Suppose $X$ is a smooth projective geometrically irreducible curve over a perfect field $k$ of positive characteristic $p$. Let $G$ be a finite group acting faithfully on $X$ over $k$ such that $G$ has non-trivial, cyclic Sylow $p$-subgroups. If $E$ is a $G$-invariant Weil divisor on $X$ with $\mathrm{deg}(E)> 2g(X)-2$, we prove that the decomposition of $\mathrm{H}^0(X,\mathcal{O}_X(E))$ into a direct sum of indecomposable $kG$-modules is uniquely determined by the class of $E$ modulo $G$-invariant principal divisors, together with the ramification data of the cover $X\to X/G$. The latter is given by the lower ramification groups and the fundamental characters of the closed points of $X$ that are ramified in the cover. As a consequence, we obtain that if $m>1$ and $g(X)\ge 2$, then the $kG$-module structure of $\mathrm{H}^0(X,\Omega_X^{\otimes m})$ is uniquely determined by the class of a canonical divisor on $X/G$ modulo principal divisors, together with the ramification data of $X\to X/G$. This extends to arbitrary $m > 1$ the $m = 1$ case treated by the first author with T. Chinburg and A. Kontogeorgis. We discuss applications to the tangent space of the global deformation functor associated to $(X,G)$ and to congruences between prime level cusp forms in characteristic $0$. In particular, we complete the description of the precise $k\mathrm{PSL}(2,\mathbb{F}_\ell)$-module structure of all prime level $\ell$ cusp forms of even weight in characteristic $p=3$.


title: Human-Centered Explainable AI (XAI): From Algorithms to User Experiences

abstract: In recent years, the field of explainable AI (XAI) has produced a vast collection of algorithms, providing a useful toolbox for researchers and practitioners to build XAI applications. With the rich application opportunities, explainability is believed to have moved beyond a demand by data scientists or researchers to comprehend the models they develop, to an essential requirement for people to trust and adopt AI deployed in numerous domains. However, explainability is an inherently human-centric property and the field is starting to embrace human-centered approaches. Human-computer interaction (HCI) research and user experience (UX) design in this area are becoming increasingly important. In this chapter, we begin with a high-level overview of the technical landscape of XAI algorithms, then selectively survey our own and other recent HCI works that take human-centered approaches to design, evaluate, and provide conceptual and methodological tools for XAI. We ask the question "what are human-centered approaches doing for XAI" and highlight three roles that they play in shaping XAI technologies by helping navigate, assess and expand the XAI toolbox: to drive technical choices by users' explainability needs, to uncover pitfalls of existing XAI methods and inform new methods, and to provide conceptual frameworks for human-compatible XAI.


title: Modeling Human-Human Collaboration: A Connection Between Inter-Personal Motor Synergy and Consensus Algorithms

abstract: Many day-to-day activities involve people working collaboratively toward reaching a desired outcome. Previous research in motor control and neuroscience have proposed inter-personal motor synergy (IPMS) as a mechanism of collaboration between people, referring to the idea of how two or more people may work together "as if they were one" to coordinate their motion. In motor control literature, uncontrolled manifold (UCM) is used for quantifying IPMS. According to this approach, coordinated motion is achieved through stabilization of a performance variable (e.g., an output in a collaborative output tracking task). We show that the UCM approach is closely related to the well-studied consensus approach in multi-agent systems that concerns processes by which a set of interacting agents agree on a shared objective. To explore the connection between these two approaches, in this paper, we provide a control-theoretic model that represents the systems-level behaviors in a collaborative task. In particular, we utilize the consensus protocol and show how the model can be systematically tuned to reproduce the behavior exhibited by human-human collaboration experiments. We discuss the association between the proposed control law and the UCM approach and validate our model using experimental results previously collected from an inter-personal finger force production task.


title: A Framework for Measures of Risk under Uncertainty

abstract: A risk analyst assesses potential financial losses based on multiple sources of information. Often, the assessment does not only depend on the specification of the loss random variable but also various economic scenarios. Motivated by this observation, we design a unified axiomatic framework for risk evaluation principles which quantifies jointly a loss random variable and a set of plausible probabilities. We call such an evaluation principle a generalized risk measure. We present a series of relevant theoretical results. The worst-case, coherent, and robust generalized risk measures are characterized via different sets of intuitive axioms. We establish the equivalence between a few natural forms of law invariance in our framework, and the technical subtlety therein reveals a sharp contrast between our framework and the traditional one. Moreover, coherence and strong law invariance are derived from a combination of other conditions, which provides additional support for coherent risk measures such as Expected Shortfall over Value-at-Risk, a relevant issue for risk management practice.


title: L\'evy-noise versus Gaussian-noise-induced Transitions in the Ghil-Sellers Energy Balance Model

abstract: We study the impact of applying stochastic forcing to the Ghil-Sellers energy balance climate model in the form of a fluctuating solar irradiance. Through numerical simulations, we explore the noise-induced transitions between the competing warm and snowball climate states. We consider multiplicative stochastic forcing driven by Gaussian and $\alpha$-stable L\'evy - $\alpha\in(0,2)$ - noise laws, examine the statistics of transition times, and estimate most probable transition paths. While the Gaussian noise case has been carefully studied in a plethora of investigations on metastable systems, much less is known about the L\'evy case, especially in the case of high- and infinite-dimensional systems. In the weak noise limit, the expected residence time in each metastable state scales in a fundamentally different way in the Gaussian vs. L\'evy noise case with respect to the intensity of the noise. In the former case, the classical Kramers-like exponential law is recovered. In the latter case, power laws are found, with the exponent equal to $-\alpha$, in apparent agreement with rigorous results obtained for additive noise in a related - yet different - reaction-diffusion equation as well as in simpler models. This can be better understood by treating the L\'evy noise as a compound Poisson process. The transition paths are studied in a projection of the state space and remarkable differences are observed between the two different types of noise. The snowball-to-warm and the warm-to-snowball most probable transition path cross at the single unstable edge state on the basin boundary. In the case of L\'evy noise, the most probable transition paths in the two directions are wholly separated, as transitions apparently take place via the closest basin boundary region to the outgoing attractor. This property can be better elucidated by considering singular perturbations to the solar irradiance.

categories: cond-mat.stat-mech physics.ao-ph physics.data-an physics.geo-ph


title: Low-overhead fault-tolerant quantum computing using long-range connectivity

abstract: Vast numbers of qubits will be needed for large-scale quantum computing due to the overheads associated with error correction. We present a scheme for low-overhead fault-tolerant quantum computation based on quantum low-density parity-check (LDPC) codes, where long-range interactions enable many logical qubits to be encoded with a modest number of physical qubits. In our approach, logic gates operate via logical Pauli measurements that preserve both the protection of the LDPC codes as well as the low overheads in terms of the required number of additional qubits. Compared with surface codes with the same code distance, we estimate order-of-magnitude improvements in the overheads for processing around one hundred logical qubits using this approach. Given the high thresholds demonstrated by LDPC codes, our estimates suggest that fault-tolerant quantum computation at this scale may be achievable with a few thousand physical qubits at comparable error rates to what is needed for current approaches.


title: Feasibility Study For Hydrogen Producing Colony on Mars

abstract: A technologically mature colony on Mars can produce and deliver at least 1 million tons of liquid hydrogen per year to one or more propellant depots at Low Earth Orbit (LEO). Production of 1 $kg$ of hydrogen at Marian colony and its delivery to LEO requires an energy expenditure of 1.4 $GJ$ on Mars. LEO propellant depot contains hydrogen produced on Mars and oxygen produced on the Moon or near-Earth asteroids. This propellant is used to deliver payload from LEO to many destinations in the Solar System including Mars. Delivery of 1 $kg$ payload from LEO to Mars requires an energy expenditure of 3.5 $GJ$ on Mars, Moon, and near-Earth asteroids. The use of liquid hydrogen produced on Mars to deliver astronauts and payload to Mars ensures exponential bootstrap growth of the Martian colony. Martian Colony and delivery of millions of tons of liquid hydrogen to LEO is the key to Colonization of Solar System. %% Martian Colony starts transporting liquid hydrogen to LEO only after it grows to significant size. It should contain about 20 million tons of steel and 3 million tons of plastic in structures and material as well as several thousand astronauts. Prior to that time, LEO hydrogen deposit will be supplied by hydrogen from Lunar poles.


title: Distribution and Moments of the Error Term in the Lattice Point Counting Problem for 3-Dimensional Cygan-Kor\'anyi Balls

abstract: We study fluctuations of the error term for the number of integer lattice points lying inside a 3-dimensional Cygan-Kor\'anyi ball of large radius. We prove that the error term, suitably normalized, has a limiting value distribution which is absolutely continuous, and we provide estimates for the decay rate of the corresponding density on the real line. In addition, we establish the existence of all moments for the normalized error term, and we prove that these are given by the moments of the corresponding density.


title: Scheduling of Graph Queries: Controlling Intra- and Inter-query Parallelism for a High System Throughput

abstract: The vast amounts of data used in social, business or traffic networks, biology and other natural sciences are often managed in graph-based data sets, consisting of a few thousand up to billions and trillions of vertices and edges, respectively. Typical applications utilizing such data either execute one or a few complex queries or many small queries at the same time interactively or as batch jobs. Furthermore, graph processing is inherently complex, as data sets can substantially differ (scale free vs. constant degree), and algorithms exhibit diverse behavior (computational intensity, local or global, push- or pull-based). This work is concerned with multi-query execution by automatically controlling the degree of parallelization, with overall objectives including high system utilization, low synchronization cost, and highly efficient concurrent execution. The underlying concept is three-fold: (1) sampling is used to determine graph statistics, (2) parallelization constraints are derived from algorithm and system properties, and (3) suitable work packages are generated based on the previous two aspects. We evaluate the proposed concept using different algorithms on synthetic and real world data sets, with up to 16 concurrent sessions (queries). The results demonstrate a robust performance in spite of these various configurations, and in particular that the performance is always close to or even slightly ahead of the performance of manually optimized implementations. Furthermore, the similar performance to manually optimized implementations under extreme configurations, which require either a full parallelization (few large queries) or complete sequential execution (many small queries), shows that the proposed concept exhibits a particularly low overhead.


title: Spin Fluctuations Yield zT Enhancement in Ferromagnets

abstract: Thermal fluctuation of local magnetization in magnetic metals intercoupled with charge carriers and phonons offers a path to enhance thermoelectric performance. The thermopower enhancement by spin fluctuations (SF) has been observed before. However, the crucial evidence for enhancing thermoelectric-figure-of-merit (zT) by SF has not been reported until now. Here we report evidence for such enhancement in the ferromagnetic CrTe. The SF leads to nearly 80% zT enhancement in ferromagnetic CrTe near and below TC~335 K. The ferromagnetism in CrTe is originated from the collective electronic and localized magnetic moments. The field-dependent transport properties demonstrate the profound impact of SF on the electrons and phonons. The SF simultaneously enhances the thermopower and reduces the thermal conductivity. Under an external magnetic field, the enhancement in thermopower is suppressed, and the thermal conductivity is enhanced, evidencing the existence of a strong SF near and below TC. The anomalous thermoelectric transport properties are analyzed based on theoretical models, and a good agreement with experimental data is found. Furthermore, the detailed analysis proves an insignificant impact from spin-wave contribution to the transport properties. This study contributes to the fundamental understanding of spin fluctuation for designing high-performance spin-driven thermoelectric materials.


title: Band Structure and Polarization Effects in Photothermoelectric Spectroscopy of a Bi2Se3 Device

abstract: Bi2Se3 is a prototypical topological insulator which has a small band gap (~0.3 eV) and topologically protected conducting surface states. This material exhibits quite strong thermoelectric effects. Here we show in a mechanically exfoliated thick (~100 nm) nanoflake device that we can measure the energy dependent optical absorption through the photothermoelectric effect. Spectral signatures are seen for a number of optical transitions between the valence and conduction bands, including a broad peak at 1.5 eV which is likely dominated by bulk band-to-band optical transitions but is at the same energy as the well known optical transition between the two topologically protected conducting surface states. We also observe a surprising linear polarization dependence in the response of the device that reflects the influence of the metal contacts.


title: Media abnormal tone, earnings announcements, and the stock market

abstract: We conduct a tone-based event study to examine the aggregate abnormal tone dynamics in media articles around earnings announcements. We test whether they convey incremental information that is useful for price discovery for nonfinancial S&P 500 firms. The relation we find between the abnormal tone and abnormal returns suggests that media articles provide incremental information relative to the information contained in earnings press releases and earnings calls.


title: Simulating Ising and Potts models at critical and cold temperatures using auxiliary Gaussian variables

abstract: Ising and Potts models are an important class of discrete probability distributions which originated from Statistical Physics and since then have found applications in several disciplines. Simulation from these models is a well known challenging problem. In this paper, we propose a class of MCMC algorithms to simulate from both Ising and Potts models, by using auxiliary Gaussian random variables. Our algorithms apply to coupling matrices with both positive and negative entries, thus including Spin Glass models such as the SK and Hopfield model. In contrast to existing methods of a similar flavor, our algorithm can take advantage of the low-rank structure of the coupling matrix, and scales linearly with the number of states in a Potts model. We compare our proposed algorithm to existing state of the art algorithms, such as the Swendsen-Wang and Wolff algorithms for Ising and Potts models on graphs, and the Heat Bath for Spin Glass models. Our comparison takes into account a wide range of coupling matrices and temperature regimes, focusing in particular on behavior at or below the critical temperature. For cold systems, augmenting our algorithm with a tempering scheme yields significant improvements.


title: A Data-Centric Optimization Framework for Machine Learning

abstract: Rapid progress in deep learning is leading to a diverse set of quickly changing models, with a dramatically growing demand for compute. However, as frameworks specialize performance optimization to patterns in popular networks, they implicitly constrain novel and diverse models that drive progress in research. We empower deep learning researchers by defining a flexible and user-customizable pipeline for optimizing training of arbitrary deep neural networks, based on data movement minimization. The pipeline begins with standard networks in PyTorch or ONNX and transforms computation through progressive lowering. We define four levels of general-purpose transformations, from local intra-operator optimizations to global data movement reduction. These operate on a data-centric graph intermediate representation that expresses computation and data movement at all levels of abstraction, including expanding basic operators such as convolutions to their underlying computations. Central to the design is the interactive and introspectable nature of the pipeline. Every part is extensible through a Python API, and can be tuned interactively using a GUI. We demonstrate competitive performance or speedups on ten different networks, with interactive optimizations discovering new opportunities in EfficientNet.


title: Propensity-scored Probabilistic Label Trees

abstract: Extreme multi-label classification (XMLC) refers to the task of tagging instances with small subsets of relevant labels coming from an extremely large set of all possible labels. Recently, XMLC has been widely applied to diverse web applications such as automatic content labeling, online advertising, or recommendation systems. In such environments, label distribution is often highly imbalanced, consisting mostly of very rare tail labels, and relevant labels can be missing. As a remedy to these problems, the propensity model has been introduced and applied within several XMLC algorithms. In this work, we focus on the problem of optimal predictions under this model for probabilistic label trees, a popular approach for XMLC problems. We introduce an inference procedure, based on the $A^*$-search algorithm, that efficiently finds the optimal solution, assuming that all probabilities and propensities are known. We demonstrate the attractiveness of this approach in a wide empirical study on popular XMLC benchmark datasets.


title: Identifiable Deep Generative Models via Sparse Decoding

abstract: We develop the sparse VAE for unsupervised representation learning on high-dimensional data. The sparse VAE learns a set of latent factors (representations) which summarize the associations in the observed data features. The underlying model is sparse in that each observed feature (i.e. each dimension of the data) depends on a small subset of the latent factors. As examples, in ratings data each movie is only described by a few genres; in text data each word is only applicable to a few topics; in genomics, each gene is active in only a few biological processes. We prove such sparse deep generative models are identifiable: with infinite data, the true model parameters can be learned. (In contrast, most deep generative models are not identifiable.) We empirically study the sparse VAE with both simulated and real data. We find that it recovers meaningful latent factors and has smaller heldout reconstruction error than related methods.


title: DVIO: Depth aided visual inertial odometry for RGBD sensors

abstract: In past few years we have observed an increase in the usage of RGBD sensors in mobile devices. These sensors provide a good estimate of the depth map for the camera frame, which can be used in numerous augmented reality applications. This paper presents a new visual inertial odometry (VIO) system, which uses measurements from a RGBD sensor and an inertial measurement unit (IMU) sensor for estimating the motion state of the mobile device. The resulting system is called the depth-aided VIO (DVIO) system. In this system we add the depth measurement as part of the nonlinear optimization process. Specifically, we propose methods to use the depth measurement using one-dimensional (1D) feature parameterization as well as three-dimensional (3D) feature parameterization. In addition, we propose to utilize the depth measurement for estimating time offset between the unsynchronized IMU and the RGBD sensors. Last but not least, we propose a novel block-based marginalization approach to speed up the marginalization processes and maintain the real-time performance of the overall system. Experimental results validate that the proposed DVIO system outperforms the other state-of-the-art VIO systems in terms of trajectory accuracy as well as processing time.


title: How are Lyman $\alpha$ absorbers in the cosmic web related to gas-rich galaxies?

abstract: We present the two-point cross-correlation function between Lyman $\alpha$ absorbers and HI galaxies in the nearby Universe ($0.01 \le z \le 0.057$). We use absorbers from 21 QSO sightlines from the Survey of the Low-Redshift Intergalactic Medium with HST/COS and the galaxy catalogs from the Arecibo Legacy Fast ALFA survey and the New York University Value-Added Galaxy Catalog. We find that Lyman $\alpha$ absorbers are strongly correlated to Lyman $\alpha$ galaxies at a projected separation of $\le$0.5 Mpc and velocity separation of $\le 50~km~s^{-1}$. Lyman $\alpha$ absorbers are 7.6 times more likely to be found near HI galaxies compared to a random distribution. The correlation decreases as the projected and/or velocity separation increase. We also find the correlation between Lyman $\alpha$ absorbers and HI galaxies to be stronger than those observed between Lyman $\alpha$ absorbers and optically selected galaxies. There is an enhancement in the number of absorbers at velocity separations of $\le30~km~s^{-1}$ from HI galaxies at distances larger than their viral radius. Combined with the fact that most of our galaxies are not driving strong outflows, we conclude that the absorbers at low-velocity separations are tracing cooler intergalactic gas around galaxies. This conclusion is consistent with the predictions from cosmological simulations where faint gas from the intergalactic medium flows into the disks of galaxies leading to galaxy growth.


title: Text-Based Person Search with Limited Data

abstract: Text-based person search (TBPS) aims at retrieving a target person from an image gallery with a descriptive text query. Solving such a fine-grained cross-modal retrieval task is challenging, which is further hampered by the lack of large-scale datasets. In this paper, we present a framework with two novel components to handle the problems brought by limited data. Firstly, to fully utilize the existing small-scale benchmarking datasets for more discriminative feature learning, we introduce a cross-modal momentum contrastive learning framework to enrich the training data for a given mini-batch. Secondly, we propose to transfer knowledge learned from existing coarse-grained large-scale datasets containing image-text pairs from drastically different problem domains to compensate for the lack of TBPS training data. A transfer learning method is designed so that useful information can be transferred despite the large domain gap. Armed with these components, our method achieves new state of the art on the CUHK-PEDES dataset with significant improvements over the prior art in terms of Rank-1 and mAP. Our code is available at https://github.com/BrandonHanx/TextReID.


title: The Arcanum Mission: Scientific Objectives and Instruments for Neptune, Triton and KBOs

abstract: The Arcanum mission is a proposed L-class spacecraft that highlights the revolutionary approach which can now be taken to future space mission design. Using the case of the SpaceX Starship vehicle and in particular the high mass and volume characteristics of this launcher, the feasible large size of future missions, even with high delta-V transfer requirements, are analysed. A demonstrator vehicle, designed to support a large and capable science platform with multiple components, is detailed, clearly showing the range and depth of science goals that will be answerable thanks to the current revolution in super heavy-lift launch vehicles.


title: Hierarchical Skills for Efficient Exploration

abstract: In reinforcement learning, pre-trained low-level skills have the potential to greatly facilitate exploration. However, prior knowledge of the downstream task is required to strike the right balance between generality (fine-grained control) and specificity (faster learning) in skill design. In previous work on continuous control, the sensitivity of methods to this trade-off has not been addressed explicitly, as locomotion provides a suitable prior for navigation tasks, which have been of foremost interest. In this work, we analyze this trade-off for low-level policy pre-training with a new benchmark suite of diverse, sparse-reward tasks for bipedal robots. We alleviate the need for prior knowledge by proposing a hierarchical skill learning framework that acquires skills of varying complexity in an unsupervised manner. For utilization on downstream tasks, we present a three-layered hierarchical learning algorithm to automatically trade off between general and specific skills as required by the respective task. In our experiments, we show that our approach performs this trade-off effectively and achieves better results than current state-of-the-art methods for end- to-end hierarchical reinforcement learning and unsupervised skill discovery. Code and videos are available at https://facebookresearch.github.io/hsd3 .


title: Stability of adsorption of Mg and Na on sulfur-functionalized MXenes

abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate adsorption of Mg and Na adatoms on five M$_{2}$CS$_{2}$ monolayers (where M= Mo, Nb, Ti, V, Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo$_{2}$CS$_{2}$ cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr$_{2}$CS$_{2}$, followed by Ti$_{2}$CS$_{2}$, Nb$_{2}$CS$_{2}$ and V$_{2}$CS$_{2}$. Using the nudged elastic band method (NEB) we computed promising diffusion barriers for both dilute and nearly-full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti$_{2}$CS$_{2}$ experience $\sim$0.47 eV and $\sim$0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti$_{2}$CS$_{2}$ experiences a $\sim$0.33 eV energy barrier, implying a concentration dependent diffusion barrier. Our molecular dynamics results indicate that three (one) layers (layer) of Mg (Na) ion on both surfaces of Ti$_{2}$CS$_{2}$ remain stable at T=300 K. While, according to voltage calculations, Zr$_{2}$CS$_{2}$ can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from Zr$_{2}$CS$_{2}$ monolayer due to weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirming the stability of an ion-electrode system - an insight that is mostly absence in previous studies.


title: HALP: Hardware-Aware Latency Pruning

abstract: Structural pruning can simplify network architecture and improve inference speed. We propose Hardware-Aware Latency Pruning (HALP) that formulates structural pruning as a global resource allocation optimization problem, aiming at maximizing the accuracy while constraining latency under a predefined budget. For filter importance ranking, HALP leverages latency lookup table to track latency reduction potential and global saliency score to gauge accuracy drop. Both metrics can be evaluated very efficiently during pruning, allowing us to reformulate global structural pruning under a reward maximization problem given target constraint. This makes the problem solvable via our augmented knapsack solver, enabling HALP to surpass prior work in pruning efficacy and accuracy-efficiency trade-off. We examine HALP on both classification and detection tasks, over varying networks, on ImageNet and VOC datasets. In particular, for ResNet-50/-101 pruning on ImageNet, HALP improves network throughput by $1.60\times$/$1.90\times$ with $+0.3\%$/$-0.2\%$ top-1 accuracy changes, respectively. For SSD pruning on VOC, HALP improves throughput by $1.94\times$ with only a $0.56$ mAP drop. HALP consistently outperforms prior art, sometimes by large margins.


title: REAL-M: Towards Speech Separation on Real Mixtures

abstract: In recent years, deep learning based source separation has achieved impressive results. Most studies, however, still evaluate separation models on synthetic datasets, while the performance of state-of-the-art techniques on in-the-wild speech data remains an open question. This paper contributes to fill this gap in two ways. First, we release the REAL-M dataset, a crowd-sourced corpus of real-life mixtures. Secondly, we address the problem of performance evaluation of real-life mixtures, where the ground truth is not available. We bypass this issue by carefully designing a blind Scale-Invariant Signal-to-Noise Ratio (SI-SNR) neural estimator. Through a user study, we show that our estimator reliably evaluates the separation performance on real mixtures. The performance predictions of the SI-SNR estimator indeed correlate well with human opinions. Moreover, we observe that the performance trends predicted by our estimator on the REAL-M dataset closely follow those achieved on synthetic benchmarks when evaluating popular speech separation models.


title: CXR-Net: An Encoder-Decoder-Encoder Multitask Deep Neural Network for Explainable and Accurate Diagnosis of COVID-19 pneumonia with Chest X-ray Images

abstract: Accurate and rapid detection of COVID-19 pneumonia is crucial for optimal patient treatment. Chest X-Ray (CXR) is the first line imaging test for COVID-19 pneumonia diagnosis as it is fast, cheap and easily accessible. Inspired by the success of deep learning (DL) in computer vision, many DL-models have been proposed to detect COVID-19 pneumonia using CXR images. Unfortunately, these deep classifiers lack the transparency in interpreting findings, which may limit their applications in clinical practice. The existing commonly used visual explanation methods are either too noisy or imprecise, with low resolution, and hence are unsuitable for diagnostic purposes. In this work, we propose a novel explainable deep learning framework (CXRNet) for accurate COVID-19 pneumonia detection with an enhanced pixel-level visual explanation from CXR images. The proposed framework is based on a new Encoder-Decoder-Encoder multitask architecture, allowing for both disease classification and visual explanation. The method has been evaluated on real world CXR datasets from both public and private data sources, including: healthy, bacterial pneumonia, viral pneumonia and COVID-19 pneumonia cases The experimental results demonstrate that the proposed method can achieve a satisfactory level of accuracy and provide fine-resolution classification activation maps for visual explanation in lung disease detection. The Average Accuracy, the Precision, Recall and F1-score of COVID-19 pneumonia reached 0.879, 0.985, 0.992 and 0.989, respectively. We have also found that using lung segmented (CXR) images can help improve the performance of the model. The proposed method can provide more detailed high resolution visual explanation for the classification decision, compared to current state-of-the-art visual explanation methods and has a great potential to be used in clinical practice for COVID-19 pneumonia diagnosis.


title: Symmetry for algebras associated to Fell bundles over groups and groupoids

abstract: To every Fell bundle $\mathscr C$ over a locally compact group ${\sf G}$ one associates a Banach $^*$-algebra $L^1({\sf G}\,\vert\,\mathscr C)$. We prove that it is symmetric whenever ${\sf G}$ with the discrete topology is rigidly symmetric. This generalizes the known case of a global action without a twist. There is also a weighted version as well as a treatment of some classes of associated integral kernels. We also deal with the case of Fell bundles over discrete groupoids. We formulate a generalization of rigid symmetry in this case and show its equivalence with an a priori stronger concept. We also study the symmetry of transformation groupoids and some permanence properties.


title: Convergence Analysis and Implicit Regularization of Feedback Alignment for Deep Linear Networks

abstract: We theoretically analyze the Feedback Alignment (FA) algorithm, an efficient alternative to backpropagation for training neural networks. We provide convergence guarantees with rates for deep linear networks for both continuous and discrete dynamics. Additionally, we study incremental learning phenomena for shallow linear networks. Interestingly, certain specific initializations imply that negligible components are learned before the principal ones, thus potentially negatively affecting the effectiveness of such a learning algorithm; a phenomenon we classify as implicit anti-regularization. We also provide initialization schemes where the components of the problem are approximately learned by decreasing order of importance, thus providing a form of implicit regularization.


title: Problems of Position Reconstruction in Silicon Microstrip Detects

abstract: The algorithms for position reconstruction in silicon micro-strip detectors are studied, and the signals of a minimum ionizing particle are simulated. The center-of-gravity distributions of the data events allow the fine tuning of the signal forms and of the strip response functions. In the sensors with floating strips the response function turns out to roughly approximate the response of a triangular function for the size of the signal involved. The simulations are extended to non orthogonal incidence. In these directions, and in general for all the asymmetric signal distributions, the standard application of the $\eta$ algorithm introduces a systematic error. A signal reconstruction theorem gives the way to implement the corrections of this error.


title: The R package sentometrics to compute, aggregate and predict with textual sentiment

abstract: We provide a hands-on introduction to optimized textual sentiment indexation using the R package sentometrics. Textual sentiment analysis is increasingly used to unlock the potential information value of textual data. The sentometrics package implements an intuitive framework to efficiently compute sentiment scores of numerous texts, to aggregate the scores into multiple time series, and to use these time series to predict other variables. The workflow of the package is illustrated with a built-in corpus of news articles from two major U.S. journals to forecast the CBOE Volatility Index.


title: Singularities of 3-parameter line congruences in $\mathbb{R}^4$

abstract: In this paper, we give the generic classification of the singularities of 3-parameter line congruences in $\mathbb{R}^4$. We also classify the generic singularities of Blaschke (affine) normal congruences.


title: Shaking the foundations: delusions in sequence models for interaction and control

abstract: The recent phenomenal success of language models has reinvigorated machine learning research, and large sequence models such as transformers are being applied to a variety of domains. One important problem class that has remained relatively elusive however is purposeful adaptive behavior. Currently there is a common perception that sequence models "lack the understanding of the cause and effect of their actions" leading them to draw incorrect inferences due to auto-suggestive delusions. In this report we explain where this mismatch originates, and show that it can be resolved by treating actions as causal interventions. Finally, we show that in supervised learning, one can teach a system to condition or intervene on data by training with factual and counterfactual error signals respectively.


title: Rigid inner forms over global function fields

abstract: We construct an fpqc gerbe $\mathcal{E}_{\dot{V}}$ over a global function field $F$ such that for a connected reductive group $G$ over $F$ with finite central subgroup $Z$, the set of $G_{\mathcal{E}_{\dot{V}}}$-torsors contains a subset $H^{1}(\mathcal{E}_{\dot{V}}, Z \to G)$ which allows one to define a global notion of ($Z$-)rigid inner forms. There is a localization map $H^{1}(\mathcal{E}_{\dot{V}}, Z \to G) \to H^{1}(\mathcal{E}_{v}, Z \to G)$, where the latter parametrizes local rigid inner forms (cf. [Kal16, Dil23]) which allows us organize local rigid inner forms across all places $v$ into coherent families. Doing so enables a construction of (conjectural) global $L$-packets and a conjectural formula for the multiplicity of an automorphic representation $\pi$ in the discrete spectrum of $G$ in terms of these $L$-packets. We also show that, for a connected reductive group $G$ over a global function field $F$, the adelic transfer factor $\Delta_{\mathbb{A}}$ for the ring of adeles $\mathbb{A}$ of $F$ serving an endoscopic datum for $G$ decomposes as the product of the normalized local transfer factors from [Dil20].


title: On the feasibility of dynamical analysis of network models of biochemical regulation

abstract: A recent article by Weidner et al. [2021] presents a method to extract graph properties that are predictive of the dynamical behavior of multivariate, discrete models of biochemical regulation. In other words, a method that uses only features from the structure of network interactions to predict which nodes are most involved in automata network dynamics. However, the authors claim that dynamical analysis of large automata network models is "not even feasible." To make sure that others are not discouraged from working on this problem, it is important to clarify that effective dynamical analysis of automata network models, to the contrary, is feasible. Unlike what is suggested in the article, graph-based analysis of static features is not the only analytical avenue for large systems biology models of regulation and signaling dynamics because there are dynamical methods that are, indeed, scalable.


title: Gauge Field Production and Schwinger Reheating in Runaway Axion Inflation

abstract: In a class of (pseudoscalar) inflation, inflationary phase is followed by a kination phase, where the Universe is dominated by the kinetic energy of the inflaton that runs away in a vanishing scalar potential. In this class of postinflationary evolution of the Universe, reheating of the Universe cannot be achieved by the inflaton particle decay, which requires its coherent oscillation in a quadratic potential. In this study, we explore the U(1) gauge field production through the Chern-Simons coupling between the pseudoscalar inflaton and the gauge field during the kination era and examine the subsequent pair-particle production induced by the amplified gauge field known as the Schwinger effect, which can lead to reheating of the Universe. We find that with a rough estimate of the Schwinger effect for the Standard Model hyper U(1) gauge field and subsequent thermalization of the pair-produced particles, a successful reheating of the Universe can be achieved by their eventual domination over the kinetic energy of the inflaton, with some reasonable parameter sets. This can be understood as a concrete realization of the "Schwinger reheating". Constraints from the later-time cosmology are also discussed.


title: Electronic confinement in quantum dots of twisted bilayer graphene

abstract: Electronic properties of quantum dots (QDs) depend sensitively on their parent materials. Therefore, confined electronic states in graphene QDs (GQDs) of monolayer and Bernal-stacked bilayer graphene are quite different. Twisted bilayer graphene (TBG) is distinct from monolayer and Bernal-stacked bilayer graphene because of the new degree of freedom: twist angle. In the past few years, numerous efforts have been made to realize the GQDs of monolayer and Bernal-stacked bilayer graphene and achieved great success. Thus far, however, strategies for realizing GQDs of TBG have been elusive. Here, we demonstrate a general approach for fabricating stationary GQDs of TBG by introducing nanoscale p-n junctions with sharp boundaries in the TBG. We verify the confinement of low-energy massless Dirac fermions via whispering-gallery modes in the GQDs of TBG. Unexpectedly, electronic states around van Hove singularities of the TBG are also strongly modified around the GQDs. Such a feature has never been reported and is attributed to spatial variation of the interlayer coupling in the TBG induced by the GQDs.


title: Dynamic Bipartite Matching Market with Arrivals and Departures

abstract: In this paper, we study a matching market model on a bipartite network where agents on each side arrive and depart stochastically by a Poisson process. For such a dynamic model, we design a mechanism that decides not only which agents to match, but also when to match them, to minimize the expected number of unmatched agents. The main contribution of this paper is to achieve theoretical bounds on the performance of local mechanisms with different timing properties. We show that an algorithm that waits to thicken the market, called the $\textit{Patient}$ algorithm, is exponentially better than the $\textit{Greedy}$ algorithm, i.e., an algorithm that matches agents greedily. This means that waiting has substantial benefits on maximizing a matching over a bipartite network. We remark that the Patient algorithm requires the planner to identify agents who are about to leave the market, and, under the requirement, the Patient algorithm is shown to be an optimal algorithm. We also show that, without the requirement, the Greedy algorithm is almost optimal. In addition, we consider the $\textit{1-sided algorithms}$ where only an agent on one side can attempt to match. This models a practical matching market such as a freight exchange market and a labor market where only agents on one side can make a decision. For this setting, we prove that the Greedy and Patient algorithms admit the same performance, that is, waiting to thicken the market is not valuable. This conclusion is in contrast to the case where agents on both sides can make a decision and the non-bipartite case by [Akbarpour et al.,$~\textit{Journal of Political Economy}$, 2020].


title: $\ell_{\infty}$-Bounds of the MLE in the BTL Model under General Comparison Graphs

abstract: The Bradley-Terry-Luce (BTL) model is a popular statistical approach for estimating the global ranking of a collection of items using pairwise comparisons. To ensure accurate ranking, it is essential to obtain precise estimates of the model parameters in the $\ell_{\infty}$-loss. The difficulty of this task depends crucially on the topology of the pairwise comparison graph over the given items. However, beyond very few well-studied cases, such as the complete and Erd\"os-R\'enyi comparison graphs, little is known about the performance of the maximum likelihood estimator MLE) of the BTL model parameters in the $\ell_{\infty}$-loss under more general graph topologies. In this paper, we derive novel, general upper bounds on the $\ell_{\infty}$ estimation error of the BTL MLE that depend explicitly on the algebraic connectivity of the comparison graph, the maximal performance gap across items and the sample complexity. We demonstrate that the derived bounds perform well and in some cases are sharper compared to known results obtained using different loss functions and more restricted assumptions and graph topologies. We carefully compare our results to Yan et al. (2012), which is closest in spirit to our work. We further provide minimax lower bounds under $\ell_{\infty}$-error that nearly match the upper bounds over a class of sufficiently regular graph topologies. Finally, we study the implications of our $\ell_{\infty}$-bounds for efficient (offline) tournament design. We illustrate and discuss our findings through various examples and simulations.


title: Ideal Mixing of Paraelectric and Ferroelectric Nematic Phases in Liquid Crystals of Distinct Molecular Species

abstract: The organic mesogens RM734 and DIO are members of separate molecular families featuring distinct molecular structures. These families, at the present time, are the only ones known to exhibit a ferroelectric nematic liquid crystal (LC) phase. Here we present an experimental investigation of the phase diagram and electro-optics of binary mixtures of RM734 and DIO. We observe paraelectric nematic (N) and ferroelectric nematic (NF) phases in both materials, each of which exhibits complete miscibility across the phase diagram, showing that the paraelectric and ferroelectric are the same phases in RM734 as in DIO. Remarkably, these molecules form ideal mixtures with respect to both the paraelectric-ferroelectric nematic phase behavior and the ferroelectric polarization density of the mixtures, the principal order parameter of the transition. Ideal mixing is also manifested in the orientational viscosity, and the onset of glassy dynamics at low temperature. This behavior is attributable in part to the similarity of their overall molecular shape and net longitudinal dipole moment (~ 11 Debye), and to a common tendency for head-to-tail molecular association. In contrast, the significant difference in molecular structures leads to poor solubility in the crystal phases, enhancing the stability of the ferroelectric nematic phase at low temperature in the mixtures and making possible room temperature electro-optic effects. In the mixtures with excess DIO, an intermediate phase appears via an ultraweak, first-order transition from the N phase, in a narrow temperature range between the paraelectric and ferroelectric nematics.


title: How to pose material design problems for flow through porous media applications?: Sensitivity of dissipation rate to medium's permeability holds the key

abstract: Recent studies have advocated using the total dissipation rate under topology optimization to realize material designs involving the flow of fluids through porous media. However, these studies decided how to pose the design problem, such as maximizing the total dissipation rate for some situations while minimizing for others, by solving one-dimensional problems and justifying their choices using numerical experiments. The rigor is lacking -- a bottleneck for further scientific advancements to computational material design. This paper provides the missing theoretical justification. We identify four classes of boundary value problems using the adjoint state method and analytically calculate the sensitivity of the total dissipation rate to the permeability field. For two of those classes in which the flow of fluids is pressure-driven, the sensitivity is positive -- the total dissipation rate increases if the medium's permeability increases. While for the other two classes, in which the flow is velocity-driven, the trend is the opposite. These sensitivities provide rigorous answers to the central question: how to pose a material design problem for flow through porous media applications. The impact of our work is multi-fold. First, this study further elevates the role of the dissipation rate in posing well-posed material design problems using topology optimization. Second, besides the theoretical significance, the results benefit computational scientists and practitioners to realize optimal designs. Third, given their simplicity yet far-reaching impact, both the approach and results possess immense pedagogical value.


title: AdamD: Improved bias-correction in Adam

abstract: Here I present a small update to the bias-correction term in the Adam optimizer that has the advantage of making smaller gradient updates in the first several steps of training. With the default bias-correction, Adam may actually make larger than requested gradient updates early in training. By only including the well-justified bias-correction of the second moment gradient estimate, $v_t$, and excluding the bias-correction on the first-order estimate, $m_t$, we attain these more desirable gradient update properties in the first series of steps. The default implementation of Adam may be as sensitive as it is to the hyperparameters $\beta_1, \beta_2$ partially due to the originally proposed bias correction procedure, and its behavior in early steps.


title: ReachBot: A Small Robot for Large Mobile Manipulation Tasks

abstract: Robots are widely deployed in space environments because of their versatility and robustness. However, adverse gravity conditions and challenging terrain geometry expose the limitations of traditional robot designs, which are often forced to sacrifice one of mobility or manipulation capabilities to attain the other. Prospective climbing operations in these environments reveals a need for small, compact robots capable of versatile mobility and manipulation. We propose a novel robotic concept called ReachBot that fills this need by combining two existing technologies: extendable booms and mobile manipulation. ReachBot leverages the reach and tensile strength of extendable booms to achieve an outsized reachable workspace and wrench capability. Through their lightweight, compactable structure, these booms also reduce mass and complexity compared to traditional rigid-link articulated-arm designs. Using these advantages, ReachBot excels in mobile manipulation missions in low gravity or that require climbing, particularly when anchor points are sparse. After introducing the ReachBot concept, we discuss modeling approaches and strategies for increasing stability and robustness. We then develop a 2D analytical model for ReachBot's dynamics inspired by grasp models for dexterous manipulators. Next, we introduce a waypoint-tracking controller for a planar ReachBot in microgravity. Our simulation results demonstrate the controller's robustness to disturbances and modeling error. Finally, we briefly discuss next steps that build on these initially promising results to realize the full potential of ReachBot.


title: Bounds for moments of modular $L$-functions to a fixed modulus

abstract: We study the $2k$-th moment of the family of twisted modular $L$-functions to a fixed prime power modulus at the central values. We establish sharp lower bounds for all real $k \geq 0$ and sharp upper bounds for $k$ in the range $0 \leq k \leq 1$.


title: Investigating the Role of Structural Water on the Electrochemical Properties of $\alpha$-V$_2$O$_5$ through Density Functional Theory

abstract: The $\alpha$ polymorph of V$_2$O$_5$ is one of the few known cathodes capable of reversibly intercalating multivalent ions such as Mg, Ca, Zn and Al, but suffers from sluggish diffusion kinetics. The role of H$_2$O within the electrolyte and between the layers of the structure in the form of a xerogel/aerogel structure, though, has been shown to lower diffusion barriers and lead to other improved electrochemical properties. This density functional theory study systematically investigates how and why the presence of structural H$_2$O within $\alpha$-V$_2$O$_5$ changes the resulting structure, voltage, and diffusion kinetics for the intercalation of Li, Na, Mg, Ca, Zn, and Al. We found that the coordination of H$_2$O molecules with the ion leads to an improvement in voltage and energy density for all ions. This voltage increase was attributed to the extra host sites for electrons present with H$_2$O, thus leading to a stronger ionization of the ion and a higher voltage. We also found that the increase in interlayer distance and a potential "charge shielding" effect drastically changes the electrostatic environment and the resulting diffusion kinetics. For Mg and Ca, this resulted in a decrease in diffusion barrier from 1.3 eV and 2.0 eV to 0.89 eV and 0.4 eV, respectively. We hope that our study motivates similar research regarding the role of water in both V$_2$O$_5$ xerogels/aerogels and other layered transition metal oxides.


title: Ensemble of Averages: Improving Model Selection and Boosting Performance in Domain Generalization

abstract: In Domain Generalization (DG) settings, models trained independently on a given set of training domains have notoriously chaotic performance on distribution shifted test domains, and stochasticity in optimization (e.g. seed) plays a big role. This makes deep learning models unreliable in real world settings. We first show that this chaotic behavior exists even along the training optimization trajectory of a single model, and propose a simple model averaging protocol that both significantly boosts domain generalization and diminishes the impact of stochasticity by improving the rank correlation between the in-domain validation accuracy and out-domain test accuracy, which is crucial for reliable early stopping. Taking advantage of our observation, we show that instead of ensembling unaveraged models (that is typical in practice), ensembling moving average models (EoA) from independent runs further boosts performance. We theoretically explain the boost in performance of ensembling and model averaging by adapting the well known Bias-Variance trade-off to the domain generalization setting. On the DomainBed benchmark, when using a pre-trained ResNet-50, this ensemble of averages achieves an average of $68.0\%$, beating vanilla ERM (w/o averaging/ensembling) by $\sim 4\%$, and when using a pre-trained RegNetY-16GF, achieves an average of $76.6\%$, beating vanilla ERM by $6\%$. Our code is available at https://github.com/salesforce/ensemble-of-averages.


title: High-resolution rainfall-runoff modeling using graph neural network

abstract: Time-series modeling has shown great promise in recent studies using the latest deep learning algorithms such as LSTM (Long Short-Term Memory). These studies primarily focused on watershed-scale rainfall-runoff modeling or streamflow forecasting, but the majority of them only considered a single watershed as a unit. Although this simplification is very effective, it does not take into account spatial information, which could result in significant errors in large watersheds. Several studies investigated the use of GNN (Graph Neural Networks) for data integration by decomposing a large watershed into multiple sub-watersheds, but each sub-watershed is still treated as a whole, and the geoinformation contained within the watershed is not fully utilized. In this paper, we propose the GNRRM (Graph Neural Rainfall-Runoff Model), a novel deep learning model that makes full use of spatial information from high-resolution precipitation data, including flow direction and geographic information. When compared to baseline models, GNRRM has less over-fitting and significantly improves model performance. Our findings support the importance of hydrological data in deep learning-based rainfall-runoff modeling, and we encourage researchers to include more domain knowledge in their models.


title: Integrating Visuospatial, Linguistic and Commonsense Structure into Story Visualization

abstract: While much research has been done in text-to-image synthesis, little work has been done to explore the usage of linguistic structure of the input text. Such information is even more important for story visualization since its inputs have an explicit narrative structure that needs to be translated into an image sequence (or visual story). Prior work in this domain has shown that there is ample room for improvement in the generated image sequence in terms of visual quality, consistency and relevance. In this paper, we first explore the use of constituency parse trees using a Transformer-based recurrent architecture for encoding structured input. Second, we augment the structured input with commonsense information and study the impact of this external knowledge on the generation of visual story. Third, we also incorporate visual structure via bounding boxes and dense captioning to provide feedback about the characters/objects in generated images within a dual learning setup. We show that off-the-shelf dense-captioning models trained on Visual Genome can improve the spatial structure of images from a different target domain without needing fine-tuning. We train the model end-to-end using intra-story contrastive loss (between words and image sub-regions) and show significant improvements in several metrics (and human evaluation) for multiple datasets. Finally, we provide an analysis of the linguistic and visuo-spatial information. Code and data: https://github.com/adymaharana/VLCStoryGan.


title: Semiclassical limit problems with concurrent use of several clocks in quantum cosmology

abstract: We revisit a recent proposal for a definition of time in quantum cosmology, to investigate the effects of having more than one possible type of clock "at the same time". We use as test tube an extension of Einstein gravity with a massless scalar field in which the gravitational coupling $G_N$ is only a constant on-shell, mimicking the procedure for $\Lambda$ in unimodular gravity. Hence we have two "simultaneous" clocks in the theory: a scalar field clock, and the conjugate of $G_N$. We find that attempts to use two coherent clocks concurrently are disastrous for recovering the classical limit. The Heisenberg relations, instead of being saturated, are always realized abundantly above their bound, with strong quantum effects expected at least in parts of the trajectory. Semi-classical states always result from situations where we effectively impose a single clock, either by making the other clock a failed clock (i.e. by choosing a state where its conjugate constant is infinitely sharp), or by choosing a basis of constants where all clocks but one are redundant, i.e. motion or change in phase space does not occur with the passing of their "times". We show how this conclusion generalizes to fluids with any equation of state. It also applies to systems where "sub-clocks" of the same type could be used, for example in mixtures of different fluids with the same equation of state.


title: Application of E-Commerce Technologies in accelerating the Success of SME Operation

abstract: Electronic commerce (e-Commerce) technologies have been increased over the past two decades in different business sectors. In particular, the technologies of B2C operations have significantly improved the productivity of online small businesses such as SMEs. Systematic literature review in this domain categorized different benefits but a limited number of studies on SME success from a view of an information systems (IS) research exists, which needs to be taken for further attention. Informing through a comprehensive analysis this study introduces a conceptual framework of the application of e-Commerce technologies in accelerating the SME operation. Content analysis methodology was adopted for generating the outcome associated with the success of the technologies in SMEs.


title: A Domain Gap Aware Generative Adversarial Network for Multi-domain Image Translation

abstract: Recent image-to-image translation models have shown great success in mapping local textures between two domains. Existing approaches rely on a cycle-consistency constraint that supervises the generators to learn an inverse mapping. However, learning the inverse mapping introduces extra trainable parameters and it is unable to learn the inverse mapping for some domains. As a result, they are ineffective in the scenarios where (i) multiple visual image domains are involved; (ii) both structure and texture transformations are required; and (iii) semantic consistency is preserved. To solve these challenges, the paper proposes a unified model to translate images across multiple domains with significant domain gaps. Unlike previous models that constrain the generators with the ubiquitous cycle-consistency constraint to achieve the content similarity, the proposed model employs a perceptual self-regularization constraint. With a single unified generator, the model can maintain consistency over the global shapes as well as the local texture information across multiple domains. Extensive qualitative and quantitative evaluations demonstrate the effectiveness and superior performance over state-of-the-art models. It is more effective in representing shape deformation in challenging mappings with significant dataset variation across multiple domains.


title: Mass and effective spin distributions of binary black hole events: Evidence for a population of hierarchical black hole mergers in active galactic nuclei disks

abstract: The origins of the coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still in debate and clues may present in the mass and effective spin ($\chi_{\rm eff}$) distributions of these merger events. Here we analyze the GWTC-1 and GWTC-2.1 BBH events, and find evidence for a primary mass dependent $\chi_{\rm eff}$ distribution. Below (above) the mass of $\sim 50M_\odot$, a nearly zero (relatively high) effective spin is typical. In the joint analysis of the primary mass and $\chi_{\rm eff}$ distributions, we identify a sub-population of BBHs characterized by both a bump-like primary mass distribution extending to $\sim 100M_\odot$ and a relatively high $\chi_{\rm eff}\sim 0.3$, consistent with the expectation for the $\geq 2$g merger events in the disks of active galactic nuclei (AGNs). Our results suggest that the merger rate in AGN disks is around $1.7~{\rm Gpc}^{-3}{\rm yr}^{-1}$.


title: Efficient, Systematic Estimation of Alloy Free Energy from Special Microscopic States

abstract: For classical discrete systems under constant composition typically refferred to substitutional alloys, we propose calculation method of Helmholtz free energy based on a set of special microscopic states. The advantage of the method is that configuration of the special states are essentially independent of energy and temperature, and they depend only on underlying lattice: The special states can be known a priori without any thermodynamic information, enabling systematic prediction of free energy for multicomponent alloys. We confirm that by comparing to conventional thermodynamic simulation, information about the special states provide reasonable predictive power above order-disorder and phase-separating transition temperature for alloys with many-body (up to 3-body) interactions.


title: An Efficient Scheme for Sampling in Constrained Domains

abstract: The creation of optimal samplers can be a challenging task, especially in the presence of constraints on the support of parameters. One way of mitigating the severity of this challenge is to work with transformed variables, where the support is more conducive to sampling. In this work, a particular transformation called inversion in a sphere is embedded within the popular Metropolis-Hastings paradigm to effectively sample in such scenarios. The method is illustrated on three domains: the standard simplex (sum-to-one constraint), a sector of an $n$-sphere, and hypercubes. The method's performance is assessed using simulation studies with comparisons to strategies from existing literature.


title: New advances in the automation of context-aware information selection and guided model assembly

abstract: The automated assembly and extension of dynamic network models using information extracted from literature are challenging due to the amount and inconsistency in published literature. Recently, efforts have been made to automatically and efficiently assemble the information extracted from literature into models. In this review, we summarize the basic concept, performance, advantages, and limitations of five automated extension methods. Each method was tested for its ability to reconstruct a model of T-cell differentiation as compared against a number of predefined system properties.


title: SMOF: Squeezing More Out of Filters Yields Hardware-Friendly CNN Pruning

abstract: For many years, the family of convolutional neural networks (CNNs) has been a workhorse in deep learning. Recently, many novel CNN structures have been designed to address increasingly challenging tasks. To make them work efficiently on edge devices, researchers have proposed various structured network pruning strategies to reduce their memory and computational cost. However, most of them only focus on reducing the number of filter channels per layer without considering the redundancy within individual filter channels. In this work, we explore pruning from another dimension, the kernel size. We develop a CNN pruning framework called SMOF, which Squeezes More Out of Filters by reducing both kernel size and the number of filter channels. Notably, SMOF is friendly to standard hardware devices without any customized low-level implementations, and the pruning effort by kernel size reduction does not suffer from the fixed-size width constraint in SIMD units of general-purpose processors. The pruned networks can be deployed effortlessly with significant running time reduction. We also support these claims via extensive experiments on various CNN structures and general-purpose processors for mobile devices.


title: Deep Learning Inter-atomic Potential for Thermal and Phonon Behaviour of Silicon Carbide with Quantum Accuracy

abstract: Silicon carbide (SiC) is an essential material for next generation semiconductors and components for nuclear plants. It's applications are strongly dependent on its thermal conductivity, which is highly sensitive to microstructures. Molecular dynamics (MD) simulation is the most used methods to address thermal transportation mechanisms in devices or microstructures of nano-meters. However, the implementation of MD is limited in SiC because of lacking accurate inter-atomic potentials. In this work, using the Deep Potential (DP) methodology, we developed two inter-atomic potentials (DP-IAPs) for SiC based on two adaptively generated datasets within the density functional approximations at the local density and the generalized gradient levels. These two DP-IAPs manifest their speed with quantum accuracy in lattice dynamics simulations as well as scattering rate analysis of phonon transportation. Combining with molecular dynamics simulations, the thermal transport and mechanical properties were systematically investigated. The presented methodology and the inter-atomic potentials pave the way for a systematic approach to model heat transport in SiC related devices using multiscale modelling.


title: Large change of interlayer vibrational coupling with stacking in Mo$_{1-x}$W$_{x}$Te$_{2}$

abstract: Stacking variations in quasi-2D materials can have an important influence on material properties, such as changing the topology of the band structure. Unfortunately, the weakness of van der Waals interactions makes it difficult to compute the stacking dependence of properties, and even in a material as simple as graphite the stacking energetics remain unclear. Mo$_{1-x}$W$_{x}$Te$_{2}$ is a material in which three differently-stacked phases are conveniently accessible by temperature changes: $1T^{\prime}$, $T^*_d$, and the reported Weyl semimetal phase $T_d$. The transitions proceed via layer sliding, and the corresponding interlayer shear mode (ISM) is relevant not just for the stacking energetics, but for understanding the relationship between the Weyl physics and structural changes. However, the interlayer interactions of Mo$_{1-x}$W$_{x}$Te$_{2}$ are not well understood, with wide variation in computed properties. We report inelastic neutron scattering of the ISM in a Mo$_{0.91}$W$_{0.09}$Te$_{2}$ crystal. The ISM energies are generally consistent with the linear chain model (LCM), as expected given the weak interlayer interaction, though there are some discrepancies from predicted intensities. However, the interlayer force constants $K_x$ in the $T^*_d$ and $1T^{\prime}$ phases are substantially weaker than that of $T_d$, at 76(3)% and 83(3)%, respectively. Considering that the relative positioning of atoms in neighboring layers is approximately the same regardless of overall stacking, our results suggest that longer-range influences, such as stacking-induced band structure changes, may be responsible for the substantial change in the interlayer vibrational coupling. These findings should elucidate the stacking energetics of Mo$_{1-x}$W$_{x}$Te$_{2}$ and other van der Waals layered materials.


title: Fast active thermal cloaking through PDE-constrained optimization and reduced-order modeling

abstract: In this paper we show how to efficiently achieve thermal cloaking from a computational standpoint in several virtual scenarios by controlling a distribution of active heat sources. We frame this problem in the setting of PDE-constrained optimization, where the reference field is the solution of the time-dependent heat equation in the absence of the object to cloak. The optimal control problem then aims at actuating the space-time control field so that the thermal field outside the obstacle is indistinguishable from the reference field. In particular, we consider multiple scenarios where material's thermal diffusivity, source intensity and obstacle's temperature are allowed to vary within a user-defined range. To tackle the thermal cloaking problem in a rapid and reliable way, we rely on a parametrized reduced order model built through the reduced basis method, thus entailing huge computational speedups compared to high-fidelity, full-order model exploiting the finite element method while dealing both with complex target shapes and disconnected control domains.


title: Merger rate density of binary black holes through isolated Population I, II, III and extremely metal-poor binary star evolution

abstract: We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitational wave (GW) observations. Notably, Pop III and EMP ($< 10^{-2}$ $Z_\odot$) binary stars yield most of the pair instability (PI) mass gap events with $m_1 = 65$--$130$ $M_\odot$. Pop III binary stars contribute more to the PI mass gap events with increasing redshift, and all the PI mass gap events have the Pop III origin at redshifts $\gtrsim 8$. Our result can be assessed by future GW observations in the following two points. First, there are no binary BHs with $m_1=100$--$130$ $M_\odot$ in our result, and thus the $m_1$ distribution should suddenly drop in the range of $m_1=100$--$130$ $M_\odot$. Second, the PI mass gap event rate should increase toward higher redshift up to $\sim 11$, since those events mainly originate from the Pop III binary stars. We find that the following three assumptions are needed to reproduce the current GW observations: a top-heavy stellar initial mass function and the presence of close binary stars for Pop III and EMP binary stars, and inefficient convective overshoot in the main-sequence phase of stellar evolution. Without any of the above, the number of PI mass gap events becomes too low to reproduce current GW observations.