title: Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

abstract: The flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear stress ratio $\mu$ are relatively insensitive to system size fo a wide range of inertial numbers $I$. Second-order effects in strain rate, such as second normal stress differences, require large system sizes to accurately extract meaningful results, notably a non-monotonic dependence in the first normal stress difference with strain rate. The first-order rheological response represented by the $\mu(I)$ relationship works well at describing the lower-order aspects of the rheology, except near the quasi-static limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems $N\le$ 1,000 near the quasi-static limit can lead to arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasi-static and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasi-static flow.


title: FairEdit: Preserving Fairness in Graph Neural Networks through Greedy Graph Editing

abstract: Graph Neural Networks (GNNs) have proven to excel in predictive modeling tasks where the underlying data is a graph. However, as GNNs are extensively used in human-centered applications, the issue of fairness has arisen. While edge deletion is a common method used to promote fairness in GNNs, it fails to consider when data is inherently missing fair connections. In this work we consider the unexplored method of edge addition, accompanied by deletion, to promote fairness. We propose two model-agnostic algorithms to perform edge editing: a brute force approach and a continuous approximation approach, FairEdit. FairEdit performs efficient edge editing by leveraging gradient information of a fairness loss to find edges that improve fairness. We find that FairEdit outperforms standard training for many data sets and GNN methods, while performing comparably to many state-of-the-art methods, demonstrating FairEdit's ability to improve fairness across many domains and models.


title: Upper bound on the rate of mixing for the Earthquake flow on moduli spaces

abstract: We prove that the earthquake flow is at most polynomially mixing with a degree bounded by a constant depending only on the topology of the surface. In particular it is not exponentially mixing.


title: Continuous gravitational waves in the lab: recovering audio signals with a table-top optical microphone

abstract: Gravitational-wave observatories around the world are searching for continuous waves: persistent signals from sources such as spinning neutron stars. These searches use sophisticated statistical techniques to look for weak signals in noisy data. In this paper, we demonstrate these techniques using a table-top model gravitational-wave detector: a Michelson interferometer where sound is used as an analog for gravitational waves. Using signal processing techniques from continuous-wave searches, we demonstrate the recovery of tones with constant and wandering frequencies. We also explore the use of the interferometer as a teaching tool for educators in physics and electrical engineering by using it as an "optical microphone" to capture music and speech. A range of filtering techniques used to recover signals from noisy data are detailed in the Supplementary Material. Here, we present highlights of our results using a combined notch plus Wiener filter and the statistical log minimum mean-square error (logMMSE) estimator. Using these techniques, we easily recover recordings of simple chords and drums, but complex music and speech are more challenging. This demonstration can be used by educators in undergraduate laboratories and can be adapted for communicating gravitational-wave and signal-processing topics to non-specialist audiences.


title: Gravitational glint: Detectable gravitational wave tails from stars and compact objects

abstract: Observations of a merging neutron star binary in both gravitational waves, by the Laser Interferometer Gravitational-wave Observatory (LIGO), and across the spectrum of electromagnetic radiation, by myriad telescopes, have been used to show that gravitational waves travel in vacuum at a speed that is indistinguishable from that of light to within one part in a quadrillion. However, it has long been expected mathematically that, when electromagnetic or gravitational waves travel through vacuum in a curved spacetime, the waves develop "tails" that travel more slowly. The associated signal has been thought to be undetectably weak. Here we demonstrate that gravitational waves are efficiently scattered by the curvature sourced by ordinary compact objects -- stars, white dwarfs, neutron stars, and planets -- and certain candidates for dark matter, populating the interior of the null cone. The resulting "gravitational glint" should imminently be detectable, and be recognizable (for all but planets) as briefly delayed echoes of the primary signal emanating from extremely near the direction of the primary source. This opens the prospect for using GRAvitational Detection And Ranging (GRADAR) to map the Universe and conduct a comprehensive census of massive compact objects, and ultimately to explore their interiors.


title: Measuring the Galactic Binary Fluxes with LISA: Metamorphoses and Disappearances of White Dwarf Binaries

abstract: The space gravitational wave detector LISA is expected to detect $\sim10^4$ of nearly monochromatic binaries, after $\sim 10$\.yr operation. We propose to measure the inspiral/outspiral binary fluxes in the frequency space, by processing tiny frequency drifts of these numerous binaries. Rich astrophysical information is encoded in the frequency dependencies of the two fluxes, and we can read the long-term evolution of white dwarf binaries, resulting in metamorphoses or disappearances. This measurement will thus help us to deepen our understanding on the strongly interacting exotic objects. Using a simplified model for the frequency drift speeds, we discuss the primary aspects of the flux measurement, including the prospects with LISA.


title: NFANet: A Novel Method for Weakly Supervised Water Extraction from High-Resolution Remote Sensing Imagery

abstract: The use of deep learning for water extraction requires precise pixel-level labels. However, it is very difficult to label high-resolution remote sensing images at the pixel level. Therefore, we study how to utilize point labels to extract water bodies and propose a novel method called the neighbor feature aggregation network (NFANet). Compared with pixellevel labels, point labels are much easier to obtain, but they will lose much information. In this paper, we take advantage of the similarity between the adjacent pixels of a local water-body, and propose a neighbor sampler to resample remote sensing images. Then, the sampled images are sent to the network for feature aggregation. In addition, we use an improved recursive training algorithm to further improve the extraction accuracy, making the water boundary more natural. Furthermore, our method utilizes neighboring features instead of global or local features to learn more representative features. The experimental results show that the proposed NFANet method not only outperforms other studied weakly supervised approaches, but also obtains similar results as the state-of-the-art ones.


title: Nonlocal signatures of hybridization between quantum dot and Andreev bound states

abstract: We investigate local and nonlocal signatures of hybridization between a quantum dot state and an extended Andreev bound state (ABS) in a gate-defined InAs nanowire with multiple side probes. When a quantum dot in one of the side probes was hybridized with an ABS in the nanowire, a characteristic spectroscopic pattern was observed both locally, i.e., in the probe with the quantum dot, and nonlocally, in the tunnel conductance of a remote probe. Nonlocal signatures of hybridization reveal the extended nature of the ABS.


title: The Hybrid Discontinuous Galerkin method for elliptic problems and applications in vertical ocean-slice modeling

abstract: Element Method. The Finite Volume Method guarantees local and global mass conservation. A property not satisfied by the Finite Volume Method. On the down side, the Finite Volume Method requires non trivial modifications to attain high order approximations unlike the Finite Volume Method. It has been contended that the Discontinuous Galerkin Method, locally conservative and high order, is a natural progression for Coastal Ocean Modeling. Consequently, as a primer we consider the vertical ocean-slice model with the inclusion of density effects. To solve these non steady Partial Differential Equations, we develop a pressure projection method for solution. We propose a Hybridized Discontinuous Galerkin solution for the required Poisson Problem in each time step. The purpose, is to reduce the computational cost of classical applications of the Discontinuous Galerkin method. The Hybridized Discontinuous Galerkin method is first presented as a general elliptic problem solver. It is shown that a high order implementation yields fast and accurate approximations on coarse meshes.


title: NuSD: A Geant4 based simulation framework for segmented anti-neutrino detectors

abstract: NuSD: Neutrino Segmented Detector is a Geant4-based user application that simulates inverse beta decay event in a variety of segmented scintillation detectors developed by different international collaborations. This simulation framework uses a combination of cross-programs and libraries including Geant4, ROOT and CLHEP developed and used by high energy physics community. It will enable the neutrino physics community to simulate and study neutrino interactions within different detector concepts using a single program. In addition to neutrino simulations in segmented detectors, this program can also be used for various research projects that use of scintillation detectors for different physics purposes.


title: Ritus functions for graphene-like systems with magnetic fields generated by first-order intertwining operators

abstract: In this work, we construct the exact propagator for Dirac fermions in graphene-like systems immersed in external static magnetic fields with non-trivial spatial dependence. Such field profiles are generated within a first-order supersymmetric framework departing from much simpler (seed) magnetic field examples. The propagator is spanned on the basis of the Ritus eigenfunctions, corresponding to the Dirac fermion asymptotic states in the non-trivial magnetic field background which nevertheless admits a simple diagonal form in momentum space. This strategy enlarges the number of magnetic field profiles in which the fermion propagator can be expressed in a closed-form. Electric charge and current densities are found directly from the corresponding propagator and compared against similar findings derived from other methods.

categories: quant-ph cond-mat.mes-hall math-ph math.MP


title: On-demand Integrated Quantum Memory for Polarization Qubits

abstract: Photonic polarization qubits are widely used in quantum computation and quantum communication due to the robustness in transmission and the easy qubit manipulation. An integrated quantum memory for polarization qubits is a fundamental building block for large-scale integrated quantum networks. However, on-demand storing polarization qubits in an integrated quantum memory is a long-standing challenge due to the anisotropic absorption of solids and the polarization-dependent features of microstructures. Here we demonstrate a reliable on-demand quantum memory for polarization qubits, using a depressed-cladding waveguide fabricated in a 151Eu3+: Y2SiO5 crystal. The site-2 151Eu3+ ions in Y2SiO5 crystal provides a near-uniform absorption for arbitrary polarization states and a new pump sequence is developed to prepare a wideband and enhanced absorption profile. A fidelity of 99.4\pm0.6% is obtained for the qubit storage process with an input of 0.32 photons per pulse, together with a storage bandwidth of 10 MHz. This reliable integrated quantum memory for polarization qubits reveals the potential for use in the construction of integrated quantum networks.


title: On-demand storage of photonic qubits at telecom wavelengths

abstract: Quantum memories at telecom wavelengths are crucial for the construction of large-scale quantum networks based on existing fiber networks. On-demand storage of telecom photonic qubits is a fundamental request for such networking applications but yet to be demonstrated. Here we demonstrate the storage and on-demand retrieval of telecom photonic qubits using a laser-written waveguide fabricated in an $^{167}$Er$^{3+}$:Y$_2$SiO$_5$ crystal. Both ends of the waveguide memory are directly connected with fiber arrays with a fiber-to-fiber efficiency of 51%. Storage fidelity of 98.3(1)% can be obtained for time-bin qubits encoded with single-photon-level coherent pulses, which is far beyond the maximal fidelity that can be achieved with a classical measure and prepare strategy. This device features high reliability, easy scalability and can be directly integrated into fiber networks, which could play an essential role in fiber-based quantum networks.


title: Perturbative expansion of irreversible works in symmetric and asymmetric processes

abstract: The systematic expansion method of the solution of the Fokker-Planck equation is developed by generalizing the formulation proposed in [J. Phys. A50, 325001 (2017)]. Using this method, we obtain a new formula to calculate the mean work perturbatively which is applicable to systems with degeneracy in the eigenvalues of the Fokker-Planck operator. This method enables us to study how the geometrical symmetry affects thermodynamic description of a Brownian particle. To illustrate the application of the derived theory, we consider the Fokker-Planck equation with a two-dimensional harmonic potential. To investigate the effect of symmetry of the potential, we study thermodynamic properties in symmetric and asymmetric deformation processes of the potential: the rotational symmetry of the harmonic potential is held in the former, but it is broken in the latter. Optimized deformations in these processes are defined by minimizing mean works. Comparing these optimized processes, we find that the difference between the symmetric and asymmetric processes is maximized when the deformation time of the potential is given by a critical time which is characterized by the relaxation time of the Fokker-Planck equation. This critical time in the mean work is smaller than that of the change of the mean energy because of the hysteresis effect in the irreversible processes.


title: The Mechanical Alignment of Dust (MAD) I: On the spin-up process of fractal grains by a gas-dust drift

abstract: Context: Aligned dust grains are commonly exploited to probe the magnetic field orientation. However, the exact physical processes that result in a coherent large-scale grain alignment are far from being constrained. Aims: In this work, we aim to investigate the impact of a gas-dust drift leading to a mechanical alignment of dust (MAD) and to dust polarization. Methods: We explore fractal dust aggregates to statistically analyze the average alignment behavior of distinct grain ensembles. The spin-up efficiencies for individual aggregates are determined utilizing MC simulations. These efficiencies are analyzed to identify stable points for the grain alignment in direction of the gas-dust drift and along the magnetic field lines. Finally, the net dust polarization is calculated per grain ensemble. Results: The mechanical spin-up within the CNM is sufficient to drive grains to a stable alignment. A likely mechanical grain alignment is parallel to the drift direction. All grains can align at subsonic conditions. Here, we predict a polarization efficiency in the order of unity for the MAD. A supersonic drift may result in a rapid rotation where dust grains may become rotationally disrupted and the polarization becomes drastically reduced. In the presence of a magnetic field, the drift required for the alignment of elongated grains is roughly one order of magnitude higher compared to the pure MAD. Here, the dust polarization efficiency is 0.6-0.9 indicating that a drift can provide the prerequisites to probe the magnetic field. The alignment is inefficient when the direction of the drift and the field lines are perpendicular. Conclusions: We find that MAD has to be taken into consideration as an alternative driving mechanism where the standard RAT alignment theory fails to account for the full spectrum of available dust polarization observations.


title: A Finite Equivariant Generalization of Motion Planning and Topological Complexity

abstract: This paper explores topological complexity in the finite equivariant setting. We first define and study an equivariant version of Tanaka's combinatorial complexity for finite topological spaces. We explore the relationships between this invariant and several others already discussed in the literature: Farber's topological complexity, Tanaka's combinatorial complexity, and Colman-Grant's equivariant Lusternik-Schnirelmann category. We find bounds for equivariant combinatorial complexity and for the necessary lengths of equivariant combinatorial motion plannings. We show that the equivariant topological complexity of any finite $G$-space is equal to its equivariant combinatorial complexity. We then adapt Gonz\`{a}lez's simplicial complexity to ordered and unordered $G$-simplicial complexes and explore its first properties. Lastly, we show that the equivariant topological complexity of the realization of any ordered $G$-simplicial complex is equal to the equivariant simplicial complexity.


title: Stratified Graph Spectra

abstract: In classic graph signal processing, given a real-valued graph signal, its graph Fourier transform is typically defined as the series of inner products between the signal and each eigenvector of the graph Laplacian. Unfortunately, this definition is not mathematically valid in the cases of vector-valued graph signals which however are typical operands in the state-of-the-art graph learning modeling and analyses. Seeking a generalized transformation decoding the magnitudes of eigencomponents from vector-valued signals is thus the main objective of this paper. Several attempts are explored, and also it is found that performing the transformation at hierarchical levels of adjacency help profile the spectral characteristics of signals more insightfully. The proposed methods are introduced as a new tool assisting on diagnosing and profiling behaviors of graph learning models.


title: Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

abstract: All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity distributions for a sample of $\sim10^5$ stars across the main sequence observed by $\textit{TESS}$, all of which exhibited power-law distributions similar to those observed in the Sun, albeit with varying slopes. Here we investigate the mechanisms required to produce such a distribution of flaring events via direct current energy deposition, in which coronal magnetic fields braid, reconnect, and produce flares. We adopt a topological model for this process which produces a power-law distribution of energetic flaring events. We expand this model to include the Coriolis effect, which we demonstrate produces a shallower distribution of flare energies in stars that rotate more rapidly (corresponding to a weaker decline in occurrence rates toward increasing flare energies). We present tentative evidence for the predicted rotation-power-law index correlation in the observations. We advocate for future observations of stellar flares that would improve our measurements of the power-law exponents, and yield key insights into the underlying dynamo mechanisms that underpin the self-similar flare intensity distributions.


title: Verified Probabilistic Policies for Deep Reinforcement Learning

abstract: Deep reinforcement learning is an increasingly popular technique for synthesising policies to control an agent's interaction with its environment. There is also growing interest in formally verifying that such policies are correct and execute safely. Progress has been made in this area by building on existing work for verification of deep neural networks and of continuous-state dynamical systems. In this paper, we tackle the problem of verifying probabilistic policies for deep reinforcement learning, which are used to, for example, tackle adversarial environments, break symmetries and manage trade-offs. We propose an abstraction approach, based on interval Markov decision processes, that yields probabilistic guarantees on a policy's execution, and present techniques to build and solve these models using abstract interpretation, mixed-integer linear programming, entropy-based refinement and probabilistic model checking. We implement our approach and illustrate its effectiveness on a selection of reinforcement learning benchmarks.


title: Inverse problems for non-linear Schr\"odinger equations with time-dependent coefficients

abstract: The paper studies the inverse problem of reconstructuring the coefficient $\beta(t,x)$ of the non-linear term and the potential $V(t,x)$ of a non-linear Schr\"odinger equation in time-domain, $( i \frac{\partial}{\partial t} + \Delta + V) u + \beta u^2 = f$ in $(0,T)\times M$, where $M\subset \mathbb{R}^n$ is a convex and compact set with smooth boundary. We consider measurements in a neighborhood $\Omega\subset M$ that is a neighborhood of the boundary of $M$ and the source-to-solution map $L_{\beta,V}$ that maps a source $f$ supported in $\Omega\times [0,T]$ to the restriction of the solution $u$ in $\Omega\times [0,T]$. We show that the map $L_{\beta,V}$ uniquely determines the time-dependent potential and the coefficient of the non-linearity, for a second-order non-linear Schr\"odinger equation and for the Gross-Pitaevskii equation, with a cubic non-linear term $\beta |u|^2u$, that is encountered in quantum physics.


title: Quantum activation functions for quantum neural networks

abstract: The field of artificial neural networks is expected to strongly benefit from recent developments of quantum computers. In particular, quantum machine learning, a class of quantum algorithms which exploit qubits for creating trainable neural networks, will provide more power to solve problems such as pattern recognition, clustering and machine learning in general. The building block of feed-forward neural networks consists of one layer of neurons connected to an output neuron that is activated according to an arbitrary activation function. The corresponding learning algorithm goes under the name of Rosenblatt perceptron. Quantum perceptrons with specific activation functions are known, but a general method to realize arbitrary activation functions on a quantum computer is still lacking. Here we fill this gap with a quantum algorithm which is capable to approximate any analytic activation functions to any given order of its power series. Unlike previous proposals providing irreversible measurement--based and simplified activation functions, here we show how to approximate any analytic function to any required accuracy without the need to measure the states encoding the information. Thanks to the generality of this construction, any feed-forward neural network may acquire the universal approximation properties according to Hornik's theorem. Our results recast the science of artificial neural networks in the architecture of gate-model quantum computers.


title: Hybrid Quantum-Classical Unit Commitment

abstract: This paper proposes a hybrid quantum-classical algorithm to solve a fundamental power system problem called unit commitment (UC). The UC problem is decomposed into a quadratic subproblem, a quadratic unconstrained binary optimization (QUBO) subproblem, and an unconstrained quadratic subproblem. A classical optimization solver solves the first and third subproblems, while the QUBO subproblem is solved by a quantum algorithm called quantum approximate optimization algorithm (QAOA). The three subproblems are then coordinated iteratively using a three-block alternating direction method of multipliers algorithm. Using Qiskit on the IBM Q system as the simulation environment, simulation results demonstrate the validity of the proposed algorithm to solve the UC problem.


title: Learning Logic Programs From Noisy Failures

abstract: Inductive Logic Programming (ILP) is a form of machine learning (ML) which in contrast to many other state of the art ML methods typically produces highly interpretable and reusable models. However, many ILP systems lack the ability to naturally learn from any noisy or partially misclassified training data. We introduce the relaxed learning from failures approach to ILP, a noise handling modification of the previously introduced learning from failures (LFF) approach which is incapable of handling noise. We additionally introduce the novel Noisy Popper ILP system which implements this relaxed approach and is a modification of the existing Popper system. Like Popper, Noisy Popper takes a generate-test-constrain loop to search its hypothesis space wherein failed hypotheses are used to construct hypothesis constraints. These constraints are used to prune the hypothesis space, making the hypothesis search more efficient. However, in the relaxed setting, constraints are generated in a more lax fashion as to avoid allowing noisy training data to lead to hypothesis constraints which prune optimal hypotheses. Constraints unique to the relaxed setting are generated via hypothesis comparison. Additional constraints are generated by weighing the accuracy of hypotheses against their sizes to avoid overfitting through an application of the minimum description length. We support this new setting through theoretical proofs as well as experimental results which suggest that Noisy Popper improves the noise handling capabilities of Popper but at the cost of overall runtime efficiency.


title: Riemann Hypothesis for Non-Abelian Zeta Functions of Curves over Finite Fields

abstract: In this paper, we develop some basic techniques towards the Riemann hypothesis for higher rank non-abelian zeta functions of an integral regular projective curve of genus $g$ over a finite field $\mathbb F_q$. As an application of the Riemann hypothesis for these genuine zeta functions, we obtain some explicit bounds on the fundamental non-abelian $\alpha$- and $\beta$-invariants of $X/\mathbb F_q$ in terms of $X$ and $n,\, q$ and $g$: $$\alpha_{X,\mathbb F_q;n}(mn) = \sum_{V}\frac{q^{h^0(X,V)}-1}{\#\mathrm{Aut}(V)} \qquad{\rm and}\qquad \beta_{X,\mathbb F_q;n}(mn ):= \sum_{V}\frac{1}{\#{\mathrm Aut}(V)}\qquad(m\in \mathbb Z)$$ where $V$ runs through all rank $n$ semi-stable $\mathbb F_q$-rational vector bundles on $X$ of degree $mn$. In particular, $$ \prod_{k=1}^{n}\frac{\ \big( \sqrt q^k-1\big)^{2g-1}\ }{(\sqrt q^k+1)}\leq q^{-\binom{n}{2}(g-1)} \beta_{X,\mathbb F_q;n}(0) \leq \prod_{k=1}^{n}\frac{\ \big( 1+\sqrt q^k\big)^{2g-1}\ }{(\sqrt q^k-1)}, $$ Finally, we demonstrate that the related bounds in lower ranks in turn play a central role in establishing the Riemann hypothesis for higher rank zetas.


title: Diffusion in multi-dimensional solids using Forman's combinatorial differential forms

abstract: The formulation of combinatorial differential forms, proposed by Forman for analysis of topological properties of discrete complexes, is extended by defining the operators required for analysis of physical processes dependent on scalar variables. The resulting description is intrinsic, different from the approach known as Discrete Exterior Calculus, because it does not assume the existence of smooth vector fields and forms extrinsic to the discrete complex. In addition, the proposed formulation provides a significant new modelling capability: physical processes may be set to operate differently on cells with different dimensions within a complex. An application of the new method to the heat/diffusion equation is presented to demonstrate how it captures the effect of changing properties of microstructural elements on the macroscopic behavior. The proposed method is applicable to a range of physical problems, including heat, mass and charge diffusion, and flow through porous media.


title: Quantum Measurement Without Collapse

abstract: It is shown that the classical book by von Neumann proposing dynamics of measured systems with "reduction (or collapse) of system's wave packet" contains also hints how to avoid this discontinuity in time evolution of the measured system (hence it is in this point not quite selfconsistent). The possibility of avoiding that collapse is a consequence of the observation that any "human observer" can observe simultaneously just mutually compatible "observables". In the paper it is shown how to describe this fact and its consequences. The proposed interpretation of quantum measurement leads also to trivial solution of the "Schr\"odinger Cat Paradox".


title: From Second-Order Differential Geometry to Stochastic Geometric Mechanics

abstract: Classical geometric mechanics, including the study of symmetries, Lagrangian and Hamiltonian mechanics, and the Hamilton-Jacobi theory, are founded on geometric structures such as jets, symplectic and contact ones. In this paper, we shall use a partly forgotten framework of second-order (or stochastic) differential geometry, developed originally by L. Schwartz and P.-A. Meyer, to construct second-order counterparts of those classical structures. These will allow us to study symmetries of stochastic differential equations (SDEs), to establish stochastic Lagrangian and Hamiltonian mechanics and their key relations with second-order Hamilton-Jacobi-Bellman (HJB) equations. Indeed, stochastic prolongation formulae will be derived to study symmetries of SDEs and mixed-order Cartan symmetries. Stochastic Hamilton's equations will follow from a second-order symplectic structure and canonical transformations will lead to the HJB equation. A stochastic variational problem on Riemannian manifolds will provide a stochastic Euler-Lagrange equation compatible with HJB one and equivalent to the Riemannian version of stochastic Hamilton's equations. A stochastic Noether's theorem will also follow. The inspirational example, along the paper, will be the rich dynamical structure of Schr\"odinger's problem in optimal transport, where the latter is also regarded as a Euclidean version of hydrodynamical interpretation of quantum mechanics.


title: Rate Distortion Theory for Descriptive Statistics

abstract: Rate distortion theory was developed for optimizing lossy compression of data, but it also has a lot of applications in statistics. In this paper we will see how rate distortion theory can be used to analyze a complicated data set involving orientations of early Islamic mosques. The analysis involves testing, identification of outliers, choice of compression rate, calculation of optimal reconstruction points, and assigning "descriptive confidence regions" to the reconstruction points. In this paper the focus will be on the methods, so the integrity of the data set and the interpretation of the results will not be discussed.


title: An Eulerian-based immersed boundary method for particle suspensions with implicit lubrication model

abstract: We describe an immersed boundary method in which the fluid-structure coupling is achieved in an Eulerian framework. The method is an improved extension of the immersed boundary method originally developed by Kajishima et al. [1], which accounts for the inertia of the fictitious fluid inside the particle volume and is thus able to reproduce the behaviour of particles both in the case of neutrally-buoyant objects and in the presence of density difference between the particles and the fluid. The method is capable to handle the presence of multiple suspended objects, i.e., a suspension, by including a soft-sphere normal collision model, while the lubrication correction typically added to similar immersed boundary methods in order to capture the sub-grid unresolved lubrication force is here treated implicitly, i.e., naturally obtained without any explicit expression, thus no additional computation is required. We show that our methodology can successfully reproduce the rheology of a particle suspension in a shear flow up to a dense regime (with a maximum particle volume fraction around 46%) without any additional correction force. The applicability of this methodology is also tested in a turbulent pressure-driven duct flow at high Reynolds number in the presence of non-negligible inertia and non-uniform shear-rate, showing good agreement with experimental measurements.


title: Pavlovian Signalling with General Value Functions in Agent-Agent Temporal Decision Making

abstract: In this paper, we contribute a multi-faceted study into Pavlovian signalling -- a process by which learned, temporally extended predictions made by one agent inform decision-making by another agent. Signalling is intimately connected to time and timing. In service of generating and receiving signals, humans and other animals are known to represent time, determine time since past events, predict the time until a future stimulus, and both recognize and generate patterns that unfold in time. We investigate how different temporal processes impact coordination and signalling between learning agents by introducing a partially observable decision-making domain we call the Frost Hollow. In this domain, a prediction learning agent and a reinforcement learning agent are coupled into a two-part decision-making system that works to acquire sparse reward while avoiding time-conditional hazards. We evaluate two domain variations: machine agents interacting in a seven-state linear walk, and human-machine interaction in a virtual-reality environment. Our results showcase the speed of learning for Pavlovian signalling, the impact that different temporal representations do (and do not) have on agent-agent coordination, and how temporal aliasing impacts agent-agent and human-agent interactions differently. As a main contribution, we establish Pavlovian signalling as a natural bridge between fixed signalling paradigms and fully adaptive communication learning between two agents. We further show how to computationally build this adaptive signalling process out of a fixed signalling process, characterized by fast continual prediction learning and minimal constraints on the nature of the agent receiving signals. Our results therefore suggest an actionable, constructivist path towards communication learning between reinforcement learning agents.


title: Online Changepoint Detection on a Budget

abstract: Changepoints are abrupt variations in the underlying distribution of data. Detecting changes in a data stream is an important problem with many applications. In this paper, we are interested in changepoint detection algorithms which operate in an online setting in the sense that both its storage requirements and worst-case computational complexity per observation are independent of the number of previous observations. We propose an online changepoint detection algorithm for both univariate and multivariate data which compares favorably with offline changepoint detection algorithms while also operating in a strictly more constrained computational model. In addition, we present a simple online hyperparameter auto tuning technique for these algorithms.


title: Robust Trajectory Tracking and Payload Delivery of a Quadrotor Under Multiple State Constraints

abstract: With quadrotors becoming immensely popular in applications such as relief operations, infrastructure maintenance etc., a key control design challenge arises when the quadrotor has to manoeuvre through constrained spaces during various operational scenarios: for example, inspecting a pipeline within predefined velocity and space, dropping relief material at a precise location under tight spaces etc., under the face of parametric uncertainties and external disturbances. To tackle such scenarios, a controller needs to ensure a predefined tracking accuracy so as not to violate the constraints while simultaneously tackling uncertainties and disturbances. However, state-of-the-art controllers dealing with constrained system motion are either not applicable for an underactuated system like quadrotor, or cannot tackle system uncertainties under full state constraints. This work attempts to fill such a gap in literature by designing Barrier Lyapunov Function (BLF) based robust controllers to satisfy multiple state-constraints while simultaneously negotiating parametric uncertainties and external disturbances. The superiority of the BLF control method over a typical unconstrained controller is demonstrated, followed by a robust control design to satisfy position and orientation constraints on quadrotor dynamics. Finally, full state-constraints on a quadrotor(i.e., constraints on the position, orientation, linear velocity and angular velocity) are satisfied with robust control. For each control design, the closed-loop system stability is studied analytically and the efficacy of the design is validated extensively either via realistic simulation scenarios or via experiments performed on a real quadrotor.


title: Optical Variability of "Light-weight" Supermassive Black Holes at a Few Percent Level from ZTF Forced-Photometry Light Curves

abstract: Large time-domain surveys provide a unique opportunity to detect and explore variability of millions of sources on timescales from days to years. Broadband photometric variability can be used as the key selection criteria for weak type-I active galactic nuclei (AGN), when other "direct" confirmation criteria like X-ray or radio emission are unavailable. However, to detect variability of rather weak AGN powered by intermediate-mass black holes, typical sensitivity provided by existing light curve databases is insufficient. Here we present an algorithm for post-processing of light curves for sources with stochastic variability, retrieved from the The Zwicky Transient Facility (ZTF) Forced Photometry service. Using our approach, we can filter out spurious data points related to data reduction artefacts and also eliminate long-term trends related to imperfect photometric calibration. We can now confidently detect the broad-band variability at the 1-3 $\%$ level which can potentially be used as a substitute for expensive X-ray follow-up observations.


title: CVSS Corpus and Massively Multilingual Speech-to-Speech Translation

abstract: We introduce CVSS, a massively multilingual-to-English speech-to-speech translation (S2ST) corpus, covering sentence-level parallel S2ST pairs from 21 languages into English. CVSS is derived from the Common Voice speech corpus and the CoVoST 2 speech-to-text translation (ST) corpus, by synthesizing the translation text from CoVoST 2 into speech using state-of-the-art TTS systems. Two versions of translation speeches are provided: 1) CVSS-C: All the translation speeches are in a single high-quality canonical voice; 2) CVSS-T: The translation speeches are in voices transferred from the corresponding source speeches. In addition, CVSS provides normalized translation text which matches the pronunciation in the translation speech. On each version of CVSS, we built baseline multilingual direct S2ST models and cascade S2ST models, verifying the effectiveness of the corpus. To build strong cascade S2ST baselines, we trained an ST model on CoVoST 2, which outperforms the previous state-of-the-art trained on the corpus without extra data by 5.8 BLEU. Nevertheless, the performance of the direct S2ST models approaches the strong cascade baselines when trained from scratch, and with only 0.1 or 0.7 BLEU difference on ASR transcribed translation when initialized from matching ST models.


title: Deeply virtual Compton scattering cross section at high Bjorken $x_B$

abstract: We report high-precision measurements of the Deeply Virtual Compton Scattering (DVCS) cross section at high values of the Bjorken variable $x_B$. DVCS is sensitive to the Generalized Parton Distributions of the nucleon, which provide a three-dimensional description of its internal constituents. Using the exact analytic expression of the DVCS cross section for all possible polarization states of the initial and final electron and nucleon, and final state photon, we present the first experimental extraction of all four helicity-conserving Compton Form Factors (CFFs) of the nucleon as a function of $x_B$, while systematically including helicity flip amplitudes. In particular, the high accuracy of the present data demonstrates sensitivity to some very poorly known CFFs.


title: An analysis of reconstruction noise from undersampled 4D flow MRI

abstract: Novel Magnetic Resonance (MR) imaging modalities can quantify hemodynamics but require long acquisition times, precluding its widespread use for early diagnosis of cardiovascular disease. To reduce the acquisition times, reconstruction methods from undersampled measurements are routinely used, that leverage representations designed to increase image compressibility. Reconstructed anatomical and hemodynamic images may present visual artifacts. Although some of these artifact are essentially reconstruction errors, and thus a consequence of undersampling, others may be due to measurement noise or the random choice of the sampled frequencies. Said otherwise, a reconstructed image becomes a random variable, and both its bias and its covariance can lead to visual artifacts; the latter leads to spatial correlations that may be misconstrued for visual information. Although the nature of the former has been studied in the literature, the latter has not received as much attention. In this study, we investigate the theoretical properties of the random perturbations arising from the reconstruction process, and perform a number of numerical experiments on simulated and MR aortic flow. Our results show that the correlation length remains limited to two to three pixels when a Gaussian undersampling pattern is combined with recovery algorithms based on $\ell_1$-norm minimization. However, the correlation length may increase significantly for other undersampling patterns, higher undersampling factors (i.e., 8x or 16x compression), and different reconstruction methods.

categories: eess.IV cs.CV physics.med-ph stat.AP stat.ME


title: Floquet-induced localization in long-range many-body systems

abstract: The fate of many-body localization in long-range interacting systems is not fully settled. For instance, the phase boundary between ergodic and many-body localized regimes is still under debate. Here, we use Floquet dynamics which can induce many-body localization in a clean long-range interacting system through spatiotemporal disorder, which are realized by regular operation of random local rotations. The phase diagram has been determined for two types of uniform and nonuniform long-range couplings. Our Floquet mechanism shows more localizing power than conventional static disorder methods as it pushes the phase boundary in favor of the localized phase. Moreover, our comprehensive long-time simulations provide strong support for obtained results based on static analysis.


title: Derivatives-based portfolio decisions. An expected utility insight

abstract: This paper challenges the use of stocks in portfolio construction, instead we demonstrate that Asian derivatives, straddles, or baskets could be more convenient substitutes. Our results are obtained under the assumptions of the Black--Scholes--Merton setting, uncovering a hidden benefit of derivatives that complements their well-known gains for hedging, risk management, and to increase utility in market incompleteness. The new insights are also transferable to more advanced stochastic settings. The analysis relies on the infinite number of optimal choices of derivatives for a maximized expected utility (EUT) agent; we propose risk exposure minimization as an additional optimization criterion inspired by regulations. Working with two assets, for simplicity, we demonstrate that only two derivatives are needed to maximize utility while minimizing risky exposure. In a comparison among one-asset options, e.g. American, European, Asian, Calls and Puts, we demonstrate that the deepest out-of-the-money Asian products available are the best choices to minimize exposure. We also explore optimal selections among straddles, which are better practical choices than out-of-the-money Calls and Puts due to liquidity and rebalancing needs. The optimality of multi-asset derivatives is also considered, establishing that a basket option could be a better choice than one-asset Asian call/put in many realistic situations.


title: Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

abstract: Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(\theta) \propto 1+\cos \theta$, where $\theta$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr


title: A Low-Cost, Highly Customizable Solution for Position Estimation in Modular Robots

abstract: Accurate position sensing is important for state estimation and control in robotics. Reliable and accurate position sensors are usually expensive and difficult to customize. Incorporating them into systems that have very tight volume constraints such as modular robots are particularly difficult. PaintPots are low-cost, reliable, and highly customizable position sensors, but their performance is highly dependent on the manufacturing and calibration process. This paper presents a Kalman filter with a simplified observation model developed to deal with the non-linearity issues that result in the use of low-cost microcontrollers. In addition, a complete solution for the use of PaintPots in a variety of sensing modalities including manufacturing, characterization, and estimation is presented for an example modular robot, SMORES-EP. This solution can be easily adapted to a wide range of applications.


title: Structure and Semantics Preserving Document Representations

abstract: Retrieving relevant documents from a corpus is typically based on the semantic similarity between the document content and query text. The inclusion of structural relationship between documents can benefit the retrieval mechanism by addressing semantic gaps. However, incorporating these relationships requires tractable mechanisms that balance structure with semantics and take advantage of the prevalent pre-train/fine-tune paradigm. We propose here a holistic approach to learning document representations by integrating intra-document content with inter-document relations. Our deep metric learning solution analyzes the complex neighborhood structure in the relationship network to efficiently sample similar/dissimilar document pairs and defines a novel quintuplet loss function that simultaneously encourages document pairs that are semantically relevant to be closer and structurally unrelated to be far apart in the representation space. Furthermore, the separation margins between the documents are varied flexibly to encode the heterogeneity in relationship strengths. The model is fully fine-tunable and natively supports query projection during inference. We demonstrate that it outperforms competing methods on multiple datasets for document retrieval tasks.


title: Physical Properties of 29 sdB+dM Eclipsing Binaries in Zwicky Transient Facility

abstract: The development of large-scale time-domain surveys provides an opportunity to study the physical properties as well as the evolutionary scenario of B-type subdwarfs (sdB) and M-type dwarfs (dM). Here, we obtained 33 sdB+dM eclipsing binaries based on the Zwicky Transient Facility (ZTF) light curves and $Gaia$ early data release 3 (EDR3) parallaxes. By using the PHOEBE code for light curve analysis, we obtain probability distributions for parameters of 29 sdB+dM. $R_1$, $R_2$, and $i$ are well determined, and the average uncertainty of mass ratio $q$ is 0.08. Our parameters are in good agreement with previous works if a typical mass of sdB is assumed. Based on parameters of 29 sdB+dM, we find that both the mass ratio $q$ and the companion's radius $R_2$ decrease with the shortening of the orbital period. For the three sdB+dMs with orbital periods less than 0.075 days, their companions are all brown dwarfs. The masses and radii of the companions satisfy the mass--radius relation for low-mass stars and brown dwarfs. Companions with radii between $0.12R_\odot$ and $0.15R_\odot$ seem to be missing in the observations. As more short-period sdB+dM eclipsing binaries are discovered and classified in the future with ZTF and $Gaia$, we will have more information to constrain the evolutionary ending of sdB+dM.


title: Recent advance in phase transition of vanadium oxide based solar reflectors and the fabrication progress

abstract: Vanadium dioxide (VO2) as a phase-change material controls the transferred heat during phase transition process between metal and insulator states. At temperature above 68C, the rutile structure VO2 keeps the heat out and increases the IR radiation reflectivity, while at the lower temperature the monoclinic structure VO2 acts as the transparent material and increase the transmission radiation. In this paper, we first present the metal-insulator phase transition (MIT) of the VO2 in high and low temperatures. Then we simulate the meta-surface VO2 of metamaterial reflector by Ansys HFSS to show the emittance tunability of the rutile and monoclinic phase of the VO2. In next section, we will review the recent progress in the deposition of thermochromic VO2 on glass and silicon substrate with modifying the pressure of sputtering gases and temperature of the substrate. Finally, we present the results of the in-situ sputtered VOx thin film on thick SiO2 substrate in different combination of oxygen and argon environment by V2O5 target at temperature higher than 300C and then, analyze it with x-ray diffraction (XRD) method. The thermochromic VO2 based metamaterial structures open a new route to the passive energy-efficient optical solar reflector in the past few years.


title: Is FRB 191001 embedded in a supernova remnant?

abstract: Fast radio burst (FRB) 191001 is localised at the spiral arm of a highly star-forming galaxy with an observed dispersion measure (DM) of 507 pc cm$^{-3}$. Subtracting the contributions of the intergalactic medium and our Milky Way Galaxy from the total DM, one gets an excess of around 200 pc cm$^{-3}$, which may have been contributed by the host galaxy of the FRB. It is found in this work that the position of FRB 191001 is consistent with the distribution of supernovae (SNe) in the spiral arm of their parent galaxies. If this event is indeed due to an SN explosion, then, from the analysis of the SN contributions to the excess DM, a core-collapse (CC) channel is preferred over a thermonuclear runaway. For the CC explosion, depending on the density of the surrounding medium, the age of the central engine that powers the radio burst is within a couple of years to a few decades. However, the observed rotation measure of FRB 191001 does not confirm the fact that the radio burst has passed through the remnant of a young SN.


title: Preparation of 3-qubit states

abstract: We will call a pure qubit state real if all its amplitudes are real numbers. We show that any real 3-qubit state can be prepared using $R_y(\theta)$ gates and at most four controlled-$Z$ gates, and we conjecture that four is optimal. We also present an algorithm -- different from the 2008 algorithm given by Znidaric, Giraud and Georgeot -- that prepares any 3-qubit state using local gates and at most three controlled-$Z$ gates. Videos showing how our method works for two- and three-qubit states can be found at https://youtu.be/LIdYSs-rE-o and https://youtu.be/Kne0Vq7gyzQ


title: Design, fabrication, and spectral characterization of TM-polarized metamaterials-based narrowband infrared filter

abstract: Hyperbolic Metamaterials, as a non-magnetic anisotropic artificial structure, show metal properties in one direction and dielectric behavior in orthogonal directions. The proposed hyperbolic metamaterial filter in this project is designed with the metal wire mesh perpendicular to the alternative layers of dielectric materials, keeps TM center wavelength unchanged for the different angle of incident light in MDIR regime. The geometric size of this nanostructure is smaller than the working wavelength and supports big wavevectors due to hyperbolic dispersion. In contrast with conventional Bragg stack, the copper fakir bed makes the transmission properties of the filter the same. For this purpose, the state-of-the-art fabrication methods are required to make such small dimensions in alternative layers of amorphous silicon and silicon dioxide. In this work, first we demonstrate the simulation of Bragg stack with RCWA and finite element methods. Then we focus on our first-time multistep lithography method used to fabricate the filter. Finally, we experimentally verify the optical characteristic of the fabricated filter using Fourier-transform infrared spectroscopy. The experimental and spectrometry data shows that transmission properties of the hyperbolic metamaterial filter remain the same for oblique TM polarized incident light.


title: Cracking the Quantum Scaling Limit with Machine Learned Electron Densities

abstract: A long-standing goal of science is to accurately solve the Schr\"odinger equation for large molecular systems. The poor scaling of current quantum chemistry algorithms on classical computers imposes an effective limit of about a few dozen atoms for which we can calculate molecular electronic structure. We present a machine learning (ML) method to break through this scaling limit and make quantum chemistry calculations of very large systems possible. We show that Euclidean Neural Networks can be trained to predict the electron density with high fidelity from limited data. Learning the electron density allows us to train a machine learning model on small systems and make accurate predictions on large ones. We show that this ML electron density model can break through the quantum scaling limit and calculate the electron density of systems of thousands of atoms with quantum accuracy.


title: Estimation and Inference with Proxy Data and its Genetic Applications

abstract: Existing high-dimensional statistical methods are largely established for analyzing individual-level data. In this work, we study estimation and inference for high-dimensional linear models where we only observe "proxy data", which include the marginal statistics and sample covariance matrix that are computed based on different sets of individuals. We develop a rate optimal method for estimation and inference for the regression coefficient vector and its linear functionals based on the proxy data. Moreover, we show the intrinsic limitations in the proxy-data based inference: the minimax optimal rate for estimation is slower than that in the conventional case where individual data are observed; the power for testing and multiple testing does not go to one as the signal strength goes to infinity. These interesting findings are illustrated through simulation studies and an analysis of a dataset concerning the genetic associations of hindlimb muscle weight in a mouse population.


title: Nonlinear evolution of the magnetorotational instability in eccentric disks

abstract: The magnetorotational instability (MRI) has been extensively studied in circular magnetized disks, and its ability to drive accretion has been demonstrated in a multitude of scenarios. There are reasons to expect eccentric magnetized disks to also exist, but the behavior of the MRI in these disks remains largely uncharted territory. Here we present the first simulations that follow the nonlinear development of the MRI in eccentric disks. We find that the MRI in eccentric disks resembles circular disks in two ways, in the overall level of saturation and in the dependence of the detailed saturated state on magnetic topology. However, in contrast with circular disks, the Maxwell stress in eccentric disks can be negative in some disk sectors, even though the integrated stress is always positive. The angular momentum flux raises the eccentricity of the inner parts of the disk and diminishes the same of the outer parts. Because material accreting onto a black hole from an eccentric orbit possesses more energy than material tracing the innermost stable circular orbit, the radiative efficiency of eccentric disks may be significantly lower than circular disks. This may resolve the "inverse energy problem" seen in many tidal disruption events.


title: Design, fabrication, and spectral characterization of temperature-dependent liquid crystal-based metamaterial to tune dielectric metasurface resonances

abstract: Tunable dielectric meta-surface nanostructures offer incredible performance in optical application due to their extraordinary tunability of the polarization and engineering the dispersion of light with low loss in infrared range. In this article, we designed and experimentally measured the tunability of all-dielectric subwavelength silicon nanoparticles with the help of the temperature-based refractive index of the liquid crystal in the telecom regime. The proposed structure composed of high dielectric nanodisk surrounded by nematic liquid crystal (NLC) is simulated with numerical software, assembled with pre-alignment material, and optically measured by Fourier-transform infrared (FTIR) spectroscopy. The simulated result is compatible with the practical measurements, shows that the tunability of 30nm is achieved. Electric and magnetic resonance modes of the high dielectric nanodisks are tailored in different rates by anisotropic temperature dependent NLC. The phase switching of anisotropic to isotropic nematic liquid crystal enables spectral tunning of the two modes of all dielectric metasurface and modifies the symmetry of the optical response of the metamaterial structure.


title: Heat Kernels and Hardy Spaces on Non-Tangentially Accessible Domains with Applications to Global Regularity of Inhomogeneous Dirichlet Problems

abstract: Let $n\ge2$ and $\Omega$ be a bounded non-tangentially accessible domain (for short, NTA domain) of $\mathbb{R}^n$. Assume that $L_D$ is a second-order divergence form elliptic operator having real-valued, bounded, measurable coefficients on $L^2(\Omega)$ with the Dirichlet boundary condition. The main aim of this article is threefold. First, the authors prove that the heat kernels $\{K_t^{L_D}\}_{t>0}$ generated by $L_D$ are H\"older continuous. Second, for any $p\in(0,1]$, the authors introduce the `geometrical' Hardy space $H^p_r(\Omega)$ by restricting any element of the Hardy space $H^p(\mathbb{R}^n)$ to $\Omega$, and show that, when $p\in(\frac{n}{n+\delta_0},1]$, $H^p_r(\Omega)=H^p(\Omega)=H^p_{L_D}(\Omega)$ with equivalent quasi-norms, where $H^p(\Omega)$ and $H^p_{L_D}(\Omega)$ respectively denote the Hardy space on $\Omega$ and the Hardy space associated with $L_D$, and $\delta_0\in(0,1]$ is the critical index of the H\"older continuity for the kernels $\{K_t^{L_D}\}_{t>0}$. Third, as applications, the authors obtain the global gradient estimates in both $L^p(\Omega)$, with $p\in(1,p_0)$, and $H^p_z(\Omega)$, with $p\in(\frac{n}{n+1},1]$, for the inhomogeneous Dirichlet problem of second-order divergence form elliptic equations on bounded NTA domains, where $p_0\in(2,\infty)$ is a constant depending only on $n$, $\Omega$, and the coefficient matrix of $L_D$. It is worth pointing out that the range $p\in(1,p_0)$ for the global gradient estimate in the scale of Lebesgue spaces $L^p(\Omega)$ is sharp and the above results are established without any additional assumptions on both the coefficient matrix of $L_D$, and the domain $\Omega$.


title: Graph Neural Network Aided Expectation Propagation Detector for MU-MIMO Systems

abstract: Multiuser massive multiple-input multiple-output (MU-MIMO) systems can be used to meet high throughput requirements of 5G and beyond networks. In an uplink MUMIMO system, a base station is serving a large number of users, leading to a strong multi-user interference (MUI). Designing a high performance detector in the presence of a strong MUI is a challenging problem. This work proposes a novel detector based on the concepts of expectation propagation (EP) and graph neural network, referred to as the GEPNet detector, addressing the limitation of the independent Gaussian approximation in EP. The simulation results show that the proposed GEPNet detector significantly outperforms the state-of-the-art MU-MIMO detectors in strong MUI scenarios with equal number of transmit and receive antennas.


title: Right mean for the $\alpha-z$ Bures-Wasserstein quantum divergence

abstract: A new quantum divergence induced from the $\alpha-z$ Renyi relative entropy, called the $\alpha-z$ Bures-Wasserstein quantum divergence, has been recently introduced. We investigate in this paper properties of the right mean, which is a unique minimizer of the weighted sum of $\alpha-z$ Bures-Wasserstein quantum divergences to each points. Many interesting operator inequalities of the right mean with the matrix power mean including the Cartan mean are presented. Moreover, we verify the trace inequality with the Wasserstein mean and provide bounds for the Hadamard product of two right means.


title: The perfectly matched layer (PML) for hyperbolic wave propagation problems: A review

abstract: It is well-known that reliable and efficient domain truncation is crucial to accurate numerical solution of most wave propagation problems. The perfectly matched layer (PML) is a method which, when stable, can provide a domain truncation scheme which is convergent with increasing layer width/damping. The difficulties in using the PML are primarily associated with stability, which can be present at the continuous level or be triggered by numerical approximations. The mathematical and numerical analysis of the PML for hyperbolic wave propagation problems has been an area of active research. It is now possible to construct stable and high order accurate numerical wave solvers by augmenting wave equations with the PML and approximating the equations using summation-by-parts finite difference methods, continuous and discontinuous Galerkin finite element methods. In this review we summarise the progress made, from mathematical, numerical and practical perspectives, point out some open problems and set the stage for future work. We also present numerical experiments of model problems corroborating the theoretical analysis, and numerical simulations of real-world wave propagation demonstrating impact. Stable and parallel implementations of the PML in the high performance computing software packages WaveQLab3D and ExaHyPE allow to sufficiently limit the computational domain of seismological problems with only a few grid points/elements around the computational boundaries where the PML is active, thus saving as much as $96\%$ of the required computational resources for a three space dimensional seismological benchmark problem.


title: Congestion Control Mechanisms for Inter-Datacenter Networks

abstract: Applications running in geographically distributed setting are becoming prevalent. Large-scale online services often share or replicate their data into multiple data centers (DCs) in different geographic regions. Driven by the data communication need of these applications, inter-datacenter network (IDN) is getting increasingly important. However, we find congestion control for inter-datacenter networks quite challenging. Firstly, the inter-datacenter communication involves both data center networks (DCNs) and wide-area networks (WANs) connecting multiple data centers. Such a network environment presents quite heterogeneous characteristics (e.g., buffer depths, RTTs). Existing congestion control mechanisms consider either DCN or WAN congestion, while not simultaneously capturing the degree of congestion for both. Secondly, to reduce evolution cost and improve flexibility, large enterprises have been building and deploying their wide-area routers based on shallow-buffered switching chips. However, with legacy congestion control mechanisms (e.g., TCP Cubic), shallow buffer can easily get overwhelmed by large BDP (bandwidth-delay product) wide-area traffic, leading to high packet losses and degraded throughput. This thesis describes my research efforts on optimizing congestion control mechanisms for the inter-datacenter networks. First, we design GEMINI - a practical congestion control mechanism that simultaneously handles congestions both in DCN andWAN. Second, we present FlashPass - a proactive congestion control mechanism that achieves near zero loss without degrading throughput under the shallow-buffered WAN. Extensive evaluation shows their superior performance over existing congestion control mechanisms.


title: A Combinatorial Proof of the Existence of Dense subsets in R without the "Steinhaus" like Property

abstract: This paper proves the existence of nonmeasurable dense sets with additional properties using combinatorial techniques.


title: Density Functional Theory Transformed into a One-electron Reduced Density Matrix Functional Theory for the Capture of Static Correlation

abstract: Density functional theory (DFT), the most widely adopted method in modern computational chemistry, fails to describe accurately the electronic structure of strongly correlated systems. Here we show that DFT can be formally and practically transformed into a one-electron reduced-density-matrix (1-RDM) functional theory, which can address the limitations of DFT while retaining favorable computational scaling compared to wavefunction-based approaches. In addition to relaxing the idempotency restriction on the 1-RDM in the kinetic energy term, we add a quadratic 1-RDM-based term to DFT's density-based exchange-correlation functional. Our approach, which we implement by quadratic semidefinite programming at DFT's computational scaling of $O(r^{3})$, yields substantial improvements over traditional DFT in the description of static correlation in chemical structures and processes such as singlet biradicals and bond dissociations.


title: D-Graph: AI-Assisted Design Concept Exploration Graph

abstract: We present an AI-assisted search tool, the "Design Concept Exploration Graph" ("D-Graph"). It assists automotive designers in creating an original design-concept phrase, that is, a combination of two adjectives that conveys product aesthetics. D-Graph retrieves adjectives from a ConceptNet knowledge graph as nodes and visualizes them in a dynamically scalable 3D graph as users explore words. The retrieval algorithm helps in finding unique words by ruling out overused words on the basis of word frequency from a large text corpus and words that are too similar between the two in a combination using the cosine similarity from ConceptNet Numberbatch word embeddings. Our experiment with participants in the automotive design field that used both the proposed D-Graph and a baseline tool for design-concept-phrase creation tasks suggested a positive difference in participants' self-evaluation on the phrases they created, though not significant. Experts' evaluations on the phrases did not show significant differences. Negative correlations between the cosine similarity of the two words in a design-concept phrase and the experts' evaluation were significant. Our qualitative analysis suggested the directions for further development of the tool that should help users in adhering to the strategy of creating compound phrases supported by computational linguistic principles.


title: Exact solutions of Kondo problems in higher-order fermions

abstract: The conformal field theory (CFT) approach to Kondo problems, originally developed by Affleck and Ludwig (AL), has greatly advanced the fundamental knowledge of Kondo physics. The CFT approach to Kondo impurities is based on a necessary approximation, i.e., the linearization of the low-lying excitations in a narrow energy window about the Fermi surface. This treatment works well in normal metal baths, but encounters fundamental difficulties in systems with Fermi points and high-order dispersion relations. Prominent examples of such systems are the recently-proposed topological semimetals with emergent higher-order fermions. Here, we develop a new CFT technique that yields exact solutions to the Kondo problems in higher-order fermion systems. Our approach does not require any linearization of the low-lying excitations, and more importantly, it rigorously bosonizes the entire energy spectrum of the higher-order fermions. Therefore, it provides a more solid theoretical base for evaluating the thermodynamic quantities at finite temperatures. Our work significantly broadens the scope of CFT techniques and brings about unprecedented applications beyond the reach of conventional methods.


title: Resonant Precession of Magnetization and Precession -- Induced DC voltages in FeGaB Thin Films

abstract: Measurements of frequency dependent ferromagnetic resonance (FMR) and spin pumping driven dc voltage (V_{dc}) are reported for amorphous films of Fe_{78}Ga_{13}B_{9} (FeGaB) alloy to address the phenomenon of self-induced inverse spin Hall effect (ISHE) in plain films of metallic ferromagnets. The V_{dc} signal, which is antisymmetric on field reversal, comprises of symmetric and asymmetric Lorentzians centered around the resonance field. Dominant role of thin film size effects is seen in setting the magnitude of static magnetization, V_{dc} and dynamics of magnetization precession in thinner films (\leq 8 nm). The film thickness dependence of magnetization parameters indicates the presence of a magnetically disordered region at the film-substrate interface, which may promote preferential flow of spins generated by the precessing magnetization towards the substrate. However, the V_{dc} signal also draws contributions from rectification effects of a \approx 0.4 \% anisotropic magnetoresistance and a large (\approx 54 n\Omega.m) anomalous Hall resistivity (AHR) of these films which ride over the effect of spin-orbit coupling driven spin-to-charge conversion near the film-substrate interface. We have addressed these data in the framework of the existing theories of electrodynamics of a ferromagnetic film subjected to radio-frequency field in a coplanar waveguide geometry. Our estimation of the self-induced ISHE for the sample with 54 n\Omega.m AHR shows that it may contribute significantly (\approx 90\%) to the measured symmetric voltage. This study is expected to be very useful for fully understanding the spin pumping induced dc voltages in metallic ferromagnets with disordered interfaces and large anomalous Hall effect.


title: A Grammar-Based Approach for Applying Visualization Taxonomies to Interaction Logs

abstract: Researchers collect large amounts of user interaction data with the goal of mapping user's workflows and behaviors to their higher-level motivations, intuitions, and goals. Although the visual analytics community has proposed numerous taxonomies to facilitate this mapping process, no formal methods exist for systematically applying these existing theories to user interaction logs. This paper seeks to bridge the gap between visualization task taxonomies and interaction log data by making the taxonomies more actionable for interaction log analysis. To achieve this, we leverage structural parallels between how people express themselves through interactions and language by reformulating existing theories as regular grammars. We represent interactions as terminals within a regular grammar, similar to the role of individual words in a language, and patterns of interactions or non-terminals as regular expressions over these terminals to capture common language patterns. To demonstrate our approach, we generate regular grammars for seven visualization taxonomies and develop code to apply them to three interaction log datasets. In analyzing our results, we find that existing taxonomies at the low-level (i.e., terminals) show mixed results in expressing multiple interaction log datasets, and taxonomies at the high-level (i.e., regular expressions) have limited expressiveness, due to primarily two challenges: inconsistencies in interaction log dataset granularity and structure, and under-expressiveness of certain terminals. Based on our findings, we suggest new research directions for the visualization community for augmenting existing taxonomies, developing new ones, and building better interaction log recording processes to facilitate the data-driven development of user behavior taxonomies.


title: Improving the minimum distance bound of Trace Goppa codes

abstract: In this article we prove that a class of Goppa codes whose Goppa polynomial is of the form $g(x) = x + x^q + \cdots + x^{q^{m-1}}$ where $m \geq 3$ (i.e. $g(x)$ is a trace polynomial from a field extension of degree $m \geq 3$) has a better minimum distance than what the Goppa bound $d \geq 2deg(g(x))+1$ implies. Our improvement is based on finding another Goppa polynomial $h$ such that $C(L,g) = C(M, h)$ but $deg(h) > deg(g)$. This is a significant improvement over Trace Goppa codes over quadratic field extensions (i.e. the case $m = 2$), as the Goppa bound for the quadratic case is sharp.


title: Explaining Predictive Uncertainty by Looking Back at Model Explanations

abstract: Predictive uncertainty estimation of pre-trained language models is an important measure of how likely people can trust their predictions. However, little is known about what makes a model prediction uncertain. Explaining predictive uncertainty is an important complement to explaining prediction labels in helping users understand model decision making and gaining their trust on model predictions, while has been largely ignored in prior works. In this work, we propose to explain the predictive uncertainty of pre-trained language models by extracting uncertain words from existing model explanations. We find the uncertain words are those identified as making negative contributions to prediction labels, while actually explaining the predictive uncertainty. Experiments show that uncertainty explanations are indispensable to explaining models and helping humans understand model prediction behavior.


title: Elastic constant of dielectric nano-thin films using three-layer resonance studied by picosecond ultrasonics

abstract: Elastic constants and sound velocities of nm-order thin films are essential for designing acoustic filters. However, it is difficult to measure them for dielectric thin films. In this study, we use a three-layer structure where a dielectric nano-thin film is sandwiched between thicker metallic films to measure the longitudinal elastic constant of the dielectric film. We propose an efficiency function to estimate the optimal thicknesses of the components. We use Pt/NiO/Pt three-layer films for confirming our proposed method. The determined elastic constant of NiO deposited at room temperature is smaller than the bulk value by $\sim$40$\%$. However, it approaches the bulk value as the deposition temperature increases. We also reveal that uncertainty of the elastic constant of the Pt film insignificantly affects the accuracy of the determined elastic constant of NiO in this structure.


title: AT2019azh: an unusually long-lived, radio-bright thermal tidal disruption event

abstract: Tidal disruption events (TDEs) occur when a star is destroyed by a supermassive black hole at the center of a galaxy, temporarily increasing the accretion rate onto the black hole and producing a bright flare across the electromagnetic spectrum. Radio observations of TDEs trace outflows and jets that may be produced. Radio detections of the outflows from TDEs are uncommon, with only about one third of TDEs discovered to date having published radio detections. Here we present over two years of comprehensive, multi-radio frequency monitoring observations of the tidal disruption event AT2019azh taken with the Very Large Array (VLA) and MeerKAT radio telescopes from approximately 10 days pre-optical peak to 810 days post-optical peak. AT2019azh shows unusual radio emission for a thermal TDE, as it brightened very slowly over two years, and showed fluctuations in the synchrotron energy index of the optically thin synchrotron emission from 450 days post-disruption. Based on the radio properties, we deduce that the outflow in this event is likely non-relativistic and could be explained by a spherical outflow arising from self-stream intersections, or a mildly collimated outflow from accretion onto the supermassive black hole. This data-set provides a significant contribution to the observational database of outflows from TDEs, including the earliest radio detection of a non-relativistic TDE to date, relative to the optical discovery.


title: Performance Bounds for Group Testing With Doubly-Regular Designs

abstract: In the group testing problem, the goal is to identify a subset of defective items within a larger set of items based on tests whose outcomes indicate whether any defective item is present. This problem is relevant in areas such as medical testing, DNA sequencing, and communications. In this paper, we study a doubly-regular design in which the number of tests-per-item and the number of items-per-test are fixed. We analyze the performance of this test design alongside the Definite Defectives (DD) decoding algorithm in several settings, namely, (i) the sub-linear regime $k=o(n)$ with exact recovery, (ii) the linear regime $k=\Theta(n)$ with approximate recovery, and (iii) the size-constrained setting, where the number of items per test is constrained. Under setting (i), we show that our design together with the DD algorithm, matches an existing achievability result for the DD algorithm with the near-constant tests-per-item design, which is known to be asymptotically optimal in broad scaling regimes. Under setting (ii), we provide novel approximate recovery bounds that complement a hardness result regarding exact recovery. Lastly, under setting (iii), we improve on the best known upper and lower bounds in scaling regimes where the maximum test size grows with the total number of items.


title: TSA-Net: Tube Self-Attention Network for Action Quality Assessment

abstract: In recent years, assessing action quality from videos has attracted growing attention in computer vision community and human computer interaction. Most existing approaches usually tackle this problem by directly migrating the model from action recognition tasks, which ignores the intrinsic differences within the feature map such as foreground and background information. To address this issue, we propose a Tube Self-Attention Network (TSA-Net) for action quality assessment (AQA). Specifically, we introduce a single object tracker into AQA and propose the Tube Self-Attention Module (TSA), which can efficiently generate rich spatio-temporal contextual information by adopting sparse feature interactions. The TSA module is embedded in existing video networks to form TSA-Net. Overall, our TSA-Net is with the following merits: 1) High computational efficiency, 2) High flexibility, and 3) The state-of-the art performance. Extensive experiments are conducted on popular action quality assessment datasets including AQA-7 and MTL-AQA. Besides, a dataset named Fall Recognition in Figure Skating (FR-FS) is proposed to explore the basic action assessment in the figure skating scene.


title: Deep Neural Network Approximation For H\"older Functions

abstract: In this work, we explore the approximation capability of deep Rectified Quadratic Unit neural networks for H\"older-regular functions, with respect to the uniform norm. We find that theoretical approximation heavily depends on the selected activation function in the neural network.


title: Ultrahigh quality infrared polaritonic resonators based on bottom-up-synthesized van der Waals nanoribbons

abstract: van der Waals nanomaterials supporting phonon polariton quasiparticles possess unprecedented light confinement capabilities, making them ideal systems for molecular sensing, thermal emission, and subwavelength imaging applications, but they require defect-free crystallinity and nanostructured form factors to fully showcase these capabilities. We introduce bottom-up-synthesized {\alpha}-MoO3 structures as nanoscale phonon polaritonic systems that feature tailorable morphologies and crystal qualities consistent with bulk single crystals. {\alpha}-MoO3 nanoribbons serve as low-loss hyperbolic Fabry-P\'erot nanoresonators, and we experimentally map hyperbolic resonances over four Reststrahlen bands spanning the far- and mid-infrared spectral range, including resonance modes beyond the tenth order. The measured quality factors are the highest from phonon polaritonic van der Waals structures to date. We anticipate that bottom-up-synthesized polaritonic van der Waals nanostructures will serve as an enabling high-performance and low-loss platform for infrared optical and optoelectronic applications.


title: Utilizing Parallelism in Smart Contracts on Decentralized Blockchains by Taming Application-Inherent Conflicts

abstract: Traditional public blockchain systems typically had very limited transaction throughput because of the bottleneck of the consensus protocol itself. With recent advances in consensus technology, the performance limit has been greatly lifted, typically to thousands of transactions per second. With this, transaction execution has become a new performance bottleneck. Exploiting parallelism in transaction execution is a clear and direct way to address this and to further increase transaction throughput. Although some recent literature introduced concurrency control mechanisms to execute smart contract transactions in parallel, the reported speedup that they can achieve is far from ideal. The main reason is that the proposed parallel execution mechanisms cannot effectively deal with the conflicts inherent in many blockchain applications. In this work, we thoroughly study the historical transaction execution traces in Ethereum. We observe that application-inherent conflicts are the major factors that limit the exploitable parallelism during execution. We propose to use partitioned counters and special commutative instructions to break up the application conflict chains in order to maximize the potential speedup. When we evaluated the maximum parallel speedup achievable, these techniques doubled this limit to an 18x overall speedup compared to serial execution, thus approaching the optimum. We also propose OCC-DA, an optimistic concurrency control scheduler with deterministic aborts, which makes it possible to use OCC scheduling in public blockchain settings.


title: Gamma-ray Burst Prompt Emission Spectrum and $E_p$ Evolution Patterns in the ICMART Model

abstract: In this paper, we simulate the gamma-ray bursts (GRBs) prompt emission light curve, spectrum and $E_p$ evolution patterns within the framework of the Internal-Collision-induced MAgnetic Reconnection and Turbulence (ICMART) model. We show that this model can produce a Band shape spectrum, whose parameters ($E_p$, $\alpha$, $\beta$) could distribute in the typical distribution from GRB observations, as long as the magnetic field and the electron acceleration process in the emission region are under appropriate conditions. On the other hand, we show that for one ICMART event, $E_p$ evolution is always a hard-to-soft pattern. However, a GRB light curve is usually composed of multiple ICMART events that are fundamentally driven by the erratic GRB central engine activity. In this case, we find that if one individual broad pulse in the GRB light curve is composed of multiple ICMART events, the overall $E_p$ evolution could be disguised as the intense-tracking pattern. Therefore, mixed $E_p$ evolution patterns can coexist in the same burst, with a variety of combined patterns. Our results support the ICMART model to be a competitive model to explain the main properties of GRB prompt emission. The possible challenges faced by the ICMART model are also discussed in details.


title: Local to global principle over number fields for higher moments

abstract: The local to global principle for densities is a very convenient tool proposed by Poonen and Stoll to compute the density of a given subset of the integers. In this paper we provide an effective criterion to find all higher moments of the density (e.g. the mean, the variance) of a subset of a finite dimensional free module over the ring of algebraic integers of a number field. More precisely, we provide a local to global principle that allows the computation of all higher moments corresponding to the density, over a general number field $K$. This work advances the understanding of local to global principles for density computations in two ways: on one hand, it extends a result of Bright, Browning and Loughran, where they provide the local to global principle for densities over number fields; on the other hand, it extends the recent result on a local to global principle for expected values over the integers to both the ring of algebraic integers and to moments higher than the expected value. To show how effective and applicable our method is, we compute the density, mean and variance of Eisenstein polynomials and shifted Eisenstein polynomials over number fields. This extends (and fully covers) results in the literature that were obtained with ad-hoc methods.


title: Numerical Relativity as a New Tool for Fundamental Cosmology

abstract: Advances in our understanding of the origin, evolution and structure of the universe have long been driven by cosmological perturbation theory, model building and effective field theory. In this review, we introduce numerical relativity as a powerful new complementary tool for fundamental cosmology. To illustrate its power, we discuss applications of numerical relativity to studying the robustness of slow contraction and inflation in homogenizing, isotropizing and flattening the universe beginning from generic unsmooth initial conditions. In particular, we describe how recent numerical relativity studies of slow contraction have revealed a novel, non-linear smoothing mechanism based on ultralocality that challenges the conventional view on what is required to explain the large-scale homogeneity and isotropy of the observable universe.


title: Discovery of Stable Titanium at the Northeastern Jet of Cassiopeia A: Need for a Weak Jet Mechanism?

abstract: The origin of the jet-like structures observed in Cassiopeia A is still unclear, although it seems to be related to its explosion mechanism. X-ray observations of the characteristic structures could provide us useful information on the explosive nucleosynthesis via the observation of elements, which is a unique approach to understand its origin. We here report the discovery of shocked stable Ti, which is produced only at the inner region of exploding stars, in the northeast jet of Cassiopeia A using the 1-Ms deep observation with the Chandra X-ray observatory. The observed Ti coexists with other intermediate-mass elements (e.g. Si, S, Ar, Ca) and Fe at the tip of the X-ray jet structure. We found that its elemental composition is well explained with the production by the incomplete Si burning regime, indicating that the formation process of the jet structure was sub-energetic at the explosion (the peak temperature during the nuclear burning must be $\lesssim$ 5$\times$10$^{9}$ K at most). Thus, we conclude that the energy source that formed the jet structure was not the primary engine for the supernova explosion. Our results are useful to limit the power of the jet-structure formation process, and a weak jet mechanism with low temperature may be needed to explain it.


title: Ab initio calculations of the hyperfine structure of 109Cd, 109Cd+ and reevaluation of the nuclear quadrupole moment Q(109Cd)

abstract: Large scale ab initio calculations of the electric contribution (i.e. the electric field gradient) to the electric quadrupole hyperfine interaction constants $B$ for the $5p~^2P_{3/2}$ state in $^{109}$Cd$^+$ ion and the $5s5p~^3P_{1,2}$ states in neutral $^{109}$Cd atom were performed. To probe the sensitivity of $B$ to different electron correlation effects, six sets of variational multi-configuration Dirac-Hartree-Fock and relativistic configuration interaction calculations employing different strategies were carried out. The calculated electric field gradients, together with experimental values of $B$, allow us to extract a new value of the nuclear electric quadrupole moment $Q$($^{109}$Cd) = 0.634(22) b which is about 5\% larger than the recommended value [Stone, At. Data Nucl. Data Tables 111-112, 1 (2016); Pyykk\"{o}, Mol. Phys. 116, 1328 (2018)]. Efforts were made to provide a realistic theoretical uncertainty for the final $Q$ value based on statistical principles and correlation with the magnetic dipole hyperfine interaction constant $A$.


title: Performance of several Lanczos eigensolvers with HISQ fermions

abstract: We investigate the state-of-the-art Lanczos eigensolvers available in the Grid and QUDA libraries. They include Implicitly Restarted Lanczos, Thick-Restart Lanczos, and Block Lanczos. We measure and analyze their performance for the Highly Improved Staggered Quark (HISQ) Dirac operator. We also discuss optimization of Chebyshev acceleration.


title: Supersolid phases of lattice dipoles tilted in three-dimensions

abstract: By means of quantum Monte Carlo simulations we study phase diagrams of dipolar bosons in a square optical lattice. The dipoles in the system are parallel to each other and their orientation can be fixed in any direction of the three-dimensional space. Starting from experimentally tunable parameters like scattering length and dipolar interaction strength, we derive the parameters entering the effective Hamiltonian. Depending on the direction of the dipoles, various types of supersolids (e.g. checkerboard, stripe) and solids (checkerboard, stripe, diagonal stripe, and an incompressible phase) can be stabilized. Remarkably, we find a cluster supersolid characterized by the formation of horizontal clusters of particles. These clusters order along a direction at an angle with the horizontal. Moreover, we find what we call a grain-boundary superfluid. In this phase, regions with solid order are separated by extended defects -- grain boundaries -- which support superfluidity. We also investigate the robustness of the stripe supersolid against thermal fluctuations. Finally, we comment on the experimental realization of the phases found.


title: New evidence and analysis of cosmological-scale asymmetry in galaxy spin directions

abstract: In the past several decades, multiple cosmological theories that are based on the contention that the Universe has a major axis have been proposed. Such theories can be based on the geometry of the Universe, or multiverse theories such as black hole cosmology. The contention of a cosmological-scale axis is supported by certain evidence such as the dipole axis formed by the CMB distribution. Here I study another form of cosmological-scale axis, based on the distribution of the spin direction of spiral galaxies. Data from four different telescopes is analyzed, showing nearly identical axis profiles when the distribution of the redshifts of the galaxies is similar.


title: Predictive Synthesis of API-Centric Code

abstract: Today's programmers, especially data science practitioners, make heavy use of data-processing libraries (APIs) such as PyTorch, Tensorflow, NumPy, Pandas, and the like. Program synthesizers can provide significant coding assistance to this community of users; however program synthesis also can be slow due to enormous search spaces. In this work, we examine ways in which machine learning can be used to accelerate enumerative program synthesis. We present a deep-learning-based model to predict the sequence of API functions that would be needed to go from a given input to a desired output, both being numeric vectors. Our work is based on two insights. First, it is possible to learn, based on a large number of input-output examples, to predict the likely API function needed in a given situation. Second, and crucially, it is also possible to learn to compose API functions into a sequence, given an input and the desired final output, without explicitly knowing the intermediate values. We show that we can speed up an enumerative program synthesizer by using predictions from our model variants. These speedups significantly outperform previous ways (e.g. DeepCoder) in which researchers have used ML models in enumerative synthesis.


title: A Communication Efficient Quasi-Newton Method for Large-scale Distributed Multi-agent Optimization

abstract: We propose a communication efficient quasi-Newton method for large-scale multi-agent convex composite optimization. We assume the setting of a network of agents that cooperatively solve a global minimization problem with strongly convex local cost functions augmented with a non-smooth convex regularizer. By introducing consensus variables, we obtain a block-diagonal Hessian and thus eliminate the need for additional communication when approximating the objective curvature information. Moreover, we reduce computational costs of existing primal-dual quasi-Newton methods from $\mathcal{O}(d^3)$ to $\mathcal{O}(cd)$ by storing $c$ pairs of vectors of dimension $d$. An asynchronous implementation is presented that removes the need for coordination. Global linear convergence rate in expectation is established, and we demonstrate the merit of our algorithm numerically with real datasets.


title: Quantum criticality and universality in the $p$-wave paired Aubry-Andr\'{e}-Harper model

abstract: We investigate the quantum criticality and universality in Aubry-Andr\'{e}-Harper (AAH) model with $p$-wave superconducting pairing $\Delta$ in terms of the generalized fidelity susceptibility (GFS). We show that the higher-order GFS is more efficient in spotlighting the critical points than lower-order ones, and thus the enhanced sensitivity is propitious for extracting the associated universal information from the finite-size scaling in quasiperiodic systems. The GFS obeys power-law scaling for localization transitions and thus scaling properties of the GFS provide compelling values of critical exponents. Specifically, we demonstrate that the fixed modulation phase $\phi=\pi$ alleviates the odd-even effect of scaling functions across the Aubry-Andr\'{e} transition with $\Delta=0$, while the scaling functions for odd and even numbers of system sizes with a finite $\Delta$ cannot coincide irrespective of the value of $\phi$. A thorough numerical analysis with odd number of system sizes reveals the correlation-length exponent $\nu \simeq 1.000$ and the dynamical exponent $z$ $\simeq$ 1.388 for transitions from the critical phase to the localized phase,suggesting the unusual universality class of localization transitions in the AAH model with a finite $p$-wave superconducting pairing lies in a different universality class from the Aubry-Andr\'{e} transition. The results may be testified in near term state-of-the-art experimental settings.